Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1310.3340 · Asia Pacific Physics Conference (APPC)

Application of Matrix Product States to the Hubbard Model in One Spatial Dimension

Yukihiro Shimizu, Koji Matsuura, Hikaru Yahagi

Published Mar 19, 2014• 1 citation
arXivHTMLPDF

Abstract

We investigate the application of matrix product states to the Hubbard model in one spatial dimension with both of open and periodic boundary conditions. We develop the variatinal method that the optimization of the variational parameters is carried out locally and sequentially in the framework of matrix product operators (MPO) by including the sign, due to the anti-commutation relation of fermion operators, in the matrix element of MPO. The numerical accuracy of the ground state energy is examined.

DOI
10.7566/jpscp.1.016011

Related papers

A new general approach for obtaining matrix product operators for spin Hamiltonians with periodic boundary conditions

Payal D. Solanki, M. S. Ramkarthik

Abstract This article concerns the matrix product operators for one-dimensional spin Hamiltonians under open and periodic boundary conditions. A novel approach called ‘the dynamics of blocks’ is proposed to obtain the matrix product operators in the case of Hamiltonians with periodic boundary conditions. This approach works universally for any Hamiltonian with two-particle as well as single-particle interactions. The main significance of the proposed approach is that, the dimensions of the matrix product operators obtained are the same for both open and periodic boundary conditions. This method gives us a generic way of constructing MPOs for any Hamiltonian with arbitrary neighboring interactions, immaterial of how the interaction strength falls. A general numerical implementation of the ‘dynamics of blocks’ method is also carried out. The DMRG algorithm was implemented, and the ground state eigenvalue was obtained using these MPOs to numerically validate the method.

Physica Scripta•Jan 24, 2023•0 citations

Efficient MPO Construction for Long-Range Hamiltonians with Periodic Boundary Conditions: Application to Many-Body Dynamics

Xin Wang, Bo Xiong

Matrix product operator (MPO) serves as a fundamental component in tensor network simulations of quantum many-body dynamics. We employ an MPO construction that introduces additional propagation channels to embed both periodic boundary conditions and finite-range couplings directly into an open boundary MPO. We apply this construction within the time-dependent variational principle (TDVP) framework to simulate quench dynamics in a spin-1/2 chain with finite-range interactions, and benchmark the results numerically against the fourth-order Runge-Kutta method, finding excellent agreement for both single-body and two-body observables. The approach offers a practical route for tensor network simulations of many-body dynamics in periodic finite-range systems.

arXiv:2609.21290•Sep 18, 2026•0 citations

Efficient MPS Algorithm for Periodic Boundary Conditions and Applications Section: Solid matter Original Author's Text: English Abstract: We present the implementation of an efficient algorithm for the calculation of the spectrum of one-dimensional quantum systems with periodic boundary conditions. This algorithm is based on a matrix product representation for quantum states (MPS) and a similar representation for Hamiltonians and other operators (MPO). It is significantly more efficient for systems of about 100 sites and more than for small quantum systems. We apply the formalism to calculate the ground state and the first excited state of a spin-1 Heisenberg ring and deduce the size of a Haldane gap. The results are compared to previous high-precision DMRG calculations. Furthermore, we study the spin-1 systems with a biquadratic nearest-neighbor interaction and show the first results of an application to a mesoscopic Hubbard ring of spinless fermions, which carries a persistent current. Key words: matrix product representation for quantum states (MPS) and Hamiltonians (MPO), spin-1 Heisenberg ring, density matrix renormalization group (DMRG). List Issues Configure

Physikalisch-Technische Bundesanstalt (Bundesallee 100, D-38116 Braunschweig, Germany), Michael Weyrauch, Mykhailo V. Rakov

We present the implementation of an efficient algorithm for the calculation of the spectrum of one-dimensional quantum systems with periodic boundary conditions. This algorithm is based on a matrix product representation for quantum states (MPS) and a similar representation for Hamiltonians and other operators (MPO). It is significantly more efficient for systems of about 100 sites and more than for small quantum systems. We apply the formalism to calculate the ground state and the first excited state of a spin-1 Heisenberg ring and deduce the size of a Haldane gap. The results are compared to previous high-precision DMRG calculations. Furthermore, we study the spin-1 systems with a biquadratic nearest-neighbor interaction and show the first results of an application to a mesoscopic Hubbard ring of spinless fermions, which carries a persistent current.

arXiv:1303.1333•Ukrainian Journal of Physics•Jul 1, 2013•12 citations

Variational Optimization of Continuous Matrix Product States

Benoît Tuybens, Jacopo De Nardis, Jutho Haegeman +1 more

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.

arXiv:2006.01801•Physical Review Letters•Jan 10, 2022•17 citations

Efficient matrix-product state method for periodic boundary conditions

Peter Pippan, Steven R. White, Hans Gerd Evertz

We introduce an efficient method to calculate the ground state of one-dimensional lattice models with periodic boundary conditions. The method works in the representation of matrix product states (MPS), related to the density matrix renormalization group (DMRG) method. It improves on a previous approach by Verstraete et al. We introduce a factorization procedure for long products of MPS matrices, which reduces the computational effort from ${m}^{5}$ to ${m}^{3}$, where $m$ is the matrix dimension, and $m\ensuremath{\simeq}100--1000$ in typical cases. We test the method on the $S=\frac{1}{2}$ and $S=1$ Heisenberg chains. It is also applicable to nontranslationally invariant cases. The method makes ground-state calculations with periodic boundary conditions about as efficient as traditional DMRG calculations for systems with open boundaries.

arXiv:0801.1947•Physical Review B•Feb 12, 2010•88 citations

Matrix Product Operator Constructions for Gauge Theories in the Thermodynamic Limit

Godfrey, Nicholas, Ian P. McCulloch

We present a general method for simulating lattice gauge theories in low dimensions using infinite matrix product states (iMPS). A central challenge in Hamiltonian formulations of gauge theories is the unbounded local Hilbert space associated with gauge degrees of freedom. In one spatial dimension, Gauss's law permits these gauge fields to be integrated out, yielding an effective Hamiltonian with long-range interactions among matter fields. We construct efficient matrix product operator (MPO) representations of these Hamiltonians directly in the thermodynamic limit. Our formulation naturally includes background fields and $θ$-terms, requiring no modifications to the standard iDMRG algorithm. This provides a broadly applicable framework for 1+1D gauge theories and can be extended to quasi-two-dimensional geometries such as infinite cylinders, where tensor-network methods remain tractable. As a benchmark, we apply our construction to the Schwinger model, reproducing expected features including confinement, string breaking, and the critical behavior at finite mass. Because the method alters only the MPO structure, it can be incorporated with little effort into a wide range of iMPS and infinite-boundary-condition algorithms, opening the way to efficient studies of both equilibrium and non-equilibrium gauge dynamics.

arXiv:2508.14145•Aug 19, 2025•0 citations
We investigate the application of matrix product states to the Hubbard model in one spatial dimension with both of open and periodic boundary conditions. We develop the variatinal method that the optimization of the variational parameters is carried out locally and sequentially in the framework of matrix product operators (MPO) by including the sign, due to the anti-commutation relation of fermion operators, in the matrix element of MPO. The numerical accuracy of the ground state energy is examined.
Abstract This article concerns the matrix product operators for one-dimensional spin Hamiltonians under open and periodic boundary conditions. A novel approach called ‘the dynamics of blocks’ is proposed to obtain the matrix product operators in the case of Hamiltonians with periodic boundary conditions. This approach works universally for any Hamiltonian with two-particle as well as single-particle interactions. The main significance of the proposed approach is that, the dimensions of the matrix product operators obtained are the same for both open and periodic boundary conditions. This method gives us a generic way of constructing MPOs for any Hamiltonian with arbitrary neighboring interactions, immaterial of how the interaction strength falls. A general numerical implementation of the ‘dynamics of blocks’ method is also carried out. The DMRG algorithm was implemented, and the ground state eigenvalue was obtained using these MPOs to numerically validate the method.
Matrix product operator (MPO) serves as a fundamental component in tensor network simulations of quantum many-body dynamics. We employ an MPO construction that introduces additional propagation channels to embed both periodic boundary conditions and finite-range couplings directly into an open boundary MPO. We apply this construction within the time-dependent variational principle (TDVP) framework to simulate quench dynamics in a spin-1/2 chain with finite-range interactions, and benchmark the results numerically against the fourth-order Runge-Kutta method, finding excellent agreement for both single-body and two-body observables. The approach offers a practical route for tensor network simulations of many-body dynamics in periodic finite-range systems.
Efficient MPS Algorithm for Periodic Boundary Conditions and Applications Section: Solid matter Original Author's Text: English Abstract: We present the implementation of an efficient algorithm for the calculation of the spectrum of one-dimensional quantum systems with periodic boundary conditions. This algorithm is based on a matrix product representation for quantum states (MPS) and a similar representation for Hamiltonians and other operators (MPO). It is significantly more efficient for systems of about 100 sites and more than for small quantum systems. We apply the formalism to calculate the ground state and the first excited state of a spin-1 Heisenberg ring and deduce the size of a Haldane gap. The results are compared to previous high-precision DMRG calculations. Furthermore, we study the spin-1 systems with a biquadratic nearest-neighbor interaction and show the first results of an application to a mesoscopic Hubbard ring of spinless fermions, which carries a persistent current. Key words: matrix product representation for quantum states (MPS) and Hamiltonians (MPO), spin-1 Heisenberg ring, density matrix renormalization group (DMRG). List Issues Configure
We present the implementation of an efficient algorithm for the calculation of the spectrum of one-dimensional quantum systems with periodic boundary conditions. This algorithm is based on a matrix product representation for quantum states (MPS) and a similar representation for Hamiltonians and other operators (MPO). It is significantly more efficient for systems of about 100 sites and more than for small quantum systems. We apply the formalism to calculate the ground state and the first excited state of a spin-1 Heisenberg ring and deduce the size of a Haldane gap. The results are compared to previous high-precision DMRG calculations. Furthermore, we study the spin-1 systems with a biquadratic nearest-neighbor interaction and show the first results of an application to a mesoscopic Hubbard ring of spinless fermions, which carries a persistent current.
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.
We introduce an efficient method to calculate the ground state of one-dimensional lattice models with periodic boundary conditions. The method works in the representation of matrix product states (MPS), related to the density matrix renormalization group (DMRG) method. It improves on a previous approach by Verstraete et al. We introduce a factorization procedure for long products of MPS matrices, which reduces the computational effort from m5{m}^{5}m5 to m3{m}^{3}m3, where mmm is the matrix dimension, and $m\ensuremath{\simeq}100--1000$ in typical cases. We test the method on the S=12S=\frac{1}{2}S=21​ and S=1S=1S=1 Heisenberg chains. It is also applicable to nontranslationally invariant cases. The method makes ground-state calculations with periodic boundary conditions about as efficient as traditional DMRG calculations for systems with open boundaries.
We present a general method for simulating lattice gauge theories in low dimensions using infinite matrix product states (iMPS). A central challenge in Hamiltonian formulations of gauge theories is the unbounded local Hilbert space associated with gauge degrees of freedom. In one spatial dimension, Gauss's law permits these gauge fields to be integrated out, yielding an effective Hamiltonian with long-range interactions among matter fields. We construct efficient matrix product operator (MPO) representations of these Hamiltonians directly in the thermodynamic limit. Our formulation naturally includes background fields and θθθ-terms, requiring no modifications to the standard iDMRG algorithm. This provides a broadly applicable framework for 1+1D gauge theories and can be extended to quasi-two-dimensional geometries such as infinite cylinders, where tensor-network methods remain tractable. As a benchmark, we apply our construction to the Schwinger model, reproducing expected features including confinement, string breaking, and the critical behavior at finite mass. Because the method alters only the MPO structure, it can be incorporated with little effort into a wide range of iMPS and infinite-boundary-condition algorithms, opening the way to efficient studies of both equilibrium and non-equilibrium gauge dynamics.