Journal of the Mechanics and Physics of Solids Nonequilibrium thermomechanics of Gaussian phase packet crystals: Application to the quasistatic quasicontinuum method

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 nonequilibrium 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. equation, which yields the equations of motion for the phase space parameters. We also show that the interatomic correlations are of fundamental importance for modeling interatomic heat flux. Since such correlations increase the degrees of freedom significantly, they are neglected in subsequent sections while retaining the phase space correlation of each atom, which we later identify as the thermal momentum . In Section 3 we discuss the importance of thermal momentum in dynamics and the implications of its vanishing limit in quasistatics. We show that, in the latter limit, the governing equations define the local thermomechanical equilibrium of the system and yield a rate-independent local thermal equation of state. To capture thermal transport, we adopt Venturini et al.’s linear Onsager kinetics-based model. We also discuss the time scale imposed by the Onsager