arXiv:1809.07732 · Journal of Chemical Theory and Computation
Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer
Fabian M. Faulstich, Mihály Máté, Andre Laestadius, Mihály A. Csirik, Libor Veis, Andrej Antalík, Jiří Brabec, Reinhold Schneider, Jiřı́ Pittner, Simen Kvaal, Örs Legeza
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO-LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.
In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO-LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.
In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schrödinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method combines the single-reference coupled-cluster (CC) approach with an approximate reference calculation in a subspace (complete active space (CAS)) of the considered Hilbert space that covers the static correlation. A one-particle spectral gap assumption is introduced, separating the CAS from the remaining Hilbert space. This replaces the nonexisting or nearly nonexisting gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital usually encountered in standard single-reference quantum chemistry. The analysis covers, in particular, CC methods tailored by tensor-network states (TNS-TCC methods). The problem is formulated in a nonlinear functional analysis framework, and, under certain conditions such as the aforementioned gap, local uniqueness and existence are proved using Zarantonello's lemma. From the Aubin--Nitsche-duality method, a quadratic error bound valid for TNS-TCC methods is derived, e.g., for linear-tensor-network TCC schemes using the density matrix renormalization group method.
In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schr\"odinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method combines the single-reference coupled-cluster (CC) approach with an approximate reference calculation in a subspace [complete active space (CAS)] of the considered Hilbert space that covers the static correlation. A one-particle spectral gap assumption is introduced, separating the CAS from the remaining Hilbert space. This replaces the nonexisting or nearly nonexisting gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital usually encountered in standard single-reference quantum chemistry. The analysis covers, in particular, CC methods tailored by tensor-network states (TNS-TCC methods). The problem is formulated in a nonlinear functional analysis framework, and, under certain conditions such as the aforementioned gap, local uniqueness and existence are proved using Zarantonello's lemma. From the Aubin--Nitsche-duality method, a quadratic error bound valid for TNS-TCC methods is derived, e.g., for linear-tensor-network TCC schemes using the density matrix renormalization group method.
In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.
We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)—we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.
Three fully variational formulations of the complete-active space coupled-cluster method are derived. The formulations include the ability to approximate the model vectors by smooth manifolds, thereby opening up the possibility for overcoming the exponential wall of scaling for model spaces of complete-active space type. In particular, model vectors of matrix-product states are considered, and it is argued that the present variational formulation allows not only favorably scaling multireference coupled-cluster calculations but also systematic correction of tailored coupled-cluster calculations and of quantum chemical density-matrix renormalization group methods, which are fast and polynomial scaling but lack the ability to properly resolve dynamical correlation at chemical accuracy. The extension of the variational formulations to the time domain is also discussed, with derivations of abstract evolution equations.
In the past decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favorable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate “post-DMRG” treatment of dynamic correlation based on the tailored coupled cluster (CC) theory in which the DMRG method is responsible for the proper description of nondynamic correlation, whereas dynamic correlation is incorporated through the framework of the CC theory. We illustrate the potential of this method on prominent multireference systems, in particular, N 2 and Cr 2 molecules and also oxo-Mn(Salen), for which we have performed the first post-DMRG computations in order to shed light on the energy ordering of the lowest spin states.