Entanglement Dynamics after a Quench in Ising Field Theory: A Branch Point Twist Field Approach

We extend the branch point twist field approach for the calculation of entanglement entropies to time-dependent problems in 1+1-dimensional massive quantum field theories. We focus on the simplest example: a mass quench in the Ising field theory from initial mass m0 to final mass m. The main analytical results are obtained from a perturbative expansion of the twist field one-point function in the post-quench quasi-particle basis. The expected linear growth of the Rényi entropies at large times mt ≫ 1 emerges from a perturbative calculation at second order. We also show that the Rényi and von Neumann entropies, in infinite volume, contain subleading oscillatory contributions of frequency 2m and amplitude proportional to (mt)−3/2. The oscillatory terms are correctly predicted by an alternative perturbation series, in the pre-quench quasi-particle basis, which we also discuss. A comparison to lattice numerical calculations carried out on an Ising chain in the scaling limit shows very good agreement with the quantum field theory predictions. We also find evidence of clustering of twist field correlators which implies that the entanglement entropies are proportional to the number of subsystem boundary points

Ballistic transport and boundary resistances in inhomogeneous quantum spin chains

Transport phenomena are central to physics, and transport in the many-body and fully-quantum regime is attracting an increasing amount of attention. It has been recently revealed that some quantum spin chains support ballistic transport of excitations at all energies. However, when joining two semi-infinite ballistic parts, such as the XX and XXZ spin-1/2 models, our understanding suddenly becomes less established. Employing a matrix-product-state ansatz of the wavefunction, we study the relaxation dynamics in this latter case. Here we show that it takes place inside a light cone, within which two qualitatively different regions coexist: an inner one with a strong tendency towards thermalization, and an outer one supporting ballistic transport. We comment on the possibility that even at infinite time the system supports stationary currents and displays a non-zero Kapitza boundary resistance. Our study paves the way to the analysis of the interplay between transport, integrability, and local defects

Hamiltonian truncation approach to quenches in the Ising field theory

In contrast to lattice systems where powerful numerical techniques such as matrix product state based methods are available to study the non-equilibrium dynamics, the non-equilibrium behaviour of continuum systems is much harder to simulate. We demonstrate here that Hamiltonian truncation methods can be efficiently applied to this problem, by studying the quantum quench dynamics of the 1+1 dimensional Ising field theory using a truncated free fermionic space approach. After benchmarking the method with integrable quenches corresponding to changing the mass in a free Majorana fermion field theory, we study the effect of an integrability breaking perturbation by the longitudinal magnetic field. In both the ferromagnetic and paramagnetic phases of the model we find persistent oscillations with frequencies set by the low-lying particle excitations not only for small, but even for moderate size quenches. In the ferromagnetic phase these particles are the various non-perturbative confined bound states of the domain wall excitations, while in the paramagnetic phase the single magnon excitation governs the dynamics, allowing us to capture the time evolution of the magnetisation using a combination of known results from perturbation theory and form factor based methods. We point out that the dominance of low lying excitations allows for the numerical or experimental determination of the mass spectra through the study of the quench dynamics

Statistical localization: From strong fragmentation to strong edge modes

Certain disorder-free Hamiltonians can be nonergodic due to a strong fragmentation of the Hilbert space into disconnected sectors. Here, we characterize such systems by introducing the notion of "statistically localized integrals of motion" (SLIOM), whose eigenvalues label the connected components of the Hilbert space. SLIOMs are not spatially localized in the operator sense, but appear localized to subextensive regions when their expectation value is taken in typical states with a finite density of particles. We illustrate this general concept on several Hamiltonians, both with and without dipole conservation. Furthermore, we demonstrate that there exist perturbations which destroy these integrals of motion in the bulk of the system while keeping them on the boundary. This results in statistically localized strong zero modes, leading to infinitely long-lived edge magnetizations along with a thermalizing bulk, constituting the first example of such strong edge modes in a nonintegrable model. We also show that in a particular example, these edge modes lead to the appearance of topological string order in a certain subset of highly excited eigenstates. Some of our suggested models can be realized in Rydberg quantum simulators