Quantum work statistics, Loschmidt echo and information scrambling

A universal relation is established between the quantum work probability distribution of an isolated driven quantum system and the Loschmidt echo dynamics of a two-mode squeezed state. When the initial density matrix is canonical, the Loschmidt echo of the purified double thermofield state provides a direct measure of information scrambling and can be related to the analytic continuation of the partition function. Information scrambling is then described by the quantum work statistics associated with the time-reversal operation on a single copy, associated with the sudden negation of the system Hamiltonian.


Introduction
Quantum thermodynamics provides a framework to unify quantum theory, statistical mechanics, information theory and thermodynamics 1 .A central object in this field is the notion of work associated with the dynamics of an isolated quantum system.At the quantum level, work becomes a stochastic variable, described by a probability distribution 2 .The analysis of the associated work statistics has guided seminal developments in stochastic thermodynamics and nonequilibrium statistical mechanics.Prominent examples include the Jarzynski equality [3][4][5] and fluctuation theorems 6,7 , with both classical and quantum counterparts 8 .Quantum thermodynamics is strongly tied to quantum dynamics and the notion of reversibility.The sensitivity of a quantum system to external perturbations is often characterized via a Loschmidt echo that measures the extent to which quantum evolution can be reversed upon an imperfect time-reversal operation [9][10][11] .In this context, the generating function of the work probability distribution associated with the driving of a pure energy eigenstate via a quantum quench is known to be identical to the Loschmidt echo of such eigenstate 12 .This observation has greatly facilitated the understanding of quantum work fluctuations in finite-time, nonequilibrium thermodynamics of many-body systems [12][13][14][15] .More recently, it has been suggested that the work probability distribution is related to the out-of-time order correlators (OTOC) 16,17 .First discussed in condensed matter physics 18 and the study of irreversible processes 19 , OTOC are currently under exhaustive investigation to diagnose quantum chaos and scrambling of information in black hole physics 20 .This Report establishes a fundamental connection between work statistics, Loschmidt echo, and information scrambling.This is done by first showing that the work statistics associated with an arbitrary driving protocol of a generic quantum state is equivalent to the Loschmidt echo dynamics of a two-mode squeezed state in an enlarged Hilbert space.For an initial thermal state, the two-mode squeezed state becomes a thermofield double state, as that used to describe eternal black holes 21 .The work probability distribution resulting from a perfect time-reversal operation on a given system is determined by the analytic continuation of the partition function, and then links to information scrambling.More generally, we show that work statistics dictates information scrambling resulting from an arbitrary Loschmidt echo.

Universal relation between quantum work statistics and Loschmidt echo dynamics
Let us consider an isolated quantum system in a Hilbert space H and described by the time-dependent Hermitian Hamiltonian Ĥs = n E s n |n s n s |, with instantaneous eigenstates |n s and eigenenergies E s n .We consider the evolution from time s = 0 to time s = τ of an arbitrary initial state with density matrix ρ, dictated by the evolution operator where T is the time-ordering operator.We include here also the case of sudden quenches, with τ → 0 + such that Û(τ) → 1, but in which the final Hamiltonian, for which we retain the notation Ĥτ , is different from the initial hamiltonian Ĥ0 .
Characterizing the work done during this driving protocol requires two projective energy measurements: one at the initial time s = 0 and another at s = τ 2, 5 .The results of both measurements, respectively E 0 n and E τ m , give the work done as W = E τ m − E 0 n .This two-energy measurement scheme prevents work from being defined as an observable in the quantum world 2 .Even so, it can be understood in terms of a generalized measurement scheme 22,23 .The corresponding work probability distribution can be written as 4,5 p(W) : where p 0 n = n 0 | ρ|n 0 is the probability that the initial state is found in the n-th eigenstate of the initial Hamiltonian, and p τ m|n ≥ 0 is the transition probability from this initial eigenstate to the m-th eigenstate of the final Hamiltonian |m τ , i.e., The integral representation of the delta function in terms of an auxiliary variable t, δ(W , allows us to write the work distribution as the Fourier transform of the characteristic function We show below that a universal relation exists between the work probability distribution p(W) in an arbitrary unitary protocol and the dynamics of a Loschmidt echo, for any initial state, including mixed states, e.g. at finite temperature.Notice that even if the initial state is pure, the first projective energy measurement in the initial eigenbasis generally leads to a (post-measurement) mixed state where P0 n denotes the projector on the n-th energy eigenstate of the initial Hamiltonian at s = 0. Using the explicit definition of the transition probability, Eq. ( 3), the characteristic function (5) can then be written as This form allows us to identify the auxiliary variable t as a second time of evolution, different from s, as was first proposed in Ref. 12 .Notice that the mixed state ρmix is stationary with respect to Ĥ0 , which entails the property χ * (t, τ) = χ(−t, τ).Whenever the initial state is an energy eigenstate of Ĥ0 , | j 0 , the initial density matrix simplifies to ρ = ρmix = | j 0 j 0 | and p 0 n = n 0 | ρ|n 0 = δ n j .In the sudden quench limit τ → 0 + , the transition probability p τ m|n further reduces to | m τ |n 0 | 2 , and the characteristic function becomes recognizable as a Loschmidt amplitude A(t), given by the survival amplitude of the initial eigenstate | j 0 first propagated forward in time with Ĥ0 and subsequently backward in time with a dynamics generated by Ĥτ 12 .Thus we are presented with the Loschmidt amplitude for a quantum quench in the t evolution, with a driving protocol of the form where Θ(t) is the Heaviside function.The survival probability L(t) = |A(t)| 2 is known as a Loschmidt echo, and can be further related to the local density of states 24 .While its exponential decay in time has often been used to characterize chaotic systems 9 , it can occur in simpler, integrable systems 25 .
In fact, the relation between the generating function and the Loschmidt echo amplitude is universal, since it holds for any generic state -pure or not-under any driving protocol.To show this, we proceed in two steps.First, we consider general initial states ρ that give rise to the post-measurement state ρmix .This state ρmix is then purified by embedding it in the extended Hilbert space H L ⊗ H R , with H L = H R = H, e.g.defining the double-copy state Secondly, we introduce an effective single-copy Hamiltonian Ĥeff τ , acting on H, We note that the dynamics associated with Ĥeff τ with an initial pure state has been recently discussed in 26,27 .
As a result, the characteristic function resulting from the evolution Û(τ) in Eq. ( 1) is equal to the Loschmidt echo amplitude A(t) = Ψ 0 |Ψ t of the purified state (10), when one of the copies evolves under the sudden quench Denoting the composition of unitary evolution operators by the echo matrix M(t, τ) ≡ e +it Ĥeff τ e −it Ĥ0 , the characteristic function (7) becomes where we have emphasized that the evolution acts exclusively on one of the copies, e.g.here, the left one.In particular, (13) corresponds to the echo dynamics of a purified initial state that evolves first forward in time under the initial Hamiltonian Ĥ0 and then backward under the effective final Hamiltonian Ĥeff τ , while the second copy is left unchanged.Explicitly, This result shows that the equivalence between the Loschmidt amplitude and the work statistics holds in a universal setting, and extends the correspondance that was first demonstrated for a system prepared in an energy eigenstate of the initial Hamiltonian Ĥ0 undergoing a sudden quench in 12 .
As direct consequence of this equivalence, the short-time decay of the Loschmidt echo is determined by the variance of the quantum work statistics.In particular, the Loschmidt echo is an even function since L(t) = |χ(t, τ)| 2 = χ(t, τ)χ(−t, τ).Under the assumption of analyticity of the generating function of the cumulants C n of the work probability density, we have ln in terms of the variance of the work fluctuations C 2 = W 2 − W 2 ≡ ∆W 2 associated with the quench (12) on the system.This universal short-time asymptotics is reminiscent of the well-known short-time quadratic decay under unitary evolution and that gives rise to the quantum Zeno effect 28 , with a Gaussian width now identified as the work fluctuations.

Quantum work statistics of time reversal and information scrambling
As the connection in Eq. ( 14) is valid for any generic state, it applies in particular to mixed states at finite temperature.
In what follows, we shall consider the initial state to be the canonical thermal state in H, with density matrix where β = (k B T ) −1 and T is the temperature.In this case, ρ = ρmix , and the purification of both is the so-called thermofield double state 29 , where Ĥ0 acts on only one of the copies and the normalization factor Z(β) = Tr H [e −β Ĥ0 ] = n e −βE 0 n is the partition function.We shall focus on the work statistics associated with the implementation of the time-reversal operation in the original system in H.When the system Hamiltonian is time-reversal, the implementation of the time-reversal operation is equivalent to the sudden negation of the Hamiltonian according to the quench In the laboratory, such negation of the Hamiltonian can be implemented making use of a control ancilla C, when the dynamics of the ancilla-system in the Hilbert space C 2 ⊗ H is generated by the Hamiltonian σ z C ⊗ Ĥ0 , where σ z C is the Pauli matrix acting on the ancilla 30,31 .
The generating function of the corresponding work statistics can be explicitly computed to be where we have used the analytically-continued partition function.The equivalence established above gives the Loschmidt echo L(t) = |χ(t, 0 + )| 2 as the fidelity between the initial thermofield double state |Ψ 0 and its time-evolution, also known as the survival probability: This observation leads to a direct connection with information scrambling in black hole physics [32][33][34][35] , where the thermofield double state represents an entangled state of two conformal field theories that is dual to an eternal black hole via the AdS/CFT correspondence 21 .Specifically, an eternal black hole, which amounts to a non-traversable wormhole between two asymptotic regions of spacetime, has been conjectured to be equivalent to a pair of entangled black holes in a disconnected space with common time 36 .In this interpretation, the thermofield double state (17) is not invariant under time-evolution.Further, it has been argued that no system scrambles information faster than a black hole 20,33 , prompting the analysis of various tools to diagnose chaos, including the fidelity decay of a thermofield double state, Eq. ( 20) [37][38][39][40] .Therefore, the work statistics associated with the implementation of the time-reversal operation in a system dictates the dynamics of information scrambling as described by the survival probability of the thermofield double state.Alternatively, information scrambling can also be related to the work probability distribution associated with a quench in the enlarged Hilbert space H ⊗ H described by the Hamiltonian with Ĥ(t) given in Eq. ( 18), as This corresponds to a local time-reversal operation acting on the left copy only, leaving the right one untouched.In fact any independent unquenched evolution for the right copy would give the same result.
The work probability distribution is the same in both interpretations and follows from ( 4) where Â β = Tr H Â ρth denotes the thermal average of Â in the canonical ensemble.Therefore, the work probability distribution for a time-reversal operation is given by the thermal average of the density of states operator, As a result, it follows that the mean work done to reverse the dynamics of one of the entangled copies in the thermofield double state |Ψ 0 ∈ H ⊗ H can be obtained from the mean thermal energy of the canonical ensemble, in H, Conversely, a number of protocols are available to measure p(W) 22,23,41 and have been successfully implemented in the laboratory [42][43][44] .The survival amplitude of the thermofield double state, A(t) = Z (β + i2t) /Z (β), can thus be accessed from a measurement of the work statistics p(W) associated with a time-reversal operation, that via inverse Fourier transform yields A(t) = χ(t, 0 + ) = dW p(W)e +itW .The Loschmidt echo simply follows as These relations result from the implementation of a perfect time-reversal operation described by the quench ( 18) or (21).More generally, one is led to consider an arbitrary quench dynamics on one of the copies, that can accommodate for imperfect time-reversal operations, driven by a quench from Ĥ0 to Ĥτ .The associated Loschmidt echo reads This quantity can generally be extracted from the work distribution function associated with the general quench in (12) via the inverse Fourier transform yielding the identity Equivalently, the work distribution function refers to the driving of the system from Ĥ0 to Ĥτ when the dynamics is described by the time-evolution operator Û(τ).

Quantum work statistics of quantum chaotic systems
We next illustrate the relation between work statistics, Loschmidt echoes and information scrambling in a driven quantum chaotic system.The Hamiltonians we consider are random N × N Hermitian matrices sampled from the Gaussian Orthogonal Ensemble (GOE) 45 , which are invariant under time reversal.A sudden random quench can be implemented by choosing the initial and final Hamiltonians, Ĥ0 = Ĥi and Ĥτ = Ĥ f , from two independent GOE 46 .The corresponding characteristic function averaged over the GOE where σ i = β + it, σ f = −it, is shown in Fig. 1.The average work probability distribution, also presented in Fig. 1, then directly follows from Fourier transformation, Eq. ( 4), and reads In turn, the Loschmidt echo associated with this driven chaotic system follows from the equivalence identified above, Eq. ( 26), and can be evaluated from the characteristic function as Figure 1 shows the characteristic function, Loschmidt echo and work probability distribution for (a) a sudden negation of the Hamiltonian Ĥ f = − Ĥi , and (b) for an arbitrary sudden quench.At infinite temperature, the work probability distribution is symmetric, representing the fact that the initial and final random Hamiltonians are drawn from ensembles with identical distributions.As the temperature is lowered, the initial thermal state samples predominantly the low-energy spectrum of the initial Hamiltonian, increasing the probability for trajectories associated with positive work, as manifested by the shift of p(W) GOE towards the positive real axis.Such a shift has a non-trivial effect on the decay of | χ(t, 0 + ) GOE | 2 whose long time behavior is characterized by higher values for increasing temperature.Considering the identified equivalence with the survival amplitude, this decay can also be interpreted in terms of information scrambling, as explicit from the Loschmidt echo L GOE (t) represented in Fig. 1.The Loschmidt echo associated with the work statistics for a perfect time reversal operation in one of the copies exhibits the following features, common to scrambling dynamics of many-body-chaotic systems [38][39][40] : it reaches a minimum value at a dip, followed by a ramp and a saturation at long values described by the long-time average.By contrast, the work statistics associated with a quench between two independent random Hamiltonians, leads to a Loschmidt echo characterized by an enhanced decay as manifested by lower values of the dip, which is also broadened.In addition, the subsequent dynamics towards the long-time asymptotics no longer exhibits a clear ramp.A comprehensive analysis of work statistics in quantum chaotic systems is presented in 47 .

Conclusion
In summary, we have shown a universal relation between the quantum work statistics and Loschmidt echo under arbitrary dynamics.Specifically, within the two projective energy measurement scheme, we have shown that the generating function of the work probability distribution of an isolated quantum system prepared in a possible mixed state can be interpreted as the Loschmidt echo amplitude of a purified density matrix in an enlarged Hilbert space, for a quench acting on one of the copies.When the initial state is thermal, the Loschmidt echo describes the evolution of a thermofield double state and is ideally suited to assess information scrambling.In particular, the work statistics associated with the time-reversal operation -the sudden negation of the system Hamiltonian -is dictated by the analytic continuation of the partition function, recently proposed to diagnose quantum chaos.As a result, our work establishes a firm connection between the finite-time thermodynamics of closed quantum systems, irreversibility, and information scrambling.

Figure 1 .
Figure 1.Quantum work statistics of time reversal and information scrambling.Absolute value of the characteristic function, χ(t, 0 + ) GOE , generating the Loschmidt echo L(t) GOE and work probability distribution p(W) GOE for a driven chaotic system implemented by sampling the initial and final Hamiltonians from the Gaussian Orthogonal Ensemble, for (a) an exact time reversal of the Hamiltonian and (b) an arbitrary sudden quench.Results are averaged over 5000 realizations of the GOE for N = 20 at infinite temperature (β = 0, red) and β = 1 (blue).