The Gutzwiller variational wavefunction (GVW) is commonly employed to capture correlation effects in condensed matter systems such as ferromagnets, ultracold bosonic gases, correlated superconductors, etc. By noticing that the grand-canonical and number-conserving Gutzwiller Ansätze are in fact the zero-order approximation of an expansion in the truncation parameter of a Matrix Product State (MPS), we argue that MPSs, and the algorithms used to operate on them, are not only flexible computational tools but also a unifying theoretical framework that can be used to generalize and improve on the GVW. In fact, we show that a number-conserving GVW is less efficient in capturing the ground state of a quantum system than a more general MPS which can be optimized with comparable computational resources. Moreover, we suggest a corrected time-dependent density matrix renormalization group algorithm that ensures the conservation of the expectation value of the number of particles when a GVW or a MPS are not explicitly number-conserving. The GVW dynamics obtained with our algorithm compares very well with the exact one in 1D. Most importantly, the algorithm works in any dimension for a GVW. We thus expect it to be of great value in the study of the dynamics of correlated quantum systems.
The Gutzwiller variational wavefunction (GVW) is commonly employed to capture correlation effects in condensed matter systems such as ferromagnets, ultracold bosonic gases, correlated superconductors, etc. By noticing that the grand-canonical and number-conserving Gutzwiller Ansätze are in fact the zero-order approximation of an expansion in the truncation parameter of a Matrix Product State (MPS), we argue that MPSs, and the algorithms used to operate on them, are not only flexible computational tools but also a unifying theoretical framework that can be used to generalize and improve on the GVW. In fact, we show that a number-conserving GVW is less efficient in capturing the ground state of a quantum system than a more general MPS which can be optimized with comparable computational resources. Moreover, we suggest a corrected time-dependent density matrix renormalization group algorithm that ensures the conservation of the expectation value of the number of particles when a GVW or a MPS are not explicitly number-conserving. The GVW dynamics obtained with our algorithm compares very well with the exact one in 1D. Most importantly, the algorithm works in any dimension for a GVW. We thus expect it to be of great value in the study of the dynamics of correlated quantum systems.
The truncation or compression of the spectrum of Schmidt values is inherent to the matrix product state (MPS) approximation of one-dimensional quantum ground states. We provide a renormalization group picture by interpreting this compression as an application of Wilson's numerical renormalization group along the imaginary time direction appearing in the path integral representation of the state. The location of the physical index is considered as an impurity in the transfer matrix and static MPS correlation functions are reinterpreted as dynamical impurity correlations. Coarse-graining the transfer matrix is performed using a hybrid variational ansatz based on matrix product operators, combining ideas of MPS and the multiscale entanglement renormalization ansatz. Through numerical comparison with conventional MPS algorithms, we explicitly verify the impurity interpretation of MPS compression, as put forward by V. Zauner et al. [New J. Phys. 17, 053002 (2015)] for the transverse-field Ising model. Additionally, we motivate the conceptual usefulness of endowing MPS with an internal layered structure by studying restricted variational subspaces to describe elementary excitations on top of the ground state, which serves to elucidate a transparent renormalization group structure ingrained in MPS descriptions of ground states.
We present a new wavefunction ansatz that combines the technique 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 |Ψ^(N,S,M)_SP−MPS⟩=P_S|Ψ(N,M)_MPS⟩, where P_S is the spin projector for total spin S and |Ψ(N,M)MPS⟩ is an MPS wavefunction 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 |Ψ(N,M)MPS⟩ in the SP-MPS uses only Abelian symmetries, one does not need the singlet embedding scheme for non-singlet states, as normally employed in spin-adapted DMRG, to achieve a single consistent variationally optimized state. (2) SP-MPS are closely connected to broken symmetry mean-field states. This allows 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 P_S 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 |Ψ(N,M)MPS⟩. 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.
During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. The underlying matrix product state (MPS) ansatz is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS controls the size of the corner of the many-body Hilbert space that can be reached. Whereas the MPS ansatz will only yield an efficient description for noncritical one-dimensional systems, it can still be used as a variational ansatz for other finite-size systems. Rather large virtual dimensions are then required. The two most important aspects to reduce the corresponding computational cost are a proper choice and ordering of the active space orbitals, and the exploitation of the symmetry group of the Hamiltonian. By taking care of both aspects, DMRG becomes an efficient replacement for exact diagonalization in quantum chemistry. DMRG and Hartree-Fock theory have an analogous structure. The former can be interpreted as a self-consistent mean-field theory in the DMRG lattice sites, and the latter in the particles. It is possible to build upon this analogy to introduce post-DMRG methods. Based on an approximate MPS, these methods provide improved ansätze for the ground state, as well as for excitations. Exponentiation of the single-particle (single-site) excitations for a Slater determinant (an MPS with open boundary conditions) leads to the Thouless theorem for Hartree-Fock theory (DMRG), an explicit nonredundant parameterization of the entire manifold of Slater determinants (MPS wavefunctions). This gives rise to the configuration interaction expansion for DMRG. The Hubbard-Stratonovich transformation lies at the basis of auxiliary field quantum Monte Carlo for Slater determinants. An analogous transformation for spin-lattice Hamiltonians allows to formulate a promising variant for MPSs.
Matrix product states (MPSs) illustrate the suitability of tensor networks for the description of interacting many-body systems: ground states of gapped 1D systems are approximable by MPSs, as shown by Hastings [M. B. Hastings, J. Stat. Mech. (2007) P08024]. By contrast, whether MPSs and more general tensor networks can accurately reproduce correlations in critical quantum systems or quantum field theories has not been established rigorously. Ample evidence exists: entropic considerations provide restrictions on the form of suitable ansatz states, and numerical studies show that certain tensor networks can indeed approximate the associated correlation functions. Here, we provide a complete positive answer to this question in the case of MPSs and 2D conformal field theory: we give quantitative estimates for the approximation error when approximating correlation functions by MPSs. Our work is constructive and yields an explicit MPS, thus providing both suitable initial values and a rigorous justification of variational methods.
Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case, however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.
Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.