Contributor info: Prof. Fabian Essler

Riccardo Senese, Jacob H. Robertson, Fabian H. L. Essler

SciPost Phys. 17, 139 (2024) · published 20 November 2024 | Toggle abstract · pdf

We consider the probability distributions of the subsystem (staggered) magnetization in ordered and disordered models of quantum magnets in D dimensions. We focus on Heisenberg antiferromagnets and long-range transverse-field Ising models as particular examples. By employing a range of self-consistent time-dependent mean-field approximations in conjunction with Holstein-Primakoff, Dyson-Maleev, Schwinger boson and modified spin-wave theory representations we obtain results in thermal equilibrium as well as during non-equilibrium evolution after quantum quenches. To extract probability distributions we derive a simple formula for the characteristic function of generic quadratic observables in any Gaussian theory of bosons.

Jacob H. Robertson, Riccardo Senese, Fabian H. L. Essler

SciPost Phys. 15, 032 (2023) · published 27 July 2023 | Toggle abstract · pdf

In a recent numerical study by Haldar et al. [Phys. Rev. X 11, 031062] it was shown that signatures of proximate quantum critical points can be observed at early and intermediate times after certain quantum quenches. Said work focused mainly on the case of the axial next-nearest neighbour Ising (ANNNI) model. Here we construct a simple time-dependent mean-field theory that allows us to obtain a quantitatively accurate description of these quenches at short times, which for reasons we explain remains a fair approximation at late times (with some caveats). Our approach provides a simple framework for understanding the reported numerical results as well as fundamental limitations on detecting quantum critical points through quench dynamics. We moreover explain the origin of the peculiar oscillatory behaviour seen in various observables as arising from the formation of a long-lived bound state.

Denis Bernard, Fabian H. L. Essler, Ludwig Hruza, Marko Medenjak

SciPost Phys. 12, 042 (2022) · published 28 January 2022 | Toggle abstract · pdf

We consider the dynamics of fluctuations in the quantum asymmetric simple exclusion process (Q-ASEP) with periodic boundary conditions. The Q-ASEP describes a chain of spinless fermions with random hoppings that are induced by a Markovian environment. We show that fluctuations of the fermionic degrees of freedom obey evolution equations of Lindblad type, and derive the corresponding Lindbladians. We identify the underlying algebraic structure by mapping them to non-Hermitian spin chains and demonstrate that the operator space fragments into exponentially many (in system size) sectors that are invariant under time evolution. At the level of quadratic fluctuations we consider the Lindbladian on the sectors that determine the late time dynamics for the particular case of the quantum symmetric simple exclusion process (Q-SSEP). We show that the corresponding blocks in some cases correspond to known Yang-Baxter integrable models and investigate the level-spacing statistics in others. We carry out a detailed analysis of the steady states and slow modes that govern the late time behaviour and show that the dynamics of fluctuations of observables is described in terms of closed sets of coupled linear differential-difference equations. The behaviour of the solutions to these equations is essentially diffusive but with relevant deviations, that at sufficiently late times and large distances can be described in terms of a continuum scaling limit which we construct. We numerically check the validity of this scaling limit over a significant range of time and space scales. These results are then applied to the study of operator spreading at large scales, focusing on out-of-time ordered correlators and operator entanglement.

Etienne Granet, Henrik Dreyer, Fabian H.L. Essler

SciPost Phys. 12, 019 (2022) · published 13 January 2022 | Toggle abstract · pdf

We consider the XY spin chain with arbitrary time-dependent magnetic field and anisotropy. We argue that a certain subclass of Gaussian states, called Coherent Ensemble (CE) following [1], provides a natural and unified framework for out-of-equilibrium physics in this model. We show that all correlation functions in the CE can be computed using form factor expansion and expressed in terms of Fredholm determinants. In particular, we present exact out-of-equilibrium expressions in the thermodynamic limit for the previously unknown order parameter 1-point function, dynamical 2-point function and equal-time 3-point function.

Etienne Granet, Fabian H. L. Essler

SciPost Phys. 11, 068 (2021) · published 27 September 2021 | Toggle abstract · pdf

We consider the time evolution of local observables after an interaction quench in the repulsive Lieb-Liniger model. The system is initialized in the ground state for vanishing interaction and then time-evolved with the Lieb-Liniger Hamiltonian for large, finite interacting strength $c$. We employ the Quench Action approach to express the full time evolution of local observables in terms of sums over energy eigenstates and then derive the leading terms of a $1/c$ expansion for several one and two-point functions as a function of time $t>0$ after the quantum quench. We observe delicate cancellations of contributions to the spectral sums that depend on the details of the choice of representative state in the Quench Action approach and our final results are independent of this choice. Our results provide a highly non-trivial confirmation of the typicality assumptions underlying the Quench Action approach.

Yuri D. van Nieuwkerk, Jörg Schmiedmayer, Fabian H.L. Essler

SciPost Phys. 10, 090 (2021) · published 26 April 2021 | Toggle abstract · pdf

We consider the non-equilibrium dynamics of a weakly interacting Bose gas tightly confined to a highly elongated double well potential. We use a self-consistent time-dependent Hartree–Fock approximation in combination with a projection of the full three-dimensional theory to several coupled one-dimensional channels. This allows us to model the time-dependent splitting and phase imprinting of a gas initially confined to a single quasi one-dimensional potential well and obtain a microscopic description of the ensuing damped Josephson oscillations.

Etienne Granet, Fabian H. L. Essler

SciPost Phys. 9, 082 (2020) · published 7 December 2020 | Toggle abstract · pdf

We introduce a framework for calculating dynamical correlations in the Lieb-Liniger model in arbitrary energy eigenstates and for all space and time, that combines a Lehmann representation with a $1/c$ expansion. The $n^{\rm th}$ term of the expansion is of order $1/c^n$ and takes into account all $\lfloor \tfrac{n}{2}\rfloor+1$ particle-hole excitations over the averaging eigenstate. Importantly, in contrast to a ``bare" $1/c$ expansion it is uniform in space and time. The framework is based on a method for taking the thermodynamic limit of sums of form factors that exhibit non integrable singularities. We expect our framework to be applicable to any local operator.\\ We determine the first three terms of this expansion and obtain an explicit expression for the density-density dynamical correlations and the dynamical structure factor at order $1/c^2$. We apply these to finite-temperature equilibrium states and non-equilibrium steady states after quantum quenches. We recover predictions of (nonlinear) Luttinger liquid theory and generalized hydrodynamics in the appropriate limits, and are able to compute sub-leading corrections to these.

Etienne Granet, Maurizio Fagotti, Fabian H. L. Essler

SciPost Phys. 9, 033 (2020) · published 4 September 2020 | Toggle abstract · pdf

We consider the problems of calculating the dynamical order parameter two-point function at finite temperatures and the one-point function after a quantum quench in the transverse field Ising chain. Both of these can be expressed in terms of form factor sums in the basis of physical excitations of the model. We develop a general framework for carrying out these sums based on a decomposition of form factors into partial fractions, which leads to a factorization of the multiple sums and permits them to be evaluated asymptotically. This naturally leads to systematic low density expansions. At late times these expansions can be summed to all orders by means of a determinant representation. Our method has a natural generalization to semi-local operators in interacting integrable models.

Yuri D van Nieuwkerk, Fabian H L Essler

SciPost Phys. 9, 025 (2020) · published 25 August 2020 | Toggle abstract · pdf

We consider a model of two tunnel-coupled one-dimensional Bose gases with hard-wall boundary conditions. Bosonizing the model and retaining only the most relevant interactions leads to a decoupled theory consisting of a quantum sine-Gordon model and a free boson, describing respectively the antisymmetric and symmetric combinations of the phase fields. We go beyond this description by retaining the perturbation with the next smallest scaling dimension. This perturbation carries conformal spin and couples the two sectors. We carry out a detailed investigation of the effects of this coupling on the non-equilibrium dynamics of the model. We focus in particular on the role played by spatial inhomogeneities in the initial state in a quantum quench setup.

Aleksandra A. Ziolkowska, Fabian H.L. Essler

SciPost Phys. 8, 044 (2020) · published 19 March 2020 | Toggle abstract · pdf

We consider Lindblad equations for one dimensional fermionic models and quantum spin chains. By employing a (graded) super-operator formalism we identify a number of Lindblad equations than can be mapped onto non-Hermitian interacting Yang-Baxter integrable models. Employing Bethe Ansatz techniques we show that the late-time dynamics of some of these models is diffusive.

Satoshi Ejima, Tomoki Yamaguchi, Fabian H. L. Essler, Florian Lange, Yukinori Ohta, Holger Fehske

SciPost Phys. 5, 059 (2018) · published 7 December 2018 | Toggle abstract · pdf

We consider the dimerized spin-1 $XXZ$ chain with single-ion anisotropy $D$. In absence of an explicit dimerization there are three phases: a large-$D$, an antiferromagnetically ordered and a Haldane phase. This phase structure persists up to a critical dimerization, above which the Haldane phase disappears. We show that for weak dimerization the phases are separated by Gaussian and Ising quantum phase transitions. One of the Ising transitions terminates in a critical point in the universality class of the dilute Ising model. We comment on the relevance of our results to experiments on quasi-one-dimensional anisotropic spin-1 quantum magnets.

Yuri D. van Nieuwkerk, Jörg Schmiedmayer, Fabian H. L. Essler

SciPost Phys. 5, 046 (2018) · published 8 November 2018 | Toggle abstract · pdf

We consider time-of-flight measurements in split one-dimensional Bose gases. It is well known that the low-energy sector of such systems can be described in terms of two compact phase fields $\hat{\phi}_{a,s}(x)$. Building on existing results in the literature we discuss how a single projective measurement of the particle density after trap release is in a certain limit related to the eigenvalues of the vertex operator $e^{i\hat{\phi}_a(x)}$. We emphasize the theoretical assumptions underlying the analysis of "single-shot" interference patterns and show that such measurements give direct access to multi-point correlation functions of $e^{i\hat{\phi}_a(x)}$ in a substantial parameter regime. For experimentally relevant situations, we derive an expression for the measured particle density after trap release in terms of convolutions of the eigenvalues of vertex operators involving both sectors of the two-component Luttinger liquid that describes the low-energy regime of the split condensate. This opens the door to accessing properties of the symmetric sector via an appropriate analysis of existing experimental data.

Stefan Groha, Fabian H. L. Essler, Pasquale Calabrese

SciPost Phys. 4, 043 (2018) · published 28 June 2018 | Toggle abstract · pdf

We consider the full probability distribution for the transverse magnetization of a finite subsystem in the transverse field Ising chain. We derive a determinant representation of the corresponding characteristic function for general Gaussian states. We consider applications to the full counting statistics in the ground state, finite temperature equilibrium states, non-equilibrium steady states and time evolution after global quantum quenches. We derive an analytical expression for the time and subsystem size dependence of the characteristic function at sufficiently late times after a quantum quench. This expression features an interesting multiple light-cone structure.

Lorenzo Piroli, Pasquale Calabrese, Fabian H. L. Essler

SciPost Phys. 1, 001 (2016) · published 14 September 2016 | Toggle abstract · pdf

We study quantum quenches to the one-dimensional Bose gas with attractive interactions in the case when the initial state is an ideal one-dimensional Bose condensate. We focus on properties of the stationary state reached at late times after the quench. This displays a finite density of multi-particle bound states, whose rapidity distribution is determined exactly by means of the quench action method. We discuss the relevance of the multi-particle bound states for the physical properties of the system, computing in particular the stationary value of the local pair correlation function $g_2$.