Search

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., Comment: 6 figures, plus on supplementary information

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.

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., Comment: 5 pages, 4 figures. (jiansun@nju.edu.cn; miguel.c.martinez@ucl.ac.uk)

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., Comment: 16 pages, 4 figures

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., Comment: 9 pages, 10 figures, 1 table

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., Comment: (5 pages, 4 figures)

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., Comment: 10 pages, 9 figures

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., Comment: 5 pages, 4 figures. (jian.sun@theochem.rub.de)

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., Comment: accepted in Phys. Rev. Lett.; 11 pages, 4 figures; includes supplementary information (9 pages, 4 figures)

First-principles density-functional-theory calculations show that compression of alkali metals stabilizes open structures with localized interstitial electrons which may exhibit a Stoner-type instability towards ferromagnetism. We find ferromagnetic phases of the lithium-IV-type, simple cubic, and simple hexagonal structures in the heavier alkali metals, which may be described as s-band ferromagnets. We predict that the most stable phases of potassium at low temperatures and pressures around 20 GPa are ferromagnets., Comment: 5 pages, 3 figures

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., Comment: 44 pages, 23 figures

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., Comment: 9 pages, 13 figures

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., Comment: 6 pages, 6 figures, 1 table

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., Comment: 7 pages, 8 figures

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has been applied to verify measured J-coupling in a silicophosphate structure, Si5O(PO4)6, Comment: 9 pages

Motivated by recent experimental work, we use first-principles density functional theory methods to conduct an extensive search for low enthalpy structures of C$_6$Ca under pressure. As well as a range of buckled structures, which are energetically competitive over an intermediate range of pressures, we show that the high pressure system ($\gtrsim 18$ GPa) is unstable towards the formation of a novel class of layered structures, with the most stable compound involving carbon sheets containing five- and eight-membered rings. As well as discussing the energetics of the different classes of low enthalpy structures, we comment on the electronic structure of the high pressure compound and its implications for superconductivity., Comment: 4 pages 4 figures

It has been suggested that hydrogen may metallise at lower pressures if it is ``precompressed''. Here we introduce a search strategy for predicting high-pressure structures and apply it to silane using first-principles electronic structure computations. It is based on relaxing randomly chosen structures, and is demonstrated to work well for unit cells containing up to at least ten atoms. We predict that silane will metallise at higher pressures than previously anticipated, but we suggest that the metallic phase might show high-temperature superconductivity at experimentally accessible pressures., Comment: 5 pages, 3 figures and 1 table

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\pi$--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., Comment: 4 Figures. Please see accompanying experimental manuscript "Superconductivity in the Intercalated Graphite Compounds C6Yb and C6Ca" by Weller et al

We show how to describe the coupling of electrons to non-uniform magnetic fields in the framework of the widely used norm-conserving pseudopotential appro ximation for electronic structure calculations. Our derivation applies to magnetic fields that are smooth on the scale of the core region. The method is validated by application to the calculation of the magnetic susceptibility of molecules. Our results are compared with high quality all electron quantum chemical results, and another recently proposed formalism., Comment: 4 pages, submitted to Physical Review Letters

We present a method for calculating the kinetic energy of localised functions represented on a regular real space grid. This method uses fast Fourier transforms applied to restricted regions commensurate with the simulation cell and is applicable to grids of any symmetry. In the limit of large systems it scales linearly with system size. Comparison with the finite difference approach shows that our method offers significant improvements in accuracy without loss of efficiency., Comment: Standard Latex article format, 16 pages, 3 figures. Published in Computer Physics Communications vol.140 (2001) pp. 315-322

A theory for the reliable prediction of the EPR g-tensor for paramagnetic defects in solids is presented. It is based on density functional theory and on the gauge including projector augmented wave (GIPAW) approach to the calculation of all-electron magnetic response. The method is validated by comparison with existing quantum chemical and experimental data for a selection of diatomic radicals. We thenperform the first prediction of EPR ${\rm g} $-tensors in the solid state and find the results to be in excellent agreement with experiment for the $E'_1$ and substitutional P defect centers in quartz., Comment: 5 pages, 4 tables

A comparison between experimental and theoretical electron Energy Loss Near Edge Structure (ELNES) of B and N K-edges in cubic boron nitride is presented. The electron energy loss spectra of cubic boron nitride particles were measured using a scanning transmission electron microscope. The theoretical calculation of the ELNES was performed within the framework of density functional theory including single particle core-hole effects. It is found that experimental and calculated ELNES of both the B and N K-edges in cubic boron nitride show excellent agreement., Comment: 9 pages, 1 figure and 1 table. Submitted to Physical Review B

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., Comment: 25 pages, 4 tables, submitted to Physical Review B