1
|
Mester D, Nagy PR, Kállay M. Basis-Set Limit CCSD(T) Energies for Large Molecules with Local Natural Orbitals and Reduced-Scaling Basis-Set Corrections. J Chem Theory Comput 2024; 20:7453-7468. [PMID: 39207805 PMCID: PMC11391584 DOI: 10.1021/acs.jctc.4c00777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The calculation of density-based basis-set correction (DBBSC), which remedies the basis-set incompleteness (BSI) error of the correlation energy, is combined with local approximations. Aiming at large-scale applications, the procedure is implemented in our efficient local natural orbital-based coupled-cluster singles and doubles with perturbative triples [LNO-CCSD(T)] scheme. To this end, the range-separation function, which characterizes the one-electron BSI in space, is decomposed into the sum of contributions from individual localized molecular orbitals (LMOs). A compact domain is constructed around each LMO, and the corresponding contributions are evaluated only within these restricted domains. Furthermore, for the calculation of the complementary auxiliary basis set (CABS) correction, which significantly improves the Hartree-Fock (HF) energy, the local density fitting approximation is utilized. The errors arising from the local approximations are examined in detail, efficient prescreening techniques are introduced to compress the numerical quadrature used for DBBSC, and conservative default thresholds are selected for the truncation parameters. The efficiency of the DBBSC-LNO-CCSD(T) method is demonstrated through representative examples of up to 1000 atoms. Based on the numerical results, we conclude that the corrections drastically reduce the BSI error using double-ζ basis sets, often to below 1 kcal/mol compared to the reliable LNO-CCSD(T) complete basis set references, while significant improvements are also achieved with triple-ζ basis sets. Considering that the calculation of the DBBSC and CABS corrections only moderately increases the wall-clock time required for the post-HF steps in practical applications, the proposed DBBSC-LNO-CCSD(T) method offers a highly efficient and robust tool for large-scale calculations.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Muegyetem rkp. 3, H-1111 Budapest, Hungary
| |
Collapse
|
2
|
Ünal A, Bozkaya U. Equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation. J Comput Chem 2024. [PMID: 39235313 DOI: 10.1002/jcc.27495] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Revised: 06/13/2024] [Accepted: 08/18/2024] [Indexed: 09/06/2024]
Abstract
Orbital-optimized coupled-cluster methods are very helpful for theoretical predictions of the molecular properties of challenging chemical systems, such as excited states. In this research, an efficient implementation of the equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting (DF) approach, denoted by DF-EOM-OCCD, is presented. The computational cost of the DF-EOM-OCCD method for excitation energies is compared with that of the conventional EOM-OCCD method. Our results demonstrate that DF-EOM-OCCD excitation energies are dramatically accelerated compared to EOM-OCCD. There are almost 17-fold reductions for theC 5 H 12 $$ {\mathrm{C}}_5{\mathrm{H}}_{12} $$ molecule in an aug-cc-pVTZ basis set with the RHF reference. This dramatic performance improvement comes from the reduced cost of integral transformation with the DF approach and the efficient evaluation of the particle-particle ladder (PPL) term, which is the most expensive term to evaluate. Further, our results show that the DF-EOM-OCCD approach is very helpful for the computation of excitation energies in open-shell molecular systems. Overall, we conclude that our new DF-EOM-OCCD implementation is very promising for the study of excited states in large-sized challenging chemical systems.
Collapse
Affiliation(s)
- Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara, Turkey
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara, Turkey
| |
Collapse
|
3
|
Kats D, Christlmaier EMC, Schraivogel T, Alavi A. Orbital optimisation in xTC transcorrelated methods. Faraday Discuss 2024. [PMID: 39072553 DOI: 10.1039/d4fd00036f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Additionally, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed.
Collapse
Affiliation(s)
- Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | | | - Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| |
Collapse
|
4
|
Zhang X, Li C, Ye HZ, Berkelbach TC, Chan GKL. Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics. J Chem Phys 2024; 161:014109. [PMID: 38949583 DOI: 10.1063/5.0212274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.
Collapse
Affiliation(s)
- Xing Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chenghan Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
5
|
Murillo-Lopez JA, Villegas-Escobar N, Vogt-Geisse S, Vöhringer-Martinez E. Molecular Environment Modulates CO 2 Liberation from Carboxy-Biotin. J Phys Chem B 2024; 128:5327-5335. [PMID: 38771940 DOI: 10.1021/acs.jpcb.4c00910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
Carboxy-biotin serves as a coenzyme in certain carboxylases, exhibiting the remarkable capability to transfer a carboxy group to specific substrates. This process is made possible by the presence of biotin, a unique molecule that consists of a sulfur-containing tetrahydrothiophene ring fused to a ureido group. It is covalently attached to the enzyme via a flexible linker, allowing for its functionality. Biotin-dependent carboxylases consist of two distinct domains. The first domain (BC) facilitates biotin carboxylation by utilizing ATP, while the second domain (CT) transfers CO2 to the substrate. The process of ATP-dependent carboxylation using bicarbonate in the biotin carboxylase domain (BC) is well-known. However, the precise mechanism by which CO2 is released in the carboxyltransferase domain (CT) is still not fully understood. We employed advanced computational chemistry methods to investigate the decarboxylation process of carboxy-biotin in various molecular environments and different protonation states. Regardless of the polarity of the molecular surroundings, decarboxylation only occurs spontaneously in the protonated form. To determine the protonation state of biotin in different environments, we established an accurate computational chemistry method for calculating the pKa value of carboxy-biotin, reaching sub-kcal/mol accuracy. Based on our findings, nonpolar environments, such as the active site of the carboxyltransferase domain, have the ability to cause the spontaneous release of CO2 from carboxy-biotin. The CO2 release takes place spontaneously from protonated carboxy-biotin, promoting the carboxylation of substrates.
Collapse
Affiliation(s)
- Juliana A Murillo-Lopez
- Departamento de Físico-Química Facultad de Ciencias Químicas, Universidad de Concepción, Concepción 4070139, Chile
| | - Nery Villegas-Escobar
- Departamento de Físico-Química Facultad de Ciencias Químicas, Universidad de Concepción, Concepción 4070139, Chile
| | - Stefan Vogt-Geisse
- Departamento de Físico-Química Facultad de Ciencias Químicas, Universidad de Concepción, Concepción 4070139, Chile
| | - Esteban Vöhringer-Martinez
- Departamento de Físico-Química Facultad de Ciencias Químicas, Universidad de Concepción, Concepción 4070139, Chile
| |
Collapse
|
6
|
Zheng Y, Ni Z, Wang Y, Li W, Li S. Analytical Gradient Using Cluster-in-Molecule RI-MP2 Method for the Geometry Optimizations of Large Systems. J Chem Theory Comput 2024; 20:3626-3636. [PMID: 38626287 DOI: 10.1021/acs.jctc.4c00087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
We present an efficient analytical energy gradient algorithm for the cluster-in-molecule resolution-of-identity second-order Møller-Plesset perturbation (CIM-RI-MP2) method based on the Lagrange multiplier method. Our algorithm independently constructs the Lagrangian formalism within each cluster, avoiding the solution of the coupled-perturbed Hartree-Fock (CPHF) equation for the whole system. Due to this feature, the computational cost of the CIM-RI-MP2 gradients is much lower than that of other local MP2 algorithms. Benchmark calculations of several molecules containing up to 312 atoms demonstrate the general applicability of our CIM-RI-MP2 gradient algorithm. The optimized structure of a 244-atom molecule using the CIM-RI-MP2 method with the cc-pVDZ basis set is in good agreement with the corresponding crystal structure. A single-point gradient calculation conducted for a molecular cage containing 972 atoms and 9612 basis functions takes 48 h on 25 nodes, utilizing a total of 600 CPU cores. The present CIM-RI-MP2 gradient program is applicable for obtaining the optimized geometries of large systems with hundreds of atoms.
Collapse
Affiliation(s)
- Yang Zheng
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, 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, Hangzhou 311121, P. R. China
| | - Yuqi Wang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
Amanollahi Z, Lampe L, Bensberg M, Neugebauer J, Feldt M. On the accuracy of orbital based multi-level approaches for closed-shell transition metal chemistry. Phys Chem Chem Phys 2023; 25:4635-4648. [PMID: 36662158 DOI: 10.1039/d2cp05056k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
In this work, we investigate the accuracy of the local molecular orbital molecular orbital (LMOMO) scheme and projection-based wave function-in-density functional theory (WF-in-DFT) embedding for the prediction of reaction energies and barriers of typical reactions involving transition metals. To analyze the dependence of the accuracy on the system partitioning, we apply a manual orbital selection for LMOMO as well as the so-called direct orbital selection (DOS) for both approaches. We benchmark these methods on 30 closed shell reactions involving 16 different transition metals. This allows us to devise guidelines for the manual selection as well as settings for the DOS that provide accurate results within an error of 2 kcal mol-1 compared to local coupled cluster. To reach this accuracy, on average 55% of the occupied orbitals have to be correlated with coupled cluster for the current test set. Furthermore, we find that LMOMO gives more reliable relative energies for small embedded regions than WF-in-DFT embedding.
Collapse
Affiliation(s)
- Zohreh Amanollahi
- Leibniz Institute for Catalysis (LIKAT), Albert-Einstein-Str. 29A, 18059 Rostock, Germany.
| | - Lukas Lampe
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Moritz Bensberg
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Milica Feldt
- Leibniz Institute for Catalysis (LIKAT), Albert-Einstein-Str. 29A, 18059 Rostock, Germany.
| |
Collapse
|
12
|
Kállay M, Horváth RA, Gyevi-Nagy L, Nagy PR. Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach. J Chem Theory Comput 2022; 19:174-189. [PMID: 36576419 PMCID: PMC9835832 DOI: 10.1021/acs.jctc.2c01031] [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/29/2022]
Abstract
Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).
Collapse
Affiliation(s)
- Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,
| | - Réka A. Horváth
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - László Gyevi-Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
13
|
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
|
14
|
Whitten JL. Correction of Residual Errors in Configuration Interaction Electronic Structure Calculations. J Chem Phys 2022; 157:124101. [DOI: 10.1063/5.0098793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Methods for correcting residual energy errors of configuration interaction (CI) calculations of molecules and other electronic systems are discussed based on the assumption that the energy defect can be mapped onto atomic regions. The methods do not consider the detailed nature of excitations, but instead define a defect energy per electron that that is unique to a specific atom. Defect energy contributions are determined from calculations on diatomic and hydride molecules and then applied to other systems. Calculated energies are compared with experimental thermodynamic and spectroscopic data for a set of forty-one mainly organic molecules representing a wide range of bonding environments. The most stringent test is based on a severely truncated virtual space in which higher spherical harmonic basis functions are removed. The errors of the initial CI calculations are large, but in each case, including defect corrections brings calculated CI energies into agreement with experimental values. The method is also applied to a NIST compilation of coupled-cluster calculations that employ a larger basis set and no truncation of the virtual space. The corrections show excellent consistency with total energies in very good agreement with experimental values. An extension of the method is applied to dmsn states of Sc, Ti, V, Mn, Cr, Fe, Co, Ni and Cu, significantly improving the agreement of calculated transition energies with spectroscopic values.
Collapse
Affiliation(s)
- Jerry L Whitten
- North Carolina State University, North Carolina State University, United States of America
| |
Collapse
|
15
|
Nandy A, Adamji H, Kastner DW, Vennelakanti V, Nazemi A, Liu M, Kulik HJ. Using Computational Chemistry To Reveal Nature’s Blueprints for Single-Site Catalysis of C–H Activation. ACS Catal 2022. [DOI: 10.1021/acscatal.2c02096] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Husain Adamji
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - David W. Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Vyshnavi Vennelakanti
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Azadeh Nazemi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Mingjie Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
16
|
Bensberg M, Neugebauer J. Orbital Pair Selection for Relative Energies in the Domain-Based Local Pair Natural Orbital Coupled-Cluster Method. J Chem Phys 2022; 157:064102. [DOI: 10.1063/5.0100010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
For the accurate computation of relative energies, domain-based local pair natural orbital coupled-cluster [DLPNO-CCSD(T0)] has become increasingly popular. Even though DLPNO-CCSD(T0) shows a formally linear scaling of the computational effort with the system size, accurate predictions of relative energies remain costly. Therefore, multi-level approaches are attractive that focus the available computational resources on a minor part of the molecular system, e.g., a reaction center, where changes in the correlation energy are expected to be the largest. We present a pair-selected multi-level DLPNO-CCSD(T0) ansatz that automatically partitions the orbital pairs according to their contribution to the overall correlation energy change in a chemical reaction. To this end, the localized orbitals are mapped between structures in the reaction; all pair energies are approximated through computationally efficient semi-canonical second-order Møller--Plesser perturbation theory, and the orbital pairs for which the pair energies change significantly are identified. This multi-level approach is significantly more robust than our previously suggested, orbital selection-based multi-level DLPNO-CCSD(T0) ansatz [ J. Chem. Phys. 2021, 155, 224102] for reactions showing only small changes in the occupied orbitals. At the same time, it is even more efficient without added input complexity or accuracy loss compared to the full DLPNO-CCSD(T0) calculation. We demonstrate the accuracy of the multi-level approach for a total of 128 chemical reactions and potential energy curves of weakly interacting complexes from the S66x8 benchmark set.
Collapse
Affiliation(s)
- Moritz Bensberg
- Westfälische Wilhelms-Universität Münster Fachbereich 12 Chemie und Pharmazie, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster Fachbereich 12 Chemie und Pharmazie, Germany
| |
Collapse
|
17
|
Schnack-Petersen AK, Koch H, Coriani S, Kjønstad EF. Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals. J Chem Phys 2022; 156:244111. [DOI: 10.1063/5.0087261] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3 O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.
Collapse
Affiliation(s)
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Eirik F. Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
18
|
Lemke Y, Kussmann J, Ochsenfeld C. Efficient Integral-Direct Methods for Self-Consistent Reduced Density Matrix Functional Theory Calculations on Central and Graphics Processing Units. J Chem Theory Comput 2022; 18:4229-4244. [DOI: 10.1021/acs.jctc.2c00231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Y. Lemke
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
| | - J. Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
| | - C. Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 5−13, D-81377 Munich, Germany
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
19
|
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
|
20
|
Ünal A, Bozkaya U. Efficient Implementation of Equation-of-Motion Coupled-Cluster Singles and Doubles Method with the Density-Fitting Approximation: An Enhanced Algorithm for the Particle-Particle Ladder Term. J Chem Theory Comput 2022; 18:1489-1500. [PMID: 35107297 PMCID: PMC8908769 DOI: 10.1021/acs.jctc.1c01000] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
An efficient implementation of the density-fitted equation-of-motion coupled-cluster singles and doubles (DF-EOM-CCSD) method is presented with an enhanced algorithm for the particle-particle ladder (PPL) term, which is the most expensive part of EOM-CCSD computations. To further improve the evaluation of the PPL term, a hybrid density-fitting/Cholesky decomposition (DF/CD) algorithm is also introduced. In the hybrid DF/CD approach, four virtual index integrals are constructed on-the-fly from the DF factors; then, their partial Cholesky decomposition is simultaneously performed. The computational cost of the DF-EOM-CCSD method for excitation energies is compared with that of the resolution of the identity EOM-CCSD (RI-EOM-CCSD) (from the Q-chem 5.3 package). Our results demonstrate that DF-EOM-CCSD excitation energies are significantly accelerated compared to RI-EOM-CCSD. There is more than a 2-fold reduction for the C8H18 molecule in the cc-pVTZ basis set with the restricted Hartree-Fock (RHF) reference. This cost savings results from the efficient evaluation of the PPL term. In the RHF based DF-EOM-CCSD method, the number of flops (NOF) is 1/4O2V4, while that of RI-EOM-CCSD was reported (Epifanovsky et al. J. Chem. Phys. 2013, 139, 134105) to be 5/8O2V4 for the PPL contraction term. Further, the NOF of VVVV-type integral transformation is 1/2V4Naux in our case, while it appears to be V4Naux for RI-EOM-CCSD. Hence, our implementation is 2.5 and 2.0 times more efficient compared to RI-EOM-CCSD for these expensive terms. For the unrestricted Hartree-Fock (UHF) reference, our implementation maintains its enhanced performance and provides a 1.8-fold reduction in the computational time compared to RI-EOM-CCSD for the C7H16 molecule. Our results indicate that our DF-EOM-CCSD implementation is 1.7 and 1.4 times more efficient compared with RI-EOM-CCSD for average computational cost per EOM-CCSD iteration. Moreover, our results show that the new hybrid DF/CD approach improves upon the DF algorithm, especially for large molecular systems. Overall, we conclude that the new hybrid DF/CD PPL algorithm is very promising for large-sized chemical systems.
Collapse
Affiliation(s)
- Aslı Ünal
- Graduate School of Science and Engineering, Hacettepe University, Ankara 06800, Turkey.,Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
21
|
Bozkaya U, Ermiş B, Alagöz Y, Ünal A, Uyar AK. MacroQC 1.0: An electronic structure theory software for large-scale applications. J Chem Phys 2022; 156:044801. [DOI: 10.1063/5.0077823] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Betül Ermiş
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Yavuz Alagöz
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Ali Kaan Uyar
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
22
|
Mullan T, Maschio L, Saalfrank P, Usvyat D. Reaction barriers on non-conducting surfaces beyond periodic local MP2: Diffusion of hydrogen on \ce{\alpha-Al2O3}(0001) as a test case. J Chem Phys 2022; 156:074109. [DOI: 10.1063/5.0082805] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Lorenzo Maschio
- Dipartimento di Chimica, Università degli Studi di Torino, Italy
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam Institut für Chemie, Germany
| | - Denis Usvyat
- Institute of Chemistry, Humboldt University of Berlin, Germany
| |
Collapse
|
23
|
Bensberg M, Neugebauer J. Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method. J Chem Phys 2021; 155:224102. [PMID: 34911318 DOI: 10.1063/5.0071347] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.
Collapse
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
24
|
Dridi N, Jellali C, Hmida F, Hui L, Kwabia Tchana F, Landsheere X, Hammami K, Rotger M, Aroui H. Line intensity parameters, He-broadening and line shift coefficients in the 2v20 and 3v21 − v21 bands of OCS. Mol Phys 2021. [DOI: 10.1080/00268976.2021.2010821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- N. Dridi
- Laboratoire de Spectroscopie et Dynamique Moléculaire, Université de Tunis, Tunis, Tunisia
- Université de Reims Champagne Ardenne, Reims, France
| | - C. Jellali
- Laboratoire de Spectroscopie et Dynamique Moléculaire, Université de Tunis, Tunis, Tunisia
| | - F. Hmida
- Laboratoire de Spectroscopie et Dynamique Moléculaire, Université de Tunis, Tunis, Tunisia
- Université de Reims Champagne Ardenne, Reims, France
| | - L. Hui
- State Key Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun, People’s Republic of China
| | | | - X. Landsheere
- Université de Paris and Univ Paris Est Creteil, Paris, France
| | - K. Hammami
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications-LSAMA, Université de Tunis El Manar, Tunis, Tunisia
| | - M. Rotger
- Université de Reims Champagne Ardenne, Reims, France
| | - H. Aroui
- Laboratoire de Spectroscopie et Dynamique Moléculaire, Université de Tunis, Tunis, Tunisia
| |
Collapse
|
25
|
Liang Q, Yang J. Third-Order Many-Body Expansion of OSV-MP2 Wave Function for Low-Order Scaling Analytical Gradient Computation. J Chem Theory Comput 2021; 17:6841-6860. [PMID: 34704757 DOI: 10.1021/acs.jctc.1c00581] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present a many-body expansion (MBE) formulation and implementation for efficient computation of analytical energy gradients from the orbital-specific-virtual second-order Møllet-Plesset perturbation theory (OSV-MP2) based on our earlier work (Zhou et al. J. Chem. Theory Comput. 2020, 16, 196-210). The third-order MBE(3) expansion of OSV-MP2 amplitudes and density matrices was developed to adopt the orbital-specific clustering and long-range termination schemes, which avoids term-by-term differentiations of the MBE energy bodies. We achieve better efficiency by exploiting the algorithmic sparsity that allows us to prune out insignificant fitting integrals and OSV relaxations. With these approximations, the present implementation is benchmarked on a range of molecules that show an economic scaling in the linear and quadratic regimes for computing MBE(3)-OSV-MP2 amplitude and gradient equations, respectively, and yields normal accuracy comparable to the original OSV-MP2 results. The MPI-3-based parallelism through shared memory one-sided communication is further developed for improving parallel scalability and memory accessibility by sorting the MBE(3) orbital clusters into independent tasks that are distributed on multiple processes across many nodes, supporting both global and local data locations in which selected MBE(3)-OSV-MP2 intermediates of different sizes are distinguished and accordingly placed. The accuracy and efficiency level of our MBE(3)-OSV-MP2 analytical gradient implementation is finally illustrated in two applications: we show that the subtle coordination structure differences of mechanically interlocked Cu-catenane complexes can be distinguished when tuning ligand lengths; and the porphycene molecular dynamics reveals the emergence of the vibrational signature arising from softened N-H stretching associated with hydrogen transfer, using an MP2 level of electron correlation and classical nuclei for the first time.
Collapse
Affiliation(s)
- Qiujiang Liang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| |
Collapse
|
26
|
Cui K, Bhattacharyya S, Nair S, Schmidt JR. Origins of Acid-Gas Stability Behavior in Zeolitic Imidazolate Frameworks: The Unique High Stability of ZIF-71. J Am Chem Soc 2021; 143:18061-18072. [PMID: 34677971 DOI: 10.1021/jacs.1c06321] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Zeolitic imidazolate frameworks (ZIFs) are promising materials for industrial process separations, but recent literature reports have highlighted their vulnerability to acid gases (e.g., SO2, CO2, NO2, H2S), often present in practical applications. While previous work has documented the widely varying stability behavior of many ZIFs under varying (humid and dry) acid gas environments, efforts to explain or correlate these experimental observations via empirical descriptors have not succeeded. A key observation is that ZIF-71 (RHO topology) is an extraordinarily stable ZIF material, retaining both structure and porosity under prolonged humid SO2 exposure whereas many other well-known ZIFs with different linkers and topologies (such as ZIF-8) were shown to degrade. Through a combination of hybrid quantum mechanics/molecular mechanics (QM/MM) based methods and statistical mechanical models, we successfully explain this important experimental observation via atomistic investigations of the reaction mechanism. Our holistic approach reveals an ∼9 times lower average defect formation rate in ZIF-71 RHO compared to ZIF-8 SOD, leading to the conclusion that the observed experimental stability of this material rises from kinetic effects. Moreover, our analysis reveals that differing stability of the two materials is determined by the distributions of acid gas molecules, which is difficult to capture using empirical descriptors. Our results suggest wider applicability of the present approach, toward identifying tuned functional groups and topologies that move the acid gas distributions away from more reactive sites and thus allow enhanced kinetic stability.
Collapse
Affiliation(s)
- Kai Cui
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Souryadeep Bhattacharyya
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, 311 Ferst Drive NW, Atlanta, Georgia 30332, United States
| | - Sankar Nair
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, 311 Ferst Drive NW, Atlanta, Georgia 30332, United States
| | - J R Schmidt
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
27
|
Ray J, Samanta SK, Tripathy T. Adsorption of toxic organophosphorus pesticides from aqueous medium using dextrin‐graft‐
poly
(2‐acrylamido‐2‐methyl propane sulfonic acid‐co‐acrylic acid) copolymer: Studies on equilibrium kinetics, isotherms, and thermodynamics of interactions. POLYM ENG SCI 2021. [DOI: 10.1002/pen.25805] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Jagabandhu Ray
- Postgraduate Division of Chemistry Midnapore College (Autonomous) Midnapore India
| | - Santu Kumar Samanta
- Postgraduate Division of Chemistry Midnapore College (Autonomous) Midnapore India
| | - Tridib Tripathy
- Postgraduate Division of Chemistry Midnapore College (Autonomous) Midnapore India
| |
Collapse
|
28
|
Agarawal V, Patra C, Maitra R. An approximate coupled cluster theory via nonlinear dynamics and synergetics: The adiabatic decoupling conditions. J Chem Phys 2021; 155:124115. [PMID: 34598582 DOI: 10.1063/5.0060087] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The coupled cluster iteration scheme is analyzed as a multivariate discrete time map using nonlinear dynamics and synergetics. The nonlinearly coupled set of equations to determine the cluster amplitudes are driven by a fraction of the entire set of cluster amplitudes. These driver amplitudes enslave all other amplitudes through a synergistic inter-relationship, where the latter class of amplitudes behave as the auxiliary variables. The driver and the auxiliary variables exhibit vastly different time scales of relaxation during the iteration process to reach the fixed points. The fast varying auxiliary amplitudes are small in magnitude, while the driver amplitudes are large, and they have a much longer time scale of relaxation. Exploiting their difference in relaxation time scale, we employ an adiabatic decoupling approximation, where each of the fast relaxing auxiliary modes is expressed as a unique function of the principal amplitudes. This results in a tremendous reduction in the independent degrees of freedom. On the other hand, only the driver amplitudes are determined accurately via exact coupled cluster equations. We will demonstrate that the iteration scheme has an order of magnitude reduction in computational scaling than the conventional scheme. With a few pilot numerical examples, we would demonstrate that this scheme can achieve very high accuracy with significant savings in computational time.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Chayan Patra
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
29
|
Alagöz Y, Ünal A, Bozkaya U. Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation. J Chem Phys 2021; 155:114104. [PMID: 34551547 DOI: 10.1063/5.0061351] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OCCD) method with the density-fitting approach, denoted by DF-OCCD(T) and DF-OCCD(T)Λ, are presented. The computational cost of the DF-OCCD(T) method is compared with that of the conventional OCCD(T). In the conventional OCCD(T) and OCCD(T)Λ methods, one needs to perform four-index integral transformations at each coupled-cluster doubles iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD(T) provides dramatically lower computational costs compared to OCCD(T), and there are more than 68-fold reductions in the computational time for the C5H12 molecule with the cc-pVTZ basis set. Our results show that the DF-OCCD(T) and DF-OCCD(T)Λ methods are very helpful for the study of single bond-breaking problems. Performances of the DF-OCCD(T) and DF-OCCD(T)Λ methods are noticeably better than that of the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for the potential energy surfaces of the molecules considered. Specifically, the DF-OCCD(T)Λ method provides dramatic improvements upon CCSD(T), and there are 8-14-fold reductions in nonparallelity errors. Overall, we conclude that the DF-OCCD(T)Λ method is very promising for the study of challenging chemical systems, where the CCSD(T) fails.
Collapse
Affiliation(s)
- Yavuz Alagöz
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
30
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
31
|
Chen J, Kato J, Harper JB, Shao Y, Ho J. On the Accuracy of QM/MM Models: A Systematic Study of Intramolecular Proton Transfer Reactions of Amino Acids in Water. J Phys Chem B 2021; 125:9304-9316. [PMID: 34355564 DOI: 10.1021/acs.jpcb.1c04876] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol-1 of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jin Kato
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jason B Harper
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
32
|
Heinen S, von Rudorff GF, von Lilienfeld OA. Toward the design of chemical reactions: Machine learning barriers of competing mechanisms in reactant space. J Chem Phys 2021; 155:064105. [PMID: 34391351 DOI: 10.1063/5.0059742] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The interplay of kinetics and thermodynamics governs reactive processes, and their control is key in synthesis efforts. While sophisticated numerical methods for studying equilibrium states have well advanced, quantitative predictions of kinetic behavior remain challenging. We introduce a reactant-to-barrier (R2B) machine learning model that rapidly and accurately infers activation energies and transition state geometries throughout the chemical compound space. R2B exhibits improving accuracy as training set sizes grow and requires as input solely the molecular graph of the reactant and the information of the reaction type. We provide numerical evidence for the applicability of R2B for two competing text-book reactions relevant to organic synthesis, E2 and SN2, trained and tested on chemically diverse quantum data from the literature. After training on 1-1.8k examples, R2B predicts activation energies on average within less than 2.5 kcal/mol with respect to the coupled-cluster singles doubles reference within milliseconds. Principal component analysis of kernel matrices reveals the hierarchy of the multiple scales underpinning reactivity in chemical space: Nucleophiles and leaving groups, substituents, and pairwise substituent combinations correspond to systematic lowering of eigenvalues. Analysis of R2B based predictions of ∼11.5k E2 and SN2 barriers in the gas-phase for previously undocumented reactants indicates that on average, E2 is favored in 75% of all cases and that SN2 becomes likely for chlorine as nucleophile/leaving group and for substituents consisting of hydrogen or electron-withdrawing groups. Experimental reaction design from first principles is enabled due to R2B, which is demonstrated by the construction of decision trees. Numerical R2B based results for interatomic distances and angles of reactant and transition state geometries suggest that Hammond's postulate is applicable to SN2, but not to E2.
Collapse
Affiliation(s)
- Stefan Heinen
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
| | | | | |
Collapse
|
33
|
Symmetrized systematic molecular fragmentation model and its application for molecular properties. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113303] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
34
|
Stoychev GL, Auer AA, Gauss J, Neese F. DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings. J Chem Phys 2021; 154:164110. [PMID: 33940835 DOI: 10.1063/5.0047125] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.
Collapse
Affiliation(s)
- Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
35
|
Ye HZ, Berkelbach TC. Fast periodic Gaussian density fitting by range separation. J Chem Phys 2021; 154:131104. [PMID: 33832275 DOI: 10.1063/5.0046617] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present an efficient implementation of periodic Gaussian density fitting (GDF) using the Coulomb metric. The three-center integrals are divided into two parts by range-separating the Coulomb kernel, with the short-range part evaluated in real space and the long-range part in reciprocal space. With a few algorithmic optimizations, we show that this new method-which we call range-separated GDF (RSGDF)-scales sublinearly to linearly with the number of k-points for small to medium-sized k-point meshes that are commonly used in periodic calculations with electron correlation. Numerical results on a few three-dimensional solids show about ten-fold speedups over the previously developed GDF with little precision loss. The error introduced by RSGDF is about 10-5Eh in the converged Hartree-Fock energy with default auxiliary basis sets and can be systematically reduced by increasing the size of the auxiliary basis with little extra work.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
36
|
Lee SJR, Husch T, Ding F, Miller TF. Analytical gradients for molecular-orbital-based machine learning. J Chem Phys 2021; 154:124120. [PMID: 33810669 DOI: 10.1063/5.0040782] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.
Collapse
Affiliation(s)
- Sebastian J R Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Tamara Husch
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Feizhi Ding
- Entos, Inc., Los Angeles, California 90027, USA
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
37
|
Feldt M, Brown A. Assessment of local coupled cluster methods for excited states of BODIPY/Aza-BODIPY families. J Comput Chem 2021; 42:144-155. [PMID: 33103817 DOI: 10.1002/jcc.26442] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 10/05/2020] [Accepted: 10/10/2020] [Indexed: 12/28/2022]
Abstract
It was previously reported that Laplace transformed local CC2 (LCC2*) provided the best agreement (MAE = 0.145 eV) when comparing vertical excitation energies to experimental λmax for a benchmark set of 17 BODIPY/Aza-BODIPY molecules. However, these energies did not agree with values obtained from canonical CC2. Here we report LCC2* computations of vertical excitation energies on the same benchmark set of molecules using a newly implemented treatment of the ground state. Comparison with resolution-of-identity approximate coupled cluster to second-order (RI-CC2) results demonstrate that the new LCC2* results agree quantitatively. Furthermore, these values can easily be corrected empirically to also provide excellent agreement with the experiment. We show that the local algebraic diagrammatic construction to second-order (LADC(2)) method exhibits the same differences between implementations as seen for LCC2. The source of the difference is traced to an improved treatment of the ground state in the local methods, which decreases agreement with the experiment (as attributed to a fortuitous cancellation of errors) but significantly improves agreement with RI-CC2. While the absolute vertical excitation energies now show larger deviations, there remains a strong linear correlation between the LCC2* results and the experiment. For the 17 BODIPY/Aza-BODIPY molecules vertical excitation energies are determined using DLPNO-STEOM-CCSD and shown to have excellent agreement with experimental λmax (MAE = 0.145 eV), which is the best of all the single-reference methods. The vertical excitation energies are determined using LCC2*, empirically corrected LCC2*, and RI-CC2 for a series of eight large BODIPYs and Aza-BODIPYs.
Collapse
Affiliation(s)
- Milica Feldt
- Department of Chemistry, KU Leuven, Leuven, Belgium.,Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Münster, Germany
| | - Alex Brown
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
38
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. J Chem Theory Comput 2021; 17:902-926. [PMID: 33405921 DOI: 10.1021/acs.jctc.0c01129] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol-1, and this maximum error is reduced to below 0.2 kcal mol-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
39
|
Tew DP. Principal domains in F12 explicitly correlated theory. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
40
|
Bozkaya U, Ünal A, Alagöz Y. Energy and analytic gradients for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation. J Chem Phys 2020; 153:244115. [PMID: 33380091 DOI: 10.1063/5.0035811] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient implementations of the orbital-optimized coupled-cluster doubles (or simply "optimized CCD," OCCD, for short) method and its analytic energy gradients with the density-fitting (DF) approach, denoted by DF-OCCD, are presented. In addition to the DF approach, the Cholesky-decomposed variant (CD-OCCD) is also implemented for energy computations. The computational cost of the DF-OCCD method (available in a plugin version of the DFOCC module of PSI4) is compared with that of the conventional OCCD (from the Q-CHEM package). The OCCD computations were performed with the Q-CHEM package in which OCCD are denoted by OD. In the conventional OCCD method, one needs to perform four-index integral transformations at each of the CCD iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD provides dramatically lower computational costs compared to OCCD, and there are almost eightfold reductions in the computational time for the C6H14 molecule with the cc-pVTZ basis set. For open-shell geometries, interaction energies, and hydrogen transfer reactions, DF-OCCD provides significant improvements upon DF-CCD. Furthermore, the performance of the DF-OCCD method is substantially better for harmonic vibrational frequencies in the case of symmetry-breaking problems. Moreover, several factors make DF-OCCD more attractive compared to CCSD: (1) for DF-OCCD, there is no need for orbital relaxation contributions in analytic gradient computations; (2) active spaces can readily be incorporated into DF-OCCD; (3) DF-OCCD provides accurate vibrational frequencies when symmetry-breaking problems are observed; (4) in its response function, DF-OCCD avoids artificial poles; hence, excited-state molecular properties can be computed via linear response theory; and (5) symmetric and asymmetric triples corrections based on DF-OCCD [DF-OCCD(T)] have a significantly better performance in near degeneracy regions.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Yavuz Alagöz
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
41
|
Whitten JL. Estimates of electron correlation based on density expansions. J Chem Phys 2020; 153:244103. [PMID: 33380085 DOI: 10.1063/5.0031279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Two methods for estimating the correlation energy of molecules and other electronic systems are discussed based on the assumption that the correlation energy can be partitioned between atomic regions. In the first method, the electron density is expanded in terms of atomic contributions using rigorous electron repulsion bounds, and in the second method, correlation contributions are associated with basis function pairs. These methods do not consider the detailed nature of localized excitations but instead define a correlation energy per electron factor that is unique to a specific atom. The correlation factors are basis function dependent and are determined by configuration interaction (CI) calculations on diatomic and hydride molecules. The correlation energy estimates are compared with the results of high-level CI calculations for a test set of 27 molecules representing a wide range of bonding environments (average error of 2.6%). An extension based on truncated CI calculations in which d-type and hydrogen p-type functions are eliminated from the virtual space combined with estimates of dynamical correlation contributions using atomic correlation factors is discussed and applied to the dissociation of several molecules.
Collapse
Affiliation(s)
- Jerry L Whitten
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
42
|
Abstract
A broad range of approaches to many-body dispersion are discussed, including empirical approaches with multiple fitted parameters, augmented density functional-based approaches, symmetry adapted perturbation theory, and a supermolecule approach based on coupled cluster theory. Differing definitions of "body" are considered, specifically atom-based vs molecule-based approaches.
Collapse
Affiliation(s)
- Peng Xu
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| | - Melisa Alkan
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| |
Collapse
|
43
|
Sorathia K, Tew DP. Basis set extrapolation in pair natural orbital theories. J Chem Phys 2020; 153:174112. [PMID: 33167642 DOI: 10.1063/5.0022077] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.
Collapse
Affiliation(s)
- Kesha Sorathia
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, GermanyUniversity of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
44
|
von Rudorff GF, Heinen SN, Bragato M, von Lilienfeld OA. Thousands of reactants and transition states for competing E2 and S$_\mathrm{N}$2 reactions. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/aba822] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
45
|
Sidorczuk D, Kozanecki M, Civalleri B, Pernal K, Prywer J. Structural and Optical Properties of Struvite. Elucidating Structure of Infrared Spectrum in High Frequency Range. J Phys Chem A 2020; 124:8668-8678. [PMID: 32972131 PMCID: PMC7586396 DOI: 10.1021/acs.jpca.0c04707] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Revised: 09/22/2020] [Indexed: 12/02/2022]
Abstract
Study of structure and optical properties of magnesium ammonium phosphate hexahydrate crystal known as struvite is presented. Experimentally determined infrared (IR) and ultraviolet-visible (UV-vis) spectra are compared with the theoretical predictions of density functional methods. Examination of the interatomic bond lengths, Mulliken atomic charges, and binding energies of water in the magnesium hexahydrate cation, together with the analysis of the hydrogen bond pattern have allowed us to explain a special feature of the IR spectrum of struvite, a blueshift of the band corresponding to the O-H stretching mode. This mode has been assigned to a "dangling" hydroxyl group in one of the water molecules in magnesium hexahydrate. Using experimentally obtained UV-vis spectrum and performing Tauc plots analysis, optical bandgap of struvite has been narrowed to a range from 5.92 to 6.06 eV.
Collapse
Affiliation(s)
- Dominik Sidorczuk
- Institute
of Physics, Lodz University of Technology, ul. Wólczańska 219, 90-924Łódź, Poland
| | - Marcin Kozanecki
- Department
of Molecular Physics, Lodz University of
Technology, ul. Żeromskiego 116, 90-924 Łódź, Poland
| | - Bartolomeo Civalleri
- Department
of Chemistry, University of Torino, Via P. Giuria 7-10125 Torino, Italy
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wólczańska 219, 90-924Łódź, Poland
| | - Jolanta Prywer
- Institute
of Physics, Lodz University of Technology, ul. Wólczańska 219, 90-924Łódź, Poland
| |
Collapse
|
46
|
Mallick S, Roy B, Kumar P. A comparison of DLPNO-CCSD(T) and CCSD(T) method for the determination of the energetics of hydrogen atom transfer reactions. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112934] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
47
|
Wang X, Lewis CA, Valeev EF. Efficient evaluation of exact exchange for periodic systems via concentric atomic density fitting. J Chem Phys 2020; 153:124116. [DOI: 10.1063/5.0016856] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiao Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Cannada A. Lewis
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
48
|
Ye HZ, Tran HK, Van Voorhis T. Bootstrap Embedding For Large Molecular Systems. J Chem Theory Comput 2020; 16:5035-5046. [PMID: 32589842 DOI: 10.1021/acs.jctc.0c00438] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N3) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
49
|
Neese F, Wennmohs F, Becker U, Riplinger C. The ORCA quantum chemistry program package. J Chem Phys 2020; 152:224108. [DOI: 10.1063/5.0004608] [Citation(s) in RCA: 697] [Impact Index Per Article: 174.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
- FAccTs GmbH, Rolandstr. 67, 50677 Köln, Germany
| | - Frank Wennmohs
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Ute Becker
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | |
Collapse
|
50
|
Lang J, Antalík A, Veis L, Brandejs J, Brabec J, Legeza Ö, Pittner J. Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T). J Chem Theory Comput 2020; 16:3028-3040. [PMID: 32275424 DOI: 10.1021/acs.jctc.0c00065] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Sciences, Charles University, Albertov 6, 128 00 Prague 2, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|