Activating critical exponent spectra with a slow drive

We uncover an aspect of the Kibble–Zurek phenomenology, according to which the spectrum of critical exponents of a classical or quantum phase transition is revealed by driving the system slowly in directions parallel to the phase boundary. This result is obtained in a renormalization group formulation of the Kibble–Zurek scenario, and based on a connection between the breaking of adiabaticity and the exiting of the critical domain via new relevant directions induced by the slow drive. The mechanism does not require fine tuning, in the sense that the scaling of subleading operators occurs as the dominant power law in an extensive regime of drive parameters. Therefore, it should be observable in quantum simulators or dynamically tunable condensed matter platforms.

We uncover an aspect of the Kibble-Zurek phenomenology, according to which the spectrum of critical exponents of a classical or quantum phase transition is revealed, by driving the system slowly in directions parallel to the phase boundary. This result is obtained in a renormalization group formulation of the Kibble-Zurek scenario, and based on a connection between the breaking of adiabaticity and the exiting of the critical domain via new relevant directions induced by the slow drive. The mechanism does not require fine tuning, in the sense that scaling originating from irrelevant operators is observable in an extensive regime of drive parameters. Therefore, it should be observable in quantum simulators or dynamically tunable condensed-matter platforms.

I. INTRODUCTION
Universality near classical and quantum second-order phase transitions builds on the fact that long-wavelength fluctuations are freed by the fine tuning to the critical point, and dominate the macroscopic physics. This is quantified by only few independent observable critical exponents, universal numbers which are independent of microscopic details and determined exclusively by the system's symmetries and dimensionality. However, universality extends beyond these so-called relevant exponents to an entire spectrum of universal exponents [1][2][3]. Yet, this fully universal information is not easily accessible in static experiments (see, e.g., [2,[4][5][6]), or even numerics, the subleading character of the associated power laws is easily overwritten by the more dominant exponents (see, e.g., [7][8][9][10][11][12][13][14]). Exceptions to this scenario are available in conformal field theories, where a relation between the scaling dimensions of operators and the energy spectrum has been established [15][16][17][18] as well as in transitions with dangerously irrelevant parameters, where specific exponents are accessible [19].
In this work, we demonstrate that irrelevant universal exponents can be turned into relevant ones by slowly driving the system in the vicinity of a second-order phase transition, classical or quantum, thus circumventing the need of fine tuning. This allows for the detection of these exponents in a robust way, and we provide a simple longitudinal drive protocol [20] to do so in dynamical experiments (see Fig. 1a).
We describe the underlying physics qualitatively in Sec. II and derive our main result [detailed below Eq. (9)] in Sec. III. In Sec. IV we choose a rather general model to illustrate our mechanism. In Sec. V we discuss how the breaking of adiabaticity is reflected in the Renormalization Group (RG) flow and how it can be seen to produce two different scales. We show how our mechanism can be observed in practice in Sec. VI, and discuss the more * smathey@thp.uni-koeln.de A generic phase diagram with two independent parameters (locally τ0 and λ0) and a continuous phase transition. In the KZ scenario, the critical line is crossed transversally (blue line, left). As the phase transition is approached, the correlation length increases as ξ ∼ |τ0| −ν until adiabaticity breaks. At this point, the usual KZ scaling with the drive amplitude v, emerges ξ ∼ v −ν/(1+zν) (left inset, log-log). z is the dynamical critical exponent. For more general drives involving a longitudinal component (red line, right), the correlation length may exhibit a different scaling ξ ∼ v −1/(z− ) , which involves an irrelevant exponent . This behavior is realized when the drive amplitude is bigger than a crossover amplitude v * (right inset), that can be made to interpolate between ∞ and 0 by changing the angle φ, at which the phase boundary is crossed: thus, an irrelevant exponent is robustly observable in drive protocols with a strong enough longitudinal component. (b): Spectrum of equilibrium (blue, left) and nonequilibrium (right, red) critical exponents resulting from a linear drive. The drive shifts the whole spectrum down by z and makes new relevant exponents out of irrelevant ones. This mechanism underlies the scaling of the correlation length with the drive amplitude described in (a) (see text).
general case of a polynomial drive in Sec. VII. Finally, we illustrate the full mechanism on an exactly solvable model in Sec. VIII.

II. BASIC PHYSICAL MECHANISM
Our result is obtained by a systematic reformulation and generalization of the Kibble-Zurek (KZ) scenario of diabatic decoupling in a time-dependent adiabatic RG language [3]. The KZ mechanism [21][22][23][24][25][26][27][28][29][30][31][32] describes the behavior of a system when slowly ramping the parameters through a second-order phase transition: At a point in parameter space close enough to the phase transition, where the driving rate becomes comparable to the system's gap, the system crosses over from adiabatic to diabatic. Beyond that point, defects lead to a cutoff for the scaling of the correlation length (and far from equilibrium dynamics [33]). A quantitative prediction of this mechanism is the scaling of the correlation length at the crossover itself; according to the KZ hypothesis, this scaling involves exclusively the exponents of the underlying equilibrium critical point.
Crucially, the quantitative predictions of the KZ mechanism rely on identifying the point where adiabaticity breaks. The key idea of our RG approach is to first formulate adiabatic flow equations, where time enters only as a parameter, and then translate the breakdown of adiabaticity into an RG language to connect to KZ. As a first result of this approach, we demonstrate that the structure of the adiabatic flow equations reduces the KZ hypothesis to the conventional scaling hypothesis: no independent information is encoded in them [27,34,35]. This is because a slow drive only affects the large-scale physics (conversely, a fast drive acts at small scales and can produce new critical exponents [36][37][38][39]). The most important new implication of our analysis is, however, based on the downshift of the spectrum of critical exponents (see Fig. 1b): Exponents that are irrelevant (positive sign) in the static problem now become relevant (negative sign). This implies that they must manifest themselves in largescale physical observables. And indeed, we show that experiments probing e.g. the correlation length can reveal irrelevant scaling exponents, provided that the slow ramp in parameter space crosses the phase boundary at a shallow enough angle (see Fig. 1a).
Our result is best understood by imagining the extreme case where the system is driven along the phase boundary. Then, the correlation length is always finite, but can get large if the system is held close enough to the phase boundary. In that case, the drive needs to be extremely slow for the dynamics to be adiabatic. With such a drive, we can imagine an alternative KZ scenario where the system is driven multiple times along the phase boundary with the drive amplitude held fixed and the distance to the phase boundary decreasing. For large distance to the phase boundary the relaxation time is short and the system is adiabatic, and the latter increases as the former decreases. There is therefore a point where adiabaticity breaks since the drive amplitude is held fixed. The smaller the drive amplitude, the closer this point is to the phase boundary. As in the usual KZ scenario, the correlation length scales with the drive amplitude at this point. The crucial difference, however, is that this longitudinal drive does not involve changing the distance to the phase boundary which is held fixed. It is instead an irrelevant parameter with a different scaling exponent that is driven. Although this parameter does not produce a diverging scale at equilibrium, it can still break adiabaticity when it is driven fast enough. It turns out that driving an irrelevant parameter in this way can actually produce a diverging length scale with its own scaling exponent. The KZ scaling is thus modified as we show in detail below.

III. ADIABATIC RG FLOW EQUATIONS AND EXPONENT SHIFT
In an adiabatic system, equilibration time scales are much shorter than those induced externally. Time then only enters through the parameters characterizing the system's partition function. The full set of nonequilibrium adiabatic flow equations -encoding the exact scaling dimensions of all operators -is then obtained from its equilibrium counterpart by promoting the static parameters to time-dependent ones. They are This equation depends on time through the parameters, which are bundled in the vectorĝ =ĝ(t ) = (ĝ 1 (t ),ĝ 2 (t ), . . . ). We work in rescaled units where the canonical and anomalous components of the scaling dimensions are accounted for by the diagonal matrices D i and the anomalous dimension η. k is a momentum RG scale, and the prime denotes a time derivative. Here, time must be rescaled like any other parameter. This produces the second term on the right-hand-side becausê g depends on the rescaled timet = k z t, where z is the dynamical critical exponent.
In principle, the above equation describes the renormalization of the full drive protocol through the entire time dependence ofĝ. We will, however, focus on linear drive protocols and simplify the problem accordingly. But first, we briefly sketch the derivation of the rescaled RG flow equation (1), starting from dimensionful adiabatic flow equations. In the adiabatic approximation, the dimensionful RG flow equations are given by k∂ k g(t) = B(g(t), Z(t), Q(t)) , k∂ k Z(t) = η(g(t), Z(t), Q(t)) Z(t) , k∂ k Q(t) = z(g(t), Z(t), Q(t)) Q(t) .
(2) Z(t) and Q(t) are the field strength and time rescaling factors respectively. At an RG fixed point they scale as Z ∼ k η and Q ∼ k z and provide the anomalous dimension and the dynamical scaling exponent respectively. These equations are obtained from the dimensionful flow equations of the equilibrium system (which are computed in a standard way [3,58]) by promoting all the parameters to time-dependent ones. This is the definition of an adiabatic time dependence. Eq. (2) is converted to its dimensionless form [Eq. (1)] by rescaling the time and the coordinates of g with appropriate powers of the RG scale k, and the rescaling factors Z and Q, We define where the rescaling (choice of D 1 and D 2 ) is made in such a way that k, Z and Q drop out on the left-hand-side. Then, we obtain Eq. (1) with the anomalous dimension and the dynamical critical exponents being functions of the rescaled parameters Again, this is a standard procedure in RG theory [3,58]. Notably however, due to time being involved as an additional parameter to be rescaled in the presence of the slow drive, an additional term in Eq. (1) (the middle term on the right-hand side) arises because the rescaled couplings now also depend on the cutoff through Q [see Eq. (3)]. Now, to fix ideas, we choose the adiabatic time dependence of the parameters to beĝ(t ) =ĝ 0 +ĝ 1t . Expanding Eq. (1) to first order int then leads to This equation describes a stationary system (with parameters given byĝ 0 ) in the presence of a small drive amplitudeĝ 1 . The explicit time dependence of Eq. (1) is traded for a doubling of the number of flowing parameters. The entire system's beta functions are obtained by joining Eqs. (6) into a single vector, k∂ k g = (k∂ kĝ0 , k∂ kĝ1 ) = β( g), with g = (ĝ 0 ,ĝ 1 ). The critical physics is then characterized by the equilibrium fixed point g fp = (g fp , 0), that satisfies β( g fp ) = 0 and the flow close to it. In particular, if we define the components of G(k) are small in the critical region, and the flow is approximated as with M is the Jacobian matrix of β( g), evaluated at the fixed point (see App. A). We see from Eq. (8) that the eigensystem of M plays an important role at criticality. Its eigenvalues are the critical exponents and its eigenvectors are used to define the corresponding parameters (see App. B). As a result of the adiabatic structure of the problem, the stability matrix M has a relatively simple structure (i) Its elements are square matrices of the same dimensionality asĝ. (ii) Its lower-left block vanishes identically because nonzero entries would signal that drive parameters (g 1 , left-hand side) can be generated from equilibrium ones (g 0 , right-hand side) alone. (iii) Its diagonal blocks M 0 and M 1 , represent the equilibrium and nonequilibrium critical physics, respectively. They are given in App. A. (iv) As a result of the adiabatic setup, the diagonal blocks are directly related: The upper-right block of M describes the part of the renormalization of the equilibrium parameters that is generated by the drive. Although we find X = 0 from Eq. (6), taking into account nonadiabatic correction leads to X = 0 (see App. A). The critical exponents (including all subleading corrections) can be identified with the eigenvalues of M [see Eq. (8)]. The upper-triangular-block structure of Eq. (9), implies that these are the eigenvalues of M 0 and M 1 , and are thus independent of X. This leads to two key observations: First, the upper-left block is the only element of M capable of producing independent exponents. This provides an RG justification of the KZ hypothesis: a slow drive does not produce any new critical exponents (see also [34,35]). Second, the relation M 1 = M 0 − z implies that the spectrum of critical exponents is doubled with the nonequilibrium copy shifted downwards by a factor z (see Fig. 1b). Intriguingly, a parameter that is irrelevant at equilibrium, and characterized by a positive exponent , can thus be made relevant if < z. This means that this exponent must be associated with a diverging length scale, i.e. it must be observable macroscopically. This will be elaborated on below.

IV. EXAMPLE
We now consider a concrete model. We show how the newly relevant operators can be observed in Sec. VI. The general picture depends, however, on the precise form of the exponent spectrum and must be evaluated separately for each universality class. In particular, different combinations z and drive protocol will lead to different scaling behaviors. We choose interacting O(N ) models. These provide common examples of well-known universality classes such as the Ising (N = 1) [40] and Bose condensation (N = 2) transitions. The spectrum of critical exponents of O(N ) models is estimated in [59] for all values of N and d = 3. For example, the first four critical exponents of the Ising phase transition are found to be . In all cases, θ 1 = −y is the only relevant exponent and θ m>1 > 0 is irrelevant. The dynamical critical exponent z, can take different values that depend on the dynamical conservation laws [3]. We choose the example of a dissipative order parameter (model A) where z ∼ = 2 + 0.36(N + 2)/(N + 8) 2 for d = 3 [3,58]. Then, a linear drive changes the relevance of a single parameter, and the four most relevant exponents, are associated with three relevant and one irrelevant parameters. The corresponding equilibrium system is controlled by the two parameters associated with the relevant and irrelevant exponents −y and , which we refer to as τ 0 and λ 0 , respectively. We denote the corresponding drive parameters as τ 1 and λ 1 , which are associated to −z − y and − z, respectively. See App. B and Eq. (B4) for a precise relation between these parameters and g.

V. ADIABATICITY BREAKING AND RG ANALYSIS
The RG approach provides a direct connection between the physical breaking of adiabaticity and the emergence of a scale resulting from the drive: adiabaticity can be quantified by taking the ratio of the relaxation time 1 T r = 1/g 1 0 , and the drive time-scale T α = |g α 0 /g α 1 | for each parameter. Under RG transformations, these quantities acquire a scale dependence Adiabaticity is then signaled by a small such dimensionless ratio, and it gets broken once | (k)| ∼ = 1.
To derive Eq. (11), we focus on linear drive protocols. 2 We write the dimensionful parameters as g = g 0 + g 1 t [see Eq. (3)]. In terms of the microscopic parameters α is given by α = g α 1 /(g 1 0 g α 0 ) . The two time scales 1 We choose the system's gap to be the first coordinate of g 0 . 2 Higher-order drives can, however, be taken into account by noting that they produce additional means of adiabaticity breaking. then acquires an additional index denoting the order of the drive α → α i ∼ g α i (k 0 ) (with i ≥ 1), and the following reasoning can be directly applied.
involved (T r and T α ) are, however, strongly renormalized in the critical regime. Fluctuations are included by using cutoff-dependent parameters and choosing k as small as possible: We have inserted the rescaled variables in the second equality, andĜ α 0 is the α th component ofĜ 0 . We consider the RG flow close to an equilibrium fixed point where g fp = (g fp , 0). For this reason, the coordinates of the distance from the fixed point can be identified witĥ g 1 in the nonequilibrium sector. Close to the fixed point we can expand the above equation to leading order in − → G and obtain Eq. (11).
We emphasize that does not depend on the equilibrium couplings τ 0 and λ 0 , at criticality. Adiabaticity can only be broken if the system is driven. This is a consequence of the block-diagonal structure of M (see App. B 1). Indeed, expanding the solution of the linearized RG flow onto the eigenvectors of M eventually provides (see App. B 2) with τ 1 and λ 1 characterizing the two components of the drive and v α the two eigenvectors of the equilibrium stability matrix [see Eq. (B10)]. τ 0 and λ 0 do not enter in the flow ofĝ 1 , which vanishes identically if the system is not driven. From the point of view of the RG, |ĝ 1 (k)| and therefore | (k)| increases as k is lowered [see Eq. (13)]. To see how this is related to the emergence of a scale, we first recapitulate how the correlation length is extracted from the RG under static equilibrium conditions. The RG flow is initialized at a large momentum scale k = Λ with the physical microscopic parameters, and evolves to effective macroscopic parameters as the RG scale k is lowered. In particular, for a near critical system, the distance from the fixed point is very small and |τ 0 | 1. As a relevant parameter, 3τ 0 (k) grows large under RG transformations, and a scale k 0 , emerges when |τ 0 (k 0 )| ≈ 1: k 0 ∼ ξ −1 provides the correlation length.
As demonstrated above, however, the system of adiabatic flow equations contains additional relevant directions emanating from the fixed point and associated to the slow drive. The three relevant directions for the case of O(N ) models are illustrated in Fig. 2. Along any of these directions, the flow will leave the critical scaling regime once the associated dimensionless parameters grow to O(1). We discuss the physics of the flow initialized close to the fixed point along the primitive axes of the coordinate system first: • The flow along theτ 0 axis corresponds to the case of an undriven system, τ 1 = λ 1 = 0. In that case, the entire flow is controlled by the equilibrium critical exponent and we have | g(k)− g fp | ∼τ 0 (k) ∼ τ 0 k −y , with τ 0 = (T − T c )/T c the reduced temperature. The correlation length -extracted from the solution of the flow equation where |τ 0 (k 0 )| ≈ 1scales as ξ ∼ τ −ν 0 with ν = 1/y. In particular we recover the Ising exponent ν ∼ = 0.65, when y ∼ = 1.54. This reproduces the usual link between the critical exponents Eq. (10), and the scaling of the correlation length.
Thus, all together there are three possible ways for a scale to emerge and two ways to break adiabaticity. We represent the crossovers from the different scaling regimes as gray surfaces in Fig. 2 with the scaling in each region determined by its encompassed axis. The system is adiabatic in the lower-left corner only, where we have ξ ∼ |τ 0 | −ν . Conversely, if the drive is strong (or if |τ 0 | is small) enough, the system is diabatic. Then, the correlation length at the crossover from adiabatic to diabatic behaves according to one of the two nonequilibrium relevant scalings. Crucially, there are two possible exponents: ν 1 = ν/(1 + zν) is the usual KZ exponent, and ν 2 = 1/(z − ) is a new one containing the irrelevant exponent . We emphasize that the difference in exponents is significant: for example, for model A at N −1 = 0, ν 1 = 1/3 and ν 2 = 1.

VI. OBSERVABILITY AND ROBUSTNESS
We now connect these RG findings to concrete observables. To this end, we start from the drive protocol illustrated in the phase diagram of Fig. 1a (in red, on the The crossover amplitude v * , depends on the drive direction (parametrized by the angle φ), and can be made arbitrarily small by choosing φ close enough to π/2. right) and parametrizê v denotes the drive amplitude and φ its direction. τ 1 and λ 1 are respectively associated with driving the system across and along the phase boundary because they are analogous to (i.e. they scale with the shifted exponent of) τ 0 and λ 0 (see Fig. 1a, green arrows and App. B 1). This provides the following interpretation: When the drive is perpendicular to the phase boundary (φ = 0), it is characterized by τ 1 and the scaling is ξ ∼ v −ν1 . In that case we recover the usual KZ scaling. When the system is driven along the phase boundary (φ = π/2), λ 1 takes over, and the irrelevant critical exponent is visible. Crucially, a fine tuning of φ is not necessary, and both scalings can be observed for an extended set of parameters. Indeed, we find that the scaling of ξ with v displays two regimes separated by a crossover drive amplitude v * . The KZ and the new exponents are visible for v smaller and bigger than v * respectively. Notably, φ provides a means to bring either one of these regimes forward because v * interpolates from infinity to zero as φ is tuned from zero to π/2 (See Fig. 3b).
The dependence of the crossover drive amplitude on φ can be estimated for angles close to zero and to π/2 [see Eq. (17)]. This is readily understood by viewing the system as depending on τ 0 with a fixed drive amplitude. For τ 0 large enough (the system is adiabatic and) the correlation length scales as ξ ∼ τ −ν 0 . We see that ξ is bigger for smaller τ 0 . It can however not be arbitrarily large since adiabaticity eventually breaks if τ 0 is decreased while v is held fixed. This happens when the correlation length reaches the smallest of the two underlying scales See [30] where a similar picture emerges. For v small enough, ξ 1 is always smaller than ξ 2 because ν/(1+zν) < 1/(z− ). This means that the usual KZ scaling is always visible when v → 0. A crossover emerges because ξ 1 is bigger than ξ 2 when v is large enough. v * is obtained by equating the two scales See Eq. (14). Its dependence on φ can be estimated for angles close to zero and π/2 (see Fig. 3b) We see that v * vanishes when φ = π/2 and is small when φ is close enough to π/2. The longitudinal scaling emerges when the system is driven at a shallow angle with the phase boundary. Finally, we connect our quasistatic findings to an experimental procedure where the system is driven across the phase transition, taking into account the instantaneous value of the distance from the phase transition τ (t) = τ 0 +v cos(φ)t (see Fig. 1a). Adiabaticity inevitably breaks since τ (t) crosses zero, and the above results are applied at this moment. In particular, v * can be made arbitrarily small by crossing the phase boundary at a shallow enough angle.
The above scaling phenomenology is not bound to the correlation length showcased so far. It will emerge in all typical KZ observables, such as the relaxation time T r ∼ ξ z ∼ v −zν , or the defect density (if these are created as the system crosses the phase boundary), which is evaluated from the correlation volume n d ∼ ξ −d ∼ v dνi (see e.g., [29]). Both carry information on the irrelevant exponent for v > v * .

VII. SPECTRUM OF IRRELEVANT EXPONENTS
In principle, the whole spectrum of irrelevant exponents is accessible through an appropriate slow drive. For example, a higher-order driveĝ =ĝ 0 +ĝ 2t 2 /2 produces a doubly shifted copy of the spectrum of critical exponents: M 2 = M 0 − 2z (see App. A). In particular the newly relevant exponent θ 3 − 2z produces a scaling regime with ξ ∼ v −1/(2z−θ3) and enables the observation of the next irrelevant exponent θ 3 by adjustingĝ 2 longitudinally to the phase boundary. In that case, however, two exponents are made relevant since the first irrelevant exponent is also shifted by 2z, and two irrelevant couplings must be taken into account. The second order drive must follow the phase boundary in a direction that favors the observation of θ 3 over (see App. B 2). As this procedure is iterated, the dimensionality of the required parameter space grows by one for each newly relevant critical exponent. The general principle is as follows: Consider an equilibrium phase diagram with r + 1 axes and an r-dimensional critical surface, i.e., with one direction crossing the transition (equilibrium relevant parameter, transversal direction) and r irrelevant (longitudinal) directions. Then, a polynomial drive of order r can activate the first r irrelevant critical exponents.
We point out that in the case of a transversal polynomial drive of order r, the most relevant operator, which is −y − rz, immediately leads to the known scaling ξ ∼ v −ν/(1+zrν) [27,[60][61][62]. Moreover, our approach can also be applied to the problem of adiabaticity restoration with a symmetry-breaking field (see, e.g., [63,64]). These fields are relevant and produce a breaking of scale invariance with an additional negative critical exponent −y 2 . In particular, driving the system through its critical point with such a field will produce a new scaling with ξ ∼ v −1/(z+y2) .

VIII. SOLVABLE MODEL
In this section we illustrate our result with a noninteracting toy model, which is not particularly realistic but exactly solvable. In particular, we recover our result without using the RG. We consider the fluctuating steady state of the Langevin dynamics in d spatial dimensions, ψ is a real space and time dependent field and ζ a Gaussian white noise with ζ(t, x)ζ(t , x ) = 2δ(t−t )δ(x−x ). We scale the field ψ such that the temperature that usually appears in the noise correlator is set to one. We start with the application of our RG analysis. The renormalization of Eq. (18) is simple because there is no interaction. The flow of m and Q are given by There is a fixed point atm =Q = 0 with one relevant (y = 4) and one irrelevant ( = 2) exponent. Furthermore, the dynamical critical exponent is given by z = 4.
In the presence of a linear drive our work predicts the equilibrium (ν = 1/4), transversal and longitudinal scaling to be respectively. The drive amplitude is given in terms of the dimensionless parameters at the beginning of the RG flow: v 2 = (m 1 /Λ 8 ) 2 + (Q 1 /Λ 2 ) 2 . We now reproduce these results from a less formal, independent argumentation. Although Eq. (18) can be solved exactly with the parameters of Eq. (20), this is actually not necessary. The equal-time correlation function, ψ(t, r)ψ(t, 0) = G(t, r; m 0 , Q 0 , m 1 , Q 1 ) , depends on space r = |(r 1 , r 2 , . . . , r d )|, time t, and the system parameters. We apply dimensional analysis and write the above equation in terms of dimensionless quantities. The different elements of Eq. (22) have the following dimension These dimensions can all be written as powers of a single length L, because the parameter in front of ∇ 4 in Eq. (18) and the temperature have both been set to 1. This dimensional analysis implies that the equation of motion in terms of the rescaled parameterŝ is identical to Eq. (18), with L not appearing explicitly.
In particular this implies that for any choice of L. We now insert the rescaled variables, choose L = r and identify G on both sides of the above equation with T r = r 4 , and G is a dimensionless function that does not depend on the couplings. We see that the four couplings m 0 , m 1 , Q 0 and Q 1 together with the ramp time t, each produce a different scale. Furthermore, the scale produced by the ramp time t/T r = [t 1/4 /r] 4 , is ξ t = t 1/4 . We choose small values of |t| because we are interested in the near critical physics. Close to the fixed point (where all the couplings vanish and |t| is small) we find that, in accordance with the RG analysis, • ξ t is small.
• Q 0 produces a scale that asymptotically vanishes as we approach the fixed point. It does not diverge and therefore is not visible on large spatial scales. Q 0 is an irrelevant parameter.
We now set Q 0 = t = 0 for simplicity and, as in the general analysis, we introduce m 1 = Λ 8 v cos(φ), Q 1 = Λ 2 v sin(φ). Then, the correlation function behaves as a power law (G ∼ r 4−d ) as long as r Min(ξ m0 , ξ m1 , ξ Q1 ) and then decays exponentially. We can therefore identify the correlation length with the smallest of the three observable scales, ξ = Min(ξ m0 , ξ m1 , ξ Q1 ). For a fixed value of φ we extract v * as the value of v where ξ m1 = ξ Q1 . It is given by [v * cos(φ)] 1/8 = [v * sin(φ)] 1/2 . Then, the three regimes that we discuss in our paper emerge naturally: • If m 0 is large enough, ξ m0 is the smallest scale and ξ = m −1/4 0 . The system is adiabatic and exhibits equilibrium scaling.
• If m 0 is small and v < v * , then ξ m1 is the smallest scale and we see the usual KZ scaling ξ ∼ v −1/8 .
• If m 0 is small and v > v * , then ξ Q1 takes over and we see the new scaling ξ ∼ v −1/2 .
We conclude with a remark concerning our choice of model. We have chosen this somewhat unusual equation (with ∇ 4 and ∇ 6 ) because it illustrates our result very cleanly. Indeed, with a more typical model (with ∇ 2 and ∇ 4 ), our RG analysis would still be applicable, but we would have to resort to a quadratic drive because z = .
In this case, the linear drive makes the irrelevant parameter marginal (vanishing critical exponent), not relevant.

IX. CONCLUSION
In an RG language, the KZ mechanism allows one to observe relevant critical exponents by driving along a relevant scaling direction, i.e., transversally to the phase boundary. We find that irrelevant exponents can be made relevant, and therefore observable, by driving longitudinally to the phase boundary. The observability is robust, persisting to the presence of weak transversal drive components. The quantitative difference between the exponents is quite significant if, as usually the case, the full critical exponents are close enough to the canonical ones. It therefore stands to reason that the mechanism uncovered here may underlie some of the difficulties in determining critical exponents in KZ experiments [52], and may help to foster progress in this direction. In this appendix we show how the structure of the stability matrix that is described in the main text emerges from the adiabatic approximation. In particular, we recover Eq. (9), provide an explicit expression for M 0 and M 1 , and identify the role of the adiabatic approximation for each element.
We now formally consider the full RG flow equations without approximation, but use a representation of the flowing parameters that is particularly suited to the adiabatic approximation, In general, the dimensionless RG flow equations take the form Every term in the above equation depends on the rescaled timet explicitly, or throughĝ(t ) and its time-derivatives. The square brackets denote a functional dependence of the form with Eq. (A2) inserted on the right-hand-side. Then, the different blocks of the stability matrix are obtained by differentiating the above equation with respect toĝ i and evaluating atĝ(t ) = g fp , with Eqs. (A2) and (A1) inserted on the right-hand-side. It is easier to evaluate Eq. (A4) if we start with the derivatives with respect toĝ i . Indeed, evaluating the parameters atĝ(t ) = g fp simplifies greatly the time dependence since g fp does not depend on time. We find that the stability matrix is block upper triangular (M ij = 0 if j < i) with the different blocks being square matrices with the dimension of the space of equilibrium parameters. They are We use Latin indices (i, j = 0, 1, . . . ) to denote the different drive sectors. We give an explicit expression for upper-left block M 00 that couples g 0 to itself (equilibrium physics): We use Greek indices [with α, β = 1, 2, . . . , dim(g)] to denote the coordinates within each block. The above equation makes the notation used in Eq. (A5) clear. We see that all the diagonal blocks [referred to as M 0 and M 1 in Eq. (9)] are determined from the equilibrium block We see from Eq. (A5), that the diagonal blocks of M are not affected by the adiabatic approximation. The off-diagonal blocks, however, are because they contain on L(ĝ)[ĝ ] and the derivatives of η[ĝ] with respect to the time-derivatives ofĝ.
In particular, Eq. (B5) becomes which reproduces Fig. 2 whereλ 0 is not represented. We see that | G(k)| increases as k decreases and there is a scale k 0 , below which Eq. (B8) is no longer applicable. Away from the fixed point (and thus on spatial scales larger than 1/k 0 ), the physics is nonuniversal. k 0 then separates the universal and nonuniversal regimes, and can be identified with the inverse correlation length ξ = 1/k 0 . It can be estimated by picking the largest among the three values of k for which the projections of G(k) along v 1 0 , v 1 1 , and v 2 1 are, respectively ±1, 1/ξ ∼ k 0 = Λ Max |τ 0 | 1/y , |τ 1 | 1/(y+z) , |λ 1 | 1/(z− ) .
This reproduces the three scaling regimes that are identified in the main text.

Transversal and longitudinal drives
We now provide some information on the eigenvectors of M and the projection Eq. (B4). We will see that this provides the interpretation of the different scaling regimes as being transversal and longitudinal.

(B9)
A α jk are vectors that depend on v α and the different blocks of M . They can be computed recursively starting from A α jj−1 = [θ α − z − M 0 ] −1 M j−1,j v α . The important element of the above equation is that the j th sub-vector of v α j is given by the equilibrium eigenvector v α .
From Eq. (B9) we find that the highest order part of the drive is a linear combination of the equilibrium eigenvectors. From Eq. (B4) we can extract α λ αr v α =ĝ r , because g fp = (g fp , 0, . . . ) only has nonvanishing components in its equilibrium (i = 0) part. Moreover, the coordinates λ αr do not depend on the lower-order driveŝ g i<r . This is the main result of this appendix. It implies that the equilibrium eigenvectors provide the basis on which to decompose the r th -order part of the drive. If g r is taken to be aligned enough with a given equilibrium eigenvector v α , then the nonequilibrium scaling will be ξ ∼ v −1/(rz−θα) . The usual KZ scaling emerges from a drive along the relevant direction and is therefore transversal. A drive that follows an irrelevant eigenvector will not cross the phase boundary and is longitudinal.