Quantum Critical Scaling under Periodic Driving

Universality is key to the theory of phase transitions, stating that the equilibrium properties of observables near a phase transition can be classified according to few critical exponents. These exponents rule an universal scaling behaviour that witnesses the irrelevance of the model’s microscopic details at criticality. Here we discuss the persistence of such a scaling in a one-dimensional quantum Ising model under sinusoidal modulation in time of its transverse magnetic field. We show that scaling of various quantities (concurrence, entanglement entropy, magnetic and fidelity susceptibility) endures up to a stroboscopic time τbd, proportional to the size of the system. This behaviour is explained by noticing that the low-energy modes, responsible for the scaling properties, are resilient to the absorption of energy. Our results suggest that relevant features of the universality do hold also when the system is brought out-of-equilibrium by a periodic driving.

follows the same scaling 3 behaviour proper of the equilibrium QCP, even though the state itself is very far from the critical ground state. This behaviour persists up to a stroboscopic time scale, τ bd , where scaling breaks down, thus setting a condition for the observation of quantum critical scaling in periodically driven many body systems.

Results
Periodically driven Ising model. We investigate the 1D quantum XY-model, driven by a periodic transverse magnetic field: where σ α (α = x y z , , ) are the Pauli matrices, sin( ) is the harmonically modulated transverse field, and γ the anisotropy parameter. For γ ∈ (0, 1], this model belongs to the Ising universality class and it exhibits a second order QPT with a critical point located at h = h c = 1, separating a ferromagnetic phase from a paramagnetic one. The XY-model with a static field is diagonalised by standard Jordan-Wigner (JW) and Bogolyubov transformations 1, 36 , enabling Eq. (1) (with Δh = 0) to be re-written as a free fermion Hamiltonian  Floquet analysis is a valuable tool to deal with time-periodic Hamiltonians, = +Ĥ t H t nT ( ) ( ) , as it allows to reduce the stroboscopic time evolution, i.e. at integer steps n of the period T, to a dynamics generated by a time-independent effective Hamiltonian 9 . Indeed, by exploiting the periodicity of Ĥ t ( ), the (stroboscopic) unitary time evolution operator Û t ( ), at times t = nT, can be written as a discrete-time quantum map . Accordingly, after n periods of the sinusoidal drive, the initial critical state (3) evolves into Ψ nT ( ) will be the stroboscopic state where we test the persistence of finite size scaling (FSS) behaviour, the characteristic trait of criticality, in several quantities and for a number of different driving conditions. Critical scaling under periodic drive. In the following, we provide evidence that the scaling behaviour of a number of physical quantities, which should strictly hold at the equilibrium QCP only, persists in fact also under periodic driving. In particular, we consider both local-in the real lattice space, quantities (e.g., nearest-neighbor concurrence and local transverse magnetic susceptibility) and non-local ones (e.g., entanglement entropy and fidelity susceptibility).
Remarkably, scaling behaviour at equilibrium has been found also for quantities that are not observables in the strict quantum mechanical sense, such as entanglement. Indeed, the concurrence 37 , a measure of bi-partite entanglement for two qubits (see Methods), has been shown to exhibit FSS at equilibrium QCP, as first demonstrated for the quantum Ising model 38 and later illustrated in other systems [39][40][41][42][43] . The scaling of concurrence was explored in details, as it bridges QPT with quantum information theory 44,45 and quantum thermodynamics, because of its close connection with the notions of irreversible work 46 and ergotropy 47 .
In the quantum Ising (QI) model at the equilibrium critical point 38 , the derivative of the concurrence between neighboring spins with respect to the transverse field, h, displays a logarithmic singularity at h = h c . An FSS analysis at equilibrium for the concurrence shows data collapse for different system sizes, consistent with the universal critical exponent ν = 1, which also governs the divergence of the correlation length of the order parameter in the QI model.
Interestingly, we find that the same scaling property still holds for a periodically driven, out of equilibrium system. The nearest-neighbor concurrence C i,i+1 (N) for different system sizes N is reported in the upper left inset of Fig. 1 both at equilibrium, nT = 0, and after time nT = 30 of the driving π = + . h t h t ( ) 0 1sin (2 ). Notice that, notwithstanding the system is driven out-of-equilibrium, the qualitative features of C i,i+1 (N) around the critical point h c are preserved. Indeed, in the main plot of Fig. 1, we report the derivative of the stroboscopic nearest-neighbor concurrence, , in the neighborhood of h c for different system sizes N and after = = n t T / 30 cycles of the driving . It is interesting to note that the logarithmic correction to the shift exponent λ = 2 is shared also by the equilibrium FSS behaviour of the half-chain entanglement entropy 48  (see the lower right inset of Fig. 1). Performing the FSS analysis for logarithmic divergences 3 , we obtain data collapse once the scaling exponent is set to its equilibrium value ν = 1, as reported in the upper right inset of Fig. 1. FFS is attained with the same value of the critical exponent ν even after changing driving amplitude, frequency, anisotropy parameter as well as number of cycles, provided the conditions outlined in the following subsection are fulfilled. This highlights that the role played by universality stretches well beyond the ground state properties, significantly affecting the system even under periodic driving. In the following we will substantiate further our claim about the persistence of FSS in periodically driven critical systems by considering additional quantities. For this purpose, we have also considered the scaling relation for the transverse magnetic susceptibility of the XY-model, . At equilibrium, it exhibits a scaling behaviour with critical exponent α = 0 49 , implying a logarithmic divergence akin to the one encountered for the concurrence. Such logarithmic divergence is preserved also under driving, and, in analogy with the scaling Ansatz for the concurrence, data collapse is obtained for different system's sizes, implying that even after the system is brought out-of-equilibrium by periodic driving, the scaling exponents keep their equilibrium values α = 0 and ν = 1. So far, we have considered single and two-site quantities. Still, physical quantities having support on a larger part of the system, e.g., the entanglement entropy, or even genuinely global, such as the Fidelity susceptibility, exhibit FSS at equilibrium in the critical QI model 48,50 . The former following a logarithmic and the latter an vs. initial magnetic field h around the QPT point for different system's sizes N. A logarithmic divergent behaviour is found close to h = 1, as inferred from the right lower inset, where the maximum of the peaks of vs N are reported in a semilogarithmic plot. Upper Left inset: Nearest-neighbor Concurrence C i,i+1 (N) vs. initial magnetic field h at equilibrium, for an infinite chain (red-dotted line) and after n = 30 cycles for N = 128, 256, 512, 1024, 2048, 4096, (green continuous lines: different lines are not distinguishable on this scale). Upper right inset: Data collapse for FSS at the logarithmic divergence, for the same N as listed above, after n = 30 cycles, according to The driving protocol is h(t) = 1 + 0.1 sin 2πt, and we took γ = 1 (cfr. Eq. 1). The chosen value of n = 30 is within the breakdown time of the shortest chain here considered (see following section). For n < 30, although the numerical values of (and C i,i+1 (N)) change, the FSS data collapse and the logarithmic divergence, upper right and left inset, respectively, is attainable with the same critical exponent. algebraic divergence, respectively. As for the Entanglement Entropy of the half chain S N/2 , defined by the von Neumann entropy based on the logarithmic law of the entanglement entropy at criticality, derived in ref. 48, holds as well under periodic driving up to times where v max is the maximum group velocity of the Floquet quasiparticles 51-53 (see Methods). Indeed t N/2 gives a time after which the quasi-particles have left the half chain and a volume law for the Entanglement Entropy is attained 51,54,55 .
Let us now turn our attention to the Fidelity susceptibility (FS). The ground state FS is defined by ) and it depends on three length scales; namely, the system size N, the correlation length ξ ∼ −  2 2 , whose scaling behaviour will be dictated by the other two length scales (for a detailed analysis of the use FS in transverse field spin models see the book by Dutta et al. in ref. 50). For the Ising model at equilibrium, it has been shown that, at criticality, where ξ  N , the FS exhibits a maximum whose height scales algebraically with the system size as . Far from criticality, on the other hand, where ξ  N , the scaling is extensive 58,59 . In the following we will investigate the former limit as the driving is around criticality. For the analysis of the stroboscopic Fidelity susceptibility in the limit where ξ  N , see the Section Methods. Notice that, contrary to the quantities considered before, χ F N does not scale logarithmically with N. Notwithstanding, FSS is still attainable, provided a new, time-dependent exponent is introduced, which takes into account the fact that the algebraic scaling gets modified. At equilibrium, one considers the susceptibility of the ground state fidelity, for δh → 0, which, by definition, is time-independent and can be related also to the irreversible work in an infinitesimal quench protocol 60 . On the other hand, in the presence of the driving, we will consider the susceptibility of the following expression for the Fidelity: , where Û and ′ U differ as they correspond to driving around h and h + δh, respectively.
reduces to the ground state Fidelity at t = 0. In Fig. 2, we report our results for the fidelity susceptibility, χ F N , obtained by taking δh = 10 −5 and the same driving parameters as in Fig. 1. Interestingly enough, we find that the algebraic divergence is preserved, according 1 32 ) is a monotonically increasing function (see left inset in Fig. 2). This behaviour can be qualitatively explained by noticing that the driving ultimately will invalidate the FSS behaviour (hence a linear scaling with N is retrieved) at later times (see following subsection). As a consequence, we modify the FSS Ansatz 61 , by introducing a time-dependent exponent r(nT) (r(0) = 1), where sgn(·) is the sign function. In the main plot of Fig. 2 we show an instance of the data collapse obtained by means of Eq. 4 after n = 15 periods, where r(nT) = 0.86. By choosing a different number of periods, the exponent r(nT) decreases following the curve reported in the right inset of Fig. 2. It is worth mentioning that our r(0) = 1 coincides with the equilibrium value, which is given by and d a = 1 are the critical adiabatic dimension and the adiabatic dimension, respectively 61 . The d's cannot be easily defined out-of-equilibrium, and, therefore, we relied on numerical tools to extract the best collapse exponent r(nT). Nevertheless, notice that the critical exponent ν equals 1 at all times where the FSS holds.
Notably, all the quantities whose FSS we have reported so far for the Ising model (γ = 1 in Eq. 1), maintain this scaling behaviour, with the same critical exponents, in the whole XY universality class 0 < γ ≤ 1.
Breakdown time of finite-size scaling. The FSS under periodic driving holds up to a characteristic time τ bd (N), which we dub breakdown time, and which depends on the driving parameters. Remarkably, we obtain for ω > 4 and ∆  h h c , a breakdown time comparable to the recurrence time τ rec m ax ), while for smaller frequencies, the breakdown time occurs before the onset of recurrences (τ bd < t rec ). Hence, for a large frequency and small amplitude driving, the breakdown of critical FSS could be arbitrarily delayed by taking systems of a large-enough-size. For the high-frequency limit, see also a recent paper by Gritsev and Polkovnikov 62 , where it has been shown that the Floquet Hamiltonian of a step-like periodically driven Ising model shares the same critical properties of the time-independent Ising model based on the structure of the Onsager algebra and the self-duality (in the limit of ω → ∞, the system freezes in its initial state and FSS is attainable for every n and δh, but this would fall us back into the equilibrium scenario).
On the other hand, the breakdown of FSS at times smaller than t rec that occurs for small frequencies (ω < 4) is mostly due to the fact that systems with small sizes loose their scaling properties. As an instance of such a behaviour, in Fig. 3 we display τ bd (N), as extracted from the transverse magnetisation susceptibility. In the upper left panel, we show that the breakdown time is comparable to the recurrence time for small driving amplitudes; conversely, by increasing Δh, the number of cycles for which FSS holds reduces significantly, as can be expected by observing that already in the first few cycles a considerable amount of energy is absorbed by the system. In the same figure, we report three instances of scaling analysis performed at increasing number of cycles. FSS holds for all cases but the one corresponding to the system of smaller size, so that all of the curves collapse near h c , but one. , which gives the probability amplitude to find the system in a state close to the initial critical ground state 63 , we find-already after a few cycles-that nT ( )  has become negligibly small. This is explicitly shown in Fig. 4(b), where an exponential decay of  is reported. A closer look, however, shows that the decay of the Loschmidt echo is essentially due to those mode k for which a quasi-degeneracy occurs in the Floquet energies 28 . These almost degenerate modes are also those responsible for energy absorption from the driving, see Fig. 4.
Actually, the Loschmidt echo can be decomposed as a product of contributions from the different modes,   = ∏ k k (see the methods section). Such a momentum resolved Loschmidt echo is shown in Fig. 4(c), where it is superimposed to the dispersion relation of the Floquet quasi energies (given by the black dashed lines on the background). One can see that the minimum of  k is found where an 'inter-band quasi-degeneracy' occurs. Moreover, this 'quasi-resonance' precisely corresponds to the peak of the energy absorbed from the driving, as shown in Fig. 4(a). There, the amount of work performed by the driving agent is reported for each mode k (see Methods). Therefore, we can conclude that the decay of the Loschmidt echo is due to these specific modes, where a Floquet resonance occurs, and which are brought significantly out of equilibrium by the absorption of energy from the driving.
On the other hand, one expects FSS to be essentially a feature of the long wavelength (and low energy) modes. These modes are much less affected by the driving as they absorb much less energy than those close to the quasi degeneracy. Correspondingly, their contribution to the overall decay of the Loschmidt echo is very small. The time dependence of the k-resolved Loschmidt echo for such small-k-modes is reported in Fig. 4(d), where we see that  nT ( ) k periodically oscillates in time. For each k, the period of such oscillations is determined by the k-eigenvalue of Ĥ F (i.e. µ / k  see Methods), which becomes larger and larger with increasing N.
It turns out that the breakdown of FSS occurs close to the time where a minimum is found for the  k corresponding to the smallest k. Therefore, bd k . After this time, in fact, it not possible anymore to obtain the collapse of curves corresponding to physical observables evaluated for different N's. This is due to the fact that systems with different sizes behave asynchronously. This means that, although the  k corresponding to the smallest k goes back to unity after a period, this revival occurs at a time that explicitly depends on the system size. As a result, since the scaling is a comparison of physical quantities for different N at the same time, it does not hold anymore because systems with different sizes have k-resolved Loschmidt echoes that are 'out of phase' from each other.
Finally, let us comment on the fact that for low-ω drivings, the FSS behaviour is lost already after a single cycle because the low-k modes absorb more energy from the drive as the Floquet resonances move towards them, (see Methods).

Discussion
In summary, we have shown that the scaling behaviour proper of an equilibrium quantum critical point retains its validity also when the system is brought out-of-equilibrium via a periodic drive. Finite-size scaling of both local and global quantities, exhibiting logarithmic as well as algebraic scaling with the system size, has been performed. We have shown that the equilibrium critical exponents are robust against the periodic perturbation up to times when the stroboscopic state is far from the quantum critical state, suggesting that the features of universality manifest themselves also under strong periodically driven settings. In addition, our claims are within reach of experimental verification, as out-of-equilibrium quantum Ising dynamics is currently under active investigation via a great variety of different physical systems, ranging from degenerate Bose gases in a optical lattice 64,65 to Rydberg atoms 66 . Our results may find applications also in the emerging field of out-of-equilibrium quantum thermodynamics, where, recently 67 , quantum Otto engines, having as working substance a many-body system at the verge of criticality, have been suggested to be able to attain the Carnot efficiency at finite power because of the validity of the FSS relations.
It would be interesting to study in the future whether the same scenario holds for non-integrable systems hosting a quantum phase transition, as, for instance, the one dimensional Bose-Hubbard model, where, however, the system is eventually driven into an infinite temperature state, and therefore the persistence of critical scaling is expected only in a temporal window delimited by the thermalization time of the system. A different scenario, on the other hand, could emerge for interacting integrable models, such as the antiferromagnetic XXZ Hamiltonian, where the Néel and the XY phase are separated by a second-order QPT and thermalisation is prevented by the integrability of the model.

Methods
Concurrence. The concurrence, C(ρ), is a measure a bipartite qubit entanglement, which can be straightforwardly computed for any two-spin− 1 2 density matrix ρ 37 . In particular, we apply it to evaluate the entanglement between two spins of the chain, residing at sites i and j. In this case, the two spin reduced density matrix, ρˆi j , , is obtained by computing the partial trace over all but the i-th and j-th spin degrees of freedom of either (i) the ground state |GS〉, or (ii) the stroboscopical state Ψ nT ( ) . Given the state, C i,j can be evaluated, once the state ρˆi j , is expressed in the logical basis of the eigenstates of the σ z operator, via the relation Loschmidt echo. The Loschmidt echo is defined as Using the explicit form of the initial and stroboscopic state, namely where the k-resolved Loschmidt echo is given by For an extensive analysis of the Loschmidt echo in periodically driven systems see ref. 63. Work. The (average) work performed up to time t by driving the system is given by the difference between the average instantaneous energy of the system and its initial value given by the ground state energy. At the discrete time instants t = nT, we have where we used the fact that =Ĥ nT H ( ) (0). Using the explicit expression of Ĥ in terms of the Bogoliubov fermions, Eq. 2, we have that the constant terms cancel out and that the work naturally decomposes into the sum of contributions arising from each mode k, where Scientific RepoRts | 7: 5672 | DOI:10.1038/s41598-017-06025-1 Notice that, since the Hamiltonian undergoes a periodic driving, and we are evaluating the work at an integer number of periods, the average work coincides in our case with both the so called irreversible work and the inner friction 25 , so that it can be used to describe also the amount of irreversibility brought into the system.
are, respectively, the positive and negative Floquet eigenvalues and eigenvectors of the Floquet Hamiltonian ĥ kF for the mode k. The evolution operator = −Û T e ( ) k iTh kF is determined by the solution of the Bogoliubov-de Gennes equations as those reported in Fig. 4 in the main text. As a consequence, resonances can occur both within the same band and between different bands, dubbed intra-band and inter-band resonances, respectively, in the main text.
Finally, as the stroboscopic evolution of the initial state is given by , the dynamics is governed by the Floquet quasiparticles energies μ k and hence the maximum velocity quasiparticles can spread out is given by = µ v max max k d dk k , that is the maximum of the "group velocity" as given by the "dispersion relation" Extensive scaling of the Fidelity susceptibility far from criticality. In the main text we have investigated the fidelity susceptibility (FS) in the limit where the order parameter correlation length ξ ∼ − ν − h h c dominates over the system size N, i.e., ξ  N . In such a regime, the FS scales as N 2 at equilibrium and FS finite-size scaling has been reported in Fig. 2  , the FS scales at equilibrium linearly with N. Here we will show that such a linear scaling is preserved also under periodic drive. In Fig. 5 we report the fidelity susceptibility both at equilibrium and after n = 25 cycles of the driving.
Furthermore, by fixing N, in the limit ξ ∼ −  N h (ln ) 1 , i.e., away from criticality, it is known that at equilibrium the fidelity susceptibility χ h ( ) z N ( ) scales as ξ 50 . In Fig. 6    showing that, for the modes near the quasi-degeneracy, now in the FSS-relevant region, the decay of the Loschmidt echo is both faster and more significant than for higher ω (see Fig. 4 in the main text). In both plots we considered a driving frequency ω = . Low-ω drivings. In the main text of the manuscript we have investigated mainly drivings at frequencies ω > 4, where the FSS behaviour is resilient under the periodic modulation of the magnetic field h(t). In this subsection we show that, in the low-frequency limit, the FSS is lost already after the first cycle. To support this claim, we report in Fig. 7 the k-resolved work and Loschmidt echo for ω = .
0 5. As already stated in ref. 28, by decreasing the frequency of the drive, the number of resonances in the Floquet spectrum increases and, more importantly, they also move towards the low k-modes. As a consequence, the energy injected into the system by the low-frequency driving is absorbed mainly by the latter and also the Loschmidt echo of the low-k modes decreases to values significantly lower than in the ω > 4 case. As a result, the FSS behaviour is not resilient to such low-ω drives.