Extensive v2DM study of the one-dimensional Hubbard model for large lattice sizes: Exploiting translational invariance and parity

Abstract

Using variational density matrix optimization with two- and three-index conditions we study the one-dimensional Hubbard model with periodic boundary conditions at various filling factors. Special attention is directed to the full exploitation of the available symmetries, more specifically the combination of translational invariance and space-inversion parity, which allows for the study of large lattice sizes. We compare the computational scaling of three different semidefinite programming algorithms with increasing lattice size, and find the boundary point method to be the most suited for this type of problem. Several physical properties, such as the two-particle correlation functions, are extracted to check the physical content of the variationally determined density matrix. It is found that the three-index conditions are needed to correctly describe the full phase diagram of the Hubbard model. We also show that even in the case of half filling, where the ground-state energy is close to the exact value, other properties such as the spin-correlation function can be flawed.

28 pages, 10 figures

References in corpus (8)

• The density-matrix renormalization group in the age of matrix product states
• Matrix Product States, Projected Entangled Pair States, and variational renormalization group methods for quantum spin systems
• The one-dimensional Hubbard model: A reminiscence
• Variational determination of the second-order density matrix for the isoelectronic series of beryllium, neon and silicon
• Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations
• Variational two-particle density matrix calculation for the Hubbard model below half filling using spin-adapted lifting conditions
• Subsystem constraints in variational second order density matrix optimization: curing the dissociative behavior
• A primal-dual semidefinite programming algorithm tailored to the variational determination of the two-body density matrix

Cited by in corpus (6)

• An efficient method for strongly correlated electrons in one dimension
• Benchmarking the variational reduced density matrix theory in the doubly-occupied configuration interaction space with integrable pairing models
• Comparison of one-dimensional and quasi-one-dimensional Hubbard models from the variational two-electron reduced-density-matrix method
• Variational optimization of the 2DM: approaching three-index accuracy using extended cluster constraints
• Geometric Constraints on Two-electron Reduced Density Matrices
• Exact ground state for the four-electron problem in a 2D finite honeycomb lattice