Search

Unless constrained by symmetry, measurement of an observable on an ensemble of identical quantum systems returns a distribution of values which are encoded in the full counting statistics. While the mean value of this distribution is important for determining certain properties of a system, the full distribution can also exhibit universal behavior. In this paper we study the full counting statistics of particle number in one dimensional interacting Bose and Fermi gases which have been quenched far from equilibrium. In particular we consider the time evolution of the Lieb-Liniger and Gaudin-Yang models quenched from a Bose-Einstein condensate initial state and calculate the full counting statistics of the particle number within a subsystem. We show that the scaled cumulants of the charge in the initial state and at long times are simply related and in particular the latter are independent of the model parameters. Using the quasi-particle picture we obtain the full time evolution of the cumulants and find that although their endpoints are fixed, the finite time dynamics depends strongly on the model parameters. We go on to construct the scaled cumulant generating functions and from this determine the limiting charge probability distributions at long time which are shown to exhibit distinct non-trivial and non-Gaussian fluctuations and large deviations., Comment: 19 pages, 8 figures, 1 appendix

We study the R\'enyi entanglement entropies along the massless renormalisation group flow that connects the tricritical and critical Ising field theories. Similarly to the massive integrable field theories, we derive a set of bootstrap equations, from which we can analytically calculate the twist field form factors in a recursive way. Additionally, we also obtain them as a non-trivial roaming limit of the sinh-Gordon theory. Then the R\'enyi entanglement entropies are obtained as expansions in terms of the form factors of these branch point twist fields. We find that the form factor expansion of the entanglement entropy along the flow organises in two different kind of terms. Those that couple particles with the same chirality, and reproduce the entropy of the infrared Ising theory, and those that couple particles with different chirality, which provide the ultraviolet contributions. The massless flow under study possesses a global $\mathbb{Z}_2$ spin-flip symmetry. We further consider the composite twist fields associated to this group, which enter in the study of the symmetry resolution of the entanglement. We derive analytical expressions for their form factors both from the bootstrap equations and from the roaming limit of the sinh-Gordon theory., Comment: 55 pages, 3 figures

Topological excitations or defects such as solitons are ubiquitous throughout physics, supporting numerous interesting phenomena like zero energy modes with exotic statistics and fractionalized charges. In this paper, we study such objects through the lens of symmetry-resolved entanglement entropy. Specifically, we compute the charge-resolved entanglement entropy for a single interval in the low-lying states of the Su-Schrieffer-Heeger model in the presence of topological defects. Using a combination of exact and asymptotic analytical techniques, backed up by numerical analysis, we find that, compared to the unresolved counterpart and to the pure system, a richer structure of entanglement emerges. This includes a redistribution between its configurational and fluctuational parts due to the presence of the defect and an interesting interplay with entanglement equipartition. In particular, in a subsystem that excludes the defect, equipartition is restricted to charge sectors of the same parity, while full equipartition is restored if the subsystem includes the defect, as long as the associated zero mode remains unoccupied. Additionally, by exciting zero modes in the presence of multiple defects, we observe a significant enhancement of entanglement in certain charge sectors, due to charge splitting on the defects. The two different scenarios featuring the breakdown of entanglement equipartition are underlied by a joint mechanism, which we unveil by relating them to degeneracies in the spectrum of the entanglement Hamiltonian. In addition, equipartition is shown to stem from an equidistant entanglement spectrum., Comment: 17 pages, 7 figures and 4 appendices

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(\theta)$ 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.

We investigate the non-equilibrium dynamics of the symmetry-resolved R\'enyi entropies in a one-dimensional gas of non-interacting spinless fermions by means of quantum generalised hydrodynamics, which recently allowed to obtain very accurate results for the total entanglement in inhomogeneous quench settings. Although our discussion is valid for any quench setting accessible with quantum generalised hydrodynamics, we focus on the case of a quantum gas initially prepared in a bipartite fashion and subsequently let evolve unitarily with a hopping Hamiltonian. For this system, we characterise the symmetry-resolved R\'enyi entropies as function of time $t$ and of the entangling position $x$ along the inhomogeneous profile. We observe an asymptotic logarithmic growth of the charged moments at half system and an asymptotic restoration of equipartition of entropy among symmetry sectors with deviations which are proportional to the square of the inverse of the total entropy., Comment: 27 pages, 10 figures; v2:minor changes

In this paper, we apply the form factor bootstrap approach to branch point twist fields in the $q$-state Potts model for $q\leq 3$. For $q=3$ this is an integrable interacting quantum field theory with an internal discrete $\mathbb{Z}_3$ symmetry and therefore provides an ideal starting point for the investigation of the symmetry resolved entanglement entropies. However, more generally, for $q\leq 3$ the standard R\'enyi and entanglement entropies are also accessible through the bootstrap programme. In our work we present form factor solutions both for the standard branch point twist field with $q\leq 3$ and for the composite (or symmetry resolved) branch point twist field with $q=3$. In both cases, the form factor equations are solved for two particles and the solutions are carefully checked via the $\Delta$-sum rule. Using our analytic predictions, we compute the leading finite-size corrections to the entanglement entropy and entanglement equipartition for a single interval in the ground state., Comment: 32 pages, 12 figures

In this paper we continue the programme initiated in Part I, that is the study of entanglement measures in the sine-Gordon model. In both parts, we have focussed on one specific technique, that is the well-known connection between branch point twist field correlators and measures of entanglement in 1+1D integrable quantum field theory. Our papers apply this technique for the first time to a non-diagonal theory with an involved particle spectrum, the sine-Gordon model. In this Part II we focus on a different entanglement measure, the symmetry resolved entanglement, and develop its associated twist field description, exploiting the underlying U(1) symmetry of the theory. In this context, conventional branch point twist fields are no longer the fields required, but instead we must work with one of their composite generalisations, which can be understood as the field resulting from the fusion of a standard branch point twist field and the sine-Gordon exponential field associated with U(1) symmetry. The resulting composite twist field has correlators which as usual admit a form factor expansion. In this paper we write the associated form factor equations and solve them for various examples in the breather sector by using the method of angular quantisation. We show that, in the attractive regime, this is the sector which provides the leading contribution to the symmetry resolved entropies, both Renyi and von Neumann. We compute the latter in the limit of a large region size and show that they satisfy the property of equipartition, that is the leading contribution to the symmetry resolved entanglement is independent of the symmetry sector.

The quantum sine-Gordon model is the simplest massive interacting integrable quantum field theory whose two-particle scattering matrix is generally non-diagonal. As such, it is a model that has been extensively studied, especially in the context of the bootstrap programme. In this paper we compute the form factors of a special local field known as the branch point twist field, whose correlation functions are building blocks for measures of entanglement. We consider the attractive regime where the theory posesses a particle spectrum consisting of a soliton, an antisoliton (of opposite $U(1)$ charges) and several (neutral) breathers. In the breather sector we exploit the fusion procedure to compute form factors of heavier breathers from those of lighter ones. We apply our results to the study of the entanglement dynamics after a small mass quench and for short times. We show that in the presence of two or more breathers the von Neumann and R\'enyi entropies display undamped oscillations in time, whose frequencies are proportional to the even breather masses and whose amplitudes are proportional to the breather's one-particle form factor., Comment: 35 pages; 9 figures; 3 tables

We generalise the form factor bootstrap approach to integrable field theories with U(1) symmetry to derive matrix elements of composite branch-point twist fields associated with symmetry resolved entanglement entropies. The bootstrap equations are solved for the free massive Dirac and complex boson theories, which are the simplest theories with U(1) symmetry. We present the exact and complete solution for the bootstrap, including vacuum expectation values and form factors involving any type and arbitrarily number of particles. The non-trivial solutions are carefully cross-checked by performing various limits and by the application of the Delta-theorem. An alternative and compact determination of the novel form factors is also presented. Based on the form factors of the U(1) composite branch-point twist fields, we re-derive earlier results showing entanglement equipartition for an interval in the ground state of the two models.

We consider the form factor bootstrap approach of integrable field theories to derive matrix elements of composite branch-point twist fields associated with symmetry resolved entanglement entropies. The bootstrap equations are determined in an intuitive way and their solution is presented for the massive Ising field theory and for the genuinely interacting sinh-Gordon model, both possessing a $\mathbb{Z}_{2}$ symmetry. The solutions are carefully cross-checked by performing various limits and by the application of the $\Delta$-theorem. The issue of symmetry resolution for discrete symmetries is also discussed. We show that entanglement equipartition is generically expected and we identify the first subleading term (in the UV cutoff) breaking it. We also present the complete computation of the symmetry resolved von Neumann entropy for an interval in the ground state of the paramagnetic phase of the Ising model. In particular, we compute the universal functions entering in the charged and symmetry resolved entanglement., Comment: 46 pages, no figure

We study Euler scale hydrodynamics of massless integrable quantum field theories interpolating between two non-trivial renormalisation group fixed points after inhomogeneous quantum quenches. Using a partitioning protocol with left and right initial thermal states and the recently developed framework of generalised hydrodynamics, we focus on current and density profiles for the energy and momentum as a function of \xi = x=t, where both x and t are sent to infinity. Studying the first few members of the An and Dn massless flows we carry out a systematic treatment of these series and generalise our results to other unitary massless models. In our analysis we find that the profiles exhibit extended plateaux and that non-trivial bounds exist for the energy and momentum densities and currents in the non-equilibrium stationary state, i.e. when \xi = 0. To quantify the magnitude of currents and densities, dynamical central charges are defined and it is shown that the dynamical central charge for the energy current satisfies a certain monotonicity property. We discuss the connection of the Landauer-Buttiker formalism of transport with our results and show that this picture can account for some of the bounds for the currents and for the monotonicity of the dynamical central charge. These properties are shown to be present not only in massless flows but also in the massive sinh-Gordon model suggesting their general validity and the correctness of the Landauer-Buttiker interpretation of transport in integrable field theories. Our results thus imply the existence of a non-equilibrium c-theorem as well, at least in integrable models. Finally we also study the interesting low energy behaviour of the A2 model that corresponds to the massless flow from the tricritical to the critical Ising field theory., Comment: 32+7 pages, 6+1 figures and 4+1 tables

Spreading on networks is influenced by a number of factors including different parts of the inter-event time distribution (IETD), the topology of the network and non-stationarity. In order to understand the role of these factors we study the SI model on temporal networks with different aggregated topologies and different IETDs. Based on analytic calculations and numerical simulations, we show that if the stationary bursty process is governed by power-law IETD, the spreading can be slowed down or accelerated as compared to a Poisson process; the speed is determined by the short time behaviour, which in our model is controlled by the exponent. We demonstrate that finite, so called "locally tree-like" networks, like the Barab\'asi-Albert networks behave very differently from real tree graphs if the IETD is strongly fat-tailed, as the lack or presence of rare alternative paths modifies the spreading. A further important result is that the non-stationarity of the dynamics has a significant effect on the spreading speed for strongly fat-tailed power-law IETDs, thus bursty processes characterized by small power-law exponents can cause slow spreading in the stationary state but also very rapid spreading heavily depending on the age of the processes.