Abstract:An optimized implementation of block-correlated coupled cluster theory based on the generalized valence bond wave function (GVB-BCCC) for the singlet ground state of strongly correlated systems is presented. The GVB-BCCC method with two-pair correlation (GVB-BCCC2b) or up to three-pair correlation (GVB-BCCC3b) will be the focus of this work. Three major techniques have been adopted to dramatically accelerate GVB-BCCC2b and GVB-BCCC3b calculations. First, the GVB-BCCC2b and GVB-BCCC3b codes are noticeably optimized by removing redundant calculations. Second, independent amplitudes are identified by constraining excited configurations to be pure singlet states and only independent amplitudes need to be solved. Third, an incremental updating scheme for the amplitudes in solving the GVB-BCCC equations is adopted. With these techniques, accurate GVB-BCCC3b calculations are now accessible for systems with relatively large active spaces (50 electrons in 50 orbitals) and GVB-BCCC2b calculations are affordable for systems with much larger active spaces. We have applied GVB-BCCC methods to investigate three typical kinds of systems: polyacenes, pentaprismane, and [Cu2O2]2+ isomers. For polyacenes, we demonstrate that GVB-BCCC3b can capture more than 94% of the total correlation energy even for 12-acene with 50 π electrons. For the potential energy curve of simultaneously stretching 15 C-C bonds in pentaprismane, our calculations show that the GVB-BCCC3b results are quite close to the results from the density matrix renormalization group (DMRG) over the whole range. For two dinuclear copper oxide isomers, their relative energy predicted by GVB-BCCC3b is also in good accord with the DMRG result. All calculations show that the inclusion of three-pair correlation in GVB-BCCC is critical for accurate descriptions of strongly correlated systems.

GUGA-based MRCI approach with Core-Valence Separation Approximation (CVS) for the calculation of the Core-Excited States of molecules

Externally-Contracted Multireference Configuration Interaction Method Using a DMRG Reference Wave Function

A DFT/MRCI Hamiltonian Parameterized Using Only Ab Initio Data: II. Core-Excited States

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry.

Theoretical study on the core-excited states of the allyl using CVS-icMRCISD method

Spin-Adapted Externally Contracted Multireference Configuration Interaction Method Based on Selected Reference Configurations

A variational Monte Carlo approach for core excitations

Efficient Implementation of Block-Correlated Coupled Cluster Theory Based on the Generalized Valence Bond Reference for Strongly Correlated Systems

Excited State Structure and Minimum Energy Conical Intersection Optimization Using DFT/MRCI

Enhancing the Accuracy of XPS Calculations: Exploring Hybrid Basis Set Schemes for CVS-EOMIP-CCSD Calculations

Complex Multireference Configuration Interaction Calculations for the K-vacancy Auger States of N(q+) (q = 2-5) Ions.

An Efficient Implementation of the "Cluster-in-molecule" Approach for Local Electron Correlation Calculations.

Excited States via Coupled Cluster Theory without Equation-of-Motion Methods: Seeking Higher Roots with Application to Doubly Excited States and Double Core Hole States

Renormalized Internally-Contracted Multireference Coupled Cluster with Perturbative Triples

Core Binding Energy Calculations: A Scalable Approach with the Quantum Embedding-Based Equation-of-Motion Coupled-Cluster Method

Optimization of highly excited matrix product states with an application to vibrational spectroscopy

Multireference Second Order Perturbation Theory with a Simplified Treatment of Dynamical Correlation.

Automatic Selection of Active Orbitals from Generalized Valence Bond Orbitals.

Comparison of some multireference electronic structure methods in illustrative applications

High-level Multireference Methods in the Quantum-Chemistry Program System COLUMBUS: Analytic MR-CISD and MR-AQCC Gradients and MR-AQCC-LRT for Excited States, GUGA Spin–orbit CI and Parallel CI Density

Density-matrix Renormalization Group Algorithm with Multi-Level Active Space.