A Caltech Library Service

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

Gupta, Prateek and Ortiz, Michael and Kochmann, Dennis M. (2021) Nonequilibrium thermomechanics of Gaussian phase packet crystals: Application to the quasistatic quasicontinuum method. Journal of the Mechanics and Physics of Solids, 153 . Art. No. 104495. ISSN 0022-5096. doi:10.1016/j.jmps.2021.104495.

[img] PDF - Published Version
Creative Commons Attribution.


Use this Persistent URL to link to this item:


The quasicontinuum (QC) method was originally introduced to bridge across length scales by coarse-graining an atomistic ensemble to significantly larger continuum scales at zero temperature, thus overcoming the crucial length-scale limitation of classical atomic-scale simulation techniques while solely relying on atomic-scale input (in the form of interatomic potentials). An associated challenge lies in bridging across time scales to overcome the time-scale limitations of atomistics at finite temperature. To address the biggest challenge, bridging across both length and time scales, only a few techniques exist, and most of those are limited to conditions of constant temperature. Here, we present a new general strategy for the space–time coarsening of an atomistic ensemble, which introduces thermomechanical coupling. Specifically, we evolve the statistics of an atomistic ensemble in phase space over time by applying the Liouville equation to an approximation of the ensemble’s probability distribution (which further admits a variational formulation). To this end, we approximate a crystalline solid as a lattice of lumped correlated Gaussian phase packets occupying atomic lattice sites, and we investigate the resulting quasistatics and dynamics of the system. By definition, phase packets account for the dynamics of crystalline lattices at finite temperature through the statistical variances of atomic momenta and positions. We show that momentum–space correlation allows for an exchange between potential and kinetic contributions to the crystal’s Hamiltonian. Consequently, local adiabatic heating due to atomic site motion is captured. Moreover, in the quasistatic limit, the governing equations reduce to the minimization of thermodynamic potentials (similar to maximum-entropy formulation previously introduced for finite-temperature QC), and they yield the local equation of state, which we derive for isothermal, isobaric, and isentropic conditions. Since our formulation without interatomic correlations precludes irreversible heat transport, we demonstrate its combination with thermal transport models to describe realistic atomic-level processes, and we discuss opportunities for capturing atomic-level thermal transport by including interatomic correlations in the Gaussian phase packet formulation. Overall, our Gaussian phase packet approach offers a promising avenue for finite-temperature non-equilibrium quasicontinuum techniques, which may be combined with thermal transport models and extended to other approximations of the probability distribution as well as to exploit the variational structure.

Item Type:Article
Related URLs:
URLURL TypeDescription
Gupta, Prateek0000-0003-3666-0257
Ortiz, Michael0000-0001-5877-4824
Kochmann, Dennis M.0000-0002-9112-6615
Additional Information:© 2021 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license. Received 10 January 2021, Revised 11 April 2021, Accepted 13 May 2021, Available online 17 May 2021. The support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 770754) is gratefully acknowledged. The authors thank Miguel Angel Spinola Fandila for help with the LAMMPS simulations used for numerical validation. CRediT authorship contribution statement: Prateek Gupta: Conceptualization, Methodology, Software, Validation, Investigation, Data curation, Writing - original draft. Michael Ortiz: Conceptualization, Methodology. Dennis M. Kochmann: Conceptualization, Methodology, Writing - review & editing, Supervision, Funding acquisition. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Funding AgencyGrant Number
European Research Council (ERC)770754
Subject Keywords:Quasicontinuum; Multiscale modeling; Atomistics; Non-equilibrium statistical mechanics; Updated-Lagrangian
Record Number:CaltechAUTHORS:20210518-103040351
Persistent URL:
Official Citation:Prateek Gupta, Michael Ortiz, Dennis M. Kochmann, Nonequilibrium thermomechanics of Gaussian phase packet crystals: Application to the quasistatic quasicontinuum method, Journal of the Mechanics and Physics of Solids, Volume 153, 2021, 104495, ISSN 0022-5096, (
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:109175
Deposited By: Tony Diaz
Deposited On:19 May 2021 18:27
Last Modified:27 May 2021 21:48

Repository Staff Only: item control page