Nonadiabatic energy fluctuations of scale-invariant quantum systems in a time-dependent trap

We consider the nonadiabatic energy fluctuations of a many-body system in a time-dependent harmonic trap. In the presence of scale-invariance, the dynamics becomes self-similar and the nondiabatic energy fluctuations can be found in terms of the initial expectation values of the second moments of the Hamiltonian, square position, and squeezing operators. Nonadiabatic features are expressed in terms of the scaling factor governing the size of the atomic cloud, which can be extracted from time-of-flight images. We apply this exact relation to a number of examples: the single-particle harmonic oscillator, the one-dimensional Calogero-Sutherland model, describing bosons with inverse-square interactions that includes the non-interacting Bose gas and the Tonks-Girdardeau gas as limiting cases, and the unitary Fermi gas. We illustrate these results for various expansion protocols involving sudden quenches of the trap frequency, linear ramps and shortcuts to adiabaticity. Our results pave the way to the experimental study of nonadiabatic energy fluctuations in driven quantum fluids.


Introduction
Nonequilibrium quantum phenomena in many-body systems are often hard to describe by numerical methods due to the growth of entanglement with the time of evolution. In some cases, exact solutions can be found by analytical methods. This is often the case when the dynamics is characterized by the presence of dynamical invariants, also known as invariants of motion. Such conserved quantities can give rise to symmetries manifested in the dynamics. A prominent example is the presence of scale-invariance in the dynamics of a broad family of quantum fluids in time-dependent traps, which is exploited with routine in ultracold atom laboratories [1,2]. In this context, time-of-flight imaging techniques rely on the scale-invariant dynamics of a suddenly released atomic cloud, that expands freely [3][4][5]. Scale invariance is then manifested in local correlation functions, such as the density profile or the density-density correlation function, that become self-similar. One can then relate the correlations measured at a given time after the expansion with those of the initial trapped system.
Scale invariance is not restricted to free dynamics and it is also common in the presence of interactions. Indeed, in harmonically trapped systems it is always realized when the interparticle potential has the same scaling dimension as the kinetic energy. As an example, Bose-Einstein condensates with s-wave contact interactions exhibit it exactly in two-spatial dimensions as a result of the Pitaevskii-Rosch symmetry [6]. Similarly, it is present in the one-dimensional (rational) Calogero-Sutherland gas

Exact Many-Body Dynamics under Scale Invariance
The description of the exact many-body dynamics under scale invariance can be simplified using a sequence of symmetry transformations. To this end, we follow closely Lohe [36] and in particular the formulation for many-body systems in [24,28]. Consider a N-particle quantum system in D spatial dimensions and let us denote by r i and p i the position and momentum of the i-th particle, respectively. We consider a confined many-particle system described by the time-dependent Hamiltonian where the interparticle interaction potential exhibits the scaling property V(λ r) = λ −2 V( r). The dynamical group for the Hamiltonian (1) is SU(1, 1) [45]. Given an arbitrary energy eigenstate at t = 0, the exact dynamics generated by the Hamiltonian (1) for any modulation of the trapping frequency ω(t) is described by a scaling law [15,16,46], where n = {n i : i = 1, . . . , DN} is a shorthand for a complete set of quantum numbers that labels the quantum state as well as the corresponding energy eigenvalue E n (0) = E n [ω(0)]. The scaling factor b = b(t) > 0 as a function of time is the solution of the Ermakov equation with the initial boundary conditions b(0) = 1 andḃ(0) = 0, such that the time-dependent state (2) reduces to the initial eigenstate at the beginning of the process t = 0. The time-dependent state (2) can be obtained by acting with a sequence of unitary transformations on the initial state. We follow Lohe [36] and Jaramillo et al. [24] to describe the time-evolution in terms of the action of two unitary transformations, elements of SU(1, 1). To this end, we introduce We note that T dil implements scaling transformations by acting on a function in the coordinate Similarly, acting on a function in the momentum representation g( p 1 , . . . , p N ), one finds T dil g( p 1 , . . . , p N )T † dil = g( p 1 b, . . . , p N b). By contrast, T r acts as a multiplicative factor in the coordinate representation while it generates translations in momentum space acting on a function g( p 1 , . . . , p N ) as T r g( p 1 , . . . , p N )T † r = g( p 1 − mḃ r 1 /b, . . . , r N − mḃ r N /b).
Denoting T = T r T dil and the Hamiltonian (1) at t = 0 by H 0 , we identify the invariant of motion as an instance of the many-body invariants of motion known under scale-invariant driving [16,17]. Equation (6) is indeed a generalization of the Lewis-Riesenfeld invariant of motion known for a single-particle driven harmonic oscillator [47], satisfies and has time-dependent eigenfunctions Φ n ( r 1 , . . . , r N , t) = TΨ n ( r 1 , . . . , r N , 0) and time-independent eigenvalues E n (0). In terms of the eigenvectors of I the time-evolving state (2) can be written as Ψ n ( r 1 , . . . , r N , t) = e iα t TΨ n ( r 1 , . . . , r N , 0) , hb(t ) 2 dt is a time-dependent dynamical phase.

Exact Nonadiabatic Mean Energy
Knowledge of the exact time-dependent dynamics under scale invariance makes possible to find the exact nonadiabatic mean energy of the system, as shown in [24]. To simplify the notation, we introduce the square position and squeezing operators and note that the time-dependent Hamiltonian (1) can be rewritten as where Using Equaion (8), the expectation value of the time-dependent Hamiltonian in the time-evolving state reads with As a result, the mean nonadiabatic energy is given by As this expression holds for an arbitrary eigenstate, the expectation value can be also carried out with respect to an arbitrary state diagonal in the energy eigenbasis, such as an equilibrium thermal state. Specifically, the nonadiabatic mean-energy following a change of ω(t) is found to be where • ≡ tr (•ρ 0 ) denotes the mean value at the time t = 0 with respect to a state ρ 0 diagonal in the eigenbasis of H(0). Similarly, • t ≡ tr (•ρ t ) denotes the mean value at the time t with respect to the state ρ t = U(t)ρ 0 U(t) † , where U(t) is the time evolution operator. Under scale invariance it is thus possible to relate the evolution of the mean nonadiabatic energy to the initial expectation values of the Hamiltonian, square position and squeezing operators. For an initial state at equilibrium, diagonal in the eigenbasis of H(0), the mean squeezing vanishes C(0) = 0. In addition, the initial expectation value of Q is proportional to the mean energy of the system [24], Using this property and the Ermakov Equation (15), the nonadiabatic mean energy reads where we have introduced the adiabatic mean energy H(t) ad = H(0) ω(t)/ω(0) as well as the nonadiabatic factor Q * (t) that accounts for the amount of energy excitations over the adiabatic dynamics Equations (17) and (18) provide the exact nonadiabatic mean energy of the family of systems governed by Hamiltonian (1). We emphasize that the nonadiabatic factor Q * (t) is simply the ratio between the nonadiabatic and adiabatic mean energies. As it can be expressed in terms of the scaling factor, that governs the size of the atomic cloud, it is an experimentally accessible quantity [23,31].
Summarizing these background results, reported in [24], it suffices to know the evolution of the scaling factor and the initial expectation values at t = 0 of the mean energy, square position and squeezing operators to determine the exact nonadiabatic dynamics of the energy fluctuations. This relation is particularly simple when the initial state is at equilibrium and the mean squeezing vanishes. As we shall see in the next section, a similar exact relation can be derived for the second order of the Hamiltonian required to describe nonadiabatic energy fluctuations, but requires the computation of expectation values at zero time which are model dependent (we have not been able to relate them to properties of the Hamiltonian via symmetries).

Nonadiabatic Energy Fluctuations
We next turn to the focus of our work, the characterization of the energy fluctuations. To this end, we consider the expectation value of the second momentum of the Hamiltonian at time t > 0 in the corresponding time-evolving state We recall that where for compactness we define the time-dependent coefficients We thus find Evaluation of the expectation value of the fourth, fifth and the sixth terms in (22) over an energy eigenstate, or any initial state ρ 0 commuting with H 0 , is straightforward and yields as Q = H 0 /mω 2 0 from equipartition theorem and C = 0 in an equilibrium eigenstate. Although the mean value of the squeezing operator vanishes, the analogue quantity for its square is generally non-zero. In this case, one needs to specify the model to compute explicitly these quantities. Nonadiabatic energy fluctuations are set by the second moment of the Hamiltonian that equals where we remind that • ≡ tr (•ρ 0 ) denotes the mean value at time t = 0. This is the main result of this work. As for the mean energy, the nonadiabatic energy fluctuations can be related to the initial expectation value of the second moment of the Hamiltonian, square position and squeezing operators. Once these quantities are known, nonadiabatic effects are simply encoded in the time-dependent scaling factor b(t), as the coefficients α(t)and β(t) are directly determined by it. At this stage we are ready to consider the adiabatic evolution, corresponding to the limiṫ ω/ω 2 1 [47]. At the level of the Ermakov equation, one can setḃ ≈ 0 andb ≈ 0 and find that the scaling factor in the adiabatic limit is set by b ad (t) = [ω(0)/ω(t)] 1/2 . As a result, α ≈ β ≈ 0 and the adiabatic evolution yields

Single-Particle Time-Dependent Quantum Harmonic Oscillator
Equation (24) determines exactly non-equilibrium energy fluctuation as a function of the initial expectation values at time t = 0. However, it is not possible in general to obtain a simple expression similar to (17) for the mean value of the non-equilibrium energy. In (23b) we use that the mean value of the squeezing operator vanishes, but this is not true for its second moment. This makes the explicit computation model dependent and challenging, in general. To illustrate this observation, let us evaluate n|C 2 |n ≡ Ψ n (0)|C 2 |Ψ n (0) for a single particle in a harmonic oscillator that corresponds to H 0 for N = 1, with Hamiltonian H 0 = p 2 /2m + mω 2 0 x 2 /2. Using the creation and annihilation operators, one finds The nonadiabatic expectation value of the second moment of the Hamiltonian then reads where Q * is given by Equation (18), E n (t) =hω(t)(n + 1 2 ) is the instantaneous energy eigenvalue at time t and H t t = Ψ n (t)|H(t)|Ψ n (t) = Q * E n (t) the corresponding nonadiabatic mean energy.
Using (24), our result reproduces the prediction by Husimi [48] and differs from that by Chen and Muga [49] (as the latter result does not vanish in the adiabatic limit, it is incorrect). We verify that in the adiabatic limit whenḃ = 0, b = b ad and Q * = 1, energy fluctuations in a single driven energy mode vanish identically, as expected.
For a canonical thermal state as initial state, with density matrix ρ 0 = e −βH(0) /Z 0 normalized in terms of the partition function

Calogero-Sutherland Gas in a Time-Dependent Harmonic Trap
We next consider the nonadiabatic energy fluctuations in a driven many-body quantum system. We focus on the (rational) Calogero-Sutherland model describing one-dimensional bosons subject to inverse-square pairwise interactions [50,51]. The system Hamiltonian is In one dimension, interactions are inextricably linked to particle statistics. The notion of generalized exclusion statistics [52,53], complementary to the familiar exchange statistics, classifies particles according to their exclusion in Hilbert space. In this context, the Calogero-Sutherland model can be described as a non-interacting gas of particles obeying generalized exclusion statistics [54]. The constant g ≥ 0 thus plays a dual role. While in the standard description it determines the strength of the inverse-square interactions, it also serves to parameterize generalized exclusion statistics. The value g = 0 corresponds to ideal bosons. The case g = 1 corresponds to hard-core bosons in the continuum, i.e., the so-called Tonks-Girardeau gas [55]. For g = 1/2 the system is sometimes referred to as composed of Haldane semions, and more generally particles are called geons for an arbitrary value g of the exclusion parameter. The many-body wavefunction of the Calogero-Sutherland ground-state at t = 0 reads where C 0 is a known normalization constant, see e.g., [35], and the corresponding energy eigenvalue reads E 0 =hω 0 N[1 + g(N − 1)]/2. Obviously, for g = 0 on recovers the wavefunction an ideal Bose gas, while for g = 1 one recognizes the wavefunction of bosons in the Tonks-Girardeau regime [56]. The time evolution of (31) exhibits scale invariance [7,35] and takes the form of Equation (2) in one spatial dimension with D = 1. As the ground state is an energy eigenstate, the expectation value of the second moment of the Hamiltonian simply reads H 2 0 = E 2 0 . For the computation of Q 2 and C 2 we note that these moments can be derived from the corresponding characteristic function, which can be written down as an overlap between ground-state wavefunctions. In particular, as the ground-state is a homogeneous function times a Gaussian term, the explicit computation of the overlap is possible [35,57].
We first proceed to the computation of Q 2 by introducing the overlap from which moments ofQ are found as Their explicit expressions are given by where for compactness we have defined c = N 2 [1 + g(N − 1)]. Similarly, we can evaluate the term C 2 . To this end, we first notice that the dilation operator acts on a many-body wavefunction as We next evaluate the overlap which can be understood as the characteristic function of the eigenvalue distribution of the squeezing operator. Thus, moments of the generator can be found using the identity As A C (b) is a real function, this implies that odd moments Ψ 0 |Ĉ k |Ψ 0 identically vanish in the ground state |Ψ 0 . By explicit evaluation, we find Thus, the second moment of the Hamiltonian is set by where Q * is given by Equation (18) and E 0 =hω 0 c. As a result, the variance simply reads and scales at most quadratically with the number of particles [35]. The exact many-body wavefunctions of the Calogero-Sutherland model are known for excited states and can be written as the product of the ground-state wavefunction and a combination of Hermite polynomials, see e.g., [58]. Due to the later, the wavefunction of an arbitrary excited state is no longer a homogeneous function times a Gaussian term, precluding the computation of the energy fluctuations via the method put forward.

Unitary Fermi Gas in a Time-Dependent Harmonic Trap
We next consider a strongly interacting degenerate Fermi gas with two spin components [59]. The unitary limit is reached when the interaction strength between atoms becomes divergent, e.g., at a Feshbach resonance. To characterize the spectrum and eigenstates of the Hamiltonian, the interactions can be replaced by a contact condition in terms of a function A such that where R ij = ( r i + r j )/2 is the center of mass position for particles i and j, and the limit of vanishing relative distance between two particles r ij = | r i − r j | → 0 is taken. Tan, Werner and Castin have shown that the low-energy wavefunctions of a unitary Fermi gas in an isotropic trap are the product of a Gaussian factor times a homogeneous function [60,61]. This result together with the scale invariance exhibited by this system [11] will make possible the exact computation of the nonadiabatic energy fluctuations.
In three spatial dimensions, we introduce the 3N vector X = ( r 1 , . . . , r N ). Similarly, we denote by P = ( p 1 , . . . , p N ) the 3N vector of the conjugated momenta. Denoting the norm of X by X = ∑ j r 2 j and the unit vectorn = X/X, we make use of the hyperspherical coordinates (X,n). Eigenstates in a trap can be related to zero-energy eigenstates in free space. We follow Werner and Castin [61] and denote by ψ 0 ν ( X) a zero-energy eigenstate in free space. The latter exhibits scale invariance where ν is the degree of homogeneity, i.e., X · Pψ ν = −ihνψ ν [60]. The corresponding trapped eigenstate reads and its energy eigenvalue is given by in terms of a non-negative integer q. Energy eigenstates are connected by raising and lowering operators associated with the SO(2, 1) Pitaevskii-Rosch symmetry. The lowest energy eigenstate of a given tower of states, that with q = 0 that is annihilated by the lowering operator, can be shown to take the form where C is a normalization constant. This completes our review of [60,61].
By scale invariance, the time-dependent state following a modulation of the trapping frequency is of the form in Equation (2). As it turns out, this suffices to compute the nonadiabatic expectation value of the first and second energy moments of the Hamiltonian of a unitary Fermi gas initially prepared in the low-energy state (49). To this end, we note that the characteristic functions of Q = X 2 and C = ( X · P + P · X)/2 are and thus whence it follows It is convenient to rewrite this expression in the form in parallel with the results of the single-particle harmonic oscillator as well as that of the Calogero-Sutherland model with the corresponding ground states as initial state. This last expression can be rewritten as follows with Q * given by Equation (18), leading to the expression of the energy variance in the form The expression of the results in terms of the ratio is particularly convenient and it also holds for the results derived in the previous subsections for the ground states of the harmonic oscillator and the Calogero-Sutherland model. As a result, we shall use it hereafter without specifying the system under consideration. Beyond these models, our results hold whenever the initial state is the product of a Gaussian term times a homogeneous function. We further note that which is but the initial size of the cloud in units of x 0 = h/(mω 0 ), and is thus an experimentally measurable quantity.

Nonadiabatic Moments of the Square Position Operator
Due to the SU(1, 1) dynamical symmetry group, we have related moments of the Hamiltonian H(t) to those of the operators H(0), Q and C at t = 0. For completeness, we note as well the nonadiabatic evolution of the moments of the square position operator Q. To this end, we introduce Noticing that As expected, the k-th momentum of the position operator scales in a trivial form, and is proportional to the corresponding moment at t = 0, in agreement with the definition of the scaling factor b(t).

Nonadiabatic Moments of the Squeezing Operator
The evolution of the squeezing operator C is more involved than that of Q, and as with the non-adiabatic energy, we shall only compute the first two moments explicitly. We first note that The first moment reads while the second one is given by after using {Q, C} = 0.

Driving Protocols
In what follows, we analyze the exact nonadiabatic energy fluctuations for specific instances of driving protocols of the time-dependent harmonic frequency.

Free Expansion
A standard protocol for time-of-flight imaging of ultracold atoms relies on suddenly releasing a cloud initially trapped in a harmonic potential with frequency ω 0 . For t > 0, ω(t) = 0. For an initial equilibrium state at t > 0, the initial conditions b(0) = 1 andḃ(0) = 0 are required for the time-dependent state (2) to reduce to the initial state. The solution of the Ermakov equation is then b(t) = 1 + ω 2 0 t 2 [9,10]. In this case there is no reference notion of adiabatic dynamics and Q * is not defined. The mean energy satisfies i.e., the mean total energy is preserved during time of flight. The same holds true for the expectation value of all the moments of the Hamiltonian. In turn, for an energy eigenstate H 2 0 = H 0 2 and thus

Linear Frequency Ramp
As an example of a generic nonadiabatic driving we consider a linear ramp of the frequency such that ω(t) = ω 0 for t < 0, ω(t) = ω 0 + (ω f − ω 0 )t/t F for 0 < t < t F and ω(t) = ω F for t ≥ t F . Such time-dependence is nonanalytic at t = {0, t F } and is thus expected to give rise to nonadiabatic dynamics [73]. The limit of t F → 0 naturally reduces to a sudden quench. The scaling factor can be obtained by solving numerically the Ermakov equation subject to b(0) = 1 andḃ(0) = 0 and is shown in Figure 1. The nonadiabatic factor increases monotonically as a function of time with the final mean-energy being four times higher than the corresponding adiabatic value. The energy variance exhibits as well a monotonic behavior, with lower values reflecting the expansion nature of the process, during which the instantaneous trapping frequency decreases. Such features are however dependent on the parameters chosen and a non-monotonic behavior of b(t), Q * and ∆H(t) 2 is observed whenever t F is of the order of the inverse of the frequencies involved in the process, e.g., ω −1 0 and ω −1 F . We note that exact closed-form solutions of the Ermakov equation can be constructed as discussed in [74].

Shortcuts to Adiabaticity by Reverse Engineering
Reverse engineering provides a convenient approach for the fast driving of a system from an initial equilibrium state towards a final nonequilibrium state, without relying on adiabatic strategies. Fixing the duration of the protocol to be t F , with initial and final frequencies ω 0 and ω F , we proceed to impose the following boundary conditions on the scaling factor: b(0) = 1 and b(t F ) = b ad (t F ) = √ ω 0 /ω F . Additional boundary conditions are required for the exact solution (2) to reduce to a stationary state. The vanishing of the derivates of the scaling factoṙ b(0) =ḃ(t F ) = 0 make sure that Q * = 1 at t = {0, t F }. It proves convenient to further imposë b(0) =b(t F ) = 0. With these conditions one can then fix an interpolating polynomial ansatz for b(t) and determine the driving protocol ω(t), using the Ermakov equation as ω(t) 2 = ω 0 /b 4 −b/b. This is the essence to design a shortcut to adiabaticity by reverse engineering the scale-invariant dynamics [14,15,28], as experimentally demonstrated in [18,19,22]. Taking the interpolating ansatz Figure 2 shows the modulation of the frequency-squared ω(t) 2 , nonadiabatic factor Q * and nonadiabatic energy fluctuations for a protocol with t F = 10/ω 0 leading to an expansion with final b(t F ) = 4. The duration time is chosen such that ω(t) 2 ≥ 0 for all 0 ≤ t ≤ t F , preventing the inversion of the trap associated with values ω(t) 2 < 0 required for shorter processes in which t F ∼ ω −1 0 . The nonadiabatic nature of the protocol is manifested in the values of Q * exceeding unity for 0 < t < t F . Nonetheless, Q * = 1 and ∆H(t) 2 = 0 for t = {0, t F } showing that at the beginning and end of the protocol excitations are cancelled out and the nonadiabatic evolution reduces to a stationary energy eigenstate. The scaling factor, frequency square, nonadiabatic factor and energy fluctuations are shown from left to right during a shortcut to adiabaticity for an expansion with final scaling factor b(t F ) = 4 and duration t F = 10/ω 0 . The latter is chosen long enough so that ω 2 (t) > 0 during the process. While values of the nonadiabatic factor Q * > 1 manifest the nonadiabatic nature of the process, the time-evolving state reduces to a stationary equilibrium state at the beginning and end of the protocol, when Q * = 1. The nonadiabatic energy variance, when normalized byhω 0 E(0), accounts for the behavior of the ground state of the harmonic oscillator, the Calogero-Sutherland gas, and the low-energy state of the unitary Fermi gas discussed. Its vanishing value at t = {0, t F } indicates that the initial and final state are energy eigenstates of the corresponding Hamiltonians, H(0) and H(t), respectively.

Local Counterdiabatic Driving
Consider an arbitrary modulation of the trapping frequency ω R (t) as a reference, interpolating between a prescheduled initial value ω R (0) = ω 0 and a final one ω R (t F ) = ω F . Such modulation is generally nonadiabatic. Shortcuts to adiabaticity engineered by local counterdiabatic driving [16,22,23,28,31] suppress nonadiabatic excitations by implementing instead the driving protocol with frequency ω(t), where Driving the system with ω(t), it follows from the Ermakov Equation (3) that the scaling factor is precisely With the boundary conditions,ω R (0) =ω It is then guaranteed that ω(0) = ω R (0) = ω 0 and ω(t F ) = ω R (t F ) = ω F , so the local counterdiabatic protocol ω(t) still interpolates between the same initial and final frequencies considered in the reference protocol ω R (t). Using Equations (69) and (70), the nonadiabatic factor (18) reads in this case This result corrects the expression found in [28,75]. The nonadiabatic energy variance takes the form and vanishes identically at t = {0, t F }. Figure 3 shows a polynomial modulation of the reference frequency-square ω 2 R (t) that decreases monotonically during the process from ω 0 to ω F . By contrast, the modulation of ω(t) 2 associated with the local counterdiabatic driving protocol is highly non-monotonic, involving a fast expansion, followed by a compression and an expansion again. This specific modulation as a function of time strongly depends on the duration of the process and is shown here for t F = 2/ω 0 and an expansion with ω F = ω 0 /2. The corresponding nonadiabatic factor Q * is shown to be greater than unity but at the beginning and end of the protocol when excitations are cancelled. The nonadiabatic energy variance grows from zero value and fluctuates during the process until its completion, when it vanishes again, as the state at t = t F is an eigenstate of the final Hamiltonian.  Figure 3. Local counterdiabatic driving of the scale-invariant dynamics. Given a reference modulation of the trapping frequency ω R (t) between the initial and final trap configurations, a shortcut to adiabaticity by local counterdiabatic driving implements a modulation of the trapping frequency ω(t) that differs from the reference one ω R (t). The reference modulation of the frequency square, the frequency square under the local counterdiabatic driving, the corresponding nonadiabatic factor and energy fluctuations are shown from left to right for an expansion with final frequency ω F = ω 0 /2 and duration t F = 2/ω 0 . The latter is chosen long enough so that ω 2 (t) > 0 during the process.

Conclusions and Discussion
Scale invariance is an important symmetry that governs the dynamics of a variety of systems, including many harmonically-trapped ultracold atomic gases. It is not only of relevance to time-of-flight imaging techniques, but it also plays an important role in the implementation of quantum control techniques and the characterization of finite-time thermodynamics. We have shown that nonadiabatic energy fluctuations in many-body systems under scale-invariant driving can be simply expressed in terms of the equilibrium properties of the initial state. In particular, knowledge of expectation value of the second moments of the initial Hamiltonian, square position and squeezing operators suffices to characterize the exact nonadiabatic dynamics in terms of the evolution of the scaling factor, which governs the size of the atomic cloud. As the later is an experimentally measurable quantity, our results provide a way to characterize nonadiabatic energy fluctuations in the laboratory using time-of-flight imaging techniques. To illustrate our results, we have considered the single-particle quantum oscillator of relevance to experiments with trapped ions and thermal atomic clouds, as well as two instances of many-body systems: the Calogero-Sutherland gas and the unitary Fermi gas. The corresponding nonadiabatic evolution has been characterized in expansion protocols involving sudden quenches and finite-time protocols, including linear frequency ramps, as well as shortcuts to adiabaticity leading to stationary states upon completion of the process.
The exact results for the nonadiabatic energy fluctuations we have provided are expected to have widespread applications in both theoretical and experimental studies of these systems. For example, they can be used to characterize far-from-equilibrium dynamics in ultracold gases, quantify the cost of quantum control such as the use of optimal protocols [68,[76][77][78] and shortcuts to adiabaticity [41,42], characterize the performance of devices and processes in quantum thermodynamics [30], and study the ultimate speed limits [37,38], quantum decay [35] and orthogonality catastrophe [79] under scale-invariant quantum dynamics.