Search

The journal of physical chemistry letters 13(36), 8447 - 8454 (2022). doi:10.1021/acs.jpclett.2c02157, Through a series of high pressure diamond anvil experiments, we report the synthesis of alkaline earth (Ca, Sr, Ba) tetrahydrides, and investigate their properties through Raman spectroscopy, X-ray diffraction, and density functional theory calculations. The tetrahydrides incorporate both atomic and quasi-molecular hydrogen, and we find that the frequency of the intramolecular stretching mode of the $H^{𝛿−}_2$ units downshifts from Ca to Sr and to Ba upon compression. The experimental results indicate that the larger the host cation, the longer the $H^{𝛿−}_2$ bond. Analysis of the electron localization function (ELF) demonstrates that the lengthening of the H–H bond is caused by the charge transfer from the metal to $H^{𝛿−}_2$ and by the steric effect of the metal host on the $H–H$ bond. This effect is most prominent for BaH$_4$, where the precompression of $H^{𝛿−}_2$ units at 50 GPa results in bond lengths comparable to that of pure $H_2$ above 275 GPa., Published by ACS, Washington, DC

Acknowledgements: H.W. is thankful to Y. Ma and M. Shiga for valuable discussions, to V. Minkov and M. Eremets for the XRD data, and to C.Wang for the EPC data. The project is supported by the National Natural Science Foundation of China (Grant No. 11974135, 11874176, 12174170, and 12074138), the Natural Sciences and Engineering Research Council of Canada, the EPSRC through grants EP/P022596/1, and EP/S021981/1, and the startup funds of the office of the Dean of SASN of Rutgers University-Newark. P. T. S. thanks the Department of Materials Science and Metallurgy at the University of Cambridge for generous funding. The work of P. T. S. is further supported through a Trinity Hall research studentship. I. E. acknowledges financial support by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 802533). We used the computing facilities at Beijing Super Cloud Computing Center., 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.

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., 20 pages, 15 figures

Funder: Korea Institute of Science and Technology (KIST); doi: https://doi.org/10.13039/501100003693, Funder: Gouvernement du Canada | Natural Sciences and Engineering Research Council of Canada (Conseil de Recherches en Sciences Naturelles et en Génie du Canada); doi: https://doi.org/10.13039/501100000038, Funder: RCUK | Science and Technology Facilities Council (STFC); doi: https://doi.org/10.13039/501100000271, Carbon dioxide capture is essential to achieve net-zero emissions. A hurdle to the design of improved capture materials is the lack of adequate tools to characterise how CO2 adsorbs. Solid-state nuclear magnetic resonance (NMR) spectroscopy is a promising probe of CO2 capture, but it remains challenging to distinguish different adsorption products. Here we perform a comprehensive computational investigation of 22 amine-functionalised metal-organic frameworks and discover that 17O NMR is a powerful probe of CO2 capture chemistry that provides excellent differentiation of ammonium carbamate and carbamic acid species. The computational findings are supported by 17O NMR experiments on a series of CO2-loaded frameworks that clearly identify ammonium carbamate chain formation and provide evidence for a mixed carbamic acid – ammonium carbamate adsorption mode. We further find that carbamic acid formation is more prevalent in this materials class than previously believed. Finally, we show that our methods are readily applicable to other adsorbents, and find support for ammonium carbamate formation in amine-grafted silicas. Our work paves the way for investigations of carbon capture chemistry that can enable materials design.

In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies)., Main text: 23 pages, 4 figures. Supplementary: 68 pages

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.

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 superconducting. Subsequent experimental studies have failed to reproduce the superconductivity, but have confirmed the existence of pressure-driven colour changes. However, these colour changes appear in different sequences and at different pressures depending on the experiment, with observed colours including blue, pink, red, violet, and orange. Given the 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 describe a full microscopic theory of colour in lutetium hydride. We find that hydrogen-deficient LuH$_2$ is the only phase which exhibits colour changes under pressure consistent with experimental reports, with a sequence blue-violet-pink-red-orange. We also find that 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 LuH$_2$ but it plays a secondary role compared to hydrogen vacancies. Therefore, we propose hydrogen-deficient LuH$_2$ 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., Comment: 3 figures

Controlling the formation and stoichiometric content of desired phases of materials has become a central interest for the study of a variety of fields, notably high temperature superconductivity under extreme pressures. The further possibility of accessing metastable states by initiating reactions by x-ray triggered mechanisms over ultra-short timescales is enabled with the development of x-ray free electron lasers (XFEL). Utilizing the exceptionally high brilliance x-ray pulses from the EuXFEL, we report the synthesis of a previously unobserved yttrium hydride under high pressure, along with non-stoichiometric changes in hydrogen content as probed at a repetition rate of 4.5\,MHz using time-resolved x-ray diffraction. Exploiting non-equilibrium pathways we synthesize and characterize a hydride with yttrium cations in an \textit{A}15 structure type at 125\,GPa, predicted using crystal structure searches, with a hydrogen content between 4.0--5.75 hydrogens per cation, that is enthalpically metastable on the convex hull. We demonstrate a tailored approach to changing hydrogen content using changes in x-ray fluence that is not accessible using conventional synthesis methods, and reveals a new paradigm in metastable chemical physics.

Dataset of 10.48550/arXiv.2304.04447, PPF and LTFE gratefully acknowledge the S ̃ao Paulo Research Foundation (FAPESP) under Grants 2020/08258-0 and 2021/13441-1. CH acknowledges the Austrian Science Fund (FWF) Project No. P 32144-N36. LB acknowledges support from Project PE0000021,"Network 4 Energy Sustainable Transition – NEST", funded by the European Union NextGenerationEU, under the National Recovery and Resilience Plan (NRRP), Mission 4 Component 2 Investment 1.3 - Call for tender No. 1561 of 11.10.2022 of Ministero dell'Universit ́a e della Ricerca (MUR). Calculations were performed on the Vienna Scientific Cluster (proj. 71754 "TEST"), on the dcluster of the Graz University of Technology, on the CINECA cluster (proj. IsC99ACME-C); Computational resources for the UK partners were provided by the UK's National Supercomputer Service through the UK Car-Parrinello consortium (EP/P022561/1) and the Cambridge Service for Data Driven Discovery (CSD3) using Tier-2 EPSRC funding (EP/T022159/1).

Density functional theory in the local or semi-local density approximation is a powerful tool for materials simulation, yet it struggles in many cases to describe collective electronic order that is driven by electronic interactions. In this work it is shown how arbitrary, fluctuation-driven electronic order may be introduced into density functional theory using the quantum order-by-disorder framework. This is a method of calculating the free energy correction due to collective spin and charge fluctuations about a state that hosts static order, in a self-consistent manner. In practical terms, the quantum order-by-disorder method is applied to the Kohn-Sham auxiliary system of density functional theory to give an order-dependent correction to the exchange-correlation functional. Calculation of fluctuation propagators within density functional theory renders the result fully first-principles. Two types of order are considered as examples -- fluctuation-driven superconductivity and spin nematic order -- and implementation schemes are presented in each case., 20 pages, 2 figures

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.

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.

Designing materials with advanced functionalities is the main focus of contemporary solid-state physics and chemistry. Research efforts worldwide are funneled into a few high-end goals, one of the oldest, and most fascinating of which is the search for an ambient temperature superconductor (A-SC). The reason is clear: superconductivity at ambient conditions implies being able to handle, measure and access a single, coherent, macroscopic quantum mechanical state without the limitations associated with cryogenics and pressurization. This would not only open exciting avenues for fundamental research, but also pave the road for a wide range of technological applications, affecting strategic areas such as energy conservation and climate change. In this roadmap we have collected contributions from many of the main actors working on superconductivity, and asked them to share their personal viewpoint on the field. The hope is that this article will serve not only as an instantaneous picture of the status of research, but also as a true roadmap defining the main long-term theoretical and experimental challenges that lie ahead. Interestingly, although the current research in superconductor design is dominated by conventional (phonon-mediated) superconductors, there seems to be a widespread consensus that achieving A-SC may require different pairing mechanisms. In memoriam, to Neil Ashcroft, who inspired us all. © 2022 The Author(s). Published by IOP Publishing Ltd. Yundi Quan has provided many useful discussions on this topic, and assistance with preparation of the figures in this paper. Giustino’s review [12] contains a wealth of references on this topic. This work was supported by US National Science Foundation Grant DMR 1607139. 3. We acknowledge Ashkan Salamat, Elliot Snider, Nathan Dasenbrock-Gammon, Raymond McBride, Mathew Debessai for useful discussions. This was supported by NSF, Grant No. DMR-1809649, and by the DOE Stockpile Stewardship Academic Alliance Program, Grant No. DE-NA0003898.5. CH acknowledges support from the Austrian Science Fund (FWF) Project No. P 32144-N36. 6. This work was supported principally by the US National Science Foundation (DMR-1933622), and by the US Department of Energy (DE-SC0020340 and DE-NA0003975). 7. I would like to thank all my collaborators and friends that contributed to this work: D M Ceperley, M Holzmann, V Gorelov, G Rillo, Y Yubo, M A Morales. This work was supported by ANR-France under the program ‘Accueil de Chercheurs de Haut Niveau 2015’ project: HyLightExtreme. 8. MA acknowledges support from the Swiss National Science Foundation (Project P4P4P2-180669). 9. This work was supported by the National Natural Science Foundation of China (Grant Nos. 52090024 and 12074138), the Science Challenge Project (Grant No. TZ2016001), the Fundamental Research Funds for the Central Universities (Jilin University, JLU), the Program for JLU Science and Technology Innovative Research Team (JLUSTIRT), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB33000000) 10. ANK acknowledges the support through the NSF Award No. DMR-1821815. 11. We thank the Russian Science Foundation (Grant 19-72-30043) for support. 12. I thank Alice Shipley and Michael Hutcheon for their careful reading of the manuscript and insightful comments. This work has been funded by the EPSRC over many years (Projects EP/G007489/1 and EP/P022596/1), and through a Royal Society Wolfson Research Merit Award. 13. We acknowledge the National Science Foundation (DMR-1827815) for financial support. 14. The author acknowledges the financial support by JSPS KAK-ENHI Grant No. 19H05825. 15. This research was supported by the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Program (Grant Agreement No. 802533). 16. We acknowledge financial support by the European Research Council Advanced Grant FACT (ERC-2017-AdG-788890). AS acknowledges hospitality of the Physics Department of La Sapienza, under the program ‘Professori Visitatori 2020’. 17. ERM acknowledges support from the National Science Foundation (Award No. OAC-1740263). 18. The work was supported by the US Department of Energy Basic Energy Sciences under Contract No. DE-SC-0020385. 19. We thank G Profeta and G Lamura for fruitful scientific discussions. We acknowledge funding from the MIUR PRIN-2017 program (Grant No. 2017Z8TS5B—‘Tuning and understanding Quantum phases in 2D materials—Quantum2D’). 20. We thank G Profeta and G Lamura for fruitful scientific discussions. We acknowledge funding from the MIUR PRIN-2017program (Grant No. 2017Z8TS5B—‘Tuning and understanding Quantum phases in 2D materials—Quantum2D’). 21. We acknowledge financial support by the Deutsche Forschungs-Gemeinschaft through Grant (DFG) for funding through TRR 288–422213477 (B05). This research was supported in part by the National Science Foundation under Grant No. NSF PHY-1748958.

Superionic states simultaneously exhibit properties of a fluid and a solid. Proton (H+) superionicity in ice, H3O, He-H2O, and He-NH3 compounds is well-studied. However, hydride (H-) superionicity in H-rich compounds is rare, being associated with instability and strongly reducing conditions. Silicon, sodium, and hydrogen are abundant elements in many astrophysical bodies. Here, we use first-principles calculations to show that, at high pressure, Na, Si, and H can form several hypervalent compounds. A previously unreported superionic state of Na2SiH6 results from unconstrained H- in the hypervalent [SiH6]2- unit. Na2SiH6 is dynamically stable at low pressure (3 GPa), becoming superionic at 5 GPa, and re-entering solid/fluid states at about 25 GPa. Our observation of H- transport opens up a new field of H- conductors. It also has implications for the formation of conducting layers at depth in exotic carbon exoplanets, potentially enhancing the habitability of such planets.

Structure prediction has become a key task of the modern atomistic sciences, and depends on the rapid and reliable computation of the energy landscape. First principles density functional based calculations are highly reliable, faithfully describing the entire energy landscape. They are, however, computationally intensive and slow compared to interatomic potentials. Great progress has been made in the development of machine learning, or data derived, potentials, which promise to describe the entire energy landscape at first principles quality. However, compared to first principles approaches, their preparation can be time consuming and delay searching. Ab initio random structure searching (AIRSS) is a straightforward and powerful approach to structure prediction, based on the stochastic generation of sensible initial structures, and their repeated local optimisation. Here, a scheme, compatible with AIRSS, for the rapid construction of disposable, or ephemeral, data derived potentials (EDDPs) is described. These potentials are constructed using a homogeneous, separable manybody environment vector, and iterative neural network fits, sparsely combined through non-negative least squares. The approach is first tested on methane, boron nitride, elemental boron and urea. In the case of boron, an EDDP generated using data from small unit cells is used to rediscover the complex $\gamma$-boron structure without recourse to symmetry or fragments. Finally, an EDDP generated for silane (SiH$_4$) at 500 GPa enables the discovery of an extremely complex, dense, structure which significantly modifies silane's high pressure phase diagram. This has implications for the theoretical exploration for high temperature superconductivity in the dense hydrides, which have so far largely depended on searches in smaller unit cells., Comment: 17 pages, 12 figures

By combining pressures up to 50 GPa and temperatures of 1200 K, we synthesize the novel barium hydride, Ba8H46, stable down to 27 GPa. We use Raman spectroscopy, X-ray diffraction, and first-principles calculations to determine that this compound adopts a highly symmetric Pm3¯n structure with an unusual 534:1 hydrogen-to-barium ratio. This singular stoichiometry corresponds to the well-defined type-I clathrate geometry. This clathrate consists of a Weaire-Phelan hydrogen structure with the barium atoms forming a topologically close-packed phase. In particular, the structure is formed by H20 and H24 clathrate cages showing substantially weakened H-H interactions. Density functional theory (DFT) demonstrates that cubic Pm3¯n Ba8H46 requires dynamical effects to stabilize the H20 and H24 clathrate cages.

Most elemental metals under ambient conditions adopt simple structures such as BCC, FCC and HCP in specific groupings across the Periodic Table, and on compression, many of these elements undergo transitions to surprisingly complex structures, including open and low-symmetry phases 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 of bonding that underlies this behavior is not yet in hand. We propose a remarkably simple theory based on large-scale high-throughput calculations of the elements over a broad range of thermodynamic conditions. The results broaden the conventional concept of metallic bonding with a new perspective that both explains the stability of different simple structures as well as lower symmetry phases arising from electron localization at interstitial sites in these structures. The success of this simple framework reveals the important role of chemical interactions in governing the structures and electronic properties of simple metals.

Helium is the second most abundant element in the universe, and together with silica, they are major components of giant planets. Exploring the reactivity and state of helium and silica under high pressure is of fundamental importance for developing and understanding of the evolution and internal structure of giant planets. Here, using first-principles calculations and crystal structure predictions, we identify four stable phases of a helium-silica compound with seven/eight-coordinated silicon atoms at pressure range of 600-4000 GPa, corresponding to the interior condition of the outer planets in the solar system. The density of HeSiO2 agrees with current structure models of the planets. This helium-silica compound exhibits a superionic-like helium diffusive state at the high pressure and high temperature conditions along the isentropes of Saturn, a metallic fluid state in Jupiter, and a solid state in the deep interiors of Uranus and Neptune. The reaction of helium and silica may lead to the erosion of the rocky core of giant planets and form a diluted core region. These results highlight the reactivity of helium under high pressure to form new compounds, and also provides evidence to help build more sophisticated interior models of giant planets., 20 pages, 5 figures

All elements that form diatomic molecules, such as H_{2}, N_{2}, O_{2}, Cl_{2}, Br_{2}, and I_{2}, are destined to become atomic solids under sufficiently high pressure. However, as revealed by many experimental and theoretical studies, these elements show very different propensity and transition paths due to the balance of reduced volume, lone pair electrons, and interatomic bonds. The study of F under pressure can illuminate this intricate behavior since F, owing to its unique position on the periodic table, can be compared with H, with N and O, and also with other halogens. Nevertheless, F remains the only element whose solid structure evolution under pressure has not been thoroughly studied. Using a large-scale crystal structure search method based on first principles calculations, we find that, before reaching an atomic phase, F solid transforms first into a structure consisting of F_{2} molecules and F polymer chains and then into a structure consisting of F polymer chains and F atoms, a distinctive evolution with pressure that has not been seen in any other elements. Both intermediate structures are found to be metallic and become superconducting, a result that adds F to the elemental superconductors.

N/A

Half-Heusler thermoelectric materials are potential candidates for high thermoelectric efficiency. We report high-pressure thermoelectric and structural property measurements, density functional theory calculations on the half-Heusler material TiNiSn, and an increase of 15% in the relative dimensionless figure of merit, ZT, around 3 GPa. Thermal and electrical properties were measured utilizing a specialized sample cell assembly designed for the Paris-Edinburgh large-volume press to a maximum pressure of 3.5 GPa. High-pressure structural measurements performed up to 50 GPa in a diamond-anvil cell indicated the emergence of a new high-pressure phase around 20 GPa. A first-principles structure search performed using an ab initio random structure search approach identified the high-pressure phase as an orthorhombic type, in good agreement with the experimental results.

Incorporating bismuth, the heaviest element stable to radioactive decay, into new materials enables the creation of emergent properties such as permanent magnetism, superconductivity, and nontrivial topology. Understanding the factors that drive Bi reactivity is critical for the realization of these properties. Using pressure as a tunable synthetic vector, we can access unexplored regions of phase space to foster reactivity between elements that do not react under ambient conditions. Furthermore, combining computational and experimental methods for materials discovery at high-pressures provides broader insight into the thermodynamic landscape than can be achieved through experiment alone, informing our understanding of the dominant chemical factors governing structure formation. Herein, we report our combined computational and experimental exploration of the Mo-Bi system, for which no binary intermetallic structures were previously known. Using the ab initio random structure searching (AIRSS) approach, we identified multiple synthetic targets between 0-50 GPa. High-pressure in situ powder X-ray diffraction experiments performed in diamond anvil cells confirmed that Mo-Bi mixtures exhibit rich chemistry upon the application of pressure, including experimental realization of the computationally predicted CuAl2-type MoBi2 structure at 35.8(5) GPa. Electronic structure and phonon dispersion calculations on MoBi2 revealed a correlation between valence electron count and bonding in high-pressure transition metal-Bi structures as well as identified two dynamically stable ambient pressure polymorphs. Our study demonstrates the power of the combined computational-experimental approach in capturing high-pressure reactivity for efficient materials discovery.

Titanium dioxide has been extensively studied in the rutile or anatase phases, 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 TiO$_2$ 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 stablities of the various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO$_2$, 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 TiO$_2$ phases that are in good agreement with DFT. We construct a pressure-temperature phase diagram of TiO$_2$ in the range $0

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., Comment: 9 pages, 4 figures in main manuscript. 19 pages, 11 figures in supplementary

Superionic states are phases of matters that can simultaneously exhibit some of the properties of a fluid and of a solid. Superionic states of ice, H$_{3}$O, He-H$_{2}$O or He-NH$_{3}$ compounds have been reported in previous works. Silicon, sodium, and hydrogen are abundant elements inside the earth. Here, we use ab initio calculations to show that, at extreme conditions inside the earth, Na, Si, and H can form many hypervalent compounds that some of them can exist every close to ambient pressure, and surprisingly a previously unknown type of superionic state of $P\overline3m1 - $Na$_{2}$SiH$_{6}$ can form as well. Our work focused on new superionic state of Na$_{2}$SiH$_{6}$, and the results also reveal several different hypervalent Si-H anions discovered, which are different from individual SiH$_{5}^{\,-}$ and octahedral SiH$_{6}^{\,2-}$ in previous research of ternary alkali hypervalent silicon hydrides. Our work provides some advice on further investigations on potential ternary hydrides inside the earth.

Two dimensionless fundamental physical constants, the fine structure constant α and the proton- to-electron mass ratio mp/me are attributed a particular importance from the point of view of nuclear synthesis, formation of heavy elements, planets, and life-supporting structures. Here, we show that a combination of these two constants results in a new dimensionless constant which provides the upper bound for the speed of sound in condensed phases, vu. We find that vu/c = α (me/2mp)^1/2 , where c is the speed of light in vacuum. We support this result by a large set of experimental data and first principles computations for atomic hydrogen. Our result expands current understanding of how fundamental constants can impose new bounds on important physical properties.

© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA 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.

Highly compressed hydrides have been at the forefront of the search for high-Tc superconductivity. The recent discovery of record-high Tc's in H3S and LaH10±x under high pressure fuels the enthusiasm for finding good superconductors in similar hydride groups. Guided by first-principles structure prediction, we successfully synthesized ZrH3 and Zr4H15 at modest pressures (30-50 GPa) in diamond anvil cells by two different reaction routes: ZrH2 + H2 at room temperature and Zr + H2 at ∼1500 K by laser heating. From the synchrotron X-ray diffraction patterns, ZrH3 is found to have a Pm3̅n structure corresponding to the familiar A15 structure, and Zr4H15 has an I4̅3d structure isostructural to Th4H15. Electrical resistance measurement and the dependence of Tc on the applied magnetic field of the sample showed the emergence of two superconducting transitions at 6.4 and 4.0 K at 40 GPa, which correspond to the two Tc's for ZrH3 and Zr4H15.

Physical review materials 4(5), 055001 (2020). doi:10.1103/PhysRevMaterials.4.055001, Carbonates are the major hosts of carbon on Earth's surface and their fate during subduction needs to be known to understand the deep carbon cycle. Magnesite (MgCO$_{3}$) is thought to be an important phase participating in deep Earth processes, but its phase stability is still a matter of debate for the conditions prevalent in the lowest part of the mantle and at the core mantle boundary. Here, we have studied the phase relations and stabilities of MgCO$_{3}$ at these P,T conditions, using Raman spectroscopy at high pressures (∼148GPa) and after heating to high temperatures (∼3600K) in laser-heated diamond anvil cell experiments. The experimental Raman experiments were supplemented by x-ray powder diffraction data, obtained at a pressure of 110 GPa. Density-functional-theory-based model calculations were used to compute Raman spectra for several MgCO$_{3}$ high-pressure polymorphs, thus allowing an unambiguous assignment of Raman modes. By combining the experimental observations with the density-functional-theory results, we constrain the phase stability field of MgCO$_{3}$ with respect to the high-pressure polymorph, MgCO$_{3}$-II. We further confirm that Fe-free MgCO$_{3}$-II is a tetracarbonate with monoclinic symmetry (space group C2/m), which is stable over the entire P,T range of the Earth's lowermost mantle geotherm., Published by APS, College Park, MD

The make-up of the outer planets, and many of their moons, are dominated by matter from the H-C-N-O chemical space, commonly assumed to originate from mixtures of hydrogen and the planetary ices H$_2$O, CH$_4$, and NH$_3$. In their interiors, these ices 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. In turn, this determines planets' interior structure, thermal history, magnetic field generation, etc. Despite its importance, the H-C-N-O space has not been surveyed systematically. Asked simply: at high-pressure conditions, what compounds emerge within this space, and what governs their stability? Here, we report on results from an unbiased crystal structure search amongst H-C-N-O compounds at 5 Mbar to answer this question.

Lanthanum hydrides, containing hydrogen framework structures under compression, display a superconducting state with a high observed critical temperature. However, this phenomenon has so far only been observed at very high pressures. Here, we computationally search for superconductors with very high critical temperatures, but at much lower pressures. We uncover two such sodalite-type clathrate hydrides, YbH6 and LuH6 (Tc =145 K at P=70 GPa for YbH6 and, especially, Tc =273 K at P=100 GPa for LuH6). These striking properties are a consequence of the strong interrelationship between the f states present at the Fermi surface, structural stability and Tc value. For example, TmH6, with unfilled 4f orbitals, is stable at 50 GPa, while it has a relatively low value of Tc of 25 K. As for the YbH6 and LuH6 compounds, they have filled f-shells and the decrease of the f-energy below the Fermi level leads to formation of the nesting regions on the Fermi surface and, as a result, to phonon "softening" and an increase in Tc., Comment: 14 pages ,4 figures

The quasiperiodic structure of high-pressure bismuth gives rise to unusually strong electron-phonon coupling., We examine the low-temperature states supported by the quasiperiodic host-guest structure of elemental bismuth at high pressure, Bi-III. Our electronic transport and magnetization experiments establish Bi-III as a rare example of type II superconductivity in an element, with a record upper critical field of ~ 2.5 T, unusually strong electron-phonon coupling, and an anomalously large, linear temperature dependence of the electrical resistivity in the normal state. These properties may be attributed to the peculiar phonon spectrum of incommensurate host-guest structures, which exhibit additional quasi-acoustic sliding modes, suggesting a pathway toward strong coupling superconductivity with the potential for enhanced transition temperatures and high critical fields.

The AIRSS method generates crystal structures for m-aminobenzoic acid; comparison is made to experimental powder X-ray diffraction and MAS NMR., 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.

Little was known about the properties of hydrogen under extreme pressure. Experiments now reveal key details about the arrangement of molecules in several of the element’s high-pressure phases.

Two novel high-pressure polymorphs of lanthanum oxyhydride have been successfully predicted and stabilized under pressure. When reacted at 3 GPa, the fluorite structure of LaHO with anion-centered tetrahedral (HLa4/OLa4) geometry is transformed to the PbCl2-type structure involving coordination number increase of H- to five (HLa5 square pyramids). Upon further application of pressure to 5 GPa, LaHO changed into the anti Fe2P-type structure. Interestingly, the 5 GPa phase contains tetrahedral HLa4 and square-pyramidal OLa5 geometry, meaning coordination switching, as confirmed by ab initio calculations. The structural analysis shows that this unprecedented phenomenon is enabled by higher compressibility of hydride anion and emphasizes its potential in the search for new high-pressure forms of hydride-based mixed anion materials.

An unambiguous crystallographic structure solution for the observed phases II-VI of high pressure hydrogen does not exist due to the failure of standard structural probes at extreme pressure. In this work we propose that nuclear magnetic resonance spectroscopy provides a complementary structural probe for high pressure hydrogen. We show that the best structural models available for phases II, III, and IV of high pressure hydrogen exhibit markedly distinct nuclear magnetic resonance spectra which could therefore be used to discriminate amongst them. As an example, we demonstrate how nuclear magnetic resonance spectroscopy could be used to establish whether phase III exhibits polymorphism. Our calculations also reveal a strong renormalisation of the nuclear magnetic resonance response in hydrogen arising from quantum fluctuations, as well as a strong isotope effect. As the experimental techniques develop, nuclear magnetic resonance spectroscopy can be expected to become a useful complementary structural probe in high pressure experiments., Comment: Accepted manuscript; 7 pages, 3 figures; contains Supplemental Material

An unambiguous crystallographic structure solution for the observed phases II-VI of high pressure hydrogen does not exist due to the failure of standard structural probes at extreme pressure. In this work we propose that nuclear magnetic resonance spectroscopy provides a complementary structural probe for high pressure hydrogen. We show that the best structural models available for phases II, III, and IV of high pressure hydrogen exhibit markedly distinct nuclear magnetic resonance spectra which could therefore be used to discriminate amongst them. As an example, we demonstrate how nuclear magnetic resonance spectroscopy could be used to establish whether phase III exhibits polymorphism. Our calculations also reveal a strong renormalisation of the nuclear magnetic resonance response in hydrogen arising from quantum fluctuations, as well as a strong isotope effect. As the experimental techniques develop, nuclear magnetic resonance spectroscopy can be expected to become a useful complementary structural probe in high pressure experiments., Accepted manuscript; 7 pages, 3 figures; contains Supplemental Material

Ab initio random structure searching is employed to generate candidate structures of hydrous wadsleyite, predicting NMR parameters for experimental comparison., 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.

The energy landscape of carbon is exceedingly complex, hosting diverse and important metastable phases, including diamond, fullerenes, nanotubes, and graphene. Searching for structures, especially those with large unit cells, in this landscape is challenging. Here we use a combined stochastic search strategy employing two algorithms (ab initio random structure search and random sampling strategy combined with space group and graph theory) to apply connectivity constraints to unit cells containing up to 100 carbon atoms. We uncover three low energy carbon polymorphs (Pbam-32, P6/mmm, and I4[over ¯]3d) with new topologies, containing 32, 36, and 94 atoms in their primitive cells, respectively. Their energies relative to diamond are 96, 131, and 112 meV/atom, respectively, which suggests potential metastability. These three carbon allotropes are mechanically and dynamically stable, insulating carbon crystals with superhard mechanical properties. The I4[over ¯]3d structure possesses a direct band gap of 7.25 eV, which is the widest gap in the carbon allotrope family. Silicon, germanium, and tin versions of Pbam-32, P6/mmm, and I4[over ¯]3d also show energetic, dynamical, and mechanical stability. The computed electronic properties show that they are potential materials for semiconductor and photovoltaic applications.

Natural specimens of the pyrochlore (A2B2O7) compounds have been found to retain foreign actinide impurities within their parent framework, undergoing metamictization to a fully amorphous state. The response to radionuclide decay identifies pyrochlore systems with having high radiation tolerance and tailored use in radioactive waste applications and radionuclide sequestration. High pressure is a powerful pathway to high density states and amorphization with parallels to radiation-induced processes. Here, La2Sn2O7 is evaluated under extreme conditions via the combination of laser heating in a diamond anvil cell with X-ray diffraction and Raman spectroscopy. The measurements are supported by ab initio random structure searching and molecular dynamics calculations. A new ground state at 70 GPa is revealed, and high temperature annealing is fundamental to access its crystalline ground state and fully determine the structure. This crystalline phase ( P21/ c) retains its structural integrity during decompression and is fully recoverable to ambient conditions. The final state of the system is shown to be highly pathway dependent due to the covalent nature of the Sn-O bonding. The Tc pyrochlore, La2Tc2O7, is analyzed for similarities in the bonding to determine the likelihood of an analogous pathway dependency to a final state.

Ice is one of the most extensively studied condensed matter systems. Yet, both experimentally and theoretically several new phases have been discovered over the last years. Here we report a large-scale density-functional-theory study of the configuration space of water ice. We geometry optimise 74,963 ice structures, which are selected and constructed from over five million tetrahedral networks listed in the databases of Treacy, Deem, and the International Zeolite Association. All prior knowledge of ice is set aside and we introduce "generalised convex hulls" to identify configurations stabilised by appropriate thermodynamic constraints. We thereby rediscover all known phases (I-XVII, i, 0 and the quartz phase) except the metastable ice IV. Crucially, we also find promising candidates for ices XVIII through LI. Using the "sketch-map" dimensionality-reduction algorithm we construct an a priori, navigable map of configuration space, which reproduces similarity relations between structures and highlights the novel candidates. By relating the known phases to the tractably small, yet structurally diverse set of synthesisable candidate structures, we provide an excellent starting point for identifying formation pathways.

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.

Superionic hydrogen was previously thought to be an exotic state predicted and confirmed only in pure H2O ice. In Earth's deep interior, H2O exists in the form of O-H groups in ultra-dense hydrous minerals, which have been proved to be stable even at the conditions of the core-mantle boundary (CMB). However, the superionic states of these hydrous minerals at high P-T have not been investigated. Using first-principles calculations, we found that pyrite structured FeO2Hx (0

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.

Crystal structure prediction methods and first-principles calculations have been used to explore low-energy structures of carbon monoxide (CO). Contrary to the standard wisdom, the most stable structure of CO at ambient pressure was found to be a polymeric structure of Pna21 symmetry rather than a molecular solid. This phase is formed from six-membered (4 Carbon + 2 Oxygen) rings connected by C=C double bonds with two double-bonded oxygen atoms attached to each ring. Interestingly, the polymeric Pna21 phase of CO has a much higher energy density than trinitrotoluene (TNT). On compression to about 7 GPa, Pna21 is found to transform into another chain-like phase of Cc symmetry which has similar ring units to Pna21. On compression to 100 GPa it is energetically favorable for CO to polymerize to form a single-bonded Cmcm phase from another structure of Cmca symmetry composed of units similar to those found in the single-bonded I212121 structure. Thermodynamic stability of these structures was verified using calculations with different density functionals, including hybrid and van der Waals corrected functionals., jiansun@nju.edu.cn

We propose here a two-dimensional material based on a single layer of violet or Hittorf’s phosphorus. Using first-principles density functional theory, we find it to be energetically very stable, comparable to other previously proposed single-layered phosphorus structures. It requires only a small energetic cost of approximately 0.04 eV/atom to be created from its bulk structure, Hittorf’s phosphorus, or a binding energy of 0.3 − 0.4 J/m^2 per layer, suggesting the possibility of exfoliation in experiments. We find single-layered Hittorf’s phosphorus to be a wide band gap semiconductor with a direct band gap of approximately 2.5 eV and our calculations show it is expected to have a high and highly anisotropic hole mobility with an upper bound lying between 3000−7000 cm^2V^−1 s^ −1 . These combined properties make single-layered Hittorf’s phosphorus a very good candidate for future applications in a wide variety of technologies, in particular for high frequency electronics, and optoelectronic devices operating in the low wavelength blue color range.

The sixth blind test of organic crystal-structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal, and a bulky flexible molecule. This blind test has seen substantial growth in the number of submissions, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and "best practices" for performing CSP calculations. All of the targets, apart from a single potentially disordered Z` = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.

We propose here a two-dimensional material based on a single layer of violet or Hittorf's phosphorus. Using first-principles density functional theory, we find it to be energetically very stable, comparable to other previously proposed single-layered phosphorus structures. It requires only a small energetic cost of approximately $0.04~\text{eV/atom}$ to be created from its bulk structure, Hittorf's phosphorus, or a binding energy of $0.3-0.4~\text{J/m}^2$ per layer, suggesting the possibility of exfoliation in experiments. We find single-layered Hittorf's phosphorus to be a wide band gap semiconductor with a direct band gap of approximately $2.5$~eV and our calculations show it is expected to have a high and highly anisotropic hole mobility with an upper bound lying between $3000-7000$~cm$^2$V$^{-1}$s$^{-1}$. These combined properties make single-layered Hittorf's phosphorus a very good candidate for future applications in a wide variety of technologies, in particular for high frequency electronics, and optoelectronic devices operating in the low wavelength blue color range., 17 pages, 3 figures

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2 GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Phosphorus has received recent attention in the context of high-capacity and high-rate anodes for lithium and sodium-ion batteries. Here, we present a first principles structure prediction study combined with NMR calculations which gives us insights into its lithiation/sodiation process. We report a variety of new phases phases found by AIRSS and the atomic species swapping methods. Of particular interest, a stable Na$_5$P$_4$-C2/m structure and locally stable structures found less than 10 meV/f.u. from the convex hull, such as Li$_4$P$_3$-P2$_1$2$_1$2$_1$, NaP$_5$-Pnma and Na$_4$P$_3$-Cmcm. The mechanical stability of Na$_5$P$_4$-C2/m and Li$_4$P$_3$-P2$_1$2$_1$2$_1$ has been studied by first principles phonon calculations . We have calculated average voltages which suggests that black phosphorus (BP) can be considered as a safe anode in lithium-ion batteries due to its high lithium insertion voltage, 1.5 V; moreover, BP exhibits a relatively low theoretical volume expansion compared with other intercalation anodes, 216\% ($\Delta V/V$). We identify that specific ranges in the calculated shielding can be associated with specific ionic arrangements, results which play an important role in the interpretation of NMR spectroscopy experiments. Since the lithium-phosphides are found to be insulating even at high lithium concentrations we show that Li-P-doped phases with aluminium have electronic states at the Fermi level suggesting that using aluminium as a dopant can improve the electrochemical performance of P anodes.