-
Polarizable Potentials For Metals: The Density Readjusting Embedded Atom Method (DR-EAM)
Authors:
Hemanta Bhattarai,
Kathie E. Newman,
J. Daniel Gezelter
Abstract:
In simulations of metallic interfaces, a critical aspect of metallic behavior is missing from the some of the most widely used classical molecular dynamics force fields. We present a modification of the embedded atom method (EAM) which allows for electronic polarization of the metal by treating the valence density around each atom as a fluctuating dynamical quantity. The densities are represented…
▽ More
In simulations of metallic interfaces, a critical aspect of metallic behavior is missing from the some of the most widely used classical molecular dynamics force fields. We present a modification of the embedded atom method (EAM) which allows for electronic polarization of the metal by treating the valence density around each atom as a fluctuating dynamical quantity. The densities are represented by a set of additional fluctuating variables (and their conjugate momenta) which are propagated along with the nuclear coordinates. This ``density readjusting EAM'' (DR-EAM) preserves nearly all of the useful qualities of traditional EAM, including bulk elastic properties and surface energies. However, it also allows valence electron density to migrate through the metal in response to external perturbations. We show that DR-EAM can successfully model polarization in response to external charges, capturing the image charge effect in atomistic simulations. DR-EAM also captures some of the behavior of metals in the presence of uniform electric fields, predicting surface charging and shielding internal to the metal. We further show that it predicts charge transfer between the constituent atoms in alloys, leading to novel predictions about unit cell geometries in layered L$1_0$ structures.
△ Less
Submitted 30 March, 2019;
originally announced April 2019.
-
Real space electrostatics for multipoles. III. Dielectric Properties
Authors:
Madan Lamichhane,
Thomas Parsons,
Kathie E. Newman,
J. Daniel Gezelter
Abstract:
In the first two papers in this series, we developed new shifted potential (SP), gradient shifted force (GSF), and Taylor shifted force (TSF) real-space methods for multipole interactions in condensed phase simulations. Here, we discuss the dielectric properties of fluids that emerge from simulations using these methods. Most electrostatic methods (including the Ewald sum) require correction to th…
▽ More
In the first two papers in this series, we developed new shifted potential (SP), gradient shifted force (GSF), and Taylor shifted force (TSF) real-space methods for multipole interactions in condensed phase simulations. Here, we discuss the dielectric properties of fluids that emerge from simulations using these methods. Most electrostatic methods (including the Ewald sum) require correction to the conducting boundary fluctuation formula for the static dielectric constants, and we discuss the derivation of these corrections for the new real space methods. For quadrupolar fluids, the analogous material property is the quadrupolar susceptibility. As in the dipolar case, the fluctuation formula for the quadrupolar susceptibility has corrections that depend on the electrostatic method being utilized. One of the most important effects measured by both the static dielectric and quadrupolar susceptibility is the ability to screen charges embedded in the fluid. We use potentials of mean force between solvated ions to discuss how geometric factors can lead to distance-dependent screening in both quadrupolar and dipolar fluids.
△ Less
Submitted 17 August, 2016;
originally announced August 2016.
-
Theoretical Study of the Stability of $\mathrm{Au_nGe_{m}}$ Clusters
Authors:
Danielle McDermott,
Kathie E. Newman
Abstract:
We investigate the properties of clusters formed from two connected $\mathrm{Ge_{m}}$ cage-like clusters, such as experimentally synthesized $\mathrm{Au_3Ge_{18}^{5-}}$, using first-principles DFT methods. We focus particularly on $\mathrm{Au_nGe_{12}^{q-}}$ formed from $\mathrm{Ge_{6}}$, where $n=0-3$ and $q=0,2$. We examine the geometries, electronic structure, and thermal excitations of these c…
▽ More
We investigate the properties of clusters formed from two connected $\mathrm{Ge_{m}}$ cage-like clusters, such as experimentally synthesized $\mathrm{Au_3Ge_{18}^{5-}}$, using first-principles DFT methods. We focus particularly on $\mathrm{Au_nGe_{12}^{q-}}$ formed from $\mathrm{Ge_{6}}$, where $n=0-3$ and $q=0,2$. We examine the geometries, electronic structure, and thermal excitations of these clusters using the SIESTA code, where we test stability using short molecular dynamics simulations. Our examination of intercluster bridges between $\mathrm{Ge_{m}}$, formed of either Ge-Ge or Au-Ge bonds, indicates that the overlap of interacting molecular orbitals of $\mathrm{Ge_{m}}$~can either bind a cluster together or tear it apart depending on the orientation of the bridging atoms with respect to the cages. We characterize the properties of stable isomers of neutrally charged $\mathrm{AuGe_{12}}$ and $\mathrm{Au_2Ge_{12}}$, seeing that radially directed molecular orbitals stabilize $\mathrm{AuGe_{12}}$ while a geometric asymmetry stabilizes $\mathrm{Au_2Ge_{12}}$ and $\mathrm{Au_3Ge_{18}}$. We examine a two-dimensional $_2^\infty[\mathrm{Au_2Ge_{6}}]$ structure, finding it more stable than other periodic $[\mathrm{Au_nGe_{6}}]$ subunits. We observe no stable neutral isomers of $\mathrm{Au_3Ge_{12}}$, but see that additional charge supports isomers of both $\mathrm{Au_2Ge_{12}}$ and $\mathrm{Au_3Ge_{12}}$.
△ Less
Submitted 5 March, 2015; v1 submitted 12 August, 2014;
originally announced August 2014.
