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). Quantum 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 inter-operable 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.Comment: 36 pages, 5 figures, resubmitted to J.Phys.: Condens. Matte

Quantum Monte Carlo for large chemical systems: Implementing efficient strategies for petascale platforms and beyond

Various strategies to implement efficiently QMC simulations for large chemical systems are presented. These include: i.) the introduction of an efficient algorithm to calculate the computationally expensive Slater matrices. This novel scheme is based on the use of the highly localized character of atomic Gaussian basis functions (not the molecular orbitals as usually done), ii.) the possibility of keeping the memory footprint minimal, iii.) the important enhancement of single-core performance when efficient optimization tools are employed, and iv.) the definition of a universal, dynamic, fault-tolerant, and load-balanced computational framework adapted to all kinds of computational platforms (massively parallel machines, clusters, or distributed grids). These strategies have been implemented in the QMC=Chem code developed at Toulouse and illustrated with numerical applications on small peptides of increasing sizes (158, 434, 1056 and 1731 electrons). Using 10k-80k computing cores of the Curie machine (GENCI-TGCC-CEA, France) QMC=Chem has been shown to be capable of running at the petascale level, thus demonstrating that for this machine a large part of the peak performance can be achieved. Implementation of large-scale QMC simulations for future exascale platforms with a comparable level of efficiency is expected to be feasible

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post–Hartree–Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension

The calculation of optical absorption spectra using linear-scaling density-functional theory

The goal of the work presented in this thesis was to develop and implement a method for calculating optical absorption spectra for large electronic systems within a linear-scaling density-functional theory (LS-DFT) formalism. The key feature of this method was the development of a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states, which are not well represented by the localized orbital basis sets used for ground state LS-DFT calculations. Three different schemes were compared for the calculation of unoccupied states using a one-dimensional “toy model” and the most promising of these, based on the use of a projection operator, was implemented in a fully-functional LS-DFT code. Using the toy model, two methods for the calculation of band structures within a localized basis set were investigated and some of the features required by localized basis sets in order to produce accurate band structures were identified. The method was tested by the application to both molecular and extended systems, with calculations of densities of states, band structures and optical absorption spectra. The results for the smaller systems were validated by comparison with a cubic-scaling plane-wave density-functional theory code, with which excellent agreement was achieved. Additionally, the method was shown to be linear-scaling for a conjugated polymer for system sizes up to 1000 atoms. The use of the projection method was shown to be crucial for calculating the above results, as was the implementation of a momentum operator based formalism for the calculation of spectra. Finally, it was shown that the method can be used to identify the transitions responsible for particular peaks in the spectra and is sensitive enough to distinguish between spectra for systems with very similar structures, demonstrating the capabilities of the method to aid the interpretation of experimental results

Feet on the potential energy surface, head in the π clouds

The landscape of a potential energy surface is marked by chemically interesting features. Hills and valleys correspond to transition states and reactive intermediates; the deepest valley gives the most stable configuration. Mapping these features for individual molecules and for the interactions between molecules is one of the goals of computational chemistry. The dispersion energy is a weak attractive force in intermolecular interactions. Dispersion energy results from a purely quantum mechanical effect, in which instantaneous multipoles on one molecule induce multipoles on another. Among neutral atoms or molecules that lack permanent multipole moments, the dispersion interaction is the principal attractive force. Dispersion also plays a significant role in the interaction between molecules with diffuse π clouds. This interaction is often difficult to capture with standard computational chemistry methods, so a comparison of the results obtained with various methods is itself important. This work presents explorations of the potential energy surface of clusters of atoms and of the interactions between molecules. First, structures of small aluminum clusters are examined and classified as ground states, transition states, or higher-order saddle points. Subsequently, the focus shifts to dispersion-dominated π-π interactions when the potential energy surfaces of benzene, substituted benzene, and pyridine dimers are explored. Because DNA nucleotide bases can be thought of as substituted heterocycles, a natural extension of the substituted benzene and pyridine investigations is to model paired nucleotide bases. Finally, the success of the dispersion studies inspires the development of an extension to the computational method used, which will enable the dispersion energy to be modeled - and the potential energy surface explored - in additional chemical systems