BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Doyl Dickel contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We develop an Fe-C-H interatomic potential based on the modified embedded-atom method (MEAM) formalism based on density functional theory to enable large-scale modular dynamics simulations of carbon steel and hydrogen.
While machine learning approaches have been successfully used to represent interatomic potentials, their speed has typically lagged behind conventional formalisms. This is often due to the complexity of the structural fingerprints used to describe the local atomic environment and the large cutoff radii and neighbor lists used in the calculation of these fingerprints. Even recent machine learned methods are at least 10 times slower than traditional formalisms. An implementation of a rapid artificial neural network (RANN) style potential in the LAMMPS molecular dynamics package is presented here which utilizes angular screening to reduce computational complexity without reducing accuracy. For the smallest neural network architectures, this formalism rivals the modified embedded atom method (MEAM) for speed and accuracy, while the networks approximately one third as fast as MEAM were capable of reproducing the training database with chemical accuracy. The numerical accuracy of the LAMMPS implementation is assessed by verifying conservation of energy and agreement between calculated forces and pressures and the observed derivatives of the energy as well as by assessing the stability of the potential in dynamic simulation. The potential style is tested using a force field for magnesium and the computational efficiency for a variety of architectures is compared to a traditional potential models as well as alternative ANN formalisms. The predictive accuracy is found to rival that of slower methods.
In a two-part publication, we propose and analyze a formal foundation for practical calculations of vibrational mode lifetimes in solids. The approach is based on a recursion method analysis of the Liouvillian. In the first part, we derived the lifetime of vibrational modes in terms of moments of the power spectrum of the Liouvillian as projected onto the relevant subspace of phase space. In practical terms, the moments are evaluated as ensemble averages of well-defined operators, meaning that the entire calculation is to be done with Monte Carlo. In this second part, we present a numerical analysis of a simple anharmonic model of lattice vibrations which exhibits two regimes of behavior, at low temperature and at high temperature. Our results show that, for this simple model, the mode lifetime as a function of temperature and wavevector can be simply approximated as a function of the shift in frequency from the harmonic limit. We next compare these calculations, obtained using both Monte Carlo and computationally intensive molecular dynamics, with those using the lowest order moment formalism from the Part I. We show that, in the high-temperature regime, the lowest order approximation gives a reliable approximation to the calculated lifetimes. The results also show that extension to at least fourth moment is required to obtain reliable results over a full range of temperatures.