Infinite boundary conditions for response functions and limit cycles in iDMRG, demonstrated for bilinear-biquadratic spin-1 chains

TL;DR

This paper introduces a method using infinite boundary conditions in iDMRG to efficiently compute response functions and limit cycles in bilinear-biquadratic spin-1 chains, reducing computational effort and enhancing resolution.

Contribution

The authors develop a novel approach employing infinite boundary conditions in iDMRG to compute response functions with fewer simulations and analyze limit cycles in the ground state of spin chains.

Findings

Reduced the number of time-evolution runs from ll to 2sqrt{ll} for finite systems.

Demonstrated high-resolution computation of dynamic structure factors in spin-1 chains.

Identified non-trivial limit cycles in iDMRG simulations depending on system parameters.

Abstract

Response functions $\langle A_x(t) B_y(0)\rangle$ for one-dimensional strongly correlated quantum many-body systems can be computed with matrix product state (MPS) techniques. Especially, when one is interested in spectral functions or dynamic structure factors of translation-invariant systems, the response for some range $|x-y|<\ell$ is needed. We demonstrate how the number of required time-evolution runs can be reduced substantially: (a) If finite-system simulations are employed, the number of time-evolution runs can be reduced from $\ell$ to $2\sqrt{\ell}$. (b) To go beyond, one can employ infinite MPS (iMPS) such that two evolution runs suffice. To this purpose, iMPS that are heterogeneous only around the causal cone of the perturbation are evolved in time, i.e., the simulation is done with infinite boundary conditions. Computing overlaps of these states, spatially shifted relative to each other, yields the response functions for all distances $|x-y|$. As a specific application, we compute the dynamic structure factor for ground states of bilinear-biquadratic spin-1 chains with very high resolution and explain the underlying low-energy physics. To determine the initial uniform iMPS for such simulations, infinite-system density-matrix renormalization group (iDMRG) can be employed. We discuss that, depending on the system and chosen bond dimension, iDMRG with a cell size $n_c$ may converge to a non-trivial limit cycle of length $m$. This then corresponds to an iMPS with an enlarged unit cell of size $m n_c$.