Ian P. McCulloch

Confinement is an intriguing phenomenon prevalent in condensed matter and high-energy physics. Exploring its effect on the far-from-equilibrium criticality of quantum many-body systems is of great interest both from a fundamental and technological point of view. Here, we employ large-scale uniform matrix product state calculations to show that a qualitative change in the type of dynamical quantum phase transitions (DQPTs) accompanies the confinement-deconfinement transition in three paradigmatic models -- the power-law interacting quantum Ising chain, the two-dimensional quantum Ising model, and the spin-$S$ $\mathrm{U}(1)$ quantum link model. By tuning a confining parameter in these models, it is found that \textit{branch} (\textit{manifold}) DQPTs arise as a signature of (de)confinement. Whereas manifold DQPTs are associated with a sign change of the order parameter, their branch counterparts are not, but rather occur even when the order parameter exhibits considerably constrained dynamics. Our conclusions can be tested in modern quantum-simulation platforms, such as ion-trap setups and cold-atom experiments of gauge theories.

In recent years, the infinite time-evolution block decimation (iTEBD) method has been demonstrated to be one of the most efficient and powerful numerical schemes for time-evolution in one-dimensional quantum many-body systems. However, a major shortcoming of the method, along with other state-of-the-art algorithms for many-body dynamics, has been their restriction to one spatial dimension. We present an algorithm based on a \textit{hybrid} extension of iTEBD where finite blocks of a chain are first locally time-evolved before an iTEBD-like method combines these processes globally. This in turn permits simulating the dynamics of many-body systems in the thermodynamic limit in $d\geq1$ dimensions including in the presence of long-range interactions. Our work paves the way for simulating the dynamics of many-body phenomena that occur exclusively in higher dimensions, and whose numerical treatments have hitherto been limited to exact diagonalization of small systems, which fundamentally limits a proper investigation of dynamical criticality. We expect the algorithm presented here to be of significant importance to validating and guiding investigations in state-of-the-art ion-trap and ultracold-atom experiments.