BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Chris J. Pickard

Chris J. Pickard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.supr-conphysics.comp-phcond-mat.otherastro-ph.EPphysics.chem-phastro-ph.SRcond-mat.mes-hallphysics.app-phphysics.geo-phquant-ph

Catalog footprint

What is connected

45works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Ephemeral data derived potentials for random structure search

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 $γ$-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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unified chemical theory of structure and bonding in elemental metals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Design Principles for High Temperature Superconductors with Hydrogen-based Alloy Backbone at Moderate Pressure

Hydrogen-based superconductors provide a route to the long-sought goal of room-temperature superconductivity, but the high pressures required to metallize these materials limit their immediate application. For example, carbonaceous sulfur hydride, the first room-temperature superconductor, can reach a critical temperature (Tc) of 288 K only at the extreme pressure of 267 GPa. The next recognized challenge is the realization of room-temperature superconductivity at significantly lower pressures. Here, we propose a strategy for the rational design of high-temperature superconductors at low pressures by alloying small-radius elements and hydrogen to form ternary hydride superconductors with alloy backbones. We identify a hitherto unknown fluorite-type backbone in compositions of the form AXH8, which exhibit high temperature superconductivity at moderate pressures. The Fm-3m phase of LaBeH8, with a fluorite-type H-Be alloy backbone, is predicted to be metastable and superconducting with a Tc ~ 191 K at 50 GPa; a substantially lower pressure than that required by the geometrically similar clathrate hydride LaH10 (170 GPa). Our approach paves the way for finding high-Tc ternary hydride superconductors at conditions close to ambient pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-pressure phase behaviors of titanium dioxide revealed by a $Δ$-learning potential

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<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 observe predominantly anatase-to-baddeleyite transformation at about 20~GPa, via a martensitic two-step mechanism with 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 TiO$_2$ phases, and our method is generally applicable to other functional materials with multiple polymorphs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Partially Diffusive Helium-Silica Compound in the Deep Interiors of Giant Planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient prediction of Nucleus Independent Chemical Shifts for polycyclic aromatic hydrocarbons

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 magnitudes faster than DFT calculations, is very promising and can be further refined in order to study more complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting the phase diagram of titanium dioxide with random search and pattern recognition

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stability and superconductivity of lanthanum and yttrium decahydrides

Rare-earth hydrides can exhibit high-temperature superconductivity under high pressure. Here, we apply a crystal structure prediction method to the current record-holding $T_c$ material, LaH$_{10}$, and a candidate for even higher $T_c$, YH$_{10}$. We find a pressure-induced phase transition from the experimentally observed cubic LaH$_{10}$ phase to a new hexagonal phase at around 420 GPa. This hexagonal phase could explain experimental observations of hcp impurities in fcc samples. We find that YH$_{10}$ forms similar structures to LaH$_{10}$ and discuss the sensitivity of superconductivity calculations to the computational parameters used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure-specific, mode-resolved phonon coherence and specularity at graphene grain boundaries

In spite of their importance for understanding phonon transport phenomena in thin films and polycrystalline solids, the effects of boundary roughness scattering on phonon specularity and coherence are poorly understood because there is no general method for predicting their dependence on phonon momentum, frequency, branch and boundary morphology. Using the recently formulated atomistic S-matrix method, we develop a theory of boundary roughness scattering to determine the mode-resolved phonon coherence and specularity parameters from the scattering amplitudes. To illustrate the theory, we apply it to phonon scattering in realistic nonsymmetric graphene grain boundary (GB) models derived from atomic structure predictions. The method is validated by comparing its predictions with frequency-resolved results from lattice dynamics-based calculations. We prove that incoherent scattering is almost perfectly diffusive. We show that phonon scattering at the graphene GB is not diffuse although coherence and specularity are significantly reduced for long-wavelength flexural acoustic phonons. Our approach can be generalized to other atomistic boundary models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Carbon Nitride Frameworks and Dense Crystalline Polymorphs

We used ab initio random structure searching (AIRSS) to investigate polymorphism in C3N4 carbon nitride as a function of pressure. Our calculations reveal new framework structures, including a particularly stable chiral polymorph of space group P43212 containing mixed sp2 and sp3-bonding, that we have produced experimentally and recovered to ambient conditions. As pressure is increased a sequence of structures with fully sp3-bonded C atoms and three-fold coordinated N atoms is predicted, culminating in a dense Pnma phase above 250 GPa. Beyond 650 GPa we find that C3N4 becomes unstable to decomposition into diamond and pyrite-structured CN2.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hexagonal structure of phase III of solid hydrogen

A hexagonal structure of solid molecular hydrogen with $P6_122$ symmetry is calculated to be more stable below about 200 GPa than the monoclinic $C2/c$ structure identified previously as the best candidate for phase III. We find that the effects of nuclear quantum and thermal vibrations play a central role in the stabilization of $P6_122$. The $P6_122$ and $C2/c$ structures are very similar and their Raman and infra-red data are in good agreement with experiment. However, our calculations show that the hexagonal $P6_122$ structure provides better agreement with the available x-ray diffraction data than the $C2/c$ structure at pressures below about 200 GPa. We suggest that two phase-III-like structures may be formed at high pressures, hexagonal $P6_122$ below about 200 GPa and monoclinic $C2/c$ at higher pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-pressure Phase Stability and Superconductivity of Pnictogen Hydrides and Chemical Trends for Compressed Hydrides

Binary hydrides formed by the pnictogens of phosphorus, arsenic and antimony are studied at high pressures using first principles methods. Stable structures are predicted and their electronic, vibrational and superconducting properties are investigated. We predict that SbH$_{4}$ and AsH$_{8}$ will be high-temperature superconductors at megabar pressures, with critical temperatures in excess of 100 K. The highly symmetric hexagonal SbH$_{4}$ phase is predicted to be stabilized above about 150 GPa, which is readily achievable in diamond anvil cell experiments. We find that all phosphorus hydrides are metastable with respect to decomposition into the elements within the pressure range studied. Trends based on our results and literature data reveal a connection between the high-pressure behaviors and ambient-pressure chemical quantities which provides insight into understanding which elements may form hydrogen-rich high-temperature superconducting phases at high pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$\textit{Ab initio}$ study of phosphorus anodes for lithium and sodium-ion batteries

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\% ($Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D ice from first principles: structures and phase transitions

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 \textit{ca.} 2 GPa at which point square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. We also find a double layer AA stacked square ice phase, which clarifies the difference between experimental observations and earlier force field simulations. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and width.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational searches for iron oxides at high pressures

We have used density-functional-theory methods and the ab initio random structure searching (AIRSS) approach to predict stable structures and stoichiometries of mixtures of iron and oxygen at high pressures. Searching was performed for 12 different stoichiometries at pressures of 100, 350 and 500 GPa, which involved relaxing more than 32,000 structures. We find that Fe$_2$O$_3$ and FeO$_2$ are the only phases stable to decomposition at 100 GPa, while at 350 and 500 GPa several stoichiometries are found to be stable or very nearly stable. We report a new structure of Fe$_2$O$_3$ with $P2_12_12_1$ symmetry which is found to be more stable than the known Rh$_2$O$_3$(II) phase at pressures above $\sim$233 GPa. We also report two new structures of FeO, with $Pnma$ and $R\bar{3}m$ symmetries, which are found to be stable within the ranges 195-285 GPa and 285-500 GPa, respectively, and two new structures of Fe$_3$O$_4$ with $Pca2_1$ and $P2_1/c$ symmetries, which are found to be stable within the ranges 100-340 GPa and 340-500 GPa, respectively. Finally, we report two new structures of Fe$_4$O$_5$ with $P4_2/n$ and $P\bar{3}m1$ symmetries, which are found to be stable within the ranges 100-231 GPa and 231-500 GPa, respectively. Our new structures of Fe$_3$O$_4$ and Fe$_4$O$_5$ are found to have lower enthalpies than their known structures within their respective stable pressure ranges.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dissociation products and structures of solid H2S at strong compression

Hydrogen sulfides have recently received a great deal of interest due to the record high superconducting temperatures of up to 203 K observed on strong compression of dihydrogen sulfide (H2S). A joint theoretical and experimental study is presented in which decomposition products and structures of compressed H2S are characterized, and their superconducting properties are calculated. In addition to the experimentally known H2S and H3S phases, our first-principles structure searches have identified several energetically competitive stoichiometries that have not been reported previously; H2S3, H3S2, and H4S3. In particular, H4S3 is predicted to be thermodynamically stable within a large pressure range of 25-113 GPa. High-pressure room-temperature X-ray diffraction measurements confirm the presence of H3S and H4S3 through decomposition of H2S that emerge at 27 GPa and coexist with residual H2S, at least up to the highest pressure studied in our experiments of 140 GPa. Electron-phonon coupling calculations show that H4S3 has a small Tc of below 2 K, and that H2S is mainly responsible for the observed superconductivity of samples prepared at low temperature (<100K).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-energy tetrahedral polymorphs of carbon, silicon, and germanium

Searches for low-energy tetrahedral polymorphs of carbon and silicon have been performed using density functional theory computations and the ab initio random structure searching (AIRSS) ap- proach. Several of the hypothetical phases obtained in our searches have enthalpies that are lower or comparable to those of other polymorphs of group 14 elements that have either been experimentally synthesized or recently proposed as the structure of unknown phases obtained in experiments, and should thus be considered as particularly interesting candidates. A structure of P bam symmetry with 24 atoms in the unit cell was found to be a low energy, low-density metastable polymorph in carbon, silicon, and germanium. In silicon, Pbam is found to have a direct band gap at the zone center with an estimated value of 1.4 eV, which suggests applications as a photovoltaic material. We have also found a low-energy chiral framework structure of P 41 21 2 symmetry with 20 atoms per cell containing fivefold spirals of atoms, whose projected topology is that of the so-called Cairo-type two- dimensional pentagonal tiling. We suggest that P41 21 2 is a likely candidate for the structure of the unknown phase XIII of silicon. We discuss Pbam and P41 21 2 in detail, contrasting their energetics and structures with those of other group 14 elements, particularly the recently proposed P42 /ncm structure, for which we also provide a detailed interpretation as a network of tilted diamond-like tetrahedra.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Metallic Icosahedron Phase of Sodium at Terapascal Pressures

Alkali metals exhibit unexpected structures and electronic behavior at high pressures. Compression of metallic sodium (Na) to 200 GPa leads to the stability of a wide-band-gap insulator with the double hexagonal hP4 structure. Post-hP4 structures remain unexplored, but they are important for addressing the question of the pressure at which Na reverts to a metal. Here we report the reentrant metallicity of Na at the very high pressure of 15.5 terapascal (TPa), predicted using first-principles structure searching simulations. Na is therefore insulating over the large pressure range of 0.2-15.5 TPa. Unusually, Na adopts an oP8 structure at pressures of 117-125 GPa, and the same oP8 structure at 1.75-15.5 TPa. Metallization of Na occurs on formation of a stable and striking body-centered cubic cI24 electride structure consisting of Na12 icosahedra, each housing at its center about one electron which is not associated with any Na ions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Hydrogen-Bond Symmetrization and High-Temperature Superconductivity in Hydrogen Sulfide

Hydrogen compounds are peculiar as the quantum nature of the proton can crucially affect their structural and physical properties. A remarkable example are the high-pressure phases of H$_2$O, where quantum proton fluctuations favor the symmetrization of the H bond and lower by 30 GPa the boundary between the asymmetric structure and the symmetric one. Here we show that an analogous quantum symmetrization occurs in the recently discovered sulfur hydride superconductor with the record superconducting critical temperature $T_c=203$ K at 155 GPa. In this system, according to classical theory, superconductivity occurs via formation of a structure of stoichiometry H$_3$S with S atoms arranged on a body-centered-cubic (bcc) lattice. For $P \gtrsim 175$ GPa, the H atoms are predicted to sit midway between two S atoms, in a structure with $Im\bar3m$ symmetry. At lower pressures the H atoms move to an off-center position forming a short H$-$S covalent bond and a longer H$\cdots$S hydrogen bond, in a structure with $R3m$ symmetry. X-ray diffraction experiments confirmed the H$_3$S stoichiometry and the S lattice sites, but were unable to discriminate between the two phases. Our present ab initio density-functional theory (DFT) calculations show that the quantum nuclear motion lowers the symmetrization pressure by 72 GPa. Consequently, we predict that the $Im\bar3m$ phase is stable over the whole pressure range within which a high $T_c$ was measured. The observed pressure-dependence of $T_c$ is closely reproduced in our calculations for the $Im\bar3m$ phase, but not for the $R3m$ phase. Thus, the quantum nature of the proton completely rules the superconducting phase diagram of H$_3$S.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures

Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is nonmetallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Backbone NxH Compounds at High Pressures

Optical and synchrotron x-ray diffraction diamond anvil cell experiments have been combined with first principles theoretical structure predictions to investigate mixed N2 and H2 up to 55 GPa. We found the formation of oligomeric NxH (x>1) compounds using mechano- and photochemistry at pressures above 47 and 10 GPa, respectively, and room temperature. These compounds can be recovered to ambient pressure at T<130 K, whereas at room temperature, they can be metastable down to 3.5 GPa. Our results suggest new pathways for synthesis of environmentally benign high energy-density materials and alternative planetary ice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predicting interface structures: From SrTiO$_3$ to graphene

We present here a fully first-principles method for predicting the atomic structure of interfaces. Our method is based on the {\it ab initio} random structure searching (AIRSS) approach, applied here to treat two dimensional defects. The method relies on repeatedly generating random structures in the vicinity of the interface and relaxing them within the framework of density functional theory (DFT). The method is simple, requiring only a small set of parameters that can be easily connected to the chemistry of the system of interest, and efficient, ideally adapted to high-throughput first-principles calculations on modern parallel architectures. Being first-principles, our method is transferable, an important requirement for a generic computational method for the determination of the structure of interfaces. Results for two structurally and chemically very different interfaces are presented here, grain boundaries in graphene and grain boundaries in strontium titanate (SrTiO$_3$). We successfully find a previously unknown low energy grain boundary structure for the graphene system, as well as recover the previously known higher energy structures. For the SrTiO$_3$ system we study both stoichiometric and non-stoichiometric compositions near the grain boundary and find previously unknown low energy structures for all stoichiometries. We predict that these low energy structures have long-range distortions to the ground state crystal structure emanating into the bulk from the interface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structures and stability of calcium and magnesium carbonates at mantle pressures

Ab initio random structure searching (AIRSS) and density functional theory methods are used to predict structures of calcium and magnesium carbonate (CaCO$_3$ and MgCO$_3$) at high pressures. We find a previously unknown CaCO$_3$ structure which is more stable than the aragonite and "post aragonite" phases in the range 32--48 GPa. At pressures from 67 GPa to well over 100 GPa the most stable phase is a previously unknown CaCO$_3$ structure of the pyroxene type with fourfold coordinated carbon atoms. We also predict a stable structure of MgCO$_3$ in the range 85--101 GPa. Our results lead to a revision of the phase diagram of CaCO$_3$ over more than half the pressure range encountered within the Earth's mantle, and smaller changes to the phase diagram of MgCO$_3$. We predict CaCO$_3$ to be more stable than MgCO$_3$ in the Earth's mantle above 100 GPa, and that CO$_2$ is not a thermodynamically stable compound under deep mantle conditions. Our results have significant implications for understanding the Earth's deep carbon cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Synthesis of Sodium Polyhydrides at High Pressures

The only known compound of sodium and hydrogen is archetypal ionic NaH. Application of high pressure is known to promote states with higher atomic coordination, but extensive searches for polyhydrides with unusual stoichiometry remain unsuccessful in spite of several theoretical predictions. Here we report the first observation of formation of polyhydrides of Na (NaH3 and NaH7) above 40 GPa and 2000 K. We combined synchrotron x-ray diffraction and Raman spectroscopy in a laser heated diamond anvil cell and theoretical Ab-Initio Random Structure search, which both agree in stable structures and compositions. Our results support the formation of multicenter bonding in a material with unusual stoichiometry. These results are applicable to the design of new energetic solids and high-temperature superconductors based on hydrogenrich materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Temperature effects in first-principles solid state calculations of the chemical shielding tensor made simple

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 beta-aspartyl-L-alanine. We therefore propose as the method of choice to incorporate the effects of temperature in solid state chemical shielding tensor calculations 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Classical and quantum ordering of protons in cold solid hydrogen under megabar pressures

A combination of state-of-the-art theoretical methods has been used to obtain an atomic-level picture of classical and quantum ordering of protons in cold high-pressure solid hydro-gen. We focus mostly on phases II and III of hydrogen, exploring the effects of quantum nuclear motion on certain features of these phases (through a number of ab initio path-integral molecular dynamics (PIMD) simulations at particular points on the phase diagram). We also examine the importance of van der Waals forces in this system by performing calculations using the optB88-vdW density functional, which accounts for non-local correlations. Our calculations reveal that the transition between phases I and II is strongly quantum in nature, resulting from a competition between anisotropic inter-molecular interactions that restrict molecular rotation and thermal plus quantum fluctuations of the nuclear positions that facilitate it. The transition from phase II to III is more classical because quantum nuclear motion plays only a secondary role and the transition is determined primarily by the underlying potential energy surface. A structure of P21c symmetry with 24 atoms in the primitive unit cell is found to be stable when anharmonic quantum nuclear vibrational motion is included at finite temperatures using the PIMD method. This structure gives a good account of the infra-red (IR) and Raman vibron frequencies of phase II. We find ad-ditional support for C2/c structure as a strong candidate for phase III, since it remains transparent up to 300 GPa, even when quantum nuclear effects are included. Finally, we find that accounting for van der Waals forces improves the agreement between experiment and theory for the parts of the phase diagram considered, when compared to previous work which employed the widely-used Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Decomposition and terapascal phases of water ice

Computational searches for stable and metastable structures of water ice and other H:O compositions at TPa pressures have led us to predict that H$_2$O decomposes into H$_2$O$_2$ and a hydrogen-rich phase at pressures of a little over 5 TPa. The hydrogen-rich phase is stable over a wide range of hydrogen contents, and it might play a role in the erosion of the icy component of the cores of gas giants as H$_2$O comes into contact with hydrogen. Metallization of H$_2$O is predicted at a higher pressure of just over 6 TPa, and therefore H$_2$O does not have a thermodynamically stable low-temperature metallic form. We have also found a new and rich mineralogy of complicated water ice phases which are more stable in the pressure range 0.8--2 TPa than any predicted previously.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron-Phonon Coupling and the Metalization of Solid Helium at Terapascal Pressures

Solid He is studied in the pressure and temperature ranges 1-40 TPa and 0-10,000 K using first-principles methods. Anharmonic vibrational properties are calculated within a self-consistent field framework, including the internal and free energies, density-pressure relation, stress tensor, thermal expansion, and the electron-phonon coupling renormalization of the electronic band gap. We find that an accurate description of electron-phonon coupling requires us to use a non-perturbative approach. The metalization pressure of 32.9 TPa at 0 K is larger than found previously. The vibrational effects are large; for example at P=30 TPa the band gap is increased by 2.8 eV by electron-phonon coupling and a further 0.1 eV by thermal expansion compared to the static value. The implications of the calculated metalization pressure for the cooling of white dwarfs are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High Energy Density Mixed Polymeric Phase From Carbon Monoxide and Nitrogen

Carbon monoxide and nitrogen are among the potentially interesting high-energy density materials. However, in spite of the physical similarities of the molecules, they behave very differently at high pressures. Using density functional theory and structural prediction methods, we examine the ability of these systems to combine their respective properties and form novel mixed crystalline phases under pressures of up to 100 GPa. Interestingly, we find that CO catalyzes the molecular dissociation of N2, which means mixed structures are favored at a relatively low pressure (below 18 GPa), and that a three-dimensional framework with Pbam symmetry becomes the most stable phase above 52 GPa, i.e., at much milder conditions than in pure solid nitrogen. This structure is dynamically stable at ambient pressure and has an energy density of approximately 2:2 kJ/g, making it a candidate for a high-energy density material, and one that could be achieved at less prohibitive experimental conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lithiation of silicon via lithium Zintl-defect complexes

An extensive search for low-energy lithium defects in crystalline silicon using density-functional-theory methods and the ab initio random structure searching (AIRSS) method shows that the four-lithium-atom substitutional point defect is exceptionally stable. This defect consists of four lithium atoms with strong ionic bonds to the four under-coordinated atoms of a silicon vacancy defect, similar to the bonding of metal ions in Zintl phases. This complex is stable over a range of silicon environments, indicating that it may aid amorphization of crystalline silicon and form upon delithiation of the silicon anode of a Li-ion rechargeable battery.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stable all nitrogen metallic salt at terapascal pressures

The phase diagram and equation of state of dense nitrogen are of interest in understanding the fundamental physics and chemistry under extreme conditions, including planetary processes, and in discovering new materials. We predict several stable phases of nitrogen at multi-TPa pressures, including a P4/nbm structure consisting of partially charged N2 pairs and N5 tetrahedra, which is stable in the range 2.5-6.8 TPa. This is followed by a modulated layered structure between 6.8 and 12.6 TPa, which also exhibits a significant charge transfer. The P4/nbm metallic nitrogen salt and the modulated structure are stable at high pressures and temperatures, and they exhibit strongly ionic features and charge density distortions, which is unexpected in an element under such extreme conditions and could represent a new class of nitrogen materials. The P-T phase diagram of nitrogen at TPa pressures is investigated using quasiharmonic phonon calculations and ab initio molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Xenon Reacts with Iron at the Conditions of the Earth's Core

Studies of the Earth's atmosphere have shown that more than 90% of xenon (Xe) is depleted compared with its abundance in chondritic meteorites. This long-standing missing Xe paradox has become the subject of considerable interest and several models for a Xe reservoir have been proposed. Whether the missing Xe is hiding in the Earth's core has remained a long unanswered question. The key to address this issue lies in the reactivity of Xe with iron (Fe, the main constituent of the Earth's core), which has been denied by earlier studies. Here we report on the first evidence of the chemical reaction of Xe and Fe at the conditions of the Earth's core, predicted through first-principles calculations and unbiased structure searching techniques. We find that Xe and Fe form a stable, inter-metallic compound of XeFe3, adopting a Cu3Au-type face-centered cubic structure above 183 GPa and at 4470 K. As the result of a Xe -> Fe charge transfer, Xe loses its chemical inertness by opening up the filled 5p electron shell and functioning as a 5p-like element, whilst Fe is unusually negatively charged, acting as an oxidant rather than a reductant as usual. Our work establishes that the Earth's core is a natural reservoir for Xe storage, and possibly provides the key to unlocking the missing Xe paradox.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Density functional theory study of phase IV of solid hydrogen

We have studied solid hydrogen up to pressures of 300 GPa and temperatures of 350 K using density functional theory methods and have found "mixed structures" that are more stable than those predicted earlier. Mixed structures consist of alternate layers of strongly bonded molecules and weakly bonded graphene-like sheets. Quasi-harmonic vibrational calculations show that mixed structures are the most stable at room temperature over the pressure range 250-295 GPa. These structures are stabilized with respect to strongly-bonded molecular phases at room temperature by the presence of lower-frequency vibrational modes arising from the graphene-like sheets. Our results for the mixed structures are consistent with the experimental Raman data [M. I. Eremets and I. A. Troyan, Nature Mater. 10 927 (2011) and R. T. Howie et al. Phys. Rev. Lett. 108 125501 (2012)]. We find that mixed phases are reasonable structural models for phase IV of hydrogen.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diamondoid Structure of Polymeric Nitrogen at High Pressures

High-pressure polymeric structures of nitrogen have attracted great attention owing to their potential application as high-energy-density materials. We report the density functional structural prediction of the unexpected stabilization of a diamondoid (or N10-cage) structure of polymeric nitrogen at high pressures. The structure adopts a highly symmetric body-centered cubic form with lattice sites occupied by N10 tetracyclic cages, each of which consists of 10 atoms and is covalently bonded with its six next-nearest N10 cages. The prediction of this diamondoid structure rules out the earlier proposed helical tunnel phase and demonstrates the high-order nature of polymeric nitrogen at extreme high pressures. Diamondoid nitrogen is a wide-gap insulator and energetically more favorable than the experimental cubic gauche and previously predicted layered Pba2 phases above 263 GPa, a pressure which is accessible to high pressure experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energetics of hydrogen/lithium complexes in silicon analyzed using the Maxwell construction

We have studied hydrogen/lithium complexes in crystalline silicon using density-functional-theory methods and the ab initio random structure searching (AIRSS) method for predicting structures. A method based on the Maxwell construction and convex hull diagrams is introduced which gives a graphical representation of the relative stabilities of point defects in a crystal and enables visualization of the changes in stability when the chemical potentials are altered. We have used this approach to study lithium and hydrogen impurities in silicon, which models aspects of the anode material in the recently-suggested lithium-ion batteries. We show that hydrogen may play a role in these anodes, finding that hydrogen atoms bind to three-atom lithium clusters in silicon, forming stable {H,3Li} and {2H,3Li} complexes, while the {H,2Li} complex is almost stable.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum simulation of low-temperature metallic liquid hydrogen

Experiments and computer simulations have shown that the melt-ing temperature of solid hydrogen drops with pressure above about 65 GPa, suggesting that a liquid state might exist at low temperatures. It has also been suggested that this low temperature liquid state might be non-molecular and metallic, although evidence for such behaviour is lacking. Here, we report results for hydrogen at high pressures using ab initio path-integral molecular dynamics methods, which include a description of the quantum motion of the protons at finite temperatures. We have determined the melting temperature as a function of pressure by direct simulation of the coexistence of the solid and liquid phases and have found an atomic solid phase from 500 to 800 GPa which melts at <200 K. Beyond this and up to pressures of 1200 GPa a metallic atomic liquid is stable at temperatures as low as 50 K. The quantum motion of the protons is critical to the low melting temperature in this system as ab initio simulations with classical nuclei lead to a considerably higher melting temperature of \sim300 K across the entire pressure range considered.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ab initio Random Structure Searching

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High pressure ionic and molecular crystals of ammonia monohydrate within density functional theory

A combination of first-principles density functional theory calculations and a search over structures predicts the stability of a proton-transfer modification of ammonia monohydrate with space group P4/nmm. The phase diagram is calculated with the 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 (AMH-I) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Persistence and eventual demise of oxygen molecules at terapascal pressures

Computational searches for structures of solid oxygen under pressures in the multi TPa range have been carried out using density-functional-theory methods. We find that molecular oxygen persists to about 1.9 TPa at which it transforms into a semiconducting square spiral-like polymeric structure (I41/acd) with a band gap of about 3.0 eV. Solid oxygen forms a metallic zig-zag chain-like structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure is increased and a structure of Fmmm symmetry forms at about 9.5 TPa in which each atom has four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with compression, although the structure becomes more symmetric. The electronic properties of oxygen have a complex evolution with pressure, swapping between insulating, semiconducting and metallic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Predicted formation of superconducting platinum-hydride crystals under pressure in the presence of molecular hydrogen

Noble metals adopt close-packed structures at ambient pressure and rarely undergo structural transformation at high pressures. Platinum (Pt), in particular, is normally considered to be unreactive and is therefore not expected to form hydrides under pressure. We predict that platinum hydride (PtH) has a lower enthalpy than its constituents solid Pt and molecular hydrogen at pressures above 21.5 GPa. We have calculated structural phase transitions from tetragonal to hexagonal close-packed or face-centered cubic (fcc) PtH between 70 and 80 GPa. Linear response calculations indicate that PtH is a superconductor at these pressures with a critical temperature of about 10--25 K. These findings help to shed light on recent observations of pressure-induced metallization and superconductivity in hydrogen-rich materials. We show that formation of fcc metal hydrides under pressure is common among noble metal hydrides and examine the possibility of superconductivity in these materials.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hydrogen/nitrogen/oxygen defect complexes in silicon from computational searches

Point defect complexes in crystalline silicon composed of hydrogen, nitrogen, and oxygen atoms are studied within density-functional theory (DFT). Ab initio Random Structure Searching (AIRSS) is used to find low-energy defect structures. We find new lowest-energy structures for several defects: the triple-oxygen defect, {3O}, triple oxygen with a nitrogen atom, {N, 3O}, triple nitrogen with an oxygen atom, {3N,O}, double hydrogen and an oxygen atom, {2H,O}, double hydrogen and oxygen atoms, {2H,2O} and four hydrogen/nitrogen/oxygen complexes, {H,N,O}, {2H,N,O}, {H,2N,O} and {H,N,2O}. We find that some defects form analogous structures when an oxygen atom is replaced by a NH group, for example, {H,N,2O} and {3O}, and {H,N} and {O}. We compare defect formation energies obtained using different oxygen chemical potentials and investigate the relative abundances of the defects.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Post-Cotunnite phase of TeO2

We have used first-principles density-functional-theory methods with a random-structure-searching technique to determine the structure of the previously unidentified post-cotunnite phase of TeO2. Our calculations indicate a transition from the cotunnite to post-cotunnite phase at 130 GPa. The predicted post-cotunnite structure has P2_1/m space group symmetry and its calculated x-ray diffraction pattern is in good agreement with the available experimental data. We find that the cotunnite phase re-enters at about 260 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Hydrogen/silicon complexes in silicon from computational searches

Defects in crystalline silicon consisting of a silicon self-interstitial atom and one, two, three, or four hydrogen atoms are studied within density-functional theory (DFT). We search for low-energy defects by starting from an ensemble of structures in which the atomic positions in the defect region have been randomized. We then relax each structure to a minimum in the energy. We find a new defect consisting of a self-interstitial and one hydrogen atom (denoted by {I,H}) which has a higher symmetry and a lower energy than previously reported structures. We recover the {I,H_2} defect found in previous studies and confirm that it is the most stable such defect. Our best {I,H_3} defect has a slightly different structure and lower energy than the one previously reported, and our lowest energy {I,H_4} defect is different to those of previous studies.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Electronic structure of superconducting graphite intercalate compounds: The role of the interlayer state

Although not an intrinsic superconductor, it has been long--known that, when intercalated with certain dopants, graphite is capable of exhibiting superconductivity. Of the family of graphite--based materials which are known to superconduct, perhaps the most well--studied are the alkali metal--graphite intercalation compounds (GIC) and, of these, the most easily fabricated is the C${}_8$K system which exhibits a transition temperature $\bm{T_c\simeq 0.14} $K. By increasing the alkali metal concentration (through high pressure fabrication techniques), the transition temperature has been shown to increase to as much as $\bm 5 $K in C${}_2$Na. Lately, in an important recent development, Weller \emph{et al.} have shown that, at ambient conditions, the intercalated compounds \cyb and \cca exhibit superconductivity with transition temperatures $\bm{T_c\simeq 6.5} $K and $\bm{11.5} $K respectively, in excess of that presently reported for other graphite--based compounds. We explore the architecture of the states near the Fermi level and identify characteristics of the electronic band structure generic to GICs. As expected, we find that charge transfer from the intercalant atoms to the graphene sheets results in the occupation of the $\bmπ$--bands. Yet, remarkably, in all those -- and only those -- compounds that superconduct, we find that an interlayer state, which is well separated from the carbon sheets, also becomes occupied. We show that the energy of the interlayer band is controlled by a combination of its occupancy and the separation between the carbon layers.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

All-electron magnetic response with pseudopotentials: NMR chemical shifts

A theory for the ab initio calculation of all-electron NMR chemical shifts in insulators using pseudopotentials is presented. It is formulated for both finite and infinitely periodic systems and is based on an extension to the Projector Augmented Wave approach of Bloechl [P. E. Bloechl, Phys. Rev. B 50, 17953 (1994)] and the method of Mauri et al [F. Mauri, B.G. Pfrommer, and S.G. Louie, Phys. Rev. Lett. 77, 5300 (1996)]. The theory is successfully validated for molecules by comparison with a selection of quantum chemical results, and in periodic systems by comparison with plane-wave all-electron results for diamond.

0 citations0 reviewsNo code linkedOpen access