Atomistic long-term simulation of heat and mass transport

We formulate a theory of non-equilibrium statistical thermodynamics for ensembles of atoms or molecules. The theory is an application of Jaynes׳ maximum entropy principle, which allows the statistical treatment of systems away from equilibrium. In particular, neither temperature nor atomic fractions are required to be uniform but instead are allowed to take different values from particle to particle. In addition, following the Coleman-Noll method of continuum thermodynamics we derive a dissipation inequality expressed in terms of discrete thermodynamic fluxes and forces. This discrete dissipation inequality effectively sets the structure for discrete kinetic potentials that couple the microscopic field rates to the corresponding driving forces, thus resulting in a closed set of equations governing the evolution of the system. We complement the general theory with a variational meanfield theory that provides a basis for the formulation of computationally tractable approximations. We present several validation cases, concerned with equilibrium properties of alloys, heat conduction in silicon nanowires and hydrogen desorption from palladium thin films, that demonstrate the range and scope of the method and assess its fidelity and predictiveness. These validation cases are characterized by the need or desirability to account for atomic-level properties while simultaneously entailing time scales much longer than those accessible to direct molecular dynamics. The ability of simple meanfield models and discrete kinetic laws to reproduce equilibrium properties and long-term behavior of complex systems is remarkable.