All Publications

Quantum impurity models (QIMs) are ubiquitous throughout physics. As simplified toy models they provide crucial insights for understanding more complicated strongly correlated systems, while in their own right are accurate descriptions of many experimental platforms. In equilibrium, their physics is well understood and have proven a testing ground for many powerful theoretical tools, both numerical and analytical, in use today. Their nonequilibrium physics is much less studied and understood. However, the recent advancements in nonequilibrium integrable quantum systems through the development of generalized hydrodynamics (GHD) coupled with the fact that many archetypal QIMs are in fact integrable presents an enticing opportunity to enhance our understanding of these systems. We take a step towards this by expanding the framework of GHD to incorporate integrable interacting QIMs. We present a set of Bethe-Boltzmann type equations which incorporate the effects of impurity scattering and discuss the new aspects which include entropy production. These impurity GHD equations are then used to study a bipartioning quench wherein a relevant backscattering impurity is included at the location of the bipartition. The density and current profiles are studied as a function of the impurity strength and expressions for the entanglement entropy and full counting statistics are derived.

Owing to its probabilistic nature, a measurement process in quantum mechanics produces a distribution of possible outcomes. This distribution - or its Fourier transform known as full counting statistics (FCS) - contains much more information than say the mean value of the measured observable, and accessing it is sometimes the only way to obtain relevant information about the system. In fact, the FCS is the limit of an even more general family of observables - the charged moments - that characterize how quantum entanglement is split in different symmetry sectors in the presence of a global symmetry. Here we consider the evolution of the FCS and of the charged moments of a U(1) charge truncated to a finite region after a global quantum quench. For large scales these quantities take a simple large-deviation form, showing two different regimes as functions of time: while for times much larger than the size of the region they approach a stationary value set by the local equilibrium state, for times shorter than region size they show a nontrivial dependence on time. We show that, whenever the initial state is also U(1) symmetric, the leading order in time of FCS and charged moments in the out-of-equilibrium regime can be determined by means of a space-time duality. Namely, it coincides with the stationary value in the system where the roles of time and space are exchanged. We use this observation to find some general properties of FCS and charged moments out of equilibrium, and to derive an exact expression for these quantities in interacting integrable models. We test this expression against exact results in the Rule 54 quantum cellular automaton and exact numerics in the XXZ spin-1/2 chain.

We consider the quantum quench in the XX spin chain starting from a tilted Néel state which explicitly breaks the U(1) symmetry of the post-quench Hamiltonian. Very surprisingly, the U(1) symmetry is not restored at large time because of the activation of a non- Abelian set of charges which all break it. The breaking of the symmetry can be effectively and quantitatively characterised by the recently introduced entanglement asymmetry. By a combination of exact calculations and quasi-particle picture arguments, we are able to exactly describe the behaviour of the asymmetry at any time after the quench. Furthermore we show that the stationary behaviour is completely captured by a non-Abelian generalised Gibbs ensemble. While our computations have been performed for a noninteracting spin chain, we expect similar results to hold for the integrable interacting case as well because of the presence of non-Abelian charges also in that case.

In one-dimensional (1D) quantum gases, the momentum distribution (MD) of the atoms is a standard experimental observable, routinely measured in various experimental setups. The MD is sensitive to correlations, and it is notoriously hard to compute theoretically for large numbers of atoms N, which often prevents direct comparison with experimental data. Here we report significant progress on this problem for the 1D Tonks-Girardeau (TG) gas in the asymptotic limit of large N, at zero temperature and driven out of equilibrium by a quench of the confining potential. We find an exact analytical formula for the one-particle density matrix (ψ †(x)ψ(x′)) of the out-of-equilibrium TG gas in the N→∞ limit, valid on distances |x-x′| much larger than the interparticle distance. By comparing with time-dependent Bose-Fermi mapping numerics, we demonstrate that our analytical formula can be used to compute the out-of-equilibrium MD with great accuracy for a wide range of momenta (except in the tails of the distribution at very large momenta). For a quench from a double-well potential to a single harmonic well, which mimics a "quantum Newton cradle"setup, our method predicts the periodic formation of peculiar, multiply peaked, momentum distributions.

We consider the ground state of two species of one-dimensional critical free theories coupled together via a conformal interface. They have an internal U(1) global symmetry and we investigate the quantum fluctuations of the total charge on one side of the interface, giving analytical predictions for the full counting statistics, the charged moments of the reduced density matrix and the symmetry resolved Rényi entropies. Our approach is based on the relation between the geometry with the defect and the homogeneous one, and it provides a way to characterize the spectral properties of the correlation functions restricted to one of the two species. Our analytical predictions are tested numerically, finding a perfect agreement.

When a free Fermi gas on a lattice is subject to the action of a linear potential it does not drift away, as one would naively expect, but it remains spatially localized. Here we revisit this phenomenon, known as Stark localization, within the recently proposed framework of generalized hydrodynamics. In particular, we consider the dynamics of an initial state in the form of a domain wall and we recover known results for the particle density and the particle current, while we derive analytical predictions for relevant observables such as the entanglement entropy and the full counting statistics. Then, we extend the analysis to generic potentials, highlighting the relationship between the occurrence of localization and the presence of peculiar closed orbits in phase space, arising from the lattice dispersion relation. We also compare our analytical predictions with numerical calculations and with the available results, finding perfect agreement. This approach paves the way for an exact treatment of the interacting case known as Stark many-body localization.

The negativity Hamiltonian, defined as the logarithm of a partially transposed density matrix, provides an operatorial characterisation of mixed-state entanglement. However, so far, it has only been studied for the mixed-state density matrices corresponding to subsystems of globally pure states. Here, we consider as a genuine example of a mixed state the one-dimensional massless Dirac fermions in a system at finite temperature and size. As subsystems, we consider an arbitrary set of disjoint intervals. The structure of the corresponding negativity Hamiltonian resembles the one for the entanglement Hamiltonian in the same geometry: in addition to a local term proportional to the stress-energy tensor, each point is non-locally coupled to an infinite but discrete set of other points. However, when the lengths of the transposed and non-transposed intervals coincide, the structure remarkably simplifies and we retrieve the mild non-locality of the ground state negativity Hamiltonian. We also conjecture an exact expression for the negativity Hamiltonian associated to the twisted partial transpose, which is a Hermitian fermionic matrix. We finally obtain the continuum limit of both the local and bi-local operators from exact numerical computations in free-fermionic chains.

The study of entanglement in the symmetry sectors of a theory has recently attracted a lot of attention since it provides better understanding of some aspects of quantum many-body systems. In this paper, we extend this analysis to the case of non-Hermitian models, in which the reduced density matrix ρA may be nonpositive definite and the entanglement entropy negative or even complex. Here we examine in detail the symmetry-resolved entanglement in the ground state of the non-Hermitian Su-Schrieffer-Heeger chain at the critical point, a model that preserves particle number and whose scaling limit is a bc-ghost nonunitary conformal field theory (CFT). By combining bosonization techniques in the field theory and exact lattice numerical calculations, we analytically derive the charged moments of ρA and |ρA|. From them, we can understand the origin of the nonpositiveness of ρA and naturally define a positive-definite reduced density matrix in each charge sector, which gives a well-defined symmetry-resolved entanglement entropy. As a by-product, we also obtain the analytical distribution of the critical entanglement spectrum.

The effect of unstable quasiparticles in the out-of-equilibrium dynamics of certain integrable systems has been the subject of several recent studies. In this paper we focus on the stationary value of the entanglement entropy density, its growth rate, and related functions, after a quantum quench. We consider several quenches, each of which is characterised by a corresponding squeezed coherent state. In the quench action approach, the coherent state amplitudes K(θ) become input data that fully characterise the large-time stationary state, thus also the corresponding Yang-Yang entropy. We find that, as function of the mass of the unstable particle, the entropy growth rate has a global minimum signalling the depletion of entropy that accompanies a slowdown of stable quasiparticles at the threshold for the formation of an unstable excitation. We also observe a separation of scales governed by the interplay between the mass of the unstable particle and the quench parameter, separating a non-interacting regime described by free fermions from an interacting regime where the unstable particle is present. This separation of scales leads to a double-plateau structure of many functions, where the relative height of the plateaux is related to the ratio of central charges of the UV fixed points associated with the two regimes, in full agreement with conformal field theory predictions. The properties of several other functions of the entropy and its growth rate are also studied in detail, both for fixed quench parameter and varying unstable particle mass and viceversa.

Whenever a system possesses a conserved charge, the density matrix splits into eigenspaces associated to the each symmetry sector and we can access the entanglement entropy in a given subspace, known as symmetry resolved entanglement (SRE). Here, we first evaluate the SRE for massless Dirac fermions in a system at finite temperature and size, i.e. on a torus. Then we add a massive term to the Dirac action and we treat it as a perturbation of the massless theory. The charge-dependent entropies turn out to be equally distributed among all the symmetry sectors at leading order. However, we find subleading corrections which depend both on the mass and on the boundary conditions along the torus. We also study the resolution of the fermionic negativity in terms of the charge imbalance between two subsystems. We show that also for this quantity, the presence of the mass alters the equipartition among the different imbalance sectors at subleading order.

We study the melting of a domain wall in a free-fermionic chain with a localised impurity. We find that the defect enhances quantum correlations in such a way that even the smallest scatterer leads to a linear growth of the entanglement entropy contrasting the logarithmic behaviour in the clean system. Exploiting the hydrodynamic approach and the quasiparticle picture, we provide exact predictions for the evolution of the entanglement entropy for arbitrary bipartitions. In particular, the steady production of pairs at the defect gives rise to non-local correlations among distant points. We also characterise the subleading logarithmic corrections, highlighting some universal features.

A gas of interacting fermions confined in a quasi one-dimensional geometry shows a BEC to BCS crossover upon slowly driving its coupling constant through a confinement-induced resonance. On one side of the crossover the fermions form tightly bound bosonic molecules behaving as a repulsive Bose gas, while on the other they form Cooper pairs, whose size is much larger than the average interparticle distance. Here we consider the situation arising when the coupling constant is varied suddenly from the BEC to the BCS value. Namely, we study a BEC-to-BCS quench. By exploiting a suitable continuum limit of recently discovered solvable quenches in the Hubbard model, we show that the local stationary state reached at large times after the quench can be determined exactly by means of the quench action approach. We provide an experimentally accessible characterization of the stationary state by computing local pair correlation function as well as the quasiparticle distribution functions. We find that the steady state is increasingly dominated by two-particle spin singlet bound states for stronger interaction strength, but that bound state formation is inhibited at larger BEC density. The bound state rapidity distribution displays quartic power-law decay suggesting a violation of Tan's contact relations.

The operator entanglement (OE) is a key quantifier of the complexity of a reduced density matrix. In out-of-equilibrium situations, e.g., after a quantum quench of a product state, it is expected to exhibit an entanglement barrier. The OE of a reduced density matrix initially grows linearly as entanglement builds up between the local degrees of freedom; it then reaches a maximum and ultimately decays to a small finite value as the reduced density matrix converges to a simple stationary state through standard thermalization mechanisms. Here, by performing a new data analysis of the published experimental results of Brydges et al. [Science 364, 260 (2019)], we obtain the first experimental estimation of the OE of a subsystem reduced density matrix in a quantum many-body system. We employ the randomized-measurements toolbox and we introduce and develop a new efficient method to postprocess experimental data in order to extract higher-order density-matrix functionals and access the OE. The OE thus obtained displays the expected barrier as long as the experimental system is large enough. For smaller systems, we observe a barrier with a double-peak structure, the origin of which can be interpreted in terms of pairs of quasiparticles being reflected at the boundary of the qubit chain. As U(1) symmetry plays a key role in our analysis, we introduce the notion of symmetry-resolved operator entanglement (SROE), in addition to the total OE. To gain further insights into the SROE, we provide a thorough theoretical analysis of this new quantity in chains of noninteracting fermions, which, in spite of their simplicity, capture most of the main features of OE and SROE. In particular, we uncover three main physical effects: the presence of a barrier in any charge sector, a time delay for the onset of the growth of SROE, and an effective equipartition between charge sectors.

We study the ground-state entanglement Hamiltonian of several disjoint intervals for the massless Dirac fermion on the half-line. Its structure consists of a local part and a bi-local term that couples each point to another one in each other interval. The bi-local operator can be either diagonal or mixed in the fermionic chiralities and it is sensitive to the boundary conditions. The knowledge of such entanglement Hamiltonian is the starting point to evaluate the negativity Hamiltonian, i.e. the logarithm of the partially transposed reduced density matrix, which is an operatorial characterisation of entanglement of subsystems in mixed states. We find that the negativity Hamiltonian inherits the structure of the corresponding entanglement Hamiltonian. We finally show how the continuum expressions for both these operators can be recovered from exact numerical computations in free-fermion chains.

The study of the entanglement dynamics plays a fundamental role in understanding the behaviour of many-body quantum systems out of equilibrium. In the presence of a globally conserved charge, further insights are provided by the knowledge of the resolution of entanglement in the various symmetry sectors. Here, we carry on the program we initiated in Parez et al (2021 Phys. Rev. B 103 L041104), for the study of the time evolution of the symmetry resolved entanglement in free fermion systems. We complete and extend our derivations also by defining and quantifying a symmetry resolved mutual information. The entanglement entropies display a time delay that depends on the charge sector that we characterise exactly. Both entanglement entropies and mutual information show effective equipartition in the scaling limit of large time and subsystem size. Furthermore, we argue that the behaviour of the charged entropies can be quantitatively understood in the framework of the quasiparticle picture for the spreading of entanglement, and hence we expect that a proper adaptation of our results should apply to a large class of integrable systems. We also find that the number entropy grows logarithmically with time before saturating to a value proportional to the logarithm of the subsystem size.

The study of non-equilibrium dynamics of many-body systems after a quantum quench received a considerable boost and a deep theoretical understanding from the path integral formulation in imaginary time. However, the celebrated problem of a quench in the Luttinger parameter of a one dimensional quantum critical system (massless quench) has so far only been solved in the real-time Heisenberg picture. In order to bridge this theoretical gap and to understand on the same ground massive and massless quenches, we study the problem of a gaussian field characterized by a coupling parameter K within a strip and a different one K0 in the remaining two semi-infinite planes. We give a fully analytical solution using the electrostatic analogy with the problem of a dielectric material within a strip surrounded by an infinite medium of different dielectric constant, and exploiting the method of charge images. After analytic continuation, this solution allows us to obtain all the correlation functions after the quench within a path integral approach in imaginary time, thus recovering and generalizing the results in real time. Furthermore, this imaginary-time approach establishes a remarkable connection between the quench and the famous problem of the conductivity of a Tomonaga-Luttinger liquid coupled to two semi-infinite leads: the two are in fact related by a rotation of the spacetime coordinates.

We consider the ground state of a theory composed by M species of massless complex bosons in one dimension coupled together via a conformal interface. We compute both the Rényi entropy and the negativity of a generic partition of wires, generalizing the approach developed in a recent work for free fermions. These entanglement measures show a logarithmic growth with the system size, and the universal prefactor depends both on the details of the interface and the bipartition. We test our analytical predictions against exact numerical results for the harmonic chain.

We study the quench dynamics of the one-dimensional Hubbard model through the quench action formalism. We introduce a class of integrable initial states—expressed as product states over two sites—for which we can provide an exact characterisation of the late-time regime. This is achieved by finding a closed-form expression for the overlaps between our states and the Bethe ansatz eigenstates, which we check explicitly in the limits of low densities and infinite repulsion. Our solution gives access to the stationary values attained by local observables (we show the explicit example of the density of doubly occupied sites) and the asymptotic entanglement dynamics directly in the thermodynamic limit. Interestingly, we find that for intermediate interaction strength Rényi entropies display a double-slope structure.

We study nonhomogeneous quantum quenches in a one-dimensional gas of repulsive spin-1/2 fermions, as described by the integrable Yang-Gaudin model. By means of generalized hydrodynamics (GHD), we analyze in detail the real-time evolution following a sudden change of the confining potential. We consider in particular release protocols and trap quenches, including a version of the quantum Newton's cradle. At zero temperature, we employ a simplified phase-space hydrodynamic picture to characterize the dynamics of the particle- and spin-density profiles. Away from zero temperatures, we perform a thorough numerical study of the GHD equations, and provide quantitative predictions for different values of the temperature, external magnetic field, and chemical potential. We highlight the qualitative features arising due to the multicomponent nature of the elementary excitations, discussing in particular effects of spin-charge separation and dynamical polarization.

The growth of Renyi entropies after the injection of energy into a correlated system provides a window upon the dynamics of its entanglement properties. We develop here a simulation scheme by which this growth can be determined in Luttinger liquids systems with arbitrary interactions, even those introducing gaps into the liquid. We apply this scheme to an experimentally relevant quench in the sine-Gordon field theory. While for short times we provide analytic expressions for the growth of the second and third Renyi entropy, to access longer times, we combine our scheme with truncated spectrum methods.