All Publications

The entanglement asymmetry measures the extent to which a symmetry is broken within a subsystem of an extended quantum system. Here, we analyse this quantity in Haar random states for arbitrary compact, semi-simple Lie groups, building on and generalising recent results obtained for the U(1) symmetric case. We find that, for any group, the average entanglement asymmetry vanishes in the thermodynamic limit when the subsystem is smaller than its complement. When the subsystem and its complement are of equal size, the entanglement asymmetry jumps to a finite value, indicating a sudden transition of the subsystem from a fully symmetric state to one devoid of any symmetry. For larger subsystem sizes, the entanglement asymmetry displays a logarithmic scaling with a coefficient fixed by the dimension of the group. We also investigate the fluctuations of the entanglement asymmetry, which tend to zero in the thermodynamic limit. We check our findings against exact numerical calculations for the SU(2) and SU(3) groups. We further discuss their implications for the thermalisation of isolated quantum systems and black hole evaporation.

We examine the behavior of the entanglement asymmetry in the ground state of a (1+1)-dimensional conformal field theory with a boundary condition that explicitly breaks a bulk symmetry. Our focus is on the asymmetry of a subsystem A originating from the symmetry-breaking boundary and extending into a semi-infinite bulk. By employing the twist field formalism, we derive a universal expression for the asymmetry, showing that the asymptotic behavior for large subsystems is approached algebraically, with an exponent which is twice the conformal dimension of a boundary condition-changing operator. As a secondary result, we also establish a similar asymptotic behavior for the string order parameter. Our exact analytical findings are validated through numerical simulations in the critical Ising and 3-state Potts models.

Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.

The entanglement Hamiltonian (EH) provides the most comprehensive characterization of bipartite entanglement in many-body quantum systems. Ground states of local Hamiltonians inherit this locality, resulting in EHs that are dominated by local, few-body terms. Unfortunately, in nonequilibrium situations, analytic results are rare and largely confined to continuous field theories, which fail to accurately describe microscopic models. To address this gap, we present an analytic result for the EH following a quantum quench in noninteracting fermionic models, valid in the ballistic scaling regime. The derivation adapts the celebrated quasiparticle picture to operators, providing detailed insights into its physical properties. The resulting analytic formula serves as a foundation for engineering EHs in quantum optics experiments.

We study the nonequilibrium dynamics of bosons in a two-dimensional optical lattice after a sudden quench from the superfluid phase to the free-boson regime. The initial superfluid state is described approximately using both the Bogoliubov theory and the Gaussian variational principle. The subsequent time evolution remains Gaussian, and we compare the results from each approximation of the initial state by examining different aspects of the dynamics. First, we analyze the entanglement entropy and observe that, in both cases, it increases linearly with time before reaching a saturation point. This behavior is attributed to the propagation of entangled pairs of quantum depletions in the superfluid state. Next, we explore the fate of particle-number symmetry, which is spontaneously broken in the superfluid phase. To do so, we use the entanglement asymmetry, a recently introduced observable that enables us to track symmetry breaking within a subsystem. We observe that its evolution varies qualitatively depending on the theory used to describe the initial state. However, in both cases, the symmetry remains broken and is never restored in the stationary state. Finally, we assess the time it takes to reach the stationary state by evaluating the quantum fidelity between the stationary reduced density matrix and the time-evolved one. Interestingly, within the Gaussian variational principle, we find that an initial state further from the stationary state can relax more quickly than one closer to it, indicating the presence of the recently discovered quantum Mpemba effect. We derive the microscopic conditions necessary for this effect to occur and demonstrate that these conditions are never met in the Bogoliubov theory.

We employ the quasiparticle picture of entanglement evolution to obtain an effective description for the out-of-equilibrium entanglement Hamiltonian at the hydrodynamical scale following quantum quenches in free fermionic systems in two or more spatial dimensions. Specifically, we begin by applying dimensional reduction techniques in cases where the geometry permits, building directly on established results from one-dimensional systems. Subsequently, we generalize the analysis to encompass a wider range of geometries. We obtain analytical expressions for the entanglement Hamiltonian valid at the ballistic scale, which reproduce the known quasiparticle picture predictions for the Renyi entropies and full counting statistics. We also numerically validate the results with excellent precision by considering quantum quenches from several initial configurations.

We introduce a general approximate method for calculating the one-body correlations and the momentum distributions of one-dimensional Bose gases at finite interaction strengths and temperatures trapped in smooth confining potentials. Our method combines asymptotic techniques for the long-distance behavior of the gas (similar to Luttinger liquid theory) with known short-distance expansions. We derive analytical results for the limiting cases of strong and weak interactions and provide a general procedure for calculating one-body correlations at any interaction strength. A step-by-step explanation of the numerical method used to compute Green's functions (needed as input to our theory) is included. We benchmark our method against exact numerical calculations and compare its predictions to recent experimental results.

The interplay between symmetries and entanglement in out-of-equilibrium quantum systems is currently at the center of an intense multidisciplinary research effort. Here we introduce a setting where these questions can be characterized exactly by considering dual-unitary circuits with an arbitrary number of U(1) charges. After providing a complete characterization of these systems we show that one can introduce a class of solvable states, which extends that of generic dual-unitary circuits, for which the nonequilibrium dynamics can be solved exactly. In contrast to the known class of solvable states, which relax to the infinite-temperature state, these states relax to a family of nontrivial generalized Gibbs ensembles. The relaxation process of these states can be simply described by a linear growth of the entanglement entropy followed by saturation to a nonmaximal value but with maximal entanglement velocity. We then move on to consider the dynamics from nonsolvable states, combining the exact results with the entanglement membrane picture we argue that the entanglement dynamics from these states is qualitatively different from that of the solvable ones. It shows two different growth regimes characterized by two distinct slopes, both corresponding to submaximal entanglement velocities. Moreover, we show that nonsolvable initial states can give rise to the quantum Mpemba effect, where less symmetric initial states restore the symmetry faster than more symmetric ones.

We study the explicit breaking of a SU(2) symmetry to a U(1) subgroup employing the entanglement asymmetry, a recently introduced observable that measures how much symmetries are broken in a part of extended quantum systems. We consider as specific model the critical XXZ spin chain, which breaks the SU(2) symmetry of spin rotations except at the isotropic point, and is described by the massless compact boson in the continuum limit. We examine the U(1) subgroup of SU(2) that is broken outside the isotropic point by applying conformal perturbation theory, which we complement with numerical simulations on the lattice. We also analyse the entanglement asymmetry of the full SU(2) group. By relying on very generic scaling arguments, we derive an asymptotic expression for it.

The entanglement asymmetry is an observable independent tool to investigate the relaxation of quantum many-body systems through the restoration of an initially broken symmetry of the dynamics. In this paper we use this to investigate the effects of interactions on quantum relaxation in a paradigmatic integrable model. Specifically, we study the dynamical restoration of the U(1) symmetry corresponding to rotations about the z-axis in the XXZ model quenched from a tilted ferromagnetic state. We find two distinct patterns of behaviour depending upon the interaction regime of the model. In the gapless regime, at roots of unity, we find that the symmetry restoration is predominantly carried out by bound states of spinons of maximal length. The velocity of these bound states is suppressed as the anisotropy is decreased toward the isotropic point leading to slower symmetry restoration. By varying the initial tilt angle, one sees that symmetry restoration is slower for an initially smaller tilt angle, signifying the presence of the quantum Mpemba effect. In the gapped regime, however, spin transport for non maximally tilted states is dominated by smaller bound states with longer bound states becoming frozen. This leads to much longer time scales for restoration compared to the gapless regime. In addition, the quantum Mpemba effect is absent in the gapped regime.

Quasicondensation in one dimension is known to occur for equilibrium systems of hard-core bosons (HCBs) at zero temperature. This phenomenon arises due to the off-diagonal long-range order in the ground state, characterized by a power-law decay of the one-particle density matrix g 1 ( x , y ) ∼ | x − y | − 1 / 2 —a well-known outcome of Luttinger liquid theory. Remarkably, HCBs, when allowed to freely expand from an initial product state (i.e. characterized by initial zero correlation), exhibit quasicondensation and demonstrate the emergence of off-diagonal long-range order during nonequilibrium dynamics. This phenomenon has been substantiated by numerical and experimental investigations in the early 2000s. In this work, we revisit the dynamical quasicondensation of HCBs, providing a fully analytical treatment of the issue. In particular, we derive an exact asymptotic formula for the equal-time one-particle density matrix by borrowing ideas from the framework of quantum Generalized Hydrodynamics. Our findings elucidate the phenomenology of quasicondensation and of dynamical fermionization occurring at different stages of the time evolution, as well as the crossover between the two.

Local relaxation after a quench in 1D quantum many-body systems is a well-known and very active problem with rich phenomenology. Except in pathological cases, the local relaxation is accompanied by the local restoration of the symmetries broken by the initial state that are preserved by unitary evolution. Recently, the entanglement asymmetry has been introduced as a probe to study the interplay between symmetry breaking and relaxation in an extended quantum system. In particular, using the entanglement asymmetry, it has been shown that the more a symmetry is initially broken, the faster it may be restored. This surprising effect, which has also been observed in trapped-ion experiments, can be seen as a quantum version of the Mpemba effect, and is manifested by the crossing at a finite time of the entanglement asymmetry curves of two different initial symmetry-breaking configurations. In this paper we show that, by tuning the initial state, the symmetry dynamics in free fermionic systems can display much richer behavior than seen previously. In particular, for certain classes of initial states, including the ground states of free fermionic models with long-range couplings, the entanglement asymmetry can exhibit multiple crossings. This illustrates that the existence of the quantum Mpemba effect can only be inferred by examining the late-time behavior of the entanglement asymmetry.

Hawking's discovery that black holes can evaporate through radiation emission has posed a number of questions that with time became fundamental hallmarks for a quantum theory of gravity. The most famous one is likely the information paradox, which finds an elegant explanation in the Page argument suggesting that a black hole and its radiation can be effectively represented by a random state of qubits. Leveraging the same assumption, we ponder the extent to which a black hole may display emergent symmetries, employing the entanglement asymmetry as a modern, information-based indicator of symmetry breaking. We find that for a random state devoid of any symmetry, a U(1) symmetry emerges and it is exact in the thermodynamic limit before the Page time. At the Page time, the entanglement asymmetry shows a finite jump to a large value. Our findings imply that the emitted radiation is symmetric up to the Page time and then undergoes a sharp transition. Conversely the black hole is symmetric only after the Page time.

The quantum Mpemba effect is the counterintuitive nonequilibrium phenomenon wherein the dynamic restoration of a broken symmetry occurs more rapidly when the initial state exhibits a higher degree of symmetry breaking. The effect has been recently discovered theoretically and observed experimentally in the framework of global quantum quenches, but so far it has only been investigated in one-dimensional systems. Here we focus on a two-dimensional free-fermion lattice employing the entanglement asymmetry as a measure of symmetry breaking. Our investigation begins with the ground-state analysis of a system featuring nearest-neighbor hoppings and superconducting pairings, the latter breaking explicitly the U(1) particle-number symmetry. We compute analytically the entanglement asymmetry of a periodic strip using dimensional reduction, an approach that allows us to adjust the extent of the transverse size, achieving a smooth crossover between one and two dimensions. Further applying the same method, we study the time evolution of the entanglement asymmetry after a quench to a Hamiltonian with only nearest-neighbor hoppings, preserving the particle-number symmetry which is restored in the stationary state. We find that the quantum Mpemba effect is strongly affected by the size of the system in the transverse dimension, with the potential to either enhance or spoil the phenomenon depending on the initial states. We establish the conditions for its occurrence based on the properties of the initial configurations, extending the criteria found in the one-dimensional case.

The highly complicated nature of far from equilibrium systems can lead to a complete breakdown of the physical intuition developed in equilibrium. A famous example of this is the Mpemba effect, which states that nonequilibrium states may relax faster when they are further from equilibrium or, put another way, hot water can freeze faster than warm water. Despite possessing a storied history, the precise criteria and mechanisms underpinning this phenomenon are still not known. Here, we study a quantum version of the Mpemba effect that takes place in closed many-body systems with a U(1) conserved charge: in certain cases a more asymmetric initial configuration relaxes and restores the symmetry faster than a more symmetric one. In contrast to the classical case, we establish the criteria for this to occur in arbitrary integrable quantum systems using the recently introduced entanglement asymmetry. We describe the quantum Mpemba effect in such systems and relate the properties of the initial state, specifically its charge fluctuations, to the criteria for its occurrence. These criteria are expounded using exact analytic and numerical techniques in several examples, a free fermion model, the Rule 54 cellular automaton, and the Lieb-Liniger model.

The nonequilibrium physics of many-body quantum systems harbors various unconventional phenomena. In this Letter, we experimentally investigate one of the most puzzling of these phenomena - the quantum Mpemba effect, where a tilted ferromagnet restores its symmetry more rapidly when it is farther from the symmetric state compared to when it is closer. We present the first experimental evidence of the occurrence of this effect in a trapped-ion quantum simulator. The symmetry breaking and restoration are monitored through entanglement asymmetry, probed via randomized measurements, and postprocessed using the classical shadows technique. Our findings are further substantiated by measuring the Frobenius distance between the experimental state and the stationary thermal symmetric theoretical state, offering direct evidence of subsystem thermalization.

Entanglement Hamiltonians provide the most comprehensive characterisation of entanglement in extended quantum systems. A key result in unitary quantum field theories is the Bisognano-Wichmann theorem, which establishes the locality of the entanglement Hamiltonian. In this work, our focus is on the non-Hermitian Su-Schrieffer-Heeger (SSH) chain. We study the entanglement Hamiltonian both in a gapped phase and at criticality. In the gapped phase we find that the lattice entanglement Hamiltonian is compatible with a lattice Bisognano-Wichmann result, with an entanglement temperature linear in the lattice index. At the critical point, we identify a new imaginary chemical potential term absent in unitary models. This operator is responsible for the negative entanglement entropy observed in the non-Hermitian SSH chain at criticality.

The entanglement asymmetry is an information based observable that quantifies the degree of symmetry breaking in a region of an extended quantum system. We investigate this measure in the ground state of one dimensional critical systems described by a CFT. Employing the correspondence between global symmetries and defects, the analysis of the entanglement asymmetry can be formulated in terms of partition functions on Riemann surfaces with multiple non-topological defect lines inserted at their branch cuts. For large subsystems, these partition functions are determined by the scaling dimension of the defects. This leads to our first main observation: at criticality, the entanglement asymmetry acquires a subleading contribution scaling as log ℓ/ℓ for large subsystem length ℓ. Then, as an illustrative example, we consider the XY spin chain, which has a critical line described by the massless Majorana fermion theory and explicitly breaks the U(1) symmetry associated with rotations about the z-axis. In this situation the corresponding defect is marginal. Leveraging conformal invariance, we relate the scaling dimension of these defects to the ground state energy of the massless Majorana fermion on a circle with equally-spaced point defects. We exploit this mapping to derive our second main result: the exact expression for the scaling dimension associated with n defects of arbitrary strengths. Our result generalizes a known formula for the n = 1 case derived in several previous works. We then use this exact scaling dimension to derive our third main result: the exact prefactor of the log ℓ/ℓ term in the asymmetry of the critical XY chain.

Conserved-charge densities are very special observables in quantum many-body systems as, by construction, they encode information about the dynamics. Therefore, their evolution is expected to be of much simpler interpretation than that of generic observables and to return universal information on the state of the system at any given time. Here, we study the dynamics of the fluctuations of conserved U(1) charges in systems that are prepared in charge-asymmetric initial states. We characterize the charge fluctuations in a given subsystem using the full-counting statistics of the truncated charge and the quantum entanglement between the subsystem and the rest resolved to the symmetry sectors of the charge. We show that, even though the initial states considered are homogeneous in space, the charge fluctuations generate an effective inhomogeneity due to the charge-asymmetric nature of the initial states. We use this observation to map the problem into that of charge fluctuations on inhomogeneous, charge-symmetric states and treat it using a recently developed space-time duality approach. Specializing the treatment to interacting integrable systems we combine the space-time duality approach with generalized hydrodynamics to find explicit predictions.

The ‘operator entanglement’ of a quantum operator O is a useful indicator of its complexity, and, in one-dimension, of its approximability by matrix product operators. Here we focus on spin chains with a global U(1) conservation law, and on operators O with a well-defined U(1) charge, for which it is possible to resolve the operator entanglement of O according to the U(1) symmetry. We employ the notion of symmetry resolved operator entanglement (SROE) introduced in Rath et al (2023 PRX Quantum 4 010318) and extend the results of the latter paper in several directions. Using a combination of conformal field theory and of exact analytical and numerical calculations in critical free fermionic chains, we study the SROE of the thermal density matrix ρ β = e − β H and of charged local operators evolving in Heisenberg picture O = e i t H O e − i t H . Our main results are: i) the SROE of ρ β obeys the operator area law; ii) for free fermions, local operators in Heisenberg picture can have a SROE that grows logarithmically in time or saturates to a constant value; iii) there is equipartition of the entanglement among all the charge sectors except for a pair of fermionic creation and annihilation operators.