Sandro Scandolo

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.

X-ray Diffraction and Computation Yield the Structure of Alkanethiols on Gold(111)

The structure of self-assembled monolayers (SAMs) of long-chain alkyl sulfides on gold(111) has been resolved by density functional theory–based molecular dynamics simulations and grazing incidence x-ray diffraction for hexanethiol and methylthiol. The analysis of molecular dynamics trajectories and the relative energies of possible SAM structures suggest a competition between SAM ordering, driven by the lateral van der Waals interaction between alkyl chains, and disordering of interfacial Au atoms, driven by the sulfur-gold interaction. We found that the sulfur atoms of the molecules bind at two distinct surface sites, and that the first gold surface layer contains gold atom vacancies (which are partially redistributed over different sites) as well as gold adatoms that are laterally bound to two sulfur atoms.

An ab initio parametrized interatomic force field for silica

We present a classical interatomic force field for liquid SiO2 which has been parametrized using the forces, stresses and energies extracted from ab initio calculations. We show how inclusion of more electronic effects in a phenomenological way and parametrization at the relevant conditions of pressure and temperature allow the creation of more accurate force fields. We compare the results of simulations with this force field both to experiment and to the results of ab initio molecular dynamics simulations and show how our procedure leads to comparisons which are greatly improved with respect to the most widely used force fields for silica.

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.

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.

How well do Car–Parrinello simulations reproduce the Born–Oppenheimer surface? Theory and examples

We derive an analytic expression for the average difference between the forces on the ions in a Car–Parrinello simulation and the forces obtained at the same ionic positions when the electrons are at their ground state. We show that for common values of the fictitious electron mass, a systematic bias may affect the Car–Parrinello forces in systems where the electron–ion coupling is large. We show that in the limit where the electronic orbitals are rigidly dragged by the ions the difference between the two dynamics amounts to a rescaling of the ionic masses, thereby leaving the thermodynamics intact. We study the examples of crystalline magnesium oxide and crystalline and molten silicon. We find that for crystalline silicon the errors are very small. For crystalline MgO the errors are very large but the dynamics can be quite well corrected within the rigid-ion model. We conclude that it is important to control the effect of the electron mass parameter on the quantities extracted from Car–Parrinello simulat...

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

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

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.

OsN2: Crystal structure and electronic properties

Osmium nitride belongs to a family of nitrides synthesized recently at high pressures from their parent elements. Here we show, based on first-principles calculations, that the crystal structure of osmium nitride is isostructural to marcasite. Excellent agreement is obtained between the authors’ results and x-ray, Raman, and compressibility measurements. In the OsN2 marcasite structure single-bonded N2 units occupy the interstitial sites of the Os close-packed lattice, giving rise to a metallic compound. A comparison between the formation energies of OsN2 and PtN2 explains the similar thermodynamic conditions of formation reported experimentally for the two compounds.

First-principle-constant pressure molecular dynamics

Abstract We present a new method for first-principles numerical simulation of solid-solid phase transformation. The method is applied to the study of pressure induced transformations in silicon and carbon.