{"id":"6ee200f7-1326-4b6a-bd17-682abb58f6c4","arxiv_id":"2512.16053","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A strength-constrained crack-growth model reproduces echelon crack fragmentation that pure energy minimization cannot.","lead":"Using computer simulations, the paper shows that adding a material strength constraint to energy-based fracture theory makes a flat crack split into the staggered 'echelon' pattern seen in experiments, whereas pure energy minimization predicts an unbroken planar crack. The result offers a single explanation for a decades-old puzzle in materials science and suggests that classical energy-only crack-growth theory is incomplete.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Echelon prediction depends on an unverified phase-field regularization of the strength constraint; no proof that Eqs. (11)-(13) faithfully implement the constrained minimization (9) in mixed-mode tearing.","rationale":"The reader identified the phase-field implementation of the strength constraint as the weakest assumption, specifically the transfer of ce from calibration problems to the tearing geometry. This is indeed the load-bearing issue: the echelon pattern is a new prediction that depends entirely on ce, and the paper provides no rigorous link between the sharp constrained minimization (9) and the regularized equations (10)-(13). An additional arithmetic inconsistency exists in the claim that PDMS has twice the compressive-to-tensile strength ratio of graphite (using Tables 1-2 and Eq. (7): graphite σcs/σts≈2.8, PDMS≈2.1), which undermines an explanatory sub-claim but is secondary. The verdict CONDITIONAL remains appropriate because the paper shows a plausible mechanism and provides reproducible code, but the central assertion of general applicability must be backed by a demonstration that the phase-field model faithfully represents the constrained theory. The proposed test (re-deriving the regularization or comparing with a direct penalty approach) would settle whether the echelon prediction is robust or an artifact.","tokens_in":12197,"tokens_out":11504,"duration_ms":114582,"concrete_test":"Independently derive the phase-field equations for the constrained minimization (9) by adding a penalty term λ∫Ω0 χ{F(S)<0} (1−z)² dX to the regularized energy (3) and taking λ→∞; compute the resulting Euler-Lagrange equations and compare with (10)-(13). If the derived driving force does not match ce in Eq. (11), implement both versions for the graphite tearing geometry of Fig. 2(a) and compare crack patterns and orientation angles. A mismatch would demonstrate that the echelon pattern is an artifact of the chosen ce rather than a prediction of the strength-constrained minimization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the strength-constrained phase-field model (10)-(13). This model is not a direct minimization of a functional representing (9); rather, it is obtained by adding a stress-based driving force ce to the classical Euler-Lagrange equations, with coefficients β1ε, β2ε, δε calibrated in [30,40] to reproduce the strength surface under homogeneous stress states (uniaxial tension, hydrostatic stress). The echelon fragmentation in Fig. 1(d) and Fig. 3 is entirely generated by this term. There is no Γ-convergence proof, nor a numerical equivalence check, showing that (10)-(13) regularizes the constrained minimization (9) in the non-uniform, mixed-mode stress field of the tearing experiment. If the functional form of ce is not the unique (or even the correct) regularized representation of (9), the predicted echelon pattern may be a numerical artifact. The paper's own caveat that the phase-field model is \"not completely connected\" (§3 summary) underscores this gap. Absent an independent derivation or a sensitivity test against an alternative strength-constraint implementation, the claim that the strength surface explains echelon fragmentation is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Using two phase-field formulations, the paper revisits Knauss's mode-III tearing experiment on a notched plate. The classical variational Griffith model (3), minimizing elastic plus fracture energy, predicts continued planar crack growth for both graphite and PDMS. The authors then introduce a strength-constrained model (10)-(13) in which a Drucker-Prager strength surface restricts where crack growth can occur; this model spontaneously produces echelon daughter cracks without stochastic disorder or geometric imperfections. From the simulations, the paper identifies the shear-to-tensile strength ratio sigma_ss/sigma_ts and the thickness-to-material-length ratio H/l_ss^ch as the controlling parameters, and argues that energy-based and stress-based crack-path criteria are reconciled within one framework. The central conclusion is that Griffith energy minimization alone is fundamentally incomplete for predicting large-crack growth and that a strength constraint is necessary.","tokens_in":12545,"tokens_out":7983,"duration_ms":87252,"significance":"If correct, these results would be a substantive advance: they offer a deterministic explanation of a classical crack morphogenesis problem and tie together competing empirical criteria. The paper's strengths include a clean binary comparison between the classical and strength-constrained models on the same geometry, use of independently measurable material parameters (with the PDMS caveat below), and distribution of an open-source FEniCS implementation. However, the main prediction depends on an imported phase-field regularization of the strength constraint whose validity in the mixed-mode tearing state is not established; and the quantitative experimental support is a single historical data point. Thus the work is promising, but at present the force of the central claim exceeds its verification.","major_comments":[{"comment":"The echelon patterns in Figs. 1(d), 3, and 5 are generated entirely by the added driving force c_e in Eq. (10). Equations (11)-(13) are taken from Refs. [30,40], where the coefficients beta1^eps, beta2^eps, delta_eps were calibrated against homogeneous uniaxial and hydrostatic stress states (poker-chip, indentation). It is not demonstrated that this analytical form enforces the sharp constrained minimization (9) in the strongly inhomogeneous, mixed-mode stress field of the tearing problem. Since the classical model (3) gives planar growth in the same geometry, the strength-surface explanation is currently as sensitive to the regularization as to the physics. The paper's own caveat in the §3 summary ('the phase-field model is not completely connected') underscores the gap. I request an eps-convergence study at fixed h/eps, an independent implementation check, or a direct comparison with t","section":"§3, Eqs. (10)-(13)"},{"comment":"For PDMS, the paper states 'We have considered lower strength values than the actual values to make the simulations computationally less expensive.' This is a significant departure from the stated goal of using standard, experimentally measured material parameters. The PDMS results in Figs. 2(d) and 3(b), and the abstract claim of 'general applicability to any brittle material', are therefore not a test of the model for actual PDMS. In addition, the full-geometry comparison with Knauss's experiment (§3, Fig. 1(d)) does not state which material parameters were used; if those are the reduced PDMS values, the approximate 40° result cannot be viewed as a quantitative validation. Please rerun at least one PDMS case at actual strength or reframe the claim as qualitative.","section":"Appendix A, Table 2"},{"comment":"The text states that 'both the number and orientation of daughter cracks were found to depend primarily on the value of l_ss^ch,' but Fig. 5 plots only the number of daughter cracks. No orientation data are shown for the H/l_ss^ch study, so the second key parameter is not supported by the presented data. Also, the error bars are not defined and no statistical procedure is described. Please either add orientation results or temper the claim.","section":"§3 and Fig. 5"},{"comment":"The sole quantitative experimental benchmark is Knauss's reported inclination of about 45° compared with the predicted approximately 40°. There is no uncertainty on either value, no discussion of how the experimental angle was measured, and no account of the differences between Solithane (experiment) and the constitutive models used here. Given that the Abstract states the model is shown 'through comparison with classical experiments' to explain echelon formation, this comparison is too thin to carry the load. Additional validation against existing hydrogel/glass data, or at least a more careful treatment of this single data point, is needed.","section":"§3, experimental comparison"}],"minor_comments":[{"comment":"The sharp constrained minimization statement is informal: 'Gamma subset V_F(t)' conflates a 2D crack set with a 3D point set, and no topology/admissibility is specified. A precise formulation would help, especially because the paper argues that (10) is its regularization.","section":"Eq. (9)"},{"comment":"Please define how theta is measured from the phase-field contours. Only two crack extensions are used to support the orientation claim, so the measurement procedure matters.","section":"§3, angle extraction"},{"comment":"There is an inconsistency in whether the sweep varied sigma_ss or sigma_ss^2 by factors of 1-5. The text near Eq. (14) says 'sigma_ss^2' and a few lines later says 'sigma_ss was varied.' Please check consistency.","section":"§3 and Fig. 5"},{"comment":"The regularization length eps is described as a free parameter, but no convergence study in eps is reported. A sentence or figure showing that the reported morphologies are stable as eps decreases, with h=eps/4, would strengthen the paper.","section":"Appendix B"},{"comment":"Some figure panels are difficult to read in the current typesetting (e.g., Fig. 2(b) axes, Fig. 4(d) labels). A cleaner presentation would improve reproducibility and readability.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to attract broad interest in the fracture mechanics community. My main concern is whether the c_e term from Refs. [30,40] can be transplanted to mode-III tearing without verification; the reduced-strength PDMS simulation and the thin experimental comparison compound that risk. I would be willing to revise my recommendation upward if the authors add an eps-convergence/equivalence check and at least one PDMS run at actual material parameters."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead this one if you care about crack path prediction. The paper's central claim is that echelon fragmentation—a planar crack splitting into stepped daughter cracks under shear—is not explained by Griffith energy minimization alone, but by a strength-constrained minimization: crack growth is allowed only where a Drucker-Prager strength surface is violated. The authors show the classical variational phase-field model predicts planar growth, while their strength-constrained version spontaneously produces echelon cracks in graphite without any disorder. That is a clean, striking result, and the two-parameter rationalization (shear-to-tensile strength ratio sets orientation, H/lch sets fragmentation) is a useful synthesis. The collapse of crack counts on H/lch is compelling, and the use of independently measured graphite strengths is a real point in their favor. They also ship code.\n\nThe soft spots are mostly at the edges, but one sits near the load-bearing wall. The echelon patterns are produced by the additional driving force ce in Eq. (11), whose coefficients were calibrated in earlier papers on uniaxial tension, hydrostatic stress, and poker-chip tests. There is no Γ-convergence or numerical equivalence proof that this term faithfully regularizes the constrained minimization (9) in the non-uniform, mixed-mode tearing field. The authors themselves say the phase-field model is 'not completely connected' to the sharp formulation. Absent a sensitivity test against an alternative implementation—or a direct derivation—the possibility that the fragmentation is an artifact of the regularization cannot be dismissed. That is a real gap, and it needs to be addressed, not just asserted away.\n\nSecond, the empirical support is thinner than the abstract implies. The experimental comparison is a single historical angle from Knauss (1970), with no uncertainty. PDMS parameters are reduced to keep the computation cheap, which weakens the 'general applicability to any brittle material' claim. And while the paper criticizes disorder-based models, there is no quantitative head-to-head comparison showing the strength model does better. The 'variational model is fundamentally incomplete' conclusion is also stronger than what a single counterexample supports, especially since the variational model's failures could conceivably be fixed by better energy splits (though the authors argue convincingly those are unphysical for shear).\n\nStill, the core idea is novel and the qualitative predictions are right. This deserves a serious referee, but the referee should push for the regularization check and a more honest boundary on the claims.\n\nRecommended: send to peer review, but the authors should expect to do substantial additional work.","headline":"Interesting and likely important claim—strength matters for crack path—but the phase-field regularization that produces the echelon pattern is imported from prior work without a proof of fidelity in this geometry; the broad claims exceed the current evidence.","tokens_in":12932,"tokens_out":3205,"would_cite":true,"duration_ms":32923,"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":"A strength-constrained energy minimization—not energy competition alone—predicts when and where a planar crack fragments into echelon daughter cracks.","keywords":["echelon cracks","brittle fracture","phase-field model","strength surface","Griffith energy","crack path prediction","mixed-mode fracture","dimensionless parameters"],"falsifier":"A tearing experiment on a brittle material with known strength ratio, varying plate thickness H while holding the material length scale fixed: if echelon cracks appear for H/lss_ch below the critical threshold (or fail to appear above it), the predicted scaling is wrong. Alternatively, varying the shear-to-tensile strength ratio by changing the material should move the crack orientation angle between about 0° and 45° as the paper's Fig. 4(b) predicts.","tokens_in":12112,"feed_emoji":"💥","tokens_out":3517,"duration_ms":37780,"temperature":0.7,"pith_summary":"The paper argues that the classical energetic theory of fracture cannot explain a basic observation: a planar crack under out-of-plane shear often breaks into staggered, disconnected daughter cracks (echelon cracks). The authors show that adding a material strength surface to the variational energy-minimization framework reproduces this morphogenesis in simulations of both hard and soft brittle materials, without any assumed disorder or defects. They identify two dimensionless parameters—shear-to-tensile strength ratio and plate-thickness-to-material-length ratio—that set whether and how strongly a crack fragments. If correct, this resolves a long-standing dispute between energy-based and stress-based crack path criteria.","feed_headline":"Strength constraint explains why cracks fragment into echelon patterns","feed_subtitle":"Energy alone fails; adding material strength predicts the staggered crack patterns in brittle solids.","key_machinery":"The strength-constrained phase-field formulation: the classical variational fracture functional is minimized subject to the constraint that the crack set lies inside the region where the stress satisfies the Drucker–Prager strength surface F(S)=0. In practice, this is implemented by adding an explicit stress-based driving force ce to the phase-field evolution equation, with coefficients calibrated on uniform stress states. Crack growth then requires both that the strength surface be exceeded (nucleation) and that energy be minimized (propagation).","core_discovery":"Echelon crack formation is a nucleation phenomenon, not an instability of a smooth crack front: the parent crack spawns disconnected daughter cracks in regions where the local stress state first violates the material's strength surface. A phase-field model that minimizes elastic plus surface energy only among crack sets confined to these strength-exceeded regions reproduces the observed angled, stepped crack patterns; the purely energetic model instead predicts continued planar growth. The orientation angle of the daughter cracks increases with the shear-to-tensile strength ratio, and fragmentation occurs only when the plate thickness exceeds a characteristic strength-based length scale. Thu","pith_inferences":["The same strength-constraint mechanism may explain other 'spontaneous fragmentation' fracture patterns, such as crack front segmentation under mixed-mode loading in geological or composite materials.","If the two-parameter scaling is robust, experiments that vary only the sample thickness should show a sharp transition from planar to echelon cracking at a critical H/lss_ch, which could be tested with a series of hydrogel or polymer plates.","The authors' claim that the phase-field regularization length ε is freely tunable suggests the model may be used to resolve process-zone effects directly, perhaps linking echelon spacing to microstructure.","The framework's extension to anisotropic strength surfaces (e.g., Mohr-Coulomb or Hoek-Brown) is straightforward, which may yield orientation preferences in rocks with anisotropic cleavage."],"forward_implications":["Predicts echelon crack formation in both soft and hard brittle materials without invoking disorder, defects, or front perturbations.","Reconciles energy-based (maximum energy release rate) and stress-based (local mode I) path criteria: which one dominates is controlled by the shear-to-tensile strength ratio.","Identifies two non-dimensional parameters, σss/σts and H/lss_ch, that govern crack orientation and fragmentation, giving experimentally testable scaling laws.","Shows the purely variational Griffith model to be fundamentally incomplete for large-crack growth under mixed tension/shear.","Because the framework uses only measurable elastic, toughness, and strength properties, it offers a direct route to prediction in new materials."],"fun_headline_variants":["Crack fragmentation needs strength, not just energy","Echelon cracks arise from nucleation, not instability","Strength-limited nucleation predicts echelon crack angles","Energy alone can't predict crack paths; strength can"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The specific analytical driving force that encodes the strength constraint, calibrated on simple axisymmetric test geometries, is assumed to remain a faithful regularization of the constrained energy minimization in the mixed-mode tearing geometry; if it does not, the predicted echelon fragmentation could be a numerical artifact.","fun_headline_variants_meta":{"raw":{"variants":["Crack fragmentation needs strength, not just energy","Echelon cracks arise from nucleation, not instability","Strength-limited nucleation predicts echelon crack angles","Energy alone can't predict crack paths; strength can"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1329,"prompt_tokens":651,"completion_tokens":678,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":395,"completion_tokens_details":{"reasoning_tokens":617}},"tokens_in":395,"tokens_out":678,"duration_ms":6992,"temperature":1.0,"reasoning_tokens":617,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T15:38:44.164288+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A tearing experiment on a brittle material with known strength ratio, varying plate thickness H while holding the material length scale fixed: if echelon cracks appear for H/lss_ch below the critical threshold (or fail to appear above it), the predicted scaling is wrong. Alternatively, varying the shear-to-tensile strength ratio by changing the material should move the crack orientation angle between about 0° and 45° as the paper's Fig. 4(b) predicts.","supporting_citations":[],"review_version":1}