{"id":"b82e59de-028d-4fa6-b427-6f52aa3646c3","arxiv_id":"2508.00193","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors combine edge-based smoothed finite elements with element splitting and a local energy-release-rate criterion to track dynamic crack growth in 2D, validating on Kalthoff plates, concrete beams, and PMMA specimens.","lead":"This paper describes a finite element technique, the Crack Element Method, for simulating how cracks spread through two-dimensional structures under fast dynamic loads. It matters for chip reliability because it aims to predict crack paths in copper and low-k interconnect stacks without expensive remeshing.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (24) defines G_G0 as σ⊥·δd/2 without deriving it from the variation of potential energy with crack area; if this local product does not equal the actual energy release rate under mixed-mode dynamic loading, the crack-advance criterion and the dissipated-energy comparisons lose their stated…","rationale":"I read the paper as proposing a practical element-splitting dynamic fracture scheme and claiming that its local criterion, Eq. (24), is a fracture energy release rate. For that claim to hold, Eq. (24) must equal -∂Π/∂A, or at least approximate it well enough that comparing with Gc controls crack advance. The manuscript never shows that. The expression is dimensionally consistent, but it multiplies a stress at a candidate point with the current stretch of the crack-tip edge; these are not conjugate for the same crack increment. The paper also concedes in Section V that the formulation is mode-I based, while two of the three benchmarks are mixed mode. This is the weakest load-bearing point because every crack-direction/timing decision and every dissipated-energy comparison in Section IV uses Eq. (24). The validation matches experiments and prior simulations with no fitted parameters, which is real supporting evidence; the coarse-mesh Kalthoff result is particularly encouraging. But good agreement cannot by itself establish that the criterion is an energy release rate. The proposed test is straightforward: compute an independent J-integral from the code's own fields at first advance and compare. If it matches, the paper needs only a derivation; if it does not, the central energetic claim is unsupported and the method should be presented as a heuristic. The reader's weakest_assumption is the same, so I agree. I do not move the verdict: the paper should remain CONDITIONAL, pending the derivation/verification, not be rejected, because the benchmark performance is credible and no parameter fitting is involved.","tokens_in":17084,"tokens_out":6351,"duration_ms":68249,"concrete_test":"Re-run the Kalthoff benchmark with the fine irregular mesh, freeze the crack at the first advance step, and compute the dynamic J-integral (or equivalent domain integral) from the same ES-FEM stress/displacement fields around the crack tip. Compare this J with Eq. (24)'s σ⊥δd/2 at the chosen candidate edge. Repeating at several early advance steps and for the PMMA compact-compression case would settle whether the local product is the true energy-release rate; a deviation consistently larger than ~20% would mean the crack-advance decisions are not energetically grounded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III presents Eq. (24), G_G0 = σ⊥_Gi δd/2, as 'a formulation of fracture energy release rate,' but no derivation from δL/δA or an equivalent crack-closure integral is given. The quantity has the right dimensions, yet conjugacy is missing: δd is the current stretch of the crack-tip edge (Eq. 23), while σ⊥_Gi is the normal projection of maximum principal stress at a candidate edge point ahead of the tip. These are not the traction and opening associated with the same virtual crack extension; a true local energy-release expression would involve an integral of σ(r)δ(Δa−r) over the would-be crack segment, not a single product of an ahead-of-tip stress with an already-existing opening. In addition, Eq. (24) uses only the normal projection of maximum principal stress and therefore contains no explicit mode-II or mode-III contribution, although the Kalthoff and PMMA benchmarks are mixed-mode; the conclusion itself states the formulation is 'based solely on mode-I fracture mechanics.' Since every crack-advance decision compares this quantity with Gc, the method's energy basis and its predicted direction/timing both hinge on Eq. (24) being a true energy release rate. If it is only a dimensionally suggestive heuristic, the good benchmark agreement may be fortuitous or mesh-tuned, and the claim of an 'energy-based formulation' is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a Crack Element Method (CEM) implemented within the Edge-based Smoothed Finite Element Method (ES-FEM) framework. The method tracks dynamic crack growth in 2D quasi-brittle solids by splitting elements along edge-quadrature-point paths and by using a locally computed fracture energy release rate (Eq. 24) to decide whether and where a crack advances. The approach is demonstrated on three classical dynamic fracture benchmarks (Kalthoff-Winkler plate, three-point bending concrete beam, compact-compression PMMA specimen) and on two Cu/Ultra Low-k interconnect case studies, with comparisons to experimental data and reference numerical solutions for crack paths, dissipated energy, and crack tip velocity.","tokens_in":17410,"tokens_out":2724,"duration_ms":28459,"significance":"If the central criterion is accepted, the paper offers a practically attractive alternative to phase-field, XFEM, and cohesive-zone methods: it avoids remeshing and damage smearing, works on coarse meshes, and uses material fracture energies taken from the literature rather than fitted to the validation benchmarks. The benchmark coverage, including the Kalthoff crack angle near 70°, the transition parameter γ_t = 0.765 versus the experimental 0.77, and the curved PMMA crack path, is a genuine strength. However, the energetic meaning of the crack advance criterion is not established, and the reported crack-tip velocities show noteworthy mesh dependence. The method's practical value is plausible, but the theoretical status of Eq. (24) must be clarified before the central claim of an 'energy-based formulation' can be accepted.","major_comments":[{"comment":"The fracture energy release rate G_G0 = σ⊥_Gi · δd / 2 is presented as 'derived' from the topology of split elements, but no derivation from the variation of potential energy with respect to crack area (e.g., δL/δA or a crack-closure integral) is given. The quantity has the correct dimensions, yet the conjugacy is missing: δd (Eq. 23) is the stretch of the current crack-tip edge, while σ⊥_Gi is the normal projection of the maximum principal stress at a candidate edge quadrature point ahead of the tip. These are not the traction and opening associated with the same virtual crack extension. Since every crack advance decision in Eqs. (25) and (26) compares this quantity with G_c, and since the dissipated-energy comparisons in Figs. 7(a) and 13(a) inherit its meaning, the method's energy basis and its predicted direction and timing both hinge on Eq. (24) being a true energy release rate. The authors should either provide a proper derivation (e.g., from the discrete variation of the total potential energy when an element is split) or explicitly state that Eq. (24) is a heuristic criterion and then treat the benchmark agreement as empirical validation of that heuristic. As written, the claim that CEM is an 'energy-based formulation' is not supported.","section":"§III, Eq. (24)"},{"comment":"Eq. (24) contains only the normal projection of the maximum principal stress and therefore has no explicit mode-II contribution. The Kalthoff-Winkler, three-point bending, and compact-compression problems are all mixed-mode to varying degrees, and the conclusion itself acknowledges that the formulation is based solely on mode-I fracture mechanics. If Eq. (24) is intended to represent the total energy release rate in mixed-mode conditions, the absence of shear traction work must be justified. If it is a mode-I-only criterion, then the good agreement with mixed-mode benchmark paths needs an explanation beyond the assertion in the conclusion, because the crack direction is selected by maximizing σ⊥ among candidate points, which is a maximum-principal-stress direction condition rather than an energy-balance direction condition.","section":"§III, Eq. (24) and benchmarks IV.A–IV.C"},{"comment":"The crack-tip velocity comparison shows a material mesh dependence: the regular fine mesh crack-tip velocity decreases after 35 µs, while the coarse and irregular meshes increase to about 1600 m/s. The authors offer a qualitative explanation (straighter path in the regular mesh), but the claim that CEM 'accurately captures' crack tip velocities requires a quantitative assessment: what is the spread relative to the reference solution, and is a 1600 m/s versus a decreasing velocity within the expected accuracy of the method? At minimum, the authors should provide the reference velocity data from [30] in the same plot and quantify the discrepancy.","section":"§IV.A, Fig. 7(b)"},{"comment":"Equation (17) is labeled as a 'diagonal mass matrix' but the displayed expression is a consistent (non-diagonal) element mass matrix, since it is the outer product of shape function vectors integrated over the element. If a diagonal/lumped mass matrix is used in the explicit time integration of Eq. (20), the lumping scheme must be described or referenced. As written, the formulation is internally inconsistent: a non-diagonal mass matrix would not allow the 'without equation solving' statement following Eq. (20).","section":"§II.B, Eq. (17)"}],"minor_comments":[{"comment":"The definition of δd uses the Heaviside function H(·) of (‖x_N2 − x_N1‖/‖X_N2 − X_N1‖ − 1). When the edge is in compression, δd = 0, so the criterion predicts no crack advance even under large compressive stretch along the edge. This should be stated and justified, as it affects the interpretation of the crack path under mixed-mode loading.","section":"§III, Eq. (23)"},{"comment":"The notation in Eq. (11) is confusing: the symbols N_gi^e, L_gi^e, M_Aj^n, and A_j are introduced in the text but the equation combines them in a way that is hard to parse. Please define all symbols directly below the equation and clarify the meaning of the quotient term (the area weight).","section":"§II.B, Eq. (11)"},{"comment":"There is a typo 'to to G3' in the description of crack pattern II in the QUAD element; also, the phrase 'from G0 to to G3' should read 'from G0 to G3.'","section":"§III, Fig. 3(b) and text"},{"comment":"The text states that the crack angle in the proposed method is 'around 65°∼70°,' while the experimental value is 70°. Given that the crack path is digitized from figures, please provide a quantitative measurement of the final crack angle for each mesh and state the precise deviation from 70°.","section":"§IV.A"},{"comment":"The qualitative comparison between the simulated crack pattern and the fab image in Fig. 15(b) is not sufficient to 'confirm the ability of the proposed CEM to model the fracturing process.' A quantitative metric (e.g., crack path length, deviation from the layer interface, or comparison with a reference simulation) would strengthen the claim.","section":"§IV.D, mechanically-induced loads"},{"comment":"The manuscript contains several typographical errors: 'Possion' for 'Poisson,' 'detedction' for 'detection,' 'the reason why more energy is dissipated... may come from discrepancy' should be 'may come from the discrepancy,' and 'making it as a strong alternative' should be 'making it a strong alternative.' A careful proofread is recommended.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper's practical contribution is clear and the benchmark results are encouraging, but the theoretical status of Eq. (24) is the main risk: it is presented as a derivation when it is at best a heuristic, and the mixed-mode benchmarks are used to validate a criterion that is explicitly mode-I only. If the authors can either provide a proper energy derivation or, failing that, explicitly reframe Eq. (24) as an empirical fracture criterion, the paper would be defensible. The mesh-dependence of the crack-tip velocity in Fig. 7(b) also needs a quantitative discussion. I do not see evidence of circularity or inappropriate fitting of parameters to the benchmarks, which is a point in the paper's favor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper reports a Crack Element Method (CEM) built on ES-FEM with an element-splitting scheme. The combination is genuinely new: crack tips live at edge quadrature points, the crack advances by splitting elements, and the advance criterion is a local comparison of a stress times a stretch against G_c. The benchmarks are relevant and the practical results are decent: the Kalthoff angle comes out near 70°, the critical notch offset is γ_t = 0.765 versus the experimental 0.77, and the curved PMMA paths look right. No parameters were fitted to these benchmarks; G_c values come from experiments and prior material data. That is real work and deserves credit.\n\nThe soft spot is the load-bearing criterion, Eq. (24). The paper calls it a derivation of the fracture energy release rate, but it is an ansatz. The product of the normal-projected maximum principal stress at a candidate edge point and the stretch of the current crack-tip edge has the right dimensions but no demonstrated conjugacy with the energy released per new crack area. The stress and the opening are not from the same virtual crack extension. And since only the normal component is used, the criterion is purely mode-I in a mixed-mode setting; the paper's own conclusion acknowledges this. This means the good benchmark agreement could be fortuitous or mesh-tuned. I do not think it is intentionally tuned—G_c is fixed and the meshes differ—but the lack of a derivation leaves a real gap.\n\nThere are also smaller issues. Eq. (17) labels a diagonal mass matrix while displaying a consistent (non-diagonal) formula; that needs fixing. The mesh-independence claims are stronger than the evidence: crack paths are confined to element edges, so orientation independence is inherently limited, and the crack-tip velocity curves in Fig. 7(b) differ noticeably across meshes. The PMMA benchmark is mostly qualitative; the reference path is compared visually. No code or data is provided, which limits reproducibility.\n\nThe paper is worth engaging with. The method is practical and the benchmarks suggest it works, but the central criterion needs either a proper derivation or a careful justification as a heuristic with a study of its limits. Without that, the energy-based claim is unsupported. A serious referee could usefully require those additions.\n\nFor peer review: send it to review, but with a request for major revision focusing on the criterion and quantitative comparisons. It deserves referee time.","headline":"A practical element-splitting crack method with good benchmark fits, but the central energy-release-rate criterion is asserted rather than derived, so the energy basis is weakly supported.","tokens_in":17899,"tokens_out":2401,"would_cite":false,"duration_ms":24016,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper introduces the Crack Element Method, an ES-FEM-based element-splitting scheme whose local fracture-energy-release-rate criterion reproduces dynamic crack angles, dissipated energy, and crack-tip velocities in quasi-brittle solids.","keywords":["dynamic fracture","Crack Element Method","ES-FEM","element-splitting algorithm","fracture energy release rate","quasi-brittle materials","Cu/ultra low-k interconnect","chip-package interaction"],"falsifier":"Compute, in one simulated split, the actual drop in total potential energy of the entire domain divided by the area of the newly created crack surface, and compare it with the value of Eq. (24) that triggered the split; if the two disagree beyond discretization error, the energy-release-rate interpretation is falsified.","tokens_in":16878,"feed_emoji":"💥","tokens_out":11874,"duration_ms":105858,"temperature":0.7,"pith_summary":"This paper puts forward the Crack Element Method, a finite-element scheme for tracking dynamic cracks in two-dimensional quasi-brittle solids. The method splits elements along their edges rather than adding new nodes or smearing damage over a region, and it decides where the crack goes by comparing a locally computed fracture energy release rate with the material's critical value. The paper demonstrates on three benchmark problems that the scheme reproduces experimental crack angles, curved crack paths, dissipated fracture energy, and crack-tip velocities, and it shows a patterned Cu/ultra low-k interconnect layout cracking under mechanical and thermal loads. The claim is that this topological, energy-based rule keeps accuracy on coarse and irregular meshes and extends naturally to three-dimensional problems.","feed_headline":"Element-splitting FEM matches dynamic crack benchmarks","feed_subtitle":"A local energy-release-rate rule predicts kinked, curved crack paths and works on chip interconnect layouts.","key_machinery":"The load-bearing mechanism is the element-splitting algorithm implemented on the Edge-based Smoothed Finite Element Method (ES-FEM), in which strains are averaged over edge-based smoothing domains. Crack tips are identified with edge quadrature points; at each step the method enumerates the two candidate next tip locations in a triangle or the three in a quadrilateral, computes the candidate energy release rates through Eq. (24), and splits the element along the winning edge. The fracture criterion is the product of the normal-projected maximum principal stress at the candidate point and half the Heaviside-weighted stretch of the current crack-tip edge, so the only bookkeeping required is the evolving topology of split elements. This topology-based rule replaces the J-integral and turns crack advance into a local comparison with $G_c$.","core_discovery":"The central discovery is that dynamic crack growth in quasi-brittle solids can be driven by a purely local, element-topology-based energy criterion instead of a path integral or a diffused damage field. When the crack tip sits at an edge quadrature point $G_0$ of an ES-FEM mesh, the method evaluates, for each candidate next edge point $G_i$, the quantity $G_{G_0} = \\sigma^{\\perp}_{G_i}\\delta_d / 2$, where $\\delta_d$ is the stretch of the current crack-tip edge and $\\sigma^{\\perp}_{G_i}$ is the component of the maximum principal stress at $G_i$ perpendicular to the candidate crack line. If $G_{G_0}$ exceeds the material's critical fracture energy release rate $G_c$, the crack advances to $G_i$; when a quadrilateral element is only partially split, it degenerates into triangle elements that continue the computation. The paper reports that this rule, without enrichment functions or a cohesive law, gives a shear-impact crack angle near 70 degrees, a critical notch-offset threshold of 0.765 versus the 0.77 measured in the concrete beam experiment, and curved PMMA crack paths matching reference simulations, with crack-tip speeds staying below roughly 60 percent of the Rayleigh wave speed.","pith_inferences":["Although Eq. (24) is built from a mode-I-style normal stress, it performs well on mixed-mode benchmarks; this suggests the opening component dominates in those cases, and a future sliding-mode term may be needed for shear-dominated crack growth.","Because the method tracks only element topology, it could be combined with adaptive refinement of the band of elements near the crack tip, reducing cost further while preserving local stress resolution.","The formula resembles a crack-closure estimate that needs no predefined crack path; if it survives broader validation, it could serve as a cheap screening tool for comparing interconnect layouts and crack-stop structures in chip-package interaction studies.","The interconnect demonstration uses a periodic, conceptual layout, which hints that the method is cheap enough for statistical studies over many layout variants."],"forward_implications":["A coarse triangle mesh with 434 nodes still produces a shear-impact crack path near 70 degrees, so the method does not require the fine meshes that smeared crack approaches typically need.","The critical notch offset 0.765 reproduces the measured 0.77 transition between pure mode-I and mixed-mode cracking in the concrete beam, without an artificially inserted mid-span notch.","In the PMMA compact compression case the curved arc-shaped crack path and the crack-tip speed history, peaking near 700 m/s and staying below about 60 percent of the Rayleigh wave speed, are reproduced with three different mesh densities.","The same formulation captures mechanically and thermally induced crack propagation in a patterned Cu/ultra low-k interconnect layout, matching the qualitative failure patterns used as reference."],"supporting_citations":[{"why":"It provides the edge-based smoothed finite element formulation whose edge smoothing domains the element-splitting algorithm is built on.","marker":"[19]"},{"why":"It supplies the shear-impact experimental dataset defining the benchmark crack angle the method matches.","marker":"[26]"},{"why":"It supplies the concrete beam impact experiment and the reported critical notch-offset value 0.77 used for comparison.","marker":"[31]"},{"why":"It supplies the PMMA compact compression experiment that defines the curved-crack benchmark.","marker":"[35]"},{"why":"It supplies the PMMA material properties used in the compact compression simulations.","marker":"[37]"},{"why":"It supplies the reference curved crack path against which the simulated PMMA crack paths are compared.","marker":"[39]"},{"why":"It supplies the reference dissipated-energy and crack-tip-velocity curves used for quantitative comparison in the shear-impact benchmark.","marker":"[30]"},{"why":"It supplies the Cu/ultra low-k interconnect crack-pattern reference used for the chip-package interaction case studies.","marker":"[5]"}],"fun_headline_variants":["Local edge energy rule predicts dynamic crack growth","Element-splitting method simulates cracks without cohesive laws","Split-element criterion drives fracture in chip interconnects","Dynamic crack paths from local energy, no path integrals","CEM: local rule for dynamic cracks in ultra low-k chips"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that Eq. (24) really measures the energy released per unit of new crack surface when the crack advances; the paper asserts this equality rather than deriving it from the variation of potential energy with crack area, and if it fails under mixed-mode dynamic loading the crack-advance decisions lose their energetic basis.","fun_headline_variants_meta":{"raw":{"variants":["Local edge energy rule predicts dynamic crack growth","Element-splitting method simulates cracks without cohesive laws","Split-element criterion drives fracture in chip interconnects","Dynamic crack paths from local energy, no path integrals","CEM: local rule for dynamic cracks in ultra low-k chips"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00171,"raw_usage":{"total_tokens":6762,"prompt_tokens":932,"completion_tokens":5830,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":5752}},"tokens_in":548,"tokens_out":5830,"duration_ms":39250,"temperature":1.0,"reasoning_tokens":5752,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T10:19:40.554190+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, in one simulated split, the actual drop in total potential energy of the entire domain divided by the area of the newly created crack surface, and compare it with the value of Eq. (24) that triggered the split; if the two disagree beyond discretization error, the energy-release-rate interpretation is falsified.","supporting_citations":[{"cited_title":"An edge-based smoothed finite element method (es-fem) for static, free and forced vibration analyses of solids,","cited_arxiv_id":null,"evidence_quote":"It provides the edge-based smoothed finite element formulation whose edge smoothing domains the element-splitting algorithm is built on."},{"cited_title":"Failure mode transition at high rates of shear loading,","cited_arxiv_id":null,"evidence_quote":"It supplies the shear-impact experimental dataset defining the benchmark crack angle the method matches."},{"cited_title":"Mixed-mode fracture of concrete subjected to impact loading,","cited_arxiv_id":null,"evidence_quote":"It supplies the concrete beam impact experiment and the reported critical notch-offset value 0.77 used for comparison."},{"cited_title":"An investigation of dynamic crack initiation in pmma,","cited_arxiv_id":null,"evidence_quote":"It supplies the PMMA compact compression experiment that defines the curved-crack benchmark."},{"cited_title":"A general mass lumping scheme for the variants of the extended finite element method,","cited_arxiv_id":null,"evidence_quote":"It supplies the PMMA material properties used in the compact compression simulations."},{"cited_title":"Efficient explicit time stepping for the extended finite element method (x-fem),","cited_arxiv_id":null,"evidence_quote":"It supplies the reference curved crack path against which the simulated PMMA crack paths are compared."},{"cited_title":"Simulation of dynamic brittle and quasi-brittle fracture: a revisited local damage approach,","cited_arxiv_id":null,"evidence_quote":"It supplies the reference dissipated-energy and crack-tip-velocity curves used for quantitative comparison in the shear-impact benchmark."},{"cited_title":"Chip-package interaction and crackstop study for cu/ultra low-k interconnects,","cited_arxiv_id":null,"evidence_quote":"It supplies the Cu/ultra low-k interconnect crack-pattern reference used for the chip-package interaction case studies."}],"review_version":1}