Pith. sign in

REVIEW 3 major objections 4 minor 52 references

Predicting temperature-dependent failure and transformation zones in 2D silica glass through quasistatic Gaussian Phase Packets

T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper claims that atoms can be represented as temperature-dependent Gaussian packets, extending quasistatic athermal simulation to finite temperature and predicting fracture sites in 2D silica glass from per-atom covariance.

desk verdict Anisotropic GPP covariance is a solid failure-site predictor; the Metropolis-GPP onset 'validation' is a two-parameter fit dressed as a rate dependence. read the letter →

arxiv 2507.13960 v1 pith:USPXZPBZ submitted 2025-07-18 cond-mat.dis-nn cond-mat.mtrl-sci

Miguel Spínola , Shashank Saxena , Franz Bamer , Dennis M. Kochmann This is my paper
classification cond-mat.dis-nncond-mat.mtrl-sci
keywords atomisticsmultiscalemodelingdisorderedsolidstress-straincurvethermalexpansionsilicaglassGaussianPhasePacketsMetropolissampling
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper is trying to establish that the athermal quasistatic method, which is limited to near-zero temperature, can be extended to finite temperature by describing each atom as a Gaussian packet in phase space and relaxing the packet means and covariances to a free-energy minimum. On 2D silica glass under uniaxial tension, it claims this Gaussian-packet picture captures thermal expansion, predicts temperature-dependent onset of fracture when combined with Metropolis sampling, and identifies the sites where bonds will break from the per-atom covariance alone. The payoff, if true, is a way to approach quasi-static loading at finite temperature without the cost of long molecular-dynamics runs, while also getting a cheap map of where rearrangement zones will form.

What carries the argument

The central object is the anisotropic Gaussian phase packet: a per-atom Gaussian probability density in phase space parametrized by a mean position $\bar{\boldsymbol{q}}_i$ and a $2\times 2$ position-covariance matrix $\boldsymbol{\Sigma}_i^{(q,q)}$, with the momentum covariance fixed by temperature through the equipartition relation $\boldsymbol{\Omega}_i = m_i k_B T$. In the quasistatic limit, the paper solves the mean-force balance $\langle \boldsymbol{f}_i\rangle = \mathbf{0}$ and the thermal-force balance (Eq. 21b) for each atom, which together define the free-energy-minimum ensemble at constant temperature. The scalar $\Sigma_i = (\det \boldsymbol{\Sigma}_i^{(q,q)})^{1/2}$ is the failure indicator, and it is available at every relaxed state at negligible extra cost. The Metropolis-GPP extension uses the relaxed GPP density as the proposal distribution and accepts or rejects trial microstates with probability $\min[1, \exp(-\Delta U/k_B T)]$, giving new starting guesses that can fall into adjacent free-energy minima.

What would settle it

Fix a target ratio $\gamma/N_p = 2\times 10^{-8}$ in the Metropolis-GPP scheme at $T^*=0.00025$ and run two sets of simulations, $\{\gamma,N_p\}=\{10^{-4},5000\}$ and $\{\gamma,N_p\}=\{10^{-3},50000\}$, on the same initial sample. If the first-inelastic-event stress and site differ by more than the observed scatter, then $N_p$ is not a physical clock and the claimed agreement with $10^4\,\mathrm{s}^{-1}$ MD is calibration rather than validation.

Watch

Extended reading notes

Core claim

The central claim is that a quasistatic ensemble of Gaussian phase packets, one per atom, with an anisotropic position covariance matrix, provides a finite-temperature extension of athermal quasistatic simulation: at each deformation step the mean positions and covariances are relaxed to a free-energy minimum at constant temperature. Applied to three 2D silica samples, the anisotropic form reproduces MD thermal expansion to within about 0.2% in box dimensions and, more importantly, the scalar $\Sigma_i = (\det \boldsymbol{\Sigma}_i^{(q,q)})^{1/2}$ marks the two oxygen sites where MD later breaks bonds, both in the undeformed state and with increasing sharpness as deformation proceeds. When the relaxed GPP density is used as the proposal distribution in a Metropolis algorithm, the resulting prototype statistics match MD for the stress and strain at the first inelastic event at temperatures $T^*=0.00025$ and $T^*=0.00075$, including the lowest MD strain rate of about $10^4\,\mathrm{s}^{-1}$. The paper interprets the matching cases as showing that $\gamma/N_p$ can be about two orders of magnitude larger than $\gamma/N_e$, making the scheme a candidate for near-quasistatic finite-temperature fracture predictions at a fraction of MD cost.

Load-bearing premise

The load-bearing premise is that a Metropolis sampling pass can be treated as taking approximately the same physical time as one molecular-dynamics timestep, so $\gamma/N_p$ can stand in for the MD strain rate $\gamma/N_e$; the paper itself notes the two are not equivalent.

Editorial extensions

If this is right

  • A finite-temperature quasistatic loading protocol for disordered solids can be run without MD timestep integration, with temperature entering through free-energy minimization instead of thermal velocities.
  • The anisotropic per-atom covariance, obtained for free at each relaxed state, provides a pre-failure map of likely fracture or rearrangement sites, and the map sharpens as deformation progresses.
  • Metropolis-GPP reproduces the temperature dependence of the first-inelastic-event stress seen in MD at about $10^4\,\mathrm{s}^{-1}$, while using $\gamma/N_p$ ratios about two orders of magnitude larger than $\gamma/N_e$.
  • Thermal expansion of the disordered network solid emerges from the same zero-stress free-energy minimization, matching MD box dimensions to within about 0.2% at the highest temperature studied.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • We infer that the covariance indicator, because it is a by-product of every relaxed state, could serve as a practical soft-spot predictor for 3D glasses and sheared metallic glasses without computing dynamical-matrix eigenmodes; the paper only demonstrates this for 2D silica.
  • We infer that the $\gamma/N_p$ equivalence is testable by varying $N_p$ while holding $\gamma/N_p$ fixed; the paper does not report such a test, so the low-strain-rate validation claim is only as strong as that untested identification.
  • We infer that the relaxed mean positions and covariances could be post-processed to estimate local activation barriers for the highest-$\Sigma$ bonds, turning the ranked failure sites into rate estimates; the paper currently stops at ranking.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper extends the athermal quasistatic (AQS) method to finite temperature by representing atoms as Gaussian phase packets (GPPs) whose mean positions and position covariances are obtained from a free-energy minimization. The method is applied to 2D silica glass under uniaxial tension. Two main results are claimed: (i) the anisotropic position covariance provides an inexpensive predictor of the sites of atomic rearrangements and incipient bond breaking, and (ii) when the GPP relaxation is augmented by a Metropolis sampling pass, the stress and strain at the onset of inelasticity match MD results at effective strain rates as low as 10^4 s^-1, including the temperature dependence between about 160 K and 485 K. A thermal-expansion benchmark against MD is also presented, with errors below 0.2%.

Significance. If the Metropolis-GPP validation were rigorous, the paper would offer a genuinely useful route toward finite-temperature quasistatic simulations of disordered solids at much lower effective strain rates than direct MD. The thermal-expansion benchmark is clean and convincing, and the observation that the anisotropic GPP covariance localizes the experimentally/MD-observed initial failure sites, both in the undeformed state and progressively during deformation, is a valuable and interesting result that appears well supported by the presented snapshots. The efficiency argument, however, is entirely dependent on an unvalidated equivalence between Metropolis passes and MD time steps, and the paper's own wall-clock comparison shows no actual speedup in the current implementation. The central validation claim is therefore weaker than the abstract suggests.

major comments (3)
  1. [Section 2.3, Section 3.2.3, Fig. 10, abstract] The claim that Metropolis-GPP 'predicts the effect of temperature on the onset of fracture' and is 'validated through MD simulations at strain rates as low as 10^4 s^-1' is not established, because the mapping between Metropolis passes and physical time is assumed rather than derived. The paper defines N_s = N_p N trial moves per deformation increment (Section 2.3) and then compares the ratio gamma/N_p directly with the MD dimensionless strain rate gamma/(L_y N_e dt) (Appendix A, Eq. 35) in Figure 10. The only justification offered is the statement, citing reference [1], that 'in some scenarios a sampling pass can be considered to be proportional to an MD step.' No calibration of this proportionality for the specific Yukawa silica potential is provided, and the constant of proportionality is not fixed. Because both gamma and N_p are free parameters, reproducing the MD onset stresses at two temperatures is not a validation: two parameters can fit two data points, and the agreement does not demonstrate that gamma/N_p controls the onset in a predictive sense. A rigorous test would fix gamma and N_p from one condition and predict the onset at other temperatures and strain rates, or would derive the pass-to-time mapping from a physical argument.
  2. [Section 3.2.3 and Conclusion] The 'two orders of magnitude' efficiency gain is based on the same unvalidated gamma/N ratio comparison and is not a statement about actual computational cost. The paper itself states that the GPP procedure 'did not result in important gains' in wall-clock time against MD in the presented implementation, and that the wall-clock comparison was performed on one core versus eight cores. The abstract's characterization of the framework as predicting onset 'efficiently without the need for expensive MD simulations' is therefore overstated relative to the evidence. The efficiency case would be strengthened by a controlled wall-clock comparison on the same hardware with a realistic accounting of the Metropolis sampling and GPP relaxation costs, or by a clear statement that only asymptotic complexity and not realized performance is being claimed.
  3. [Section 3.2.3] The matching between Metropolis-GPP and MD onset is only shown for two temperatures and a single sample (sample (a) in Figure 1). The paper states that 'the first inelastic onset in the different MD runs ... are close to those predicted by the Metropolis-GPP procedure' for particular parameter pairs, but the scatter among the 25 runs and the differences among the parameter pairs are not quantified. Given the underdetermination described above, a more extensive test across more temperatures, more samples, and a fixed parameter choice is necessary to support the general claim of temperature-dependent onset prediction.
minor comments (4)
  1. [Section 2.2 (after Eq. 21)] The word 'anisotorpic' appears as a typo; it should be 'anisotropic.'
  2. [Section 3.2.3, Figure 10] The top and bottom axes of Figure 10 compare gamma/N_p with the MD dimensionless strain rate, but the axes are not labeled with the assumed proportionality; adding an explicit statement that the comparison presumes one Metropolis pass equals one MD timestep would make the figure less misleading.
  3. [Section 3.2.2] The failure-site prediction using the anisotropic covariance is demonstrated for sample (a) only. Applying the same analysis to samples (b) and (c) would strengthen the claim that this is a general predictor rather than a selected example.
  4. [Section 3.2.3] The paper notes that Metropolis trials and MD steps are not equivalent, but this caveat appears only after the comparison has been made in Figure 10; moving this caveat earlier and stating it as a formal limitation of the method would improve clarity.

Circularity Check

1 steps flagged · score 6.0 of 10

Metropolis-GPP low-rate onset validation reduces to a post-hoc match of the free parameters {γ, N_p}; thermal expansion and failure-site predictions remain independent.

  1. fitted input called prediction [Section 3.2.3 (Metropolis-Gaussian Phase Packets), Figs. 9-10; abstract and conclusion]
    "Specifically, at both temperatures shown in Figure 9 the first inelastic onset in the different MD runs with \dot{\epsilon}^*=2.4\cdot10^{-10} (\gamma=10^{-6}, N_e=5\cdot10^3) and \dot{\epsilon}^*=1.2\cdot10^{-10} (\gamma=10^{-7}, N_e=5\cdot10^3) are close to those predicted by the Metropolis-GPP procedure using {\gamma,N_p}={10^{-4},5\cdot10^3} and {\gamma,N_p}={10^{-4},10^4}, respectively."

    The two MD onset stresses being reproduced are exactly the two temperature/rate points used for the validation, while the Metropolis-GPP parameters {\gamma,N_p} are free and are selected post hoc to produce the agreement. The paper establishes no calibrated relation between a Metropolis pass and MD time; it only states that 'in some scenarios a sampling pass can be considered to be proportional to an MD step [1]', with [1] being a generic Monte-Carlo time-scale reference not calibrated to this Yukawa silica potential. Consequently, comparing \gamma/N_p with \gamma/N_e (Fig.

full rationale

The thermal-expansion benchmark (Section 3.1) and the anisotropic-covariance failure-site predictor (Section 3.2.2) are outputs of the GPP free-energy minimization and are checked against independent MD data, so those components are not circular. The GPP equations are restated in Section 2.2, and the prior GPP citations [15,39,41] are derivational background rather than a self-citation chain that forbids alternatives. The circularity is confined to the Metropolis-GPP validation: N_p is a free sampling-count parameter with no established physical-time meaning, and the paper explicitly acknowledges that 'Metropolis trials and MD steps are not equivalent'. Choosing the GPP parameters so that selected runs land near selected low-rate MD onsets, then plotting \gamma/N_p against \gamma/N_e and calling the result 'validated at 10^4 s^{-1}', is a fitted-input-called-prediction pattern rather than a first-principles rate prediction. The paper's own efficiency caveat (that the wall-clock gain was not realized in the current implementation) further weakens the advertised advantage but is not itself circular. Overall, the central claim is partially circular because the low-rate onset agreement is constructed by parameter selection, while the independent failure-site and thermal-expansion results give the paper substantial non-circular content.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The GPP method rests on a Gaussian ansatz with independent atomic packets, equipartition closure, and ergodicity, plus a heuristic mapping from Metropolis passes to MD time. The main hand-chosen parameters are the Metropolis pass count Np and strain increment gamma, which are adjusted to reproduce MD onset stresses.

free parameters (2)
  • Number of Metropolis passes between deformation increments (Np) = 5000 to 50000 (e.g., Np = 5000, 10000, 50000)
    Controls how far the Markov chain explores before free-energy relaxation; values are chosen so that onset stresses approach low-rate MD results, with no first-principles calibration to physical time.
  • GPP strain increment gamma = 1e-3 and 1e-4 (dimensionless)
    Chosen large enough for efficiency but small enough to access intermediate free-energy minima; the paper states that the initial guess determines accessible minima, so gamma is a hand-chosen loading parameter.
assumptions (4)
  • ad hoc to paper The ensemble probability density is a product of independent atomic Gaussian packets (Eq. 13).
    Central approximation; eliminates interatomic correlations in the covariance matrix; acknowledged in Section 3.1 as a source of error.
  • domain assumption Isothermal conditions fix Omega = m_i k_B T delta via equipartition.
    Standard canonical ensemble result used to close Eqs. (17) and (21).
  • domain assumption Time averages of equilibrated MD equal phase averages of the stationary GPP density (ergodicity).
    Stated in Section 2.2; needed to compare GPP predictions with MD reference data.
  • ad hoc to paper A Metropolis pass can be treated as proportional to an MD timestep, allowing gamma/Np to be compared with gamma/Ne.
    Load-bearing for the low-strain-rate validation and efficiency claim; based on a heuristic analogy cited to [1], not derived.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Predicting temperature-dependent failure and transformation zones in 2D silica glass through quasistatic Gaussian Phase Packets." pith.science (2026). https://pith.science/paper/USPXZPBZ

@misc{pith2026250713960,
  author       = {Pith},
  title        = {Pith review of: Predicting temperature-dependent failure and transformation zones in 2D silica glass through quasistatic Gaussian Phase Packets},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/USPXZPBZ}},
  note         = {Machine review of arXiv:2507.13960}
}
abstract

The athermal quasistatic (AQS) method is a powerful technique to study the mechanical behavior of disordered systems. However, its applicability is limited to temperatures near zero, where thermal activation is unlikely. In this work, we extend the AQS method to finite temperatures, based on a formulation that describes atoms as temperature-dependent Gaussian packets (GPPs) in phase space under quasistatic conditions, thus equivalent to minimum free energy conditions. This framework is used to study the effect of temperature on the onset of inelasticity and fracture in amorphous two-dimensional silica glass approaching quasistatic conditions under uniaxial tensile loading. An important characteristic of this formulation is the directional dependence of the variance of each Gaussian packet in configuration space, making this formulation an inexpensive and accurate predictor of zones prone to atomic-scale rearrangements, both in the undeformed state and (with increasing accuracy) as the deformation progresses. This method is also shown to accurately capture the thermal expansion of the disordered material. Furthermore, combining the GPP description with Metropolis sampling predicts the effect of temperature on the onset of fracture of the material, which is validated through MD simulations at strain rates as low as $10^{4}$s$^{-1}$. The presented framework therefore provides a valuable technique for studying the nonlinear mechanics of disordered materials at finite temperature and for predicting local rearrangement zones in disordered solids efficiently without the need for expensive MD simulations.

Figures

Figures reproduced from arXiv: 2507.13960 by the authors.

Figure 1
Figure 1. shows three samples with increasing disorder. (a) (b) (c) [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Dimensions and average potential energy vs. (dimensionless) temperature for the three different 2D silica samples of [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Stress-strain curves obtained from MD simulations for different values of {𝛾, 𝑁𝑒 } for sample (a) in [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Stresses at the onset of inelasticity for the MD simulations at different strain rates. The markers correspond to those in [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: (a) was observed. For example, at 𝑇 ∗ = 0.00025 for the case {𝛾, 𝑁𝑒 } = {10−6 , 5 ⋅ 103} only in one of the 25 simulations the case of [PITH_FULL_IMAGE:figures/full_fig_p013_5.png]
Figure 6
Figure 6. Figure 6: Stress-strain curves and a magnification of the section of the onset of inelasticity (highlighted on the left) for sample (a) of [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]
Figure 7
Figure 7. Figure 7: Different stages during the uniaxial deformation response shown of the cases corresponding to 𝑇 ∗ = 0.00025 and 𝛾 = 10−3 in [PITH_FULL_IMAGE:figures/full_fig_p015_7.png]
Figure 8
Figure 8. Figure 8: Snapshots at different strains of sample (a) in [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: Stress-strain responses obtained from the Metropolis-GPP framework, using the anisotropic Gaussian approximation, compared to selected MD results (from [PITH_FULL_IMAGE:figures/full_fig_p017_9.png]
Figure 10
Figure 10. Figure 10: Stresses at the onset of inelasticity in simulations using the Metropolis-GPP framework (red) as a function of the ratio between the incremental deformation parameter 𝛾 and the number of Metropolis passes 𝑁𝑝 (top axis). For comparison, the stresses at the onset of ine…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

52 extracted references · 40 canonical work pages

  1. [1]

    On the time scale associated with monte carlo simulations

    Bal, K.M., Neyts, E.C., 2014. On the time scale associated with monte carlo simulations. The Journal of chemical physics 141

    work page 2014

  2. [2]

    Molecular Mechanics of Disordered Solids

    Bamer, F., Ebrahem, F., Markert, B., Stamm, B., 2023. Molecular Mechanics of Disordered Solids. Archives of Computational Methods in Engineering 30, 2105–2180. URL: https://link.springer.com/10.1007/s11831-022-09861-1, doi:10.1007/ s11831-022-09861-1

    doi:10.1007/s11831-022-09861-1 2023

  3. [3]

    Structural relaxation made simple

    Bitzek, E., Koskinen, P., Gähler, F., Moseler, M., Gumbsch, P., 2006. Structural relaxation made simple. Physical review letters 97, 170201

    work page 2006

  4. [4]

    Defectsinbilayersilicaandgraphene:commontrendsindiversehexagonaltwo-dimensionalsystems

    Björkman, T., Kurasch, S., Lehtinen, O., Kotakoski, J., Yazyev, O.V., Srivastava, A., Skakalova, V., Smet, J.H., Kaiser, U., Krasheninnikov, A.V.,2013. Defectsinbilayersilicaandgraphene:commontrendsindiversehexagonaltwo-dimensionalsystems. ScientificReports3,3482. URL: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3863822/, doi:10.1038/srep03482

    work page doi:10.1038/srep03482 2013

  5. [5]

    Ultrathin silica films: The atomic structure of two-dimensional crystals and glasses

    Büchner, C., Lichtenstein, L., Yu, X., Boscoboinik, J.A., Yang, B., Kaden, W.E., Heyde, M., Shaikhutdinov, S.K., Włodarczyk, R., Sierka, M., et al., 2014. Ultrathin silica films: The atomic structure of two-dimensional crystals and glasses. Chemistry–A European Journal 20, 9176–9183

    work page 2014

  6. [6]

    Molecular Dynamics, in: Comprehensive Nuclear Materials

    Cai, W., Li, J., Yip, S., 2012. Molecular Dynamics, in: Comprehensive Nuclear Materials. Elsevier, pp. 249–265. URL: https: //linkinghub.elsevier.com/retrieve/pii/B9780080560335001282, doi:10.1016/B978-0-08-056033-5.00128-2

    work page doi:10.1016/b978-0-08-056033-5.00128-2 2012

  7. [7]

    A general construction for parallelizing metropolis-hastings algorithms

    Calderhead, B., 2014. A general construction for parallelizing metropolis-hastings algorithms. Proceedings of the National Academy of Sciences111,17408–17413. URL: https://www.pnas.org/doi/abs/10.1073/pnas.1408184111,doi: 10.1073/pnas.1408184111, arXiv:https://www.pnas.org/doi/pdf/10.1073/pnas.1408184111

    work page doi:10.1073/pnas.1408184111 2014

  8. [8]

    Molecular dynamics simulations of the structure and mechanical properties of silica glass using ReaxFF

    Chowdhury, S.C., Haque, B.Z.G., Gillespie, J.W., 2016. Molecular dynamics simulations of the structure and mechanical properties of silica glass using ReaxFF. Journal of Materials Science 51, 10139–10159. URL:https://doi.org/10.1007/s10853-016-0242-8, doi:10.1007/s10853-016-0242-8

    work page doi:10.1007/s10853-016-0242-8 2016

Show all 52 references
  1. [9]

    Strain rate and orientation dependencies of the strength of single crystalline copper under compression

    Dupont, V., Germann, T.C., 2012. Strain rate and orientation dependencies of the strength of single crystalline copper under compression. Physical Review B 86, 134111. URL:https://link.aps.org/doi/10.1103/PhysRevB.86.134111, doi:10.1103/PhysRevB.86. 134111. publisher: American...

    2012 doi

  2. [10]

    Vitreous 2d silica under tension: From brittle to ductile behaviour

    Ebrahem, F., Bamer, F., Markert, B., 2020. Vitreous 2d silica under tension: From brittle to ductile behaviour. Materials Science and Engineering:A780,139189. URL: https://www.sciencedirect.com/science/article/pii/S0921509320302756,doi: 10.1016/ j.msea.2020.139189

    2020

  3. [11]

    Dynamics of viscoplastic deformation in amorphous solids

    Falk, M.L., Langer, J.S., 1998. Dynamics of viscoplastic deformation in amorphous solids. Physical Review E 57, 7192–7205. URL: https://link.aps.org/doi/10.1103/PhysRevE.57.7192,doi: 10.1103/PhysRevE.57.7192.publisher:AmericanPhysicalSociety

    1998 doi

  4. [12]

    mdpi.com/1996-1944/17/11/2507, doi:10.3390/ma17112507

    Fan,B.,Li,M.,2024.Effectofstrainrateonmechanicaldeformationbehaviorincuzrmetallicglass.Materials17,2507.URL: https://www. mdpi.com/1996-1944/17/11/2507, doi:10.3390/ma17112507. number: 11 Publisher: Multidisciplinary Digital Publishing Institute

    2024 doi

  5. [13]

    Connecting Structural Relaxation with the Low Frequency Modes in a Hard-Sphere Colloidal Glass

    Ghosh, A., Chikkadi, V., Schall, P., Bonn, D., 2011. Connecting Structural Relaxation with the Low Frequency Modes in a Hard-Sphere Colloidal Glass. Physical Review Letters 107, 188303. URL: https://link.aps.org/doi/10.1103/PhysRevLett.107.188303, doi:10.1103/PhysRevLett.107.1...

    2011 doi

  6. [14]

    Elementary principles in statistical mechanics: developed with especial reference to the rational foundations of thermodynamics

    Gibbs, J.W., 1902. Elementary principles in statistical mechanics: developed with especial reference to the rational foundations of thermodynamics. C. Scribner’s sons

    1902

  7. [15]

    Nonequilibrium thermomechanics of gaussian phase packet crystals: Application to the quasistatic quasicontinuum method

    Gupta, P., Ortiz, M., Kochmann, D.M., 2021. Nonequilibrium thermomechanics of gaussian phase packet crystals: Application to the quasistatic quasicontinuum method. Journal of the Mechanics and Physics of Solids 153, 104495

    2021

  8. [16]

    Spínola et al.: Preprint submitted to Elsevier Page 20 of 22 Predicting temperature-dependent failure and transformation zones in 2D silica

    Huang,P.Y.,Kurasch,S.,Srivastava,A.,Skakalova,V.,Kotakoski,J.,Krasheninnikov,A.V.,Hovden,R.,Mao,Q.,Meyer,J.C.,Smet,J.,Muller, D.A.,Kaiser,U.,2012.Directimagingofatwo-dimensionalsilicaglassongraphene.NanoLetters12,1081–1086.doi: 10.1021/nl204423x. Spínola et al.: Preprint submi...

    2012 doi

  9. [17]

    Statistical physics of the yielding transition in amorphous solids

    Karmakar, S., Lerner, E., Procaccia, I., 2010. Statistical physics of the yielding transition in amorphous solids. Physical Review E 82, 055103. URL: https://link.aps.org/doi/10.1103/PhysRevE.82.055103, doi:10.1103/PhysRevE.82.055103. publisher: American Physical Society

    2010 doi

  10. [18]

    A variational approach to coarse graining of equilibrium and non-equilibrium atomistic description at finite temperature

    Kulkarni, Y., Knap, J., Ortiz, M., 2008. A variational approach to coarse graining of equilibrium and non-equilibrium atomistic description at finite temperature. Journal of the Mechanics and Physics of Solids 56, 1417–1449

    2008

  11. [19]

    Localityandnonlocalityinelastoplasticresponsesofamorphoussolids

    Lerner,E.,Procaccia,I.,2009. Localityandnonlocalityinelastoplasticresponsesofamorphoussolids. PhysicalReviewE79,066109. URL: https://link.aps.org/doi/10.1103/PhysRevE.79.066109, doi:10.1103/PhysRevE.79.066109. publisher: American Physical Society

    2009 doi

  12. [20]

    Finite-temperature defect properties from free-energy minimization

    LeSar, R., Najafabadi, R., Srolovitz, D., 1989. Finite-temperature defect properties from free-energy minimization. Physical Review Letters 63, 624

    1989

  13. [21]

    Temperature rise at shear bands in metallic glasses

    Lewandowski, J.J., Greer, A.L., 2006. Temperature rise at shear bands in metallic glasses. Nature Materials 5, 15–18. URL:https: //www.nature.com/articles/nmat1536, doi:10.1038/nmat1536. publisher: Nature Publishing Group

    2006 doi

  14. [22]

    Diffusive molecular dynamics and its application to nanoindentation and sintering

    Li, J., Sarkar, S., Cox, W.T., Lenosky, T.J., Bitzek, E., Wang, Y., 2011. Diffusive molecular dynamics and its application to nanoindentation and sintering. Physical Review B—Condensed Matter and Materials Physics 84, 054103

    2011

  15. [23]

    Crystalline-vitreous interface intwo dimensionalsilica

    Lichtenstein, L.,Heyde, M.,Freund, H.J.,2012. Crystalline-vitreous interface intwo dimensionalsilica. Phys. Rev.Lett. 109,106101. URL: https://link.aps.org/doi/10.1103/PhysRevLett.109.106101, doi:10.1103/PhysRevLett.109.106101

    2012 doi

  16. [24]

    Universal Breakdown of Elasticity at the Onset of Material Failure

    Maloney, C., Lemaître, A., 2004. Universal Breakdown of Elasticity at the Onset of Material Failure. Physical Review Letters 93, 195501. URL: https://link.aps.org/doi/10.1103/PhysRevLett.93.195501, doi:10.1103/PhysRevLett.93.195501

    2004 doi

  17. [25]

    Amorphous systems in athermal, quasistatic shear

    Maloney, C.E., Lemaître, A., 2006. Amorphous systems in athermal, quasistatic shear. Physical Review E 74, 016118. URL:https: //link.aps.org/doi/10.1103/PhysRevE.74.016118, doi:10.1103/PhysRevE.74.016118

    2006 doi

  18. [26]

    Vibrational Modes Identify Soft Spots in a Sheared Disordered Packing

    Manning, M.L., Liu, A.J., 2011. Vibrational Modes Identify Soft Spots in a Sheared Disordered Packing. Physical Review Letters 107, 108302. URL: https://link.aps.org/doi/10.1103/PhysRevLett.107.108302, doi:10.1103/PhysRevLett.107.108302. publisher: American Physical Society

    2011 doi

  19. [27]

    Equationofstatecalculationsbyfastcomputingmachines

    Metropolis,N.,Rosenbluth,A.W.,Rosenbluth,M.N.,Teller,A.H.,Teller,E.,1953. Equationofstatecalculationsbyfastcomputingmachines. The journal of chemical physics 21, 1087–1092

    1953

  20. [28]

    Soft Modes and Nonaffine Rearrangements in the Inherent Structures of SupercooledLiquids

    Mosayebi, M., Ilg, P., Widmer-Cooper, A., Del Gado, E., 2014. Soft Modes and Nonaffine Rearrangements in the Inherent Structures of SupercooledLiquids. PhysicalReviewLetters112,105503. URL: https://link.aps.org/doi/10.1103/PhysRevLett.112.105503, doi:10.1103/PhysRevLett.112.10...

    2014 doi

  21. [29]

    Springer International Publishing, Cham, pp

    Mozgawa,W.,Sitarz,M.,Król,M.,2019.SpectroscopicCharacterizationofSilicateAmorphousMaterials,in:Koleżyński,A.,Król,M.(Eds.), Molecular Spectroscopy—Experiment and Theory: From Molecules to Functional Materials. Springer International Publishing, Cham, pp. 457–481. doi:10.1007/9...

    2019 doi

  22. [30]

    Molecular dynamics studies of brittle failure in silica: effect of thermal vibrations

    Ochoa, R., Swiler, T.P., Simmons, J.H., 1991. Molecular dynamics studies of brittle failure in silica: effect of thermal vibrations. JournalofNon-CrystallineSolids128,57–68. URL: https://www.sciencedirect.com/science/article/pii/0022309391907763, doi:10.1016/0022-3093(91)90776-3

    1991

  23. [31]

    Connecting Local Yield Stresses with Plastic Activity in Amorphous Solids

    Patinet, S., Vandembroucq, D., Falk, M.L., 2016. Connecting Local Yield Stresses with Plastic Activity in Amorphous Solids. Physical Review Letters 117, 045501. URL:https://link.aps.org/doi/10.1103/PhysRevLett.117.045501, doi:10.1103/PhysRevLett. 117.045501

    2016 doi

  24. [32]

    Moleculardynamicsstudiesofstress-strainbehaviorofsilicaglass under a tensile load

    Pedone,A.,Malavasi,G.,Menziani,M.C.,Segre,U.,Cormack,A.N.,2008. Moleculardynamicsstudiesofstress-strainbehaviorofsilicaglass under a tensile load. Chemistry of Materials 20, 4356–4366. URL:https://doi.org/10.1021/cm800413v, doi:10.1021/cm800413v. publisher: American Chemical Society

    2008 doi

  25. [33]

    Predicting plasticity in disordered solids fromstructuralindicators

    Richard,D.,Ozawa,M.,Patinet,S.,Stanifer,E.,Shang,B.,Ridout,S.A.,Xu,B.,Zhang,G.,Morse,P.K.,Barrat,J.L.,Berthier,L.,Falk,M.L., Guan, P., Liu, A.J., Martens, K., Sastry, S., Vandembroucq, D., Lerner, E., Manning, M.L., 2020. Predicting plasticity in disordered solids fromstructur...

    2020 doi

  26. [34]

    Modeling the mechanics of amorphous solids at different length scale and time scale

    Rodney, D., Tanguy, A., Vandembroucq, D., 2011. Modeling the mechanics of amorphous solids at different length scale and time scale. Modelling and Simulation in Materials Science and Engineering 19, 083001. URL:https://dx.doi.org/10.1088/0965-0393/19/8/ 083001, doi:10.1088/096...

    2011 doi

  27. [35]

    Extended molecular dynamics: Seamless temporal coarse-graining and transition between deterministic and probabilistic paradigms

    Romero, I., Ortiz, M., 2023. Extended molecular dynamics: Seamless temporal coarse-graining and transition between deterministic and probabilistic paradigms. European Journal of Mechanics - A/Solids 97, 104858. URL:https://www.sciencedirect.com/science/ article/pii/S0997753822...

    2023

  28. [36]

    Modelling the atomic arrangement of amorphous 2d silica: a network analysis

    Roy, P.K., Heyde, M., Heuer, A., 2018. Modelling the atomic arrangement of amorphous 2d silica: a network analysis. Physical Chemistry Chemical Physics 20, 14725–14739

    2018

  29. [37]

    Avalanches in the Athermal Quasistatic Limit of Sheared Amorphous Solids: An Atomistic Perspective

    Ruscher, C., Rottler, J., 2021. Avalanches in the Athermal Quasistatic Limit of Sheared Amorphous Solids: An Atomistic Perspective. Tribology Letters 69, 64. URL:https://doi.org/10.1007/s11249-021-01439-5, doi:10.1007/s11249-021-01439-5

    2021 doi

  30. [38]

    Gnn-assisted phase space integration with application to atomistics

    Saxena, S., Bastek, J.H., Spinola, M., Gupta, P., Kochmann, D.M., 2023. Gnn-assisted phase space integration with application to atomistics. Mechanics of Materials 182, 104681

    2023

  31. [39]

    A fast atomistic approach to finite-temperature surface elasticity of crystalline solids

    Saxena, S., Spinola, M., Gupta, P., Kochmann, D.M., 2022. A fast atomistic approach to finite-temperature surface elasticity of crystalline solids. Computational Materials Science 211, 111511

    2022

  32. [40]

    Criticality in the fracture of silica glass: Insights from molecular mechanics

    Shekh Alshabab, S., Markert, B., Bamer, F., 2024. Criticality in the fracture of silica glass: Insights from molecular mechanics. Physical ReviewE109,034110. URL: https://link.aps.org/doi/10.1103/PhysRevE.109.034110,doi: 10.1103/PhysRevE.109.034110. publisher: American Physica...

    2024 doi

  33. [41]

    Finite-temperature grain boundary properties from quasistatic atomistics

    Spínola, M., Saxena, S., Gupta, P., Runnels, B., Kochmann, D.M., 2024. Finite-temperature grain boundary properties from quasistatic atomistics. Computational Materials Science 244, 113270. URL: https://www.sciencedirect.com/science/article/pii/ Spínola et al.: Preprint submit...

    2024

  34. [42]

    Approximate calculation of multiple integrals

    Stroud, A.H., 1971. Approximate calculation of multiple integrals. Englewood Cliffs, N.J.: Prentice-Hall. URL:http://archive.org/ details/approximatecalcu0000stro_b8j7

    1971

  35. [43]

    Modeling materials: continuum, atomistic and multiscale techniques

    Tadmor, E.B., Miller, R.E., 2011. Modeling materials: continuum, atomistic and multiscale techniques. Cambridge university press

    2011

  36. [44]

    Plastic response of a 2D Lennard-Jones amorphous solid: Detailed analysis of the local rearrangements at very slow strain rate

    Tanguy, A., Leonforte, F., Barrat, J.L., 2006. Plastic response of a 2D Lennard-Jones amorphous solid: Detailed analysis of the local rearrangements at very slow strain rate. The European Physical Journal E 20, 355–364. URL: https://doi.org/10.1140/epje/ i2006-10024-2, doi:10....

    2006 doi

  37. [45]

    LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales

    Thompson, A.P., Aktulga, H.M., Berger, R., Bolintineanu, D.S., Brown, W.M., Crozier, P.S., in ’t Veld, P.J., Kohlmeyer, A., Moore, S.G., Nguyen, T.D., Shan, R., Stevens, M.J., Tranchida, J., Trott, C., Plimpton, S.J., 2022. LAMMPS - a flexible simulation tool for particle-base...

    2022

  38. [46]

    Local elasticity map and plasticity in a model Lennard-Jones glass

    Tsamados, M., Tanguy, A., Goldenberg, C., Barrat, J.L., 2009. Local elasticity map and plasticity in a model Lennard-Jones glass. Physical Review E 80, 026112. URL: https://link.aps.org/doi/10.1103/PhysRevE.80.026112, doi:10.1103/PhysRevE.80.026112. publisher: American Physica...

    2009 doi

  39. [47]

    Statistical mechanics: theory and molecular simulation

    Tuckerman, M.E., 2023. Statistical mechanics: theory and molecular simulation. Oxford university press

    2023

  40. [48]

    Frequency factors and isotope effects in solid state rate processes

    Vineyard, G.H., 1957. Frequency factors and isotope effects in solid state rate processes. Journal of Physics and Chemistry of Solids 3, 121–127

    1957

  41. [49]

    Anomalous temperature dependence of elastic limit in metallic glasses

    Wang, Y., Liu, J., Jiang, J.Z., Cai, W., 2024. Anomalous temperature dependence of elastic limit in metallic glasses. Nature Communications 15,171. URL: https://www.nature.com/articles/s41467-023-44048-7,doi: 10.1038/s41467-023-44048-7.publisher:Nature Publishing Group

    2024 doi

  42. [50]

    Localized soft modes and the supercooled liquid’s irreversible passage through its configuration space

    Widmer-Cooper, A., Perry, H., Harrowell, P., Reichman, D.R., 2009. Localized soft modes and the supercooled liquid’s irreversible passage through its configuration space. The Journal of Chemical Physics 131, 194508. URL:https://doi.org/10.1063/1.3265983, doi:10.1063/1.3265983

    2009 doi

  43. [51]

    Topologyofvibrationalmodespredictsplasticeventsinglasses

    Wu,Z.W.,Chen,Y.,Wang,W.H.,Kob,W.,Xu,L.,2023. Topologyofvibrationalmodespredictsplasticeventsinglasses. NatureCommunica- tions 14, 2955. URL:https://www.nature.com/articles/s41467-023-38547-w, doi:10.1038/s41467-023-38547-w. publisher: Nature Publishing Group

    2023 doi

  44. [52]

    Two-dimensional ultrathin silica films

    Zhong, J.Q., Freund, H.J., 2022. Two-dimensional ultrathin silica films. Chemical Reviews 122, 11172–11246. Spínola et al.: Preprint submitted to Elsevier Page 22 of 22

    2022

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.