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]

Fluorite-perovskite heterointerfaces garner great interest for enhanced ionic conductivity for application in electronic and energy devices. However, the origin of observed enhanced ionic conductivity as well as the details of the atomic structure at these interfaces remain elusive. Here, systematic, multi-stoichiometry computational searches and experimental investigations are performed to obtain stable and exact atomic structures of interfaces between CeO2 and SrTiO3—two archetypes of the corresponding structural families. Local reconstructions take place at the interface because of mismatched lattices. TiO2 terminated SrTiO3 causes a buckled rock salt CeO interface layer to emerge. In contrast, SrO terminated SrTiO3 maintains the fluorite structure at the interface compensated by a partially occupied anion lattice. Moderate enhancement in oxygen diffusion is found along the interface by simulations, yet evidence to support further significant enhancement is lacking. Our findings demonstrate the control of interface termination as an effective pathway to achieve desired device performance. [ABSTRACT FROM AUTHOR]

The recent claim of room temperature superconductivity in a copper-doped lead apatite compound, called LK-99, has sparked remarkable interest and controversy. Subsequent experiments have largely failed to reproduce the claimed superconductivity, while theoretical works have identified multiple key features including strong electronic correlation, structural instabilities, and dopability constraints. A puzzling claim of several recent theoretical studies is that both parent and copper-doped lead apatite structures are dynamically unstable at the harmonic level, questioning decades of experimental reports of the parent compound structures and the recently proposed copper-doped structures. In this work, we demonstrate that both parent and copper-doped lead apatite structures are dynamically stable at room temperature. Anharmonic phonon–phonon interactions play a key role in stabilizing some copper-doped phases, while most phases are largely stable even at the harmonic level. We also show that dynamical stability depends on both volume and correlation strength, suggesting controllable ways of exploring the copper-doped lead apatite structural phase diagram. Our results fully reconcile the theoretical description of the structures of both parent and copper-doped lead apatite with the experiment. [ABSTRACT FROM AUTHOR]

Using first-principles calculations and crystal structure search methods, we found that many covalently bonded molecules such as H2, N2, CO2, NH3, H2O and CH4 may react with NaCl, a prototype ionic solid, and form stable compounds under pressure while retaining their molecular structure. These molecules, despite whether they are homonuclear or heteronuclear, polar or non-polar, small or large, do not show strong chemical interactions with surrounding Na and Cl ions. In contrast, the most stable molecule among all examples, N2, is found to transform into cyclo-N5− anions while reacting with NaCl under high pressures. It provides a new route to synthesize pentazolates, which are promising green energy materials with high energy density. Our work demonstrates a unique and universal hybridization propensity of covalently bonded molecules and solid compounds under pressure. This surprising miscibility suggests possible mixing regions between the molecular and rock layers in the interiors of large planets. [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]

Nitrogen-doped lutetium hydride has recently been proposed as a near-ambient-conditions superconductor. Interestingly, the sample transforms from blue to pink to red as a function of pressure, but only the pink phase is claimed to be superconducting. Subsequent experimental studies have failed to reproduce the superconductivity, but have observed pressure-driven colour changes including blue, pink, red, violet, and orange. However, discrepancies exist among these experiments regarding the sequence and pressure at which these colour changes occur. Given the claimed relationship between colour and superconductivity, understanding colour changes in nitrogen-doped lutetium hydride may hold the key to clarifying the possible superconductivity in this compound. Here, we present a full microscopic theory of colour in lutetium hydride, revealing that hydrogen-deficient LuH2 is the only phase which exhibits colour changes under pressure consistent with experimental reports, with a sequence blue-violet-pink-red-orange. The concentration of hydrogen vacancies controls the precise sequence and pressure of colour changes, rationalising seemingly contradictory experiments. Nitrogen doping also modifies the colour of LuH2 but it plays a secondary role compared to hydrogen vacancies. Therefore, we propose hydrogen-deficient LuH2 as the key phase for exploring the superconductivity claim in the lutetium-hydrogen system. Finally, we find no phonon-mediated superconductivity near room temperature in the pink phase. Nitrogen-doped lutetium hydride, recently proposed as a superconductor at near-ambient conditions, features distinct color changes from blue to pink to red as a function of pressure. Using theoretical calculations, the authors identify the pink phase as hydrogen-deficient LuH2 and find that this phase is not a phonon-mediated superconductor near room temperature. Further, the color is controlled by the concentration of hydrogen vacancies. [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]

Motivated by the recent report of room-temperature superconductivity at near-ambient pressure in N-doped lutetium hydride, we performed a comprehensive, detailed study of the phase diagram of the Lu–N–H system, looking for superconducting phases. We combined ab initio crystal structure prediction with ephemeral data-derived interatomic potentials to sample over 200,000 different structures. Out of the more than 150 structures predicted to be metastable within ~50 meV from the convex hull we identify 52 viable candidates for conventional superconductivity, for which we computed their superconducting properties from Density Functional Perturbation Theory. Although for some of these structures we do predict a finite superconducting Tc, none is even remotely compatible with room-temperature superconductivity as reported by Dasenbrock et al. Our work joins the broader community effort that has followed the report of near-ambient superconductivity, confirming beyond reasonable doubt that no conventional mechanism can explain the reported Tc in Lu–N–H. Superconductivity was recently reported experimentally in nitrogen-doped lutetium hydride with Tc = 294 K at a pressure of 1 GPa. Here, via theoretical calculations, the authors find no structures capable of supporting conventional superconductivity in the Lu-N-H system at ambient pressure. [ABSTRACT FROM AUTHOR]

The discovery of 250-kelvin superconducting lanthanum polyhydride under high pressure marked a significant advance toward the realization of a room‐temperature superconductor. X-ray diffraction (XRD) studies reveal a nonstoichiometric LaH9.6 or LaH10±δ polyhydride responsible for the superconductivity, which in the literature is commonly treated as LaH10 without accounting for stoichiometric defects. Here, we discover significant nuclear quantum effects (NQE) in this polyhydride, and demonstrate that a minor amount of stoichiometric defects will cause quantum proton diffusion in the otherwise rigid lanthanum lattice in the ground state. The diffusion coefficient reaches ~10−7 cm2/s in LaH9.63 at 150 gigapascals and 240 kelvin, approaching the upper bound value of interstitial hydrides at comparable temperatures. A puzzling phenomenon observed in previous experiments, the positive pressure dependence of the superconducting critical temperature Tc below 150 gigapascals, is explained by a modulation of the electronic structure due to a premature distortion of the hydrogen lattice in this quantum fluxional structure upon decompression, and resulting changes of the electron-phonon coupling. This finding suggests the coexistence of the quantum proton fluxion and hydrogen-induced superconductivity in this lanthanum polyhydride, and leads to an understanding of the structural nature and superconductivity of nonstoichiomectric hydrogen-rich materials. The role of stoichiometric defects in the superconducting polyhydride LaH10±δ has received little attention so far. Here, the authors use molecular-dynamics simulations to show that a small amount of stoichiometric defects will cause quantum proton diffusion in the otherwise rigid lanthanum lattice. [ABSTRACT FROM AUTHOR]

Most metals adopt simple structures such as body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal close-packed (HCP) structures in specific groupings across the periodic table, and many undergo transitions to surprisingly complex structures on compression, not expected from conventional free-electron-based theories of metals. First-principles calculations have been able to reproduce many observed structures and transitions, but a unified, predictive theory that underlies this behavior is not yet in hand. Discovered by analyzing the electronic properties of metals in various lattices over a broad range of sizes and geometries, a remarkably simple theory shows that the stability of metal structures is governed by electrons occupying local interstitial orbitals and their strong chemical interactions. The theory provides a basis for understanding and predicting structures in solid compounds and alloys over a broad range of conditions. [ABSTRACT FROM AUTHOR]

KSiH3 exhibits 4.1 wt% experimental hydrogen storage capacity and shows reversibility under moderate conditions, which provides fresh impetus to the search for other complex hydrides in the K–Si–H system. Here, we reproduce the stable Fm3¯m phase of K2SiH6 and uncover two denser phases, space groups P3¯m1 and P63mc at ambient pressure, by means of first-principles structure searches. We note that P3¯m1-K2SiH6 has a high hydrogen content of 5.4 wt% and a volumetric density of 88.3 g L−1. Further calculations suggest a favorable dehydrogenation temperature Tdes of −20.1/55.8 °C with decomposition into KSi + K + H2. The higher hydrogen density and appropriate dehydrogenation temperature indicate that K2SiH6 is a promising hydrogen storage material, and our results provide helpful and clear guidance for further experimental studies. We found three further potential hydrogen storage materials stable at high pressure: K2SiH8, KSiH7 and KSiH8. These results suggest the need for further investigations into hydrogen storage materials among such ternary hydrides at high pressure. [ABSTRACT FROM AUTHOR]

The choice of cathode material in Li-ion batteries underpins their overall performance. Discovering new cathode materials is a slow process, and all major commercial cathode materials are still based on those identified in the 1990s. Discovery of materials using high-throughput calculations has attracted great research interest; however, reliance on databases of existing materials begs the question of whether these approaches are applicable for finding truly novel materials. In this work, we demonstrate that ab initio random structure searching (AIRSS), a first-principles structure prediction method that does not rely on any pre-existing data, can locate low energy structures of complex cathode materials efficiently based only on chemical composition. We use AIRSS to explore three Fe-containing polyanion compounds as low-cost cathodes. Using known quaternary LiFePO4 and quinary LiFeSO4F cathodes as examples, we easily reproduce the known polymorphs, in addition to predicting other, hitherto unknown, low energy polymorphs and even finding a new polymorph of LiFeSO4F that is more stable than the known ones. We then explore the phase space for Fe-containing fluoroxalates, predicting a range of redox-active phases that are yet to be experimentally synthesized, demonstrating the suitability of AIRSS as a tool for accelerating the discovery of novel cathode materials. [ABSTRACT FROM AUTHOR]

Sulfur and lanthanum hydrides under compression display superconducting states with high observed critical temperatures. It has been recently demonstrated that carbonaceous sulfur hydride displays room temperature superconductivity. However, this phenomenon has been observed only at very high pressure. Here, we theoretically search for superconductors with very high critical temperatures, but at much lower pressures. We describe two of such sodalite-type clathrate hydrides, YbH6 and LuH6. These hydrides are metastable and are predicted to superconduct with T c ∼ 145 K at 70 GPa and T c ∼ 273 K at 100 GPa, respectively. This striking result is a consequence of the strong interrelationship between the f states present at the Fermi level, structural stability, and the final T c value. For example, TmH6, with unfilled 4 f orbitals, is stable at 50 GPa, but has a relatively low value of T c of 25 K. The YbH6 and LuH6 compounds, with their filled f -shells, exhibit prominent phonon “softening”, which leads to a strong electron-phonon coupling, and as a result, an increase in T c. [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]

Using spin-assisted ab initio random structure searches, we explore an exhaustive quantum phase diagram of archetypal interfaced Mott insulators, i.e. lanthanum-iron and lanthanum-titanium oxides. In particular, we report that the charge transfer induced by the interfacial electronic reconstruction stabilises a high-spin ferrous Fe2+ state. We provide a pathway to control the strength of correlation in this electronic state by tuning the epitaxial strain, yielding a manifold of quantum electronic phases, i.e. Mott-Hubbard, charge transfer and Slater insulating states. Furthermore, we report that the electronic correlations are closely related to the structural oxygen octahedral rotations, whose control is able to stabilise the low-spin state of Fe2+ at low pressure previously observed only under the extreme high pressure conditions in the Earth's lower mantle. Thus, we provide avenues for magnetic switching via THz radiations which have crucial implications for next generation of spintronics technologies. [ABSTRACT FROM AUTHOR]

Double perovskite oxides, with generalized formula A2BB ′ O6, attract wide interest due to their multiferroic and charge transfer properties. They offer a wide range of potential applications such as spintronics and electrically tunable devices. However, great practical limitations are encountered, since a spontaneous order of the B-site cations is notoriously hard to achieve. In this joint experimental-theoretical work, we focused on the characterization of double perovskites La2TiFeO6 and La2VCuO6 films grown by pulsed laser deposition and interpretation of the observed B-site disorder and partial charge transfer between the B-site ions. A random structure sampling method was used to show that several phases compete due to their corresponding configurational entropy. In order to capture a representative picture of the most relevant competing microstates in realistic experimental conditions, this search included the potential formation of non-stoichiometric phases as well, which could also be directly related to the observed partial charge transfer. We optimized the information encapsulated in the potential energy landscape, captured via structure sampling, by evaluating both enthalpic and entropic terms. These terms were employed as a metric for the competition of different phases. This approach, applied herein specifically to La2TiFeO6, highlights the presence of highly entropic phases above the ground state which can explain the disorder observed frequently in the broader class of double perovskite oxides. [ABSTRACT FROM AUTHOR]

The solar system's outer planets, and many of their moons, are dominated by matter from the H-C-N-O chemical space, based on solar system abundances of hydrogen and the planetary ices H2O, CH4, and NH3. In the planetary interiors, these ices will experience extreme pressure conditions, around 5 Mbar at the Neptune mantle-core boundary, and it is expected that they undergo phase transitions, decompose, and form entirely new compounds. While temperature will dictate the formation of compounds, groundstate density functional theory allows us to probe the chemical effects resulting from pressure alone. These structural developments in turn determine the planets' interior structures, thermal evolution, and magnetic field generation, among others. Despite its importance, the H-C-N-O system has not been surveyed systematically to explore which compounds emerge at high-pressure conditions, and what governs their stability. Here, we report on and analyze an unbiased crystal structure search among H-C-N-O compounds between 1 and 5 Mbar. We demonstrate that simple chemical rules drive stability in this composition space, which explains why the simplest possible quaternary mixture HCNO--isoelectronic to diamond--emerges as a stable compound and discuss dominant decomposition products of planetary ice mixtures. [ABSTRACT FROM AUTHOR]

The Na–P and Li–P chemical systems were studied under pressure using synchrotron X-ray diffraction in a diamond anvil cell up to 20 GPa, combined with the AIRSS ab initio random structure searching technique. The results reveal an enhanced reactivity of both alkali metals with phosphorous at slightly elevated pressures. This enables the synthesis of Li3P and Na3P at room temperature (RT) starting from element precursors, bypassing the established chemical synthesis methods. Both compounds undergo a pressure-induced phase transition from the hexagonal Na3As-type structure (stable at ambient conditions) towards a Fm3¯m (FCC) structure that remains stable up to 20 GPa. Attempts to synthesize compounds with higher alkali metal content (such as Li5P) using high-temperature and -pressure conditions (up to 2000+ K and 30 GPa), inspired by recent theoretical predictions, were not successful. [ABSTRACT FROM AUTHOR]

A combined experimental and computational approach was used to distinguish between different polymorphs of the pharmaceutical drug aspirin. This method involves the use of abinitio random structure searching (AIRSS), a density functional theory (DFT)‐based crystal structure prediction method for the high‐accuracy prediction of polymorphic structures, with DFT calculations of nuclear magnetic resonance (NMR) parameters and solid‐state NMR experiments at natural abundance. AIRSS was used to predict the crystal structures of form‐I and form‐II of aspirin. The root‐mean‐square deviation between experimental and calculated 1H chemical shifts was used to identify form‐I as the polymorph present in the experimental sample, the selection being successful despite the large similarities between the molecular environments in the crystals of the two polymorphs. [ABSTRACT FROM AUTHOR]

Helium and methane are major components of giant icy planets and are abundant in the universe. However, helium is the most inert element in the periodic table and methane is one of the most hydrophobic molecules, thus whether they can react with each other is of fundamental importance. Here, our crystal structure searches and first-principles calculations predict that a He3CH4 compound is stable over a wide range of pressures from 55 to 155 GPa and a HeCH4 compound becomes stable around 105 GPa. As nice examples of pure van der Waals crystals, the insertion of helium atoms changes the original packing of pure methane molecules and also largely hinders the polymerization of methane at higher pressures. After analyzing the diffusive properties during the melting of He3CH4 at high pressure and high temperature, in addition to a plastic methane phase, we have discovered an unusual phase which exhibits coexistence of diffusive helium and plastic methane. In addition, the range of the diffusive behavior within the helium-methane phase diagram is found to be much narrower compared to that of previously predicted helium-water compounds. This may be due to the weaker van der Waals interactions between methane molecules compared to those in helium-water compounds, and that the helium-methane compound melts more easily. [ABSTRACT FROM AUTHOR]

Hydrogen, the simplest and most abundant element in the Universe, develops a remarkably complex behaviour upon compression1. Since Wigner predicted the dissociation and metallization of solid hydrogen at megabar pressures almost a century ago2, several efforts have been made to explain the many unusual properties of dense hydrogen, including a rich and poorly understood solid polymorphism1,3–5, an anomalous melting line6 and the possible transition to a superconducting state7. Experiments at such extreme conditions are challenging and often lead to hard-to-interpret and controversial observations, whereas theoretical investigations are constrained by the huge computational cost of sufficiently accurate quantum mechanical calculations. Here we present a theoretical study of the phase diagram of dense hydrogen that uses machine learning to 'learn' potential-energy surfaces and interatomic forces from reference calculations and then predict them at low computational cost, overcoming length- and timescale limitations. We reproduce both the re-entrant melting behaviour and the polymorphism of the solid phase. Simulations using our machine-learning-based potentials provide evidence for a continuous molecular-to-atomic transition in the liquid, with no first-order transition observed above the melting line. This suggests a smooth transition between insulating and metallic layers in giant gas planets, and reconciles existing discrepancies between experiments as a manifestation of supercritical behaviour. Simulations using machine-learning-based interatomic potentials in dense hydrogen overcome system size and timescale limitations, providing evidence of a supercritical behaviour of high-pressure liquid hydrogen and reconciling theoretical and experimental discrepancies. [ABSTRACT FROM AUTHOR]

The discovery of materials is increasingly guided by quantum‐mechanical crystal‐structure prediction, but the structural complexity in bulk and nanoscale materials remains a bottleneck. Here we demonstrate how data‐driven approaches can vastly accelerate the search for complex structures, combining a machine‐learning (ML) model for the potential‐energy surface with efficient, fragment‐based searching. We use the characteristic building units observed in Hittorf's and fibrous phosphorus to seed stochastic ("random") structure searches over hundreds of thousands of runs. Our study identifies a family of hierarchically structured allotropes based on a P8 cage as principal building unit, including one‐dimensional (1D) single and double helix structures, nanowires, and two‐dimensional (2D) phosphorene allotropes with square‐lattice and kagome topologies. These findings yield new insight into the intriguingly diverse structural chemistry of phosphorus, and they provide an example for how ML methods may, in the long run, be expected to accelerate the discovery of hierarchical nanostructures. [ABSTRACT FROM AUTHOR]

The discovery of materials is increasingly guided by quantum‐mechanical crystal‐structure prediction, but the structural complexity in bulk and nanoscale materials remains a bottleneck. Here we demonstrate how data‐driven approaches can vastly accelerate the search for complex structures, combining a machine‐learning (ML) model for the potential‐energy surface with efficient, fragment‐based searching. We use the characteristic building units observed in Hittorf's and fibrous phosphorus to seed stochastic ("random") structure searches over hundreds of thousands of runs. Our study identifies a family of hierarchically structured allotropes based on a P8 cage as principal building unit, including one‐dimensional (1D) single and double helix structures, nanowires, and two‐dimensional (2D) phosphorene allotropes with square‐lattice and kagome topologies. These findings yield new insight into the intriguingly diverse structural chemistry of phosphorus, and they provide an example for how ML methods may, in the long run, be expected to accelerate the discovery of hierarchical nanostructures. [ABSTRACT FROM AUTHOR]

Nuclear Magnetic Resonance (NMR) is one of the most powerful experimental techniques to characterize the structure of molecules and confined liquids. Nevertheless, the complexity of the systems under investigation usually requires complementary computational studies to interpret the NMR results. In this work we focus on polycyclic aromatic hydrocarbons (PAHs), an important class of organic molecules which have been commonly used as simple analogues for the spectroscopic properties of more complex systems, such as porous disordered carbons. We use Density Functional Theory (DFT) to calculate 13C chemical shifts and Nucleus Independent Chemical Shifts (NICS) for 34 PAHs. The results show a clear molecular size dependence of the two quantities, as well as the convergence of the 13C NMR shifts towards the values observed for graphene. We then present two computationally cheap models for the prediction of NICS in simple PAHs. We show that while a simple dipolar model fails to produce accurate values, a perturbative tight-binding approach can be successfully applied for the prediction of NICS in this series of molecules, including some non-planar ones containing 5- and 7-membered rings. This model, one to two orders of magnitude faster than DFT calculations, is very promising and can be further refined in order to study more complex systems. [ABSTRACT FROM AUTHOR]

Predicting phase stabilities of crystal polymorphs is central to computational materials science and chemistry. Such predictions are challenging because they first require searching for potential energy minima and then performing arduous free-energy calculations to account for entropic effects at finite temperatures. Here, we develop a framework that facilitates such predictions by exploiting all the information obtained from random searches of crystal structures. This framework combines automated clustering, classification and visualisation of crystal structures with machine-learning estimation of their enthalpy and entropy. We demonstrate the framework on the technologically important system of TiO2, which has many polymorphs, without relying on prior knowledge of known phases. We find a number of new phases and predict the phase diagram and metastabilities of crystal polymorphs at 1600 K, benchmarking the results against full free-energy calculations. [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]

Vertically aligned nanocomposite (VAN) films have self-assembled pillar-matrix nanostructures. Owing to their large area-to-volume ratios, interfaces in VAN films are expected to play key roles in inducing functional properties, but our understanding is hindered by limited knowledge about their structures. Motivated by the lack of definitive explanation for the experimentally found enhanced ionic conductivity in Sm-doped-CeO2/SrTiO3 VAN films, we determine the structure at vertical interfaces using random structure searching and explore how it can affect ionic conduction. Interatomic potentials are used to perform the initial searching, followed by first-principles calculations for refinement. Previously unknown structures are found, with lower energy than that of an optimized hand-built model. We find a strongly distorted oxygen sublattice which gives a complex landscape of vacancy energies. The cation lattice remains similar to the bulk phase, but has a localized strain field. The excess energy of the interface is similar to that of high angle grain boundaries in SrTiO3. [ABSTRACT FROM AUTHOR]

We present computational discoveries of new structural phases of the B2O compound exhibiting novel bonding networks and electronic states at ambient and elevated pressures. Our advanced crystal structure searches in conjunction with density functional theory calculations have identified an orthorhombic phase of B2O that is energetically stable at ambient pressure and contains an intriguing bonding network of icosahedral B12 clusters bridged by oxygen atoms. As pressure increases above 1.9 GPa, a structural transformation takes the orthorhombic B2O into a pseudo-layered trigonal phase. We have performed extensive studies to investigate the evolution of chemical bonds and electronic states associated with the B12 icosahedral unit in the orthorhombic phase and the covalent B–O bonds in the trigonal phase. We have also examined the nature of the charge carriers and their coupling to the lattice vibrations in the newly identified B2O crystals. Interestingly, our results indicate that both B2O phases become superconducting at low temperatures, with transition temperatures of 6.4 K and 5.9 K, respectively, in the ambient and high-pressure phase. The present findings establish new B2O phases and characterize their structural and electronic properties, which offer insights and guidance for exploration toward further fundamental understanding and potential synthesis and application. [ABSTRACT FROM AUTHOR]

Abstract: The application of pressure allows systematic tuning of the charge density of a material cleanly, that is, without changes to the chemical composition via dopants, and exploratory high‐pressure experiments can inform the design of bulk syntheses of materials that benefit from their properties under compression. The electronic and structural response of semiconducting tin nitride Sn3N4 under compression is now reported. A continuous opening of the optical band gap was observed from 1.3 eV to 3.0 eV over a range of 100 GPa, a 540 nm blue‐shift spanning the entire visible spectrum. The pressure‐mediated band gap opening is general to this material across numerous high‐density polymorphs, implicating the predominant ionic bonding in the material as the cause. The rate of decompression to ambient conditions permits access to recoverable metastable states with varying band gaps energies, opening the possibility of pressure‐tuneable electronic properties for future applications. [ABSTRACT FROM AUTHOR]

Abstract: The application of pressure allows systematic tuning of the charge density of a material cleanly, that is, without changes to the chemical composition via dopants, and exploratory high‐pressure experiments can inform the design of bulk syntheses of materials that benefit from their properties under compression. The electronic and structural response of semiconducting tin nitride Sn3N4 under compression is now reported. A continuous opening of the optical band gap was observed from 1.3 eV to 3.0 eV over a range of 100 GPa, a 540 nm blue‐shift spanning the entire visible spectrum. The pressure‐mediated band gap opening is general to this material across numerous high‐density polymorphs, implicating the predominant ionic bonding in the material as the cause. The rate of decompression to ambient conditions permits access to recoverable metastable states with varying band gaps energies, opening the possibility of pressure‐tuneable electronic properties for future applications. [ABSTRACT FROM AUTHOR]

Crystal structure prediction algorithms, including ab initio random structure searching (AIRSS), are intrinsically limited by the huge computational cost of the underlying quantum-mechanical methods. We have recently shown that a novel class of machine learning (ML) based interatomic potentials can provide a way out: by performing a high-dimensional fit to the ab initio energy landscape, these potentials reach comparable accuracy but are orders of magnitude faster. In this paper, we develop our approach, dubbed Gaussian approximation potential-based random structure searching (GAP-RSS), towards a more general tool for exploring configuration spaces and predicting structures. We present a GAP-RSS interatomic potential model for elemental phosphorus, which identifies and correctly “learns” the orthorhombic black phosphorus (A17) structure without prior knowledge of any crystalline allotropes. Using the tubular structure of fibrous phosphorus as an example, we then discuss the limits of free searching, and discuss a possible way forward that combines a recently proposed fragment analysis with GAP-RSS. Examples of possible tubular (1D) and extended (3D) hypothetical allotropes of phosphorus as found by GAP-RSS are discussed. We believe that in the future, ML potentials could become versatile and routine computational tools for materials discovery and design. [ABSTRACT FROM AUTHOR]

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated. [ABSTRACT FROM AUTHOR]

We present a systematic study of the stability of nineteen different periodic structures using the finite range Lennard-Jones potential model discussing the effects of pressure, potential truncation, cutoff distance and Lennard-Jones exponents. The structures considered are the hexagonal close packed (hcp), face centred cubic (fcc) and seventeen other polytype stacking sequences, such as dhcp and 9R. We found that at certain pressure and cutoff distance values, neither fcc nor hcp is the ground state structure as previously documented, but different polytypic sequences. This behaviour shows a strong dependence on the way the tail of the potential is truncated. [ABSTRACT FROM AUTHOR]

Elemental phosphorus displays an impressive number of allotropes with highly diverse chemical and physical properties. White phosphorus has now been filled into single-wall carbon nanotubes (SWCNTs) from the liquid and thereby stabilized against the highly exothermic reaction with atmospheric oxygen. The encapsulated tetraphosphorus molecules were visualized with transmission electron microscopy, but found to convert readily into chain structures inside the SWCNT 'nanoreactors'. The energies of the possible chain structures were determined computationally, highlighting a delicate balance between the extent of polymerization and the SWCNT diameter. Experimentally, a single-stranded zig-zag chain of phosphorus atoms was observed, which is the lowest energy structure at small confinement diameters. These one-dimensional chains provide a glimpse into the very first steps of the transformation from white to red phosphorus. [ABSTRACT FROM AUTHOR]

Elemental phosphorus displays an impressive number of allotropes with highly diverse chemical and physical properties. White phosphorus has now been filled into single-wall carbon nanotubes (SWCNTs) from the liquid and thereby stabilized against the highly exothermic reaction with atmospheric oxygen. The encapsulated tetraphosphorus molecules were visualized with transmission electron microscopy, but found to convert readily into chain structures inside the SWCNT 'nanoreactors'. The energies of the possible chain structures were determined computationally, highlighting a delicate balance between the extent of polymerization and the SWCNT diameter. Experimentally, a single-stranded zig-zag chain of phosphorus atoms was observed, which is the lowest energy structure at small confinement diameters. These one-dimensional chains provide a glimpse into the very first steps of the transformation from white to red phosphorus. [ABSTRACT FROM AUTHOR]

Solid solution perovskite Pb(Zr1−xTix)O3 (PZT) is an industrially important material. Despite the long history of experimental and theoretical studies, the structure of this material is still under intensive discussion. In this work, we have applied structure searching coupled with density functional theory methods to provide a multiphase description of this material at x = 0.4. We demonstrate that the permutational freedom of B-site cations leads to the stabilisation of a variety of local phases reflecting a relatively flat energy landscape of PZT. Using a set of predicted local phases we reproduce the experimental pair distribution function (PDF) profile with high accuracy. We introduce a complex multiphase picture of the structure of PZT and show that additional monoclinic and rhombohedral phases account for a better description of the experimental PDF profile. We propose that such a multiphase picture reflects the entropy reached in the sample during the preparation process. [ABSTRACT FROM AUTHOR]

Materials informatics owes much to bioinformatics and the Materials Genome Initiative has been inspired by the Human Genome Project. But there is more to bioinformatics than genomes, and the same is true for materials informatics. Here we describe the rapidly expanding role of searching for structures of materials using firstprinciples electronic-structure methods. Structure searching has played an important part in unraveling structures of dense hydrogen and in identifying the record-hightemperature superconducting component in hydrogen sulfide at high pressures. We suggest that first-principles structure searching has already demonstrated its ability to determine structures of a wide range of materials and that it will play a central and increasing part in materials discovery and design. [ABSTRACT FROM AUTHOR]

The structural chemistry of materials containing low levels of nonstoichiometric hydrogen is difficult to determine, and producing structural models is challenging where hydrogen has no fixed crystallographic site. Here we demonstrate a computational approach employing ab initio random structure searching (AIRSS) to generate a series of candidate structures for hydrous wadsleyite (β-Mg2SiO4 with 1.6 wt% H2O), a high-pressure mineral proposed as a repository for water in the Earth's transition zone. Aligning with previous experimental work, we solely consider models with Mg3 (over Mg1, Mg2 or Si) vacancies. We adapt the AIRSS method by starting with anhydrous wadsleyite, removing a single Mg2+ and randomly placing two H+ in a unit cell model, generating 819 candidate structures. 103 geometries were then subjected to more accurate optimisation under periodic DFT. Using this approach, we find the most favourable hydration mechanism involves protonation of two O1 sites around the Mg3 vacancy. The formation of silanol groups on O3 or O4 sites (with loss of stable O1–H hydroxyls) coincides with an increase in total enthalpy. Importantly, the approach we employ allows observables such as NMR parameters to be computed for each structure. We consider hydrous wadsleyite (∼1.6 wt%) to be dominated by protonated O1 sites, with O3/O4–H silanol groups present as defects, a model that maps well onto experimental studies at higher levels of hydration (J. M. Griffin et al., Chem. Sci., 2013, 4, 1523). The AIRSS approach adopted herein provides the crucial link between atomic-scale structure and experimental studies. [ABSTRACT FROM AUTHOR]

Selenium-rich Ge–Te–Se glasses have been synthesized along the GeSe4–GeTe4 pseudo-composition line and acquired by 77Se Hahn echo magic-angle spinning NMR. The comparison with the GeSe4 spectrum shows a drastic modification of the typical double-resonance lineshape even at low Te concentrations (＜10%). In order to rationalize this feature and to understand the effect of Te on the structure of our glasses, first-principles molecular dynamics simulations and gauge including projector augmented wave NMR parameter calculations have been performed. The distribution of the tellurium atoms in the selenium phase was shown to be mainly responsible for the 77Se lineshape changes. Another possible factor related to the perturbation of the δiso value due to Te proximity appears to be much more limited in the bulk, while the results obtained using molecular models suggest shifts of several hundreds of ppm. [ABSTRACT FROM AUTHOR]