We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C 2 and Cr 2, which are notorious for their multireference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.
We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C 2 and Cr 2, which are notorious for their multireference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.
We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C2 and Cr2, which are notorious for their multi-reference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.
We describe a low cost alternative to the standard variational DMRG (density matrix renormalization group) algorithm that is analogous to the combination of the selected configuration interaction plus perturbation theory (SCI+PT). We denote the resulting method p-DMRG (perturbative DMRG) to distinguish it from the standard variational DMRG. p-DMRG is expected to be useful for systems with very large active spaces, for which variational DMRG becomes too expensive. Similar to SCI+PT, in p-DMRG, a zeroth-order wave function is first obtained by a standard DMRG calculation but with a small bond dimension. Then, the residual correlation is recovered by a second-order perturbative treatment. We discuss the choice of partitioning for perturbation theory, which is crucial for its accuracy and robustness. To circumvent the problem of a large bond dimension in the first-order wave function, we use a sum of matrix product states to expand the first-order wave function, yielding substantial savings in computational cost and memory. We also propose extrapolation schemes to reduce the errors in the zeroth- and first-order wave functions. Numerical results for Cr 2 with a (28e, 76o) active space and 1,3-butadiene with an (22e, 82o) active space reveal that p-DMRG provides ground state energies of a similar quality to variational DMRG with very large bond dimensions but at a significantly lower computational cost. This suggests that p-DMRG will be an efficient tool for benchmark studies in the future.
This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.
We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.
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.
This paper describes an interface between the density matrix renormalization group (DMRG) method and the complete active-space self-consistent field (CASSCF) method and its analytical gradient, as well as an extension to the second-order perturbation theory (CASPT2) method. This interfacing allows large active-space multi-reference computations to be easily performed. The interface and its extension are both implemented in terms of reduced density matrices (RDMs) which can be efficiently computed via the DMRG sweep algorithm. We also present benchmark results showing that, in practice, the DMRG-CASSCF calculations scale with active-space size in a polynomial manner in the case of quasi-1D systems. Geometry optimization of a binuclear iron-sulfur cluster using the DMRG-CASSCF analytical gradient is demonstrated, indicating that the inclusion of the valence p-orbitals of sulfur and double-shell d-orbitals of iron lead to non-negligible changes in the geometry compared to the results of small active-space calculations. With the exception of the selection of M values, many computational settings in these practical DMRG calculations have been tuned and black-boxed in our interface, and so the resulting DMRG-CASSCF and DMRG-CASPT2 calculations are now available to novice users as a common tool to compute strongly correlated electronic wavefunctions.