1
|
Wei Y, Debnath S, Weber JL, Mahajan A, Reichman DR, Friesner RA. Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry. J Phys Chem A 2024. [PMID: 38970826 DOI: 10.1021/acs.jpca.4c02853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/08/2024]
Abstract
This study evaluates the precision of widely recognized quantum chemical methodologies, CCSD(T), DLPNO-CCSD(T), and localized ph-AFQMC, for determining the thermochemistry of main group elements. DLPNO-CCSD(T) and localized ph-AFQMC, which offer greater scalability compared to canonical CCSD(T), have emerged over the past decade as pivotal in producing precise benchmark chemical data. Our investigation includes closed-shell, neutral molecules, focusing on their heat of formation and atomization energy sourced from four specific small molecule data sets. First, we selected molecules from the G2 and G3 data sets, noted for their reliable experimental heat of formation data. Additionally, we incorporate molecules from the W4-11 and W4-17 sets, which provide high-level theoretical reference values for atomization energy at 0 K. Our findings reveal that both DLPNO-CCSD(T) and ph-AFQMC methods are capable of achieving a root-mean-square deviation of less than 1 kcal/mol across the combined data set, aligning with the threshold for chemical accuracy. Moreover, we make efforts to confine the maximum deviations within 2 kcal/mol, a degree of precision that significantly broadens the applicability of these methods in fields such as biology and materials science.
Collapse
Affiliation(s)
- Yujing Wei
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Sibali Debnath
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Ankit Mahajan
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
2
|
Otlyotov AA, Moshchenkov AD, Minenkov Y. Ni, Cu, Zn, Pd, Ag and Cd Tetraphenylporphyrin Ab Initio Thermochemistry: Enthalpy of Formation of ZnTPP Revisited. Inorg Chem 2024; 63:10230-10239. [PMID: 38780084 DOI: 10.1021/acs.inorgchem.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
Groups 10-12 metalloporphyrins have been recognized for their numerous properties essential for the development of new sensing materials. In this work, accurate gas-phase enthalpies of formation, ΔfHm0(g,298.15), are predicted for the series of Ni, Cu, Zn, Pd, Ag, and Cd tetraphenylporphyrins (MTPPs) on the basis of the reaction-based Feller-Peterson-Dixon approach and high-level ab initio DLPNO-CCSD(T) calculations. Our recently developed automatic generator of the balanced chemical reactions was employed to reduce the bias of the theoretical ΔfHm0(g,298.15) toward a particular reaction. Theoretical ΔfHm0(g,298.15) for ZnTPP (227.0 ± 3.4 kcal mol-1) does not support the previously reported experimental value of 132 ± 2 kcal mol-1. The origin of the discrepancy probably lies in the experimental solid-state ΔfHm0(ZnTPP, cr,298.15) as it stems from our theoretical evaluations of the ΔfHm0(cr,298.15) values for the entire set of transition metal TPP complexes. The large discrepancy between experiment and theory also holds when different DFT functionals (ωB97M-V, PBE0-D4, and B3LYP-D4) paired with quadruple-ζ quality basis sets are used for the theoretical calculations. Experimental revisiting of the solid-state enthalpy of formation of ZnTPP and analogue measurements for other transition metal TPPs are needed to resolve the observed discrepancy. Based on the predicted enthalpies of formation of MTPPs, the relative energies of the metal-ligand bonding are evaluated and the trends are compared to those for the complexes of the unsubstituted porphyrin with the same set of metals derived in [Can. J. Chem., 2009, 87, 1063]. According to both studies, Pd complexes exhibit the strongest bonding, while the Cd species are the least stable metallocomplexes within the considered series.
Collapse
Affiliation(s)
- Arseniy A Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| | - Andrey D Moshchenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation
| |
Collapse
|
3
|
Gómez S, Flórez E, Acelas N, Cappelli C, Hadad C, Restrepo A. Encapsulation of charged halogens by the 5 12 water cage. Phys Chem Chem Phys 2024; 26:15426-15436. [PMID: 38747303 DOI: 10.1039/d4cp01340a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Elizabeth Flórez
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Nancy Acelas
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Cacier Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| |
Collapse
|
4
|
Beletsan OB, Gordiy I, Lunkov SS, Kalinin MA, Alkhimova LE, Nosach EA, Ilin EA, Bespalov AV, Dallakyan OL, Chamkin AA, Prolomov IV, Zaripov RA, Pershin AA, Protsenko BO, Rusalev YV, Oganov RA, Kovaleva DK, Mironov VA, Dotsenko VV, Genaev AM, Sharapa DI, Tikhonov DS. From a humorous post to a detailed quantum-chemical study: isocyanate synthesis revisited. Phys Chem Chem Phys 2024; 26:13850-13861. [PMID: 38656824 DOI: 10.1039/d3cp04654k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Isocyanates play an essential role in modern manufacturing processes, especially in polyurethane production. There are numerous synthesis strategies for isocyanates both under industrial and laboratory conditions, which do not prevent searching for alternative highly efficient synthetic protocols. Here, we report a detailed theoretical investigation of the mechanism of sulfur dioxide-catalyzed rearrangement of phenylnitrile oxide into phenyl isocyanate, which was first reported in 1977. The DLPNO-CCSD(T) method and up-to-date DFT protocols were used to perform a highly accurate quantum-chemical study of the rearrangement mechanism. An overview of various organic and inorganic catalysts has revealed other potential catalysts, such as sulfur trioxide and selenium dioxide. Furthermore, the present study elucidated how substituents in phenylnitrile oxide influence reaction kinetics. This study was performed by a self-organized collaboration of scientists initiated by a humorous post on the VK social network.
Collapse
Affiliation(s)
- Oleg B Beletsan
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Igor Gordiy
- ChemU Corporation Ltd, 17 17 Gr. Xenopoulou St., 3106 Limassol, Cyprus
| | - Sergey S Lunkov
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
- N. D. Zelinsky Institute of Organic Chemistry, Leninsky prosp. 47, 119991 Moscow, Russia
| | - Mikhail A Kalinin
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Larisa E Alkhimova
- Center for Nature-Inspired Engineering, University of Tyumen, 625003 Tyumen, Russia
- School of Natural Sciences, University of Tyumen, 625003 Tyumen, Russia
| | - Egor A Nosach
- Department of Fundamental Physical and Chemical Engineering, Lomonosov Moscow State University, 119234 Moscow, Russia
| | - Egor A Ilin
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
- N. D. Zelinsky Institute of Organic Chemistry, Leninsky prosp. 47, 119991 Moscow, Russia
| | - Alexandr V Bespalov
- Department of Chemistry and High Technologies, Kuban State University, 149 Stavropolskaya St., 350040 Krasnodar, Russia
| | - Olgert L Dallakyan
- Computational Material Science Laboratory, Department of Physics, Yerevan State University, 0025 Yerevan, Armenia
| | - Aleksandr A Chamkin
- A.N. Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, 119334 Moscow, Russia
| | - Ilya V Prolomov
- D. Mendeleev University of Chemical Technology of Russia, 125047 Moscow, Russia
- N. D. Zelinsky Institute of Organic Chemistry, Leninsky prosp. 47, 119991 Moscow, Russia
| | - Radion A Zaripov
- Skolkovo Institute of Science and Technology, 121205 Moscow, Russia
| | - Andrey A Pershin
- Samara Branch of Lebedev Physical Institute, 443011 Samara, Russia
- Department of Physics, Samara University, 443086 Samara, Russia
| | - Bogdan O Protsenko
- The Smart Materials Research Institute, Southern Federal University, 344090 Rostov-on-Don, Russia
| | - Yury V Rusalev
- The Smart Materials Research Institute, Southern Federal University, 344090 Rostov-on-Don, Russia
| | - Ruslan A Oganov
- Department of Biochemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Diana K Kovaleva
- Department of Biochemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Vladimir A Mironov
- A. M. Butlerov Chemistry Institute, Kazan Federal University, 420008 Kazan, Russia
| | - Victor V Dotsenko
- Department of Chemistry and High Technologies, Kuban State University, 149 Stavropolskaya St., 350040 Krasnodar, Russia
- Faculty of Chemistry and Pharmacy, North-Caucasus Federal University, 355017 Stavropol, Russia
| | - Alexandr M Genaev
- N.N. Vorozhtsov Institute of Organic Chemistry, 630090 Novosibirsk, Russia.
| | | | | |
Collapse
|
5
|
Bistoni G, Altun A, Wang Z, Neese F. Local Energy Decomposition Analysis of London Dispersion Effects: From Simple Model Dimers to Complex Biomolecular Assemblies. Acc Chem Res 2024; 57:1411-1420. [PMID: 38602396 PMCID: PMC11080063 DOI: 10.1021/acs.accounts.4c00085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Revised: 04/04/2024] [Accepted: 04/04/2024] [Indexed: 04/12/2024]
Abstract
ConspectusLondon dispersion (LD) forces are ubiquitous in chemistry, playing a pivotal role in a wide range of chemical processes. For example, they influence the structure of molecular crystals, the selectivity of organocatalytic transformations, and the formation of biomolecular assemblies. Harnessing these forces for chemical applications requires consistent quantification of the LD energy across a broad and diverse spectrum of chemical scenarios. Despite the great progress made in recent years in the development of experimental strategies for LD quantification, quantum chemical methods remain one of the most useful tools in the hand of chemists for the study of these weak interactions. Unfortunately, the accurate quantification of LD effects in complex systems poses many challenges for electronic structure theories. One of the problems stems from the fact that LD forces originate from long-range electronic dynamic correlation, and hence, their rigorous description requires the use of complex, highly correlated wave function-based methods. These methods typically feature a steep scaling with the system size, limiting their applicability to small model systems. Another core challenge lies in disentangling short-range from long-range dynamic correlation, which from a rigorous quantum mechanical perspective is not possible.In this Account, we describe our research endeavors in the development of broadly applicable computational methods for LD quantification in molecular chemistry as well as challenging applications of these schemes in various domains of chemical research. Our strategy lies in the use of local correlation theories to reduce the computational cost associated with complex electronic structure methods while providing at the same time a simple means of decomposition of dynamic correlation into its long-range and short-range components. In particular, the local energy decomposition (LED) scheme at the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) level has emerged as a powerful tool in our research, offering a clear-cut quantitative definition of the LD energy that remains valid across a plethora of different chemical scenarios. Typical applications of this scheme are examined, encompassing protein-ligand interactions and reactivity studies involving many fragments and complex electronic structures. In addition, our research also involves the development of novel cost-effective methodologies, which exploit the LED definition of the LD energy, for LD energy quantification that are, in principle, applicable to systems with thousands of atoms. The Hartree-Fock plus London Dispersion (HFLD) scheme, correcting the HF interaction energy using an approximate CCSD(T)-based LD energy, is a useful, parameter-free electronic structure method for the study of LD effects in systems with hundreds of molecular fragments. However, the usefulness of the LED scheme reaches beyond providing an interpretation of the calculated DLPNO-CCSD(T) or DLPNO-MP2 interaction energies. For example, the dispersion energies obtained from the LED can be fruitfully used in order to parametrize semiempirical dispersion models. We will demonstrate this in the context of an emerging semiempirical method, namely, the Natural Orbital Tied Constructed Hamiltonian (NOTCH) method. NOTCH incorporates LED-derived LD energies and shows promising accuracy at a minimum amount of empiricism. Thus, it holds substantial promise for large and complex systems.
Collapse
Affiliation(s)
- Giovanni Bistoni
- Department
of Chemistry, Biology and Biotechnology, University of Perugia Via Elce di Sotto, 8, 06123 Perugia, Italy
| | - Ahmet Altun
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Zikuan Wang
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
6
|
Xi ZK, Ding YH, Tian X. Building a New Platform for Significantly Improving Performance of Hartree-Fock and CCSD(T) Correlation Energy Based on Two-Point Complete Basis Set Extrapolation Schemes. J Phys Chem A 2024. [PMID: 38686765 DOI: 10.1021/acs.jpca.4c01712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2024]
Abstract
The leading cause of high expense in gold standard coupled cluster theory is that calculations of electronic energies converge exceedingly slowly with an increased basis set size. Extrapolation principally allows for achieving higher-quality outcomes at reduced costs. Numerous extrapolation formulas have been developed, with attempts to predict energies up to the complete basis set limit. Unfortunately, since the intricate shape of the function hinges on the molecular properties with the highest angular momentum of the basis set, the accuracy of the extrapolated energies highly depends on the fitted empirical parameters, which rely on the quality of the data sets for fitting. In this work, to overcome the extrapolation deficiency caused by the very limited data sets and smaller basis sets in the early stages, we constructed a new benchmark platform that includes a broader data set of 183 species (containing open-shell, closed-shell, ionic, and neutral species) and a larger basis set up to aug-cc-pV6Z. The newly optimized parameters can significantly improve the energy-predictive abilities of ten published formulas. Notably, all ten formulas perform quite similarly under the new platform with the reoptimized parameters. Finally, we built an online calculator for researchers to use for these extrapolation schemes. Our work would reignite the interest and applications of the underestimated formulas.
Collapse
Affiliation(s)
- Zhao-Kai Xi
- †Key Laboratory of Carbon Materials of Zhejiang Province, Wenzhou Key Lab of Advanced Energy Storage and Conversion, Zhejiang Province Key Lab of Leather Engineering, College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou 325035, P. R. China
| | - Yi-Hong Ding
- Key Laboratory of Carbon Materials of Zhejiang Province, Wenzhou Key Lab of Advanced Energy Storage and Conversion, Zhejiang Province Key Lab of Leather Engineering, College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou 325035, P. R. China
- Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Xiao Tian
- School of Mathematics and Science, Hebei GEO University, Shijiazhuang 050031, P. R. China
| |
Collapse
|
7
|
Drabik G, Radoń M. Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes. J Chem Theory Comput 2024; 20:3199-3217. [PMID: 38574194 PMCID: PMC11044276 DOI: 10.1021/acs.jctc.4c00092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 03/19/2024] [Accepted: 03/20/2024] [Indexed: 04/06/2024]
Abstract
Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.
Collapse
Affiliation(s)
- Gabriela Drabik
- Jagiellonian
University, Doctoral School
of Exact and Natural Sciences, Łojasiewicza 11, 30-348 Kraków, Poland
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| | - Mariusz Radoń
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| |
Collapse
|
8
|
Sorathia K, Frantzov D, Tew DP. Improved CPS and CBS Extrapolation of PNO-CCSD(T) Energies: The MOBH35 and ISOL24 Data Sets. J Chem Theory Comput 2024; 20:2740-2750. [PMID: 38513261 PMCID: PMC11008106 DOI: 10.1021/acs.jctc.3c00974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 03/06/2024] [Accepted: 03/08/2024] [Indexed: 03/23/2024]
Abstract
Computation of heats of reaction of large molecules is now feasible using the domain-based pair natural orbital (PNO)-coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] theory. However, to obtain agreement within 1 kcal/mol of experiment, it is necessary to eliminate basis set incompleteness error, which comprises both the AO basis set error and the PNO truncation error. Our investigation into the convergence to the canonical limit of PNO-CCSD(T) energies with the PNO truncation threshold T shows that errors follow the model E ( T ) = E + A T 1 / 2 . Therefore, PNO truncation errors can be eliminated using a simple two-point CPS extrapolation to the canonical limit so that subsequent CBS extrapolation is not limited by the residual PNO truncation error. Using the ISOL24 and MOBH35 data sets, we find that PNO truncation errors are larger for molecules with significant static correlation and that it is necessary to use very tight thresholds of T = 10 - 8 to ensure that errors do not exceed 1 kcal/mol. We present a lower-cost extrapolation scheme that uses information from small basis sets to estimate the PNO truncation errors for larger basis sets. In this way, the canonical limit of CCSD(T) calculations on sizable molecules with large basis sets can be reliably estimated in a practical way. Using this approach, we report near complete basis set (CBS)-CCSD(T) reaction energies for the full ISOL24 and MOBH35 data sets.
Collapse
Affiliation(s)
- Kesha Sorathia
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| | - Damyan Frantzov
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| | - David P. Tew
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| |
Collapse
|
9
|
Wappett DA, Goerigk L. Exploring CPS-Extrapolated DLPNO-CCSD(T 1) Reference Values for Benchmarking DFT Methods on Enzymatically Catalyzed Reactions. J Phys Chem A 2024; 128:62-72. [PMID: 38124376 DOI: 10.1021/acs.jpca.3c05086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Domain-based local pair natural orbital coupled-cluster singles doubles with perturbative triples [DLPNO-CCSD(T)] is regularly used to calculate reliable benchmark reference values at a computational cost significantly lower than that of canonical CCSD(T). Recent work has shown that even greater accuracy can be obtained at only a small additional cost through extrapolation to the complete PNO space (CPS) limit. Herein, we test two levels of CPS extrapolation, CPS(5,6), which approximates the accuracy of standard TightPNO, and CPS(6,7), which surpasses it, as benchmark values to test density functional approximations (DFAs) on a small set of organic and transition-metal-dependent enzyme active site models. Between the different reference levels of theory, there are changes in the magnitudes of the absolute deviations for all functionals, but these are small and there is minimal impact on the relative rankings of the tested DFAs. The differences are more significant for the metalloenzymes than the organic enzymes, so we repeat the tests on our entire ENZYMES22 set of organic enzyme active site models [Wappett, D. A.; Goerigk, L. J. Phys. Chem. A 2019, 123, 7057-7074] to confirm that using the CPS extrapolations for the reference values has negligible impact on the benchmarking outcomes. This means that we can particularly recommend CPS(5,6) as an alternative to standard TightPNO settings for calculating reference values, increasing the applicability of DLPNO-CCSD(T) in benchmarking reaction energies and barrier heights of larger models of organic enzymes. DLPNO-CCSD(T1)/CPS(6,7) energies for ENZYMES22 are finally presented as updated reference values for the set, reflecting the recent improvements in the method.
Collapse
Affiliation(s)
- Dominique A Wappett
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
10
|
Stylianakis I, Zervos N, Lii JH, Pantazis DA, Kolocouris A. Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory. J Comput Aided Mol Des 2023; 37:607-656. [PMID: 37597063 PMCID: PMC10618395 DOI: 10.1007/s10822-023-00513-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Accepted: 06/03/2023] [Indexed: 08/21/2023]
Abstract
We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.
Collapse
Affiliation(s)
- Ioannis Stylianakis
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Nikolaos Zervos
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Jenn-Huei Lii
- Department of Chemistry, National Changhua University of Education, Changhua City, Taiwan
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| | - Antonios Kolocouris
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece.
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece.
| |
Collapse
|
11
|
Wang Y, Guo Y, Neese F, Valeev EF, Li W, Li S. Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules. J Chem Theory Comput 2023; 19:8076-8089. [PMID: 37920973 DOI: 10.1021/acs.jctc.3c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
12
|
Wappett D, Goerigk L. Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set. J Chem Theory Comput 2023; 19:8365-8383. [PMID: 37943578 PMCID: PMC10688432 DOI: 10.1021/acs.jctc.3c00558] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 10/13/2023] [Accepted: 10/16/2023] [Indexed: 11/10/2023]
Abstract
We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.
Collapse
Affiliation(s)
- Dominique
A. Wappett
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| |
Collapse
|
13
|
Neugebauer H, Pinski P, Grimme S, Neese F, Bursch M. Assessment of DLPNO-MP2 Approximations in Double-Hybrid DFT. J Chem Theory Comput 2023; 19:7695-7703. [PMID: 37862406 PMCID: PMC10653103 DOI: 10.1021/acs.jctc.3c00896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Indexed: 10/22/2023]
Abstract
The unfavorable scaling (N5) of the conventional second-order Møller-Plesset theory (MP2) typically prevents the application of double-hybrid (DH) density functionals to large systems with more than 100 atoms. A prominent approach to reduce the computational demand of electron correlation methods is the domain-based local pair natural orbital (DLPNO) approximation that is successfully used in the framework of DLPNO-CCSD(T). Its extension to MP2 [Pinski P.; Riplinger, C.; Valeev, E. F.; Neese, F. J. Chem. Phys. 2015, 143, 034108.] paved the way for DLPNO-based DH (DLPNO-DH) methods. In this work, we assess the accuracy of the DLPNO-DH approximation compared to conventional DHs on a large number of 7925 data points for thermochemistry and 239 data points for structural features, including main-group and transition-metal systems. It is shown that DLPNO-DH-DFT can be applied successfully to perform energy calculations and geometry optimizations for large molecules at a drastically reduced computational cost. Furthermore, PNO space extrapolation is shown to be applicable, similar to its DLPNO-CCSD(T) counterpart, to reduce the remaining error.
Collapse
Affiliation(s)
- Hagen Neugebauer
- Mulliken
Center for Theoretical Chemistry, Clausius Institute for Physical
and Theoretical Chemistry, University of
Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Peter Pinski
- HQS
Quantum Simulations GmbH, Rintheimer Straße 23, D-76131 Karlsruhe, Germany
| | - Stefan Grimme
- Mulliken
Center for Theoretical Chemistry, Clausius Institute for Physical
and Theoretical Chemistry, University of
Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Markus Bursch
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
14
|
Debnath S, Neufeld VA, Jacobson LD, Rudshteyn B, Weber JL, Berkelbach TC, Friesner RA. Accurate Quantum Chemical Reaction Energies for Lithium-Mediated Electrolyte Decomposition and Evaluation of Density Functional Approximations. J Phys Chem A 2023; 127:9178-9184. [PMID: 37878768 PMCID: PMC10795021 DOI: 10.1021/acs.jpca.3c04369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
An important concern related to the performance of Li-ion batteries is the formation of a solid electrolyte interphase on the surface of the anode. This film is formed from the decomposition of electrolytes and can have important effects on the stability and performance. Here, we evaluate the decomposition pathway of ethylene carbonate and related organic electrolyte molecules using a series of density functional approximations and correlated wave function (WF) methods, including the coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] and auxiliary-field quantum Monte Carlo (AFQMC). We find that the transition state barrier associated with ring opening varies widely across different functionals, ranging from 3.01 to 17.15 kcal/mol, which can be compared to the value of 12.84 kcal/mol predicted by CCSD(T). This large variation underscores the importance of benchmarking against accurate WF methods. A performance comparison of all of the density functionals used in this study reveals that the M06-2X-D3 (a meta-hybrid GGA), CAM-B3LYP-D3 (a range-separated hybrid), and B2GP-PLYP-D3 (a double hybrid) perform the best, with average errors of about 1.50-1.60 kcal/mol compared to CCSD(T). We also compared the performance of the WF methods that are more scalable than CCSD(T), finding that DLPNO-CCSD(T) and phaseless AFQMC with a DFT trial wave function exhibit average errors of 1.38 and 1.74 kcal/mol, respectively.
Collapse
Affiliation(s)
- Sibali Debnath
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Verena A Neufeld
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Benjamin Rudshteyn
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| | - John L Weber
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| |
Collapse
|
15
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
16
|
Pogrebetsky J, Siklitskaya A, Kubas A. MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations. J Chem Theory Comput 2023. [PMID: 37338422 DOI: 10.1021/acs.jctc.3c00444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.
Collapse
Affiliation(s)
- James Pogrebetsky
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Alexandra Siklitskaya
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| |
Collapse
|
17
|
Shee J, Weber JL, Reichman DR, Friesner RA, Zhang S. On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond. J Chem Phys 2023; 158:140901. [PMID: 37061483 PMCID: PMC10089686 DOI: 10.1063/5.0134009] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.
Collapse
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
| |
Collapse
|
18
|
Altun A, Riplinger C, Neese F, Bistoni G. Exploring the Accuracy Limits of PNO-Based Local Coupled-Cluster Calculations for Transition-Metal Complexes. J Chem Theory Comput 2023; 19:2039-2047. [PMID: 36917767 PMCID: PMC10100528 DOI: 10.1021/acs.jctc.3c00087] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/15/2023]
Abstract
While the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples (DLPNO-CCSD(T)) has proven instrumental for computing energies and properties of large and complex systems accurately, calculations on first-row transition metals with a complex electronic structure remain challenging. In this work, we identify and address the two main error sources that influence the DLPNO-CCSD(T) accuracy in this context, namely, (i) correlation effects from the 3s and 3p semicore orbitals and (ii) dynamic correlation-induced orbital relaxation (DCIOR) effects that are not described by the local MP2 guess. We present a computational strategy that allows us to completely eliminate the DLPNO error associated with semicore correlation effects, while increasing, at the same time, the efficiency of the method. As regards the DCIOR effects, we introduce a diagnostic for estimating the deviation between DLPNO-CCSD(T) and canonical CCSD(T) for systems with significant orbital relaxation.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Department of Chemistry, Biology, and Biotechnology, University of Perugia, 06123 Perugia, Italy
| |
Collapse
|
19
|
Aðalsteinsson HM, Bjornsson R. Ionization energies of metallocenes: a coupled cluster study of cobaltocene. Phys Chem Chem Phys 2023; 25:4570-4587. [PMID: 36723003 DOI: 10.1039/d2cp04715b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T1) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T1) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.
Collapse
Affiliation(s)
| | - Ragnar Bjornsson
- Science Institute, University of Iceland, 107 Reykjavik, Iceland.,Univ Grenoble Alpes, CNRS, CEA, IRIG, Laboratoire de Chimie et Biologie des Métaux, 17 Rue des Martyrs, F-38054 Grenoble Cedex, France.
| |
Collapse
|
20
|
Schümann JM, Ochmann L, Becker J, Altun A, Harden I, Bistoni G, Schreiner PR. Exploring the Limits of Intramolecular London Dispersion Stabilization with Bulky Dispersion Energy Donors in Alkane Solution. J Am Chem Soc 2023; 145:2093-2097. [PMID: 36688409 DOI: 10.1021/jacs.2c13301] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
We present an experimental study of a cyclooctatetraene-based molecular balance disubstituted with increasingly bulky tert-butyl (tBu), adamantyl (Ad), and diamantyl (Dia) substituents in the 1,4-/1,6-positions for which we determined the valence-bond shift equilibrium in n-hexane (hex), n-octane (oct), and n-dodecane (dod). Computations including implicit and explicit solvation support our temperature-dependent NMR equilibrium measurements indicating that the more sterically crowded 1,6-isomer is always favored, irrespective of solvent, and that the free energy is quite insensitive to substituent size.
Collapse
Affiliation(s)
- Jan M Schümann
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Lukas Ochmann
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Jonathan Becker
- Institute of Analytical and Inorganic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Ingolf Harden
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Dipartmento di Chimica, Biologia e Biotechnologie, Università Degli Studi Di Perugia, Via Elce di sotto, 8, 06123 Perugia, Italy
| | - Peter R Schreiner
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| |
Collapse
|
21
|
Nakai H, Kobayashi M, Yoshikawa T, Seino J, Ikabata Y, Nishimura Y. Divide-and-Conquer Linear-Scaling Quantum Chemical Computations. J Phys Chem A 2023; 127:589-618. [PMID: 36630608 DOI: 10.1021/acs.jpca.2c06965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.
Collapse
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo, Hokkaido060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo, Hokkaido001-0021, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba274-8510, Japan
| | - Junji Seino
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Yasuhiro Ikabata
- Information and Media Center, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan.,Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| |
Collapse
|
22
|
de Oliveira MT, Alves JMA, Vrech NL, Braga AAC, Barboza CA. A comprehensive benchmark investigation of quantum chemical methods for carbocations. Phys Chem Chem Phys 2023; 25:1903-1922. [PMID: 36541431 DOI: 10.1039/d2cp04603b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.
Collapse
Affiliation(s)
- Marcelo T de Oliveira
- Department of Chemistry and Physics, La Trobe Institute of Molecular Sciences, La Trobe University, Melbourne, Victoria 3086, Australia. .,Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Júlia M A Alves
- Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Natália L Vrech
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Ataualpa A C Braga
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Cristina A Barboza
- Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland.,Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, Warsaw, Poland
| |
Collapse
|
23
|
Santra G, Martin JM. Performance of Localized-Orbital Coupled-Cluster Approaches for the Conformational Energies of Longer n-Alkane Chains. J Phys Chem A 2022; 126:9375-9391. [PMID: 36508714 PMCID: PMC9791657 DOI: 10.1021/acs.jpca.2c06407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We report an update and enhancement of the ACONFL (conformer energies of large alkanes [J. Phys. Chem. A2022,126, 3521-3535]) dataset. For the ACONF12 (n-dodecane) subset, we report basis set limit canonical coupled-cluster with singles, doubles, and perturbative triples [i.e., CCSD(T)] reference data obtained from the MP2-F12/cc-pV{T,Q}Z-F12 extrapolation, [CCSD(F12*)-MP2-F12]/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/aug-cc-pV{D,T}Z calculations. Then, we explored the performance of a variety of single and composite localized-orbital CCSD(T) approximations, ultimately finding an affordable localized natural orbital CCSD(T) [LNO-CCSD(T)]-based post-MP2 correction that agrees to 0.006 kcal/mol mean absolute deviation with the revised canonical reference data. In tandem with canonical MP2-F12 complete basis set extrapolation, this was then used to re-evaluate the ACONF16 and ACONF20 subsets for n-hexadecane and n-icosane, respectively. Combining those with the revised canonical reference data for the dodecane conformers (i.e., ACONF12 subset), a revised ACONFL set was obtained. It was then used to assess the performance of different localized-orbital coupled-cluster approaches, such as pair natural orbital localized CCSD(T) [PNO-LCCSD(T)] as implemented in MOLPRO, DLPNO-CCSD(T0) and DLPNO-CCSD(T1) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, at their respective "Normal", "Tight", "vTight", and "vvTight" accuracy settings. For a given accuracy threshold and basis set, DLPNO-CCSD(T1) and DLPNO-CCSD(T0) perform comparably. With "VeryTightPNO" cutoffs, explicitly correlated DLPNO-CCSD(T1)-F12/VDZ-F12 is the best pick among all the DLPNO-based methods tested. To isolate basis set incompleteness from localized-orbital-related truncation errors (domain, LNOs), we have also compared the localized coupled-cluster approaches with canonical DF-CCSD(T)/aug-cc-pVTZ for the ACONF12 set. We found that gradually tightening the cutoffs improves the performance of LNO-CCSD(T), and using a composite scheme such as vTight + 0.50[vTight - Tight] improves things further. For DLPNO-CCSD(T1), "TightPNO" and "VeryTightPNO" offer a statistically similar accuracy, which gets slightly better when TCutPNO is extrapolated to the complete PNO space limit. Similar to Brauer et al.'s [Phys. Chem. Chem. Phys.2016,18 (31), 20905-20925] previous report for the S66x8 noncovalent interactions, the dispersion-corrected direct random phase approximation (dRPA)-based double hybrids perform remarkably well for the ACONFL set. While the revised reference data do not affect any conclusions on the less accurate methods, they may upend orderings for more accurate methods with error statistics on the same order as the difference between reference datasets.
Collapse
|
24
|
Reimann M, Kaupp M. Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes. J Chem Theory Comput 2022; 18:7442-7456. [PMID: 36417564 DOI: 10.1021/acs.jctc.2c00924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = Equintet - Esinglet) of [FeIIL6]2+ complexes (L = CNH, CO, NCH, NH3, H2O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH3, H2O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.
Collapse
Affiliation(s)
- Marc Reimann
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
25
|
A Cost Effective Scheme for the Highly Accurate Description of Intermolecular Binding in Large Complexes. Int J Mol Sci 2022; 23:ijms232415773. [PMID: 36555413 PMCID: PMC9780852 DOI: 10.3390/ijms232415773] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Revised: 11/23/2022] [Accepted: 12/07/2022] [Indexed: 12/15/2022] Open
Abstract
There has been a growing interest in quantitative predictions of the intermolecular binding energy of large complexes. One of the most important quantum chemical techniques capable of such predictions is the domain-based local pair natural orbital (DLPNO) scheme for the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)], whose results are extrapolated to the complete basis set (CBS) limit. Here, the DLPNO-based focal-point method is devised with the aim of obtaining CBS-extrapolated values that are very close to their canonical CCSD(T)/CBS counterparts, and thus may serve for routinely checking a performance of less expensive computational methods, for example, those based on the density-functional theory (DFT). The efficacy of this method is demonstrated for several sets of noncovalent complexes with varying amounts of the electrostatics, induction, and dispersion contributions to binding (as revealed by accurate DFT-based symmetry-adapted perturbation theory (SAPT) calculations). It is shown that when applied to dimeric models of poly(3-hydroxybutyrate) chains in its two polymorphic forms, the DLPNO-CCSD(T) and DFT-SAPT computational schemes agree to within about 2 kJ/mol of an absolute value of the interaction energy. These computational schemes thus should be useful for a reliable description of factors leading to the enthalpic stabilization of extended systems.
Collapse
|
26
|
Bready CJ, Fowler VR, Juechter LA, Kurfman LA, Mazaleski GE, Shields GC. The driving effects of common atmospheric molecules for formation of prenucleation clusters: the case of sulfuric acid, formic acid, nitric acid, ammonia, and dimethyl amine. ENVIRONMENTAL SCIENCE: ATMOSPHERES 2022; 2:1469-1486. [PMID: 36561556 PMCID: PMC9648633 DOI: 10.1039/d2ea00087c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Accepted: 09/30/2022] [Indexed: 11/12/2022]
Abstract
How secondary aerosols form is critical as aerosols' impact on Earth's climate is one of the main sources of uncertainty for understanding global warming. The beginning stages for formation of prenucleation complexes, that lead to larger aerosols, are difficult to decipher experimentally. We present a computational chemistry study of the interactions between three different acid molecules and two different bases. By combining a comprehensive search routine covering many thousands of configurations at the semiempirical level with high level quantum chemical calculations of approximately 1000 clusters for every possible combination of clusters containing a sulfuric acid molecule, a formic acid molecule, a nitric acid molecule, an ammonia molecule, a dimethylamine molecule, and 0-5 water molecules, we have completed an exhaustive search of the DLPNO-CCSD(T)/CBS//ωB97X-D/6-31++G** Gibbs free energy surface for this system. We find that the detailed geometries of each minimum free energy cluster are often more important than traditional acid or base strength. Addition of a water molecule to a dry cluster can enhance stabilization, and we find that the (SA)(NA)(A)(DMA)(W) cluster has special stability. Equilibrium calculations of SA, FA, NA, A, DMA, and water using our quantum chemical ΔG° values for cluster formation and realistic estimates of the concentrations of these monomers in the atmosphere reveals that nitric acid can drive early stages of particle formation just as efficiently as sulfuric acid. Our results lead us to believe that particle formation in the atmosphere results from the combination of many different molecules that are able to form highly stable complexes with acid molecules such as SA, NA, and FA.
Collapse
Affiliation(s)
- Conor J. Bready
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| | - Vance R. Fowler
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| | - Leah A. Juechter
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| | - Luke A. Kurfman
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| | - Grace E. Mazaleski
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| | - George C. Shields
- Department of Chemistry, Furman UniversityGreenvilleSouth Carolina 29613USA
| |
Collapse
|
27
|
Wang Y, Ni Z, Neese F, Li W, Guo Y, Li S. Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems. J Chem Theory Comput 2022; 18:6510-6521. [PMID: 36240189 DOI: 10.1021/acs.jctc.2c00412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H2O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou311121, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der RuhrD-45470, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong266237, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| |
Collapse
|
28
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
29
|
Chan B, Karton A. Assessment of DLPNO-CCSD(T)-F12 and its use for the formulation of the low-cost and reliable L-W1X composite method. J Comput Chem 2022; 43:1394-1402. [PMID: 35709311 DOI: 10.1002/jcc.26892] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 03/28/2022] [Accepted: 04/28/2022] [Indexed: 02/03/2023]
Abstract
In the present study, we have investigated the performance of RIJCOSX DLPNO-CCSD(T)-F12 methods for a wide range of systems. Calculations with a high-accuracy option ["DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12"] extrapolated to the complete-basis-set limit using the maug-cc-pV[D+d,T+d]Z basis sets provides fairly good agreements with the canonical CCSD(T)/CBS reference for a diverse set of thermochemical and kinetic properties [with mean absolute deviations (MADs) of ~1-2 kJ mol-1 except for atomization energies]. On the other hand, the low-cost "RIJCOSX DLPNO-CCSD(T)-F12D" option leads to substantial deviations for certain properties, notably atomization energies (MADs of up to tens of kJ mol-1 ). With the high-accuracy CBS approach, we have formulated the L-W1X method, which further includes a low-cost core-valence plus scalar-relativistic term. It shows generally good accuracy. For improved accuracies in specific cases, we advise replacing maug-cc-pV(n+d)Z with jun-cc-pV(n+d)Z for the calculation of electron affinities, and using well-constructed isodesmic-type reactions to obtain atomization energies. For medium-sized systems, DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 calculations are several times faster than the corresponding canonical computation; the use of the local approximations (RIJCOSX and DLPNO) leads to a better scaling than that for the canonical calculation (from ~6-7 down to ~2-4 for our test systems). Thus, the DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 method, and the L-W1X protocol that based on it, represent a useful means for obtaining accurate thermochemical quantities for larger systems.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Nagasaki, Japan
| | - Amir Karton
- School of Molecular Sciences, University of Western Australia, Perth, WA, Australia
| |
Collapse
|
30
|
Pliego JR. Copper-Catalyzed Aromatic Fluorination of 2-(2-bromophenyl)pyridine via Cu(I)/Cu(III) Mechanism in Acetonitrile Solvent: Cluster-Continuum Free Energy Profile and Microkinetic Analysis. MOLECULAR CATALYSIS 2022. [DOI: 10.1016/j.mcat.2022.112560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
31
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
32
|
Bruno G, de Souza B, Neese F, Bistoni G. Can domain-based local pair natural orbitals approaches accurately predict phosphorescence energies? Phys Chem Chem Phys 2022; 24:14228-14241. [PMID: 35649286 DOI: 10.1039/d2cp01623k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Since the discovery of the peculiar conducting and optical properties of aromatics, many efforts have been made to characterize and predict their phosphorescence. This physical process is exploited in modern Organic Emitting Light Diodes (OLEDs), and it is also one of the processes decreasing the efficiency of Dye-sensitized solar cells (DSSCs). Herein, we propose a computational strategy for the accurate calculation of singlet-triplet gaps of aromatic compounds, which provides results that are in excellent agreement with available experimental data. Our approach relies on the domain-based local pair natural orbital (DLPNO) variant of the "gold standard" CCSD(T) method. The convergence of our results with respect to the key technical parameters of the calculation, such as the basis set used, the approximations employed in the perturbative triples correction, and the dimension of the PNOs space, was thoroughly discussed.
Collapse
Affiliation(s)
- Giovanna Bruno
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany.,Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Perugia, Italy.
| |
Collapse
|
33
|
Drosou M, Mitsopoulou CA, Pantazis DA. Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics. J Chem Theory Comput 2022; 18:3538-3548. [PMID: 35582788 PMCID: PMC9202354 DOI: 10.1021/acs.jctc.2c00265] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
Spin-state energetics
of transition metal complexes remain one
of the most challenging targets for electronic structure methods.
Among single-reference wave function approaches, local correlation
approximations to coupled cluster theory, most notably the domain-based
local pair natural orbital (DLPNO) approach, hold the promise of bringing
the accuracy of coupled cluster theory with single, double, and perturbative
triple excitations, CCSD(T), to molecular systems of realistic size
with acceptable computational cost. However, recent studies on spin-state
energetics of iron-containing systems raised doubts about the ability
of the DLPNO approach to adequately and systematically approximate
energetics obtained by the reference-quality complete active space
second-order perturbation theory with coupled-cluster semicore correlation,
CASPT2/CC. Here, we revisit this problem using a diverse set of iron
complexes and examine several aspects of the application of the DLPNO
approach. We show that DLPNO-CCSD(T) can accurately reproduce both
CASPT2/CC and canonical CCSD(T) results if two basic principles are
followed. These include the consistent use of the improved iterative
(T1) versus the semicanonical perturbative triple corrections
and, most importantly, a simple two-point extrapolation to the PNO
space limit. The latter practically eliminates errors arising from
the default truncation of electron-pair correlation spaces and should
be viewed as standard practice in applications of the method to transition
metal spin-state energetics. Our results show that reference-quality
results can be readily achieved with DLPNO-CCSD(T) if these principles
are followed. This is important also in view of the applicability
of the method to larger single-reference systems and multinuclear
clusters, whose treatment of dynamic correlation would be challenging
for multireference-based approaches.
Collapse
Affiliation(s)
- Maria Drosou
- Inorganic Chemistry Laboratory, National and Kapodistrian University of Athens, Panepistimiopolis, Zografou 15771, Greece
| | - Christiana A Mitsopoulou
- Inorganic Chemistry Laboratory, National and Kapodistrian University of Athens, Panepistimiopolis, Zografou 15771, Greece
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
34
|
Altun A, Neese F, Bistoni G. Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies. J Chem Theory Comput 2022; 18:2292-2307. [PMID: 35167304 PMCID: PMC9009084 DOI: 10.1021/acs.jctc.1c01295] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The London dispersion
(LD)-corrected Hartree–Fock (HF) method
(HFLD) is an ab initio approach for the quantification
and analysis of noncovalent interactions (NCIs) in large systems that
is based on the domain-based local pair natural orbital coupled-cluster
(DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation,
accuracy, and efficiency of its closed-shell variant. Herein, an extension
of this method to open-shell molecular systems is presented. Its accuracy
is tested on challenging benchmark sets for NCIs, using CCSD(T) energies
at the estimated complete basis set limit as reference. The HFLD scheme
was found to be as accurate as the best-performing dispersion-corrected
exchange-correlation functionals, while being nonempirical and equally
efficient. In addition, it can be combined with the well-established
local energy decomposition (LED) for the analysis of NCIs, thus yielding
additional physical insights.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany.,Department of Chemistry, Biology and Biotechnology, University of Perugia, 06123 Perugia, Italy
| |
Collapse
|
35
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- 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
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth 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
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
36
|
Neese F. Software update: The
ORCA
program system—Version 5.0. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1606] [Citation(s) in RCA: 49] [Impact Index Per Article: 24.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Frank Neese
- Max Planck Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|
37
|
Affiliation(s)
- Milica Feldt
- Leibniz Institute for Catalysis: Leibniz-Institut fur Katalyse eV Theory & Catalysis Albert-Einstein-Str 29A 18059 Rostock GERMANY
| | - Quan Manh Phung
- Nagoya University: Nagoya Daigaku Department of Chemistry JAPAN
| |
Collapse
|
38
|
Bhattacharjee S, Isegawa M, Garcia-Ratés M, Neese F, Pantazis DA. Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches. J Chem Theory Comput 2022; 18:1619-1632. [PMID: 35191695 PMCID: PMC8908766 DOI: 10.1021/acs.jctc.1c01267] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
![]()
Hydrated transition
metal ions are prototypical systems that can
be used to model properties of transition metals in complex chemical
environments. These seemingly simple systems present challenges for
computational chemistry and are thus crucial in evaluations of quantum
chemical methods for spin-state and redox energetics. In this work,
we explore the applicability of the domain-based pair natural orbital
implementation of coupled cluster (DLPNO-CC) theory to the calculation
of ionization energies and redox potentials for hydrated ions of all
first transition row (3d) metals in the 2+/3+ oxidation states, in
connection with various solvation approaches. In terms of model definition,
we investigate the construction of a minimally explicitly hydrated
quantum cluster with a first and second hydration layer. We report
on the convergence with respect to the coupled cluster expansion and
the PNO space, as well as on the role of perturbative triple excitations.
A recent implementation of the conductor-like polarizable continuum
model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent
redox potentials at the coupled cluster level. Our results establish
conditions for the convergence of DLPNO-CCSD(T) energetics and stress
the absolute necessity to explicitly consider the second solvation
sphere even when CPCM is used. The achievable accuracy for redox potentials
of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore,
multilayer approaches that combine a higher-level DLPNO-CCSD(T) description
of the first solvation sphere with a lower-level description of the
second solvation layer are investigated. The present work establishes
optimal and transferable methodological choices for employing DLPNO-based
coupled cluster theory, the associated CPCM implementation, and cost-efficient
multilayer derivatives of the approach for open-shell transition metal
systems in complex environments.
Collapse
Affiliation(s)
- Sinjini Bhattacharjee
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Miho Isegawa
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Miquel Garcia-Ratés
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
39
|
Semidalas E, Martin JM. The MOBH35 Metal–Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes. J Chem Theory Comput 2022; 18:883-898. [PMID: 35045709 PMCID: PMC8830049 DOI: 10.1021/acs.jctc.1c01126] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We have revisited
the MOBH35 (Metal–Organic Barrier Heights,
35 reactions) benchmark [Iron, Janes, J.
Phys. Chem. A, 2019, 123 ( (17), ), 3761−378130973722; ibid. 2019, 123, 6379–6380] for realistic organometallic catalytic reactions, using both canonical
CCSD(T) and localized orbital approximations to it. For low levels
of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T)
perform well; for moderately strong levels of static correlation,
DLPNO-CCSD(T) and (T1) may break down catastrophically,
and PNO-LCCSD(T) is vulnerable as well. In contrast, LNO-CCSD(T) converges
smoothly to the canonical CCSD(T) answer with increasingly tight convergence
settings. The only two reactions for which our revised MOBH35 reference
values differ substantially from the original ones are reaction 9
and to a lesser extent 8, both involving iron. For the purpose of
evaluating density functional theory (DFT) methods for MOBH35, it
would be best to remove reaction 9 entirely as its severe level of
static correlation makes it just too demanding for a test. The magnitude
of the difference between DLPNO-CCSD(T) and DLPNO-CCSD(T1) is a reasonably good predictor for errors in DLPNO-CCSD(T1) compared to canonical CCSD(T); otherwise, monitoring all of T1, D1, max|tiA|, and 1/(εLUMO – εHOMO) should provide adequate warning
for potential problems. Our conclusions are not specific to the def2-SVP
basis set but are largely conserved for the larger def2-TZVPP, as
they are for the smaller def2-SV(P): the latter may be an economical
choice for calibrating against canonical CCSD(T). Finally, diagnostics
for static correlation are statistically clustered into groups corresponding
to (1) importance of single excitations in the wavefunction; (2a)
the small band gap, weakly separated from (2b) correlation entropy;
and (3) thermochemical importance of correlation energy, as well as
the slope of the DFT reaction energy with respect to the percentage
of HF exchange. Finally, a variable reduction analysis reveals that
much information on the multireference character is provided by T1, IND/Itot, and the exchange-based diagnostic A100[TPSS].
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| | - Jan M.L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| |
Collapse
|
40
|
Santra G, Semidalas E, Mehta N, Karton A, Martin JML. S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods. Phys Chem Chem Phys 2022; 24:25555-25570. [DOI: 10.1039/d2cp03938a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The S66x8 noncovalent interactions benchmark has been re-evaluated at the “sterling silver” level. Against this, a selection of computationally more economical alternatives has been assayed, ranging from localized CC to double hybrids and SAPT(DFT).
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
- School of Science and Technology, University of New England, Armidale, NSW 2351, Australia
| | - Jan M. L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
41
|
Altun A, Ghosh S, Riplinger C, Neese F, Bistoni G. Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations. J Phys Chem A 2021; 125:9932-9939. [PMID: 34730360 PMCID: PMC8607505 DOI: 10.1021/acs.jpca.1c09106] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Over the last two decades, the local approximation has been successfully used to extend the range of applicability of the "gold standard" singles and doubles coupled-cluster method with perturbative triples CCSD(T) to systems with hundreds of atoms. The local approximation error grows in absolute value with the increasing system size, i.e., by increasing the number of electron pairs in the system. In this study, we demonstrate that the recently introduced two-point extrapolation scheme for approaching the complete pair natural orbital (PNOs) space limit in domain-based pair natural orbital CCSD(T) calculations drastically reduces the dependence of the error on the system size, thus opening up unprecedented opportunities for the calculation of benchmark quality relative energies for large systems.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Soumen Ghosh
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
42
|
Drosou M, Mitsopoulou CA, Pantazis DA. Spin-state energetics of manganese spin crossover complexes: Comparison of single-reference and multi-reference ab initio approaches. Polyhedron 2021. [DOI: 10.1016/j.poly.2021.115399] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
43
|
Maurer LR, Bursch M, Grimme S, Hansen A. Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions. J Chem Theory Comput 2021; 17:6134-6151. [PMID: 34546754 DOI: 10.1021/acs.jctc.1c00659] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r2SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.
Collapse
Affiliation(s)
- Leonard R Maurer
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Markus Bursch
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
44
|
Drosou M, Pantazis DA. Redox Isomerism in the S 3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory. Chemistry 2021; 27:12815-12825. [PMID: 34288176 PMCID: PMC8518824 DOI: 10.1002/chem.202101567] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2021] [Indexed: 01/19/2023]
Abstract
The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.
Collapse
Affiliation(s)
- Maria Drosou
- Inorganic Chemistry LaboratoryNational and Kapodistrian University of AthensPanepistimiopolisZografou15771Greece
| | - Dimitrios A. Pantazis
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an derRuhrGermany
| |
Collapse
|
45
|
Affiliation(s)
- Frank Jensen
- Department of Chemistry Aarhus University DK-8000 Aarhus Denmark
| |
Collapse
|
46
|
Ghosh S, Neese F, Izsák R, Bistoni G. Fragment-Based Local Coupled Cluster Embedding Approach for the Quantification and Analysis of Noncovalent Interactions: Exploring the Many-Body Expansion of the Local Coupled Cluster Energy. J Chem Theory Comput 2021; 17:3348-3359. [PMID: 34037397 PMCID: PMC8190956 DOI: 10.1021/acs.jctc.1c00005] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Indexed: 11/30/2022]
Abstract
Herein, we introduce a fragment-based local coupled cluster embedding approach for the accurate quantification and analysis of noncovalent interactions in molecular aggregates. Our scheme combines two different expansions of the domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) energy: the many-body expansion (MBE) and the local energy decomposition (LED). The low-order terms in the MBE are initially computed in the presence of an environment that is treated at a low level of theory. Then, LED is used to decompose the energy of each term in the embedded MBE into additive fragment and fragment-pairwise contributions. This information is used to quantify the total energy of the system while providing at the same time in-depth insights into the nature and cooperativity of noncovalent interactions. Two different approaches are introduced and tested, in which the environment is treated at different levels of theory: the local coupled cluster in the Hartree-Fock (LCC-in-HF) method, in which the environment is treated at the HF level; and the electrostatically embedded local coupled cluster method (LCC-in-EE), in which the environment is replaced by point charges. Both schemes are designed to preserve as much as possible the accuracy of the parent local coupled cluster method for total energies, while being embarrassingly parallel and less memory intensive. These schemes appear to be particularly promising for the study of large and complex molecular aggregates at the coupled cluster level, such as condensed phase systems and protein-ligand interactions.
Collapse
Affiliation(s)
- Soumen Ghosh
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|