{"id":"804613dc-8f0d-4ccc-aa53-160fc04250f7","arxiv_id":"2506.10867","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Orientation patterns in growing rod colonies follow the shear rate of an isotropic expansion flow, and n-sided polygonal boundaries are predicted to produce a total topological defect charge of 1 - n/2.","lead":"Growing rod-shaped cells align in patterns that depend strongly on the shape of the container, and this paper shows how to predict those patterns from geometry alone. The same framework lets researchers design container shapes that produce chosen alignment patterns and topological defects.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Isotropic-growth flow field is never directly verified in polygonal simulations; the intermediate link between simulated velocity and predicted uST is missing.","rationale":"The reader correctly identifies the isotropic-growth/no-feedback assumption as the weakest link. I agree this is the load-bearing premise. However, I want to sharpen the concern: even in the regime where the paper claims validation (lmax=2, steady state), the paper never directly tests the operative assumption at the level of the velocity gradient. It computes the theoretical uST from the Poisson equation and compares the resulting orientation to the simulated orientation; it also visually compares simulated and theoretical velocity for one triangle. But it does not extract the actual shear-rate tensor from the simulated velocity field for the polygons that yield Eq. (7). This leaves open the possibility that the orientation agreement is not due to the isotropic-growth uST mechanism. The proposed test--measuring omega and comparing uST in a pentagonal simulation--would close this gap. If the test passes, the 'geometry alone' claim is mechanistically supported; if it fails, the central explanation would need to be revised, even though Eq. (7) might still hold empirically. I therefore keep the verdict at CONDITIONAL (unchanged), but with a sharper condition.","tokens_in":15275,"tokens_out":10318,"duration_ms":119800,"concrete_test":"Run an agent-based simulation of a pentagonal (or hexagonal) open domain with the same parameters as Fig. 2 (lmax=2), time-average the velocity field over at least 10 generations, and compute on a grid the vorticity omega=(grad v - (grad v)^T)/2 and the traceless symmetric shear uST_num. Compare uST_num's principal direction to the theoretical uST from the Poisson solution (Eq. 3 with alpha and absorbing boundaries) using the area-weighted mean absolute angle difference in regions where |uST| is non-negligible. Also compute the ratio ||omega||/||uST_num||. If the angle difference is consistently >10 degrees in the ordered outer layer, or if ||omega||/||uST_num|| is O(1), the isotropic-growth assumption fails and the predicted-orientation agreement is not explained by the proposed mechanism; if the angle difference is small and the vorticity ratio is negligible, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the director is slaved to the traceless shear rate uST of a curl-free, isotropic-growth flow (Eqs. 2-3, 6). This requires that in the agent-based model the actual time-averaged velocity field is (i) irrotational (omega=0) and (ii) its symmetric shear rate matches the predicted uST. The paper verifies (i)-(ii) only qualitatively and only for the triangle: Fig. 2e/f shows a visual match between simulated and theoretical velocity fields. For the full set of polygons used to establish Eq. (7) (n=4 to 17), neither the vorticity nor the simulated uST is ever measured. Consequently, the agreement between the predicted (from the theoretical isotropic-growth uST) and simulated orientation fields could in principle arise from a different effective shear mechanism, leaving the proposed 'geometry-alone' mechanism unvalidated at its operative step. The paper's own Discussion acknowledges that beta and mu vary with geometry and that incompressibility is approximate, which further suggests the velocity field may deviate from the ideal solution. This is not an objection to Eq. (7) as an empirical law, but to the claim that it is explained by the isotropic-growth shear tensor.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that in dense colonies of growing rod-shaped cells, steady-state nematic orientation patterns can be predicted from boundary geometry alone. The key assumption is that the expansion flow is generated by isotropic growth, v = -ζ∇p, so the flow is curl-free and its traceless shear rate tensor uST can be computed from geometry and boundary conditions; the nematic director is then assumed to align with the principal axis of uST. The framework is first applied to radially symmetric geometries, recovering known channel, radial, and inward-growth results, and then to regular polygonal domains where it predicts a total topological defect charge s = 1 - n/2 (Eq. 7). This prediction is compared with agent-based simulations for polygons with n = 3 to n = 17. The paper then extends the model to a quantitative advection-decay equation, v·∇Q = βuST - μQ, whose solution is fitted to radial alignment profiles in ring and wedge geometries and to a channel-like geometry with excess growth in the center. The central claim is that orientation patterns, and hence defect charge, can be forward-engineered through domain geometry without invoking active-stress feedback.","tokens_in":15490,"tokens_out":4199,"duration_ms":54578,"significance":"If the geometry-only mechanism is correct, this is a valuable unifying framework: it links a wide class of previously studied growing-colony geometries to a single scalar Poisson problem and produces an explicit, falsifiable prediction for defect charge in polygons. The paper gives credit to several concrete strengths: Eq. (7) is tested against independent agent-based simulations for a wide range of polygon side counts at aspect ratio 2; the qualitative match between simulated and predicted velocity fields is shown for the triangle (Fig. 2e-f); the advection-decay model reproduces the radial alignment profiles in ring and excess-growth geometries; and the authors openly acknowledge the limitations of their assumptions, including approximate incompressibility and the geometry dependence of fitted parameters. The main significance hinges, however, on the unverified intermediate step for the polygon systems: the actual simulated velocity field is never compared with the isotropic-growth solution for n ≥ 4, so the observed defect-charge agreement is not yet causally tied to the proposed uST mechanism.","major_comments":[{"comment":"The central claim that geometry determines orientation through the isotropic-growth shear tensor requires that the simulated velocity field be, at least approximately, curl-free and equal to the predicted gradient-pressure solution. This is verified only for the equilateral triangle (Fig. 2e-f). For all other polygons used to establish Eq. (7), no simulated velocity field, no vorticity measurement, and no comparison between simulated and predicted uST is reported. The defect-charge agreement could therefore in principle arise from a different shear mechanism with the same boundary-driven topology. Please compute, for at least a representative subset of polygons (e.g., n = 4, 6, 8, 12), the time-averaged agent-based velocity field, its vorticity, and the symmetric shear-rate tensor uST, and compare them with the numerical solutions of Eq. (3) and the resulting uST. This is a load-bearing check for the mechanistic claim.","section":"§IIC, Figs. 2e-f and 3"},{"comment":"The quantitative advection-decay model is fitted with two free parameters β and μ, but the authors report that these parameters vary with ring curvature and with growth excess, and that the fit landscape is degenerate (Supp. Fig. S4). Consequently, the statement that the framework can be extended to 'quantitatively capture alignment strength' and to enable cross-prediction is stronger than the evidence supports: the only transferable parameter combination is qc = αβ/μ, and Fig. 4f shows that this quantity, computed from fits, is lower than the measured channel value at low curvature, with no reported uncertainty. To substantiate the quantitative claim, please either provide a microscopic or independent determination of β and μ, or perform an explicit out-of-sample transfer test (e.g., predict the channel alignment from ring fits and compare with simulation) with confidence intervals. Without this, the model remains a two-parameter interpolation of the radial profiles.","section":"§IIIB, Eqs. (10)-(11), Fig. 4d-g, Supp. Fig. S4"},{"comment":"The authors correctly note that incompressibility is only approximate in the simulations and that density profiles are parabolic, and they restrict their fits to the inner region where density gradients are smaller. However, the polygon defect-charge comparison in Fig. 3 is not accompanied by any estimate of how much the density variation or finite cell hardness affects the validity of the isotropic-growth solution over the region used to compute the charge. Since Eq. (7) is a topological prediction derived from an ideal incompressible flow, please quantify the sensitivity of the computed charge to the observed density profile, for example by repeating the charge computation on a subset of polygons with a stiffness or growth-rate variation, or by showing that the measured velocity deviations remain small over the charge-integration path.","section":"§IV and Supp. Fig. S5"}],"minor_comments":[{"comment":"The phrase 'strongly effected by the system geometry' should read 'strongly affected'.","section":"Introduction"},{"comment":"The sentence 'any anisotropic redistribution induced by the boundary conditions results in nonzero uST = 0' contains a typo; the intended statement should be 'uST ≠ 0'.","section":"§IV Discussion"},{"comment":"The text 'with decay exponent varying with decreasing with increasing R0' is ambiguous; please restate the monotonic trend of the fitted exponent with R0.","section":"§IIIB"},{"comment":"The heading 'TANGENTIAL ALIGNMENT: EXCESS GROWTH' contains a typo; it should read 'TANGENTIAL'.","section":"§IIIC heading"},{"comment":"Example code is stated to be made available 'alongside the final publication'; please clarify the current availability of the simulation and analysis code used to produce the figures, since the model library link is given but the exact analysis scripts are not.","section":"Code availability"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the qualitative defect-charge result is promising, but the missing velocity-field verification for polygons is a genuine load-bearing gap that can be addressed by reanalyzing existing simulations. The quantitative advection-decay model is honestly framed as parameter-dependent, but the wording in the abstract and discussion overstates its predictive power. I therefore recommend major revision rather than rejection: the fixes are concrete and feasible within the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new and testable result here is the defect-charge law s = 1 - n/2 for polygonal domains, confirmed in agent-based simulations up to n = 17 at aspect ratio 2. That is a clean, falsifiable design rule. The isotropic-growth shear-alignment mechanism itself is not new—Langeslay and Juarez already used the shear component of the expansion flow—but the polygon charge law and the advection-decay extension are new. The paper is clearly written, the derivations from the Poisson problem are transparent, and the authors are candid about their model's limitations.\n\nThe soft spots, in proportion. First, the stress-test note is right: the central mechanism is not directly verified where it matters most. For the triangle, they show a qualitative match between simulated and theoretical velocity fields, but for the polygons n = 4 to 17 that establish Eq. (7), they never measure the actual velocity gradient, its vorticity, or the simulated uST. The predicted orientation fields (or at least their defect charge) match, but whether they arise from the isotropic-growth shear tensor is not demonstrated. This is a missing link, not a fatal flaw. Eq. (7) stands as an empirical law, but the \"explained by uST\" part is under-supported. The remedy is straightforward: compute the time-averaged velocity gradient in the simulations and compare with the theoretical curl-free field for all polygons, not just the triangle.\n\nSecond, the advection-decay model is a two-parameter fit with parameters that vary with ring curvature and growth excess. The authors admit this and rescue cross-prediction via the qc ratio, which does tie to the channel measurement. Still, the fit range is restricted after observing mismatches, and the model's predictive content is modest. This is a moderate weakness, largely acknowledged.\n\nThird, the code and data for the specific results are not yet released, so independent reproducibility is hard to assess. The link to the simulation library is there, but the example code is promised \"alongside final publication.\" For a preprint, acceptable, but worth noting.\n\nWho gets value: experimentalists working with microfluidic chambers or growing colonies, and modelers looking for a simple design rule for topological defect charge. The paper deserves a serious referee. My recommendation would be \"major revision,\" with the main request being the velocity-field verification across polygon geometries, plus release of the simulation analysis code.\n\nBest.","headline":"Defect-charge law for polygons is a clean, useful result, but the paper doesn't directly verify that the isotropic-growth shear tensor is the mechanism driving it.","tokens_in":16035,"tokens_out":2866,"would_cite":true,"duration_ms":30786,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In dense growing colonies of rod-shaped cells, the container geometry alone sets the steady-state nematic orientation pattern.","keywords":["growing active matter","nematic order","isotropic growth","shear alignment","topological defect charge","polygonal confinement","agent-based simulation","bacterial colonies"],"falsifier":"Measure the steady velocity field in a growing colony inside a triangular microfluidic chamber and compare it with the curl-free solution of $\\nabla^2 p = -\\alpha/\\zeta$; a significant $\\nabla\\times\\mathbf{v}$, or a time-averaged director that does not rotate by the interior angle at each corner (charge $-1/2$), would refute the central claim. A cheaper test is to simulate division aspect ratio 3 or higher, where the paper already finds the charge law breaking at $n=12$, and check whether the breakpoint tracks the emergence of vorticity in the flow.","tokens_in":15055,"feed_emoji":"🧫","tokens_out":11322,"duration_ms":114849,"temperature":0.7,"pith_summary":"In a dense two-dimensional colony of growing rod-shaped cells, this paper argues, the steady-state pattern of nematic orientation is set by the geometry of the colony boundary rather than by active-stress feedback. The core assumption is that the expansion flow can be approximated as the gradient of a scalar pressure, so geometry and boundary conditions alone determine the velocity field; the traceless shear rate tensor of that field then gives the local preferred orientation. This reproduces previously reported alignment in channels, radial expansion, inward growth, and curved-strain geometries, and it predicts a new, simple law for polygonal domains: in an $n$-sided polygon the director rotates by the interior angle at each corner, giving a net topological defect charge of $s = 1 - n/2$. Agent-based simulations confirm this law for polygons from triangles to 17-gons. A minimal extension that adds advection and a linear decay of order turns the directional prediction into quantitative alignment-strength profiles, including tangential alignment when the inner growth rate is raised.","feed_headline":"Polygonal walls set defect charge in growing colonies","feed_subtitle":"Rods align with shear of isotropic growth flow; every extra corner adds -1/2 net topological charge.","key_machinery":"The load-bearing object is the isotropic-growth approximation combined with the traceless shear rate tensor $\\mathbf{u}^{ST}$. The approximation replaces the full active-nematic feedback loop with $\\mathbf{v} = -\\zeta\\nabla p$ and $\\nabla^2 p = -\\alpha/\\zeta$, so the velocity field is a gradient field with no vorticity and is determined entirely by the growth rate $\\alpha$ and the boundary conditions. The preferred orientation is then read from the principal axis of $\\mathbf{u}^{ST}$, with $\\mathbf{Q}\\approx\\beta\\mathbf{u}^{ST}$ at low alignment; in radial geometries this gives $\\mathbf{u}^{ST} = \\alpha (r/R_0)^{-2}(\\hat{\\mathbf{r}}\\otimes\\hat{\\mathbf{r}} - \\mathbf{I}/2)$. For polygonal domains the same Poisson problem is solved analytically for the equilateral triangle and numerically otherwise, producing the corner-rotation rule and the charge law $s=1-n/2$. The quantitative extension replaces instantaneous alignment with the advection-decay equation $(\\mathbf{v}\\cdot\\nabla)\\mathbf{Q} = \\beta\\mathbf{u}^{ST} - \\mu\\mathbf{Q}$, whose radial solution is governed by the ratio $\\mu/\\alpha$ and explains why the measured alignment decay is slower than the $r^{-2}$ shear profile.","core_discovery":"The central claim is that geometry-induced nematic order in growing rods is slaved to the traceless shear rate tensor $\\mathbf{u}^{ST}$ of an isotropic-growth expansion flow. Writing the flow as $\\mathbf{v} = -\\zeta \\nabla p$ with $\\nabla \\cdot \\mathbf{v} = \\alpha$ makes the velocity field curl-free and computable from the domain shape and boundary conditions alone; at low alignment strength the nematic tensor is taken proportional to $\\mathbf{u}^{ST}$. For regular polygons with absorbing boundaries, solving the Poisson problem yields a director that rotates by the interior angle at each corner, so the total topological charge is $s = 1 - n/2$. The paper verifies this charge law in agent-based simulations of dividing rods for $n = 3$ through $17$, and shows that the same framework, extended to the steady advection-decay equation $(\\mathbf{v}\\cdot\\nabla)\\mathbf{Q} = \\beta \\mathbf{u}^{ST} - \\mu \\mathbf{Q}$, quantitatively reproduces the measured radial alignment profiles in ring and wedge geometries, with deviations only near the outer boundary.","pith_inferences":["Editorial inference: if the corner-rotation law holds in experiments, sequences of polygon corners could be used like optical elements for the director field, writing arbitrary rotations and defect charges into a growing tissue.","Editorial inference: the curl-free flow assumption is the least protected part of the argument; a direct test is to measure vorticity in a growing colony and ask whether the director still tracks $\\mathbf{u}^{ST}$ wherever $\\nabla\\times\\mathbf{v}\\neq 0$.","Editorial inference: the framework suggests that in dense short-rod colonies, orientation is a passive recorder of the expansion geometry; if so, the same Poisson calculation should predict order in smoothly curved or irregular tissue domains without any new physics.","Editorial inference: because the model suppresses feedback, it predicts that average orientation patterns should be identical for colonies with the same geometry but different noise levels, while fluctuations around the average grow with noise; varying division-rate randomization in simulations could test this."],"forward_implications":["Polygonal chambers become a design tool: since the net defect charge is $s = 1 - n/2$, choosing the number of sides prescribes the topological charge of the growing monolayer.","Alignment patterns reported for channels, free radial expansion, inward growth, and curvature-induced strain follow from a single mechanism, the anisotropic part of the growth flow, without invoking active-stress feedback.","The alignment-strength profile decays as a power law whose exponent is controlled by $\\mu/\\alpha$, and since $\\mu$ is found to scale with the growth rate $\\alpha$, the profile is independent of growth rate in the incompressible limit.","Local modifications of channel width can flip the sign of the shear rate tensor and rotate the nematic director by $90^\\circ$ in the expanded region, turning static alignment patterns into prescribed rotational dynamics.","At division aspect ratio 2, the charge law holds on time-averaged fields over the full tested range $n=3$ to $17$; at aspect ratio 3 it starts to deviate at $n=12$, marking the limit of the isotropic-growth picture."],"supporting_citations":[{"why":"Supplies the prior strain-rate alignment model for growing monolayers that this framework generalizes and contrasts.","marker":"[25]"},{"why":"Provides the free-expansion case whose isotropy the theory recovers through $\\mathbf{u}^{ST}=0$.","marker":"[21]"},{"why":"Presents the confinement-induced channel ordering whose interpretation the paper revises via the traceless shear argument.","marker":"[19]"},{"why":"Documents inward-growth orientational order that the radial solution reproduces.","marker":"[24]"},{"why":"Gives the biomechanical channel-ordering model that the advection-decay extension is contrasted with.","marker":"[13]"},{"why":"Reports the reorientation-cascade geometry that the radial family covers.","marker":"[12]"},{"why":"Provides the agent-based model and coarse-grained binning procedure used for the numerical validation.","marker":"[20]"},{"why":"Supplies the analytic solution of the Poisson equation on equilateral triangles used for the triangular flow field.","marker":"[44]"}],"fun_headline_variants":["Polygon shape sets nematic order in growing rods","Geometry dictates alignment patterns in growing colonies","Shear flow from isotropic growth predicts nematic order","Corner count controls topological defect charge","Analytic map from domain shape to rod orientation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The flow field is assumed to come from isotropic growth alone, with no feedback from local cell orientation to the velocity, and the director is assumed to follow the shear rate tensor of that flow; if active-stress feedback or vorticity is significant, the predicted patterns need not hold.","fun_headline_variants_meta":{"raw":{"variants":["Polygon shape sets nematic order in growing rods","Geometry dictates alignment patterns in growing colonies","Shear flow from isotropic growth predicts nematic order","Corner count controls topological defect charge","Analytic map from domain shape to rod orientation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000212,"raw_usage":{"total_tokens":1420,"prompt_tokens":951,"completion_tokens":469,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":400}},"tokens_in":567,"tokens_out":469,"duration_ms":5818,"temperature":1.0,"reasoning_tokens":400,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:15:51.544838+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the steady velocity field in a growing colony inside a triangular microfluidic chamber and compare it with the curl-free solution of $\\nabla^2 p = -\\alpha/\\zeta$; a significant $\\nabla\\times\\mathbf{v}$, or a time-averaged director that does not rotate by the interior angle at each corner (charge $-1/2$), would refute the central claim. A cheaper test is to simulate division aspect ratio 3 or higher, where the paper already finds the charge law breaking at $n=12$, and check whether the breakpoint tracks the emergence of vorticity in the flow.","supporting_citations":[{"cited_title":"Strain rate controls alignment in growing bacterial monolayers","cited_arxiv_id":"2406.09615","evidence_quote":"Supplies the prior strain-rate alignment model for growing monolayers that this framework generalizes and contrasts."},{"cited_title":"Dell’Arciprete, M","cited_arxiv_id":null,"evidence_quote":"Provides the free-expansion case whose isotropy the theory recovers through $\\mathbf{u}^{ST}=0$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the confinement-induced channel ordering whose interpretation the paper revises via the traceless shear argument."},{"cited_title":"Basaran, Y","cited_arxiv_id":null,"evidence_quote":"Documents inward-growth orientational order that the radial solution reproduces."},{"cited_title":"Volfson, S","cited_arxiv_id":null,"evidence_quote":"Gives the biomechanical channel-ordering model that the advection-decay extension is contrasted with."},{"cited_title":"Nijjer, C","cited_arxiv_id":null,"evidence_quote":"Reports the reorientation-cascade geometry that the radial family covers."},{"cited_title":"Isensee, L","cited_arxiv_id":null,"evidence_quote":"Provides the agent-based model and coarse-grained binning procedure used for the numerical validation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the analytic solution of the Poisson equation on equilateral triangles used for the triangular flow field."}],"review_version":1}