QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

Superionic and metallic states of water and ammonia at giant planet conditions

The phase diagrams of water and ammonia were determined by constant pressure ab initio molecular dynamic simulations at pressures (30 to 300 gigapascal) and temperatures (300 to 7000 kelvin) of relevance for the middle ice layers of the giant planets Neptune and Uranus. Along the planetary isentrope water and ammonia behave as fully dissociated ionic, electronically insulating fluid phases, which turn metallic at temperatures exceeding 7000 kelvin for water and 5500 kelvin for ammonia. At lower temperatures, the phase diagrams of water and ammonia exhibit a superionic solid phase between the solid and the ionic liquid. These simulations improve our understanding of the properties of the middle ice layers of Neptune and Uranus.

Physics of iron at Earth's core conditions

The bulk properties of iron at the pressure and temperature conditions of Earth's core were determined by a method that combines first-principles and classical molecular dynamic simulations. The theory indicates that (i) the iron melting temperature at inner-core boundary (ICB) pressure (330 gigapascals) is 5400 (+/-400) kelvin; (ii) liquid iron at ICB conditions is about 6% denser than Earth's outer core; and (iii) the shear modulus of solid iron close to its melting line is 140 gigapascals, consistent with the seismic value for the inner core. These results reconcile melting temperature estimates based on sound velocity shock wave data with those based on diamond anvil cell experiments.

Liquid-liquid phase transition in compressed hydrogen from first-principles simulations

The properties of compressed liquid hydrogen, the most abundant fluid in the universe, have been investigated by means of first-principles molecular dynamics at pressures between 75 and 175 GPa and temperatures closer to the freezing line than so far reported in shock-wave experiments. Evidence for a liquid-liquid transition between a molecular and a dissociated phase is provided. The transition is accompanied by a 6% increase in density and by metallization. This finding has important implications for our understanding of the interiors of giant planets and supports predictions of a quantum fluid state at low temperatures.

Synthesis of novel transition metal nitrides IrN2 and OsN2

Two new transition metal nitrides, IrN2 and OsN2, were synthesized at high pressures and temperatures using laser-heated diamond-anvil cell techniques. Synchrotron x-ray diffraction was used to determine the structures of novel nitrides and the equations of states of both the parent metals as well as the newly synthesized materials. The compounds have bulk moduli comparable with those of the traditional superhard materials. For IrN2, the measured bulk modulus [K0 = 428(12) GPa] is second only to that of diamond (K0 = 440 GPa). Ab initio calculations indicate that both compounds have a metal:nitrogen stoichiometry of 1:2 and that nitrogen intercalates in the lattice of the parent metal in the form of single-bonded N-N units.

Amorphous silica-like carbon dioxide

Among the group IV elements, only carbon forms stable double bonds with oxygen at ambient conditions. At variance with silica and germania, the non-molecular single-bonded crystalline form of carbon dioxide, phase V, only exists at high pressure. The amorphous forms of silica (a-SiO2) and germania (a-GeO2) are well known at ambient conditions; however, the amorphous, non-molecular form of CO2 has so far been described only as a result of first-principles simulations. Here we report the synthesis of an amorphous, silica-like form of carbon dioxide, a-CO2, which we call 'a-carbonia'. The compression of the molecular phase III of CO2 between 40 and 48 GPa at room temperature initiated the transformation to the non-molecular amorphous phase. Infrared spectra measured at temperatures up to 680 K show the progressive formation of C-O single bonds and the simultaneous disappearance of all molecular signatures. Furthermore, state-of-the-art Raman and synchrotron X-ray diffraction measurements on temperature-quenched samples confirm the amorphous character of the material. Comparison with vibrational and diffraction data for a-SiO2 and a-GeO2, as well as with the structure factor calculated for the a-CO2 sample obtained by first-principles molecular dynamics, shows that a-CO2 is structurally homologous to the other group IV dioxide glasses. We therefore conclude that the class of archetypal network-forming disordered systems, including a-SiO2, a-GeO2 and water, must be extended to include a-CO2.

First-principles codes for computational crystallography in the Quantum-ESPRESSO package

The Quantum-ESPRESSO package is a multi-purpose and multi-platform software for ab-initio calculations of condensed matter (periodic and disordered) systems. Codes in the package are based on density functional theory and on a plane wave/pseudopotential description of the electronic ground state and are ideally suited for structural optimizations (both at zero and at finite temperature), linear response calculations (phonons, elastic constants, dielectric and Raman tensors, etc.) and high-temperature molecular dynamics. Examples of applications of the codes included in the package are briefly discussed.

Dissociation of methane into hydrocarbons at extreme (planetary) pressure and temperature

Constant-pressure, first-principles molecular dynamic simulations were used to investigate the behavior of methane at high pressure and temperature. Contrary to the current interpretation of shock-wave experiments, the simulations suggest that, below 100 gigapascals, methane dissociates into a mixture of hydrocarbons, and it separates into hydrogen and carbon only above 300 gigapascals. The simulation conditions (100 to 300 gigapascals; 4000 to 5000 kelvin) were chosen to follow the isentrope in the middle ice layers of Neptune and Uranus. Implications on the physics of these planets are discussed.

Pressure-induced solid carbonates from molecular CO2 by computer simulation

A combination of ab initio molecular dynamic simulations and fully relaxed total energy calculations is used to predict that molecular CO2 should transform to nonmolecular carbonate phases based on CO4 tetrahedra at pressures in the range of 35 to 60 gigapascals. The simulation suggests a variety of competing phases, with a more facile transformation of the molecular phase at high temperatures. Thermodynamically, the most stable carbonate phase at high pressure is predicted to be isostructural to SiO2 alpha-quartz (low quartz). A class of carbonates, involving special arrangements of CO4 tetrahedra, is found to be more stable than all the other silica-like polymorphs.

High-pressure polymeric phases of carbon dioxide

Understanding the structural transformations of solid CO(2) from a molecular solid characterized by weak intermolecular bonding to a 3-dimensional network solid at high pressure has challenged researchers for the past decade. We employ the recently developed metadynamics method combined with ab initio calculations to provide fundamental insight into recent experimental reports on carbon dioxide in the 60-80 GPa pressure region. Pressure-induced polymeric phases and their transformation mechanisms are found. Metadynamics simulations starting from the CO(2)-II (P4(2)/mnm) at 60 GPa and 600 K proceed via an intermediate, partially polymerized phase, and finally yield a fully tetrahedral, layered structure (P-4m2). Based on the agreement between calculated and experimental Raman and X-ray patterns, the recently identified phase VI [Iota V, et al. (2007) Sixfold coordinated carbon dioxide VI. Nature Mat 6:34-38], assumed to be disordered stishovite-like, is instead interpreted as the result of an incomplete transformation of the molecular phase into a final layered structure. In addition, an alpha-cristobalite-like structure (P4(1)2(1)2), is predicted to be formed from CO(2)-III (Cmca) via an intermediate Pbca structure at 80 GPa and low temperatures (<300 K). Defects in the crystals are frequently observed in the calculations at 300 K whereas at 500 to 700 K, CO(2)-III transforms to an amorphous form, consistent with experiment [Santoro M, et al. (2006) Amorphous silica-like carbon dioxide. Nature 441:857-860], but the simulation yields additional structural details for this disordered solid.

Thermal conductivity of solid argon from molecular dynamics simulations

The thermal conductivity of solid argon in the classical limit has been calculated by equilibrium molecular dynamic simulations using the Green-Kubo formalism and a Lennard-Jones interatomic potential. Contrary to previous theoretical reports, we find that the computed thermal conductivities are in good agreement with experimental data. The computed values are also in agreement with the high-temperature limit of the three-phonon scattering contribution to the thermal conductivity. We find that finite-size effects are negligible and that phonon lifetimes have two characteristic time scales, so that agreement with kinetic theory is obtained only after appropriate averaging of the calculated phonon lifetimes.

Structure of a CH3S monolayer on Au(111) solved by the interplay between molecular dynamics calculations and diffraction measurements

We have investigated the controversy surrounding the (sqrt[3] x sqrt[3]) R30 degrees structure of self-assembled monolayers of methylthiolate on Au(111) by first principles molecular dynamics simulations, energy and angle resolved photoelectron diffraction, and grazing incidence x-ray diffraction. Our simulations find a dynamic equilibrium between bridge site adsorption and a novel structure where 2 CH3S radicals are bound to an Au adatom that has been lifted from the gold substrate. As a result, the interface is characterized by a large atomic roughness with both adatoms and vacancies. This result is confirmed by extensive photoelectron and grazing incidence x-ray diffraction measurements.

Carbon phase diagram from ab initio molecular dynamics

We compute the free energy of solid and liquid diamond from first-principles electronic structure theory using efficient thermodynamic integration techniques. Our calculated melting curve is in excellent agreement with the experimental estimate of the graphite-diamond-liquid triple point and is consistent with shock wave experiments. We predict the phase diagram of diamond at pressures and temperatures that are difficult to access experimentally. We confirm early speculations on the presence of a reentrant point in the diamond melting line but find no evidence for a first order liquid-liquid phase transition near the reentrant point.

Metal work-function changes induced by organic adsorbates: a combined experimental and theoretical study

The role of molecular dipole moment, charge transfer, and Pauli repulsion in determining the work-function change (Deltaphi) at organic-metal interfaces has been elucidated by a combined experimental and theoretical study of (CH(3)S)(2)/Au(111) and CH(3)S/Au(111). Comparison between experiment and theory allows us to determine the origin of the interface dipole layer for both phases. For CH(3)S/Au(111), Deltaphi can be ascribed almost entirely to the dipole moment of the CH(3)S layer. For (CH(3)S)(2)/Au(111), a Pauli repulsion mechanism occurs. The implications of these results on the interpretation of Deltaphi in the presence of strongly and weakly adsorbed molecules is discussed.

Infrared and Raman spectra of silica polymorphs from anab initioparametrized polarizable force field

The general aim of this study is to test the reliability of polarizable model potentials for the prediction of vibrational (infrared and Raman) spectra in highly anharmonic systems such as high temperature crystalline phases. By using an ab initio parametrized interatomic potential for SiO2 and molecular dynamics simulations, we calculate the infrared and Raman spectra for quartz, cristobalite, and stishovite at various thermodynamic conditions. The model is found to perform very well in the prediction of infrared spectra. Raman peak positions are also reproduced very well by the model; however, Raman intensities calculated by explicitly taking the derivative of the polarizability with respect to the atomic displacements are found to be in poorer agreement than intensities calculated using a parametrized "bond polarizability" model. Calculated spectra for the high temperature beta phases, where the role of dynamical disorder and anharmonicities is predominant, are found to be in excellent agreement with experiments. For the octahedral phases, our simulations are able to reproduce changes in the Raman spectra across the rutile-to-CaCl2 transition around 50 GPa, including the observed phonon softening.

Structure and Dynamics of Low-Density and High-Density Liquid Water at High Pressure

Liquid water has a primary role in ruling life on Earth in a wide temperature and pressure range as well as a plethora of chemical, physical, geological, and environmental processes. Nevertheless, a full understanding of its dynamical and structural properties is still lacking. Water molecules are associated through hydrogen bonds, with the resulting extended network characterized by a local tetrahedral arrangement. Two different local structures of the liquid, called low-density (LDW) and high-density (HDW) water, have been identified to potentially affect many different chemical, biological, and physical processes. By combining diamond anvil cell technology, ultrafast pump-probe infrared spectroscopy, and classical molecular dynamics simulations, we show that the liquid structure and orientational dynamics are intimately connected, identifying the P-T range of the LDW and HDW regimes. The latter are defined in terms of the speeding up of the orientational dynamics, caused by the increasing probability of breaking and reforming the hydrogen bonds.

Trapping of excitons at chemical defects in polyethylene

In a previous paper we studied an injected electron-hole pair in crystalline polyethylene (PE) and found that the exciton becomes weakly self-trapped in a narrow interchain pocket comprised between two gauche defects. Despite the large energy stored in the trapped excitation, there did not appear to be a direct nonradiative channel for electron-hole recombination. Actual polyethylene systems of practical use are, however, neither crystalline nor pure. To understand the fate of an electron-hole pair in the impure case, we studied by ab initio simulations the evolution of an exciton trapped on three common chemical defects found in polyethylene: a grafted carbonyl (C=O); an intrachain vinyl group (C=C); a grafted carboxyl (COOH). Ab initio simulations lead to predict three different outcomes: trapping, nonradiative recombination, and homolitic bond-breaking, respectively. This suggests that extrinsic self-trapping of electron-hole pairs over chemical defects inside the quasicrystalline fraction of PE could be relevant for electrical damage in high-voltage cables.

Ab initio theory of planetary materials

Ab initio simulations play an increasingly important role in the studies of deep planetary interiors. Here we review the current state of this field, concentrating on studies of the materials of the Earth’s deep interior (MgO—SiO2—FeO—Al2O3, Fe—Si—S—O) and of the interiors of giant planets (H—He system, H2O—CH4—NH3 system). In particular, novel phases and phase diagrams, insights into structural and electronic phase transitions, melting curves, thermoelasticity and the effects of impurities on physical properties of planet-forming materials are discussed.