Canonical transformation (CT) theory describes dynamic correlation in multireference systems with large active spaces. Here we discuss CT theory's intruder state problem and why our previous approach of overlap matrix truncation becomes infeasible for sufficiently large active spaces. We propose the use of strongly and weakly contracted excitation operators as alternatives for dealing with intruder states in CT theory. The performance of these operators is evaluated for the H(2)O, N(2), and NiO molecules, with comparisons made to complete active space second order perturbation theory and Davidson-corrected multireference configuration interaction theory. Finally, using a combination of strongly contracted CT theory and orbital-optimized density matrix renormalization group theory, we evaluate the singlet-triplet gap of free base porphin using an active space containing all 24 out-of-plane 2p orbitals. Modeling dynamic correlation with an active space of this size is currently only possible using CT theory.
Canonical transformation (CT) theory describes dynamic correlation in multireference systems with large active spaces. Here we discuss CT theory's intruder state problem and why our previous approach of overlap matrix truncation becomes infeasible for sufficiently large active spaces. We propose the use of strongly and weakly contracted excitation operators as alternatives for dealing with intruder states in CT theory. The performance of these operators is evaluated for the H(2)O, N(2), and NiO molecules, with comparisons made to complete active space second order perturbation theory and Davidson-corrected multireference configuration interaction theory. Finally, using a combination of strongly contracted CT theory and orbital-optimized density matrix renormalization group theory, we evaluate the singlet-triplet gap of free base porphin using an active space containing all 24 out-of-plane 2p orbitals. Modeling dynamic correlation with an active space of this size is currently only possible using CT theory.
We describe the joint application of the density matrix renormalization group and canonical transformation theory to multireference quantum chemistry. The density matrix renormalization group provides the ability to describe static correlation in large active spaces, while the canonical transformation theory provides a high-order description of the dynamic correlation effects. We demonstrate the joint theory in two benchmark systems designed to test the dynamic and static correlation capabilities of the methods, namely, (i) total correlation energies in long polyenes and (ii) the isomerization curve of the [Cu(2)O(2)](2+) core. The largest complete active spaces and atomic orbital basis sets treated by the joint DMRG-CT theory in these systems correspond to a (24e,24o) active space and 268 atomic orbitals in the polyenes and a (28e,32o) active space and 278 atomic orbitals in [Cu(2)O(2)](2+).
The canonical transformation (CT) theory has been developed as a multireference electronic structure method to compute high-level dynamic correlation on top of a large active space reference treated with the ab initio density matrix renormalization group method. This article describes a parallelized algorithm and implementation of the CT theory to handle large computational demands of the CT calculation, which has the same scaling as the coupled cluster singles and doubles theory. To stabilize the iterative solution of the CT method, a modification to the CT amplitude equation is introduced with the inclusion of a level shift parameter. The level-shifted condition has been found to effectively remove a type of intruder state that arises in the linear equations of CT and to address the discontinuity problems in the potential energy curves observed in the previous CT studies.
We present two efficient and intruder-free methods for treating dynamic correlation on top of general multiconfiguration reference wave functions — including such as obtained by the density matrix renormalization group (DMRG) with large active spaces. The new methods are the second order variant of the recently proposed multireference linearized coupled cluster method (MRLCC) [Sharma, S.; Alavi, A. J. Chem. Phys. 2015, 143, 102815] and of N-electron valence perturbation theory (NEVPT2), with expected accuracies similar to MRCI+Q and (at least) CASPT2, respectively. Great efficiency gains are realized by representing the first order wave function with a combination of internal contraction (IC) and matrix product state perturbation theory (MPSPT). With this combination, only third order reduced density matrices (RDMs) are required. Thus, we obviate the need for calculating (or estimating) RDMs of fourth or higher order; these had so far posed a severe bottleneck for dynamic correlation treatments involving the large active spaces accessible to DMRG. Using several benchmark systems, including first and second row containing small molecules, Cr 2, pentacene, and oxo-Mn(Salen), we show that active spaces containing at least 30 orbitals can be treated using this method. On a single node, MRLCC2 and NEVPT2 calculations can be performed with over 550 and 1100 virtual orbitals, respectively. We also critically examine the errors incurred due to the three sources of errors introduced in the present implementation – calculating second order instead of third order energy corrections, use of internal contraction, and approximations made in the reference wave function due to DMRG.
We report development of the multireference configuration interaction (MRCI) method that can use active space scalable to much larger size references than has previously been possible. The recent development of the density matrix renormalization group (DMRG) method in multireference quantum chemistry offers the ability to describe static correlation in a large active space. The present MRCI method provides a critical correction to the DMRG reference by including high-level dynamic correlation through the CI treatment. When the DMRG and MRCI theories are combined (DMRG-MRCI), the full internal contraction of the reference in the MRCI ansatz, including contraction of semi-internal states, plays a central role. However, it is thought to involve formidable complexity because of the presence of the five-particle rank reduced-density matrix (RDM) in the Hamiltonian matrix elements. To address this complexity, we express the Hamiltonian matrix using commutators, which allows the five-particle rank RDM to be canceled out without any approximation. Then we introduce an approximation to the four-particle rank RDM by using a cumulant reconstruction from lower-particle rank RDMs. A computer-aided approach is employed to derive the exceedingly complex equations of the MRCI in tensor-contracted form and to implement them into an efficient parallel computer code. This approach extends to the size-consistency-corrected variants of MRCI, such as the MRCI+Q, MR-ACPF, and MR-AQCC methods. We demonstrate the capability of the DMRG-MRCI method in several benchmark applications, including the evaluation of single-triplet gap of free-base porphyrin using 24 active orbitals.
The recent development of the density matrix renormalization group (DMRG) method in multireference quantum chemistry makes it practical to evaluate static correlation in a large active space, while dynamic correlation provides a critical correction to the DMRG reference for strong-correlated systems and is usually obtained using multireference perturbation (MRPT) or configuration interaction (MRCI) methods with internal contraction (ic) approximation. These methods can use an active space scalable to relatively larger size references than has previously been possible. However, they are still hardly applicable to systems with an active space larger than 30 orbitals and/or a large basis set because of high computation and storage costs of high-order reduced density matrices (RDMs) and the crucial dependence of the MRCI Hamiltonian dimension on the number of virtual orbitals. In this work, we propose a new effective implementation of DMRG-MRCI, in which we use reconstructed CASCI-type configurations from DMRG wave function via the entropy-driving genetic algorithm (EDGA) [Luo et al. J. Chem. Theory Comput. 2017, 13, 4699−4710.] and integrate it with MRCI by an external contraction (ec) scheme. This bypasses the bottleneck of computing high-order RDMs in traditional DMRG dynamic correlation methods with ic approximation, and the number of MRCI configurations is not dependent on the number of virtual orbitals. Therefore, the DMRG-ec-MRCI method is promising for dealing with a larger active space than 30 orbitals and large basis sets. We demonstrate the capability of our DMRG-ec-MRCI method in several benchmark applications, including the evaluation of the potential energy curve of Cr 2, single-triplet gaps of higher n-acene molecules, and the energy of the Eu-BTBP(NO 3 ) 3 complex.
In this thesis, we will describe several extensions of the Density Matrix Renormalization Group (DMRG) and Canonical Transformation Theory (CT). In the first part we describe a new way to solve for excited states in DMRG. To overcome the limitations of the traditional state-averaging approaches in excited state calculations, where one solves for all states between the ground state and excited state of interest, we have investigated a number of new excited state algorithms. Building on the work of van der Vorst and Sleijpen (SIAM J. Matrix Anal. Appl., 17, 401 (1996)), we have implemented Harmonic Davidson and State-Averaged Harmonic Davidson algorithms within the context of DMRG. We have assessed their accuracy and stability of convergence in complete active space DMRG calculations on the low-lying excited states in the acenes ranging from naphthalene to pentacene. We find that both algorithms offer increased accuracy over the traditional State-Averaged Davidson approach, and in particular, the State-Averaged Harmonic Davidson algorithm offers an optimal combination of accuracy and stability in convergence. In the second part, we propose an analytic response theory for DMRG whereby response properties correspond to analytic derivatives of DMRG observables with respect to the applied perturbations. Both static and frequencydependent response theories are formulated and implemented. We evaluate our pilot implementation by calculating static and frequency-dependent polar- isabilities of short oligo-di-acetylenes. The analytic response theory is competitive with dynamical DMRG methods and yields significantly improved accuracies when using a small number of density matrix renormalisation group states. Strengths and weaknesses of the analytic approach are discussed. In the third part, we describe how to calculate density matrices in CT theory. Density matrices are useful not only for computing observables but also for characterizing the nature of individual states. We demonstrate this with a preliminary application of the CT theory to understand the low-lying excited states of oligo-phenylvinylenes. We finish by presenting the theory and equations for calculating the response (e.g. to an external field) in the CT and DMRG-CT theories.