Source author record

E. R. Smith

E. R. Smith appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2026arXiv

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.

preprint2022arXiv

The Intrinsic Fragility of the Liquid-Vapor Interface: A Stress Network Perspective

The evolution of the liquid-vapour interface of a Lennard-Jones fluid is examined with molecular dynamics simulations using the intrinsic sampling method. Results suggest, in agreement with capillary wave theory, clear damping of the density profiles as the temperature is increased. We identify a linear variation of the space-filling nature (fractal dimension) of the stress-clusters at the intrinsic surface with increasing surface tension, or equivalently, with decreasing temperature. A percolation analysis of these stress networks indicates that the stress field is more disjointed at higher temperatures. This leads to more fragile interfaces that result in a reduction in surface tension at higher temperature.

preprint2015arXiv

A Molecular Dynamics Simulation of the Turbulent Couette Minimal Flow Unit

A molecular dynamics (MD) simulation of planar Couette flow is presented for the minimal channel in which turbulence structures can be sustained. Evolution over a single breakdown and regeneration cycle is compared to computational fluid dynamics (CFD) simulations. Qualitative similar structures are observed and turbulent statistics show excellent quantitative agreement. The molecular scale law of the wall is presented in which stick-slip molecular wall-fluid interactions replace the no-slip conditions. The impact of grid resolution is explored and the observed structures are seen to be dependant on averaging time and length scales. The kinetic energy spectra show a range of scales are present in the molecular system and that spectral content is dependent on the grid resolution employed. The subgrid velocity of the molecules is compared to spatial averaged velocity using joint probability density functions. Molecular trajectories, diffusions and Lagrangian statistics are presented. The importance of sub-grid scales, relevance of the Kolmogorov lengthscale and implications of molecular turbulence are discussed.

preprint2012arXiv

Charge Transport in one Dimension:Dissipative and Non-Dissipative Space-Charge Limited Currents

We consider charge transport in nanopores where the dielectric constant inside the nanopore is much greater than in the surrounding material, so that the flux of the electric fields due to the charges is almost entirely confined to the nanopore. That means that we may model the electric fields due to charge densities in the nanopore in terms of average properties across the nanopore as solutions of one dimensional Poisson equations. We develop basic equations for an M component system using equations of continuity to relate concentrations to currents, and flux equations relating currents to concentration gradients and conductivities. We then derive simplified scaled versions of the equations. We develop exact solutions for the one component case in a variety of boundary conditions using a Hopf-Cole transformation, Fourier series, and periodic solutions of the Burgers equation. These are compared with a simpler model in which the scaled diffusivity is zero so that all charge motion is driven by the electric field. In this non-dissipative case, recourse to an admissibility condition is utilised to obtain the physically relevant weak solution of a Riemann problem concerning the electric field. It is shown that the admissibility condition is Poynting's theorem.

preprint2012arXiv

Control-volume representation of molecular dynamics

A Molecular Dynamics (MD) parallel to the Control Volume (CV) formulation of fluid mechanics is developed by integrating the formulas of Irving and Kirkwood, J. Chem. Phys. 18, 817 (1950) over a finite cubic volume of molecular dimensions. The Lagrangian molecular system is expressed in terms of an Eulerian CV, which yields an equivalent to Reynolds' Transport Theorem for the discrete system. This approach casts the dynamics of the molecular system into a form that can be readily compared to the continuum equations. The MD equations of motion are reinterpreted in terms of a Lagrangian-to-Control-Volume (\CV) conversion function $\vartheta_{i}$, for each molecule $i$. The \CV function and its spatial derivatives are used to express fluxes and relevant forces across the control surfaces. The relationship between the local pressures computed using the Volume Average (VA, Lutsko, J. Appl. Phys 64, 1152 (1988)) techniques and the Method of Planes (MOP, Todd et al, Phys. Rev. E 52, 1627 (1995)) emerges naturally from the treatment. Numerical experiments using the MD CV method are reported for equilibrium and non-equilibrium (start-up Couette flow) model liquids, which demonstrate the advantages of the formulation. The CV formulation of the MD is shown to be exactly conservative, and is therefore ideally suited to obtain macroscopic properties from a discrete system.