1
|
Jiang A, Glick ZL, Poole D, Turney JM, Sherrill CD, Schaefer HF. Accurate and efficient open-source implementation of domain-based local pair natural orbital (DLPNO) coupled-cluster theory using a t1-transformed Hamiltonian. J Chem Phys 2024; 161:082502. [PMID: 39171707 DOI: 10.1063/5.0219963] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2024] [Accepted: 08/05/2024] [Indexed: 08/23/2024] Open
Abstract
We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol-1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol-1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms).
Collapse
Affiliation(s)
- Andy Jiang
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - David Poole
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Justin M Turney
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
2
|
Aliverti-Piuri D, Chatterjee K, Ding L, Liao K, Liebert J, Schilling C. What can quantum information theory offer to quantum chemistry? Faraday Discuss 2024. [PMID: 39099471 DOI: 10.1039/d4fd00059e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
It is the ultimate goal of this work to foster synergy between quantum chemistry and the flourishing field of quantum information theory. For this, we first translate quantum information concepts, such as entanglement and correlation, into the context of quantum chemical systems. In particular, we establish two conceptually distinct perspectives on 'electron correlation', leading to a notion of orbital and particle correlation. We then demonstrate that particle correlation equals total orbital correlation minimized over all orbital bases. Accordingly, particle correlation resembles the minimal, thus intrinsic, complexity of many-electron wave functions, while orbital correlation quantifies their complexity relative to a basis. We illustrate these concepts of intrinsic and extrinsic correlation complexity in molecular systems, which also manifests the crucial link between the two correlation pictures. Our results provide theoretical justification for the long-favored natural orbitals for simplifying electronic structures, and open new pathways for developing more efficient approaches towards the electron correlation problem.
Collapse
Affiliation(s)
- Damiano Aliverti-Piuri
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Kaustav Chatterjee
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Lexin Ding
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Ke Liao
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Julia Liebert
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Christian Schilling
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| |
Collapse
|
3
|
Ye HZ, Berkelbach TC. Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory. Faraday Discuss 2024. [PMID: 39049598 DOI: 10.1039/d4fd00041b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, NY 10027, USA.
- Initiative for Computational Catalysis, Flatiron Institute, New York, NY 10010, USA
| |
Collapse
|
4
|
Wei Y, Debnath S, Weber JL, Mahajan A, Reichman DR, Friesner RA. Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry. J Phys Chem A 2024; 128:5796-5807. [PMID: 38970826 DOI: 10.1021/acs.jpca.4c02853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/08/2024]
Abstract
This study evaluates the precision of widely recognized quantum chemical methodologies, CCSD(T), DLPNO-CCSD(T), and localized ph-AFQMC, for determining the thermochemistry of main group elements. DLPNO-CCSD(T) and localized ph-AFQMC, which offer greater scalability compared to canonical CCSD(T), have emerged over the past decade as pivotal in producing precise benchmark chemical data. Our investigation includes closed-shell, neutral molecules, focusing on their heat of formation and atomization energy sourced from four specific small molecule data sets. First, we selected molecules from the G2 and G3 data sets, noted for their reliable experimental heat of formation data. Additionally, we incorporate molecules from the W4-11 and W4-17 sets, which provide high-level theoretical reference values for atomization energy at 0 K. Our findings reveal that both DLPNO-CCSD(T) and ph-AFQMC methods are capable of achieving a root-mean-square deviation of less than 1 kcal/mol across the combined data set, aligning with the threshold for chemical accuracy. Moreover, we make efforts to confine the maximum deviations within 2 kcal/mol, a degree of precision that significantly broadens the applicability of these methods in fields such as biology and materials science.
Collapse
Affiliation(s)
- Yujing Wei
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Sibali Debnath
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Ankit Mahajan
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
5
|
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
|
6
|
Liao K, Ding L, Schilling C. Quantum Information Orbitals (QIO): Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality. J Phys Chem Lett 2024; 15:6782-6790. [PMID: 38913404 DOI: 10.1021/acs.jpclett.4c01105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C2 and Cr2 molecule.
Collapse
Affiliation(s)
- Ke Liao
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Lexin Ding
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Christian Schilling
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| |
Collapse
|
7
|
Ibrahim A, Roy PN. A neural network-based four-body potential energy surface for parahydrogen. J Chem Phys 2024; 160:244308. [PMID: 38916269 DOI: 10.1063/5.0214495] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Accepted: 06/06/2024] [Indexed: 06/26/2024] Open
Abstract
We present an isotropic ab initio (para-H2)4 four-body interaction potential energy surface (PES). The electronic structure calculations are performed at the correlated coupled-cluster theory level, with single, double, and perturbative triple excitations. They use an atom-centered augmented correlation-consistent double zeta basis set, supplemented by a (3s3p2d) midbond function. We use a multilayer perceptron to construct the PES. We apply a rescaling transformation to the output energies during training to improve the prediction of weaker energies in the sample data. At long distances, the interaction energies are adjusted to match the empirically derived four-body dispersion interaction. The four-body interaction energy at short intermolecular separations is net repulsive. The use of this four-body PES, in combination with a first principles pair potential for para-H2 [J. Chem. Phys. 119, 12551 (2015)] and an isotropic ab initio three-body potential for para-H2 [J. Chem. Phys. 156, 044301 (2022)], is expected to provide closer agreement with experimental results.
Collapse
Affiliation(s)
- Alexander Ibrahim
- Department of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
8
|
Türkmen I, Dolg M. Linear Scaling Incremental Scheme for Correlation Energies with Embedding Generated Virtuals. J Chem Theory Comput 2024; 20:3154-3168. [PMID: 38588492 DOI: 10.1021/acs.jctc.3c01377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
A novel incremental scheme is presented including an incremental expansion of the virtual space for the calculation of electron correlation energies, which is compatible with any size-extensive correlation method and scales asymptotically linear for large molecules. The performance is studied for organic molecules, water clusters, and a La(III)-water complex, where the compatibility with pseudopotentials is also examined. The computational requirements are already reduced tremendously for medium-sized water clusters and hydrocarbons with respect to the canonical CCSD as well as the ordinary incremental scheme references. Correlation energies within chemical accuracy have been observed for all studied systems. The novelty of the method is that relatively small virtual spaces are used in combination with tuples of localized occupied spaces. The corresponding orthonormal occupied and virtual orbitals are obtained from QM/QM embedding calculations and can thus be used with standard quantum chemistry codes for correlation calculations. It is presented how relevant virtual spaces are selected and the correlation energies are linked in the new virtual space expansion.
Collapse
Affiliation(s)
- Ilyas Türkmen
- Institute for Theoretical Chemistry, University of Cologne, Greinstr. 4, D-50939 Cologne, Germany
| | - Michael Dolg
- Institute for Theoretical Chemistry, University of Cologne, Greinstr. 4, D-50939 Cologne, Germany
| |
Collapse
|
9
|
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
|
10
|
Kurian JS, Ye HZ, Mahajan A, Berkelbach TC, Sharma S. Toward Linear Scaling Auxiliary-Field Quantum Monte Carlo with Local Natural Orbitals. J Chem Theory Comput 2024; 20:134-142. [PMID: 38113195 DOI: 10.1021/acs.jctc.3c01122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
We develop a local correlation variant of auxiliary-field quantum Monte Carlo (AFQMC) based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with the system size for a target accuracy, the method has linear scaling. Applying LNO-AFQMC to molecular problems containing a few hundred to a thousand orbitals, we demonstrate convergence of total energies with significantly reduced costs. The savings are more significant for larger systems and larger basis sets. However, even for our smallest system studied, we find that LNO-AFQMC is cheaper than canonical AFQMC, in contrast with many other reduced-scaling methods. Perhaps most significantly, we show that energy differences converge much more quickly than total energies, making the method ideal for applications in chemistry and material science. Our work paves the way for linear scaling AFQMC calculations of strongly correlated systems, which would have a transformative effect on ab initio quantum chemistry.
Collapse
Affiliation(s)
- Jo S Kurian
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, United States
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Ankit Mahajan
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, United States
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, United States
| |
Collapse
|
11
|
Li W, Wang Y, Ni Z, Li S. Cluster-in-Molecule Local Correlation Method for Dispersion Interactions in Large Systems and Periodic Systems. Acc Chem Res 2023; 56:3462-3474. [PMID: 37991873 DOI: 10.1021/acs.accounts.3c00538] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2023]
Abstract
ConspectusThe noncovalent interactions, including dispersion interactions, control the structures and stabilities of complex chemical systems, including host-guest complexes and the adsorption process of molecules on the solid surfaces. The density functional theory (DFT) with empirical dispersion correction is now the working horse in many areas of applications. Post-Hartree-Fock (post-HF) methods have been well recognized to provide more accurate descriptions in a systematic way. However, traditional post-HF methods are mainly limited to small- or medium-sized systems, and their applications to periodic condensed phase systems are still very limited due to their expensive computational costs.To extend post-HF calculations to large molecules, the cluster-in-molecule (CIM) local correlation approach has been established, allowing highly accurate electron correlation calculations that are routinely available for very large systems. In the CIM approach, the electron correlation energy of a large molecule could be obtained from electron correlation calculations on a series of clusters, each of which contains a subset of occupied and virtual localized molecular orbitals. The CIM method could be massively and efficiently parallelized on general computer clusters. The CIM method has been implemented at various electron correlation levels, including second-order Mo̷ller-Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), CCSD with perturbative triples correction [CCSD(T)], etc. The CIM-MP2 energy gradient algorithm was developed and applied to the geometry optimizations of large systems. The CIM method has also been extended to condensed-phase systems under periodic boundary conditions (PBC-CIM). For periodic systems, the correlation energy per unit cell could be evaluated with correlation energy contributions from a series of clusters that are built with localized Wannier functions.CIM-based electron correlation calculations have been employed to investigate a number of chemical problems in which the dispersion interaction is important. CIM-based post-HF methods including CIM domain-based local pair natural orbital (DLPNO) CCSD(T) are applied to compute the relative or binding energies of biological systems or supramolecular complexes, the reaction barrier in a relatively complex chemical reaction. The CIM-MP2 method is used to obtain the optimized geometry of large systems. CIM-based post-HF calculations have also been used to compute the cohesive energies of molecular crystals and adsorption energies of molecules on the solid surfaces. The CIM and its PBC variant are expected to become a powerful theoretical tool for accurate calculations of the energies and structures for a broad range of large systems and condensed-phase systems with significant dispersion interactions.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of 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 210023, People's Republic of China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Hangzhou Normal University, Hangzhou 311121, People's Republic of 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 210023, People's Republic of China
| |
Collapse
|
12
|
Wang W, Whitfield JD. Basis Set Generation and Optimization in the NISQ Era with Quiqbox.jl. J Chem Theory Comput 2023; 19:8032-8052. [PMID: 37924295 DOI: 10.1021/acs.jctc.3c00011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2023]
Abstract
In the noisy intermediate-scale quantum era, ab initio computation of electronic structure problems has become one of the major benchmarks for identifying the boundary between classical and quantum computational power. Basis sets play a key role in the electronic structure methods implemented on both classical and quantum devices. To investigate the consequences of single-particle basis sets, we propose a framework for more customizable basis set generation and optimization. This framework allows composite basis sets to go beyond typical basis set frameworks, such as atomic basis sets, by introducing the concept of mixed-contracted Gaussian-type orbitals. These basis set generations set the stage for more flexible variational optimization of basis set parameters. To realize this framework, we have developed an open-source software package named "Quiqbox" in the Julia programming language. We demonstrate various examples of using Quiqbox for basis set optimization and generation, ranging from optimizing atomic basis sets on the Hartree-Fock level, preparing the initial state for variational quantum eigensolver computation, and constructing basis sets with completely delocalized orbitals. We also include various benchmarks of Quiqbox for basis set optimization and ab initial electronic structure computation.
Collapse
Affiliation(s)
- Weishi Wang
- Department of Physics and Astronomy, Dartmouth College, Hanover, New Hampshire 03755, United States
| | - James D Whitfield
- Department of Physics and Astronomy, Dartmouth College, Hanover, New Hampshire 03755, United States
- AWS Center for Quantum Computing, Pasadena, California 91106, United States
| |
Collapse
|
13
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Local Natural Orbital CCSD(T) Approach for Open-Shell Systems: Algorithms, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2023; 19:8166-8188. [PMID: 37921429 PMCID: PMC10687875 DOI: 10.1021/acs.jctc.3c00881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 10/05/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023]
Abstract
The extension of the highly optimized local natural orbital (LNO) coupled cluster (CC) with single-, double-, and perturbative triple excitations [LNO-CCSD(T)] method is presented for high-spin open-shell molecules based on restricted open-shell references. The techniques enabling the outstanding efficiency of the closed-shell LNO-CCSD(T) variant are adopted, including the iteration- and redundancy-free second-order Møller-Plesset and (T) formulations as well as the integral-direct, memory- and disk use-economic, and OpenMP-parallel algorithms. For large molecules, the efficiency of our open-shell LNO-CCSD(T) method approaches that of its closed-shell parent method due to the application of restricted orbital sets for demanding integral transformations and a novel approximation for higher-order long-range spin-polarization effects. The accuracy of open-shell LNO-CCSD(T) is extensively tested for radicals and reactions thereof, ionization processes, as well as spin-state splittings, and transition-metal compounds. At the size range where the canonical CCSD(T) reference is accessible (up to 20-30 atoms), the average open-shell LNO-CCSD(T) correlation energies are found to be 99.9 to 99.95% accurate, which translates into average absolute deviations of a few tenths of kcal/mol in the investigated energy differences already with the default settings. For more extensive molecules, the local errors may grow, but they can be estimated and decreased via affordable systematic convergence studies. This enables the accurate modeling of large systems with complex electronic structures, as illustrated on open-shell organic radicals and transition-metal complexes of up to 179 atoms as well as on challenging biochemical systems, including up to 601 atoms and 11,000 basis functions. While the protein models involve difficulties for local approximations, such as the spin states of a bounded iron ion or an extremely delocalized singly occupied orbital, the corresponding single-node LNO-CCSD(T) computations were feasible in a matter of days with 10s to 100 GB of memory use. Therefore, the new LNO-CCSD(T) implementation enables highly accurate computations for open-shell systems of unprecedented size and complexity with widely accessible hardware.
Collapse
Affiliation(s)
- P. Bernát Szabó
- 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
| | - József Csóka
- 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
- HUN-REN-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
| | - 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
- HUN-REN-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
- HUN-REN-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
|
14
|
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
|
15
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
16
|
Semidalas E, Martin JML. Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals. J Chem Theory Comput 2023; 19:5806-5820. [PMID: 37540641 PMCID: PMC10500978 DOI: 10.1021/acs.jctc.3c00506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Indexed: 08/06/2023]
Abstract
We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VnZ(-PP)-F12-wis (n = D,T), for the d-block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3d-block, while its ECP counterpart for the 4d and 5d-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [9s8p6d3f2g] for the 3d-block elements and [7s7p6d3f2g] for the 4d and 5d-block elements. Our VnZ(-PP)-F12-wis basis sets are evaluated on challenging test sets for metal-organic barrier heights (MOBH35) and group-11 metal clusters (CUAGAU-2). In F12 calculations, they are found to be about as close to the complete basis set limit as the combination of standard cc-pVnZ-F12 on main-group elements with the standard aug-cc-pV(n+1)Z(-PP) basis sets on the transition metal(s). While our basis sets are somewhat more compact than aug-cc-pV(n+1)Z(-PP), the CPU time benefit is negligible for catalytic complexes that contain only one or two transition metals among dozens of main-group elements; however, it is somewhat more significant for metal clusters.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
17
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1256510. [PMID: 37654900 PMCID: PMC10466216 DOI: 10.3389/fchem.2023.1256510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 07/11/2023] [Indexed: 09/02/2023] Open
Abstract
[This corrects the article DOI: 10.3389/fchem.2023.1154526.].
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
18
|
Vysotskiy VP, Torbjörnsson M, Jiang H, Larsson ED, Cao L, Ryde U, Zhai H, Lee S, Chan GKL. Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art ab initio methods. J Chem Phys 2023; 159:044106. [PMID: 37486046 DOI: 10.1063/5.0152611] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 07/03/2023] [Indexed: 07/25/2023] Open
Abstract
We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.
Collapse
Affiliation(s)
- Victor P Vysotskiy
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Magne Torbjörnsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ernst D Larsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Lili Cao
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
19
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1154526. [PMID: 37388945 PMCID: PMC10303140 DOI: 10.3389/fchem.2023.1154526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 05/11/2023] [Indexed: 07/01/2023] Open
Abstract
This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
20
|
Monzel L, Holzer C, Klopper W. Natural virtual orbitals for the GW method in the random-phase approximation and beyond. J Chem Phys 2023; 158:144102. [PMID: 37061489 DOI: 10.1063/5.0144469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G0W0 and evGW0 as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.
Collapse
Affiliation(s)
- Laurenz Monzel
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| |
Collapse
|
21
|
Demel O, Lecours MJ, Nooijen M. Further investigations into a Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals: Multipole correction, OSV extrapolation, and critical assessment. J Chem Phys 2023; 158:114120. [PMID: 36948803 DOI: 10.1063/5.0135113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
We report further investigations to aid the development of a Laplace MP2 (second-order Møller Plesset) method with a range separated Coulomb potential partitioned into short- and long-range parts. The implementation of the method extensively uses sparse matrix algebra, density fitting techniques for the short-range part, and a Fourier transformation in spherical coordinates for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas virtual space is described by orbital specific virtual orbitals (OSVs) associated with localized molecular orbitals. The Fourier transform is deficient for very large distances between localized occupied orbitals, and a multipole expansion for widely separated pairs is introduced for the direct MP2 contribution, which is applicable also to non-Coulombic potentials that do not satisfy the Laplace equation. For the exchange contribution, an efficient screening of contributing localized occupied pairs is employed, which is discussed more completely here. To mitigate errors due to the truncation of OSVs, a simple and efficient extrapolation procedure is used to obtain results close to MP2 for the full basis set of atomic orbitals Using a suitable set of default parameters, the accuracy of the approach is demonstrated. The current implementation of the approach is not very efficient, and the aim of this paper is to introduce and critically discuss ideas that can have more general applicability beyond MP2 calculations for large molecules.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
22
|
Jiang A, Turney JM, Schaefer HF. Tensor Hypercontraction Form of the Perturbative Triples Energy in Coupled-Cluster Theory. J Chem Theory Comput 2023; 19:1476-1486. [PMID: 36802552 PMCID: PMC10018738 DOI: 10.1021/acs.jctc.2c00996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Abstract
We present the working equations for a reduced-scaling method of evaluating the perturbative triples (T) energy in coupled-cluster theory, through the tensor hypercontraction (THC) of the triples amplitudes (tijkabc). Through our method, we can reduce the scaling of the (T) energy from the traditional O(N7) to a more modest O(N5). We also discuss implementation details to aid future research, development, and software realization of this method. Additionally, we show that this method yields submillihartree (mEh) differences from CCSD(T) when evaluating absolute energies and sub-0.1 kcal/mol energy differences when evaluating relative energies. Finally, we demonstrate that this method converges to the true CCSD(T) energy through the systematic increasing of the rank or eigenvalue tolerance of the orthogonal projector, as well as exhibiting sublinear to linear error growth with respect to system size.
Collapse
Affiliation(s)
- Andy Jiang
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Justin M Turney
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Henry F Schaefer
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
23
|
Naim C, Besalú-Sala P, Zaleśny R, Luis JM, Castet F, Matito E. Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties? J Chem Theory Comput 2023; 19:1753-1764. [PMID: 36862983 DOI: 10.1021/acs.jctc.2c01212] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/04/2023]
Abstract
Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for β and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.
Collapse
Affiliation(s)
- Carmelo Naim
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France.,Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, 20080 Donostia, Euskadi, Spain
| | - Pau Besalú-Sala
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia, Spain
| | - Robert Zaleśny
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Josep M Luis
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia, Spain
| | - Frédéric Castet
- Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France
| | - Eduard Matito
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Ikerbasque Foundation for Science, 48011 Bilbao, Euskadi, Spain
| |
Collapse
|
24
|
Khire SS, Gadre SR. Development and testing of an algorithm for efficient MP2/CCSD(T) energy estimation of molecular clusters with the 2-body approach. J Comput Chem 2023; 44:261-267. [PMID: 35514315 DOI: 10.1002/jcc.26881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2022] [Revised: 04/08/2022] [Accepted: 04/13/2022] [Indexed: 12/31/2022]
Abstract
This work reports the development and testing of an automated algorithm for estimating the energies of weakly bound molecular clusters employing correlated theory. Firstly, the monomers and dimers of (homo/hetero) clusters are identified, and the sum of one-body and two-body contributions to correlation energy is calculated. The addition of this contribution to the Hartree-Fock full calculation (FC) energies provides a good estimate of the total energies at Møller-Plesset second-order perturbation theory (MP2)/coupled-cluster method with singles and doubles (CCSD) (T)-level theory using augmented Dunning basis sets. The estimated energies for several test clusters show an excellent agreement with their FC counterparts, with a substantial wall-clock time saving employing off-the-shelf hardware. Furthermore, the complete basis set (CBS) limit for MP2 energy computed using the two-body approach also agrees with its CBS energy with its FC counterpart.
Collapse
Affiliation(s)
- Subodh S Khire
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune, India
| |
Collapse
|
25
|
Győri T, Czakó G. A comprehensive benchmark ab initio survey of the stationary points and products of the OH· + CH 3OH system. J Chem Phys 2023; 158:034301. [PMID: 36681627 DOI: 10.1063/5.0133978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Reactions between methanol and the hydroxyl radical are of significant interest for combustion-, atmospheric-, and astrochemistry. While the two primary product channels (the formation of H2O with either CH3O· or ·CH2OH) have been the subject of numerous studies, the possibility of other products has seen little attention. Here, we present a comprehensive thermochemical survey of the stationary points and plausible products of the reaction, featuring 29 geometries optimized at the UCCSD(T)-F12b/aug-cc-pVTZ level, followed by accurate composite ab initio computations for all stationary points (including ·CH2OH dissociation and isomerization) and five product channels, with a detailed evaluation of basis set convergence and efficiency. The computations reveal that the formation of methanediol and the hydroxymethoxy radical is thermodynamically favorable and the endothermicity of formaldehyde formation is low enough to be a plausible product channel. We also observe unexpectedly large energy deviations between the partially-spin-adapted ROHF-RCCSD(T) method and ROHF-UCCSD(T) as well as between UHF-UCCSDT(Q) and ROHF-UCCSDT(Q) results.
Collapse
Affiliation(s)
- Tibor Győri
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
26
|
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
|
27
|
Zhou X, Cao Z, Wang F, Wang Z. Barrier heights, reaction energies and bond dissociation energies for RH + HO 2 reactions with coupled-cluster theory, density functional theory and diffusion quantum Monte Carlo methods. Phys Chem Chem Phys 2022; 25:341-350. [PMID: 36477176 DOI: 10.1039/d2cp04463c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Hydrogen abstraction reactions by the HO2 radical from hydrocarbon molecules are an important class of reactions in the autoignition of hydrocarbon fuels. Performance of DLPNO-CC and DFT methods using three hybrids and four double hybrids as well as FN-DMC with the single-Slater-Jastrow trial wavefunction on barrier heights and reaction energies of RH + HO2 reactions as well as bond dissociation energies of the involved X-H molecules is evaluated by comparison with the highly accurate CCSD(T)-F12b/CBS results in this study. Our results show that the DLPNO-CCSD(T)-F12 method can achieve highly accurate barrier heights, reaction energies and X-H bond energies for RH + HO2 reactions at a relatively low computational cost, and it is applicable to the H-abstraction reactions of larger molecules. Among all DFAs, MN15 and the employed double hybrids can achieve accurate barrier heights and reaction energies with MADs of less than or around 2 kJ mol-1, but their error on X-H bond energies is more pronounced. Only DSD-BLYP and DSD-PBEB95 can provide X-H bond energies with MADs less than 4 kJ mol-1. Considering dispersion correction in DFT calculations does not improve these barrier heights and reaction energies. The error of FN-DMC on barrier heights and reaction energies is slightly larger than that of MN15 and those of double hybrids, but it can achieve results within chemical accuracy for these reactions and the X-H bond energies.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Department of Physics, Shaanxi University of Science and Technology, Xi'an, 710021, P. R. China.
| | - Zhanli Cao
- School of Science, Xi'an University of Posts and Telecommunications, Xi'an, 710121, P. R. China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, P. R. China
| | - Zhifan Wang
- School of Electronic Engineering, Chengdu Technological University, Chengdu, P. R. China
| |
Collapse
|
28
|
Karton A. Tightening the Screws: The Importance of Tight d Functions in Coupled-Cluster Calculations up to the CCSDT(Q) Level. J Phys Chem A 2022; 126:8544-8555. [DOI: 10.1021/acs.jpca.2c06522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Amir Karton
- School of Science and Technology, University of New England, Armidale, New South Wales2351, Australia
| |
Collapse
|
29
|
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
|
30
|
Demir Gİ, Demir S, Tekin A. 2D‐FFCASP—A New Approach for 2D Structure Prediction Applied to Self‐Assemblies of DNA Bases. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Gözde İniş Demir
- Informatics Institute Istanbul Technical University Maslak Istanbul 34469 Turkey
| | - Samet Demir
- Informatics Institute Istanbul Technical University Maslak Istanbul 34469 Turkey
- TÜBİTAK Research Institute for Fundamental Sciences Gebze Kocaeli 41470 Turkey
| | - Adem Tekin
- Informatics Institute Istanbul Technical University Maslak Istanbul 34469 Turkey
- TÜBİTAK Research Institute for Fundamental Sciences Gebze Kocaeli 41470 Turkey
| |
Collapse
|
31
|
Bozkaya U, Ermiş B. Linear-Scaling Systematic Molecular Fragmentation Approach for Perturbation Theory and Coupled-Cluster Methods. J Chem Theory Comput 2022; 18:5349-5359. [PMID: 35972734 PMCID: PMC9476663 DOI: 10.1021/acs.jctc.2c00587] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The coupled-cluster (CC) singles and doubles with perturbative
triples [CCSD(T)] method is frequently referred to as the “gold
standard” of modern computational chemistry. However, the high
computational cost of CCSD(T) [O(N7)], where N is the number of basis functions,
limits its applications to small-sized chemical systems. To address
this problem, efficient implementations of linear-scaling coupled-cluster
methods, which employ the systematic molecular fragmentation (SMF)
approach, are reported. In this study, we aim to do the following:
(1) To achieve exact linear scaling and to obtain a pure ab
initio approach, we revise the handling of nonbonded interactions
in the SMF approach, denoted by LSSMF. (2) A new fragmentation algorithm,
which yields smaller-sized fragments, that better fits high-level
CC methods is introduced. (3) A modified nonbonded fragmentation scheme
is proposed to enhance the existent algorithm. Performances of the
LSSMF-CC approaches, such as LSSMF-CCSD(T), are compared with their
canonical versions for a set of alkane molecules, CnH2n+2 (n = 6–10),
which includes 142 molecules. Our results demonstrate that the LSSMF
approach introduces negligible errors compared with the canonical
methods; mean absolute errors (MAEs) are between 0.20 and 0.59 kcal
mol–1 for LSSMF(3,1)-CCSD(T). For a larger alkanes
set (L12), CnH2n+2 (n = 50–70), the performance of
LSSMF for the second-order perturbation theory (MP2) is investigated.
For the L12 set, various bonded and nonbonded levels are considered.
Our results demonstrate that the combination of bonded level 6 with
nonbonded level 2, LSSMF(6,2), provides very accurate results for
the MP2 method with a MAE value of 0.32 kcal mol–1. The LSSMF(6,2) approach yields more than a 26-fold reduction in
errors compared with LSSMF(3,1). Hence, we obtain substantial improvements
over the original SMF approach. To illustrate the efficiency and applicability
of the LSSMF-CCSD(T) approach, we consider an alkane molecule with
10,004 atoms. For this molecule, the LSSMF(3,1)-CCSD(T)/cc-pVTZ energy
computation, on a Linux cluster with 100 nodes, 4 cores, and 5 GB
of memory provided to each node, is performed just in ∼24 h.
As a second test, we consider a biomolecular complex (PDB code: 1GLA), which includes
10,488 atoms, to assess the efficiency of the LSSMF approach. The
LSSMF(3,1)-FNO–CCSD(T)/cc-pVTZ energy computation is completed
in ∼7 days for the biomolecular complex. Hence, our results
demonstrate that the LSSMF-CC approaches are very efficient. Overall,
we conclude the following: (1) The LSSMF(m, n)-CCSD(T) methods can be reliably used for large-scale
chemical systems, where the canonical methods are not computationally
affordable. (2) The accuracy of bonded level 3 is not satisfactory
for large chemical systems. (3) For high-accuracy studies, bonded
level 5 (or higher) and nonbonded level 2 should be employed.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Betül Ermiş
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
32
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
33
|
Barcza B, Szirmai ÁB, Szántó KJ, Tajti A, Szalay PG. Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles. J Comput Chem 2022; 43:1079-1093. [PMID: 35478353 PMCID: PMC9321956 DOI: 10.1002/jcc.26866] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/26/2022] [Accepted: 03/29/2022] [Indexed: 01/15/2023]
Abstract
The ground state intermolecular potential of bimolecular complexes of N-heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short- and long-distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine-uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.
Collapse
Affiliation(s)
- Bónis Barcza
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Ádám B. Szirmai
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Katalin J. Szántó
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Péter G. Szalay
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| |
Collapse
|
34
|
Peng L, Peng D, Gu FL, Yang W. Regularized Localized Molecular Orbitals in a Divide-and-Conquer Approach for Linear Scaling Calculations. J Chem Theory Comput 2022; 18:2975-2982. [PMID: 35416665 PMCID: PMC9972215 DOI: 10.1021/acs.jctc.2c00142] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Non-orthogonal localized molecular orbitals (NOLMOs) have been employed as building blocks for the divide-and-conquer (DC) linear scaling method. The NOLMOs are calculated from subsystems and used for constructing the density matrix (DM) of the entire system, instead of the subsystem DM in the original DC approach. Also, unlike the original DC method, the inverse electronic temperature parameter β is not needed anymore. Furthermore, a new regularized localization approach for NOLMOs has been developed, in which the localization cost function is a sum of the spatial spread function, as in the Boys method, and the kinetic energy, as a regularization measure to limit the oscillation of the NOLMOs. The optimal weight of the kinetic energy can be determined by optimization with analytical gradients. The resulting regularized NOLMOs have enhanced smoothness and better transferability because of reduced kinetic energies. Compared with the original DC, while NOLMO-DC has a similar computational linear scaling cost, the accuracy of NOLMO-DC is better by several orders of magnitude for large conjugated systems and by about 1 order of magnitude for other systems. The NOLMO-DC method is thus a promising development of the DC approach for linear scaling calculations.
Collapse
Affiliation(s)
- Liang Peng
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Daoling Peng
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708-0346, United States
| |
Collapse
|
35
|
Vitale E, Li Manni G, Alavi A, Kats D. FCIQMC-Tailored Distinguishable Cluster Approach: Open-Shell Systems. J Chem Theory Comput 2022; 18:3427-3437. [PMID: 35522217 PMCID: PMC9202306 DOI: 10.1021/acs.jctc.2c00059] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A recently proposed
tailored approach based on the distinguishable
cluster method and the stochastic FCI solver, FCIQMC [J. Chem.
Theory Comput. 2020, 16, 5621], is extended to open-shell
molecular systems. The method is employed to calculate spin gaps of
various Fe(II) complexes, including a Fe(II) porphyrin model system.
Both distinguishable cluster and fully relaxed CASSCF natural orbitals
were used in this work as reference for the subsequent tailored distinguishable
cluster calculations. The distinguishable cluster natural orbitals
occupation numbers were also used as an aid to the selection of the
active space. The effect of the active space sizes and of the explicit
correlation correction (F12) onto the predicted spin gaps is investigated.
The tailored distinguishable cluster with singles and doubles yields
consistently more accurate results compared to the tailored coupled
cluster with singles and doubles.
Collapse
Affiliation(s)
- Eugenio Vitale
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Giovanni Li Manni
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
36
|
Karton A, Thimmakondu VS. From Molecules with a Planar Tetracoordinate Carbon to an Astronomically Known C 5H 2 Carbene. J Phys Chem A 2022; 126:2561-2568. [PMID: 35426667 PMCID: PMC9442649 DOI: 10.1021/acs.jpca.2c01261] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Ethynylcyclopropenylidene
(2), an isomer of C5H2, is a known
molecule in the laboratory and has
recently been identified in Taurus Molecular Cloud-1 (TMC-1). Using
high-level coupled-cluster methods up to the CCSDT(Q)/CBS level of
theory, it is shown that two isomers of C5H2 with a planar tetracoordinate carbon (ptC) atom, (SP-4)-spiro[2.2]pent-1,4-dien-1,4-diyl
(11) and (SP-4)-spiro[2.2]pent-1,4-dien-1,5-diyl (13), serve as the reactive intermediates for the formation
of 2. Here, a theoretical connection has been established
between molecules containing ptC atoms (11 and 13) and a molecule (2) that is present nearly
430 light years away, thus providing evidence for the existence of
ptC species in the interstellar medium. The reaction pathways connecting
the transition states and the reactants and products have been confirmed
by intrinsic reaction coordinate calculations at the CCSDT(Q)/CBS//B3LYP-D3BJ/cc-pVTZ
level. While isomer 11 is non-polar (μ = 0), isomers 2 and 13 are polar, with dipole moment values
of 3.52 and 5.17 Debye at the CCSD(T)/cc-pVTZ level. Therefore, 13 is also a suitable candidate for both laboratory and radioastronomical
studies.
Collapse
Affiliation(s)
- Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, Western Australia 6009, Australia
| | - Venkatesan S Thimmakondu
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182-1030, USA
| |
Collapse
|
37
|
Nguyen TL, Perera A. Reaction of Methylidyne with Ethane: The C-C Insertion Is Unimportant. J Phys Chem A 2022; 126:1966-1972. [PMID: 35302775 DOI: 10.1021/acs.jpca.2c00735] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-accuracy coupled-cluster calculations in combination with the E,J-resolved master-equation analysis are used to study the reaction mechanism and kinetics of methylidyne with ethane. This reaction plays an important role in the combustion of hydrocarbon fuels and in interstellar chemistry. Two distinct mechanisms, the C-C and the C-H insertions of CH in C2H6, are characterized. The C-C insertion pathway is identified to have a large barrier of 34.5 kcal mol-1 and hence plays no significant role in kinetics. The C-H insertion pathway is found to have no barrier, leading to a highly vibrationally excited n-C3H7 radical, which rapidly dissociates (within 50 ps) to yield CH3 + C2H4 and H + C3H6 in a roughly 7:3 ratio. These findings are in good agreement with an experimental result that indicates that about 20% of the reaction goes to H + C3H6. The reaction of the electronically excited quartet state of the CH radical with C2H6 is examined for the first time and found to proceed as a direct H-abstraction via a small barrier of 0.4 kcal mol-1 to yield triplet CH2 and C2H5. The reaction on the quartet state surface is negligibly slow at low temperatures characteristic of interstellar environments but becomes important at high combustion temperatures.
Collapse
Affiliation(s)
- Thanh Lam Nguyen
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, United States
| | - Ajith Perera
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
38
|
Shi BX, Kapil V, Zen A, Chen J, Alavi A, Michaelides A. General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides. J Chem Phys 2022; 156:124704. [DOI: 10.1063/5.0087031] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
The O vacancy (Ov) formation energy, EOv, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g., the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute EOv accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods, such as local coupled cluster with single, double, and perturbative triple excitations. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rock salt MgO, producing the first accurate and well-converged determinations of EOv with this method. These reference values are used to benchmark exchange–correlation functionals in DFT, and we find that all the studied functionals underestimate EOv, with the average error decreasing along the rungs of Jacob’s ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.
Collapse
Affiliation(s)
- Benjamin X. Shi
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Venkat Kapil
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Churchill College, University of Cambridge, Storey’s Way, Cambridge CB3 0DS, United Kingdom
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Ji Chen
- School of Physics, Peking University, Beijing 100871, China
| | - Ali Alavi
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| |
Collapse
|
39
|
Xia Y, Long B, Lin S, Teng C, Bao JL, Truhlar DG. Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure. J Am Chem Soc 2022; 144:4828-4838. [PMID: 35262353 DOI: 10.1021/jacs.1c12324] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.
Collapse
Affiliation(s)
- Yu Xia
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Bo Long
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Shiru Lin
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
40
|
Ibrahim A, Roy PN. Three-body potential energy surface for para-hydrogen. J Chem Phys 2022; 156:044301. [PMID: 35105099 DOI: 10.1063/5.0076494] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Alexander Ibrahim
- Department of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
41
|
Lesiuk M. Quintic-scaling rank-reduced coupled cluster theory with single and double excitations. J Chem Phys 2022; 156:064103. [DOI: 10.1063/5.0071916] [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
|
42
|
Vasiliu M, Peterson KA, Marshall M, Zhu Z, Tufekci BA, Bowen KH, Dixon DA. Interaction of Th with H 0/-/+: Combined Experimental and Theoretical Thermodynamic Properties. J Phys Chem A 2022; 126:198-210. [PMID: 34989579 DOI: 10.1021/acs.jpca.1c07598] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
High-level electronic structure calculations of the low-lying energy electronic states for ThH, ThH-, and ThH+ are reported and compared to experimental measurements. The inclusion of spin-orbit coupling is critical to predict the ground-state ordering as inclusion of spin-orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH+. At the multireference spin-orbit SO-CASPT2 level, the ground states of ThH, ThH-, and ThH+ are predicted to be the 2Δ3/2, 3Φ2, and 3Δ1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH- photoelectron spectrum has many transitions, which approximately correlate with excitations of Th+ and/or Th. The adiabatic ionization energy of ThH including spin-orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th+H- ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH-, and ThH+ using the Feller-Peterson-Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.
Collapse
Affiliation(s)
- Monica Vasiliu
- Department of Chemistry and Biochemistry, The University of Alabama, Tuscaloosa, Alabama 35401, Unites States
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164, Unites States
| | - Mary Marshall
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, Unites States
| | - Zhaoguo Zhu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, Unites States
| | - Burak A Tufekci
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, Unites States
| | - Kit H Bowen
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, Unites States
| | - David A Dixon
- Department of Chemistry and Biochemistry, The University of Alabama, Tuscaloosa, Alabama 35401, Unites States
| |
Collapse
|
43
|
Semidalas E, Martin JM. The MOBH35 Metal–Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes. J Chem Theory Comput 2022; 18:883-898. [PMID: 35045709 PMCID: PMC8830049 DOI: 10.1021/acs.jctc.1c01126] [Citation(s) in RCA: 29] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We have revisited
the MOBH35 (Metal–Organic Barrier Heights,
35 reactions) benchmark [Iron, Janes, J.
Phys. Chem. A, 2019, 123 ( (17), ), 3761−378130973722; ibid. 2019, 123, 6379–6380] for realistic organometallic catalytic reactions, using both canonical
CCSD(T) and localized orbital approximations to it. For low levels
of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T)
perform well; for moderately strong levels of static correlation,
DLPNO-CCSD(T) and (T1) may break down catastrophically,
and PNO-LCCSD(T) is vulnerable as well. In contrast, LNO-CCSD(T) converges
smoothly to the canonical CCSD(T) answer with increasingly tight convergence
settings. The only two reactions for which our revised MOBH35 reference
values differ substantially from the original ones are reaction 9
and to a lesser extent 8, both involving iron. For the purpose of
evaluating density functional theory (DFT) methods for MOBH35, it
would be best to remove reaction 9 entirely as its severe level of
static correlation makes it just too demanding for a test. The magnitude
of the difference between DLPNO-CCSD(T) and DLPNO-CCSD(T1) is a reasonably good predictor for errors in DLPNO-CCSD(T1) compared to canonical CCSD(T); otherwise, monitoring all of T1, D1, max|tiA|, and 1/(εLUMO – εHOMO) should provide adequate warning
for potential problems. Our conclusions are not specific to the def2-SVP
basis set but are largely conserved for the larger def2-TZVPP, as
they are for the smaller def2-SV(P): the latter may be an economical
choice for calibrating against canonical CCSD(T). Finally, diagnostics
for static correlation are statistically clustered into groups corresponding
to (1) importance of single excitations in the wavefunction; (2a)
the small band gap, weakly separated from (2b) correlation entropy;
and (3) thermochemical importance of correlation energy, as well as
the slope of the DFT reaction energy with respect to the percentage
of HF exchange. Finally, a variable reduction analysis reveals that
much information on the multireference character is provided by T1, IND/Itot, and the exchange-based diagnostic A100[TPSS].
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| | - Jan M.L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| |
Collapse
|
44
|
Nagy PR, Gyevi-Nagy L, Kállay M. Basis set truncation corrections for improved frozen natural orbital CCSD(T) energies. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1963495] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| |
Collapse
|
45
|
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
|
46
|
Margócsy Á, Szabados Á. Dressing of Vertices by Cumulants in Multi-Reference Coupled Cluster. J Chem Theory Comput 2021; 17:6947-6964. [PMID: 34643380 DOI: 10.1021/acs.jctc.1c00730] [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/28/2022]
Abstract
A new scheme is introduced in Multi-Reference (MR) Coupled Cluster (CC) based on the MR Generalized Normal Ordering (MRGNO) and the corresponding MR Generalized Wick Theorem (MRGWT) of Kutzelnigg and Mukherjee. The key element is the identification of a structure in MRGWT generated terms, facilitated by Goldstone diagram techniques. This allows for bundling the many terms of the MRGWT expansion and introduces a hierarchy in the equations that can be harnessed in devising approximations. One- and two-particle interaction vertices are found to be uniformly substituted for their counterpart dressed by density cumulants. This allows for a straightforward rewriting of the ordinary energy expression of CC with interaction dressed (id) one- and two-particle terms and reveals the presence of three- and higher-rank dressed interaction vertices too. Cumulants appearing out of dressed interaction vertices contribute to the amplitude equations and can be interpreted to have an amplitude dressing role. Dressing of one- and two-particle interaction vertices is the most straightforward to implement and does not hinder computational feasibility, provided that the reference function involves strictly limited active space sizes. The Generalized Valence Bond wave function, underlying pilot numerical tests, fulfills this criterion. Results on multiple bond breaking scenarios point to the need of stepping beyond one- and two-particle id. An extremely simple version of incorporating amplitude dressing in addition to one- and two-particle id is seen to cure the potential energy curves remarkably, stimulating further investigations along this line.
Collapse
Affiliation(s)
- Ádám Margócsy
- ELTE Eötvös Loránd University, Faculty of Science, Institute of Chemistry Laboratory of Theoretical Chemistry Hevesy György PhD School of Chemistry, P.O. Box 32, Budapest 1518, Hungary
| | - Ágnes Szabados
- ELTE Eötvös Loránd University, Faculty of Science, Institute of Chemistry Laboratory of Theoretical Chemistry, P.O. Box 32, Budapest 1518, Hungary
| |
Collapse
|
47
|
Schraivogel T, Kats D. Accuracy of the distinguishable cluster approximation for triple excitations for open-shell molecules and excited states. J Chem Phys 2021; 155:064101. [PMID: 34391360 DOI: 10.1063/5.0059181] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The distinguishable cluster approximation for triple excitations has been applied to calculate thermochemical properties and excited states involving closed-shell and open-shell species, such as small molecules, 3d transition metal atoms, ozone, and an iron-porphyrin model. Excitation energies have been computed using the ΔCC approach by directly optimizing the excited states. A fixed-reference technique has been introduced to target selected spin-states for open-shell molecular systems. The distinguishable cluster approximation consistently improves coupled cluster with singles doubles and triples results for absolute and relative energies. For excited states dominated by a single configuration state function, the fixed-reference approach combined with high-level coupled-cluster methods has a comparable accuracy to the corresponding equation-of-motion coupled-cluster methods with a negligible amount of spin contamination.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
48
|
Nguyen TL, Peeters J. The CH(X 2Π) + H 2O reaction: two transition state kinetics. Phys Chem Chem Phys 2021; 23:16142-16149. [PMID: 34296725 DOI: 10.1039/d1cp02234b] [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/21/2022]
Abstract
The reaction of ground state methylidyne (CH) with water vapor (H2O) is theoretically re-investigated using high-level coupled cluster computations in combination with semi-classical transition state theory (SCTST) and two-dimensional master equation simulations. Insertion of CH into a H-O bond of H2O over a submerged barrier via a well-skipping mechanism yielding solely H and CH2O is characterized. The reaction kinetics is effectively determined by the formation of a pre-reaction van der Waals complex (PRC, HC-OH2) and its subsequent isomerization to activated CH2OH in competition with PRC re-dissociation. The tunneling effects are found to be minor, while variational effects in the PRC → CH2OH step are negligible. The calculated rate coefficient k(T) is nearly pressure-independent, but strongly depends on temperature with pronounced down-up behavior: a high value of 2 × 10-10 cm3 s-1 at 50 K, followed by a fairly steep decrease down to 8 × 10-12 cm3 s-1 at 900 K, but increasing again to 5 × 10-11 cm3 s-1 at 3500 K. Over the T-range of this work, k(T) can be expressed as: k(T, P = 0) = 2.31 × 10-11 (T/300 K)-1.615 exp(-38.45/T) cm3 s-1 for T = 50-400 K k(T, P = 0) = 1.15 × 10-12 (T/300 K)0.8637 exp(892.6/T) cm3 s-1 for T = 400-1000 K k(T, P = 0) = 4.57 × 10-15 (T/300 K)3.375 exp(3477.4/T) cm3 s-1 for T = 1000-3500 K.
Collapse
Affiliation(s)
- Thanh Lam Nguyen
- Quantum Theory Project, Department of Chemistry and Physics, University of Florida, Gainesville, FL 32611, USA.
| | | |
Collapse
|
49
|
Accurate estimation of singlet-triplet gap of strongly correlated systems by CCSD(T) method using improved orbitals. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113326] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
50
|
Ciborowski SM, Liu G, Blankenhorn M, Harris RM, Marshall MA, Zhu Z, Bowen KH, Peterson KA. The electron affinity of the uranium atom. J Chem Phys 2021; 154:224307. [PMID: 34241193 DOI: 10.1063/5.0046315] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U-. Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin-orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th-, the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both cases, the ground state of the anion corresponds to electron attachment to the 6d orbital.
Collapse
Affiliation(s)
- Sandra M Ciborowski
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Gaoxiang Liu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Moritz Blankenhorn
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Rachel M Harris
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Mary A Marshall
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Zhaoguo Zhu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Kit H Bowen
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99162, USA
| |
Collapse
|