Zhendong Li

Matrix product operators, matrix product states, and ab initio density matrix renormalization group algorithms

Current descriptions of the ab initio density matrix renormalization group (DMRG) algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. In this work, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms.

Spin-Projected Matrix Product States: Versatile Tool for Strongly Correlated Systems

We present a new wave function ansatz that combines the strengths of spin projection with the language of matrix product states (MPS) and matrix product operators (MPO) as used in the density matrix renormalization group (DMRG). Specifically, spin-projected matrix product states (SP-MPS) are constructed as [Formula: see text], where [Formula: see text] is the spin projector for total spin S and |ΨMPS(N,M)⟩ is an MPS wave function with a given particle number N and spin projection M. This new ansatz possesses several attractive features: (1) It provides a much simpler route to achieve spin adaptation (i.e., to create eigenfunctions of Ŝ2) compared to explicitly incorporating the non-Abelian SU(2) symmetry into the MPS. In particular, since the underlying state |ΨMPS(N,M)⟩ in the SP-MPS uses only Abelian symmetries, one does not need the singlet embedding scheme for nonsinglet states, as normally employed in spin-adapted DMRG, to achieve a single consistent variationally optimized state. (2) Due to the use of |ΨMPS(N,M)⟩ as its underlying state, the SP-MPS can be closely connected to broken-symmetry mean-field states. This allows one to straightforwardly generate the large number of broken-symmetry guesses needed to explore complex electronic landscapes in magnetic systems. Further, this connection can be exploited in the future development of quantum embedding theories for open-shell systems. (3) The sum of MPOs representation for the Hamiltonian and spin projector [Formula: see text] naturally leads to an embarrassingly parallel algorithm for computing expectation values and optimizing SP-MPS. (4) Optimizing SP-MPS belongs to the variation-after-projection (VAP) class of spin-projected theories. Unlike usual spin-projected theories based on determinants, the SP-MPS ansatz can be made essentially exact simply by increasing the bond dimensions in |ΨMPS(N,M)⟩. Computing excited states is also simple by imposing orthogonality constraints, which are simple to implement with MPS. To illustrate the versatility of SP-MPS, we formulate algorithms for the optimization of ground and excited states, develop perturbation theory based on SP-MPS, and describe how to evaluate spin-independent and spin-dependent properties such as the reduced density matrices. We demonstrate the numerical performance of SP-MPS with applications to several models typical of strong correlation, including the Hubbard model, and [2Fe-2S] and [4Fe-4S] model complexes.

Time-Step Targeting Time-Dependent and Dynamical Density Matrix Renormalization Group Algorithms with ab Initio Hamiltonians

We study the dynamical density matrix renormalization group (DDMRG) and time-dependent density matrix renormalization group (td-DMRG) algorithms in the ab initio context to compute dynamical correlation functions of correlated systems. We analyze the strengths and weaknesses of the two methods in small model problems and propose two simple improved formulations, DDMRG++ and td-DMRG++, that give increased accuracy at the same bond dimension at a nominal increase in cost. We apply DDMRG++ to obtain the oxygen core-excitation energy in the water molecule in a quadruple-zeta quality basis, which allows us to estimate the remaining correlation error in existing coupled cluster results. Further, we use DDMRG++ to compute the local density of states and gaps and td-DMRG++ to compute the complex polarization function, in linear hydrogen chains with up to 50 H atoms, to study metallicity and delocalization as a function of bond length.