A new algorithm allows for the extremely efficient calculation of thermally averaged quantities in one dimension, in conjunction with the density matrix renormalization group method. The key is the judicious selection of a few representative states.
Abstract. – We describe here the extension of the density matrix renormalization group algorithm to the case where the Hamiltonian has a non-Abelian global symmetry group. The block states transform as irreducible representations of the non-Abelian group. Since the representations are multi-dimensional, a single block state in the new representation corresponds to multiple states of the original density matrix renormalization group basis. We demonstrate the usefulness of the construction via the the one-dimensional Hubbard model as the symmetry group is enlarged from U(1) × U(1), up to SU(2) × SU(2). In past years, the density matrix renormalization group (DMRG) method [1] has been extensively used to study one and two dimensional strongly correlated electron systems [2]. This method became very popular when it was realized that it enabled a level of numerical accuracy for one dimensional systems that was not possible using other methods [3]. One major drawback of DMRG is that calculations are performed in a subspace of purely Abelian symmetries, such as the U(1) symmetries of total particle number and the z component of the total spin. Thus one can only obtain a few states in different total particle
Systems of Y-junctions are interesting both from a fundamental viewpoint and because of their potential use in nanoscale devices. These systems can be studied numerically with the density matrix renormalization group (DMRG), but existing algorithms are inefficient. Here, we introduce a much more efficient DMRG algorithm for Y-junction systems. As an example of the use of this method, we study S=21 bound states in Heisenberg S=1 junctions with two geometries, one where the junction consists of a single site, and the other where it consists of a triangle of three sites.
We have applied the momentum space version of the density-matrix renormalization-group method (k-DMRG) in quantum chemistry in order to study the accuracy of the algorithm in this new context. We have shown numerically that it is possible to determine the desired accuracy of the method in advance of the calculations by dynamically controlling the truncation error and the number of block states using a novel protocol that we dubbed dynamical block state selection protocol. The relationship between the real error and truncation error has been studied as a function of the number of orbitals and the fraction of filled orbitals. We have calculated the ground state of the molecules CH2,H2O, and F2 as well as the first excited state of CH2. Our largest calculations were carried out with 57 orbitals, the largest number of block states was 1500--2000, and the largest dimensions of the Hilbert space of the superblock configuration was 800 000--1 200 000.