We formulated and implemented a perturbative triples correction for the state-specific multireference coupled cluster approach with singles and doubles suggested by Mukherjee and co-workers, Mk-MRCCSD [Mol. Phys. 94, 157 (1998)]. Our derivation of the energy correction [Mk-MRCCSD(T)] is based on a constrained search for stationary points of the Mk-MRCC energy functional together with a perturbative expansion with respect to the appearing triples cluster operator. The Lambda-Mk-MRCCSD(T) approach derived in this way consists in (1) a correction to the off-diagonal matrix elements of the effective Hamiltonian which is unique to coupled cluster methods based on the Jeziorski-Monkhorst ansatz, and (2) an asymmetric energy correction to the diagonal elements of the effective Hamiltonian. The Mk-MRCCSD(T) correction is obtained from the Lambda-Mk-MRCCSD(T) method by approximating the singles and doubles Lagrange multipliers with the corresponding cluster amplitudes. We investigate the performance of the Mk-MRCCSD(T) method by applying it to the potential energy curve of the BeH(2) model and F(2) and the geometry and harmonic vibrational frequencies of ozone. Computation of the energy difference between the mono- and bicyclic forms of the 2,6-pyridyne diradical illustrates the potential of Mk-MRCCSD(T) as a tool for the study of realistic chemical problems requiring multireference zeroth-order wave functions.
The general theory of analytic energy gradients is presented for the state-specific multireference coupled cluster method introduced by Mukherjee and co-workers [Mol. Phys. 94, 157 (1998)], together with an implementation within the singles and doubles approximation, restricted to two closed-shell determinants and Hartree-Fock orbitals. Expressions for the energy gradient are derived based on a Lagrangian formalism and cast in a density-matrix notation suitable for implementation in standard quantum-chemical program packages. In the present implementation, we exploit a decomposition of the multireference coupled cluster gradient expressions, i.e., lambda equations and the corresponding density matrices, into a so-called single-reference part for each reference determinant and a coupling term. Our implementation exhibits the proper scaling, i.e., O(dN6) with d as the number of reference determinants and N as the number of orbitals, and it is thus suitable for large-scale applications. The applicability of our multireference coupled cluster gradients is illustrated by computations for the equilibrium geometry of the 2,6-isomers of pyridyne and the pyridynium cation. The results are compared to those from single-reference coupled cluster calculations and are discussed with respect to the future perspectives of multireference coupled cluster theory.
High-level ab initio benchmark calculations of the (15)N and (31)P NMR chemical shielding constants for a representative set of molecules are presented. The computations have been carried out at the Hartree-Fock self-consistent field (HF-SCF), density functional theory (DFT) (B-P86 and B3-LYP), second-order Moller-Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), and CCSD augmented by a perturbative treatment of triple excitations [CCSD(T)] level of theory using basis sets of triple zeta quality or better. The influence of the geometry, the treatment of electron correlation, as well as basis set and zero-point vibrational effects on the shielding constants are discussed and the results are compared to gas-phase experimental shifts. As for the first time a study using high-level post-HF methods is carried out for a second-row element, we also propose a family of basis sets suitable for the computation of (31)P shielding constants. The mean deviations observed for (15)N and (31)P are 0.9 [CCSD(T)/13s9p4d3f] and -3.3 ppm [CCSD(T)/15s12p4d3f2g], respectively, when corrected for zero-point vibrational effects. Results obtained at the DFT level of theory are of comparable accuracy to MP2 for (15)N and of comparable accuracy to HF-SCF for (31)P. However, they are not improved by inclusion of zero-point vibrational effects. The PN molecule is an especially interesting case with exceptionally large electron correlation effects on shielding constants beyond MP2 which, therefore, represents an excellent example for further benchmark studies.
Analytic gradients for the state-specific multireference coupled-cluster method suggested by Mahapatra et al. [Mol. Phys. 94, 157 (1998)] (Mk-MRCC) are reported within the singles and doubles approximation using two-configurational self-consistent field (TCSCF) orbitals. The present implementation extends our previous work on Mk-MRCC gradients [E. Prochnow et al., J. Chem. Phys. 131, 064109 (2009)] which is based on restricted Hartree-Fock orbitals and consequently the main focus of the present paper is on the treatment of orbital relaxation at the TCSCF level using coupled-perturbed TCSCF theory. Geometry optimizations on m-arynes and nitrenes are presented to illustrate the influence of the orbitals on the computed equilibrium structures. The results are compared to those obtained at the single-reference coupled-cluster singles and doubles and at the Mk-MRCC singles and doubles level of theory when using restricted Hartree-Fock orbitals.
A scheme for the parallel calculation of energies at the coupled-cluster singles, doubles, and triples (CCSDT) level of theory, several approximate iterative CCSDT schemes (CCSDT-1a, CCSDT-1b, CCSDT-2, CCSDT-3, and CC3), and for the state-specific multireference coupled-cluster ansatz suggested by Mukherjee with a full treatment of triple excitations (Mk-MRCCSDT) is presented. The proposed scheme is based on the adaptation of a highly efficient serial coupled-cluster code leading to a communication-minimized implementation by parallelizing the time-determining steps. The parallel algorithm is tailored for affordable cluster architectures connected by standard communication networks such as Gigabit Ethernet. In this way, CCSDT and Mk-MRCCSDT computations become feasible even for larger molecular systems and basis sets. An analysis of the time-determining steps for CCSDT and Mk-MRCCSDT, namely the computation of the triple-excitation amplitudes and their individual contributions, is carried out. Benchmark calculations are presented for the N2O, ozone, and benzene molecules, proving that the parallelization of these steps is sufficient to obtain an efficient parallel scheme. A first application to the case of 2,6-pyridyne using a triple-ζ quality basis (222 basis functions) is presented demonstrating the efficiency of the current implementation.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.