The Dicke model as the contraction limit of a pseudo-deformed Richardson-Gaudin model

The Dicke model is derived in the contraction limit of a pseudo-deformation of the quasispin algebra in the su(2)-based Richardson-Gaudin models. Likewise, the integrability of the Dicke model is established by constructing the full set of conserved charges, the form of the Bethe Ansatz state, and the associated Richardson-Gaudin equations. Thanks to the formulation in terms of the pseudo-deformation, the connection from the su(2)-based Richardson-Gaudin model towards the Dicke model can be performed adiabatically.

Linear Response Theory for the Density Matrix Renormalization Group: Efficient Algorithms for Strongly Correlated Excited States

Linear response theory for the density matrix renormalization group (DMRG-LRT) was first presented in terms of the DMRG renormalization projectors [Dorando et al., J. Chem. Phys. 130, 184111 (2009)]. Later, with an understanding of the manifold structure of the matrix product state (MPS) ansatz, which lies at the basis of the DMRG algorithm, a way was found to construct the linear response space for general choices of the MPS gauge in terms of the tangent space vectors [Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)]. These two developments led to the formulation of the Tamm-Dancoff and random phase approximations (TDA and RPA) for MPS. This work describes how these LRTs may be efficiently implemented through minor modifications of the DMRG sweep algorithm, at a computational cost which scales the same as the ground-state DMRG algorithm. In fact, the mixed canonical MPS form implicit to the DMRG sweep is essential for efficient implementation of the RPA, due to the structure of the second-order tangent space. We present ab initio DMRG-TDA results for excited states of polyenes, the water molecule, and a [2Fe-2S] iron-sulfur cluster.

The Thouless theorem for matrix product states and subsequent post-density matrix renormalization group methods

The similarities between Hartree-Fock (HF) theory and the density-matrix renormalization group (DMRG) are explored. Both methods can be formulated as the variational optimization of a wave-function ansatz. Linearization of the time-dependent variational principle near a variational minimum allows to derive the random phase approximation (RPA). We show that the non-redundant parametrization of the matrix product state (MPS) tangent space [J. Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)] leads to the Thouless theorem for MPS, i.e. an explicit non-redundant parametrization of the entire MPS manifold, starting from a specific MPS reference. Excitation operators are identified, which extends the analogy between HF and DMRG to the Tamm-Dancoff approximation (TDA), the configuration interaction (CI) expansion, and coupled cluster theory. For a small one-dimensional Hubbard chain, we use a CI-MPS ansatz with single and double excitations to improve on the ground state and to calculate low-lying excitation energies. For a symmetry-broken ground state of this model, we show that RPA-MPS allows to retrieve the Goldstone mode. We also discuss calculations of the RPA-MPS correlation energy. With the long-range quantum chemical Pariser-Parr-Pople Hamiltonian, low-lying TDA-MPS and RPA-MPS excitation energies for polyenes are obtained.