{"id":"8c181b52-c1e5-42fa-9b73-7019d53cec73","arxiv_id":"2607.13734","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The stress-strain response of metallic and polymer glasses—overshoot, yielding, flow, and hardening—is modeled by adding irreversible relaxation and chain extensibility to nonaffine elasticity.","lead":"A new model combines nonaffine elasticity with slow structural relaxation to predict full stress-strain curves for metallic and polymer glasses, including the yield peak and later strain hardening. If it holds up, it offers a common microscopic language for deformation in amorphous materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Yield overshoot position and height are controlled by an imported exponential coordination-decay law (Eq. 1) that is never independently measured or derived here; a direct nb(gamma) test is needed.","rationale":"The reader and I identify the same weakest assumption: the exponential coordination-decay law. This is the single point on which the microscopic claim rests. In the paper's favor, the model has sensible limiting behavior for the elastic term alone, and beta is extracted independently for the metallic glass; the viscous and hardening terms are standard forms. But those do not test the origin of the overshoot. The yield peak appears because the elastic stress is sigma_el proportional to gamma * exp(-A*gamma) * (2 - A*gamma); the maximum is at gamma = 1/A. Any other strain dependence of nb(gamma) changes the predicted position and height, and the paper provides no measurement of nb(gamma). A direct simulation check is feasible and would settle the question. Because this concern is the same one the reader used to justify a conditional verdict, my read does not change the verdict.","tokens_in":10354,"tokens_out":9843,"duration_ms":106058,"concrete_test":"Perform constant-strain-rate MD or athermal quasistatic shear simulations of a model glass (e.g., 80/20 binary Lennard-Jones or CuZr) and measure the mean number of mechanically active nearest-neighbor contacts nb as a function of shear strain gamma. Fit the measured decay to nb(gamma) = n0/2 * (1 + exp(-A*gamma)) and compare the fitted A with the model's A = Tg/T + 1/(eps_dot*tau_c). If the decay is not exponential, or if the fitted A does not follow the stated strain-rate dependence, Eq. (12) and the predicted overshoot epsilon_y ~ 1/A are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The predicted overshoot is generated entirely by the factor gamma * exp(-A*gamma) * (2 - A*gamma) in Eq. (1) (uniaxial form Eq. (12)); the yield strain is epsilon_y ~ 1/A with A = Tg/T + sqrt(3)/(2(1+nu)*eps_dot*tau_c) (Eq. 11). Hence the central claim—that overshoot is a nonaffine elastic instability from strain-driven connectivity loss—is mathematically equivalent to the assumed exponential decay nb(gamma) = n0/2 * [1 + exp(-A*gamma)], imported from Ref. [36] but not derived or measured in this paper. If the true coordination-loss kinetics under strain is not exponential (e.g., sigmoidal or cooperative), the predicted overshoot position and height change. Since tau_c is the main adjustable timescale, the agreement with the four materials shown in Figs. 2-5 may reflect fitting the assumed functional form rather than supporting the microscopic mechanism. This is the load-bearing premise: without a direct check of nb(gamma), the connection between relaxation dynamics and the yield peak is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a closed-form constitutive model for amorphous solids in which the total stress is a sum of a nonaffine elastic term, a stretched-exponential viscous term, and (for polymers) a Langevin-based strain-hardening term. The elastic term uses a strain-dependent coordination number nb(γ) that decays exponentially with strain, which produces a stress maximum interpreted as yielding/overshoot. The authors fit the model to experimental stress–strain curves for a metallic glass, polycarbonate, PMMA, and epoxy over a range of strain rates and report good agreement. The central claim is that stress overshoot originates from the competition between a nonaffine elastic instability, induced by strain-driven loss of mechanical connectivity, and viscous dissipation, and that a single microscopic framework captures elasticity, yielding, flow, and hardening.","tokens_in":10675,"tokens_out":7686,"duration_ms":80179,"significance":"If established, the framework would provide a rare unified, analytically tractable description of nonlinear deformation in amorphous solids, connecting relaxation dynamics to the entire stress–strain curve. The paper has clear strengths: explicit equations, recovery of expected limits (with the caveat noted below), a demonstrated independent extraction of β for the metallic glass, and comparisons against data for four material classes spanning many strain-rate decades. However, the central mechanistic conclusion is not yet independently supported: the overshoot is generated by an imported exponential coordination-decay law, and the validation relies on per-material, per-rate fitting parameters. The paper would be significant if the coordination-decay law were derived or directly verified, and if the independent-input claim for polymers were substantiated.","major_comments":[{"comment":"Eq. (7) with β<1 produces a non-analytic viscous stress: for small ε, σvisc(ε) ≈ η εdot^{1−β} τv^{−β} ε^β. Since β<1, the tangent modulus dσvisc/dε diverges as ε→0. Therefore Eq. (10) does not reduce to linear elasticity as claimed; the total stress has a cusp at the origin. This affects the interpretation of the fitted E and the stated limiting behavior. Please repair the viscous term (e.g., regularize it) or explicitly justify that the fitting range avoids the singular regime and that the claim is only asymptotic for β=1.","section":"§2, Eq. (7) and text after Eq. (12)"},{"comment":"The stress overshoot is produced by the assumed exponential coordination-decay law nb(γ)=n0/2[1+e^{−Aγ}], imported from Ref. [36]. This functional form is not derived from microscopic dynamics in the present paper, nor is nb(γ) measured or simulated. Because the yield strain is estimated as εy≈1/A and A contains the fitted structural relaxation time τc, the predicted overshoot position and height are largely controlled by this assumed form and by the fitted τc. A direct test of nb(γ) under strain, a microscopic derivation of this law, or a clear statement that the overshoot mechanism is an assumption rather than an established prediction is needed to support the central claim.","section":"§2, Eqs. (1), (11), (12)"},{"comment":"The abstract and §3 state that the stretched-exponential exponent β is obtained independently from relaxation measurements and provides the primary dynamical input. Fig. 1 documents this only for the metallic glass. For PC, PMMA, and epoxy, β=0.85 (and 0.80 for one PC rate) is used without showing corresponding relaxation data or citing independent measurements. Please provide the relaxation fits for the polymers or explicitly identify β as a fit parameter; otherwise the claim of independent dynamical input is overstated.","section":"§3, Table 1 and Fig. 1"}],"minor_comments":[{"comment":"The matrix for the strain tensor appears garbled: the entries read 'ε0 0; 0−νε0; 0 0−νε'. The intended structure (ε on the diagonal, with lateral contraction −νε) should be typeset clearly.","section":"§2, Eq. (4)"},{"comment":"The absolute-value notation ε≡|ε|, σ≡|σ| is introduced after Eq. (5), but Eq. (6) and Eq. (12) use ε in exponential and prefactor without consistently indicating that these are magnitudes. Please make the convention explicit throughout.","section":"§2, Eqs. (5)–(6)"},{"comment":"The text says β changes only slightly 'at and above the glass transition temperature.' Since Tg=625 K and the two temperatures are 643 K and 623 K, one measurement is above and one is below Tg. Please rephrase.","section":"§3, Fig. 1 and text"},{"comment":"η is described both as a fitting parameter and as 'extracted from the linear viscous-flow region of the stress–strain curves.' Clarify whether η is fitted or directly read off the data, since the distinction affects the count of free parameters.","section":"Table 1 caption"},{"comment":"The parameter A is introduced twice: first in the text as A=Δ/kBT+1/(γdot τc) and then in Eq. (11) in uniaxial form. It would help to define a single symbol and consistently state its relation to γdot versus εdot.","section":"§2, Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The paper's central mechanistic claim is plausible but not yet fully supported: the overshoot mechanism rests on an imported exponential coordination-decay law that is neither derived nor directly tested, and the validation is largely curve fitting with several adjustable parameters per material and strain rate. The viscous-stress cusp for β<1 is a separate technical inconsistency that should be fixed. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper offers a single constitutive equation — nonaffine elastic term plus stretched-exponential viscous term plus Arruda-Boyce hardening — and shows it can reproduce stress-strain curves for a metallic glass, PC, PMMA, and epoxy across a wide range of strain rates. That combination, applied to uniaxial compression of these materials, is new in this form. The limiting behavior is correct (linear elasticity at small strain, Newtonian flow at large strain), the parameter values in Table 1 are physically plausible, and the stretching exponent for the metallic glass is taken from independent stress-relaxation data rather than fit to the stress-strain curves. The polymer hardening term is standard but handled sensibly with an invariant stretch measure. Credit where due: the synthesis is genuinely useful and clearly presented.\n\nThe soft spots are real but not fatal. The validation is dominated by curve fitting. Each material and strain rate has its own E, η, τc, τv, and for polymers ER and β; τc directly sets the overshoot position through A, and the predicted yield strain is essentially 1/A. More importantly, the overshoot peak comes from the assumed exponential form of the coordination decay, nb(γ) = (n0/2)(1+e^{−Aγ}), imported from Ref. [36]. That functional form is not derived here and no direct measurement of nb(γ) is offered. If the real connectivity-loss kinetics is sigmoidal or cooperative — as one might expect from STZ-like or strain-localization mechanisms — the predicted overshoot position and height would change. So the central mechanistic claim, that overshoot is a nonaffine elastic instability competing with viscous relaxation, is not independently established; it is contingent on an ansatz. For the polymers, β is set to 0.85 without independent support, and no baseline models are compared, so we cannot tell how much the fit quality comes from the physics versus from having enough adjustable pieces. These are limitations a careful referee should probe, not reasons to dismiss the paper.\n\nWho this is for: people working on constitutive modeling of amorphous solids, especially those who want a compact analytical expression that captures elastic loading, yield, flow, and hardening in one place. It deserves a serious referee. The referee should ask for a derivation or direct test of the coordination-decay law, a comparison against a simpler phenomenological model with the same number of parameters, and at least one genuinely predictive check where parameters are fixed by one dataset and used to predict another. I would not block publication over the curve fitting, but the 'microscopic' language should be softened until the load-bearing ansatz is actually tested.","headline":"A readable, useful three-term constitutive law that tracks a lot of data, but its yield peak is essentially installed by an imported exponential coordination-decay ansatz, so the 'microscopic' story is not as supported as the fits suggest.","tokens_in":11163,"tokens_out":1549,"would_cite":true,"duration_ms":19266,"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":"This paper claims that the full stress–strain curve of an amorphous solid—linear loading, stress overshoot, yielding, flow, and strain hardening—follows from a single constitutive equation built on nonaffine elasticity and structural relaxa","keywords":["amorphous solids","stress overshoot","yielding","strain hardening","nonaffine elasticity","constitutive model","metallic glasses","polymer glasses"],"falsifier":"Measure the yield strain εy as a function of temperature and strain rate for one material over a wide range; if εy does not track 1/A = 1/(Tg/T + 1/(γ̇τc)) with a single τc, or if the shape of the overshoot deviates systematically from the predicted exponential form, the central mechanism is falsified. A cleaner check: at fixed strain rate, the theory predicts a specific shift of the overshoot with temperature; a dataset where yield strain decreases as temperature rises (as A decreases) would support it, while a non-monotonic dependence would contradict it.","tokens_in":1400,"feed_emoji":"⚙️","tokens_out":1436,"duration_ms":43409,"temperature":0.7,"pith_summary":"The paper tries to establish that stress overshoot, yielding, and strain hardening in amorphous solids are not separate phenomena needing separate models: all follow from the competition between strain-driven loss of mechanical connectivity and viscous relaxation. The authors write the elastic stress as a function of a strain-dependent coordination number that decays exponentially with strain; this term alone creates an elastic instability whose peak marks yielding. A stretched-exponential viscous stress, with its exponent taken independently from relaxation measurements, supplies the dissipative plateau, and an inverse-Langevin term accounts for polymer strain hardening. The same three-term constitutive law is shown to reproduce experimental stress–strain curves for a metallic glass, polycarbonate, PMMA, and epoxy across strain rates spanning several decades. If correct, this unifies relaxation dynamics and mechanical yielding in amorphous materials.","feed_headline":"Yield point in glasses traced to loss of atomic connectivity","feed_subtitle":"A three-term constitutive law reproduces stress-strain curves for metallic and polymer glasses across many strain rates.","key_machinery":"The central object is the strain-dependent coordination number nb(γ)=n0/2(1+e^{-Aγ}), imported from earlier work, which describes how mechanical connectivity is lost as strain grows. It enters the elastic free energy as Fel = (1/2)K[nb(γ)-n_c]γ², where n_c is the isostatic (marginal-stability) coordination, so the elastic stress acquires a maximum when connectivity loss outruns strain. A generalized Maxwell viscous stress with stretched-exponential relaxation, σ_visc=ηε̇(1-e^{-(ε/(ε̇τ_v))^β}), supplies the dissipative branch, and an inverse-Langevin chain-stretch term adds strain hardening. Together these three terms form the constitutive law whose yield strain is set by A=Tg/T+√3/(2(1+ν)ε̇τ","core_discovery":"At its core, the paper claims that stress overshoot is not primarily a plastic or damage event but an elastic instability: as strain grows, the mean number of mechanically active neighbors per atom falls, the elastic modulus softens, and the stress–strain curve reaches a maximum at a yield strain of order 1/A, where A combines Tg/T with a strain-rate term 1/(γ̇τc). The same instability, opposed by a growing viscous stress, produces the overshoot-and-drop shape seen in experiments. For polymers, finite chain extensibility adds hardening at large strain. The authors show fits to data for metallic glass, PC, PMMA, and epoxy using mostly independently estimated parameters, with τc as the main ad","pith_inferences":["A direct extension would test the model on start-up shear of colloidal glasses, where strain rate and packing fraction are independently controllable, to see whether εy≈1/A holds without adjustable τc.","If the exponential connectivity-loss law is replaced by cooperative dynamics, such as shear-transformation-zone activation, the predicted overshoot position and height would change; comparing the two functional forms against low-temperature or high-rate data could discriminate between them.","The framework implies that apparent yield stress is not a material constant but a rate- and temperature-dependent outcome of the elastic-viscous competition, which could guide design of tougher amorphous polymers by tuning τc or the connectivity-decay rate.","The success of the model hinges on the borrowed exponential form of nb(γ); a first-principles derivation of that law from cage-breaking kinetics would be the natural next step."],"forward_implications":["If overshoot is an elastic instability from connectivity loss, the yield strain should scale as εy≈1/A across temperatures and strain rates; the fits provide a direct test of this prediction.","The same timescale τc and stretching exponent β extracted from linear relaxation measurements feed the nonlinear response without additional dynamical input.","Metallic glasses need only the elastic plus viscous terms, while hardening is a polymer-specific chain-extensibility effect.","The model predicts a viscosity that decreases as a power law with strain rate (shear thinning), with distinct exponents for metallic versus polymeric glasses.","Because the elastic term vanishes at large strain, the theory naturally recovers steady viscous flow in the large-strain limit."],"fun_headline_variants":["Glassy yield traced to connectivity loss, not damage","Stress overshoot in glasses from elastic instability","Atomic connectivity loss drives glass yielding","Unified theory predicts full stress-strain of glasses","Why glasses overshoot: elastic instability wins over plasticity"],"cache_read_input_tokens":12416,"weakest_assumption_plain":"The load-bearing premise is the assumed exponential decay of mechanical coordination with strain, nb(γ)=n0/2(1+e^{-Aγ}); the yield-strain value εy≈1/A and the overshoot height follow directly from this functional form, which is imported from earlier work rather than derived from the relaxation dynamics in this paper.","fun_headline_variants_meta":{"raw":{"variants":["Glassy yield traced to connectivity loss, not damage","Stress overshoot in glasses from elastic instability","Atomic connectivity loss drives glass yielding","Unified theory predicts full stress-strain of glasses","Why glasses overshoot: elastic instability wins over plasticity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1161,"prompt_tokens":709,"completion_tokens":452,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":453,"completion_tokens_details":{"reasoning_tokens":382}},"tokens_in":453,"tokens_out":452,"duration_ms":5462,"temperature":1.0,"reasoning_tokens":382,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T03:52:56.580962+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the yield strain εy as a function of temperature and strain rate for one material over a wide range; if εy does not track 1/A = 1/(Tg/T + 1/(γ̇τc)) with a single τc, or if the shape of the overshoot deviates systematically from the predicted exponential form, the central mechanism is falsified. A cleaner check: at fixed strain rate, the theory predicts a specific shift of the overshoot with temperature; a dataset where yield strain decreases as temperature rises (as A decreases) would support it, while a non-monotonic dependence would contradict it.","supporting_citations":[],"review_version":1}