{"id":"ec282d62-6a79-41ad-971a-a509f133f663","arxiv_id":"2412.16021","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A curvature-penalized variational model shows that homeostatic stress in growing tissue emerges from minimizing geometric incompatibility, and yields a maximum-size bound for positively curved spheroids.","lead":"This paper proposes that biological tissues grow by balancing the chemical cost of adding mass against the geometric frustration caused by non-uniform growth. It derives a new formula for the target stress tissues seek at rest, which may help predict the final size and internal stress of structures like tumor spheroids.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed equivalence between the second minimization of S(ĝ_B) and the boundary condition R'(ĝ_B)=0 is asserted but not proven; finite-λ size selection rests on it.","rationale":"The reader's weakest_assumption concerned the global-minimum correspondence and open existence/uniqueness. My review narrows this to a specific unproved step: the second-minimization equivalence. The paper's own Section 8 admits the general problem is open, and the numerical demonstration in Fig. 2 is not a proof, because it starts from solutions that already impose R'(ĝ_B)=0. The analytic upper bound (84) is derived in the λ→∞ limit and does not depend on this equivalence, so the central formula for T* (Eq. 3) and the curvature-dominated sphere solutions remain on solid footing. Thus the reader's CONDITIONAL verdict remains appropriate; no change is needed.","tokens_in":24798,"tokens_out":31180,"duration_ms":270079,"concrete_test":"Solve the non-dimensional boundary-value problem (77) for fixed Ŵ*=0.02, λ̂=0.01, k=1, once with the boundary condition R'(ĝ_B)=0 and once with W*=W(ĝ_B) replacing it (with the same other boundary conditions). Compare the resulting ĝ_B and r̂(ĝ_B). If the two procedures give distinct ĝ_B beyond a tight numerical tolerance (e.g., 0.1%), the claimed equivalence fails and the size-selection result for finite λ is unsupported. Repeat for λ̂=0.005 and λ̂=500 to see whether the equivalence only holds in the curvature-dominated limit, where the analytic bound (84) is independent of this issue.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Secs. 6.2 and 7.1 the paper argues that minimizing the reduced action S(ĝ_B) (Eq. 78) over the boundary value ĝ_B is equivalent to imposing R'(ĝ_B)=0, and hence that the homeostatic size corresponds to a global minimum of the growth action. The argument uses S'(ĝ_B)=0 and the fundamental theorem of calculus to conclude W*=W at g=ĝ_B, and then asserts this is equivalent to R'(ĝ_B)=0. But S(ĝ_B) is a shape functional: the integrand contains Γ and r which themselves depend on ĝ_B through the boundary conditions R(ĝ_B)=k and R'(ĝ_B)=0. Therefore dS/dĝ_B is not simply the integrand at the endpoint; it acquires contributions from the variation of the solution with respect to ĝ_B. These contributions vanish only if the ODE system (77) is the full Euler-Lagrange system with natural boundary conditions, which is not demonstrated. Fig. 2 checks the equivalence only for the particular parameters shown and for solutions that already satisfy R'(ĝ_B)=0, so it cannot settle the general claim. Section 8 concedes that existence and uniqueness of minimizers is open and that gauge fixing in lower symmetry is unresolved. Because the finite-λ size predictions of Fig. 3 and the associated claim of a unique homeostatic size depend on this equivalence, this is the load-bearing weak point.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a variational framework for biological growth in which the homeostatic Eshelby stress is not prescribed by hand but derived from a 'growth action' that penalizes deviations of the Ricci curvature of the growth metric from a target value. The central result is the explicit formula T*_AB = W*G_AB + 2λ[(R−R*)R_AB − ¼(R−R*)²G_AB − ∇_A∇_B R + G_AB ΔR] (Eq. 3), obtained as the variational derivative of the action S = ∫_B √G [(W−W*) + (λ/2)(R−R*)²] d³X. The paper then restricts to spherical symmetry, performs a gauge fix to simplify the equilibrium equations, and shows numerically that the moving-boundary problem yields a unique homeostatic size for given W*, R*, and λ. In the limit λ→∞ it recovers the authors' previous uniform-curvature theory and derives an upper bound r ≤ 3.26/√R* for positive target curvature. The paper also discusses why a linear curvature coupling is unsuitable and frames the quadratic coupling as the simplest viable choice.","tokens_in":25137,"tokens_out":17584,"duration_ms":161158,"significance":"If the framework and its spherical reduction are accepted, the paper offers a genuinely useful reformulation: it replaces the six components of an arbitrary, inhomogeneous homeostatic stress tensor with two scalar targets (W*, R*), connects growth mechanics to f(R) gravity and Ricci-flow-like dynamics, and produces falsifiable predictions about size regulation and non-uniform residual-stress distributions that could be tested in cutting experiments. Concrete strengths include the explicit variational derivation of T*_AB, the clean recovery of the λ→∞ limit, and the availability of the Mathematica notebooks used to generate all figures, which supports reproducibility. The main caveat is that the quadratic curvature coupling is a constitutive postulate rather than a derivation from cell-level mechanisms, so the 'geometric nature' of homeostatic stress is established within the model rather than from independent microscale data; nevertheless, the model's predictions are specific enough to be tested against experiments.","major_comments":[{"comment":"The claimed equivalence between the second minimization of S(gB) and the boundary condition R'(gB)=0 is not established. The statement that S'(gB)=0 'by the fundamental theorem of calculus' gives W* = W at g=gB ignores that S(gB) is a shape functional: the fields r, σ, and Γ depend on gB through the ODE system (77) and its boundary conditions, so dS/dgB contains terms involving ∂y/∂gB and boundary variations in addition to the endpoint integrand. The authors should either provide the full endpoint-variation calculation (including the transversality condition and the boundary terms from δr, δσ, δΓ) or present a numerical check that is not circular: compute S(gB) for solutions that do not impose R'(gB)=0, locate its minimum, and compare with the solution of the moving-boundary system that imposes R'(gB)=0, over a range of λ and W*. As presented, Fig. 2 is ambiguous about whether R'(gB)=0 was imposed while generating the plotted solutions; if it was, the agreement is built in, and if it was not, the evidence covers only three λ values, W*=0.02, and two signs of curvature. Because the finite-λ size predictions in Fig. 3 rest on this equivalence, and because Section 8 concedes that existence and uniqueness of minimizers remain open, this is a load-bearing gap.","section":"Sec. 3.3, divergence-free claim"},{"comment":"The assertion that the target stress tensor is divergence-free, ∇^B T*_AB = 0, is stated for the general theory, but it is only true if the chemical potential W* is covariantly constant. Since ∇^A(W* G_AB) = ∇_B W*, an inhomogeneous W* would spoil the conservation law, and the conclusion that 'T*_AB satisfies momentum balance by construction' would fail. The spherical calculations use constant W* and R*, and Section 1 says the targets are 'not required' to be inhomogeneous, but the general claim as written in Section 3.3 is too strong. Please either state the assumption that W* is constant in the general framework or provide the additional terms and consistency conditions that arise for inhomogeneous targets.","section":"Sec. 6.2, Eq. (74), Fig. 2"},{"comment":"The paper's stated conclusions that 'the homeostatic state corresponds to the minimisation of the growth action within the class of spherically symmetric solutions' and that equilibrium yields 'a well-defined equilibrium size' are stronger than what is actually demonstrated. The spherical result is numerical for finite λ, the analytical equivalence behind it is the unproven one discussed above, and Section 8 explicitly leaves open the existence and uniqueness of minimizers as well as gauge fixing in lower symmetry. The abstract and introduction should be tempered to say that the framework offers a new variational route and that, in spherical symmetry, numerical evidence supports unique size selection, rather than presenting these as established general results.","section":"Sec. 8 and Introduction"}],"minor_comments":[{"comment":"The axis labels 'g = g R*+/6' are dimensionally inconsistent and should read 'ĝ = g sqrt(R*+/6)' (or 'g = g sqrt(R*+/6)'), matching the scaling in Eq. (76).","section":"Fig. 2 and Fig. 4 captions"},{"comment":"The caption says the plot shows 'dimensionless equilibrium size r(gB) as a function of the dimensionless compatible growth metric g', but the horizontal axis is the dimensionless chemical potential W*; please correct the wording.","section":"Fig. 3 caption"},{"comment":"The layout of the reduced system (77) is very difficult to parse: boundary conditions are interleaved with the algebraic T_RR = (T*)_RR condition, and it is not immediately clear how many equations, unknowns, and boundary conditions are being counted. Reformatting this as a numbered system with a clean list of boundary conditions would substantially improve readability.","section":"Eq. (77)"},{"comment":"The displayed expression for (T*)_θθ contains a problematic line break and a leading minus sign that make the formula hard to verify; please simplify it, move the full component to an appendix, or use a symbolic shorthand.","section":"Eq. (54)"},{"comment":"Reference [45] (Riccobelli) is incomplete: it lacks year, journal, volume, and pages. Please update it.","section":"References"},{"comment":"The logical relation between the general boundary conditions (37), which already give R = R* and f''(R)∂_C R = 0, and the spherical boundary conditions in (72) should be clarified: in the quadratic case, (37) yields both R(gB)=R* and R'(gB)=0, so the 'second minimization' over gB appears to be an alternative characterization rather than the source of the extra boundary condition. A sentence explaining this relation would prevent confusion.","section":"Sec. 6.2, Eqs. (37) and (72)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the JMPS readership, and the availability of the computational notebooks is a clear asset. The main issue is not the variational derivation itself, which is standard and clean, but the proof or numerical validation of the free-boundary/minimization equivalence that underwrites the finite-λ size predictions. This is fixable within the scope of the manuscript, so I do not recommend rejection. The divergence-free claim for T*_AB with inhomogeneous W* should also be corrected. One additional editorial concern: the validation relies partly on the authors' own previous paper [21], which is appropriate here because the current work is a direct extension, but the new paper would be stronger with an independent experimental or numerical benchmark for the finite-λ regime."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Alexander, here's my read on the Erlich–Zurlo paper.\n\nThe genuinely new thing is a variational principle for growth in which the homeostatic Eshelby stress is not prescribed but derived: penalize deviations of the Ricci scalar of the growth metric from a target R*, and T*_AB emerges as the Euler–Lagrange derivative. This reduces the six arbitrary components of a homeostatic stress field to two scalars, W* and R*, both with a mechanical and geometric meaning. The earlier paper [21] only treated the λ→∞ uniform-curvature limit; the finite-λ extension and the size bound (84) are new and explicit.\n\nThe derivation is mostly sound. The variation of the f(R)-type action is standard and correctly adapted, with boundary terms handled in Appendix A. The spherical reduction and gauge fixing are transparent, the code is provided, and the λ→∞ limit recovers [21] as promised. Section 8 is honest about the open existence and uniqueness question and the unresolved gauge issue in lower symmetry.\n\nThe soft spot is the claimed equivalence between the second minimization of S(g_B) and the boundary condition R'(g_B)=0. The paper says that by the fundamental theorem of calculus S'(g_B)=0 implies W*=W at g_B, which \"is equivalent\" to R'(g_B)=0. That is too quick. S(g_B) is a shape functional: r and Γ inside the integral depend on g_B through the boundary conditions, so dS/dg_B has endpoint terms in addition to the integrand. The stress-test note is right that the paper does not prove those endpoint terms vanish. Appendix A checks only the curvature part, not the elastic and incompressibility contributions. So the global-minimum interpretation of homeostatic equilibrium is not established, and Fig. 2 is a consistency check, not a proof. That said, the finite-λ size predictions in Fig. 3 are obtained by solving the BVP with R'(g_B)=0 directly, so those numerical results do not rest on the unproven equivalence. The conceptual claim about minimization needs either a proof or a softer statement.\n\nAlso worth flagging: the quadratic action is postulated, not derived from cell-level mechanics. That might be fine for a first paper, but the microscopic bridge is still missing. And the experimental support leans on the authors' own [21]; the new finite-λ predictions (e.g., surface-peaked incompatibility) are not yet tested.\n\nOverall: a serious, useful paper for growth mechanics and applied geometry. The derivation is clean, the predictions are concrete, and the limitations are acknowledged. The boundary-equivalence gap should be fixed or explicitly qualified, but it does not sink the paper. I'd send it to a good referee.","headline":"A genuinely new variational route to homeostatic stress in growth mechanics, with a clean reduction to two scalars, but the paper overclaims the proof that equilibrium size equals the action's global minimum.","tokens_in":25585,"tokens_out":5943,"would_cite":true,"duration_ms":54353,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74L15","74B20","53Z05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Homeostatic growth stress is curvature, not a prescribed tensor","keywords":["morphogenesis","nonlinear elasticity","growth laws","incompatibility","homeostatic stress","metric curvature","variational methods","residual stress"],"falsifier":"Grow a spheroid with positive target curvature under increasingly abundant nutrients and measure its terminal radius: if it can exceed $3.26/\\sqrt{R^*}$, the predicted upper bound is wrong. Alternatively, cut a spheroid radially at several distances from the centre: the finite-$\\lambda$ theory predicts opening that grows toward the periphery, whereas a spatially uniform opening would contradict the predicted non-uniform Ricci profile.","tokens_in":24610,"feed_emoji":"🧫","tokens_out":8157,"duration_ms":66468,"temperature":0.7,"pith_summary":"This paper claims that the target stress state which growing tissue tries to reach—the homeostatic Eshelby stress—is not an arbitrary input but the variational derivative of a 'growth action' that penalises deviations from a target amount of geometric frustration. The action combines a chemical potential for adding mass with a quadratic penalty on the Ricci scalar curvature of the growth metric, so that reaching homeostasis means establishing a physiological level of incompatibility rather than hitting a hand-picked stress tensor. The practical gain is that the six independent components of the target stress are replaced by two scalar fields, the chemical potential and the target curvature, and in spherical symmetry the minimisation selects a unique equilibrium size. This matters because it gives a mechanistic and experimentally testable account of how organs regulate their size and store residual stress.","feed_headline":"Homeostatic growth stress is curvature, not a prescribed tensor","feed_subtitle":"A variational action penalising Ricci curvature yields the homeostatic Eshelby stress and a universal upper bound on spheroid size.","key_machinery":"The load-bearing object is the growth action, a volume integral whose density is elastic energy minus chemical potential plus a quadratic penalty on the Ricci scalar $R$ of the growth metric $G_{AB}$. Ricci curvature is the paper's measure of growth incompatibility, the geometric frustration left when locally grown pieces cannot be assembled without residual stress, and the quadratic penalty makes $R=R^*$ the preferred state while avoiding the ill-posed boundary conditions that a linear Einstein–Hilbert coupling would produce. The derivation uses the standard variation of $f(R)$ gravity, including the Laplace–Beltrami terms, and a gauge-fixing diffeomorphism that trades the reference coordinate for the compatible part of the growth metric, reducing the spherical equilibrium to a moving-boundary problem whose extra boundary condition determines the final size.","core_discovery":"The central claim is that the homeostatic Eshelby stress is determined by the growth action $S = \\int_B \\sqrt{G}\\,[(W-W^*) + \\frac{\\lambda}{2}(R-R^*)^2]\\,d^3X$, whose variation yields $T^*_{AB} = W^* G_{AB} + 2\\lambda[(R-R^*)R_{AB} - \\frac{1}{4}(R-R^*)^2 G_{AB} - \\nabla_A\\nabla_B R + G_{AB}\\Delta R]$. Here $W^*$ is the chemical potential of the nutrient bath, $R$ is the Ricci scalar of the growth metric, $R^*$ is the target curvature, and $\\lambda$ is a material modulus. The paper argues that this expression replaces the arbitrary prescription of an inhomogeneous anisotropic homeostatic stress tensor with two scalar fields that carry clear physical meaning, and that the earlier uniform-curvature theory reappears as the $\\lambda\\to\\infty$ limit. For spherical symmetry, solving the equilibrium system with the boundary conditions $R'(0)=R'(g_B)=0$ selects a unique minimising body size, and positive target curvature imposes the explicit upper bound $r \\le 3.26/\\sqrt{R^*}$.","pith_inferences":["If the variational selection of size carries over to lower symmetry, then local cell-level regulation of curvature (for example through area and perimeter targets or cell rearrangements) could robustly set organ size without any global positional information; the paper leaves this connection implicit.","The bound $r \\le 3.26/\\sqrt{R^*}$ is directly testable: measuring the terminal radius of spheroids with different target curvatures would calibrate $R^*$ and could falsify the theory if the bound is violated.","The predicted non-uniform Ricci profile suggests an optimal-cutting protocol: incisions at several radii could reconstruct the curvature field from opening shapes, giving a way to image incompatibility in living tissues.","Because existence and uniqueness of minimisers remain open outside spherical symmetry, the same variational principle may need extra selection criteria in ellipsoidal or layered geometries; numerical experiments there would probe how robust the size-selection mechanism is."],"forward_implications":["To model growth of a tissue, one no longer needs to prescribe six unknown components of target stress; two scalars, $W^*$ and $R^*$, determine the homeostatic state.","Homeostatic size is set by local parameters: positive target curvature gives a hard upper bound on radius, so nutrient supply alone cannot make a spheroid arbitrarily large.","The growth law $\\dot G_{AB} = k(T^*_{AB}-T_{AB})$ becomes a generalised Ricci flow, providing a thermodynamically grounded evolution equation for the growth metric.","Spatially resolved cutting experiments should distinguish the finite-$\\lambda$ regime, where curvature is concentrated near the periphery, from the uniform-curvature $\\lambda\\to\\infty$ regime.","The prior uniform-curvature theory for discs and spheres is recovered as the rigid-curvature limit, so the framework extends that theory rather than replacing it."],"supporting_citations":[{"why":"First recognised compatibility as the source of residual stress in volumetric growth, providing the geometric premise the action formalises.","marker":"[48]"},{"why":"Established the uniform-curvature disc/sphere theory, including size control and cutting patterns, which this paper recovers as $\\lambda\\to\\infty$ and extends.","marker":"[21]"},{"why":"Supplies the variation of $f(R)$ gravity used to derive the bulk target-stress term (28).","marker":"[50]"},{"why":"Supplies the boundary-term calculation in metric $f(R)$ gravity used for the natural boundary conditions.","marker":"[29]"},{"why":"Thermodynamic growth laws with the homeostatic condition $T_{AB}=T^*_{AB}$, the starting point the paper modifies.","marker":"[18]"},{"why":"A chemo-mechanical spheroid model with chemical-potential tensor anisotropy, one of the residual-stress models the new framework is compared with.","marker":"[20]"},{"why":"A poroelastic chemo-mechanical tumour growth model producing residual stress via a diffusing concentration field, a baseline for the size-and-stress predictions.","marker":"[6]"},{"why":"Introduced Ricci flow, which the paper identifies as the form taken by its growth law in the geometric setting.","marker":"[31]"},{"why":"Documents the ill-posed variational problem of linear Einstein–Hilbert action, motivating the quadratic curvature coupling.","marker":"[10]"}],"fun_headline_variants":["Homeostatic growth stress is curvature, not a tensor","Ricci curvature defines homeostatic stress in growth","Growth action penalizes curvature, sets stress naturally","Geometric growth: target curvature bounds spheroid size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a tissue's homeostatic state is the global minimum of the growth action over all admissible growth metrics, with the spherical boundary conditions $R'(0)=R'(g_B)=0$ selecting that minimum uniquely; the paper verifies this numerically for spheres and states in Section 8 that existence and uniqueness of minimisers remains entirely open in general.","fun_headline_variants_meta":{"raw":{"variants":["Homeostatic growth stress is curvature, not a tensor","Ricci curvature defines homeostatic stress in growth","Growth action penalizes curvature, sets stress naturally","Geometric growth: target curvature bounds spheroid size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1588,"prompt_tokens":1110,"completion_tokens":478,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":726,"completion_tokens_details":{"reasoning_tokens":416}},"tokens_in":726,"tokens_out":478,"duration_ms":4646,"temperature":1.0,"reasoning_tokens":416,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:52:56.634942+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Grow a spheroid with positive target curvature under increasingly abundant nutrients and measure its terminal radius: if it can exceed $3.26/\\sqrt{R^*}$, the predicted upper bound is wrong. Alternatively, cut a spheroid radially at several distances from the centre: the finite-$\\lambda$ theory predicts opening that grows toward the periphery, whereas a spatially uniform opening would contradict the predicted non-uniform Ricci profile.","supporting_citations":[{"cited_title":"Compatibility and the genesis of residual stress by volumetric growth","cited_arxiv_id":null,"evidence_quote":"First recognised compatibility as the source of residual stress in volumetric growth, providing the geometric premise the action formalises."},{"cited_title":"Incompatibility-driven growth and size control during development","cited_arxiv_id":null,"evidence_quote":"Established the uniform-curvature disc/sphere theory, including size control and cutting patterns, which this paper recovers as $\\lambda\\to\\infty$ and extends."},{"cited_title":"f (r) theories of gravity","cited_arxiv_id":null,"evidence_quote":"Supplies the variation of $f(R)$ gravity used to derive the bulk target-stress term (28)."},{"cited_title":"Boundary term in metric f (r) gravity: field equations in the metric formalism","cited_arxiv_id":null,"evidence_quote":"Supplies the boundary-term calculation in metric $f(R)$ gravity used for the natural boundary conditions."},{"cited_title":"Thermomechanics of volumetric growth in uniform bodies","cited_arxiv_id":null,"evidence_quote":"Thermodynamic growth laws with the homeostatic condition $T_{AB}=T^*_{AB}$, the starting point the paper modifies."},{"cited_title":"Mechanical feedback in regulating the size of growing multicellular spheroids","cited_arxiv_id":null,"evidence_quote":"A chemo-mechanical spheroid model with chemical-potential tensor anisotropy, one of the residual-stress models the new framework is compared with."},{"cited_title":"Solid tumors are poroelastic solids with a chemo-mechanical feedback on growth","cited_arxiv_id":null,"evidence_quote":"A poroelastic chemo-mechanical tumour growth model producing residual stress via a diffusing concentration field, a baseline for the size-and-stress predictions."},{"cited_title":"Three-manifolds with positive ricci curvature","cited_arxiv_id":null,"evidence_quote":"Introduced Ricci flow, which the paper identifies as the form taken by its growth law in the geometric setting."},{"cited_title":"Boundary terms of the einstein–hilbert action, in: Gravity and the Quantum: Pedagogical Essays on Cosmology, Astrophysics, and Quantum Gravity","cited_arxiv_id":null,"evidence_quote":"Documents the ill-posed variational problem of linear Einstein–Hilbert action, motivating the quadratic curvature coupling."}],"review_version":1}