1
|
Ganoe B, Shee J. On the notion of strong correlation in electronic structure theory. Faraday Discuss 2024. [PMID: 39072670 DOI: 10.1039/d4fd00066h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Strong correlation has been said to have many faces, and appears to have many synonyms of questionable suitability. In this work we aim not to define the term once and for all, but to highlight one possibility that is both rigorously defined and physically transparent, and remains so in reference to molecules and quantum lattice models. We survey both molecular examples - hydrogen systems (Hn, n = 2, 4, 6), Be2, H-He-H, and benzene - and the half-filled Hubbard model over a range of correlation regimes. Various quantities are examined including the extent of spin symmetry breaking in correlated single-reference wave functions, energetic ratios inspired by the Hubbard model and the Virial theorem, and metrics derived from the one- and two-electron reduced density matrices (RDMs). The trace and the square norm of the cumulant of the two-electron reduced density matrix capture what may well be defined as strong correlation. Accordingly, strong correlation is understood as a statistical dependence between two electrons, and is distinct from the concepts of "correlation energy" and more general than entanglement quantities that require a partitioning of a quantum system into distinguishable subspaces. This work enables us to build a bridge between a rigorous and quantifiable regime of strong electron correlation and more familiar chemical concepts such as anti-aromaticity in the context of Baird's rule.
Collapse
Affiliation(s)
- Brad Ganoe
- Department of Chemistry, Rice University, Houston, TX, 77005, USA.
| | - James Shee
- Department of Chemistry, Rice University, Houston, TX, 77005, USA.
| |
Collapse
|
2
|
Jana SK, Som NN, Jha PK. Size-Dependent Fullerenes for Enhanced Interaction of l-Leucine: A Combined DFT and MD Simulations Approach. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:13844-13859. [PMID: 38916256 DOI: 10.1021/acs.langmuir.4c00806] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
Fullerene-based biosensors have received great attention due to their unique electronic properties that allow them to transduce electrical signals by accepting electrons from amino acids. Babies with MSUD (maple syrup urine disease) are unable to break down amino acids such as l-leucine, and excess levels of the l-leucine are harmful. Therefore, sensing of l-leucine is foremost required. We aim to investigate the interaction tendencies of size-variable fullerenes (CX; X = 24, 36, 50, and 70) toward l-leucine (LEU) using density functional theory (DFT-D3) and classical molecular dynamics (MD) simulation. The C24 fullerene shows the highest affinity of the LEU biomolecule in the gas phase. Smaller fullerenes (C24 and C36) show stronger interactions with leucine due to their higher curvature in water environments. Moreover, recovery times in the ranges of 1010 and 104 s make it a viable candidate for the isolation application of LEU from the biological system. Further, the interaction between LEU and fullerenes is in line with the natural bond order (NBO) analysis, Mulliken charge analysis, quantum theory atom in molecule (QTAIM) analysis, and reduced density gradient (RDG) analysis. At 310 K, employing the explicit water model in classical MD simulations, fullerenes C24 and C36 demonstrate notably elevated binding free energies (-24.946 kJ/mol) in relation to LEU, showcasing their potential as sensors for l-leucine. Here, we demonstrate that the smaller fullerene exhibits a higher potential for l-leucine sensors than the larger fullerene.
Collapse
Affiliation(s)
- Sourav Kanti Jana
- Department of Physics, Faculty of Science, The Maharaja Sayajirao University of Baroda, Vadodara, Gujarat 39002, India
| | - Narayan N Som
- Institute of High-Pressure Physics, Polish Academy of Sciences, Sokolowska 29/37, 01-142 Warsaw, Poland
| | - Prafulla K Jha
- Department of Physics, Faculty of Science, The Maharaja Sayajirao University of Baroda, Vadodara, Gujarat 39002, India
| |
Collapse
|
3
|
Neugebauer H, Vuong HT, Weber JL, Friesner RA, Shee J, Hansen A. Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC. J Chem Theory Comput 2023; 19:6208-6225. [PMID: 37655473 DOI: 10.1021/acs.jctc.3c00617] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.
Collapse
Affiliation(s)
- Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Hung T Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
4
|
Marie A, Burton HGA. Excited States, Symmetry Breaking, and Unphysical Solutions in State-Specific CASSCF Theory. J Phys Chem A 2023; 127:4538-4552. [PMID: 37141564 DOI: 10.1021/acs.jpca.3c00603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
State-specific electronic structure theory provides a route toward balanced excited-state wave functions by exploiting higher-energy stationary points of the electronic energy. Multiconfigurational wave function approximations can describe both closed- and open-shell excited states and avoid the issues associated with state-averaged approaches. We investigate the existence of higher-energy solutions in complete active space self-consistent field (CASSCF) theory and characterize their topological properties. We demonstrate that state-specific approximations can provide accurate higher-energy excited states in H2 (6-31G) with more compact active spaces than would be required in a state-averaged formalism. We then elucidate the unphysical stationary points, demonstrating that they arise from redundant orbitals when the active space is too large or symmetry breaking when the active space is too small. Furthermore, we investigate the singlet-triplet crossing in CH2 (6-31G) and the avoided crossing in LiF (6-31G), revealing the severity of root flipping and demonstrating that state-specific solutions can behave quasi-diabatically or adiabatically. These results elucidate the complexity of the CASSCF energy landscape, highlighting the advantages and challenges of practical state-specific calculations.
Collapse
Affiliation(s)
- Antoine Marie
- Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, U.K
| | - Hugh G A Burton
- Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, U.K
| |
Collapse
|
5
|
Wong J, Ganoe B, Liu X, Neudecker T, Lee J, Liang J, Wang Z, Li J, Rettig A, Head-Gordon T, Head-Gordon M. An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3. J Chem Phys 2023; 158:164116. [PMID: 37114707 PMCID: PMC10148725 DOI: 10.1063/5.0145130] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 04/10/2023] [Indexed: 04/29/2023] Open
Abstract
We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.
Collapse
Affiliation(s)
- Jonathan Wong
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Brad Ganoe
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Xiao Liu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tim Neudecker
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joonho Lee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhe Wang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jie Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | |
Collapse
|
6
|
Lee J, Pham HQ, Reichman DR. Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking. J Chem Theory Comput 2022; 18:7024-7042. [PMID: 36255074 DOI: 10.1021/acs.jctc.2c00802] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) approach from a computational quantum chemistry perspective and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single-determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multireference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach.
Collapse
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Q Pham
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
7
|
Rana R, Thorpe JH, Stanton JF. An interesting case of symmetry breaking in the absence of symmetry: the bicarbonate radical (HCO 3). Mol Phys 2022. [DOI: 10.1080/00268976.2022.2144518] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Affiliation(s)
- Rohit Rana
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, FL, USA
| | - James H. Thorpe
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, FL, USA
| | - John F. Stanton
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, FL, USA
| |
Collapse
|
8
|
Rettig A, Shee J, Lee J, Head-Gordon M. Revisiting the Orbital Energy-Dependent Regularization of Orbital-Optimized Second-Order Møller-Plesset Theory. J Chem Theory Comput 2022; 18:5382-5392. [PMID: 36050889 DOI: 10.1021/acs.jctc.2c00641] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Optimizing orbitals in the presence of electron correlation, as in orbital-optimized second-order Møller-Plesset perturbation theory (OOMP2), can remove artifacts associated with mean-field orbitals such as spin contamination and artificial symmetry-breaking. However, OOMP2 is known to suffer from divergent correlation energies in regimes of small orbital energy gaps. To address this issue, several approaches to amplitude regularization have been explored, with those featuring energy-gap-dependent regularizers appearing to be most transferable and physically justifiable. For instance, κ-OOMP2 was shown to address the energy divergence issue in, for example, bond-breaking processes while offering a significant improvement in accuracy for the W4-11 thermochemistry data set, and a parameter of κ = 1.45 was recommended. A more recent investigation of regularized MP2 with Hartree-Fock orbitals revealed that stronger regularization (i.e., smaller values of κ) than what had previously been recommended for κ-OOMP2 may offer huge improvements in certain cases such as noncovalent interactions while retaining a high level of accuracy for main-group thermochemistry data sets. In this study, we investigate the transferability of those findings to κ-OOMP2 and assess the implications of stronger regularization on the ability of κ-OOMP2 to diagnose strong static correlation. We found similar results using κ-OOMP2 for several main-group thermochemistry, barrier height, and noncovalent interaction data sets including both closed shell and open shell species. However, stronger regularization yielded substantially higher accuracy for open-shell transition-metal (TM) thermochemistry and is necessary to provide qualitatively correct spin symmetry breaking behavior for several large and electrochemically relevant TM systems. We therefore find a single κ value insufficient to treat all systems using κ-OOMP2.
Collapse
Affiliation(s)
- Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Abstract
The computational modeling of fullerenes plays a fundamental role in designing low-dimension carbon nanostructures. Nevertheless, the relative energies of fullerenes larger than C20 and C24 have not been comprehensively examined by means of highly accurate ab initio methods, for example, the CCSD(T) method. Here we report such an investigation for a diverse set of 29 C40 isomers. We calculate the energies of the C40 fullerenes using the G4(MP2) composite ab initio method, which approximates the CCSD(T) energy in conjunction with a triple-ζ-quality basis set (CCSD(T)/TZ). The CCSD(T)/TZ isomerization energies span 43.1-763.3 kJ mol-1. We find a linear correlation (R2 = 0.96) between the CCSD(T)/TZ isomerization energies and the fullerene pentagon signatures (P1 index), which reflect the strain associated with fused pentagon-pentagon rings. Using the reference CCSD(T)/TZ isomerization energies, we examine the relationship between the percentage of exact Hartree-Fock (HF) exchange in hybrid density functional theory (DFT) methods and the pentagon-pentagon strain energies. We find that the performance of hybrid DFT methods deteriorates with the pentagon-pentagon strain energy. This deterioration in performance becomes more pronounced with the inclusion of high amounts of HF exchange. For example, for B3LYP (20% HF exchange), the root-mean-square deviation (RMSD) relative to G4(MP2) increases from 8.9 kJ mol-1 for the low-strain isomers (P1 = 11) to 18.0 kJ mol-1 for the high-strain isomers (P1 > 13). However, for BH&HLYP (50% HF exchange) the RMSD increases from 23.0 (P1 = 11) to 113.2 (P1 > 13) kJ mol-1. A similar trend is observed for the M06/M06-2X pair of functionals. Namely, for M06 (27% HF exchange) the RMSD increases from 0.8 (P1 = 11) to 21.0 (P1 > 13) kJ mol-1, whereas for M06-2X (54% HF exchange) the RMSD increases from 16.7 (P1 = 11) to 77.7 (P1 > 13) kJ mol-1. Overall, we find that the strain associated with pentagon adjacency is an inherently challenging problem for hybrid DFT methods involving high amounts of HF exchange and that there is an inverse relationship between the optimal percentage of HF exchange and the pentagon-pentagon strain energy. For example, for BLYP the optimal percentages of HF exchange are 13% (P1 = 11), 10% (P1 = 12), 7.5% (P1 = 13), and 6% (P1 > 13).
Collapse
Affiliation(s)
- Amir Karton
- School of Science and Technology, University of New England, Armidale, NSW 2351, Australia
| |
Collapse
|
10
|
Abstract
State-specific approximations can provide a more accurate representation of challenging electronic excitations by enabling relaxation of the electron density. While state-specific wave functions are known to be local minima or saddle points of the approximate energy, the global structure of the exact electronic energy remains largely unexplored. In this contribution, a geometric perspective on the exact electronic energy landscape is introduced. On the exact energy landscape, ground and excited states form stationary points constrained to the surface of a hypersphere, and the corresponding Hessian index increases at each excitation level. The connectivity between exact stationary points is investigated, and the square-magnitude of the exact energy gradient is shown to be directly proportional to the Hamiltonian variance. The minimal basis Hartree-Fock and excited-state mean-field representations of singlet H2 (STO-3G) are then used to explore how the exact energy landscape controls the existence and properties of state-specific approximations. In particular, approximate excited states correspond to constrained stationary points on the exact energy landscape, and their Hessian index also increases for higher energies. Finally, the properties of the exact energy are used to derive the structure of the variance optimization landscape and elucidate the challenges faced by variance optimization algorithms, including the presence of unphysical saddle points or maxima of the variance.
Collapse
Affiliation(s)
- Hugh G A Burton
- Physical and Theoretical Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
11
|
Schattenberg CJ, Lehmann M, Bühl M, Kaupp M. Systematic Evaluation of Modern Density Functional Methods for the Computation of NMR Shifts of 3d Transition-Metal Nuclei. J Chem Theory Comput 2021; 18:273-292. [PMID: 34968062 DOI: 10.1021/acs.jctc.1c00964] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A wide range of density functionals from all rungs of Jacob's ladder have been evaluated systematically for a set of experimental 3d transition-metal NMR shifts of 70 complexes encompassing 12 × 49Ti, 10 × 51V, 10 × 53Cr, 11 × 55Mn, 9 × 57Fe, 9 × 59Co, and 9 × 61Ni shift values, as well as a diverse range of electronic structure characteristics. The overall 39 functionals evaluated include one LDA, eight GGAs, seven meta-GGAs (including their current-density-functional─CDFT─versions), nine global hybrids, four range-separated hybrids, eight local hybrids, and two double hybrids, and we also include Hartree-Fock and MP2 calculations. While recent evaluations of the same functionals for a very large coupled-cluster-based benchmark of main-group shieldings and shifts achieved in some cases aggregate percentage mean absolute errors clearly below 2%, the best results for the present 3d-nuclei set are in the range between 4 and 5%. Strikingly, the overall best-performing functionals are the recently implemented CDFT versions of two meta-GGAs, namely cM06-L (4.0%) and cVSXC (4.3%), followed by cLH14t-calPBE (4.9%), B3LYP (5.0%), and cLH07t-SVWN (5.1%), i.e., the previously best-performing global hybrid and two local hybrids. A number of further functionals achieve aggregate deviations in the range 5-6%. Range-separated hybrids offer no particular advantage over global hybrids. Due to the overall poor performance of Hartree-Fock theory for all systems except the titanium complexes, MP2 and double-hybrid functionals are unsuitable for these 3d-nucleus shifts and provide large errors. Global hybrid functionals with larger EXX admixtures, such as BHLYP or M06-2X, also perform poorly, and some other highly parametrized global hybrids also are unsuitable. For many functionals depending on local kinetic energy τ, their CDFT variants perform much better than their "non-CDFT" versions. This holds notably also for the above-mentioned M06-L and VSXC, while the effect is small for τ-dependent local hybrids and can even be somewhat detrimental to the agreement with experiment for a few other cases. The separation between well-performing and more poorly performing functionals is mainly determined by their results for the most critical nuclei 55Mn, 57Fe, and 59Co. Here either moderate exact-exchange admixtures or CDFT versions of meta-GGAs are beneficial for the accuracy. The overall deviations of the better-performing global or local hybrids are then typically dominated by the 53Cr shifts, where triplet instabilities appear to disfavor exact-exchange admixture. Further detailed analyses help to pinpoint specific nuclei and specific types of complexes that are challenges for a given functional.
Collapse
Affiliation(s)
- Caspar Jonas Schattenberg
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Morten Lehmann
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Michael Bühl
- School of Chemistry, University of St. Andrews, Purdie Building, North Haugh, St Andrews KY16 9ST, Fife, U.K
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie Sekretariat C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623, Berlin, Germany
| |
Collapse
|
12
|
Shee J, Loipersberger M, Rettig A, Lee J, Head-Gordon M. Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost. J Phys Chem Lett 2021; 12:12084-12097. [PMID: 34910484 PMCID: PMC10037552 DOI: 10.1021/acs.jpclett.1c03468] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ2 regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.
Collapse
Affiliation(s)
- James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthias Loipersberger
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Stauch T, Ganoe B, Wong J, Lee J, Rettig A, Liang J, Li J, Epifanovsky E, Head-Gordon T, Head-Gordon M. Molecular magnetisabilities computed via finite fields: assessing alternatives to MP2 and revisiting magnetic exaltations in aromatic and antiaromatic species. Mol Phys 2021; 119. [DOI: 10.1080/00268976.2021.1990426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Tim Stauch
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Brad Ganoe
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jonathan Wong
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Joonho Lee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jie Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | | | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
14
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
15
|
Epifanovsky E, Gilbert ATB, Feng X, Lee J, Mao Y, Mardirossian N, Pokhilko P, White AF, Coons MP, Dempwolff AL, Gan Z, Hait D, Horn PR, Jacobson LD, Kaliman I, Kussmann J, Lange AW, Lao KU, Levine DS, Liu J, McKenzie SC, Morrison AF, Nanda KD, Plasser F, Rehn DR, Vidal ML, You ZQ, Zhu Y, Alam B, Albrecht BJ, Aldossary A, Alguire E, Andersen JH, Athavale V, Barton D, Begam K, Behn A, Bellonzi N, Bernard YA, Berquist EJ, Burton HGA, Carreras A, Carter-Fenk K, Chakraborty R, Chien AD, Closser KD, Cofer-Shabica V, Dasgupta S, de Wergifosse M, Deng J, Diedenhofen M, Do H, Ehlert S, Fang PT, Fatehi S, Feng Q, Friedhoff T, Gayvert J, Ge Q, Gidofalvi G, Goldey M, Gomes J, González-Espinoza CE, Gulania S, Gunina AO, Hanson-Heine MWD, Harbach PHP, Hauser A, Herbst MF, Hernández Vera M, Hodecker M, Holden ZC, Houck S, Huang X, Hui K, Huynh BC, Ivanov M, Jász Á, Ji H, Jiang H, Kaduk B, Kähler S, Khistyaev K, Kim J, Kis G, Klunzinger P, Koczor-Benda Z, Koh JH, Kosenkov D, Koulias L, Kowalczyk T, Krauter CM, Kue K, Kunitsa A, Kus T, Ladjánszki I, Landau A, Lawler KV, Lefrancois D, Lehtola S, Li RR, Li YP, Liang J, Liebenthal M, Lin HH, Lin YS, Liu F, Liu KY, Loipersberger M, Luenser A, Manjanath A, Manohar P, Mansoor E, Manzer SF, Mao SP, Marenich AV, Markovich T, Mason S, Maurer SA, McLaughlin PF, Menger MFSJ, Mewes JM, Mewes SA, Morgante P, Mullinax JW, Oosterbaan KJ, Paran G, Paul AC, Paul SK, Pavošević F, Pei Z, Prager S, Proynov EI, Rák Á, Ramos-Cordoba E, Rana B, Rask AE, Rettig A, Richard RM, Rob F, Rossomme E, Scheele T, Scheurer M, Schneider M, Sergueev N, Sharada SM, Skomorowski W, Small DW, Stein CJ, Su YC, Sundstrom EJ, Tao Z, Thirman J, Tornai GJ, Tsuchimochi T, Tubman NM, Veccham SP, Vydrov O, Wenzel J, Witte J, Yamada A, Yao K, Yeganeh S, Yost SR, Zech A, Zhang IY, Zhang X, Zhang Y, Zuev D, Aspuru-Guzik A, Bell AT, Besley NA, Bravaya KB, Brooks BR, Casanova D, Chai JD, Coriani S, Cramer CJ, Cserey G, DePrince AE, DiStasio RA, Dreuw A, Dunietz BD, Furlani TR, Goddard WA, Hammes-Schiffer S, Head-Gordon T, Hehre WJ, Hsu CP, Jagau TC, Jung Y, Klamt A, Kong J, Lambrecht DS, Liang W, Mayhall NJ, McCurdy CW, Neaton JB, Ochsenfeld C, Parkhill JA, Peverati R, Rassolov VA, Shao Y, Slipchenko LV, Stauch T, Steele RP, Subotnik JE, Thom AJW, Tkatchenko A, Truhlar DG, Van Voorhis T, Wesolowski TA, Whaley KB, Woodcock HL, Zimmerman PM, Faraji S, Gill PMW, Head-Gordon M, Herbert JM, Krylov AI. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J Chem Phys 2021; 155:084801. [PMID: 34470363 PMCID: PMC9984241 DOI: 10.1063/5.0055522] [Citation(s) in RCA: 466] [Impact Index Per Article: 155.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
Collapse
Affiliation(s)
- Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alec F. White
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zhengting Gan
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Leif D. Jacobson
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Jörg Kussmann
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Daniel S. Levine
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Simon C. McKenzie
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Marta L. Vidal
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | - Ying Zhu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Benjamin J. Albrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | - Ethan Alguire
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Josefine H. Andersen
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | - Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Dennis Barton
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Khadiza Begam
- Department of Physics, Kent State University, Kent, Ohio 44242, USA
| | - Andrew Behn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yves A. Bernard
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Abel Carreras
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Alan D. Chien
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | - Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Marc de Wergifosse
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jia Deng
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Hainam Do
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Beringstr. 4, 53115 Bonn, Germany
| | - Po-Tung Fang
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Qingguo Feng
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Triet Friedhoff
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - James Gayvert
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Qinghui Ge
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Matthew Goldey
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joe Gomes
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anastasia O. Gunina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Phillip H. P. Harbach
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Hauser
- Institute of Experimental Physics, Graz University of Technology, Graz, Austria
| | | | - Mario Hernández Vera
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Shannon Houck
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xunkun Huang
- Department of Chemistry, Xiamen University, Xiamen 361005, China
| | - Kerwin Hui
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Bang C. Huynh
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Maxim Ivanov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Ádám Jász
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Hyunjun Ji
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Sven Kähler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Kirill Khistyaev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jaehoon Kim
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Gergely Kis
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | | | - Zsuzsanna Koczor-Benda
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Joong Hoon Koh
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Dimitri Kosenkov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Laura Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | - Caroline M. Krauter
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Karl Kue
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Alexander Kunitsa
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Thomas Kus
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Arie Landau
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Keith V. Lawler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Run R. Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Yi-Pei Li
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marcus Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Fenglai Liu
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Arne Luenser
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Aaditya Manjanath
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Prashant Manohar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Erum Mansoor
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Sam F. Manzer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Thomas Markovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Stephen Mason
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Simon A. Maurer
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Peter F. McLaughlin
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | - Jan-Michael Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Pierpaolo Morgante
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J. Wayne Mullinax
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | | | | | - Alexander C. Paul
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Suranjan K. Paul
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Zheng Pei
- School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Stefan Prager
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Emil I. Proynov
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Ádám Rák
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Eloy Ramos-Cordoba
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alan E. Rask
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fazle Rob
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Elliot Rossomme
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tarek Scheele
- Institute for Physical and Theoretical Chemistry, University of Bremen, Bremen, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Matthias Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Nickolai Sergueev
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Shaama M. Sharada
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Wojciech Skomorowski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yu-Chuan Su
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Jonathan Thirman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gábor J. Tornai
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Takashi Tsuchimochi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Norm M. Tubman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Oleg Vydrov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jan Wenzel
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jon Witte
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Atsushi Yamada
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Sina Yeganeh
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Shane R. Yost
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - Igor Ying Zhang
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xing Zhang
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Yu Zhang
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Dmitry Zuev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Alexis T. Bell
- Department of Chemical Engineering, University of California, Berkeley, California 94720, USA
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biophysics, National Institute of Health, Bethesda, Maryland 20892, USA
| | - David Casanova
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Thomas R. Furlani
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - William A. Goddard
- Materials and Process Simulation Center, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | | | | - Yousung Jung
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
| | - Jing Kong
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Daniel S. Lambrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | - C. William McCurdy
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - John A. Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Roberto Peverati
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | | | | | | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, 9774AG Groningen, The Netherlands
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA,Author to whom correspondence should be addressed:
| |
Collapse
|
16
|
Cunha LA, Lee J, Hait D, McCurdy CW, Head-Gordon M. Exploring spin symmetry-breaking effects for static field ionization of atoms: Is there an analog to the Coulson-Fischer point in bond dissociation? J Chem Phys 2021; 155:014309. [PMID: 34241406 DOI: 10.1063/5.0054387] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree-Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single Ms (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree-Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis.
Collapse
Affiliation(s)
- Leonardo A Cunha
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - C William McCurdy
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
17
|
Shee J, Loipersberger M, Hait D, Lee J, Head-Gordon M. Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition. J Chem Phys 2021; 154:194109. [PMID: 34240907 DOI: 10.1063/5.0047386] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.
Collapse
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Matthias Loipersberger
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
18
|
Bertels LW, Lee J, Head-Gordon M. Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction. J Chem Theory Comput 2021; 17:742-755. [PMID: 33404238 DOI: 10.1021/acs.jctc.0c00746] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm-1, and 8.50 cm-1, and 8.75 cm-1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).
Collapse
Affiliation(s)
- Luke W Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
19
|
Tran HK, Ye HZ, Van Voorhis T. Bootstrap embedding with an unrestricted mean-field bath. J Chem Phys 2020; 153:214101. [PMID: 33291897 DOI: 10.1063/5.0029092] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
A suite of quantum embedding methods have recently been developed where the Schmidt decomposition is applied to the full system wavefunction to derive basis states that preserve the entanglement between the fragment and the bath. The quality of these methods can depend heavily on the quality of the initial full system wavefunction. Most of these methods, including bootstrap embedding (BE) [M. Welborn et al; J. Chem. Phys. 145, 074102 (2016)], start from a spin-restricted mean-field wavefunction [call this restricted BE (RBE)]. Given that spin-unrestricted wavefunctions can capture a significant amount of strong correlation at the mean-field level, we suspect that starting from a spin-unrestricted mean-field wavefunction will improve these embedding methods for strongly correlated systems. In this work, BE is generalized to an unrestricted Hartree-Fock bath [call this unrestricted BE (UBE)], and UBE is applied to model hydrogen ring systems. UBE's improved versatility over RBE is utilized to calculate high spin symmetry states that were previously unattainable with RBE. Ionization potentials, electron affinities, and spin-splittings are computed using UBE with accuracy on par with spin-unrestricted coupled cluster singles and doubles. Even for cases where RBE is viable, UBE converges more reliably. We discuss the limitations or weaknesses of each calculation and how improvements to RBE and density matrix embedding theory these past few years can also improve UBE.
Collapse
Affiliation(s)
- Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
20
|
Eken Y, Almeida NMS, Wang C, Wilson AK. SAMPL7: Host-guest binding prediction by molecular dynamics and quantum mechanics. J Comput Aided Mol Des 2020; 35:63-77. [PMID: 33150463 DOI: 10.1007/s10822-020-00357-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Accepted: 10/28/2020] [Indexed: 01/11/2023]
Abstract
Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenges provide routes to compare chemical quantities determined using computational chemistry approaches to experimental measurements that are shared after the competition. For this effort, several computational methods have been used to calculate the binding energies of Octa Acid (OA) and exo-Octa Acid (exoOA) host-guest systems for SAMPL7. The initial poses for molecular dynamics (MD) were generated by molecular docking. Binding free energy calculations were performed using molecular mechanics combined with Poisson-Boltzmann or generalized Born surface area solvation (MMPBSA/MMGBSA) approaches. The factors that affect the utility of the MMPBSA/MMGBSA approaches including solvation, partial charge, and solute entropy models were also analyzed. In addition to MD calculations, quantum mechanics (QM) calculations were performed using several different density functional theory (DFT) approaches. From SAMPL6 results, B3PW91-D3 was found to overestimate binding energies though it was effective for geometry optimizations, so it was considered for the DFT geometry optimizations in the current study, with single-point energy calculations carried out with B2PLYP-D3 with double-, triple-, and quadruple-ζ level basis sets. Accounting for dispersion effects, and solvation models was deemed essential for the predictions. MMGBSA and MMPBSA correlated better to experiment when used in conjunction with an empirical/linear correction.
Collapse
Affiliation(s)
- Yiğitcan Eken
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Nuno M S Almeida
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Cong Wang
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, 48864, USA.
- Department of Chemistry, University of North Texas, Denton, TX, 76201, USA.
| |
Collapse
|
21
|
Chan B. Fullerene Thermochemical Stability: Accurate Heats of Formation for Small Fullerenes, the Importance of Structural Deformation on Reactivity, and the Special Stability of C 60. J Phys Chem A 2020; 124:6688-6698. [PMID: 32786665 DOI: 10.1021/acs.jpca.0c04732] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
We have used quantum chemistry computations, in conjunction with isodesmic-type reactions, to obtain accurate heats of formation (HoFs) for the small fullerenes C20 (2358.2 ± 8.0 kJ mol-1), C24 (2566.2 ± 7.6), and the lowest-energy isomers of C32 (2461.1 ± 15.4), C42 (2629.0 ± 20.5), and C54 (2686.2 ± 25.3). As part of this endeavor, we have also obtained accurate HoFs for several medium-sized molecules, namely 216.6 ± 1.4 for fulvene, 375.5 ± 1.5 for pentalene, 670.8 ± 2.9 for acepentalene, and 262.7 ± 2.5 for acenaphthylene. We combine the energies of the small fullerenes and previously obtained energies for larger fullerenes (from C60 to C6000) into a full picture of fullerene thermochemical stability. In general, the per-carbon energies can be reasonably approximated by the "R+D" model that we have previously developed [Chan et al. J. Chem. Theory Comput. 2019, 15, 1255-1264], which takes into account Resonance and structural Deformation factors. In a case study on C54, we find that most of the high-deformation-energy atoms correspond to the sites of the C-Cl bond in the experimentally captured C54Cl8. In another case study, we find that C60 has the lowest value for the maximum local-deformation energy when compared with similar-sized fullerenes, which is consistent with its "special stability". These results are indicative of structural deformation playing an important role in the reactivity of fullerenes.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki-shi, Nagasaki 852-8521, Japan
| |
Collapse
|
22
|
Ye HZ, Tran HK, Van Voorhis T. Bootstrap Embedding For Large Molecular Systems. J Chem Theory Comput 2020; 16:5035-5046. [PMID: 32589842 DOI: 10.1021/acs.jctc.0c00438] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N3) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
23
|
Haggag OS, Malakar P, Pokhilko P, Stanton JF, Krylov AI, Ruhman S. The elusive dynamics of aqueous permanganate photochemistry. Phys Chem Chem Phys 2020; 22:10043-10055. [PMID: 32338267 DOI: 10.1039/c9cp07028a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Despite decades of investigation, mechanistic details of aqueous permanganate photo-decomposition remain unclear. Here we follow photoinduced dynamics of aqueous permanganate with femtosecond spectroscopy. Photoexcitation of KMnO4(aq) in the visible unleashes a sub-picosecond cascade of non-radiative transitions, leading to a distinct species which relaxes to S0 with a lifetime of 16 ps. Tuning excitation to the UV shows increasing formation of a metastable intermediate, which outlives our ∼1 ns window of detection. Guided by electronic structure calculations and observations from three pulse excitation experiments, we assign the 16 ps species as the lowest Jahn-Teller component of the 3T1 triplet state and suggest a plausible sequence of radiationless transitions, which rapidly populate it. In conjunction with photodecomposition quantum yields obtained from the literature, these results demonstrate that aqueous permanganate photo-decomposition proceeds through a long-lived intermediate which is formed in parallel to the triplet in less than one ps upon UV absorption. The possibility that this is the postulated highly oxidative peroxo species, a fraction of which leads to the stable (MnO2- + O2) fragments, is discussed. Finally, periodic modulations detected in the pump-probe signal are assigned to ground-state vibrational coherences excited by impulsive Raman. Their wavelength-dependent absolute phases outline the borders between adjacent electronic transitions in the linear spectrum of permanganate.
Collapse
Affiliation(s)
- Omer S Haggag
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel.
| | | | | | | | | | | |
Collapse
|
24
|
Rudshteyn B, Coskun D, Weber JL, Arthur EJ, Zhang S, Reichman DR, Friesner RA, Shee J. Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:3041-3054. [DOI: 10.1021/acs.jctc.0c00070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L. Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
25
|
Lee J, Malone FD, Morales MA. Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C36 Fullerene and an Iron Porphyrin Model Complex. J Chem Theory Comput 2020; 16:3019-3027. [DOI: 10.1021/acs.jctc.0c00055] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| |
Collapse
|
26
|
Yaghoobi F, Salehzadeh S, Maddah M. Quantum mechanics and molecular dynamics studies on the C…H interaction between small fullerenes (C36 and C24) and [M(H2O)6]2+ (M = Ca2+, Zn2+) cations. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.112339] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
27
|
Pokhilko P, Izmodenov D, Krylov AI. Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets. J Chem Phys 2020; 152:034105. [DOI: 10.1063/1.5138643] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Daniil Izmodenov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
28
|
Abstract
This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.
Collapse
Affiliation(s)
- David Casanova
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain. and IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Euskadi, Spain
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
29
|
Lee J, Small DW, Head-Gordon M. Excited states via coupled cluster theory without equation-of-motion methods: Seeking higher roots with application to doubly excited states and double core hole states. J Chem Phys 2019; 151:214103. [DOI: 10.1063/1.5128795] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
30
|
Lee J, Lin L, Head-Gordon M. Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory. J Chem Theory Comput 2019; 16:243-263. [DOI: 10.1021/acs.jctc.9b00820] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, California 94720, United States
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
31
|
Veccham SP, Lee J, Head-Gordon M. Making many-body interactions nearly pairwise additive: The polarized many-body expansion approach. J Chem Phys 2019; 151:194101. [DOI: 10.1063/1.5125802] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Srimukh Prasad Veccham
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
32
|
Loipersberger M, Lee J, Mao Y, Das AK, Ikeda K, Thirman J, Head-Gordon T, Head-Gordon M. Energy Decomposition Analysis for Interactions of Radicals: Theory and Implementation at the MP2 Level with Application to Hydration of Halogenated Benzene Cations and Complexes between CO2–· and Pyridine and Imidazole. J Phys Chem A 2019; 123:9621-9633. [DOI: 10.1021/acs.jpca.9b08586] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
33
|
Lee J, Malone FD, Morales MA. An auxiliary-Field quantum Monte Carlo perspective on the ground state of the dense uniform electron gas: An investigation with Hartree-Fock trial wavefunctions. J Chem Phys 2019. [DOI: 10.1063/1.5109572] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| |
Collapse
|
34
|
Bertels LW, Lee J, Head-Gordon M. Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics, and Intermolecular Interactions. J Phys Chem Lett 2019; 10:4170-4176. [PMID: 31259560 DOI: 10.1021/acs.jpclett.9b01641] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We develop and test methods that include second- and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, κ-OOMP2, denoted as MP3:κ-OOMP2. Testing MP3:κ-OOMP2 shows RMS errors that are 1.7-5 times smaller than those of MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled (c2 = 1.0; c3 = 0.8) denoted as MP2.8:κ-OOMP2. MP2.8:κ-OOMP2 yields significant additional improvement over MP3:κ-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no O(N6) iterations).
Collapse
Affiliation(s)
- Luke W Bertels
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Joonho Lee
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
35
|
Lee J, Head-Gordon M. Two single-reference approaches to singlet biradicaloid problems: Complex, restricted orbitals and approximate spin-projection combined with regularized orbital-optimized Møller-Plesset perturbation theory. J Chem Phys 2019; 150:244106. [DOI: 10.1063/1.5097613] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
36
|
Stein CJ, Reiher M. autoCAS: A Program for Fully Automated Multiconfigurational Calculations. J Comput Chem 2019; 40:2216-2226. [DOI: 10.1002/jcc.25869] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2019] [Revised: 05/11/2019] [Accepted: 05/15/2019] [Indexed: 12/11/2022]
Affiliation(s)
- Christopher J. Stein
- ETH Zürich, Laboratorium für Physikalische Chemie Vladimir‐Prelog‐Weg 2, 8093 Zürich Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie Vladimir‐Prelog‐Weg 2, 8093 Zürich Switzerland
| |
Collapse
|
37
|
Oosterbaan KJ, White AF, Head-Gordon M. Non-Orthogonal Configuration Interaction with Single Substitutions for Core-Excited States: An Extension to Doublet Radicals. J Chem Theory Comput 2019; 15:2966-2973. [DOI: 10.1021/acs.jctc.8b01259] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Katherine J. Oosterbaan
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Alec F. White
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|