{"id":"39c57a7d-75b7-47ff-b577-032101100c5f","arxiv_id":"1908.06722","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Sheared glassy Kob-Andersen mixtures show nearly identical reduced-unit dynamics along isomorphs, including the exponential stress-drop tails that signal avalanches.","lead":"This paper tests whether a theory called isomorph invariance, which predicts that certain density and temperature pairs behave identically, holds for glassy materials being sheared. The authors show that simulated glasses at paired densities and temperatures produce nearly identical stress fluctuations, particle motion, and avalanche-like stress drops.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The glass isomorphs are generated from NVT fluctuations that assume negligible aging; for the T=0.55 isomorph this is untested and potentially violated, so the compared state points may not lie on a true isomorph.","rationale":"The reader's weakest-assumption analysis identifies the same concern: the isomorphs are generated by applying equilibrium fluctuation theory inside a glass, where aging is neglected. This is genuinely load-bearing because every subsequent comparison depends on those curves being true isomorphs. The paper is transparent about the assumption and flags it as needing future work, which supports a conditional rather than a reject verdict. The observed RDF and dynamics collapse is nontrivial evidence that the assumption may be approximately valid, but it does not by itself rule out a protocol-dependent gamma. The systematic trend in epsilon_c is a warning sign that the isomorphs may be slightly misidentified, and the paper explicitly entertains this possibility. A direct test of gamma stationarity, followed by regeneration of the isomorph with a corrected gamma, would settle whether the central claim is robust. The reader's CONDITIONAL verdict remains appropriate, so no change is needed.","tokens_in":17259,"tokens_out":6410,"duration_ms":76251,"concrete_test":"At the first point of the high-T isomorph (rho=1.265, T=0.55), run NVT for 10^6, 10^7, and 10^8 steps from the same starting configuration, and compute gamma = <ΔUΔW>/<ΔU^2> in successive time windows with block errors. If gamma changes by more than the block error, Eq. (5) is not well defined. Then regenerate the isomorph using the late-time gamma and repeat the stress-change histogram collapse and epsilon_c trend; if collapse degrades or epsilon_c drift disappears, the original invariance claim was contingent on the aging-biased isomorph.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV integrates Eq. (4) by Euler stepping with gamma from NVT fluctuations of U and W in the glass (Eq. (5)), explicitly assuming 'no significant aging occurs during the simulation.' This is the sole construction of the isomorphs. For the high-temperature isomorph, the starting point T=0.55 corresponds to about 0.42 at rho=1.2 with a relaxation time of 3820 LJ units, while the NVT runs are 10^7 steps (approximately 25,000 LJ units), so the system can partially relax during the measurement. If gamma drifts with waiting time, the integrated path is not a well-defined configurational adiabat, and the 'isomorph' is protocol dependent. Section VII C(1) acknowledges that this assumption 'needs to be developed and evaluated theoretically,' and the Conclusion concedes the assumption is 'less clear for the high temperature one.' The systematic decrease of epsilon_c with density in Fig. 6 is exactly the kind of residual trend expected from an off-isomorph path; the paper itself lists an incorrect isomorph as one possible explanation. Because the central claim is invariance along these particular curves, an untested drift in gamma is the load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper tests whether isomorph theory, developed for equilibrium liquids, applies to a sheared Kob-Andersen binary Lennard-Jones glass below the glass transition. The authors generate two glassy isomorphs by numerically integrating the configurational-adiabat equation, Eq. (4), using the potential-energy/virial fluctuation slope gamma from NVT simulations, starting from configurations cooled at constant pressure to T=0.55 and T=0.10. They then shear the glasses with SLLOD and Lees-Edwards boundary conditions at fixed reduced strain rates and compare, in reduced units, the radial distribution function, steady-state flow stress, stress fluctuations and their autocorrelation, stress-change histograms, transverse intermediate scattering functions, and mean-squared displacements. They report good collapse of most observables along each isomorph, including an exponential negative tail in the stress-change distribution that is interpreted as an avalanche signature and is claimed to be isomorph invariant. The paper also discusses how isomorph invariance constrains expressions for flow stress and proposes an alternative zero-temperature compatible reduced-unit system based on a density-scaling function h(rho).","tokens_in":17517,"tokens_out":4125,"duration_ms":45501,"significance":"If the main claim holds, the paper provides a nontrivial extension of isomorph theory into the non-equilibrium glassy regime, with practical consequences: structure and sheared dynamics, including avalanche statistics, would depend on density and temperature only through the isomorph and reduced strain rate, not separately. The test is genuinely predictive: the isomorphs are constructed from unsheared NVT fluctuations, while the tested observables come from independent sheared simulations, and the main collapse is not a fit. The paper is also refreshingly explicit about its limitations, flagging the aging assumption and the imperfect collapse of the stress autocorrelation. However, the evidence for the central claim is weakened by an unverified assumption in the construction of the high-temperature isomorph, a systematic unexplained drift in one key fit parameter, and the absence of any quantitative measure of collapse. With those points addressed, the result would be a solid contribution to the glassy-dynamics and isomorph-theory literature.","major_comments":[{"comment":"","section":"Sec. IV, Eq. (5); Sec. VII C(1); Conclusion"},{"comment":"","section":"Sec. V, Fig. 6"},{"comment":"","section":"Secs. V and VI, Figs. 4, 5, 7, 8, 11, 12"}],"minor_comments":[{"comment":"","section":"Sec. II A"},{"comment":"","section":"Sec. VII A, Eq. (9)"},{"comment":"","section":"Sec. III, Fig. 1 and Sec. IV"},{"comment":"","section":"Sec. V, Fig. 4 caption"},{"comment":"","section":"References"}],"recommendation":"major_revision","confidential_remarks":"I agree with the reader's assessment that the aging assumption in the construction of the high-temperature isomorph is the main correctness risk. The paper is honest about this risk, and the central idea is sound and testable, so I would not reject. However, the requested checks (convergence of gamma with run length, quantitative collapse metrics, and resolution of the epsilon_c drift) are necessary before the claim of isomorph invariance of sheared glassy dynamics is fully supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this is the first test I know of for isomorph invariance of finite-temperature deformation of a glass, and the core empirical result looks solid. For the Kob-Andersen system below Tg, state points linked by the standard U/W-fluctuation isomorph construction collapse well in reduced units for RDFs, steady-state flow stress, stress-change histograms, transverse MSDs, and ISFs. The genuinely new bit is the avalanche signature: the exponential tail on the negative side of the stress-change distribution is isomorph invariant. That is a useful constraint for plasticity models. The paper is also honest in the right places — it flags the aging assumption, the eps_c drift, and the provisional strain-scale hierarchy.\n\nThe soft spot the stress-test identifies is real. Isomorphs are generated from NVT fluctuations assuming 'no significant aging' during the run. For the T=0.55 isomorph, the reference relaxation time is about 3820 LJ units and the NVT runs are about 25,000 LJ units, so partial relaxation is possible. If gamma drifts with waiting time, the integrated path may be protocol-dependent, and the compared points are not strictly the equilibrium-theory isomorphs. That does not sink the demonstration — the collapse still holds along the generated curves — but it does mean the paper's central interpretation is not fully nailed. The authors admit this in Sec. VII C(1) and the conclusion. The linear decrease of eps_c with density in Fig. 6 is a residual trend with three candidate explanations (statistics, off-isomorph path, or genuine limit of invariance), and they don't discriminate among them. That is a real gap, not a fatal one.\n\nTwo smaller concerns. Collapse quality is assessed visually; a quantitative metric would tighten the claim, especially since 'good collapse' is the paper's main evidence. And no data or code are released, though RUMD itself is open source. Neither is disqualifying.\n\nThe audience is specialists in isomorph theory and in shear rheology/plasticity of glasses. The Lerner–Procaccia correspondence and the alternative reduced-unit discussion are useful beyond the simulations. I would send this to peer review and ask for quantitative collapse metrics and a sharper treatment of the aging assumption. Worth citing, worth discussing in the group.","headline":"First finite-temperature sheared-glass isomorph test: real empirical result, honest paper, but the aging assumption in generating the glass isomorphs is the load-bearing weakness.","tokens_in":18018,"tokens_out":2420,"would_cite":true,"duration_ms":26127,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that isomorphs generated from potential-energy and virial fluctuations remain valid inside the glassy state: state points along a glassy isomorph have identical reduced-unit structure and steady-state shearing dynamics…","keywords":["isomorph theory","hidden scale invariance","glassy dynamics","shear deformation","avalanche statistics","stress fluctuations","binary Lennard-Jones glass","reduced units"],"falsifier":"Generate the same glassy isomorph two ways — from the NVT potential-energy/virial fluctuations used here and from fluctuations measured during steady-state shearing at the same reduced strain rate — and compare the density–temperature pairs; systematic disagreement would show the equilibrium-based isomorphs are history-dependent. As a second check, run the shearing at each state point for several times more total strain and see whether the characteristic stress-decay strain continues to drift monotonically with density; a persistent drift beyond statistical error would mark a genuine limit of isomorph invariance.","tokens_in":17081,"feed_emoji":"🧊","tokens_out":7966,"duration_ms":78166,"temperature":0.7,"pith_summary":"The paper argues that isomorph theory — the idea that certain density–temperature pairs have identical structure and dynamics once lengths, times, and energies are expressed in reduced units — continues to hold inside the glassy state, not just in equilibrium liquids. Using the Kob-Andersen binary Lennard-Jones model below its glass transition, the authors generate two glassy isomorphs from potential-energy and virial fluctuations, then shear the glasses with the SLLOD algorithm and Lees-Edwards boundary conditions. They find that the steady-state flow stress, stress fluctuations, stress-change histograms, transverse mean-squared displacement, and intermediate scattering function all collapse along the isomorphs when the reduced strain rate is held fixed. The clearest signature of avalanches, an exponential tail on the negative side of the stress-change distribution at low temperature and low strain rate, is also isomorph invariant. If correct, this means glassy plasticity depends on density and temperature only through the isomorph coordinate, simplifying the phase diagram for out-of-equilibrium amorphous solids.","feed_headline":"Isomorphs hold below the glass transition, even for avalanches","feed_subtitle":"Density–temperature pairs on the same isomorph share flow stress, stress statistics, and the exponential tail marking avalanches.","key_machinery":"The central object is the isomorph: a curve in the density–temperature plane along which reduced-unit structure and dynamics are invariant. Isomorphs are identified through hidden scale invariance, quantified by the correlation coefficient R between potential energy U and virial W fluctuations; the slope gamma = <$\\Delta$ U $\\Delta$ W> / <($\\Delta$ U)^2> gives the configurational adiabat d ln T / d ln rho = gamma, integrated numerically to step between state points. Dynamics are probed by Couette shear via the SLLOD algorithm with Lees-Edwards boundary conditions, with the reduced strain rate gamma_dot_tilde = gamma_dot (T/m)^-1/2 $rho^{-1}$/3 held fixed along each isomorph. The avalanche analysis rests on histograms of reduced stress changes $\\Delta$ $\\sigma$ / (rho k_B T) over strain intervals, whose negative exponential tail marks correlated plastic events.","core_discovery":"On the paper's own terms, the central discovery is that isomorphs generated in the glassy phase from equilibrium fluctuation formulas carry the dynamics of steady-state shearing: state points connected by the integration of d ln T / d ln rho = gamma, with gamma from potential-energy–virial correlations, have the same reduced-unit structure and the same statistical dynamics under shear, including the avalanche signature. The authors demonstrate this for two glassy isomorphs, one just below the glass transition and one deep in the glass, at nominal strain rates from $10^{-2}$ down to $10^{-5}$. The collapse holds for the mean flow stress, its standard deviation, the distribution of stress changes over strain intervals, and the transverse particle dynamics; the stress autocorrelation shows the poorest collapse, with a systematic decrease of the characteristic decay strain along the isomorph. The paper also identifies a hierarchy of strain scales — thermal/vibrational, avalanche, maximum-skewness, and correlation-decay — and shows that this structure is invariant along the isomorph.","pith_inferences":["A testable extension is to apply the same protocol to glass formers with weaker potential-energy/virial correlations: the quality of the stress-drop collapse would then serve as a quantitative predictor of how much density dependence a plasticity theory must contain.","The hierarchy of strain scales the paper identifies — thermal/vibrational, avalanche, maximum-skewness, and correlation-decay — suggests that avalanche initiation and stress relaxation are governed by separate strain intervals roughly an order of magnitude apart, a separation that mesoscale plasticity models could target.","Because the reduced flow stress differs by almost a factor of ten between the two isomorphs, the choice of reduced units is not neutral for comparing different isomorphs; the paper's alternative h(rho)-based units may be the more physical ones near zero temperature, a point worth testing on the same data."],"forward_implications":["Steady-state shearing of a glass is controlled by two variables, the isomorph coordinate and the reduced strain rate, rather than by density and temperature separately.","Flow stress and its fluctuations are isomorph invariant in reduced units, so theories of glassy plasticity that treat temperature and density separately are incomplete.","Avalanche statistics, read from the negative exponential tail of stress-change distributions, are unchanged along an isomorph at fixed reduced strain rate.","Density dependence in existing flow-stress formulas can be absorbed into isomorph-invariant forms using the density-scaling function, and an alternative reduced-unit choice gives a well-behaved zero-temperature limit.","Transverse particle diffusion is largely strain-rate controlled at low temperature, with thermal activation entering only on high-temperature isomorphs."],"supporting_citations":[{"why":"Defines isomorphs and the configurational-adiabat relation d ln T / d ln rho = gamma used to generate them.","marker":"[5]"},{"why":"States hidden scale invariance as preservation of potential-energy order, the basis for extracting gamma from U-W fluctuations.","marker":"[23]"},{"why":"Shows isomorph invariance of steady-state shearing dynamics in the non-viscous fluid regime, the protocol this paper extends below the glass transition.","marker":"[11]"},{"why":"Establishes that the equations of motion are identical in reduced units, justifying reduced-unit comparison of dynamics.","marker":"[4]"},{"why":"Introduces the binary Lennard-Jones glass former that all simulations use.","marker":"[18]"},{"why":"Supplies the strain-rate criterion that places the runs in the glassy deformation regime, and the flow-stress expression later rewritten in isomorph-invariant form.","marker":"[30]"},{"why":"Extends isomorph theory to athermal shearing of glasses, the zero-temperature precedent for this finite-temperature study.","marker":"[14]"},{"why":"Provides the analytical isomorph construction for the Lennard-Jones reference liquid and the density-scaling function h(rho) used in the discussion of alternative reduced units.","marker":"[9]"}],"fun_headline_variants":["Isomorphs survive shearing in glasses, avalanches included","Shear dynamics collapse along glass isomorphs, avalanches too","Glassy isomorphs unify shear stress stats and avalanches","Isomorph invariance holds for sheared glasses, even avalanches","Isomorphs predict avalanche dynamics in sheared glasses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that isomorphs can be generated inside the glass by ordinary NVT fluctuation formulas: the paper assumes aging is negligible during those runs, a point it explicitly flags in the section on future improvements, and if aging biases the potential-energy and virial fluctuations the resulting density–temperature pairs are not true isomorphs and the observed collapse could be approximate or coincidental.","fun_headline_variants_meta":{"raw":{"variants":["Isomorphs survive shearing in glasses, avalanches included","Shear dynamics collapse along glass isomorphs, avalanches too","Glassy isomorphs unify shear stress stats and avalanches","Isomorph invariance holds for sheared glasses, even avalanches","Isomorphs predict avalanche dynamics in sheared glasses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00057,"raw_usage":{"total_tokens":2684,"prompt_tokens":918,"completion_tokens":1766,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":1680}},"tokens_in":534,"tokens_out":1766,"duration_ms":13150,"temperature":1.0,"reasoning_tokens":1680,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:35:54.914666+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Generate the same glassy isomorph two ways — from the NVT potential-energy/virial fluctuations used here and from fluctuations measured during steady-state shearing at the same reduced strain rate — and compare the density–temperature pairs; systematic disagreement would show the equilibrium-based isomorphs are history-dependent. As a second check, run the shearing at each state point for several times more total strain and see whether the characteristic stress-decay strain continues to drift monotonically with density; a persistent drift beyond statistical error would mark a genuine limit of isomorph invariance.","supporting_citations":[{"cited_title":"high temperature","cited_arxiv_id":null,"evidence_quote":"Defines isomorphs and the configurational-adiabat relation d ln T / d ln rho = gamma used to generate them."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows isomorph invariance of steady-state shearing dynamics in the non-viscous fluid regime, the protocol this paper extends below the glass transition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that the equations of motion are identical in reduced units, justifying reduced-unit comparison of dynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the binary Lennard-Jones glass former that all simulations use."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the strain-rate criterion that places the runs in the glassy deformation regime, and the flow-stress expression later rewritten in isomorph-invariant form."},{"cited_title":"Bøhling, T","cited_arxiv_id":null,"evidence_quote":"Extends isomorph theory to athermal shearing of glasses, the zero-temperature precedent for this finite-temperature study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the analytical isomorph construction for the Lennard-Jones reference liquid and the density-scaling function h(rho) used in the discussion of alternative reduced units."}],"review_version":1}