{"id":"e813784c-f279-4d51-8aa6-d2146385413f","arxiv_id":"2411.15918","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Mechanical stability of arbitrarily loaded crystals of any symmetry reduces to checking that all Kelvin moduli of a stress-modified stiffness matrix in Mandel notation are positive.","lead":"A methods paper unifies mechanical stability checks for any 3D or 2D crystal under arbitrary load: build a stress-modified stiffness matrix in Mandel notation and test whether all its eigenvalues, the Kelvin moduli, are positive. The recipe includes explicit formulas and a Mathematica notebook, giving computational materials scientists a single algorithmic check instead of case-specific inequalities.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central criterion L=C+H may not give the correct tangent modulus for all loadings; the paper's Appendix B example undercuts its own recommendation by showing C and L predict different instabilities, yet Section 5 recommends checking only L.","rationale":"The reader's weakest assumption (that L=C+H is the correct stability criterion) is exactly the load-bearing premise I identify. I agree with the reader's CONDITIONAL verdict because the paper convincingly assembles the Mandel-notation spectral machinery and explicit stability tables, which are internally consistent where spot-checked, but it does not prove that the L criterion applies to all arbitrary loads; it merely cites prior work and selects L because it is stringent. Appendix B even concedes the need to check both C and L, weakening the universal recipe stated in Sections 4 and 5. My proposed test is concrete and decisive: reproduce the critical strain with an independent energy Hessian or phonon calculation using the same interatomic potential, and then compare to the C/L/A/Z predictions. This would settle whether the Kelvin-modulus test on L̃ is physically correct or merely one of several plausible criteria. No ad hominem is involved; the concern is about the argument's foundation, not the authors' competence.","tokens_in":19775,"tokens_out":1657,"duration_ms":14568,"concrete_test":"Re-execute the paper's Mathematica notebook for the Appendix B biaxial compression case and, in parallel, run molecular statics in LAMMPS with the same EAM potential (Pun-Mishin) to compute the phonon spectrum or the energy Hessian along the loading path; record the α at which a real phonon frequency or a negative curvature of the energy appears. Compare that onset α to the predictions of the C criterion (α=0.7985), the L criterion (α=0.9087), and the other criteria (A, Z). If the true instability occurs at or near α=0.9087, the L-based Kelvin recipe is vindicated; if it occurs near α=0.7985 or elsewhere, the Section 5 recipe is checking the wrong stability condition and the paper's central claim fails.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim (Sections 4-5) is that mechanical stability under arbitrary load is decided by positive definiteness of L = C + H (Eq. 17), i.e., positivity of the Kelvin moduli of the Mandel image L̃. The paper adopts the L criterion (Morris-Krenn [3]) because it is 'among the most stringent,' but competing criteria C, A, Z exist and the text explicitly notes in Appendix B that for biaxial compression L̃ becomes singular at α=0.9087 while C̃ alone becomes singular at α=0.7985, concluding 'it is recommended to check the condition on C̃ as well as on L̃.' That recommendation contradicts the Section 5 recipe, which states only that positivity of the Kelvin moduli of L̃ indicates stability. No rigorous derivation or validation shows that Eq. (17) is the physically relevant tangent modulus for arbitrary loads, arbitrary symmetry, and arbitrary loading paths, especially finite deformation states reached via different strain paths where C is configuration-dependent and H depends on the Cauchy stress. The paper does not prove L positive definiteness is necessary and sufficient for stability, nor does it specify when C (based on Green strain) is the correct criterion instead of L. Absent such proof, the algorithmic stability check may be testing the wrong quadratic form in exactly the finite-deformation, stressed regimes the paper targets. A direct check would recover the Appendix B numbers and, crucially, evaluate instability against phonon spectra or molecular-statics energy Hessians at the predicted critical strains, and recompute the critical α using each competing criterion (C, L, A, Z) to see which one matches the true instability onset.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified algorithmic procedure for testing the mechanical stability of 3D and 2D crystals under arbitrary load. It advocates the use of Mandel (orthonormal) notation, mapping fourth-order stiffness tensors to symmetric second-order tensors and checking the positivity of their eigenvalues, the Kelvin moduli. For stress-free crystals, it recovers explicit stability conditions for all symmetry classes (Appendices C and D). For deformed and stressed crystals, it adopts the incremental tangent modulus L = C + H (Eq. 17) and asserts that positivity of the Kelvin moduli of its Mandel image L tilde is the stability criterion. The manuscript also supplies a Mathematica notebook and an application to NiAl under biaxial deformation (Appendix B).","tokens_in":20048,"tokens_out":7711,"duration_ms":67176,"significance":"If the central claim is correct, the paper provides a practical, coordinate-free stability check that replaces lengthy lists of minorant conditions for stressed crystals of arbitrary symmetry with a single eigenvalue computation. The manuscript is strong in reproducibility: the explicit matrix representations of H tilde in Eqs. (21), (23), (26), and (28) are checkable, a Mathematica notebook is provided, and the stress-free Kelvin-modulus conditions are standard and well established. The main unresolved question is whether the chosen modulus L is the physically relevant tangent modulus for all finite-deformation regimes; the paper does not demonstrate this, and its own example in Appendix B shows that competing criteria differ. The contribution is therefore useful and promising, but the central criterion needs justification or careful reframing.","major_comments":[{"comment":"The paper's central claim that positivity of the Kelvin moduli of L tilde decides mechanical stability under arbitrary load rests on the unproven premise that L = C + H is the physically correct tangent modulus. The paper itself lists competing criteria C, A, and Z (Section 1) and justifies L only by stating that it is 'among the most stringent' and reproduces observed instabilities, citing Refs. [3,6]. Because the paper explicitly targets finite-deformation states where these criteria differ, as Appendix B demonstrates, this premise is load-bearing. A concrete validation would be to compute phonon spectra or molecular-statics energy second variations along the biaxial paths of Appendix B and to compare the predicted instability thresholds with the C tilde and L tilde criteria. Without such a comparison, the proposed algorithmic check may be evaluating a quadratic form that is not the relevant one for the stated scope.","section":"Section 4, Eq. (17), and Section 5"},{"comment":"The recommendation in Appendix B is inconsistent with the recipe in Section 5. For biaxial compression the text reports that L tilde becomes singular at alpha = 0.9087 while C tilde alone is not singular, and then recommends 'to check the condition on C tilde as well as on L tilde.' This contradicts the Section 5 statement that positivity of the Kelvin moduli of L tilde indicates mechanical stability. The paper must either adopt a two-criterion check in the conclusions or demonstrate, against a reference such as phonon spectra, that L is the correct criterion and C is not. As written, the reader cannot determine which criterion the paper actually endorses.","section":"Appendix B, Eqs. (B.6)-(B.9), and Section 5"},{"comment":"The assertion 'the crystal still remains mechanically stable' at alpha = 0.9087 is not justified. Since L tilde has a zero eigenvalue at this alpha, the Section 5 criterion would place the crystal at the boundary of stability. If the assertion is based on C tilde, it contradicts Section 5; if it is based on phonon or energy data, those data are not shown. This statement needs to be clarified or removed.","section":"Appendix B, text after Eq. (B.9)"}],"minor_comments":[{"comment":"The characters 'alpha > 1' and '0 < alpha < 1' are typeset with non-ASCII symbols ('alpha ¿ 1' and '0 ¡ alpha ¡ 1'); these should be fixed.","section":"Appendix B, Eq. (B.1)"},{"comment":"The denominator for the Reuss bulk modulus is ambiguous; explicit brackets should be added to clarify the division.","section":"Appendix E, Eq. (E.3)"},{"comment":"The norm equivalence for the elasticity tensor is stated without proof; citing the relevant result from Ref. [21] would aid readers who want to verify the isometry.","section":"Section 2, Eq. (14)"},{"comment":"The principal-minor conditions contain superscript notation such as 'C2_12' that is difficult to read; the typesetting of powers and indices should be corrected.","section":"Appendix C, orthotropic case"},{"comment":"The supplementary material is identified only through Reference [45]; giving the DOI or URL directly in the text would make the notebook easier to locate.","section":"Section 6"},{"comment":"The abstract says the paper unifies stability conditions 'under arbitrary load,' but the explicit conditions in Appendices C and D are for stress-free crystals; the wording should distinguish the general Section 5 recipe from the explicit stress-free lists.","section":"Abstract and Section 5"}],"recommendation":"major_revision","confidential_remarks":"This is essentially a methods and tutorial paper with a useful supporting notebook. The referee suggests the editor consider whether the unresolved criterion-selection issue is acceptable for the journal; it is fixable but currently undermines the advertised unification. The novelty relative to existing literature on Kelvin moduli and stressed stability lies mainly in the Mandel representation of the stress correction H and the automated procedure, which is modest but potentially useful. The internal inconsistency between Appendix B and Section 5 should be resolved before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful compilation that puts the L-criterion stability check in Mandel notation for 3D and 2D, with the explicit H tilde matrices as the genuinely new contribution. But the paper's headline recipe—check only the Kelvin moduli of L = C + H—contradicts its own Appendix B, which recommends checking C as well.\n\nWhat's new and good: Equations 23 and 28 give the stress-dependent tangent modulus in Mandel form explicitly; I don't know of another source that writes these out. Spot-checking against Eq. 17, the entries are consistent. The compiled conditions for all symmetry classes (Appendices C, D) are a handy reference, and the orientation-invariance demonstration in Appendix A is a nice illustration of why Mandel notation beats Voigt for this. The Mathematica notebook is a plus.\n\nSoft spots: First, the internal inconsistency. Section 5 states that positivity of the L-tilde eigenvalues indicates stability. Appendix B shows that for biaxial tension the C tilde singularity occurs at alpha=1.15365 versus alpha=1.15454 for L tilde, and for compression L tilde fails first (0.9087 vs 0.7985). The appendix then recommends checking both C and L. That recommendation is sensible, but it undercuts the unified \"check L only\" recipe. Second, the paper does not justify why L is the physically relevant tangent modulus for all loads and symmetries; it cites Morris-Krenn and appeals to stringency, but competing criteria exist (C, A, Z) and the paper's own example shows L isn't always the most stringent. Third, the predicted critical stretches are never validated against phonon spectra or energy Hessians, so the unified criterion remains plausible but untested in the finite-deformation regime it targets.\n\nWho it's for: computational materials scientists doing atomistic or ab initio work on loaded crystals. It's a practical tool paper, not a foundational one. I'd send it to review: the explicit matrices and compiled conditions are worth publishing. Revisions should address the criterion question and add at least one validation of a predicted instability.\n\nRecommendation: engage with it, but expect revisions.","headline":"Useful Mandel-notation stability compilation, but the paper's own Appendix B contradicts its L-only recipe.","tokens_in":20640,"tokens_out":4178,"would_cite":true,"duration_ms":37601,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74B05","74E10"],"pacs":["62.20.de","61.50.Ah"],"model":"deepseek-v4-flash","headline":"For any 3D or 2D crystal of any symmetry under any load, mechanical stability is decided by the positivity of the Kelvin moduli of the load-corrected stiffness tensor in Mandel notation.","keywords":["mechanical stability","Born stability","Kelvin moduli","orthonormal (Mandel) notation","2D materials","stressed crystals","elastic stability criteria"],"falsifier":"Drive a well-characterized crystal, say biaxially compressed B2 NiAl, through the predicted threshold alpha of about 0.9087 with a finite-strain molecular-dynamics or phonon calculation and watch for the first instability: if a softening mode or bifurcation appears while all Kelvin moduli of L-tilde are still positive, or fails to appear after one has turned negative, the claimed criterion is not the physically operative one. The comparison is concrete because the paper already supplies the stress state and stiffness at the critical stretch.","tokens_in":19538,"feed_emoji":"⚖️","tokens_out":11992,"duration_ms":95985,"temperature":0.7,"pith_summary":"This paper claims that one algorithmic recipe settles mechanical stability for every crystal: 3D or 2D, of any point-group symmetry, stress-free or under arbitrary load. The recipe is to write the fourth-order stiffness tensor together with a stress-dependent correction as a symmetric second-order tensor in orthonormal (Mandel) notation, then check that all its eigenvalues, the Kelvin moduli, are positive. This replaces the current practice of laborious, orientation-dependent lists of principal-minor conditions that differ for each symmetry class and cannot easily handle loaded crystals. If the claim is right, stability analysis becomes a standard numerical eigenproblem that any atomistic or continuum user can run directly from computed elastic constants and the Cauchy stress.","feed_headline":"One eigenproblem settles stability of any stressed crystal","feed_subtitle":"The same load-corrected stiffness tensor with positive Kelvin moduli works for every symmetry class, 3D and 2D.","key_machinery":"The load-bearing object is the symmetric second-order tensor L-tilde = C-tilde + H-tilde, assembled from the Mandel-notation images of the stiffness tensor C and of the stress correction H of Eq. (17). Mandel (orthonormal) notation is the map that rewrites a fourth-order stiffness tensor as a real symmetric 6x6 matrix in 3D or 3x3 matrix in 2D, with sqrt(2) factors on the shear entries, so that the matrix entries transform as a genuine second-order tensor; this is what makes the eigenvalues, the Kelvin moduli, orientation-independent, unlike Voigt-notation principal-minor conditions. H carries the explicit stress terms of the Morris-Krenn symmetrized tangent modulus conjugated to Cauchy stress. The argument runs: positive definiteness of the quadratic form delta-epsilon L delta-epsilon is the internal stability condition; L-tilde is symmetric, so Sylvester's criterion is replaced by the equivalent and simpler demand that all its eigenvalues be positive; and for 2D problems those eigenvalues can even be written in closed form via Cardano formulas.","core_discovery":"The paper's central claim is that internal mechanical stability of a crystal under arbitrary load reduces to one numerical test: positive definiteness of the symmetrized incremental tangent modulus L = C + H of Eq. (17), where C is the stiffness tensor in the current deformed configuration and H is the stress-dependent correction. In Mandel notation this is the statement that all eigenvalues of the 6x6 (or 3x3 in 2D) symmetric matrix L-tilde = C-tilde + H-tilde, called the Kelvin moduli, are positive. Because the Mandel representation maps fourth-order stiffness tensors to genuine second-order tensors while preserving norms, the eigenvalue test is objective: the Kelvin moduli do not depend on the crystal orientation in which the constants were computed, as the three differently oriented NiAl representations in Appendix A demonstrate. The worked analysis of B2 NiAl in Appendix B shows the test in action, with L-tilde becoming singular at stretch alpha of about 1.15454 in biaxial tension and 0.9087 in biaxial compression, while the bare tensor C-tilde alone turns singular respectively slightly earlier and far later, which is why the paper recommends checking both.","pith_inferences":["A next step the paper leaves implicit is to pair the L-tilde Kelvin-modulus test with finite-strain phonon calculations on the same deformed configuration; wherever the two thresholds disagree, that discrepancy itself maps where Born-type mechanical stability and dynamical stability part ways.","Because the minimum Kelvin modulus of L-tilde is a smooth scalar function of the deformation, it could double as an order parameter in high-throughput searches or strain-engineering protocols: maximizing the minimum eigenvalue over loading paths is a well-defined target for designing metastable states.","The paper's own example hints at a load-direction asymmetry, with C-tilde and L-tilde thresholds nearly coinciding in biaxial tension but differing strongly in compression, suggesting that the choice among the four competing criteria (C, A, Z, L) may itself depend on the loading path; a systematic scan over stress states for a single crystal could test whether the L criterion is always the governi","The recipe is zero-temperature and static; extending it to finite temperature through a temperature-dependent effective stiffness tensor would make the Kelvin-modulus criterion a candidate for a unified finite-temperature stability check."],"forward_implications":["Stability checking becomes a single numerical eigenproblem: feed the computed elastic constants and the Cauchy stress of the current configuration into L-tilde = C-tilde + H-tilde and test the sign of the smallest eigenvalue, for any material symmetry and any load.","The result is objective: because C-tilde and H-tilde are genuine tensors in six- or three-dimensional space, the Kelvin moduli are invariant under rotation of the crystal axes, so calculations in non-conventional or rotated computational cells need no special handling.","The same recipe covers 2D crystals with a 3x3 eigenproblem whose eigenvalues can be given in closed form, so two-dimensional materials are handled with the same code path as bulk crystals.","Explicit stability inequalities for every 3D symmetry class and every 2D class follow as special cases, unifying the scattered literature conditions, and the accompanying computational notebook turns the check into a routine procedure.","For stressed crystals the criterion automatically includes load-induced contributions through H, so pressure-, tension-, and shear-driven instabilities are detected by the same test that governs stress-free crystals."],"supporting_citations":[{"why":"Defines the symmetrized incremental tangent modulus L = C + H, the tensor whose positive definiteness the whole test checks.","marker":"[3]"},{"why":"Establishes the orthonormal (Mandel) notation and the eigentensor decomposition that make the Kelvin-moduli test tensorially well posed.","marker":"[21]"},{"why":"Original source of the term Kelvin moduli for the eigenvalues whose positivity is the stability criterion.","marker":"[28]"},{"why":"Shows the eigenvalue problem for Hooke's tensor is well posed, grounding the spectral-decomposition route.","marker":"[25]"},{"why":"The classical stability problem for unstressed crystal lattices that this work generalizes to arbitrary symmetry and load.","marker":"[1]"},{"why":"Supplies the standard principal-minor stability conditions in Voigt notation that the new eigenproblem replaces.","marker":"[2]"},{"why":"Basis for choosing the L criterion: it is among the two most stringent proposed criteria and reproduces observed instabilities.","marker":"[6]"},{"why":"Earlier application of the L-type criterion to homogeneous crystals under finite strain, precursor of Eq. (17).","marker":"[40]"},{"why":"Source of the B2 NiAl stiffness data used in the orientation-invariance and biaxial-deformation worked examples.","marker":"[36]"}],"fun_headline_variants":["One matrix test decides stability of any crystal under load","Kelvin moduli unify stability for all 3D and 2D crystals","Stability under arbitrary load: positive Kelvin moduli suffice","Crystal stability, any symmetry, any stress: one eigenproblem","From 2D to 3D, one condition checks crystal stability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The test stands or falls on treating positive definiteness of the symmetrized tangent modulus L = C + H as the correct definition of internal stability under load; competing criteria C, A, and Z are also in circulation, and for a stressed crystal they need not agree.","fun_headline_variants_meta":{"raw":{"variants":["One matrix test decides stability of any crystal under load","Kelvin moduli unify stability for all 3D and 2D crystals","Stability under arbitrary load: positive Kelvin moduli suffice","Crystal stability, any symmetry, any stress: one eigenproblem","From 2D to 3D, one condition checks crystal stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000759,"raw_usage":{"total_tokens":3340,"prompt_tokens":882,"completion_tokens":2458,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":498,"completion_tokens_details":{"reasoning_tokens":2367}},"tokens_in":498,"tokens_out":2458,"duration_ms":17903,"temperature":1.0,"reasoning_tokens":2367,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:45:51.397719+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Drive a well-characterized crystal, say biaxially compressed B2 NiAl, through the predicted threshold alpha of about 0.9087 with a finite-strain molecular-dynamics or phonon calculation and watch for the first instability: if a softening mode or bifurcation appears while all Kelvin moduli of L-tilde are still positive, or fails to appear after one has turned negative, the claimed criterion is not the physically operative one. The comparison is concrete because the paper already supplies the stress state and stiffness at the critical stretch.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the symmetrized incremental tangent modulus L = C + H, the tensor whose positive definiteness the whole test checks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the orthonormal (Mandel) notation and the eigentensor decomposition that make the Kelvin-moduli test tensorially well posed."},{"cited_title":"Thomson, XXI","cited_arxiv_id":null,"evidence_quote":"Original source of the term Kelvin moduli for the eigenvalues whose positivity is the stability criterion."},{"cited_title":"Rychlewski, On Hooke’s law, Journal of Applied Mathematics and Me- chanics 48 (3) (1984) 303–314","cited_arxiv_id":null,"evidence_quote":"Shows the eigenvalue problem for Hooke's tensor is well posed, grounding the spectral-decomposition route."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier application of the L-type criterion to homogeneous crystals under finite strain, precursor of Eq. (17)."}],"review_version":1}