{"id":"6dcf329e-bf36-4f63-8687-9d3ffe76d491","arxiv_id":"2602.09737","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The Slenderness Ratio of a triangular beam lattice controls how much Weibull-distributed strut strength disorder affects crack path and toughness, yielding three failure regimes.","lead":"This paper shows that in 2D triangular beam lattices with randomly varying strut strengths, a geometric parameter called the Slenderness Ratio determines whether disorder is suppressed, bends the crack locally, or triggers widespread pre-failure. It offers a design rule for 'disorder-tolerant' architected materials and challenges the common idea that disorder toughens materials simply by making cracks longer.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Frozen stress-ratio and independent-event approximations break down at low n, exactly where the non-monotonic toughening peak and the diffuse-failure regime live; the central phase diagram is therefore not established in its most distinctive regime.","rationale":"The reader's weakest assumption is the same one I judge most load-bearing: the model computes P_s and P_d from stress ratios in the pristine lattice and then assumes these apply for every successive failure event. The numerical data validate this in the weak-to-moderate disorder range, which is a genuine and well-supported contribution. However, the central claim includes a full phase diagram and a decoupling statement, and both are least secure where the approximation fails—low n, where diffuse damage accumulates and alters the stress field. The non-monotonic toughening maximum is in this same low-n region, so the paper has not yet provided a mechanistic account of the most distinctive part of the toughening curve. This does not overturn the paper; it says the unconditional version of the claim is too strong. Since the reader already returned CONDITIONAL on essentially these grounds, my stress test leaves the verdict unchanged. The proposed test is deliberately targeted: re-evaluate κ after each failure event in the damaged lattice at a representative low-n point. That single experiment separates the two possible causes of the Fig. 5 breakdown—frozen stress ratios versus genuinely correlated/damage-dependent failure statistics—and would show whether a repair of the model is possible.","tokens_in":18217,"tokens_out":12017,"duration_ms":129531,"concrete_test":"Fix λ=10 and n=5 (a point where Fig. 5(a) deviates from Eq. 18) and re-run the open-source simulator, but after every crack-tip failure event recompute the six crack-tip stress ratios κ_ij from the instantaneous stress field in the damaged configuration rather than from the undamaged lattice. Compare the cumulative scattering count and P_d^0 against Eq. (18) using (a) frozen κ and (b) updated κ. If the updated-κ model restores agreement, the frozen-κ/independence assumption is the identified soft spot; if it does not, the independent-event approximation itself fails in the strong-disorder regime and the phase diagram needs an additional state variable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central predictions—the three-regime phase diagram and the disorder-induced toughening law G/G_u = 1+B/n—rest on the assumption that the crack-tip stress ratios κ_ij and the bulk-to-tip Weibull-stress ratio entering P_d are frozen at their undamaged-lattice values and that successive failure events are statistically independent (Eqs. 8, 14, 18–19; Sec. V). This assumption is not derived; it is an approximation. The paper's own Fig. 5(a) shows that the resulting P_s prediction agrees with simulation only for Weibull moduli n≳7–9 and then systematically under-predicts scattering at lower n, coinciding with the onset of diffuse damage that changes the stress field the crack sees. Because the non-monotonic peak in G/G_u (Fig. 6 inset) occurs at low n, the paper's claim that toughening is decoupled from the amount of damage or tortuosity is weakest precisely where the supporting theoretical model is acknowledged to break down. The authors flag this breakdown as a 'notable feature', but the phase-diagram boundaries in regime (iii) and the strong-disorder part of the central claim remain empirical observations rather than predictions of the framework.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a statistical framework for fracture of pre-cracked triangular beam lattices with Weibull-distributed beam strengths. The central idea is that the slenderness ratio λ = a/t controls the relative crack-tip stresses, and therefore controls whether disorder is expressed along the main crack ('scattering') or in the bulk ('diffuse failure'). Using the Weibull survival function, the authors derive closed-form probabilities P_s (Eqs. 18–19) and P_d^(0) (Eqs. 14–15), identify three regimes in Fig. 4, and validate against simulations in Fig. 5. They further report that normalized fracture energy G/G_u collapses to 1 + B/n for n ≳ 10 (Eq. 21), while absolute G is non-monotonic in n (Fig. 6b). They argue that disorder-induced toughening is not explained by crack tortuosity or damage count. The paper is clearly written, and the authors explicitly flag the regime where their model breaks down, but the abstract and final remarks present the three-regime picture as more broadly established than the supporting analysis warrants.","tokens_in":18604,"tokens_out":6357,"duration_ms":66970,"significance":"If correct, the paper establishes geometry as a design axis for controlling disorder expression and provides a parameter-free prediction of crack-path statistics at moderate disorder. Its concrete strengths include the overlaid theory lines in Fig. 5, which contain no fitted parameters, and the public availability of code and data. The paper also cleanly separates the geometric contribution G_u from the statistical contribution G/G_u. However, the framework's predictive range is limited: the assumption of independent events with frozen crack-tip stress ratios is explicitly shown to fail for Weibull moduli n ≲ 7–9, exactly where diffuse failure and the non-monotonic toughening maximum occur. The phase diagram and the strong-disorder portion of the central claim therefore rest on simulation data rather than on the statistical model, and the manuscript overstates the degree to which the theory explains them.","major_comments":[{"comment":"The scattering probability is constructed from crack-tip stress ratios κ_i(λ) computed in the undamaged lattice and assumes statistically independent failure events via Eq. (3). The paper's own Fig. 5(a) shows this prediction systematically under-predicts excess broken bonds for Weibull moduli n ≲ 7–9, the range where diffuse damage appears and the stress field seen by the crack changes. Because regime (iii) in Fig. 4 and the peak of G in Fig. 6(b) lie in this range, the phase-diagram boundaries and the strong-disorder portion of the central claim are not predictions of the framework; they are empirical observations. The authors note the breakdown, but the three-regime picture is presented as established in the abstract and final remarks. Please either restrict the theoretical claims to the validated regime or extend the model (e.g., by including the effect of pre-existing diffuse damage","section":"Sec. V, Eqs. (18)–(19); Fig. 5(a)"},{"comment":"The master curve G/G_u = 1 + B/n is obtained by fitting B to data normalized by the simulated G_u; B is not predicted from the model. The explanation via crack-arrest scaling with the standard deviation of Weibull strengths (∼1/n) is an analogy, not a derivation from the lattice model. Moreover, the claim that toughening is decoupled from the amount of damage/tortuosity is not demonstrated by a direct comparison: the paper shows that G/G_u collapses while P_s varies with λ, but it does not plot G/G_u against the measured excess-damage fraction or effective crack path. To make this central interpretative claim quantitative, include such a plot or correlation test, or soften the wording to say the data are not directly correlated rather than 'cannot be connected' to damage.","section":"Sec. VI B, Eq. (21); Fig. 6(b)"},{"comment":"The comparison in Fig. 5(b) is less direct than stated. P_d^(0) is exact only for the first failure event (as noted in the text), yet the simulation values are y-intercepts of fits of Eq. (17), which contains free parameters β and γ. The theoretical line for P_d^(0) also depends on a representative bulk stress σ̄ and the effective volume V/V0, whose values are not specified in the main text. Please provide either a direct measurement of the first-failure probability from simulations or a clear statement of how σ̄ and V/V0 are obtained, so the reader can assess whether the agreement in Fig. 5(b) is a parameter-free test.","section":"Sec. IV C, Eqs. (14)–(15), (17); Fig. 5(b)"}],"minor_comments":[{"comment":"Typographical issues: 'theSlenderness Ratio' in the abstract, and Sec. II's 'mode axial and bending failure stresses' should likely be 'moduli' or 'ultimate axial and bending failure stresses'. Please define λ at first use.","section":"Abstract and Sec. II"},{"comment":"The hazard-rate construction uses W(σ_i; σ_0, n) and S(σ_j; σ_0, n); please state explicitly that S is the survival probability of element j and that the two failure thresholds are assumed independent and identically distributed.","section":"Eq. (5)"},{"comment":"The nonlinear correction in Appendix D is useful, but it is not incorporated into the main predictions. State explicitly that the main-text results assume constant κ_ij and that the nonlinear correction is a separate, non-validated extension.","section":"Appendix D"},{"comment":"The hierarchical agglomerative clustering cutoff of 1.1a is central to extracting the main-crack cluster, but no sensitivity analysis is given. A brief robustness check (varying the cutoff) would increase confidence in the reported N_f and excess-damage ratios.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and, with revision, could be a useful contribution. The main risk is overclaiming from a model that is only self-consistent at moderate disorder; the authors should be pushed to either extend the model or explicitly delimit the regime of validity. No concerns about citation practices or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis one is worth a look. The genuinely new piece: geometry—specifically the slenderness ratio λ—is treated as a control variable for whether Weibull disorder in a triangular lattice expresses itself in the crack path, and the paper derives closed-form probabilities for scattering and diffuse failure from crack-tip stress ratios and Weibull order statistics. The theory lines have no free parameters and match simulations for n ≳ 7–9 across all tested λ. That is a real result, and the code/data are public.\n\nWhat the paper does well. The three-element discretization is a sensible extension of the random beam model and is experimentally anchored to prior PMMA work. The derivation of Ps (Eq. 18) and Pd(0) (Eq. 14) is clean, and the renormalization via Weibull stress is a nice touch. The reproduction of disorder-induced toughening and the collapse of G/G_u onto 1 + B/n for n ≳ 10 is convincing evidence for the crack-arrest interpretation over simple tortuosity arguments. The authors are also unusually honest: they flag the breakdown at low n as “a notable feature” and speculate about the cause rather than sweeping it under the rug.\n\nSoft spots, in proportion. The central weakness is exactly where the reader placed it: the frozen-stress-ratio and independent-event approximation is not derived, and the paper’s own Fig. 5 shows it under-predicts scattering for n ≲ 7–9, coinciding with where diffuse damage sets in. So regime (iii) and the strong-disorder part of the phase diagram are empirical observations, not predictions of the framework. The non-monotonic toughening peak lives in that same regime and is unexplained, and the 1/n law is fitted in the collapsed regime only (with B a free parameter). Also, the claim that crack-path morphology is system-size-independent is asserted but never tested by varying system size for disordered lattices; only Gu convergence in uniform lattices is tested. Those are addressable rather than fatal, but they should be acknowledged more carefully.\n\nWho it’s for: anyone working on fracture in architected lattices or on disorder-controlled damage. The weak-to-moderate disorder regime is established, and the separation of morphology from energetics is worth arguing with.\n\nI’d send it to a serious referee. My own verdict would be conditional: keep the strong-disorder claims descriptive, test the size-independence claim, and the paper will be solid.","headline":"A clean, largely parameter-free theory for how slenderness ratio controls crack-path disorder in the weak-to-moderate disorder regime; the strong-disorder part of the phase diagram rests on simulation, but the core result stands.","tokens_in":19035,"tokens_out":2700,"would_cite":true,"duration_ms":28391,"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":"The slenderness ratio of a triangular beam lattice controls how Weibull-strength disorder is expressed, yielding three fracture regimes in one phase diagram.","keywords":["beam lattice","Weibull disorder","slenderness ratio","fracture toughness","crack path","diffuse failure","damage evolution","triangular lattice"],"falsifier":"The independent-event assumption predicts that the scattering probability at a given λ and n is unaffected by pre-existing diffuse damage; therefore, a lattice in which a few random bulk bonds are removed before loading should show the same P_s as an undamaged one, and the G/G_u master curve should not shift. If either changes, the fixed-κ picture is wrong.","tokens_in":18131,"feed_emoji":"📐","tokens_out":10770,"duration_ms":94584,"temperature":0.7,"pith_summary":"This paper argues that in a pre-cracked triangular beam lattice with Weibull-distributed failure stresses, the slenderness ratio (unit-cell size divided by beam thickness) is the geometric control parameter that determines how disorder expresses itself: high slenderness suppresses disorder, intermediate values cause local crack-tip scattering that roughens the crack, and low values allow initially diffuse bulk failure. The authors derive closed-form probabilities for these outcomes from pairwise Weibull failure statistics together with the crack-tip stress ratios, then test them against quasistatic lattice simulations. They also reproduce the disorder-induced increase in apparent fracture energy and find that, for Weibull moduli n at or above about 10, it collapses to G/G_u = 1 + B/n, independent of slenderness, while the absolute fracture energy depends non-monotonically on disorder. Their key negative result is that this toughening does not scale with the number of excess broken bonds or crack tortuosity, challenging the common explanation of disorder toughening by diffuse damage. The paper is explicit that it does not yet model the non-monotonic peak at strong disorder, and its own simulations show the independent-event approximation degrades for n below about 7 to 9, where diffuse damage is significant.","feed_headline":"Slenderness ratio controls crack path in disordered lattices","feed_subtitle":"A single geometric parameter decides whether disorder is muted, local, or diffuse — and toughening obeys a 1/n law.","key_machinery":"The central object is the Slenderness Ratio λ ≡ a/t, the ratio of the triangular unit-cell size to the beam thickness, which sets the relative axial and bending stress content of the six crack-tip beams and therefore the stress ratios κ_ij among them. The central identity is the pairwise Weibull failure probability P(s_i > s_j) = κ_ij^n / (1 + κ_ij^n), obtained by absorbing κ_ij into a rescaling of the reference failure stress; this identity makes the statistics functions of λ and n only. The paper extends this pairwise rule to groups using the Weibull stress, the L^n norm of the element stresses, which lets it compare the whole crack-tip set against the bulk and thus predict diffuse failure","core_discovery":"The central claim is that the slenderness ratio λ = a/t, the ratio of the triangular unit-cell size to the in-plane beam thickness, acts as a geometric control parameter that selects among three fracture regimes in a disordered beam lattice: disorder suppressed, local crack-tip scattering, and initially diffuse failure. The mechanism is the reordering of the six crack-tip stress ratios with λ, because the Weibull survival function is scale-invariant, the probability that element i fails before element j is exactly the ratio κ_ij^n/(1 + κ_ij^n), where κ_ij is the fixed stress ratio. Summing these over the crack-tip elements and comparing the crack tip to the bulk through the Weibull stress yi","pith_inferences":["A natural extension is to pattern slenderness spatially across a lattice: gradients in λ should steer the crack through different regimes along a single specimen, effectively writing the crack path by design rather than by disorder realization.","The identity P = κ^n/(1 + κ^n) suggests an inverse-problem use: measuring the scattering probability as a function of n at fixed λ would recover the effective crack-tip stress ratio κ_s, making disorder a quantitative probe of crack-tip micromechanics.","The 1/n collapse hints that weak Weibull disorder acts like a spatially fluctuating local toughness field; a testable consequence is that any weak disorder distribution whose failure thresholds have variance proportional to 1/n should produce the same master curve.","The low-n breakdown could itself be modeled by letting κ_ij depend on accumulated diffuse damage, effectively renormalizing the crack-tip stress ratios as the bulk weakens; such an extension would move the phase diagram's strong-disorder boundary from an assumption into a prediction."],"forward_implications":["Choosing the slenderness ratio at fixed disorder strength lets a designer select the fracture regime: cracks can be forced straight, made locally tortuous, or driven to fail diffusely, independent of the statistical spread of failure strengths.","In the weak-to-moderate disorder regime (n ≳ 10), the disorder-induced toughening relative to the uniform lattice collapses onto a single 1/n curve, so apparent fracture energy enhancement can be predicted without fitting slenderness-dependent parameters.","The amount of excess damage along the crack path is not a proxy for fracture energy: lattices with very different degrees of tortuosity can show the same normalized toughening.","Below a Weibull modulus of about 7 to 9, diffuse damage accumulates before the main crack advances, the independent-event approximation breaks down, and geometric control of the crack path weakens.","Because the uniform reference fracture energy G_u depends on slenderness while the relative enhancement G/G_u does not, geometry and statistical disorder contribute separately to the fracture energy in the collapsed regime."],"fun_headline_variants":["Geometry dictates how disorder breaks lattices","Slenderness ratio: master switch for fracture regimes","One geometric knob tunes disorder in crack paths","Slenderness ratio selects fracture regime","Geometry decides: disorder muted, local, or diffuse"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the crack-tip stress ratios κ_ij and the bulk-to-tip stress ratio that enters the diffuse-failure probability stay fixed at their undamaged-lattice values, so each failure event sees the same stress hierarchy and successive events are statistically independent; the paper's own simulations show this premise degrades once the Weibull modulus falls below about 7 to 9, where diffuse damage alters the stress field.","fun_headline_variants_meta":{"raw":{"variants":["Geometry dictates how disorder breaks lattices","Slenderness ratio: master switch for fracture regimes","One geometric knob tunes disorder in crack paths","Slenderness ratio selects fracture regime","Geometry decides: disorder muted, local, or diffuse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000662,"raw_usage":{"total_tokens":2859,"prompt_tokens":741,"completion_tokens":2118,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":2050}},"tokens_in":485,"tokens_out":2118,"duration_ms":14159,"temperature":1.0,"reasoning_tokens":2050,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T02:41:44.982172+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The independent-event assumption predicts that the scattering probability at a given λ and n is unaffected by pre-existing diffuse damage; therefore, a lattice in which a few random bulk bonds are removed before loading should show the same P_s as an undamaged one, and the G/G_u master curve should not shift. If either changes, the fixed-κ picture is wrong.","supporting_citations":[],"review_version":1}