1
|
Chakravarti D, Sen S, Mukherjee D. A systematic improvement to UGA-SSMRCCSD equations and its implication for potential energy curves. J Chem Phys 2023; 159:134102. [PMID: 37787129 DOI: 10.1063/5.0168941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 09/14/2023] [Indexed: 10/04/2023] Open
Abstract
The Unitary Group Adaptation (UGA) offers a very compact and efficient spin adaptation strategy for any spin-free Hamiltonian in a many body framework. Our use of UGA in the context of state-specific (SS) Jeziorski-Monkhorst Ansatz based multireference coupled cluster (MRCC) theory obviates the non-commutativity between the spin-free cluster operators via a normal ordered exponential parametrization in the wave operator. A previous formulation of UGA-SSMRCC by us [R. Maitra, D. Sinha, and D. Mukherjee, J. Chem. Phys. 137, 024105 (2012)], using the same ansatz, employed certain sufficiency conditions to reach the final working equations, which cannot be improved systematically. In this article, we will present a more rigorous formulation that follows from an exact factorization of the unlinked terms of the Bloch equation, resulting in equations on which a hierarchy of approximations can be systematically performed on the emergent additional terms. This derivation was shown in our recent article [D. Chakravarti, S. Sen, and D. Mukherjee, Mol. Phys. 119, e1979676 (2021)] in the context of a single open shell CC formalism and was applied to spectroscopic energy differences where the contribution of the new terms was found to be of the order of ∼0.001 eV for ionization potential, electron affinity, and excitation energy. In the current work, we will present a comparison between the earlier and current formulations via both a theoretical analysis and a numerical demonstration of the dramatic effect of the additional terms brought in by the factorization on potential energy curves. The contribution of such terms was found to gain importance with an increase in the number of singly occupied active orbitals in the model space functions.
Collapse
Affiliation(s)
- Dibyajyoti Chakravarti
- Centre for Quantum Engineering, Research, and Education (CQuERE), TCG CREST, Kolkata, India
| | - Sangita Sen
- Department of Chemical Sciences, Indian Institute of Science, Education and Research, Kolkata, India
| | - Debashis Mukherjee
- Centre for Quantum Engineering, Research, and Education (CQuERE), TCG CREST, Kolkata, India
| |
Collapse
|
2
|
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
|
3
|
Margócsy Á, Szabados Á. Ring coupled cluster doubles at the multireference level. J Chem Phys 2020; 152:204114. [PMID: 32486660 DOI: 10.1063/5.0005075] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A ring approximation within an internally contracted multireference (MR) Coupled Cluster (CC) framework is worked out and tested. Derivation of equations utilizes MR based, generalized normal ordering and the corresponding generalized Wick-theorem (MR-GWT). Contractions among cluster operators are avoided by adopting a normal ordered exponential ansatz. The original version of the MR ring CC doubles (MR-rCCD) equations [Á. Szabados and Á. Margócsy, Mol. Phys. 115, 2731 (2017)] is rectified in two aspects. On the one hand, over-completeness of double excitations is treated by relying on the concept of frames. On the other hand, restriction on the maximal cumulant rank is lifted from two to four. This is found essential for obtaining reliable correlation corrections to the energy. The MR function underlying the approach is provided by the Generalized Valence Bond (GVB) model. The pair structure of the reference ensures a fragment structure of GVB cumulants. This represents a benefit when evaluating cumulant contractions appearing as a consequence of MR-GWT. In particular, cumulant involving terms remain less expensive than their traditional, pair-contracted counterpart, facilitating an O(N6) eventual scaling of the proposed MR-rCCD method. Pilot applications are presented for covalent bond breaking, deprotonation energies, and torsional potentials.
Collapse
Affiliation(s)
- Á Margócsy
- Laboratory of Theoretical Chemistry, Institute of Chemistry, Faculty of Science, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Á Szabados
- Laboratory of Theoretical Chemistry, Institute of Chemistry, Faculty of Science, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
4
|
Wu CH, Magers DB, Harding LB, Klippenstein SJ, Allen WD. Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory. J Chem Theory Comput 2020; 16:1511-1525. [PMID: 32073856 DOI: 10.1021/acs.jctc.9b00966] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H2+ Be, H + NH2 → H2 + NH, CH3 + C2H5 → CH4 + C2H4, H + C2H5 → H2 + C2H4, and H + HCO → H2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH2, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH3 + C2H5, H + C2H5, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.
Collapse
Affiliation(s)
- Chia-Hua Wu
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - D Brandon Magers
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States.,Department of Chemistry and Physics, Belhaven University, Jackson, Mississippi 39202, United States
| | - Lawrence B Harding
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Stephen J Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Wesley D Allen
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
5
|
Affiliation(s)
- Ágnes Szabados
- Laboratory of Theoretical Chemistry, Loránd Eötvös University, Budapest, Hungary
| | - Ádám Margócsy
- Laboratory of Theoretical Chemistry, Loránd Eötvös University, Budapest, Hungary
| |
Collapse
|
6
|
Sen A, Sen S, Mukherjee D. Aspects of Size-Consistency of Orbitally Noninvariant Size-Extensive Multireference Perturbation Theories: A Case Study Using UGA-SSMRPT2 as a Prototype. J Chem Theory Comput 2015; 11:4129-45. [PMID: 26575908 DOI: 10.1021/acs.jctc.5b00457] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Profiling a potential energy surface (PES), all the way to dissociate a molecular state into particular fragments and to display real or avoided crossings, requires a multireference description and the maintenance of size-consistency. The many body methods, which suit this purpose, should thus be size-extensive. Size-extensive theories, which are invariant with respect to transformation among active orbitals are, in principle, size-consistent. Relatively cheaper size-extensive theories, which do not possess this invariance, can still be size-consistent if the active orbitals are localized on the asymptotic fragments. Such methods, if perturbative in nature, require the use of an unperturbed Hamiltonian, which has orbital invariance with respect to the transformation within active, core, and virtual orbitals. The principal focus of this paper is to numerically realize size-consistency with localized active orbitals using our recently developed orbitally noninvariant Unitary Group Adapted State Specific Multireference second order Perturbation Theory (UGA-SSMRPT2) as a prototype method. Our findings expose certain generic potential pitfalls of size-extensive but orbitally noninvariant MRPT theories, which are mostly related to the inability of reaching proper localized active orbitals in the fragments due to the artifacts of the orbital generation procedure. Despite the invariance of the zeroth order CAS function, lack of invariance of the MRPT itself then leads to size-inconsistency. In particular, reaching symmetry broken fragment active orbitals is an issue of concern where suitable state-averaging might ameliorate the problem, but then one has to abandon full orbital optimization. Additionally, there can be situations where the orbitals of the fragment reached as an asymptote of the supermolecule are not the same as those obtained from the optimization of the fragments individually and will require additional transformation. Moreover, for a certain PES, one may either abandon the use of optimized orbitals for that state to preserve proper symmetry and degeneracy in the fragment orbitals or be satisfied with the use of optimized orbitals, which generate broken symmetric orbitals in the fragmentation limit. All these pathologies are illustrated using the PES of various electronic states of multiply bonded systems like N2, C2H2, HCN, C2, and O2. Subject to such proviso, the UGA-SSMRPT2 turns out to be an excellent theory for studying the PES leading to fragmentation of strongly correlated systems satisfying the requirements of size-consistency with localized active orbitals. An unexpected spin-off of our studies is the realization that the size-inextensive MRMP2, which bears a close structural similarity with our theory, might under certain situations display size-consistency. We analyze this feature concretely in our paper. Our studies may serve as a benchmark for monitoring numerically the size-consistency of any state specific multireference theory which is size-extensive but not invariant with respect to transformation of active orbitals.
Collapse
Affiliation(s)
- Avijit Sen
- Raman Center for Atomic, Molecular and Optical Sciences, Indian Association for the Cultivation of Science , Kolkata 700 032, India
| | - Sangita Sen
- Raman Center for Atomic, Molecular and Optical Sciences, Indian Association for the Cultivation of Science , Kolkata 700 032, India
| | - Debashis Mukherjee
- Raman Center for Atomic, Molecular and Optical Sciences, Indian Association for the Cultivation of Science , Kolkata 700 032, India
| |
Collapse
|
7
|
Sen A, Sen S, Samanta PK, Mukherjee D. Unitary group adapted state specific multireference perturbation theory: Formulation and pilot applications. J Comput Chem 2015; 36:670-88. [DOI: 10.1002/jcc.23851] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2014] [Revised: 01/01/2015] [Accepted: 01/05/2015] [Indexed: 11/11/2022]
Affiliation(s)
- Avijit Sen
- Raman Centre for Atomic, Molecular and Optical Sciences; Indian Association for the Cultivation of Sciences; Jadavpur Kolkata India
| | - Sangita Sen
- Raman Centre for Atomic, Molecular and Optical Sciences; Indian Association for the Cultivation of Sciences; Jadavpur Kolkata India
| | - Pradipta Kumar Samanta
- Raman Centre for Atomic, Molecular and Optical Sciences; Indian Association for the Cultivation of Sciences; Jadavpur Kolkata India
| | - Debashis Mukherjee
- Raman Centre for Atomic, Molecular and Optical Sciences; Indian Association for the Cultivation of Sciences; Jadavpur Kolkata India
| |
Collapse
|