Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLinear Response Theory for the Density Matrix Renormalization Group: Efficient Algorithms for Strongly Correlated Excited States
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
The purpose of this paper is (i) to present a generic and fully functional implementation of the density-matrix renormalization group (DMRG) algorithm, and (ii) to describe how to write additional strongly-correlated electron models and geometries by using templated classes. Besides considering general models and geometries, the code implements Hamiltonian symmetries in a generic way and parallelization over symmetry-related matrix blocks.
In a Physical Review B paper Chandross and Hicks claim that an analysis of the density-density correlation function in the dimerised Hubbard model of polyacetylene indicates that the optical exciton is bound, and that a previous study by Boman and Bursill that concluded otherwise was incorrect due to numerical innacuracy. We show that the method used in our original paper was numerically sound and well established in the literature. We also show that, when the scaling with lattice size is analysed, the interpretation of the density-density correlation function adopted by Chandross and Hicks in fact implies that the optical exciton is unbound.
We present a linear-response formulation of density cumulant theory (DCT) that provides a balanced and accurate description of many electronic states simultaneously. In the original DCT formulation, only information about a single electronic state (usually, the ground state) is obtained. We discuss the derivation of linear-response DCT, present its implementation for the ODC-12 method (LR-ODC-12), and benchmark its performance for excitation energies in small molecules (N$_2$, CO, HCN, HNC, C$_2$H$_2$, and H$_2$CO), as well as challenging excited states in ethylene, butadiene, and hexatriene. For small molecules, LR-ODC-12 shows smaller mean absolute errors in excitation energies than equation-of-motion coupled cluster theory with single and double excitations (EOM-CCSD), relative to the reference data from EOM-CCSDT. In a study of butadiene and hexatriene, LR-ODC-12 correctly describes the relative energies of the singly-excited $1^1mathrm{B_{u}}$ and the doubly-excited $2^1mathrm{A_{g}}$ states, in excellent agreement with highly accurate semistochastic heat-bath configuration interaction results, while EOM-CCSD overestimates the energy of the $2^1mathrm{A_{g}}$ state by almost 1 eV. Our results demonstrate that linear-response DCT is a promising theoretical approach for excited states of molecules.
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.
A distributed-memory parallelization strategy for the density matrix renormalization group is proposed for cases where correlation functions are required. This new strategy has substantial improvements with respect to previous works. A scalability analysis shows an overall serial fraction of 9.4% and an efficiency of around 60% considering up to eight nodes. Sources of possible parallel slowdown are pointed out and solutions to circumvent these issues are brought forward in order to achieve a better performance.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
