Spin-adapted density matrix renormalization group algorithms for quantum chemistry.

We extend the spin-adapted density matrix renormalization group (DMRG) algorithm of McCulloch and Gulacsi [Europhys. Lett. 57, 852 (2002)] to quantum chemical Hamiltonians. This involves using a quasi-density matrix, to ensure that the renormalized DMRG states are eigenfunctions of Ŝ(2), and the Wigner-Eckart theorem, to reduce overall storage and computational costs. We argue that the spin-adapted DMRG algorithm is most advantageous for low spin states. Consequently, we also implement a singlet-embedding strategy due to Tatsuaki [Phys. Rev. E 61, 3199 (2000)] where we target high spin states as a component of a larger fictitious singlet system. Finally, we present an efficient algorithm to calculate one- and two-body reduced density matrices from the spin-adapted wavefunctions. We evaluate our developments with benchmark calculations on transition metal system active space models. These include the Fe(2)S(2), [Fe(2)S(2)(SCH(3))(4)](2-), and Cr(2) systems. In the case of Fe(2)S(2), the spin-ladder spacing is on the microHartree scale, and here we show that we can target such very closely spaced states. In [Fe(2)S(2)(SCH(3))(4)](2-), we calculate particle and spin correlation functions, to examine the role of sulfur bridging orbitals in the electronic structure. In Cr(2) we demonstrate that spin-adaptation with the Wigner-Eckart theorem and using singlet embedding can yield up to an order of magnitude increase in computational efficiency. Overall, these calculations demonstrate the potential of using spin-adaptation to extend the range of DMRG calculations in complex transition metal problems.

[1]  Takeshi Yanai,et al.  Second-order perturbation theory with a density matrix renormalization group self-consistent field reference function: theory and application to the study of chromium dimer. , 2011, The Journal of chemical physics.

[2]  Sandeep Sharma,et al.  The density matrix renormalization group in quantum chemistry. , 2011, Annual review of physical chemistry.

[3]  U. Schollwoeck.  The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[4]  Garnet Kin-Lic Chan,et al.  A review of canonical transformation theory , 2010 .

[5]  Garnet Kin-Lic Chan,et al.  Multireference quantum chemistry through a joint density matrix renormalization group and canonical transformation theory. , 2010, The Journal of chemical physics.

[6]  T. Yanai,et al.  High-performance ab initio density matrix renormalization group method: applicability to large-scale multireference problems for metal compounds. , 2009, The Journal of chemical physics.

[7]  Marcel Nooijen,et al.  Obtaining the two-body density matrix in the density matrix renormalization group method. , 2008, The Journal of chemical physics.

[8]  Marcel Nooijen,et al.  On the spin and symmetry adaptation of the density matrix renormalization group method. , 2008, The Journal of chemical physics.

[9]  M. Reiher,et al.  Density matrix renormalization group calculations on relative energies of transition metal complexes and clusters. , 2008, The Journal of chemical physics.

[10]  Debashree Ghosh,et al.  Orbital optimization in the density matrix renormalization group, with applications to polyenes and beta-carotene. , 2007, The Journal of chemical physics.

[11]  M. Reiher,et al.  Decomposition of density matrix renormalization group states into a Slater determinant basis. , 2007, The Journal of chemical physics.

[12]  Markus Reiher,et al.  Construction of environment states in quantum-chemical density-matrix renormalization group calculations. , 2006, The Journal of chemical physics.

[13]  S. White,et al.  Measuring orbital interaction using quantum information theory , 2005, cond-mat/0508524.

[14]  Markus Reiher,et al.  Relativistic DMRG calculations on the curve crossing of cesium hydride. , 2005, The Journal of chemical physics.

[15]  Markus Reiher,et al.  Convergence behavior of the density-matrix renormalization group algorithm for optimized orbital orderings. , 2005, The Journal of chemical physics.

[16]  U. Schollwoeck.  The density-matrix renormalization group , 2004, cond-mat/0409292.

[17]  P. Knowles,et al.  The CIPT2 method: Coupling of multi-reference configuration interaction and multi-reference perturbation theory. Application to the chromium dimer , 2004 .

[18]  Garnet Kin-Lic Chan,et al.  An algorithm for large scale density matrix renormalization group calculations. , 2004, The Journal of chemical physics.

[19]  J. Sólyom,et al.  Optimizing the density-matrix renormalization group method using quantum information entropy , 2003 .

[20]  Guido Fano,et al.  Quantum chemistry using the density matrix renormalization group II , 2003 .

[21]  B. A. Hess,et al.  Controlling the accuracy of the density-matrix renormalization-group method: The dynamical block state selection approach , 2002, cond-mat/0204602.

[22]  B. Roos.  The Ground State Potential for the Chromium Dimer Revisited , 2003 .

[23]  J. Sauer,et al.  Structure and thermochemistry of Fe2S2???/0/+ gas phase clusters and their fragments. B3LYP calculations , 2002 .

[24]  B. A. Hess,et al.  QC-DMRG study of the ionic-neutral curve crossing of LiF , 2002, cond-mat/0208187.

[25]  M. Head‐Gordon,et al.  Highly correlated calculations with a polynomial cost algorithm: A study of the density matrix renormalization group , 2002 .

[26]  I. McCulloch,et al.  The non-Abelian density matrix renormalization group algorithm , 2000, cond-mat/0012319.

[27]  Guido Fano,et al.  Quantum chemistry using the density matrix renormalization group , 2001 .

[28]  I. McCulloch,et al.  Total spin in the density matrix renormalization group algorithm , 2001 .

[29]  I. McCulloch,et al.  Density matrix renormalisation group method and symmetries of the Hamiltonian , 2000 .

[30]  Wada Tatsuaki.  Interaction-round-a-face density-matrix renormalization-group method applied to rotational-invariant quantum spin chains , 2000 .

[31]  I. Ciofini,et al.  Full‐CI quantum chemistry using the density matrix renormalization group , 1999, cond-mat/9912348.

[32]  S. White,et al.  Ab initio quantum chemistry using the density matrix renormalization group , 1998, cond-mat/9808118.

[33]  P. Palmieri,et al.  Epstein-Nesbet second-order perturbation treatment of dynamical electron correlation and ground state potential energy curve of Cr2 , 1997 .

[34]  T. Nishino,et al.  The density matrix renormalization group method applied to interaction round a face Hamiltonians , 1996, cond-mat/9610221.

[35]  I. Bertini,et al.  Individual Reduction Potentials of the Iron Ions in Fe(2)S(2) and High-Potential Fe(4)S(4) Ferredoxins. , 1996, Inorganic chemistry.

[36]  Xiang,et al.  Density-matrix renormalization-group method in momentum space. , 1996, Physical review. B, Condensed matter.

[37]  Östlund,et al.  Thermodynamic limit of density matrix renormalization. , 1995, Physical review letters.

[38]  K. Andersson,et al.  THE CR2 POTENTIAL ENERGY CURVE STUDIED WITH MULTICONFIGURATIONAL SECOND-ORDER PERTURBATION THEORY , 1994 .

[39]  White,et al.  Density-matrix algorithms for quantum renormalization groups. , 1993, Physical review. B, Condensed matter.

[40]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[41]  White,et al.  Density matrix formulation for quantum renormalization groups. , 1992, Physical review letters.

[42]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

[43]  Peter Pulay,et al.  The unrestricted natural orbital–complete active space (UNO–CAS) method: An inexpensive alternative to the complete active space–self‐consistent‐field (CAS–SCF) method , 1989 .

[44]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[45]  B. Roos,et al.  A simple method for the evaluation of the second-order-perturbation energy from external double-excitations with a CASSCF reference wavefunction , 1982 .

[46]  Hans-Joachim Werner,et al.  The self‐consistent electron pairs method for multiconfiguration reference state functions , 1982 .

[47]  B. Brooks,et al.  The graphical unitary group approach to the electron correlation problem. Methods and preliminary applications , 1979 .

[48]  I. Shavitt.  Matrix element evaluation in the unitary group approach to the electron correlation problem , 1978 .

[49]  S. Denmark,et al.  Synthetic Analogs of the Active Sites of Iron-Sulfur Proteins. XI. Synthesis and Properties of Complexes Containing the Fe2S2 Core and the Structures of Bis[o-xylyl-α,α′-dithiolato-μ-sulfido-ferrate(III)] and Bis[p-tolylthiolato-μ-sulfido-ferrate(III)] Dianions , 1975 .

[50]  Klaus Ruedenberg,et al.  EXPECTATION VALUES OF MANY-FERMION SPIN EIGENSTATES. , 1971 .

[51]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .