{"id":"4c452e0b-5625-4853-9a9b-3c0e92f2cdcb","arxiv_id":"1908.01197","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Critical points of a discrete network energy functional satisfy a generalized Murray's law with exponent (gamma+1)/2, and formal continuum analogues hold under extra regularity assumptions.","lead":"This paper proves a generalized version of Murray's law, a scaling rule for branch sizes in transport networks, for a discrete optimization model, and derives formal analogues for two continuum models. It matters because it connects an old empirical law to modern variational models of biological network formation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The continuum Murray laws (31), (38), and (50) are formal: they require a locally Lipschitz flux, a regularity property that is neither proved nor implied by the H1/W1,2γ/(γ-1) solutions constructed in Lemma 3 and Lemma 5.","rationale":"The reader's weakest assumption correctly identifies the missing local Lipschitz regularity of the flux as the load-bearing gap in the continuum claims. The paper itself acknowledges this in Sections 3 and 4, stating that the calculations are formal and that higher regularity is postponed to future work. The discrete theorem is independent and rigorous, so the paper should not be rejected outright; the continuum claims should remain conditional on an unproved regularity statement, exactly as the reader's CONDITIONAL verdict prescribes. No data fitting or circularity is involved, and the self-acknowledged limitation is the decisive factor. My independent reading therefore agrees with the reader's assessment and does not warrant a change of verdict.","tokens_in":13022,"tokens_out":9776,"duration_ms":108215,"concrete_test":"Settle the regularity question analytically for the variational solution of (33): prove or disprove that q_k = (r + nu^{-1/(gamma-1)}|partial_{x_k} p|^{2/(gamma-1)}) partial_{x_k} p belongs to W^{1,infty}_{loc} for smooth data. A tractable first case is d=2, Omega=(0,1)^2, r=1, gamma=2, S=1: derive the explicit PDE -div((1 + nu^{-1}(partial_{x_1}p)^2, 1 + nu^{-1}(partial_{x_2}p)^2) grad p)=1 and check boundedness of all second derivatives of p (or of the flux components) via elliptic estimates. If a counterexample with S in L^2 but q not locally Lipschitz exists, then (31), (38), and (50) fail as theorems in the generality stated; if q is always Lipschitz for smooth S, the remaining gap is to prove the corresponding regularity theorem, without which the formal derivation cannot be accepted as a proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The discrete generalized Murray law (15) is sound, since it follows directly from the derivative identity (11) and the critical-point condition (13). The load-bearing weakness is in the continuum sections. The derivations of (31), (38), and (50) integrate the PDE over an arbitrary C^{0,1} subset Lambda and invoke the De Giorgi-Federer formula, which requires the flux q = -(rI+c)grad p (or q = -(rI+m tensor m)grad p) to be locally Lipschitz continuous on Lambda. The paper explicitly states this requirement in Sections 3 and 4 and says in Section 3.1 that the proof of higher regularity is beyond its scope. However, the weak solutions whose existence is proved in Lemma 3 and Lemma 5 are only in H1 (resp. W1,2γ/(γ-1)), and no C^{1,alpha} or W^{2,infty} estimate is established for p or m. Since the flux is a nonlinear function of grad p (and, in the diffusive variants, of Delta m), weak compactness does not yield well-defined traces of n·q on ∂Lambda. Thus the boundary-flux identities (31), (38), and (50) are not proved as theorems; the abstract's phrase 'prove an analogue' overstates the rigorous content. This is an acknowledged limitation rather than an internal inconsistency, but it is exactly the condition on which the central continuum claims rest. The Section 4 result additionally depends on a cited linear-stability theorem to justify the assumption meas(A)=meas(Omega), and Remark 1 shows that (50) is not the correct law when that assumption fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives Murray's law-type scaling relations from variational models of biological transportation networks. In the discrete setting (Section 2), for a graph with conductivities C_ij and energy E[C] = sum_{edges} (Q_ij^2/C_ij + nu/gamma C_ij^gamma) L_ij subject to Kirchhoff's law, the authors prove, via the derivative identity in Lemma 1, that every critical point satisfies Q_ij^2 = nu C_ij^(gamma+1), and hence the generalized Murray law (15): at each node, sqrt(nu) times the sum of incoming C_ij^((gamma+1)/2) plus the source S_i equals sqrt(nu) times the sum of outgoing C_ij^((gamma+1)/2); for gamma = 1/2 this reduces to the classical 3/4-law. In Section 3, for the continuum limit model (19)--(23), the paper formally derives the Euler-Lagrange relation (24) and, using the De Giorgi-Federer Green formula, the boundary-flux Murray law (31), together with a diffusive variant (38). In Section 4, for the phenomenological model (42)--(44), the paper derives the Euler-Lagrange equation (45) and, assuming linear stability so that m is nonzero almost everywhere, the Murray law (50). The manuscript explicitly labels the continuum derivations as formal and states that the required flux regularity is beyond its scope.","tokens_in":13360,"tokens_out":14052,"duration_ms":142299,"significance":"The discrete part is the strongest contribution: Lemma 1 is proved in detail, and the resulting identity (15) is a clean and nontrivial generalization of Murray's law to arbitrary graphs with multiple branching nodes and a general metabolic exponent gamma. The continuum sections are less rigorous: equations (31), (38), and (50) are conditional identities that hold provided the flux is locally Lipschitz on the integration set, and the weak solutions constructed in Lemmas 3 and 5 do not provide this regularity. As formal derivations the material is useful, and the manuscript is transparent about the gap, but the abstract's wording 'prove an analogue' overstates the rigorous content. The paper also relies on external results for the continuum limit and for the linear-stability criterion; these are supporting ingredients rather than circular assumptions, so no circularity concern arises.","major_comments":[{"comment":"The continuum Murray laws are not established as theorems. The derivations integrate the PDE over an arbitrary C^{0,1} subset Lambda and apply the De Giorgi-Federer formula to the flux q = -(rI+c) grad p (or q = -(rI+m tensor m) grad p), which requires q to be locally Lipschitz continuous on Lambda. Lemma 3 provides only p in H^1, and Lemma 5 provides p in H^1 intersect W^{1,2gamma/(gamma-1)} on A_+ union A_-; neither regularity implies the existence of well-defined traces of n dot q on partial Lambda. The paper acknowledges this in Sections 3 and 4 and in the introduction, but the abstract's claim to 'prove an analogue' is stronger than the content. The authors should either add the missing regularity proof or explicitly present (31), (38), and (50) as conditional formal identities and adjust the abstract accordingly. This is the main load-bearing gap.","section":"Sections 3, 3.1, 4, 4.1; Eqs. (29)--(31), (49)--(50)"},{"comment":"The proof of Lemma 3 asserts that the functional F in (34) is uniformly convex on H^1 and obtains a minimizer p in H^1. However, the term integral |grad p|^{2gamma/(gamma-1)} is not finite for every H^1 function when d >= 3 and 2gamma/(gamma-1) > 2, so F is not defined on all of H^1 and cannot be uniformly convex there. The natural solution space is W^{1,2gamma/(gamma-1)}(Omega), or H^1 intersect W^{1,2gamma/(gamma-1)}, and Lemma 3 should be restated in that space. This is a fixable but nontrivial correction.","section":"Section 3.1, Lemma 3; Eqs. (33)--(34)"},{"comment":"The sentence 'On the set A we have (grad p tensor grad p)m = 0, which is equivalent to m dot grad p = 0' contradicts the Euler-Lagrange equation (grad p tensor grad p)m = -D^2 Delta m + nu |m|^{2(gamma-1)}m and the definition of A. On the set A the scalar product of the right-hand side with m is strictly positive, so (grad p dot m)^2 > 0 and hence (grad p tensor grad p)m cannot vanish. The statement should presumably refer to the complement A^c, where indeed m dot grad p = 0. As written this is an internal inconsistency in the derivation of the diffusive Murray law and needs correction.","section":"Section 4.2, around Eq. (58)"}],"minor_comments":[{"comment":"The text reads 'the formal continuum limit of of the discrete model'; 'of of' should be 'of'.","section":"Section 3, first paragraph"},{"comment":"The phrase 'with the flux := -(rI + m tensor m) grad p' is missing the symbol q before ':='.","section":"Section 4, before Eq. (50)"},{"comment":"In the displayed formula of Remark 1, the last integral is printed as + int_{A^c cap partial Lambda_+} r |grad p dot n| ds, but to match the outflow term it should presumably be over A^c cap partial Lambda_-; otherwise the same boundary segment appears twice with the same sign.","section":"Remark 1"},{"comment":"In equation (38), the characteristic functions chi_{A_+} and chi_{A_-} should carry the subscript k to match the sets A_k^+ and A_k^- defined in the previous paragraph, since the sets depend on the coordinate direction k.","section":"Section 3.2, Eq. (38)"},{"comment":"The symbol 'c2' in equations (53), (54), and (56) is used without definition; it appears to encode the constraint |grad p|^2 <= nu (for instance c2 = 1/nu), and it should be defined explicitly or removed in favor of nu.","section":"Section 4.1, Eqs. (53)--(56)"},{"comment":"The expression (nu |m|^{2gamma} - D^2 m dot Delta m)^{-1/2} is only meaningful where the radicand is positive; the paper should state explicitly that this is precisely the definition of A and explain how the formula is interpreted on A^c, as it later does in the final boundary integral.","section":"Section 4.2, Eq. (58)"},{"comment":"The introduction's statement that the continuum calculations are formal is useful and should be mirrored in the abstract, which currently says the paper 'prove[s] an analogue' for the continuum systems.","section":"Abstract and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The discrete part of the paper is rigorous and publishable, and the formal continuum derivations are clearly presented with an honest statement of the missing regularity. The main reason for major revision rather than rejection is that the gap is explicitly acknowledged and can be addressed by reframing the continuum claims as conditional formal identities, or by adding the missing regularity proof, rather than by discarding the paper's contribution. I would advise the editor that the abstract must be brought in line with the actual theorems before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the discrete half of this paper is solid and worth knowing; the continuum half is a formal derivation that the authors themselves flag, and the abstract oversells it. The generalized Murray law for arbitrary metabolic exponent gamma and multiple branching nodes, equation (15), is a clean consequence of Lemma 1's sensitivity identity. That identity is restated from the authors' earlier work but the proof is given in the paper, so the discrete derivation is self-contained. For gamma=1/2 it reproduces the classical 3/4 law. That is a real, modest contribution.\n\nThe continuum sections are where I part company with the abstract. Equations (31), (38), (50) are obtained by integrating the PDE over a C^{0,1} subset and applying the De Giorgi-Federer formula. That requires the flux to be locally Lipschitz. The solutions constructed in Lemmas 3 and 5 are only H^1 (or W^{1,2gamma/(gamma-1)}), and the paper says in Section 3.1 that higher regularity is beyond its scope. So the continuum 'Murray laws' are conditional on regularity that is not established. To be fair, the paper states this explicitly in Sections 3 and 4.1; it is an acknowledged limitation, not a hidden one. But the abstract's phrase 'prove an analogue' overstates the rigorous content. The Section 4 result also leans on a cited linear stability theorem to assume meas(A)=meas(Omega), and Remark 1 shows the law changes if that fails, so that part is doubly conditional.\n\nI don't see a circularity problem: they start from the energy functional and derive Murray's law as a first-order condition; the target law is not baked into the assumptions. The reliance on self-citations for Lemma 1 and the continuum limit is fine because those are supporting results, not the target conclusion.\n\nThe paper is worth refereeing. The discrete theorem should be accepted; the continuum claims need either a regularity proof or honest rewording to 'formal analogue.' A serious referee could reasonably ask for the abstract to match the body. I'd be happy to cite the discrete result if I worked on network optimization.","headline":"Solid discrete generalization of Murray's law, with continuum analogues that are formally derived and explicitly conditional on unproved flux regularity; the abstract oversells the rigorous content.","tokens_in":13862,"tokens_out":1656,"would_cite":true,"duration_ms":16523,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C35","05C21","76S05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that Murray's law—the classical scaling relation between branch conductivity and metabolic cost—holds at every critical point of a discrete transport-network energy and carries over to continuum and diffusion-augmented…","keywords":["Murray's law","biological transportation networks","continuum limit","energy minimization","conductivity scaling","Kirchhoff law","metabolic cost","3/4-law"],"falsifier":"One concrete test: construct, for $\\gamma>1$ and a smooth source $S\\in L^2(\\Omega)$, the unique minimizer $p$ of (34) (or of (52)) on a rectangle, and compute the flux $q=-(rI+c)\\nabla p$; if $q$ has an interior discontinuity or cusp where the pressure gradient changes sign, the Green-formula step behind (31), (38), and (50) is not justified for that solution, so the continuum Murray law would not follow for the constructed solution.","tokens_in":12776,"feed_emoji":"🌿","tokens_out":16724,"duration_ms":122670,"temperature":0.7,"pith_summary":"This paper shows that Murray's law—the rule that branch radii or conductivities in a transport network scale to balance pumping cost against metabolic upkeep—follows from energy minimization in a general discrete network model and survives passage to continuum descriptions. In the discrete setting, every critical point of the total energy satisfies a generalized relation: at each node, the sum of incoming conductivities raised to the power $(\\gamma+1)/2$, plus the scaled source strength, equals the same sum of outgoing conductivities. For blood vessels ($\\gamma=1/2$) this reduces to the classical $3/4$-power law. The same balance is derived, at least formally, for two families of continuum models—one obtained as the continuum limit of the discrete model on a rectangular mesh, and one phenomenological model of network formation under Darcy flow—and for their diffusion-augmented versions. The continuum results require the flux to be locally Lipschitz continuous, a regularity property the paper does not prove, leaving that step as the main open gap.","feed_headline":"Murray's law holds in discrete and continuum network models","feed_subtitle":"The generalized 3/4-law ties branch conductivities to metabolic cost and holds at every critical point.","key_machinery":"The load-bearing object is the derivative identity (11): $\\partial_{C_{kl}}\\sum_{(i,j)\\in E}(Q_{ij}^2/C_{ij})L_{ij}=-Q_{kl}^2/C_{kl}^2 L_{kl}$, which holds because the Kirchhoff-constrained pressures respond to conductivity changes so that all cross terms cancel. It turns the stationarity condition into the pointwise relation $Q^2=\\nu C^{\\gamma+1}$, and Kirchhoff's law then converts the flux balance at each node into the conductivity-power identity (15). In the continuum analogues, the same role is played by the Euler–Lagrange equations (24) and (45), which tie pressure gradients to powers of the conductivity field, and by the De Giorgi–Federer Green formula, which transplants the node balance into a boundary-flux balance over arbitrary regions $\\Lambda$ (equations (31), (38), (50)).","core_discovery":"The central claim is that Murray's law is not a special feature of single tubes or symmetric trees but a generic consequence of minimizing the energy $E[C]=\\sum_{(i,j)\\in E}(Q_{ij}^2/C_{ij} + \\nu C_{ij}^\\gamma/\\gamma)L_{ij}$ under Kirchhoff's law. At any critical point one has $Q_{ij}^2=\\nu C_{ij}^{\\gamma+1}$ on every edge, and Kirchhoff balance then yields the generalized Murray law (15), $\\sum_{j\\in N^+(i)}C_{ij}^{(\\gamma+1)/2}+S_i/\\sqrt{\\nu}=\\sum_{j\\in N^-(i)}C_{ij}^{(\\gamma+1)/2}$. The paper further claims the analogous balance holds in the continuum: for the diagonal-tensor model (19)–(23) the Euler–Lagrange relation $(\\partial_{x_k}p)^2=\\nu|c_k|^{\\gamma-1}$ converts the integrated Poisson equation over any region $\\Lambda$ into the boundary-flux identity (31), and the same structure yields (50) for the $m\\otimes m$ model when restricted to linearly stable steady states. These derivations are formal, justified only under the assumption that the flux is locally Lipschitz continuous.","pith_inferences":["The same variational argument likely yields Murray-type relations for other cost exponents and graph ensembles, since only (11), stationarity, and Kirchhoff's law are used; one testable extension is to random graphs with random source placements, where (15) should hold at every node individually.","The regularity gap suggests a concrete numerical experiment: compute minimizers of (34) and measure the Hölder exponent of the flux near points where $\\partial_{x_k}p$ changes sign; if the exponent is strictly below 1, the Green-formula step fails and a different proof or a counterexample would be needed.","For $\\gamma=1$, the free-boundary structure of (54)–(56) means the Murray law might serve as a sharp interface condition: the set where $|\\nabla p|^2=\\nu$ is the network, and (50) becomes a transmission condition across that interface—an interpretation the paper does not develop.","The discrete law (15) has a diagnostic corollary the authors do not state: given measured conductivities and sources in a real vascular network, the residual at each node quantifies how far that network is from optimality, offering a statistical test of Murray's law in vivo."],"forward_implications":["In any network that minimizes (10) under Kirchhoff balance, the generalized Murray law (15) holds at every node, so optimal conductivities obey a universal power-law relation regardless of graph topology.","Setting $\\gamma=1/2$ recovers the classical $3/4$-law $C_0^{3/4}=\\sum_i C_i^{3/4}$ for laminar blood flow, placing the classical observation inside a general variational framework.","The continuum limit of the discrete model preserves the Murray balance in the form (31): for any region $\\Lambda$, the difference of in- and out-fluxes through $\\partial\\Lambda$ equals the integrated source strength.","For the phenomenological $m\\otimes m$ model, only linearly stable steady states are admissible, and on those states the Murray balance (50) holds; unstable states with $m=0$ on a positive-measure set require the modified law in Remark 1.","Adding diffusion to either continuum model does not destroy the Murray balance; it only replaces the power $|c_k|^{\\gamma-1}$ (or the analogous $m$-expression) by $\\nu|c_k|^{\\gamma-1}-D^2\\Delta c_k$, so the balance persists in (38) and (58)."],"supporting_citations":[{"why":"It derives the derivative identity (11) and the continuum-limit energy (16) that this paper restates and builds on.","marker":"[12]"},{"why":"It supplies the rigorous continuum limit of the discrete model, giving the rectangular-mesh PDE model (19)–(23).","marker":"[13]"},{"why":"It is Murray's original minimum-work principle, the source of the cost balance (2) and the classical radius-cubed law.","marker":"[14]"},{"why":"It is Murray's companion paper on oxygen exchange, the second origin of the classical 3/4-law.","marker":"[15]"},{"why":"It gives the laminar-flow relation that conductivity is proportional to radius to the fourth power, used to state the law in terms of $C^{3/4}$.","marker":"[17]"},{"why":"It provides the De Giorgi–Federer Green formula that converts the integrated PDE into the boundary-flux Murray laws (31), (38), (50).","marker":"[7]"},{"why":"It establishes the linear-stability criterion meas(A)=meas(Ω) and constructs weak solutions of (52) by direct methods.","marker":"[8]"},{"why":"It develops the γ=1 free-boundary reformulation and the cutoff construction for 1/2≤γ<1 stationary solutions.","marker":"[9]"},{"why":"It derives the phenomenological continuum model (42)–(44) and its Darcy/Kirchhoff structure from discrete-network considerations.","marker":"[2]"}],"fun_headline_variants":["Murray's law holds across discrete and continuous models","Generalized Murray's law derived for biological networks","Branching networks follow Murray's law at every critical point","From discrete to continuum: Murray's law verified"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The continuum versions of Murray's law depend on the flux through the network being locally Lipschitz continuous so that integrating the pressure equation over arbitrary regions and applying the boundary-flux formula is legitimate; the paper constructs only less regular weak solutions and leaves this smoothness unproved.","fun_headline_variants_meta":{"raw":{"variants":["Murray's law holds across discrete and continuous models","Generalized Murray's law derived for biological networks","Branching networks follow Murray's law at every critical point","From discrete to continuum: Murray's law verified"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000603,"raw_usage":{"total_tokens":2789,"prompt_tokens":892,"completion_tokens":1897,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":508,"completion_tokens_details":{"reasoning_tokens":1836}},"tokens_in":508,"tokens_out":1897,"duration_ms":15879,"temperature":1.0,"reasoning_tokens":1836,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:20:04.238459+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete test: construct, for $\\gamma>1$ and a smooth source $S\\in L^2(\\Omega)$, the unique minimizer $p$ of (34) (or of (52)) on a rectangle, and compute the flux $q=-(rI+c)\\nabla p$; if $q$ has an interior discontinuity or cusp where the pressure gradient changes sign, the Green-formula step behind (31), (38), and (50) is not justified for that solution, so the continuum Murray law would not follow for the constructed solution.","supporting_citations":[{"cited_title":"ODE and PDE based modeling of biological transportation networks","cited_arxiv_id":"1805.08526","evidence_quote":"It derives the derivative identity (11) and the continuum-limit energy (16) that this paper restates and builds on."},{"cited_title":"Haskovec, L","cited_arxiv_id":null,"evidence_quote":"It supplies the rigorous continuum limit of the discrete model, giving the rectangular-mesh PDE model (19)–(23)."},{"cited_title":"Murray: The Physiological Principle of Minimum Work: II","cited_arxiv_id":null,"evidence_quote":"It is Murray's companion paper on oxygen exchange, the second origin of the classical 3/4-law."},{"cited_title":"Rogers: Laminar flow analysis","cited_arxiv_id":null,"evidence_quote":"It gives the laminar-flow relation that conductivity is proportional to radius to the fourth power, used to state the law in terms of $C^{3/4}$."},{"cited_title":"Evans and R.F","cited_arxiv_id":null,"evidence_quote":"It provides the De Giorgi–Federer Green formula that converts the integrated PDE into the boundary-flux Murray laws (31), (38), (50)."},{"cited_title":"Haskovec, P","cited_arxiv_id":null,"evidence_quote":"It establishes the linear-stability criterion meas(A)=meas(Ω) and constructs weak solutions of (52) by direct methods."},{"cited_title":"Haskovec, P","cited_arxiv_id":null,"evidence_quote":"It develops the γ=1 free-boundary reformulation and the cutoff construction for 1/2≤γ<1 stationary solutions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It derives the phenomenological continuum model (42)–(44) and its Darcy/Kirchhoff structure from discrete-network considerations."}],"review_version":1}