Search

Machine-learned interatomic potentials are fast becoming an indispensable tool in computational materials science. One approach is the ephemeral data-derived potential (EDDP), which was designed to accelerate atomistic structure prediction. The EDDP is simple and cost-efficient. It relies on training data generated in small unit cells and is fit using a lightweight neural network, leading to smooth interactions which exhibit the robust transferability essential for structure prediction. Here, we present a variety of applications of EDDPs, enabled by recent developments of the open-source EDDP software. New features include interfaces to phonon and molecular dynamics codes, as well as deployment of the ensemble deviation for estimating the confidence in EDDP predictions. Through case studies ranging from elemental carbon and lead to the binary scandium hydride and the ternary zinc cyanide, we demonstrate that EDDPs can be trained to cover wide ranges of pressures and stoichiometries, and used to evaluate phonons, phase diagrams, superionicity, and thermal expansion. These developments complement continued success in accelerated structure prediction. [ABSTRACT FROM AUTHOR]

Titanium dioxide has been extensively studied in the rutile or anatase phase, while its high-pressure phases are less well-understood, despite that many are thought to have interesting optical, mechanical, and electrochemical properties. First-principles methods, such as density functional theory (DFT), are often used to compute the enthalpies of TiO2 phases at 0 K, but they are expensive and, thus, impractical for long time scale and large system-size simulations at finite temperatures. On the other hand, cheap empirical potentials fail to capture the relative stabilities of various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO2, we design an empirical model as a baseline and then train a machine learning potential based on the difference between the DFT data and the empirical model. This so-called Δ-learning potential contains long-range electrostatic interactions and predicts the 0 K enthalpies of stable TiO2 phases that are in good agreement with DFT. We construct a pressure–temperature phase diagram of TiO2 in the range 0 < P < 70 GPa and 100 < T < 1500 K. We then simulate dynamic phase transition processes by compressing anatase at different temperatures. At 300 K, we predominantly observe an anatase-to-baddeleyite transformation at about 20 GPa via a martensitic two-step mechanism with a highly ordered and collective atomic motion. At 2000 K, anatase can transform into cotunnite around 45–55 GPa in a thermally activated and probabilistic manner, accompanied by diffusive movement of oxygen atoms. The pressures computed for these transitions show good agreement with experiments. Our results shed light on how to synthesize and stabilize high-pressure TiO2 phases, and our method is generally applicable to other functional materials with multiple polymorphs. [ABSTRACT FROM AUTHOR]

Through a series of high-pressure x-ray diffraction experiments combined with in situ laser heating, we explore the pressure–temperature phase diagram of germanium (Ge) at pressures up to 110 GPa and temperatures exceeding 3000 K. In the pressure range of 64–90 GPa, we observe orthorhombic Ge-IV transforming above 1500 K to a previously unobserved high-temperature phase, which we denote as Ge-VIII. This high-temperature phase is characterized by a tetragonal crystal structure, space group I4/mmm. Density functional theory simulations confirm that Ge-IV becomes unstable at high temperatures and that Ge-VIII is highly competitive and dynamically stable at these conditions. The existence of Ge-VIII has profound implications for the pressure–temperature phase diagram, with melting conditions increasing to much higher temperatures than previous extrapolations would imply. [ABSTRACT FROM AUTHOR]

Cathodes are critical components of rechargeable batteries. Conventionally, the search for cathode materials relies on experimental trial-and-error and a traversing of existing computational/experimental databases. While these methods have led to the discovery of several commercially viable cathode materials, the chemical space explored so far is limited and many phases will have been overlooked, in particular, those that are metastable. We describe a computational framework for battery cathode exploration based on ab initio random structure searching (AIRSS), an approach that samples local minima on the potential energy surface to identify new crystal structures. We show that by delimiting the search space using a number of constraints, including chemically aware minimum interatomic separations, cell volumes, and space group symmetries, AIRSS can efficiently predict both thermodynamically stable and metastable cathode materials. Specifically, we investigate LiCoO2, LiFePO4, and LixCuyFz to demonstrate the efficiency of the method by rediscovering the known crystal structures of these cathode materials. The effect of parameters, such as minimum separations and symmetries, on the efficiency of the sampling is discussed in detail. The adaptation of the minimum interatomic distances on a species-pair basis, from low-energy optimized structures to efficiently capture the local coordination environment of atoms, is explored. A family of novel cathode materials based on the transition-metal oxalates is proposed. They demonstrate superb energy density, oxygen-redox stability, and lithium diffusion properties. This article serves both as an introduction to the computational framework and as a guide to battery cathode material discovery using AIRSS. [ABSTRACT FROM AUTHOR]

Optical and synchrotron x-ray diffraction diamond anvil cell experiments have been combined with first-principles theoretical structure predictions to investigate mixtures of N2 and H2 up to 55 GPa. Our experiments show the formation of structurally complex van der Waals compounds [see also D. K. Spaulding et al., Nat. Commun. 5, 5739 (2014)] above 10 GPa. However, we found that these NxH (0.5 < x < 1.5) compounds transform abruptly to new oligomeric materials through barochemistry above 47 GPa and photochemistry at pressures as low as 10 GPa. These oligomeric compounds can be recovered to ambient pressure at T < 130 K, whereas at room temperature, they can be metastable on pressure release down to 3.5 GPa. Extensive theoretical calculations show that such oligomeric materials become thermodynamically more stable in comparison to mixtures of N2, H2, and NH3 above approximately 40 GPa. Our results suggest new pathways for synthesis of environmentally benign high energy-density materials. These materials could also exist as alternative planetary ices. [ABSTRACT FROM AUTHOR]

We study the effects of atomic vibrations on the solid-state chemical shielding tensor using first principles density functional theory calculations. At the harmonic level, we use a Monte Carlo method and a perturbative expansion. The Monte Carlo method is accurate but computationally expensive, while the perturbative method is computationally more efficient, but approximate. We find excellent agreement between the two methods for both the isotropic shift and the shielding anisotropy. The effects of zero-point quantum mechanical nuclear motion are important up to relatively high temperatures: at 500 K they still represent about half of the overall vibrational contribution. We also investigate the effects of anharmonic vibrations, finding that their contribution to the zero-point correction to the chemical shielding tensor is small. We exemplify these ideas using magnesium oxide and the molecular crystals L-alanine and β-aspartyl-L-alanine. We therefore propose as the method of choice to incorporate the effects of temperature in solid state chemical shielding tensor calculations using the perturbative expansion within the harmonic approximation. This approach is accurate and requires a computational effort that is about an order of magnitude smaller than that of dynamical or Monte Carlo approaches, so these effects might be routinely accounted for. [ABSTRACT FROM AUTHOR]

A combination of first-principles density functional theory calculations and a search over structures is used to predict the stability of a proton-transfer modification of ammonia monohydrate with space group P4/nmm. The phase diagram is calculated with the Perdew-Burke-Ernzerhof (PBE) density functional, and the effects of a semi-empirical dispersion correction, zero point motion, and finite temperature are investigated. Comparison with MP2 and coupled cluster calculations shows that the PBE functional over-stabilizes proton transfer phases because too much electronic charge moves with the proton. This over-binding is partially corrected by using the PBE0 hybrid exchange-correlation functional, which increases the enthalpy of P4/nmm by about 0.6 eV per formula unit relative to phase I of ammonia monohydrate and shifts the transition to the proton transfer phase from the PBE pressure of 2.8 GPa to about 10 GPa. This is consistent with experiment as proton transfer phases have not been observed at pressures up to ∼9 GPa, while higher pressures have not yet been explored experimentally. [ABSTRACT FROM AUTHOR]

NMR chemical shifts were calculated from first principles for well defined crystalline organic solids. These density functional theory calculations were carried out within the plane-wave pseudopotential framework, in which truly extended systems are implicitly considered. The influence of motional effects was assessed by averaging over vibrational modes or over snapshots taken from ab initio molecular dynamics simulations. It is observed that the zero-point correction to chemical shifts can be significant, and that thermal effects are particularly noticeable for shielding anisotropies and for a temperature-dependent chemical shift. This study provides insight into the development of highly accurate first principles calculations of chemical shifts in solids, highlighting the role of motional effects on well defined systems. [ABSTRACT FROM AUTHOR]

We have combined a computational search strategy with first-principles density-functional-theory calculations to identify metastable phases of H2O under pressure. The most stable structures consist of water molecules, while the most energetic metastable phases consist of oxygen and hydrogen molecules. In between lie many other metastable phases, consisting of various combinations of a few small molecules. It may be possible to synthesize some of these metastable phases, and we use our results to understand the nature of the crystalline metastable phase of H2O recently synthesized by Mao et al. [Science 314, 636 (2006)]. [ABSTRACT FROM AUTHOR]

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has also been applied to the silicophosphate, Si5O(PO4)6, giving 31P–29Si-couplings which are in excellent agreement with experiment. [ABSTRACT FROM AUTHOR]

The use of the spin-dependent pseudopotentials has been shown to markedly enhance the transferability of the commonly used spin-neutral pseudopotential method for the study of the structural and magnetic properties of transition-metal-containing materials. Unfortunately, because the method was based on the rather expensive norm-conserving pseudopotential formalism, the method was limited to the study of fairly small systems. Here we present an extension of the spin-dependent pseudopotential method for the far more computationally advantageous ultrasoft formalism and show that it is very easy to add such a feature to any preexisting computer code. We benchmark our new method by comparing to previously published results and then apply it to the study of several relevant test cases: bulk Ni, Fe, and Co, as well as a Pd atomic wire. [ABSTRACT FROM AUTHOR]

We present a method for the first principles density functional calculation of relativistic all-electron nuclear magnetic resonance chemical shifts using pseudopotentials. The method is based on the gauge including projector augmented wave approach of Pickard and Mauri [Phys. Rev. B 63, 245101 (2001)]. Relativistic effects are included at the level of the scalar-relativistic zeroth-order regular approximation. The method allows chemical shifts of large, low symmetry structures containing heavy elements to be calculated efficiently. We demonstrate its success for a range of Se and Te containing molecules. [ABSTRACT FROM AUTHOR]

We have used density-functional-theory (DFT) methods together with a structure searching algorithm to make an experimentally constrained prediction of the structure of ammonia dihydrate II (ADH-II). The DFT structure is in good agreement with neutron diffraction data and verifies the prediction. The structure consists of the same basic structural elements as ADH-I, with a modest alteration to the packing, but a considerable reduction in volume. The phase diagram of the known ADH and ammonia monohydrate + water-ice structures is calculated with the Perdew-Burke-Ernzerhof density functional, and the effects of a semi-empirical dispersion corrected functional are investigated. The results of our DFT calculations of the finite-pressure elastic constants of ADH-II are compared with the available experimental data for the elastic strain coefficients. [ABSTRACT FROM AUTHOR]

Differentiable programming has facilitated numerous methodological advances in scientific computing. Physics engines supporting automatic differentiation have simpler code, accelerating the development process and reducing the maintenance burden. Furthermore, fully differentiable simulation tools enable direct evaluation of challenging derivatives-including those directly related to properties measurable by experiment-that are conventionally computed with finite difference methods. Here, we investigate automatic differentiation in the context of orbital-free density functional theory (OFDFT) simulations of materials, introducing PROFESS-AD. Its automatic evaluation of properties derived from first derivatives, including functional potentials, forces, and stresses, facilitates the development and testing of new density functionals, while its direct evaluation of properties requiring higher-order derivatives, such as bulk moduli, elastic constants, and force constants, offers more concise implementations than conventional finite difference methods. For these reasons, PROFESS-AD serves as an excellent prototyping tool and provides new opportunities for OFDFT.

Optical and synchrotron x-ray diffraction diamond anvil cell experiments have been combined with first-principles theoretical structure predictions to investigate mixtures of N2 and H2 up to 55 GPa. Our experiments show the formation of structurally complex van der Waals compounds [see also D. K. Spaulding et al., Nat. Commun. 5, 5739 (2014)] above 10 GPa. However, we found that these NxH (0.5 < x < 1.5) compounds transform abruptly to new oligomeric materials through barochemistry above 47 GPa and photochemistry at pressures as low as 10 GPa. These oligomeric compounds can be recovered to ambient pressure at T < 130 K, whereas at room temperature, they can be metastable on pressure release down to 3.5 GPa. Extensive theoretical calculations show that such oligomeric materials become thermodynamically more stable in comparison to mixtures of N2, H2, and NH3 above approximately 40 GPa. Our results suggest new pathways for synthesis of environmentally benign high energy-density materials. These materials could also exist as alternative planetary ices.

We have combined a computational search strategy with first-principles density-functional-theory calculations to identify metastable phases of H(2)O under pressure. The most stable structures consist of water molecules, while the most energetic metastable phases consist of oxygen and hydrogen molecules. In between lie many other metastable phases, consisting of various combinations of a few small molecules. It may be possible to synthesize some of these metastable phases, and we use our results to understand the nature of the crystalline metastable phase of H(2)O recently synthesized by Mao et al. [Science 314, 636 (2006)].