From plane waves to local Gaussians for the simulation of correlated periodic systems

Can Spin-Component Scaled MP2 Achieve kJ/mol Accuracy for Cohesive Energies of Molecular Crystals?

Attaining kJ/mol accuracy in cohesive energy for molecular crystals is a persistent challenge in computational materials science. In this study, we evaluate second-order Møller-Plesset perturbation theory (MP2) and its spin-component scaled models for calculating cohesive energies for 23 molecular crystals (X23 data set). Using periodic boundary conditions and Brillouin zone sampling, we converge results to the thermodynamic and complete basis set limits, achieving an accuracy of about 2 kJ/mol (0.5 kcal/mol), which is rarely achieved in previous MP2 calculations for molecular crystals. When compared to experimental data, our results have a mean absolute error of 12.9 kJ/mol, comparable to Density Functional Theory with the PBE functional and TS dispersion correction. By separately scaling the opposite-spin and same-spin correlation energy components, using predetermined parameters, we reduce the mean absolute error to 9.5 kJ/mol. Further fine-tuning of these scaling parameters specifically for the X23 data set brings the mean absolute error down to 7.5 kJ/mol.

Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2

We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin-scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators.

Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems

The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H₂O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.

Ground-State Properties of Metallic Solids from Ab Initio Coupled-Cluster Theory

Metallic solids are an enormously important class of materials, but they are a challenging target for accurate wave function-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energy and compare them to experimental values, finding accuracy comparable to common density functionals. Because the quantum chemical "gold standard" CCSD(T) (CCSD with perturbative triple excitations) is inapplicable to metals in the thermodynamic limit, we test two approximate improvements to CCSD, which are found to improve the results.

Integral-Direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal

We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems 1 order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of -72.8 kJ/mol, which is about 10-15 kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals.

Full Configuration Interaction Quantum Monte Carlo treatment of fragments embedded in a periodic mean field

We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach we investigate the performanceof full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be ratherlarge. As a test case we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.

Correlation-Consistent Gaussian Basis Sets for Solids Made Simple

The rapidly growing interest in simulating condensed-phase materials using quantum chemistry methods calls for a library of high-quality Gaussian basis sets suitable for periodic calculations. Unfortunately, most standard Gaussian basis sets commonly used in molecular simulation show significant linear dependencies when used in close-packed solids, leading to severe numerical issues that hamper the convergence to the complete basis set (CBS) limit, especially in correlated calculations. In this work, we revisit Dunning's strategy for construction of correlation-consistent basis sets and examine the relationship between accuracy and numerical stability in periodic settings. We find that limiting the number of primitive functions avoids the appearance of problematic small exponents while still providing smooth convergence to the CBS limit. As an example, we generate double-, triple-, and quadruple-ζ correlation-consistent Gaussian basis sets for periodic calculations with Goedecker-Teter-Hutter (GTH) pseudopotentials. Our basis sets cover the main-group elements from the first three rows of the periodic table. Especially for atoms on the left side of the periodic table, our basis sets are less diffuse than those used in molecular calculations. We verify the fast and reliable convergence to the CBS limit in both Hartree-Fock and post-Hartree-Fock (MP2) calculations, using a diverse test set of 19 semiconductors and insulators.

Optical and Electronic Properties of Organic NIR-II Fluorophores by Time-Dependent Density Functional Theory and Many-Body Perturbation Theory: GW-BSE Approaches

Organic-molecule fluorophores with emission wavelengths in the second near-infrared window (NIR-II, 1000-1700 nm) have attracted substantial attention in the life sciences and in biomedical applications because of their excellent resolution and sensitivity. However, adequate theoretical levels to provide efficient and accurate estimations of the optical and electronic properties of organic NIR-II fluorophores are lacking. The standard approach for these calculations has been time-dependent density functional theory (TDDFT). However, the size and large excitonic energies of these compounds pose challenges with respect to computational cost and time. In this study, we used the GW approximation combined with the Bethe-Salpeter equation (GW-BSE) implemented in many-body perturbation theory approaches based on density functional theory. This method was used to perform calculations of the excited states of two NIR molecular fluorophores (BTC980 and BTC1070), going beyond TDDFT. In this study, the optical absorption spectra and frontier molecular orbitals of these compounds were compared using TDDFT and GW-BSE calculations. The GW-BSE estimates showed excellent agreement with previously reported experimental results.

Efficient Evaluation of Two-Center Gaussian Integrals in Periodic Systems

By using Poisson's summation formula, we calculate periodic integrals over Gaussian basis functions by partitioning the lattice summations between the real and reciprocal space, where both sums converge exponentially fast with a large exponent. We demonstrate that the summation can be performed efficiently to calculate two-center Gaussian integrals over various kernels including overlap, kinetic, and Coulomb. The summation in real space is performed using an efficient flavor of the McMurchie-Davidson recurrence relation. The expressions for performing summation in the reciprocal space are also derived and implemented. The algorithm for reciprocal space summation allows us to reuse several terms and leads to a significant improvement in efficiency when highly contracted basis functions with large exponents are used. We find that the resulting algorithm is only between a factor of 5 and 15 slower than that for molecular integrals, indicating the very small number of terms needed in both the real and reciprocal space summations. An outline of the algorithm for calculating three-center Coulomb integrals is also provided.

DMRG on Top of Plane-Wave Kohn-Sham Orbitals: A Case Study of Defected Boron Nitride

In this paper, we analyze the numerical aspects of the inherent multireference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

Regional Embedding Enables High-Level Quantum Chemistry for Surface Science

Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, this paper introduces the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal-point corrections can provide converged CCSD(T) (coupled-cluster) adsorption energies with very small fragment sizes.

Electronic Structure Modeling of Metal-Organic Frameworks

Owing to their molecular building blocks, yet highly crystalline nature, metal-organic frameworks (MOFs) sit at the interface between molecule and material. Their diverse structures and compositions enable them to be useful materials as catalysts in heterogeneous reactions, electrical conductors in energy storage and transfer applications, chromophores in photoenabled chemical transformations, and beyond. In all cases, density functional theory (DFT) and higher-level methods for electronic structure determination provide valuable quantitative information about the electronic properties that underpin the functions of these frameworks. However, there are only two general modeling approaches in conventional electronic structure software packages: those that treat materials as extended, periodic solids, and those that treat materials as discrete molecules. Each approach has features and benefits; both have been widely employed to understand the emergent chemistry that arises from the formation of the metal-organic interface. This Review canvases these approaches to date, with emphasis placed on the application of electronic structure theory to explore reactivity and electron transfer using periodic, molecular, and embedded models. This includes (i) computational chemistry considerations such as how functional, k-grid, and other model variables are selected to enable insights into MOF properties, (ii) extended solid models that treat MOFs as materials rather than molecules, (iii) the mechanics of cluster extraction and subsequent chemistry enabled by these molecular models, (iv) catalytic studies using both solids and clusters thereof, and (v) embedded, mixed-method approaches, which simulate a fraction of the material using one level of theory and the remainder of the material using another dissimilar theoretical implementation.

Featured properties of Li+-based battery anode: Li4Ti5O12

3D ternary Li4Ti5O12, a Li+-based battery anode, presents an unusual lattice symmetry (triclinic crystal), band structure, charge density, and density of states under first-principles calculations. It is a large direct-gap semiconductor with E d g ∼ 2.98 eV. The atom-dominated valence and conduction bands, the spatial charge distribution and the atom- and orbital-decomposed van Hove singularities are available for delicate identifications of multi-orbital hybridizations in Li-O and Ti-O bonds. The extremely non-uniform chemical environment, which induces very complicated hopping integrals, directly arises from the large bonding fluctuations and the highly anisotropic configurations. Also, the developed theoretical framework is very useful for fully understanding cathodes and electrolytes of oxide compounds.

Pseudodiagonalization-based wavefunction optimization with contracted planewave basis functions

Electronic structure calculations representing the molecular orbitals (MOs) with contracted planewave basis functions (CPWBFs) have been reported recently. CPWBFs are Fourier-series representations of atom-centered basis functions. The mathematical features of CPWBFs permit the construction of matrix-vector products, FC _o , involving the application of the Fock matrix, F, to the set of occupied MOs, C _o , without the explicit evaluation of F. This approach offers a theoretical speed-up of M/n over F-based methods, where M and n are the number of basis functions and occupied MOs, respectively. The present study reports methodological advances that permit FC _o -based optimization of wavefunction formed from CPWBFs. In particular, a technique is reported for optimizing wavefunctions by combining pseudodiagonalization techniques based on an exact representation of FC _o , approximate information regarding the virtual orbital energies, and direct inversion of the iterative subspace optimization schemes to guide the wavefunction to a converged solution. This method is found to speed-up wavefunction optimizations by factors of up to ~6 - 8 over F-based optimization methods while providing identical results. Further, the computational cost of this technique is relatively insensitive to basis set size, thus providing further benefits in calculations using large CPWBF basis sets. The results of density functional theory calculations show that this method permits the use of hybrid exchange-correlation (XC) functionals with a small increase in effort over analogous calculations using generalized gradient approximation XC functionals. © 2019 Wiley Periodicals, Inc.

Duality of Ring and Ladder Diagrams and Its Importance for Many-Electron Perturbation Theories

We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the contribution of the ladder diagrams to the electronic correlation energy can be approximated in an effective manner using second-order perturbation theory. We employ the latter approximation to reduce the computational cost of coupled cluster theory calculations for insulators and semiconductors by 2 orders of magnitude without compromising accuracy.

An optimized twist angle to find the twist-averaged correlation energy applied to the uniform electron gas

We explore an alternative to twist averaging in order to obtain more cost-effective and accurate extrapolations to the thermodynamic limit (TDL) for coupled cluster doubles (CCD) calculations. We seek a single twist angle to perform calculations at, instead of integrating over many random points or a grid. We introduce the concept of connectivity, a quantity derived from the nonzero four-index integrals in an MP2 calculation. This allows us to find a special twist angle that provides appropriate connectivity in the energy equation, which yields results comparable to full twist averaging. This special twist angle effectively makes the finite electron number CCD calculation represent the TDL more accurately, reducing the cost of twist-averaged CCD over N_s twist angles from N_s CCD calculations to N_s MP2 calculations plus one CCD calculation.

Physisorption of Water on Graphene: Subchemical Accuracy from Many-Body Electronic Structure Methods

Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.

Understanding non-covalent interactions in larger molecular complexes from first principles

Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.

Direct inversion of the iterative subspace with contracted planewave basis functions

Ways to reduce the computational cost of periodic electronic structure calculations by using basis functions corresponding to linear combinations of planewaves have been examined recently. These contracted planewave (CPW) basis functions correspond to Fourier series representations of atom-centered basis functions, and thus provide access to some beneficial properties of planewave (PW) and localized basis functions. This study reports the development and assessment of a direct inversion of the iterative subspace (DIIS) method that employs unique properties of CPW basis functions to efficiently converge electronic wavefunctions. This method relies on access to a PW-based representation of the electronic structure to provide a means of efficiently evaluating matrix-vector products involving the application of the Fock matrix to the occupied molecular orbitals. These matrix-vector products are transformed into a form permitting the use of direct diagonalization techniques and DIIS methods typically employed with atom-centered basis sets. The abilities of this method are assessed through periodic Hartree-Fock calculations of a range of molecules and solid-state systems. The results show that the method reported in this study is approximately five times faster than CPW-based calculations in which the entire Fock matrix is calculated. This method is also found to be weakly dependent upon the size of the basis set, thus permitting the use of larger CPW basis sets to increase variational flexibility with a minor impact on computational performance. © 2018 Wiley Periodicals, Inc.

Communication: A novel implementation to compute MP2 correlation energies without basis set superposition errors and complete basis set extrapolation

By using a formulation based on the dynamical polarizability, we propose a novel implementation of second-order Møller-Plesset perturbation (MP2) theory within a plane wave (PW) basis set. Because of the intrinsic properties of PWs, this method is not affected by basis set superposition errors. Additionally, results are converged without relying on complete basis set extrapolation techniques; this is achieved by using the eigenvectors of the static polarizability as an auxiliary basis set to compactly and accurately represent the response functions involved in the MP2 equations. Summations over the large number of virtual states are avoided by using a formalism inspired by density functional perturbation theory, and the Lanczos algorithm is used to include dynamical effects. To demonstrate this method, applications to three weakly interacting dimers are presented.

A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface

We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.

Projector Quantum Monte Carlo Method for Nonlinear Wave Functions

We reformulate the projected imaginary-time evolution of the full configuration interaction quantum Monte Carlo method in terms of a Lagrangian minimization. This naturally leads to the admission of polynomial complex wave function parametrizations, circumventing the exponential scaling of the approach. While previously these functions have traditionally inhabited the domain of variational Monte Carlo approaches, we consider recent developments for the identification of deep-learning neural networks to optimize this Lagrangian, which can be written as a modification of the propagator for the wave function dynamics. We demonstrate this approach with a form of tensor network state, and use it to find solutions to the strongly correlated Hubbard model, as well as its application to a fully periodic ab initio graphene sheet. The number of variables which can be simultaneously optimized greatly exceeds alternative formulations of variational Monte Carlo methods, allowing for systematic improvability of the wave function flexibility towards exactness for a number of different forms, while blurring the line between traditional variational and projector quantum Monte Carlo approaches.

Exact exchange with non-orthogonal generalized Wannier functions

The evaluation of exact exchange (EXX) is an important component of quantum chemical calculations performed with ab initio and hybrid density functional methods. While evaluating exact exchange is routine in molecular quantum chemical calculations performed with localized basis sets, the non-local nature of the exchange operator presents a major impediment to the efficient use of exact exchange in calculations that employ planewave basis sets. Non-orthogonal generalized Wannier functions (NGWFs) corresponding to planewave expansions of localized basis functions are an alternative form of basis set that can be used in quantum chemical calculations. The periodic nature of these functions renders them suitable for calculations of periodic systems, while the contraction of sets of planewaves into individual basis functions reduces the number of variational parameters, permitting the construction and direct diagonalization of the Fock matrix. The present study examines how NGWFs corresponding to Fourier series representations of conventional atom-centered basis sets can be used to evaluate exact exchange in periodic systems. Specifically, an approach for constructing the exchange operator with NGWFs is presented and used to perform Hartree-Fock calculations with a series of molecules in periodically repeated simulation cells. The results demonstrate that the NGWF approach is significantly faster than the EXX method, which is a standard approach for evaluating exact exchange in periodic systems.