We have compared different strategies for ab initio quantum chemistry density matrix renormalization group treatments. The two starting orbital blocks include all valence and active orbitals of the reference complete active space self consistent field wave function. To generate the remaining blocks, we propose following the order of the contributions to the correlation energy: a posteriori using approximate occupation numbers or a priori, using a Møller–Plesset type of arguments, by explicit evaluation of second-order interactions. We have compared two different schemes for orbital localization to identify the important and less important orbital interactions and simplify the generation of the orbital blocks. To truncate the expansion we have compared two approaches, keeping constant the number m of components or the threshold λ to fix the residue of the expansion at each step. The extrapolation of the energies is found to provide accurate estimates of the full configuration interaction energy, making the expansion independent on the actual values of the two parameters m and λ. We propose to generate the factors for the two blocks from ground and excited eigenvectors of the Hamiltonian matrix.
We have compared different strategies for ab initio quantum chemistry density matrix renormalization group treatments. The two starting orbital blocks include all valence and active orbitals of the reference complete active space self consistent field wave function. To generate the remaining blocks, we propose following the order of the contributions to the correlation energy: a posteriori using approximate occupation numbers or a priori, using a Møller–Plesset type of arguments, by explicit evaluation of second-order interactions. We have compared two different schemes for orbital localization to identify the important and less important orbital interactions and simplify the generation of the orbital blocks. To truncate the expansion we have compared two approaches, keeping constant the number m of components or the threshold λ to fix the residue of the expansion at each step. The extrapolation of the energies is found to provide accurate estimates of the full configuration interaction energy, making the expansion independent on the actual values of the two parameters m and λ. We propose to generate the factors for the two blocks from ground and excited eigenvectors of the Hamiltonian matrix.
Despite the success of modern quantum chemistry in predicting properties of organic molecules, the treatment of inorganic systems, which have many close lying states, remains out of quantitative reach for current methods. To treat non-dynamic correlation, we take advantage of the density matrix renormalization group (DMRG) method that has become very successful in the field of solid state physics. We present a detailed study of the DMRG method, and we pay special attention to the evolution of the understanding behind the mathematical structure of the DMRG wave function. Our primary goal is to develop a density matrix renormalization group self--consistent--field (DMRG-SCF) approach, analogous to the complete active space self--consistent field (CASSCF) method, but dealing with large active spaces that are too demanding for the full configuration interaction (FCI) method. \nAs a first step towards such a DMRG-SCF procedure, we present a spin-adapted DMRG algorithm designed to target spin- and spatial-symmetry states that are hard to obtain while using an unrestricted algorithm. \nOur next step is a modification of the DMRG algorithm to obtain decreasing energy at every step during the sweep. This monotonically convergent DMRG scheme lets us obtain the two-body density matrix as a by--product of the existing procedure without any additional cost in storage. Additionally, the two-body density matrix produced at convergence using this scheme is free from the N-representability problem that is present when the two--body density matrix is produced with the two-site DMRG scheme without additional storage cost. Finally, taking advantage of the modifications developed herein, we present results obtained from our DMRG-SCF method. Lastly, we discuss possible ways of incorporating dynamical correlation into the DMRG scheme, in order to obtain a modern multireference approach.
An implementation of the single-site density matrix renormalization group (DMRG) approach for the quantum-chemical application to calculate approximate as well as exact electronic ground-state energies and wave functions of atoms and molecules is presented.The algorithm is formulated in terms of the variational and sizeconsistent class of wave functions, matrix product states, and their operator analogues, matrix product operators.A procedure to construct quantum-chemical Hamiltonian operators with long-range Coulomb interactions as matrix product operators is developed.The many-electron Hamiltonian operators are represented as sums and products of one-electron operators.For the efficient matrix product operator construction of quantum-chemical Hamiltonians, different basis spaces including discrete position and momentum space as well as molecular and atomic orbitals are investigated.Discrete position space yields the best matrix product operator size scaling, but introduces artificial long-range couplings in the kinetic energy term of the Hamiltonian which are a consequence of the DMRG algorithm in terms of matrix product states.The huge amount of discrete space points needed for an accurate description poses a problem because of the electron-electron interaction.Therefore, a momentum space or Fourier series expansion of the electron-electron interaction is shown which makes it possible to compress its matrix product operator representation.DMRG applied to quantum-chemical Hamiltonians in matrix product operator form can serve as means to directly (without mean-field treatment in advance) investigate quantum-chemical problems and find the analytically exact ground-state solution of small atoms and molecules.The next basis space analyzed in detail is the atomic orbital basis.The matrix product operator construction and the direct application of the DMRG method in this basis is examined.A physical compression scheme based on a singular value decomposition of the electron-electron interaction in the atomic orbital basis is demonstrated.Furthermore, DMRG calculations in the atomic orbital basis for certain systems where mean-field theory yields qualitatively wrong results and where higher-order approximative correction methods like truncated configuration interaction and coupled cluster cannot fully capture strong electron correlation effects are compared to results obtained in the molecular orbital basis.The comparison shows that atomic orbitals describe the physical situation well and thus lead to a tremendous reduction of the matrix product state dimension and the computational effort.i
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 present the density matrix renormalization group self-consistent field (DMRG-SCF) approach that is analogous to the complete active space self-consisted field (CASSCF) method but instead of using for the description of the active space the full configuration interaction (FCI) method, the DMRG-SCF uses the density matrix renormalization group (DMRG) method. The DMRG-SCF approach, similarly to CASSCF, properly describes the multiconfigurational character of the wave function but avoids the exponential scaling of the FCI method and replaces it with a polynomial scaling. Hence, calculations for a larger number of orbitals and electrons in the active space are possible since the DMRG method provides an efficient tool to automatically select from the full Hilbert space the many-body contracted basis states that are the most important for the description of the wave function.
A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.
During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. The underlying matrix product state (MPS) ansatz is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS controls the size of the corner of the many-body Hilbert space that can be reached. Whereas the MPS ansatz will only yield an efficient description for noncritical one-dimensional systems, it can still be used as a variational ansatz for other finite-size systems. Rather large virtual dimensions are then required. The two most important aspects to reduce the corresponding computational cost are a proper choice and ordering of the active space orbitals, and the exploitation of the symmetry group of the Hamiltonian. By taking care of both aspects, DMRG becomes an efficient replacement for exact diagonalization in quantum chemistry. DMRG and Hartree-Fock theory have an analogous structure. The former can be interpreted as a self-consistent mean-field theory in the DMRG lattice sites, and the latter in the particles. It is possible to build upon this analogy to introduce post-DMRG methods. Based on an approximate MPS, these methods provide improved ansätze for the ground state, as well as for excitations. Exponentiation of the single-particle (single-site) excitations for a Slater determinant (an MPS with open boundary conditions) leads to the Thouless theorem for Hartree-Fock theory (DMRG), an explicit nonredundant parameterization of the entire manifold of Slater determinants (MPS wavefunctions). This gives rise to the configuration interaction expansion for DMRG. The Hubbard-Stratonovich transformation lies at the basis of auxiliary field quantum Monte Carlo for Slater determinants. An analogous transformation for spin-lattice Hamiltonians allows to formulate a promising variant for MPSs.