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

Richard J. Needs

Richard J. Needs 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-phphysics.chem-phcond-mat.othercond-mat.str-elquant-phastro-ph.EPcond-mat.mes-hallcond-mat.quant-gasphysics.atm-clusphysics.atom-ph

Catalog footprint

What is connected

31works
12topics
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

31 published item(s)

preprint2020arXiv

Predicting novel superconducting hydrides using machine learning approaches

Searching for superconducting hydrides has so far largely focused on finding materials exhibiting the highest possible critical temperatures ($T_c$). This has led to a bias towards materials stabilised at very high pressures, which introduces a number of technical difficulties in experiment. Here we apply machine learning methods in an effort to identify superconducting hydrides which can operate closer to ambient conditions. The output of these models informs structure searches, from which we identify and screen stable candidates before performing electron-phonon calculations to obtain $T_c$. Hydrides of alkali and alkaline earth metals are identified as particularly promising; a $T_c$ of up to 115 K is calculated for RbH$_{12}$ at 50 GPa and a $T_c$ of up to 90 K is calculated for CsH$_7$ at 100 GPa.

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

Structures and finite-temperature abundances of defects in In$_2$O$_3$-II from first-principles calculations

We have studied intrinsic defect complexes in In$_2$O$_3$ using ab initio random structure searching (AIRSS). Our first-principles density-functional-theory calculations predict the thermodynamic stability of several novel defect structures. We combine the static lattice energy and harmonic vibrational energy with the often-neglected configurational entropy to construct the free energy, which is minimised to predict defect abundances at finite temperatures. We predict that some of our new defect structures - in particular our {In,2V$_\mathrm{O}$} and {2In,3V$_\mathrm{O}$} defects - can exist in significant abundances at finite temperatures, and their densities of electronic states indicate that they could play an important role in the unexpectedly high density of n-type charge carriers observed in In$_2$O$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Vibrational renormalisation of the electronic band gap in hexagonal and cubic ice

Electron-phonon coupling in hexagonal and cubic water ice is studied using first-principles quantum mechanical methods. We consider 29 distinct hexagonal and cubic ice proton-orderings with up to 192 molecules in the simulation cell to account for proton-disorder. We find quantum zero-point vibrational corrections to the minimum electronic band gaps ranging from -1.5 to -1.7 eV, which leads to improved agreement between calculated and experimental band gaps. Anharmonic nuclear vibrations play a negligible role in determining the gaps. Deuterated ice has a smaller band-gap correction at zero-temperature of -1.2 to -1.4eV. Vibrations reduce the differences between the electronic band gaps of different proton-orderings from around 0.17 eV to less than 0.05 eV, so that the electronic band gaps of hexagonal and cubic ice are almost independent of the proton-ordering when quantum nuclear vibrations are taken into account. The comparatively small reduction in the band gap over the temperature range 0-240 K of around 0.1 eV does not depend on the proton ordering, or whether the ice is protiated or deuterated, or hexagonal or cubic. We explain this in terms of the atomistic origin of the strong electron-phonon coupling in ice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice

We use extensive first-principles quantum mechanical calculations to show that, although the static lattice and harmonic vibrational energies are almost identical, the anharmonic vibrational energy of hexagonal ice is significantly lower than that of cubic ice. This difference in anharmonicity is crucial, stabilising hexagonal ice compared with cubic ice by at least 1.4 meV/H2O, in agreement with experimental estimates. The difference in anharmonicity arises predominantly from molecular O-H bond stretching vibrational modes and is related to the different stacking of atomic layers.

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

Giant electron-phonon interactions in molecular crystals and the importance of non-quadratic coupling

We investigate electron-phonon coupling in the molecular crystals CH$_4$, NH$_3$, H$_2$O, and HF, using first-principles quantum mechanical calculations. We find vibrational corrections to the electronic band gaps at zero temperature of -1.97 eV, -1.01 eV, -1.52 eV, and -1.62 eV, respectively, which are comparable in magnitude to those from electron-electron correlation effects. Microscopically, the strong electron-phonon coupling arises in roughly equal measure from the almost dispersionless high-frequency molecular modes and from the lower frequency lattice modes. We also highlight the limitations of the widely used Allen-Heine-Cardona theory, which gives significant discrepancies compared to our more accurate treatment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrogen sulphide at high pressure: a strongly-anharmonic phonon-mediated superconductor

We use first principles calculations to study structural, vibrational and superconducting properties of H$_2$S at pressures $P\ge 200$ GPa. The inclusion of zero point energy leads to two different possible dissociations of H$_2$S, namely 3H$_2$S $\to$ 2H$_3$S + S and 5H$_2$S $\to$ 3H$_3$S + HS$_2$, where both H$_3$S and HS$_2$ are metallic. For H$_3$S, we perform non-perturbative calculations of anharmonic effects within the self-consistent harmonic approximation and show that the harmonic approximation strongly overestimates the electron-phonon interaction ($λ\approx 2.64$ at 200 GPa) and T$_c$. Anharmonicity hardens HS bond-stretching modes and softens H--S bond-bending modes. As a result, the electron-phonon coupling is suppressed by $30\%$ ($λ\approx 1.84$ at 200 GPa). Moreover, while at the harmonic level T$_c$ decreases with increasing pressure, the inclusion of anharmonicity leads to a T$_c$ that is almost independent of pressure. High pressure hydrogen sulfide is a strongly anharmonic superconductor.

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
preprint2014arXiv

Anomalous non-additive dispersion interactions in systems of three one-dimensional wires

The non-additive dispersion contribution to the binding energy of three one-dimensional (1D) wires is investigated using wires modelled by (i) chains of hydrogen atoms and (ii) homogeneous electron gases. We demonstrate that the non-additive dispersion contribution to the binding energy is significantly enhanced compared with that expected from Axilrod-Teller-Muto-type triple-dipole summations and follows a different power-law decay with separation. The triwire non-additive dispersion for 1D electron gases scales according to the power law $d^{-β}$, where $d$ is the wire separation, with exponents $β(r_s)$ smaller than 3 and slightly increasing with $r_s$ from 2.4 at $r_s = 1$ to 2.9 at $r_s=10$, where $r_s$ is the density parameter of the 1D electron gas. This is in good agreement with the exponent $β=3$ suggested by the leading-order charge-flow contribution to the triwire non-additivity, and is a significantly slower decay than the $\sim d^{-7}$ behaviour that would be expected from triple-dipole summations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of Smooth Hartree-Fock Pseudopotentials

The accuracy of two widely used scalar-relativistic Hartree-Fock pseudopotentials, the Trail-Needs-Dirac-Fock (TNDF) and the Burkatzki-Filippi-Dolg (BFD) pseudopotentials, is assessed. The performance of the pseudopotentials is tested for a chemically representative set of 34 first-row molecules. All comparisons are made at the Hartree-Fock level of theory, and both sets of pseudopotentials give good results. The all-electron equilibrium geometries, molecular dissociation energies, and zero-point vibrational energies are reproduced a little more accurately by the TNDF pseudopotentials than the BFD ones.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact and approximate Kohn-Sham potentials in ensemble density-functional theory

We construct exact Kohn-Sham potentials for the ensemble density-functional theory (EDFT) from the ground and excited states of helium. The exchange-correlation (XC) potential is compared with the quasi-local-density approximation and both single determinant and symmetry eigenstate ghost-corrected exact exchange approximations. Symmetry eigenstate Hartree-exchange recovers distinctive features of the exact XC potential and is used to calculate the correlation potential. Unlike the exact case, excitation energies calculated from these approximations depend on ensemble weight, and it is shown that only the symmetry eigenstate method produces an ensemble derivative discontinuity. Differences in asymptotic and near-ground-state behavior of exact and approximate XC potentials are discussed in the context of producing accurate optical gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excitations and benchmark ensemble density functional theory for two electrons

A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange (SEHX), is derived. Exact conditions that are proven include the signs of the correlation energy components, the virial theorem for both exchange and correlation, and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-fidelity contact pseudopotentials and p-wave superconductivity

We develop ultratransferable pseudopotentials for the contact interaction that are 100 times more accurate than contemporary approximations. The pseudopotential offers scattering properties very similar to the contact potential, has a smooth profile to accelerate numerics by a factor of up to 4,000, and, for positive scattering lengths, does not support an unwanted bound state. We demonstrate these advantages in a Diffusion Monte Carlo study of fermions with repulsive interactions, delivering the first numerical evidence for the formation of a p-wave superconducting state.

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

Compression algorithm for multi-determinant wave functions

A compression algorithm is introduced for multi-determinant wave functions which can greatly reduce the number of determinants that need to be evaluated in quantum Monte Carlo calculations. We have devised an algorithm with three levels of compression, the least costly of which yields excellent results in polynomial time. We demonstrate the usefulness of the compression algorithm for evaluating multi-determinant wave functions in quantum Monte Carlo calculations, whose computational cost is reduced by factors of between 1.885(3) and 25.23(4) for the examples studied. We have found evidence of sub-linear scaling of quantum Monte Carlo calculations with the number of determinants when the compression algorithm is used.

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

Excitons and biexcitons in symmetric electron-hole bilayers

Symmetric electron-hole bilayer systems have been studied at zero temperature using the diffusion quantum Monte Carlo method. A flexible trial wave function is used that can describe fluid, excitonic and biexcitonic phases. We calculate condensate fractions and pair correlation functions for a large number of densities rs and layer separations d. At small d we find a one-component fluid phase, an excitonic fluid phase, and a biexcitonic fluid phase, and the transitions among them appear to be continuous. At d = 0, excitons appear to survive down to about rs = 0.5 a.u., and biexcitons form at rs > 2.5 a.u.

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

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

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