The Quantum Method of Planes -- Local Pressure Definitions for Machine Learning Potentials
Stress, or pressure, is a central quantity in engineering and remains vital in molecular modelling. However, the commonly used virial stress tensor is invalid for an inhomogeneous fluid, which is essential in fluid dynamics and non-equilibrium molecular dynamics (NEMD) simulation. This is solved by using the method of planes (MoP), a mechanical form of pressure, simply interpreted as the force divided by area, yet is derived from the firm foundations of statistical mechanics. We present an extension of MoP stress1 to the MACE potential, a particular form of machine learning (ML) potentials allowing the incorporation of quantum mechanical (QM) physics into classical simulation. We present the derivation of this local stress for the MACE potential using the theoretical framework set out by Irving and Kirkwood 2 . For the test case of an interface between water and Zirconium Oxide, we show that the MoP measures the correct force balance while the virial form fails. Further, we demonstrate that this planar definition of stress is valid arbitrarily far from equilibrium, showing exact conservation every timestep in a control volume bounded by MoP planes. This links the stress directly to the conservation equations and demonstrates the validity in non equilibrium molecular dynamics (NEMD) systems. All code to reproduce these validations for any MACE system, together with ASE accelerated code to calculate the MoP, is provided as open source. This work helps build the foundation to extend the ML revolution in materials to NEMD and molecular fluid dynamics modelling.