It is essential to know the arrangement of the atoms in a material in order to compute and understand its properties. Searching for stable structures of materials using first-principles electronic structure methods, such as density-functional-theory (DFT), is a rapidly growing field. Here we describe our simple, elegant and powerful approach to searching for structures with DFT, which we call ab initio random structure searching (AIRSS). Applications to discovering the structures of solids, point defects, surfaces, and clusters are reviewed. New results for iron clusters on graphene, silicon clusters, polymeric nitrogen, hydrogen-rich lithium hydrides, and boron are presented.

https://iopscience.iop.org/article/10.1088/0953-8984/23/5/053201/pdf

Ab initio random structure searching

A systematic structure search in the La–H and Y–H systems under pressure reveals some hydrogen-rich structures with intriguing electronic properties. For example, LaH10 is found to adopt a sodalite-like face-centered cubic (fcc) structure, stable above 200 GPa, and LaH8 a C2/m space group structure. Phonon calculations indicate both are dynamically stable; electron phonon calculations coupled to Bardeen–Cooper–Schrieffer (BCS) arguments indicate they might be high-Tc superconductors. In particular, the superconducting transition temperature Tc calculated for LaH10 is 274–286 K at 210 GPa. Similar calculations for the Y–H system predict stability of the sodalite-like fcc YH10 and a Tc above room temperature, reaching 305–326 K at 250 GPa. The study suggests that dense hydrides consisting of these and related hydrogen polyhedral networks may represent new classes of potential very high-temperature superconductors.

Potential high-Tc superconducting lanthanum and yttrium hydrides at high pressure.

Pressure is a fundamental thermodynamic variable that can be used to control the properties of materials, because it reduces interatomic distances and profoundly modifies electronic orbitals and bonding patterns. It is thus a versatile tool for the creation of exotic materials not accessible at ambient conditions. Recently developed static and dynamic high-pressure experimental techniques have led to the synthesis of many functional materials with excellent performance: for example, superconductors, superhard materials and high-energy-density materials. Some of these advances have been aided and accelerated by first-principles crystal-structure searching simulations. In this Review, we discuss recent progress in high-pressure materials discovery, placing particular emphasis on the record high-temperature superconductivity in hydrogen sulfide and on nanotwinned cubic boron nitride and diamond, the hardest known materials. Energy materials and exotic chemical materials obtained under high pressures are also discussed. The main drawback of high-pressure materials is their destabilization after pressure release; this problem and its possible solutions are surveyed in the conclusions, which also provide an outlook on the future developments in the field. High pressure offers a unique degree of freedom for the creation of new materials, leading to new superconductors, superhard materials, high-energy-density materials and exotic chemical materials with unprecedented properties. This Review discusses these materials, along with recently developed theoretical and experimental methods for materials discovery at high pressures.

Materials discovery at high pressures

We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and nonmetals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation PBEsol GGA, the second-order GGA SOGGA, and the Armiento-Mattsson 2005 AM05 GGA. For completeness, we also test more standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties lattice constant and bulk modulus. For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zeropoint phonon and finite-temperature effects ignored by many workers. We show how GAUSSIAN basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline-earth metal, and alkali halide crystals where the maximum value of the reduced density gradient is about 2, but perform very similarly for most of the other solids where it is often about 1. Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.

https://dspace.mit.edu/openaccess-disseminate/1721.1/51869

Assessing the performance of recent density functionals for bulk solids

https://pure.mpg.de/pubman/item/item_3314332_3/component/file_3324953/1-s2.0-S0370157320300363-main.pdf

A perspective on conventional high-temperature superconductors at high pressure: Methods and materials

The cesium polyhydrides (CsHn, n > 1) are predicted to become stable, with respect to decomposition into CsH and H2, at pressures as low as 2 GPa. The CsH3 stoichiometry is found to have the lowest enthalpy of formation from CsH and H2 between 30 and 200 GPa. Evolutionary algorithms predict five distinct, mechanically stable, nearly isoenthalpic CsH3 phases consisting of H3– molecules and Cs+ atoms. The H3– sublattices in two of these adopt a hexagonal three-connected net; in the other three the net is twisted, like the silicon sublattice in the α-ThSi2 structure. The former emerge as being metallic below 100 GPa in our screened hybrid density functional theory calculations, whereas the latter remain insulating up to pressures greater than 250 GPa. The Cs+ cations in the most-stable I41/amd CsH3 phase adopt the positions of the Cs atoms in Cs-IV, and the H3– molecules are found in the (interstitial) regions which display a maximum in the electron density.

Compressed Cesium Polyhydrides: Cs+ Sublattices and H3– Three-Connected Nets

XTALOPT version r9: An open-source evolutionary algorithm for crystal structure prediction

Evolutionary structure searches predict three new phases of iodine polyhydrides stable under pressure. Insulating P1-H5I, consisting of zigzag chains of (HI)δ+ and H2δ− molecules, is stable between 30 and 90 GPa. Cmcm-H2I and P6/mmm-H4I are found on the 100, 150, and 200 GPa convex hulls. These two phases are good metals, even at 1 atm, because they consist of monatomic lattices of iodine. At 100 GPa the superconducting transition temperature, Tc, of H2I and H4I is estimated to be 7.8 and 17.5 K, respectively. The increase in Tc relative to elemental iodine results from a larger ωlog from the light mass of hydrogen and an enhanced λ from modes containing H/I and H/H vibrations.

https://pubs.acs.org/doi/pdf/10.1021/acs.jpclett.5b01839

Superconducting High-Pressure Phases Composed of Hydrogen and Iodine.

At 80 GPa, phases with the PH2 stoichiometry, which are composed of simple cubic like phosphorus layers capped with hydrogen atoms and layers of H2 molecules, are predicted to be important species contributing to the recently observed superconductivity in compressed phosphine. The electron–phonon coupling in these phases results from the motions of the phosphorus atoms and the hydrogen atoms bound to them. The role of the mobile H2 layers is to decrease the Coulomb repulsion between the negatively charged hydrogen atoms capping the phosphorus layers. An insulating PH5 phase, the structure and bonding of which is reminiscent of diborane, is also predicted to be metastable at this pressure.

/pdf/superconducting-phases-of-phosphorus-hydride-under-pressure-4m0rxk7wh0.pdf

Superconducting Phases of Phosphorus Hydride Under Pressure: Stabilization by Mobile Molecular Hydrogen

Abstract A priori crystal structure prediction techniques have been used to explore the phase diagrams of hydrides of main group elements under pressure. A number of novel phases with the chemical formulas MHn, n > 1 and M = Li, Na, K, Rb, Cs; MHn, n > 2 and M= Mg, Ca, Sr, Ba; HnI with n > 1 and PH, PH2, PH3 have been predicted to be stable at pressures achievable in diamond anvil cells. The hydrogenic lattices within these phases display a number of structural motifs including H2δ− , H−, H−3 , as well as one-dimensional and three-dimensional extended structures. A wide range of superconducting critical temperatures, Tcs, are predicted for these hydrides. The mechanism of metallization and the propensity for superconductivity are dependent upon the structural motifs present in these phases, and in particular on their hydrogenic sublattices. Phases that are thermodynamically unstable, but dynamically stable, are accessible experimentally. The observed trends provide insight on how to design hydrides that are superconducting at high temperatures.

/pdf/superconductivity-in-hydrides-doped-with-main-group-elements-4ntijh31hf.pdf

Andrew Shamp

Papers

Compressed Cesium Polyhydrides: Cs+ Sublattices and H3– Three-Connected Nets

XTALOPT version r9: An open-source evolutionary algorithm for crystal structure prediction

Superconducting High-Pressure Phases Composed of Hydrogen and Iodine.

Superconducting Phases of Phosphorus Hydride Under Pressure: Stabilization by Mobile Molecular Hydrogen

Superconductivity in Hydrides Doped with Main Group Elements Under Pressure