Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost

We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.

Benchmark accuracy of predicted NMR observables for quadrupolar 14 N and 17 O nuclei in molecular crystals

Nuclear quadrupole resonances for 14 N and ¹⁷ O nuclei are exquisitely sensitive to interactions with surrounding atoms. As a result, nitrogen and oxygen solid-state nuclear magnetic resonance (ssNMR) provides a powerful tool for investigating structure and dynamics in complex systems. First-principles calculations are increasingly used to facilitate spectral assignment and to evaluate and adjust crystal structures. Recent work combining the strengths of planewave density functional theory (DFT) calculations with a single molecule correction obtained using a higher level of theory has proven successful in improving the accuracy of predicted chemical shielding (CS) tensors and ¹⁷ O quadrupolar coupling constants ( C q ). Here we extend this work by examining the accuracy of predicted ¹⁴ N and ¹⁷ O electric field gradient (EFG) tensor components obtained using alternative planewave-corrections involving cluster and two-body fragment-based calculations. We benchmark the accuracy of CS and EFG tensor predictions for both nitrogen and oxygen using planewave, two-body fragment, and enhanced planewave-corrected techniques. Combining planewave and two-body fragment calculations reduces the error in predicted ¹⁷ O C q values by 35% relative to traditional planewave calculations. These enhanced planewave-correction methods improve the accuracy of ¹⁷ O and ¹⁴ N EFG tensor components by 15% relative to planewave DFT but yield minimal improvement relative to a simple molecular correction. However, in structural environments involving either high symmetry or strong intermolecular interactions, enhanced planewave-corrected methods provide a distinct advantage.

Linear fractional charge behavior in density functional theory through dielectric tuning of conductor-like polarizable continuum model

A property of exact density functional theory is linear fractional charge behavior as electrons are added or removed from a molecule. Typical density functional approximations (DFAs) exhibit delocalization error, which overstabilizes this fractional charge. Conversely, solvent corrections have been shown to erroneously destabilize this fractional charge. This work will show that an implicit solvent correction with a tuned dielectric can be used as an ad hoc correction to offset the delocalizing character of DFAs and achieve linear fractional charge behavior. While desirable, in principle, we find that this linear charge behavior degrades the vertical ionization energies reported by DFAs. Our results reveal that the localizing character of the solvent correction and the Hartree-Fock (HF) exchange offset each other. This helps explain the decreased ratios of HF exchange to DFA exchange in long-range hybrid tuning studies that use a solvent correction.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles

In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H₂O)₁₂₈) and ≈6.8 s ((H₂O)₂₅₆) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.

Crystal Structure Prediction Methods for Organic Molecules: State of the Art

The prediction of the crystal structures that a given organic molecule is likely to form is an important theoretical problem of significant interest for the pharmaceutical and agrochemical industries, among others. As evidenced by a series of six blind tests organized over the past 2 decades, methodologies for crystal structure prediction (CSP) have witnessed substantial progress and have now reached a stage of development where they can begin to be applied to systems of practical significance. This article reviews the state of the art in general-purpose methodologies for CSP, placing them within a common framework that highlights both their similarities and their differences. The review discusses specific areas that constitute the main focus of current research efforts toward improving the reliability and widening applicability of these methodologies, and offers some perspectives for the evolution of this technology over the next decade.

Measuring Density-Driven Errors Using Kohn-Sham Inversion

Kohn-Sham (KS) inversion, that is, the finding of the exact KS potential for a given density, is difficult in localized basis sets. We study the precision and reliability of several inversion schemes, finding estimates of density-driven errors at a useful level of accuracy. In typical cases of substantial density-driven errors, Hartree-Fock density functional theory (HF-DFT) is almost as accurate as DFT evaluated on CCSD(T) densities. A simple approximation in practical HF-DFT also makes errors much smaller than the density-driven errors being calculated. Two paradigm examples, stretched NaCl and the HO·Cl^- radical, illustrate just how accurate HF-DFT is.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

Cost-Effective Quantum Mechanical Approach for Predicting Thermodynamic and Mechanical Stability of Pure-Silica Zeolites

Several computational techniques for solid-state applications have recently been proposed to enlarge the scope of computer simulations of large molecular systems. In this contribution, we focused on two of these, namely, HF-3c and PBEh-3c. They were recently proposed by the Grimme's group, as "low-cost" ab initio-based techniques for electronic structure calculation of large systems and were proved to be effective essentially for organic molecules. HF-3c is based on a Hartree-Fock Hamiltonian with a minimal Gaussian quality basis set, whereas PBEh-3c is a density functional theory (DFT) based method with a hybrid functional and a medium-quality basis set. Both HF-3c and PBEh-3c account for dispersion (London) interactions and are free from the basis set superposition error due to limited basis set size, through several pairwise semiempirical corrections. To the best of our knowledge, despite the promising results on the cost-accuracy side of molecular simulations of organic molecules, these methods have been used only in few cases for solid-state applications. In this contribution, we studied the performance of HF-3c and PBEh-3c for predicting the properties of inorganic crystals to enlarge the applicability of these cheap and fast methodologies. As a testing ground, we have chosen a well-known class of material, e.g., microporous all-silica zeolites. We benchmarked geometries, formation energies, vibrational features, and mechanical properties by comparing the results with literature data from both experiment and computer simulation. For structures, HF-3c is extremely accurate in predicting the zeolites cell volume, albeit we do not include any vibrational contribution, neither zero point nor thermal, on the computed volumes, which may introduce small variations in the predicted values. For the energetic, the relative stability of the zeolites using the DFT//HF-3c approach allows predictions within the experimental error for most of the cases taken into consideration when the experimental enthalpies were corrected back to electronic energies by using the HF-3c thermodynamic contributions computed in the harmonic approximation. This strategy is particularly convenient, as the slow step (geometry optimization) is carried out with the cheapest HF-3c method, whereas the fast step (single point energy evaluation) is carried out with costly DFT methods. In this sense, the use of the DFT//HF-3c approach results to be a promising one to predict the stability and structure of microporous materials. Finally, the HF-3c method predicts the mechanical properties of the zeolite set in reasonable agreement with respect to those computed with the state-of-the-art DFT simulations, indicating the HF-3c method as a possible technique for the mechanical stability screenings of microporous materials.

Accuracy and reproducibility in crystal structure prediction: the curious case of ROY

Because of excessive electron delocalization, the polymorphs of ROY constitute a surprisingly challenging system for crystal structure prediction.