1
|
Graf D, Thom AJW. Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost. J Chem Phys 2023; 159:174106. [PMID: 37921249 DOI: 10.1063/5.0168569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| |
Collapse
|
2
|
Hartman JD, Spock LE, Harper JK. Benchmark accuracy of predicted NMR observables for quadrupolar 14 N and 17 O nuclei in molecular crystals. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2023; 61:253-267. [PMID: 36567433 DOI: 10.1002/mrc.5328] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 12/21/2022] [Accepted: 12/22/2022] [Indexed: 06/17/2023]
Abstract
Nuclear quadrupole resonances for 14 N and 17 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 17 O quadrupolar coupling constants ( C q ). Here we extend this work by examining the accuracy of predicted 14 N and 17 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 17 O C q values by 35% relative to traditional planewave calculations. These enhanced planewave-correction methods improve the accuracy of 17 O and 14 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.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside, Riverside, California, USA
| | - Lilian E Spock
- Department of Chemistry, University of California, Riverside, Riverside, California, USA
| | - James K Harper
- Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah, USA
| |
Collapse
|
3
|
Hemmingsen LO, Hervir OAJ, Dale SG. Linear fractional charge behavior in density functional theory through dielectric tuning of conductor-like polarizable continuum model. J Chem Phys 2022; 156:014106. [PMID: 34998325 DOI: 10.1063/5.0067685] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Luke O Hemmingsen
- Research School of Chemistry, Australian National University, Acton 2601, Australia
| | - Oliver A J Hervir
- Research School of Chemistry, Australian National University, Acton 2601, Australia
| | - Stephen G Dale
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan 4111, Australia
| |
Collapse
|
4
|
Ko HY, Santra B, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles. J Chem Theory Comput 2021; 17:7789-7813. [PMID: 34775753 DOI: 10.1021/acs.jctc.0c01194] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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 ((H2O)128) and ≈6.8 s ((H2O)256) 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.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
5
|
Watts TA, Price LS, Price SL, Niederberger SM, Bertke JA, Swift JA. The Crystal Structure of 5‐Aminouracil and the Ambiguity of Alternative Polymorphs. Isr J Chem 2021. [DOI: 10.1002/ijch.202100062] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Taylor A. Watts
- Department of Chemistry Georgetown University 37th and O Sts NW Washington DC 20057–1227 USA
| | - Louise S. Price
- Department of Chemistry University College London 20 Gordon Street London WC1H 0AJ UK
| | - Sarah L. Price
- Department of Chemistry University College London 20 Gordon Street London WC1H 0AJ UK
| | - Sara M. Niederberger
- Department of Chemistry Georgetown University 37th and O Sts NW Washington DC 20057–1227 USA
| | - Jeffery A. Bertke
- Department of Chemistry Georgetown University 37th and O Sts NW Washington DC 20057–1227 USA
| | - Jennifer A. Swift
- Department of Chemistry Georgetown University 37th and O Sts NW Washington DC 20057–1227 USA
| |
Collapse
|
6
|
Bowskill DH, Sugden IJ, Konstantinopoulos S, Adjiman CS, Pantelides CC. Crystal Structure Prediction Methods for Organic Molecules: State of the Art. Annu Rev Chem Biomol Eng 2021; 12:593-623. [PMID: 33770462 DOI: 10.1146/annurev-chembioeng-060718-030256] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
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.
Collapse
Affiliation(s)
- David H Bowskill
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Isaac J Sugden
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Stefanos Konstantinopoulos
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Claire S Adjiman
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Constantinos C Pantelides
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| |
Collapse
|
7
|
Nam S, Song S, Sim E, Burke K. Measuring Density-Driven Errors Using Kohn-Sham Inversion. J Chem Theory Comput 2020; 16:5014-5023. [PMID: 32667787 DOI: 10.1021/acs.jctc.0c00391] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Seungsoo Nam
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
8
|
Ko HY, Jia J, Santra B, Wu X, Car R, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance. J Chem Theory Comput 2020; 16:3757-3785. [PMID: 32045232 DOI: 10.1021/acs.jctc.9b01167] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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., (H2O)256) 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.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Princeton University, Princeton, New Jersey 08544, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
9
|
Cutini M, Civalleri B, Ugliengo P. Cost-Effective Quantum Mechanical Approach for Predicting Thermodynamic and Mechanical Stability of Pure-Silica Zeolites. ACS OMEGA 2019; 4:1838-1846. [PMID: 31459438 PMCID: PMC6648702 DOI: 10.1021/acsomega.8b03135] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/09/2018] [Accepted: 01/10/2019] [Indexed: 06/10/2023]
Abstract
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.
Collapse
|
10
|
Nyman J, Yu L, Reutzel-Edens SM. Accuracy and reproducibility in crystal structure prediction: the curious case of ROY. CrystEngComm 2019. [DOI: 10.1039/c8ce01902a] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Because of excessive electron delocalization, the polymorphs of ROY constitute a surprisingly challenging system for crystal structure prediction.
Collapse
Affiliation(s)
- Jonas Nyman
- School of Pharmacy
- University of Wisconsin – Madison
- Madison
- USA
- Small Molecule Design & Development
| | - Lian Yu
- School of Pharmacy
- University of Wisconsin – Madison
- Madison
- USA
| | | |
Collapse
|