BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Doyl Dickel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher 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.
A recent proposal for practical calculation of vibrational mode lifetimes is tested on simple, low-dimensional anharmonic models. The proposed scheme approximates the mode lifetime in terms of ensemble averages of specific functions in phase-space; various levels of approximation correspond to ensemble moments of the Liouvillian. It is shown that, for systems where the vibrational density of states is well-approximated by a single broadened peak, the fourth-moment approximation works well over the full range of temperature.
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.