{"id":"dd24b4cb-37c5-4a7f-990e-0a40f23ad920","arxiv_id":"2507.12247","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Configurational force computations on idealized crack geometries show that facet tilting, local binormal shear, and facet spacing govern echelon crack growth and coalescence in soft materials.","lead":"This paper uses a computational method called configurational force analysis to study how echelon cracks, the tilted zigzag fracture patterns seen in gels, rocks, and glass, form in soft materials. It finds that facet spacing and local shear control whether crack segments grow, merge, or shield each other, offering a new lens on mixed-mode fracture.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central mechanistic claim that local binormal shear F_BN controls the deviation of the force angle β rests on a single hand-chosen facet tip profile, with no propagation simulation or variation of the tip shape; a different realistic profile could change the trend and undermine the explanation.","rationale":"The reader's weakest_assumption identifies exactly the issue I find most load-bearing: the idealized facet tip profile and the unvalidated causal link between F_BN and the deviation of β. My stress-test agrees with that assessment and with the conditional verdict. The concern does not require a change to the reader's verdict—it reinforces it. I considered whether the circularity of using echelon geometries to explain echelon cracks is an even more fundamental objection, but the paper's focus is on interactions and propagation tendencies of already-segmented cracks, which is a legitimate but more modest claim. The title's 'explain' is stronger than the evidence, but the underlying computations are genuine and the qualitative trends are internally consistent. The conditional verdict correctly captures that the paper is a plausible computational contribution whose central mechanistic conclusions need further validation, most directly by testing robustness to the tip-shape parameterization and by comparing with measured propagation directions. No independent evidence, such as machine-checked proofs or reproducible code, is provided, so the burden falls on these additional checks. I therefore see no reason to move away from the reader's conditional assessment.","tokens_in":20705,"tokens_out":7572,"duration_ms":93095,"concrete_test":"Repeat the configurational-force computation of Sections 4.1–4.3 using a different facet tip profile, e.g., an elliptical profile or the circular-crack approximation used in Ortellado et al. [21], keeping all other mesh, material, and loading parameters fixed. If the F_BN versus φ trend, the β versus φ curves, or the sign reversal of β for φ ≤ 20° under torsion changes qualitatively, then the conclusion that F_BN causes the deviation from LEFM scaling (Section 4.3) is not robust to the choice of idealized tip geometry. If the trends persist across profiles, the concern is mitigated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that configurational forces explain echelon cracks in soft materials, including the direction of propagation and the role of local shear. The evidence for this claim has a critical gap: the facet geometry is an input, not a prediction. Section 2.2 prescribes a specific tip profile f(y) = ±sqrt(1-2y+2y^3-y^4)/2, and the tilt angle φ is varied as a parameter, but the model never predicts why a planar crack should segment into tilted facets. The observed tilt angle is put into the model, making the title-level 'explain' claim circular for the phenomenon of facet formation itself. More specifically, Section 4.3 attributes the deviation of the force angle β from LEFM scaling to the binormal shear component F_BN. This conclusion is supported only by showing that F_BN and β both vary with φ and loading; no independent variation of the tip shape is performed. If a different realistic facet profile, such as the elliptical or circular-crack approximations used in the prior literature [21, 11], changed the F_BN trend or the sign reversal of β for φ ≤ 20° under torsion, the proposed mechanistic explanation would be undermined. The absence of a propagation simulation also means that the computed configurational force direction is not shown to correspond to actual crack growth directions; the paper itself notes that crack propagation is not directly simulated. Together, these omissions make the central explanatory claim—especially the F_BN-driven deviation mechanism—insufficiently supported for the strong conclusions drawn in Sections 4.3 and 7.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies mixed-mode I+III fracture in a soft hydrogel by computing configurational forces with the Configurational Force Method (CFM) implemented as a finite-element post-processor. The authors construct idealized three-dimensional crack geometries—a planar parent crack with tilted finger-like facets, and a variant with a coalesced type B crack—informed by X-ray tomography observations of echelon cracks, and subject these geometries to simulated uniaxial tension followed by torsion. They report how the magnitude and orientation of the computed configurational force per unit crack-front length vary with facet tilt angle, facet spacing, and loading, and they interpret the results as revealing the roles of facet tilting, local binormal shear, facet spacing, and coalescence in determining propagation direction and driving force. The paper explicitly notes that crack propagation is not directly simulated.","tokens_in":20948,"tokens_out":5295,"duration_ms":66245,"significance":"If the central interpretation were fully supported, the work would offer a useful nonlinear, finite-strain alternative to LEFM-based analyses of echelon cracks, with the CFM being generally applicable to arbitrary material behavior and complex crack-front geometry. The paper has clear strengths: the FE and CFM procedures are standard and carefully described, the neo-Hookean material parameters are calibrated against experiments, the mesh-generation pipeline is documented in detail, parameter sweeps over tilt angle and spacing are systematic, and the authors include appendices addressing facet-spacing and finite-size effects on the direction angle β. However, the title-level claim that configurational forces \"explain echelon cracks\" is not yet supported because the segmented facet geometry is prescribed as an input rather than predicted, and because the mechanistic attribution of the deviations from LEFM scaling to a single local shear component rests on a single hand-chosen tip profile. These are load-bearing gaps, but they are addressable by reframing the claims and by adding targeted numerical tests.","major_comments":[{"comment":"The idealized crack geometry is an input, not a prediction: the tilt angle φ, the angular spacing Λ, and the tip profile f(y)=±√(1−2y+2y^3−y^4)/2 are all taken from experimental observations of echelon cracks (Figs. 2a–b, Eq. B.1), and φ is then varied as a parameter. The simulations therefore characterize the configurational forces on an already-segmented front; they do not explain why a planar front segments into tilted facets or why the observed tilt and spacing are selected. Since the paper itself states in §7 that \"crack propagation is not directly simulated,\" and no comparison with measured propagation directions is provided, the title-level claim that the results \"explain echelon cracks\" overreaches the evidence. The authors should either reframe the claims as a characterization of propagation tendencies of prescribed geometries or add a test that predicts the observed morphology, for example by scanning φ and Λ and showing that the energy-release-rate-maximizing configuration matches the experimental values.","section":"§2.2, Appendix B, and title/abstract"},{"comment":"The central mechanistic claim that local binormal shear F_BN causes the deviation of the force angle β from the LEFM-based scaling is supported only by a single prescribed tip profile and by the observation that the maximum value of F_BN and the angle β both vary with φ and loading. No independent variation of the tip shape is performed: if a different realistic facet profile, such as the circular or elliptical crack approximations used in earlier work [21, 11], changed the F_BN trend or the sign reversal of β for φ≤20°, the explanation would be undermined. In addition, F_BN is a raw component of the deformation gradient, which includes rigid rotation, not a shear strain measure; the claim that it quantifies \"shear deformation\" should be verified using a frame-invariant measure such as the Green-Lagrange strain component E_BN or a polar-decomposition-based shear angle. I recommend adding a sensitivity study over at least two alternative tip profiles and reporting the corresponding F_BN/E_BN and β trends.","section":"§4.3, Fig. 5"},{"comment":"The attribution of the discrepancy between simulation and LEFM scaling to F_BN is underdetermined. The scaling model neglects non-singular T-stress terms, facet–facet interactions, finite specimen size, and nonlinear material behavior, and Appendix D itself shows a finite-size effect on β for torsion angles α>36°. The observation that β<φ under pure tension and β<0 for φ≤20° under torsion could therefore be caused by any of these omitted effects or by combinations of them. A direct causal test is needed: for example, isolating a single facet in an otherwise identical geometry, or comparing the simulation against an LEFM calculation that includes the T-stress, would help establish that the F_BN contribution—rather than facet interactions or non-singular fields—is the operative mechanism.","section":"§4.2, Eqs. (12)–(14), and Appendix D"}],"minor_comments":[{"comment":"The caption sentence \"Angle β ... calculated along the loading of the cylinder in (f)\" is unclear; the definition of β and the meaning of panel (f) should be stated explicitly.","section":"Fig. 4 caption"},{"comment":"The construction of the facet profile from f(y)=±√(1−2y+2y^3−y^4)/2 on y∈[−1,0.8] plus the vertical line x=0.8 should be explained more fully: it is not obvious how the coordinates are normalized, why the endpoint is 0.8, and how the resulting width-to-length ratio compares with the experimental facet measurements quoted in §2.1.","section":"Appendix B, Eq. (B.1)"},{"comment":"The rheology specimen height is given as 3500 nm (3.5 μm), which appears to be a typo; the intended value is likely 3.5 mm. Please correct the unit.","section":"Appendix A"},{"comment":"The legend label \"planar\" is used both for the facet-free penny-shaped reference crack and for the planar segments of the parent crack; this dual use is confusing and should be disambiguated.","section":"§5.1"},{"comment":"The paragraph beginning \"To assess potential finite-size effects, we repeated the analysis...\" is duplicated verbatim; one copy should be removed.","section":"Appendix D"},{"comment":"The statement that data and software \"will be available upon reasonable request\" is vague for a computational study; depositing the meshes, scripts, and post-processing routines in a public repository would materially improve reproducibility.","section":"Data and Software Availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of cond-mat.soft and computational fracture mechanics, and the CFM-based analysis is novel in this context. My main concern is that the manuscript's central explanatory claim is broader than what the computations actually establish: the segmented geometry is prescribed from experiments and propagation is not simulated, so the results characterize driving forces on assumed morphologies rather than explaining their formation. The F_BN mechanism in particular needs a sensitivity study over tip profiles and a proper shear-strain measure. These issues are fixable by reframing and adding targeted numerical tests, so I recommend major revision rather than rejection. I would also encourage the editor to ask for the missing convergence data for the configurational-force integration along the crack front, since the quantitative comparisons in Figs. 4 and 6–9 depend on it."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis is a solid computational study, but the title oversells it. The paper computes configurational forces on idealized 3D echelon crack geometries in a neo-Hookean gel, and that alone is new relative to prior phase-field and LEFM work. It shows clear trends: facet tilt controls driving force magnitude, spacing modulates amplification versus shielding, and coalescence redistributes force asymmetrically. The methods are standard and transparent, the material parameters come from their own tensile and rheology calibration, and the molecular-dynamics mesh relaxation for the coalescence geometry is a clever technical fix.\n\nThe soft spots are real. The facet geometry is prescribed from experiment, not predicted, so the title's 'explain' is too strong for facet formation itself. More importantly, the central mechanistic claim—that binormal shear F_BN causes the deviation of the force angle β from LEFM scaling—rests on a single hand-chosen tip profile and no propagation simulation. The authors themselves note propagation is not simulated, and they do check spacing and finite-size effects, which helps. But correlation between F_BN and β across tilt angles is not causation; a different tip shape could change that trend. This is a gap, not a fatal flaw.\n\nThe circularity concern is fair but moderate. The paper isn't claiming to predict that segmentation happens; it's asking what configurational forces act on observed geometries. That is a legitimate question. The weaker part is presenting the result as an 'explanation' of echelon cracks when the geometry is an input.\n\nThe audience is researchers in fracture mechanics, especially those studying mixed-mode I+III. The paper deserves a serious referee. For a public version, I'd ask for either a propagation simulation or at least a second facet tip profile, and a softened title and abstract. The qualitative results on amplification and shielding are likely useful even if the F_BN mechanism turns out not to be the whole story.","headline":"Solid computational study of configurational forces on idealized echelon cracks, but the title overclaims: geometry is input and the F_BN mechanism rests on a single tip profile.","tokens_in":21578,"tokens_out":2875,"would_cite":true,"duration_ms":33461,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74R10","74B20","74S05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that configurational forces computed on idealized echelon crack geometries explain both the magnitude and direction of crack growth in soft materials under mixed-mode I+III loading, via facet tilting, local binormal…","keywords":["echelon cracks","mixed-mode I+III fracture","configurational forces","soft fracture","hydrogel","facet coalescence","finite element method","crack segmentation"],"falsifier":"Compute the configurational forces on the same echelon geometry while varying the facet tip profile (e.g., sharp versus blunted tips, different $f(y)$) and check whether the evolution of $F_{BN}$ and the angle $\\beta$ tracks the tip-shape changes; if $\\beta$'s deviation from LEFM scaling persists without a corresponding $F_{BN}$ trend, the claimed mechanism is false. Alternatively, measure the actual propagation directions of echelon facets in a hydrogel under mixed-mode I+III loading using time-resolved X-ray tomography and compare them with the configurational-force directions predicted here, including the predicted sign reversal of $\\beta$ for small facet tilting angles.","tokens_in":20423,"feed_emoji":"💧","tokens_out":7488,"duration_ms":81051,"temperature":0.7,"pith_summary":"The paper aims to explain why echelon cracks—the segmented, tilted facet patterns that form when a tensile crack is loaded with out-of-plane shear—grow the way they do in soft materials, using configurational forces rather than linear elastic fracture mechanics. It constructs idealized finite-element models of echelon crack geometries informed by X-ray tomography of hydrogels and computes the configurational force (the energetic driving force for crack advance) along the crack front before propagation. The central claim is that this single field captures both the magnitude and direction of facet propagation: facet tilting controls the driving force, local binormal shear $F_{BN}$ causes the propagation direction to deviate from LEFM-based scaling, facet spacing modulates whether neighboring facets amplify or shield each other, and facet coalescence redistributes forces asymmetrically. If correct, this gives a material-model-independent way to predict echelon crack morphology in soft, quasi-brittle solids, where LEFM assumptions fail.","feed_headline":"Configurational forces explain echelon crack growth in soft solids","feed_subtitle":"Facet tilting and local shear set the crack direction; spacing decides whether facets amplify or shield growth.","key_machinery":"The Configurational Force Method: the Eshelby stress tensor $\\Sigma = \\Psi I - F^T \\partial \\Psi/\\partial F$ is evaluated from a finite-element solution of the finite-strain neo-Hookean boundary value problem, and its divergence is integrated against nodal basis functions to yield nodal configurational forces $F^{CNF}_A$; summing these forces along the crack front and dividing by arc length gives the energetic driving force per unit length whose magnitude and direction indicate propagation tendency and preferred direction. The idealized echelon geometry (a parent crack with tilted finger-like facets of profile $f(y)=\\pm\\sqrt{1-2y+2y^3-y^4}/2$, symmetric angular spacing $\\Lambda$, and tilting angle $\\varphi$) and a local coordinate frame $(\\hat{T},\\hat{N},\\hat{B})$ at the crack front are the other central objects; the off-diagonal deformation-gradient component $F_{BN}$ is the local shear measure invoked to explain deviations from LEFM scaling. A molecular-dynamics-based mesh relaxation smooths the coalescence junction between facets.","core_discovery":"The paper's central discovery is that the configurational force vector at the crack tip, computed as a nodal force from the Eshelby stress on an idealized echelon crack geometry, characterizes both the tendency and the preferred direction of propagation for maximal energy release rate. Under pure tension, the force magnitude is higher for smaller facet tilting angles $\\varphi$, and the binormal component dominates while the tangent component is negligible, meaning facets propagate while keeping their tilt. A non-zero normal component produces an angle $\\beta$ between the force and the facet plane that increases under tension and decreases under torsion, even reversing sign for $\\varphi \\leq 20^\\circ$ at large torsion; the paper attributes this deviation from the LEFM scaling to local shear deformation $F_{BN}$ in the binormal direction. The paper further shows that facet spacing modulates the driving force on both facets and planar segments, with a critical spacing below which shielding dominates, and that the growth of a type B crack connecting two facets shields nearby facet regions while amplifying forces farther away.","pith_inferences":["We infer that the facet tip profile is not a passive detail: because $F_{BN}$ is invoked as the mechanism, changing the tip shape should measurably change $\\beta$ and hence the crack direction, giving a testable knob for steering cracks in soft adhesives and gels.","We infer that the amplification-shielding crossover with facet spacing implies a design principle: patterning initial flaws or tuning the mode III fraction could set facet spacing to maximize toughness, since closer spacing first amplifies and then shields the driving force.","We infer that applying the same nodal configurational force calculation directly to X-ray tomography meshes of real cracks, rather than idealized geometries, would test the $F_{BN}$ mechanism without the idealization; the paper itself notes this as a future direction.","We infer that if the predicted sign reversal of $\\beta$ for small facet tilting angles is confirmed experimentally, local-symmetry-type criteria for mixed-mode I+III fracture in soft materials would need a correction term accounting for local binormal shear."],"forward_implications":["If the central claim holds, planar segments between facets carry a configurational force up to about four times that of the facets, indicating that these planar regions reach the propagation threshold first and so type B cracks form before further facet growth.","Facet tilting angle directly sets the driving force magnitude: lower $\\varphi$ aligns the facet with the Mode I opening stress, maximizing the configurational force, while the dominant binormal component means facets advance without changing their tilt.","The non-zero normal component $F_{CNF,\\hat{N}}$ implies that echelon facets can undergo slight rotations in their propagation direction even without in-plane shear, offering a mechanism for the qualitative facet waviness seen in experiments.","Facet spacing acts as a control variable: intermediate spacings amplify the driving force through elastic interactions, but below a critical spacing shielding dominates and the force drops, so spacing governs the stability of segmented crack growth.","Facet coalescence via a type B crack breaks the symmetry of the force distribution, shielding facet regions near the junction while amplifying more distant regions, which directs subsequent crack growth toward the unshielded parts of the coalescing front."],"supporting_citations":[{"why":"Provides the experimental hydrogel echelon crack morphologies and the principle-of-local-symmetry context that the idealized geometry is built to replicate.","marker":"[21]"},{"why":"Supplies the X-ray tomography reconstructions of mixed-mode I+III cracks whose finger-like facet shapes and spacings inform the simplified geometrical model.","marker":"[53]"},{"why":"Establishes the Configurational Force Method as an accurate estimator of the J-integral in soft solids, the methodological foundation for the paper's force computations.","marker":"[48]"},{"why":"Provides the phase-field and coupled-criterion modeling of echelon cracks that motivates the configurational-mechanics approach and supplies the $K_{III}\\sim\\sin 2\\varphi$ scaling used in the LEFM comparison.","marker":"[30]"},{"why":"Documents the experimental growth of type A and type B cracks and facet interactions that the paper's spacing and coalescence results are compared against.","marker":"[11]"},{"why":"Justifies including spurious nodal configurational forces near the crack tip to improve the approximation of the total driving force.","marker":"[60]"},{"why":"Provides the stress intensity factor relations for inclined cracks used to build the LEFM scaling estimate of the propagation angle $\\beta$.","marker":"[62]"},{"why":"Shows how spatial arrangement of nearby cracks amplifies or shields stress intensity factors, supporting the paper's interpretation of facet spacing effects.","marker":"[64]"}],"fun_headline_variants":["Configurational forces map echelon crack growth in soft solids","Shear and tilt steer echelon cracks via configurational forces","Configurational forces reveal facet spacing rules for echelon cracks","Echelon crack morphology emerges from configurational forces","Configurational forces dictate echelon crack paths in soft solids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The idealized facet tip profile $f(y)=\\pm\\sqrt{1-2y+2y^3-y^4}/2$ and the symmetric equal-spacing arrangement are assumed to represent real echelon cracks, and the conclusion that local binormal shear $F_{BN}$ causes the deviation of $\\beta$ from LEFM scaling is drawn without varying the tip shape or comparing against measured propagation directions.","fun_headline_variants_meta":{"raw":{"variants":["Configurational forces map echelon crack growth in soft solids","Shear and tilt steer echelon cracks via configurational forces","Configurational forces reveal facet spacing rules for echelon cracks","Echelon crack morphology emerges from configurational forces","Configurational forces dictate echelon crack paths in soft solids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000967,"raw_usage":{"total_tokens":4126,"prompt_tokens":967,"completion_tokens":3159,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":3069}},"tokens_in":583,"tokens_out":3159,"duration_ms":28815,"temperature":1.0,"reasoning_tokens":3069,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:50:52.173719+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the configurational forces on the same echelon geometry while varying the facet tip profile (e.g., sharp versus blunted tips, different $f(y)$) and check whether the evolution of $F_{BN}$ and the angle $\\beta$ tracks the tip-shape changes; if $\\beta$'s deviation from LEFM scaling persists without a corresponding $F_{BN}$ trend, the claimed mechanism is false. Alternatively, measure the actual propagation directions of echelon facets in a hydrogel under mixed-mode I+III loading using time-resolved X-ray tomography and compare them with the configurational-force directions predicted here, including the predicted sign reversal of $\\beta$ for small facet tilting angles.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental hydrogel echelon crack morphologies and the principle-of-local-symmetry context that the idealized geometry is built to replicate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the X-ray tomography reconstructions of mixed-mode I+III cracks whose finger-like facet shapes and spacings inform the simplified geometrical model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Configurational Force Method as an accurate estimator of the J-integral in soft solids, the methodological foundation for the paper's force computations."},{"cited_title":"& Lazarus, V","cited_arxiv_id":null,"evidence_quote":"Provides the phase-field and coupled-criterion modeling of echelon cracks that motivates the configurational-mechanics approach and supplies the $K_{III}\\sim\\sin 2\\varphi$ scaling used in the LEFM comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the experimental growth of type A and type B cracks and facet interactions that the paper's spacing and coalescence results are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies including spurious nodal configurational forces near the crack tip to improve the approximation of the total driving force."},{"cited_title":"Analysis of stress intensity factors of modes i, ii and iii for inclined surface cracks of arbitrary shape","cited_arxiv_id":null,"evidence_quote":"Provides the stress intensity factor relations for inclined cracks used to build the LEFM scaling estimate of the propagation angle $\\beta$."},{"cited_title":"N., Paluszny, A","cited_arxiv_id":null,"evidence_quote":"Shows how spatial arrangement of nearby cracks amplifies or shields stress intensity factors, supporting the paper's interpretation of facet spacing effects."}],"review_version":1}