Reactive Molecular Dynamics at Constant Pressure via Nonreactive Force Fields: Extending the Empirical Valence Bond Method to the Isothermal-Isobaric Ensemble.

The Empirical Valence Bond (EVB) method offers a suitable framework to obtain reactive potentials through the coupling of nonreactive force fields. In this formalism, most of the implemented coupling terms are built using functional forms that depend on spatial coordinates, while parameters are fitted against reference data to model the change of chemistry between the participating nonreactive states. In this work, we demonstrate that the use of such coupling terms precludes the computation of the stress tensor for condensed phase systems and prevents the possibility to carry out EVB molecular dynamics in the isothermal-isobaric (NPT) ensemble. Alternatively, we make use of coupling terms that depend on the energy gaps, defined as the energy differences between the participating nonreactive force fields, and derive a general expression for the EVB stress tensor suitable for computation. Implementation of this new methodology is tested for a model of a single reactive malonaldehyde solvated in nonreactive water. Mass densities and probability distributions for the values of the energy gaps computed in the NPT ensemble reveal a negligible role of the reactive potential in the limit of low concentrated solutions, thus corroborating for the first time the validity of approximations based on the canonical NVT ensemble, customarily adopted for EVB simulations. The presented formalism also aims to contribute to future implementations and extensions of the EVB method to research the limit of highly concentrated solutions.

Charge-State-Dependent Diffusion of Individual Gold Adatoms on Ionic Thin NaCl Films.

It is known that individual metal atoms on insulating ionic films can occur in several different (meta)stable charge states, which can be reversibly switched in a controlled fashion. Here we show that the diffusion of gold adatoms on NaCl thin films depends critically on their charge state. Surprisingly, the anionic species has a lower diffusion barrier than the neutral one. Furthermore, for the former we observe that the diffusion atop a bilayer of NaCl is strongly influenced by the interface between NaCl and the underlying copper substrate. This effect disappears for a trilayer of NaCl. These observations open the prospect of controlling the diffusion properties of individual metal atoms on thin insulating films.

Combined Role of Biaxial Strain and Nonstoichiometry for the Electronic, Magnetic, and Redox Properties of Lithiated Metal-Oxide Films: The LiMn2O4 Case.

Understanding the interplay between strain and nonstoichiometry for the electronic, magnetic, and redox properties of LiMn2O4 films is essential for their development as Li-ion battery (LIB) cathodes, photoelectrodes, and systems for sustainable spintronics applications as well as for emerging applications that combine these technologies. Here, density functional theory (DFT) simulations suggest that compressive strain increases the reduction drive of (111) LiMn2O4 films by inducing >1 eV upshift of the valence band edge. The DFT results indicate that, regardless of the crystallographic orientation for the LiMn2O4 film, biaxial expansion increases the magnetic moments of the Mn atoms. Conversely, biaxial compression reduces them. For ferromagnetic films, these changes can be substantial and as large as over 4 Bohr magnetons per unit cell over the simulated range of strain (from -6 to +3%). The DFT simulations also uncover a compensation mechanism whereby strain induces opposite changes in the magnetic moment of the Mn and O atoms, leading to an overall constant magnetic moment for the ferromagnetic films. The calculated strain-induced changes in atomic magnetic moments reflect modifications in the local electronic hybridization of both the Mn and O atoms, which in turn suggests strain-tunable, local chemical, and electrochemical reactivity. Several energy-favored (110) and (111) ferromagnetic surfaces turn out to be half-metallic with minority-spin band gaps as large as 3.2 eV and compatible with spin-dependent electron-transport and possible spin-dependent electrochemical and electrocatalytic properties. The resilience of the ferromagnetic, half-metallic states to surface nonstoichiometry and compositional changes invites exploration of the potential of LiMn2O4 thin films for sustainable spintronic applications beyond state-of-the-art, rare-earth metal-based, ferromagnetic half-metallic oxides.

Reorganization energy upon charging a single molecule on an insulator measured by atomic force microscopy.

Intermolecular single-electron transfer on electrically insulating films is a key process in molecular electronics1-4 and an important example of a redox reaction5,6. Electron-transfer rates in molecular systems depend on a few fundamental parameters, such as interadsorbate distance, temperature and, in particular, the Marcus reorganization energy 7 . This crucial parameter is the energy gain that results from the distortion of the equilibrium nuclear geometry in the molecule and its environment on charging8,9. The substrate, especially ionic films 10 , can have an important influence on the reorganization energy11,12. Reorganization energies are measured in electrochemistry 13 as well as with optical14,15 and photoemission spectroscopies16,17, but not at the single-molecule limit and nor on insulating surfaces. Atomic force microscopy (AFM), with single-charge sensitivity18-22, atomic-scale spatial resolution 20 and operable on insulating films, overcomes these challenges. Here, we investigate redox reactions of single naphthalocyanine (NPc) molecules on multilayered NaCl films. Employing the atomic force microscope as an ultralow current meter allows us to measure the differential conductance related to transitions between two charge states in both directions. Thereby, the reorganization energy of NPc on NaCl is determined as (0.8 ± 0.2) eV, and density functional theory (DFT) calculations provide the atomistic picture of the nuclear relaxations on charging. Our approach presents a route to perform tunnelling spectroscopy of single adsorbates on insulating substrates and provides insight into single-electron intermolecular transport.

Frontier molecular orbitals of a single molecule adsorbed on thin insulating films supported by a metal substrate: electron and hole attachment energies.

We present calculations of vertical electron and hole attachment energies to the frontier orbitals of a pentacene molecule absorbed on multi-layer sodium chloride films supported by a copper substrate using a simplified density functional theory (DFT) method. The adsorbate and the film are treated fully within DFT, whereas the metal is treated implicitly by a perfect conductor model. We find that the computed energy gap between the highest and lowest unoccupied molecular orbitals-HOMO and LUMO -from the vertical attachment energies increases with the thickness of the insulating film, in agreement with experiments. This increase of the gap can be rationalised in a simple dielectric model with parameters determined from DFT calculations and is found to be dominated by the image interaction with the metal. We find, however, that this simplified model overestimates the downward shift of the energy gap in the limit of an infinitely thick film.

Chemically Selective Alternatives to Photoferroelectrics for Polarization-Enhanced Photocatalysis: The Untapped Potential of Hybrid Inorganic Nanotubes.

Linear-scaling density functional theory simulation of methylated imogolite nanotubes (NTs) elucidates the interplay between wall-polarization, bands separation, charge-transfer excitation, and tunable electrostatics inside and outside the NT-cavity. The results suggest that integration of polarization-enhanced selective photocatalysis and chemical separation into one overall dipole-free material should be possible. Strategies are proposed to increase the NT polarization for maximally enhanced electron-hole separation.

A simplified density functional theory method for investigating charged adsorbates on an ultrathin, insulating film supported by a metal substrate.

A simplified density functional theory (DFT) method for investigating charged adsorbates on an ultrathin, insulating film supported by a metal substrate is developed and presented. This new method is based on a previous DFT development that uses a perfect conductor (PC) model to approximate the electrostatic response of the metal substrate, while the film and the adsorbate are both treated fully within DFT (Scivetti and Persson 2013 J. Phys.: Condens. Matter 25 355006). The missing interactions between the metal substrate and the insulating film in the PC approximation are modelled by a simple force field (FF). The parameters of the PC model and the force field are obtained from DFT calculations of the film and the substrate, here shown explicitly for a NaCl bilayer supported by a Cu(100) surface. In order to obtain some of these parameters and the polarizability of the force field, we have to include an external, uniformly charged plane in the DFT calculations, which has required the development of a periodic DFT formalism to include such a charged plane in the presence of a metal substrate. This extension and implementation should be of more general interest and applicable to other challenging problems, for instance, in electrochemistry. As illustrated for the gold atom on the NaCl bilayer supported by a Cu(100) surface, our new DFT-PC-FF method allows us to handle different charge states of adsorbates in a controlled and accurate manner with a considerable reduction of the computational time. In addition, it is now possible to calculate vertical transition and reorganization energies for the charging and discharging of adsorbates that cannot be obtained by current DFT methodologies that include the metal substrate. We find that the computed vertical transition energy for charging of the gold adatom is in good agreement with experiments.

The electrostatic interaction of an external charged system with a metal surface: a simplified density functional theory approach.

As a first step to meet the challenge to calculate the electronic structure and total energy of charged states of atoms and molecules adsorbed on ultrathin insulating films supported by a metallic substrate using density functional theory (DFT), we have developed a simplified new DFT scheme that only describes the electrostatic interaction of an external charged system with a metal surface. This purely electrostatic interaction is obtained from the assumption that the electron densities of the two fragments (charged system and metal surface) are non-overlapping and by neglecting non-local exchange-correlation effects such as the van der Waals interactions between the two fragments. In addition, the response of the metal surface to the electrostatic potential from the charged system is treated to linear order, whereas the charged system is treated fully within DFT. In particular, we consider the classical perfect conductor model for the metal response, although our formalism is not limited to this approximation. To test the computational implementation of this new scheme, we have considered the case of a Na(+) cation interacting with a perfect conductor. The application of this new methodology to realistic problems involving charged systems adsorbed on insulating films supported by a metal surface are deferred to a separate following publication.