Advanced capabilities for materials modelling with Q <scp>uantum</scp> ESPRESSO

P Giannozzi; O Andreussi; T Brumme; O Bunau; M Buongiorno Nardelli; M Calandra; R Car; C Cavazzoni; D Ceresoli; M Cococcioni; N Colonna; I Carnimeo; A Dal Corso; S de Gironcoli; P Delugas; R A DiStasio; A Ferretti; A Floris; G Fratesi; G Fugallo; R Gebauer; U Gerstmann; F Giustino; T Gorni; J Jia; M Kawamura; H-Y Ko; A Kokalj; E Küçükbenli; M Lazzeri; M Marsili; N Marzari; F Mauri; N L Nguyen; H-V Nguyen; A Otero-de-la-Roza; L Paulatto; S Poncé; D Rocca; R Sabatini; B Santra; M Schlipf; A P Seitsonen; A Smogunov; I Timrov; T Thonhauser; P Umari; N Vast; X Wu; S Baroni

doi:10.1088/1361-648x/aa8f79

Abstract

Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

Research topics

Advanced Chemical Physics Studies Machine Learning in Materials Science Physics of Superconductivity and Magnetism

Full text

Read the research paper

PDF provided by the publishing journal.

View PDF Download PDF

Continue discovering

Related research

Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen

In the past, basis sets for use in correlated molecular calculations have largely been taken from single configuration calculations. Recently, Almlöf, Taylor, and co-workers have found that basis sets of natural orbitals derived from correlated atomic calculations (ANOs) provide an excellent description of molecular correlation effects. We report here a careful study of correlation effects in the oxygen atom, establishing that compact sets of primitive Gaussian functions effectively and efficiently describe correlation effects if the exponents of the functions are optimized in atomic correlated calculations, although the primitive (sp) functions for describing correlation effects can be taken from atomic Hartree–Fock calculations if the appropriate primitive set is used. Test calculations on oxygen-containing molecules indicate that these primitive basis sets describe molecular correlation effects as well as the ANO sets of Almlöf and Taylor. Guided by the calculations on oxygen, basis sets for use in correlated atomic and molecular calculations were developed for all of the first row atoms from boron through neon and for hydrogen. As in the oxygen atom calculations, it was found that the incremental energy lowerings due to the addition of correlating functions fall into distinct groups. This leads to the concept of correlation consistent basis sets, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation consistent sets are given for all of the atoms considered. The most accurate sets determined in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding ANO sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estimated that this set yields 94%–97% of the total (HF+1+2) correlation energy for the atoms neon through boron.

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions.

Accurate and simple analytic representation of the electron-gas correlation energy

We propose a simple analytic representation of the correlation energy ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$ for a uniform electron gas, as a function of density parameter ${\mathit{r}}_{\mathit{s}}$ and relative spin polarization \ensuremath{\zeta}. Within the random-phase approximation (RPA), this representation allows for the ${\mathit{r}}_{\mathit{s}}^{\mathrm{\ensuremath{-}}3/4}$ behavior as ${\mathit{r}}_{\mathit{s}}$\ensuremath{\rightarrow}\ensuremath{\infty}. Close agreement with numerical RPA values for ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$,0), ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$,1), and the spin stiffness ${\mathrm{\ensuremath{\alpha}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$)=${\mathrm{\ensuremath{\partial}}}^{2}$${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$, \ensuremath{\zeta}=0)/\ensuremath{\delta}${\mathrm{\ensuremath{\zeta}}}^{2}$, and recovery of the correct ${\mathit{r}}_{\mathit{s}}$ln${\mathit{r}}_{\mathit{s}}$ term for ${\mathit{r}}_{\mathit{s}}$\ensuremath{\rightarrow}0, indicate the appropriateness of the chosen analytic form. Beyond RPA, different parameters for the same analytic form are found by fitting to the Green's-function Monte Carlo data of Ceperley and Alder [Phys. Rev. Lett. 45, 566 (1980)], taking into account data uncertainties that have been ignored in earlier fits by Vosko, Wilk, and Nusair (VWN) [Can. J. Phys. 58, 1200 (1980)] or by Perdew and Zunger (PZ) [Phys. Rev. B 23, 5048 (1981)]. While we confirm the practical accuracy of the VWN and PZ representations, we eliminate some minor problems with these forms. We study the \ensuremath{\zeta}-dependent coefficients in the high- and low-density expansions, and the ${\mathit{r}}_{\mathit{s}}$-dependent spin susceptibility. We also present a conjecture for the exact low-density limit. The correlation potential ${\mathrm{\ensuremath{\mu}}}_{\mathit{c}}^{\mathrm{\ensuremath{\sigma}}}$(${\mathit{r}}_{\mathit{s}}$,\ensuremath{\zeta}) is evaluated for use in self-consistent density-functional calculations.