Variational optimization of continuous matrix product states

TL;DR

This paper introduces a novel piecewise linear parameterization method for optimizing continuous matrix product states (cMPS), enabling efficient and exact variational calculations for inhomogeneous quantum field systems.

Contribution

It develops a new parameterization technique that allows exact density matrix calculations and efficient optimization of cMPS in systems with external potentials.

Findings

Successfully finds ground states of interacting bosons in external potentials.

Calculates boundary and Casimir energy corrections in continuous many-body systems.

Achieves a computational cost scaling as the cube of the bond dimension.

Abstract

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.