Finite-Time Universality in Nonequilibrium CFT

Recently, remarkably simple exact results were presented about the dynamics of heat transport in the local Luttinger model for nonequilibrium initial states defined by position-dependent temperature profiles. We present mathematical details on how these results were obtained. We also give an alternative derivation using only algebraic relations involving the energy-momentum tensor which hold true in any unitary conformal field theory (CFT). This establishes a simple universal correspondence between initial temperature profiles and the resulting heat-wave propagation in CFT. We extend these results to larger classes of nonequilibrium states. It is proposed that such universal CFT relations provide benchmarks to identify nonuniversal properties of nonequilibrium dynamics in other models.

In [33] two of us (E.L. and P.M.) together with Joel L. Lebowitz and Vieri Mastropietro studied in the Luttinger model the dynamics of heat transport starting from a particular class of nonequilibrium initial states. These states were given by position-dependent temperature profiles 1/β(x) > 0, and the time evolution was determined by the standard translation invariant Hamiltonian at temperature 1/β 0 , and it is therefore natural to interpret 1/β(x) as a position-dependent temperature profile. We stick to this interpretation throughout this paper, although a more common definition of the local temperature would link it directly to E(x, t) neq . Other interpretations of β(x) are possible, in particular if the states in (1.3) arise as equilibria for dynamics defined by inhomogeneous Hamiltonians, as will be briefly discussed at the end of the paper. The setup considered in [33] resembled that of inhomogeneous quantum quenches [34] except that the evolution was studied after quenches from mixed states of (1.3) and it was analyzed directly in real rather than imaginary time.
To study heat transport we were particularly interested in kink-like profiles 1/β(x) interpolating between temperatures 1/β L to the far left and 1/β R to the far right, see Fig. 1. The smooth temperature profile protocol described above allows one to analytically compute the nontrivial behavior of the energy density and the heat current around the location of the kink at early times and their subsequent development into heat waves moving ballistically to the right and left. This should be contrasted with the results of the CFT description of > 0 and b the full system with periodic boundary conditions. Note that, in addition to the kink at x = 0, there is an opposite one at x = ±L/2, which is necessary to have a smooth periodic function. As explained in Sect. 2, the effect of this second kink is eliminated in the thermodynamic limit L → ∞ the dynamics in the partitioning protocol employed in similar previous studies, see [27] and references therein. In such a description, argued to be valid after a transient time, the ballistic heat waves are compressed to simple jumps (shocks) without internal structure moving away from the contact point. In the smooth initial states that we consider, this happens only in the limit when t and x are sent to infinity at the same rate, as such a limit wipes out the nontrivial internal structure of the heat waves. The evolution of the energy density and the heat current obtained in [33] permits then to better understand the shortcomings of the partitioning protocol. It also sheds a new interesting light on transport in integrable systems and, in particular, on how its universal features [35] emerge at long times for a large class of nonequilibrium initial states, see [36] and also Sect. 4.3 below for a related discussion of charge transport in the Luttinger model. As a representative example, we plot E(x, t) neq and J (x, t) neq in Fig. 2 at four times for the Luttinger model with local interactions (defined in more detail below) starting from the kink-like temperature profile in Fig. 1. Note a peak and a dip in the energy density at time t = 0 in the region where the temperature changes and how this local shape evolves into heat waves. This is accompanied by a universal heat current building up in the region between the two heat-wave fronts. We note that, for local interactions, the wave fronts preserve their shapes in time. For the Luttinger model with nonlocal interactions, there are additional dispersive effects, which, however, eventually become unobservable in any finite region as the wave fronts leave such regions in finite time, see Fig. 1 in [33]. In the remainder of this paper we restrict our discussion to the local case.
]/J and the heat current for different times in the nonequilibrium state with the inverse-temperature profile given by (2.1) and (2.3). The parameters are β L = 19.9, β R = 20.1, δ/ = 0.06, and v/ = 0.025 In the local Luttinger model, the energy density is given formally by 2 (1.5) with fermionic field operators ψ ± (x) (with antiperiodic boundary conditions) obeying the usual canonical anticommutation relations {ψ r (x), ψ † r (y)} = δ r,r δ(x − y), etc., :· · ·: denoting Wick (normal) ordering, the bare Fermi velocity v F > 0, and the coupling strength λ > −πv F /2. The (diverging) constant E 0 subtracts the ground-state energy density up to the finite contribution −πv/(6L 2 ) left for later convenience. The heat current is given by [36,39]. We found in [33] the following exact expectation values of the energy density and the heat current in the thermodynamic limit marked by the superscript ∞ on the expectations: 3 where x ± = x ± vt are the light-cone coordinates and the function is determined by the temperature profile. We also observed that F(x) can be written in terms of the Schwarzian derivative of the function for some constant β 0 > 0. Noting that Schwarzian derivatives appear in CFT [40] and that the local Luttinger model is a CFT, we also argued that it should be possible to obtain the result in (1.7) and (1.8) in a simpler way using conformal transformations, and that it may be possible to also obtain expectation values of other observables in that way [33]. Our main result in this paper is to show that this is indeed the case. The method used in [33] was perturbative in a small parameter measuring the distance to equilibrium in the initial state (i.e., the case = 0 corresponds to the Gibbs state). This method is general, but generically one can only obtain useful low-order results. In the special case of local interactions, however, we were able to push the computations to all orders in , and, summing the resulting infinite series, we obtained the results in (1.7) and (1.8). In this paper we give a simpler derivation of these results extending them to all unitary CFT models, including models at finite L, and to other observables.
Our analysis is based on the Minkowskian version of CFT, and it generalizes to other classes of nonequilibrium states. One such generalization is (1.3) but where, in addition to the temperature profile 1/β(x), we also allow for a "velocity" profile ν(x) taking (1.12) for |ν(x)| < v (this condition ensures that G is positive). Note that for constant profiles this is a generalization [41] of the Gibbs state in (1.4) where β 0 H is replaced by β 0 (H + ν 0 P) with the momentum operator This may be also be viewed in CFT as a Lorentz boost of the ordinary Gibbs state [27]. As will be shown, the corresponding generalization of the result in (1.7) is obtained by replacing F(x ∓ ) by F ± (x ∓ ) given by (1.8) with the right-hand side multiplied by the central charge c of CFT (which is equal to 1 for the local Luttinger model) and with β(x) replaced by (1.14) In particular, in the long-time limit where β ±,L and β ±,R are the asymptotic values of β ± (x) to the left and to the right, respectively. 4 Finally, for CFTs with a double U(1) current algebra (e.g., the Luttinger model itself), we may also handle chemical-potential profiles μ ± (x) in addition to temperature profiles 1/β ± (x), possibly different for right and left movers. In this case, the functions F ± (x) pick up an additional term proportional to μ ± (x) 2 . This is the largest class of nonequilibrium states that we consider. They are defined as in (1.3) but with G in (4.40), see Sect. 4.3 for details.
The above results show that in CFT there is a universal relation between the initial temperature, velocity, and chemical-potential profiles and the resulting heat-and density-wave propagation even at finite times. 5 We believe that this provides a useful benchmark for other models as follows. Typically, finite-time results are model-dependent, and more universal behavior is only obtained at long times [27]. As an example, we mention the nonlocal Luttinger model which exhibits finite-time dispersion effects depending on short-distance details of the interaction potential [33]. However, these effects have some qualitative features that are always present. We postulate that the CFT results have to be subtracted in order to identify the effects that come from the microscopic details in the propagation of the heat or density waves emanating from the inhomogeneities of the initial state. In addition, such a subtraction should allow to identify the time scales when such model-dependent effects are important and when not. Moreover, it was argued in [19] that integrable systems come in two kinds: those that are purely ballistic and those with a nonzero diffusive contribution. For heat transport, in particular, one way to view this is through the thermal conductivity in the frequency domain where a nonzero thermal Drude weight D th indicates the presence of a ballistic contribution and a finite nonzero value of the regular part Re κ reg th (ω) at ω = 0 signals the presence of a diffusive component in the heat transport [9]. It is known that "pure" CFT captures the ballistic part [27], with Re κ reg th = 0. This also follows from our results. The cause of dispersive and diffusive effects thus must come from short-distance details, randomness, or other relevant perturbations, see, e.g., [42] for a recent discussion of that issue within CFT.
The plan of the rest of this paper is as follows. In Sect. 2, we sketch the original derivation of the result in (1.7) and (1.8) and explain the physical significance of the limit L → ∞ since this is also relevant for our more direct CFT argument. We also present special integrals whose exact evaluation was the key to this result and which, as we believe, are interesting in their own right. The reader may skip the second half of Sect. 2 without loss of continuity. The CFT derivation is given in Sect. 3. After collecting the results about Minkowski-space CFT that are needed in Sect. 3.1, we show in Sect. 3.2 how to use conformal transformations to straighten out position-dependent temperature profiles 1/β(x) on a periodic interval. This allows one to exactly map the nonequilibrium expectations in CFT to the corresponding equilibrium ones. We make this mapping explicit for products of the components of the energy-momentum tensor. In Sect. 3.3, we study the thermodynamic limit of the finite-volume relations which allows one to treat temperature profiles on the infinite line with different asymptotic values on the left and right sides. As a byproduct, we calculate the thermal Drude weight. Section 4 is devoted to various generalizations. In Sect. 4.1, we briefly discuss the correlators of primary fields. In Sect. 4.2, we consider states with different temperature profiles for the right and left movers. They form a class of nonequilibrium states preserved by the Schrödinger-picture dynamics that lead to simple examples of generalized Gibbs states at long times. In Sect. 4.3, we extend the analysis to CFTs with a U(1) current algebra and states with temperature and chemical-potential profiles. For the Luttinger model we discuss how this implies universality of conductance for both the charge and axial currents, generalizing previous results in [36], see also [43]. Finally, in Sect. 5, we make contact with [29][30][31][32], discussing the dynamics preserving states in (1.3) and the related Euclidian CFT description. We end with conclusions and directions for further developments in Sect. 6. The Appendix contains some mathematical details on the special integrals mentioned above.

Perturbative Derivation and Remarkable Integrals
The perturbative computation method used in [33] is based on introducing an expansion parameter measuring the deviation of the temperature profile from constant temperature 1/β as follows: with W (x) a function defined by this relation.
In [33] we were mainly interested in kink-like functions where W (x) becomes (say) 1/2 and − 1/2 to the far left and right, respectively, so thatβ = (β L +β R )/2 and = (β L −β R )/β in terms of the asymptotic values. On the other hand, for technical reasons, we used a model on a circle with circumference L < ∞, which at first sight seems incompatible: a smooth function W (x) on the circle with a single kink is not possible, and there has to be at least one other opposite one. As an example, consider the periodic function which is kink-like in the vicinity of x = 0 as desired but has an opposite kink in the vicinity of x = ±L/2 and leads (for negative ) to the temperature profile of Fig. 1. The effect of the additional step around ±L/2 can be eliminated by computing results O(t) neq for L < ∞ at finite times t, and then taking the limit L → ∞ [33] in which (2.2) turns into In this way technical problems with L = ∞ on the level of quantum field theory are avoided and, at the same time, boundary conditions have no influence on the final results. The physical interpretation is that the kink at x = ±L/2 is "behind the moon" and does not affect the physics in any finite region at times significantly smaller than L/v, which is a time scale that becomes infinite in the limit L → ∞.
In the rest of this section we sketch the perturbative derivation of the result in (1.7) and (1.8), concentrating on remarkable integrals which were the key to this results. The readers mainly interested in our CFT derivation may pass directly to Sect. 3 without loss of continuity.
With the inverse-temperature profile (2.1), one can use the Dyson series to obtain an for any local observable O(x, t), such as E(x, t) and J (x, t). The leading term in (2.4) is the equilibrium expectation value (which is time independent, i.e., O(x, t) β = O(x) β ), and the n-th order term is an (n + 1)-point correlation function for n = 1, 2, . . . [33]. This method works, in principle, for any model but, in practice, it is difficult to go beyond leading order n = 1. For a quasi-free bosonic model, to which the Luttinger model reduces, one can use general mathematical results [44][45][46] to replace the many-body computation by a much simpler one-particle one, and this makes it possible to extend the calculation to all orders [33]. In particular, for the local Luttinger model, after taking the limit L → ∞, this method gives (1.7) with It is interesting to note that the limit L → ∞ is not only useful to eliminate the effect of boundary conditions (as explained) but also computationally: this limit turns Riemann sums into integrals and eliminates zero-mode contributions which would be more difficult to handle [33]. The integrals in (2.7) are certainly nontrivial, but we found that they all can be computed exactly, giving the result in Lemma 2.1, and this was the key that led to (1.8).

Lemma 2.1
For all n = 1, 2, . . ., A proof can be found in the Appendix.
It is remarkable that the result are even second order polynomials in the variables q j . It is clear from (2.6) that the constant term leads to contributions to F n (x) which are proportional to W (x) n , whereas the terms with q 2 j and q j q k , j = k, lead to W (x) n−2 W (x) and W (x) n−2 W (x) 2 , respectively. Thus, the special form of the integrals in Lemma 2.1 implies that we get at most terms involving second derivatives of W (x). The explicit expression of these integrals allows one to compute F n (x) exactly, and the result is simple enough to analytically sum the series in (2.5), which gives the result in (1.8) [33].
We describe this computation above since it allows one to interpret the argument in the next section as a partial proof of Lemma 2.1. Such a proof is only partial since all terms with q 2 j are identified with (say) q 2 1 , and all terms with q j q k for j = k are identified with q 1 q 2 . This identification is also useful in order to explicitly derive (1.8) from Lemma 2.1, see Eq.
(A4) in [33], which is implied by (2.9). Thus, the exact integrals in Lemma 2.1 contain more information than the result in (1.8). Since nontrivial integrals that can be computed exactly are rare and often not only have one but several applications in physics, we prove (2.9) in this paper. Moreover, since our derivation of (1.7) and (1.8) in the next section works even for finite L, it suggests interesting Riemann sum generalizations of the exact integrals in Lemma 2.1. We believe it would be interesting to work them out, but this is left for a future study.
We finally mention that our argument in the next section allows one to interpret the computation described above as a derivation of the conformal anomaly in CFT by a direct computation, see (3.7).

CFT in Minkowski Space
We consider the Minkowskian version of a unitary two-dimensional CFT where space is a circle S 1 parameterized by the periodic coordinate x with the basic range −L/2 ≤ x ≤ L/2 and where t ∈ R is time. We keep the propagation speed v in our equations to clearly indicate how it effects the (otherwise) universal law relating the temperature profiles to the heat-wave dynamics.
The basic objects of such a CFT are the periodic light-cone components T ± (x ∓ ) = T ± (x ∓ + L) of the energy-momentum tensor 6 where, as before, x ± = x ± vt. They are distributions with values in the self-adjoint operators on the Hilbert space of states of the theory that satisfy the equal-time commutation relations where δ(x) stands for the L-periodized δ-function. The real number c is the central charge of the theory. In terms of the Fourier modes, the commutation relations in (3.1) and (3.2) reduce to those of the Virasoro algebra, Technically, we assume that the Hilbert space of the theory is a (possibly infinite) direct sum of unitary highest-weight representations of two commuting copies of the Virasoro algebra. The local Luttinger model is an example with c = 1 of such a theory where with the fermion densities ρ ± (x) = :ψ † ± (x)ψ ± (x): and tanh 2ϕ = −λ/(πv F +λ). A related quantity describing the interactions is the Luttinger parameter 7 K = e 2ϕ . The effective densitiesρ ± act in a direct sum of bosonic Fock spaces that contains the interacting vacuum | and the Wick ordering in (3.5) is with respect to | [37,38]. In the following arguments, the explicit form of the operators T ± (x) is not used.
Let Diff + (S 1 ) denote the covering group of the group Diff + (S 1 ) of orientation-preserving diffeomorphisms of the circle. The elements of Diff + (S 1 ) are represented by smooth func- and f (x) + nL corresponding to the same diffeomorphism in Diff + (S 1 ). The operatorvalued distributions T ± generate two commuting projective unitary representations U ± of Diff + (S 1 ) on the Hilbert space of the theory such that for infinitesimal diffeomorphisms and under the adjoint action with the Schwarzian derivative (S f )(x) given by (1.9). This was proven in [47] for the unitary highest-weight representations of the Virasoro algebra and carries over to the present context.
The energy-momentum tensor determines the Hamiltonian, and under the Heisenberg picture evolution, T ± (x, t) = T ± (x ∓ ) as claimed above. The energy and momentum density operators satisfy the continuity equations As we shall see, the finite-volume nonequilibrium expectation values in (1.
are well defined for any such CFT provided the inverse-temperature profile is periodic.

Relating Nonequilibrium to Equilibrium Expectations
From the above, it is clear that the calculation of the time evolution of the nonequilibrium expectation values of the energy density and current operators is equivalent to computing The key observation is that it is possible to find an f such that for some constant β 0 > 0. In order to see this, take the function f given by (1.10) with the constant β 0 determined by The above choice of β 0 ensures that f (x + L) = f (x) + L and thus that f defines an element in Diff + (S 1 ). Using this function f , it follows from (1.3), (3.7), and (3.9) that (3.13) Upon using the relation f (x) = β 0 /β(x) for the derivative of f , this reduces to where the last equality follows by the change of variables y = f (x). This establishes (3.11).
In short, the conjugation with U ( f ) straightens out the temperature profile. From the definition in (1.3) and the relation in (3.11), we infer that This translates the nonequilibrium expectations to the equilibrium ones defined by (1.4). For (with noncoincident points), using again (3.7), we obtain the relation for f given by (1.10) and (3.12). In particular, The Gibbs state is translation invariant so that the expectations T ± (y) β 0 do not depend on y, but they depend, in general, on vβ 0 and L. By scaling, however, (vβ 0 ) 2 T ± (y) β 0 depends only on vβ 0 /L but in a way dependent on the representation content of the CFT. As we shall see, what is universal, depending only on the central charge, is the L → ∞ limit of (vβ 0 ) 2 T ± (y) β 0 .

Thermodynamic Limit
Let us consider the limit L → ∞ of the nonequilibrium expectations in (3.17). For a large class of kink-like β(x) profiles (not necessarily symmetric) with an antikink around ±L/2, where β L and β R are the asymptotic values of the plateau to the left and to the right of the kink. 8 Similarly, for fixed x, β(x) will stabilize up to O(L −1 ) terms with any trace of the antikink gradually wiped out, and so does the function f given by (1.10). The question about the large-L limit of the nonequilibrium expectations in (3.17) then boils down to the one for the equilibrium expectations j T r j (x j ) β 0 , where, by rescaling, β 0 may be set to its asymptotic value and the insertion points are allowed to have O(L −1 ) variations.
In CFT the control of the thermodynamic limit of the equilibrium expectations is an easy exercise. The thermal expectations like those mentioned above may be viewed as the ones in the Euclidean theory on the torus S 1 × S 1 where the circles have circumferences L and vβ 0 , respectively. In a modular invariant CFT [40], they also have a dual representation (3.20) as the equilibrium expectations with inverse temperature L/v in the theory on the circle with circumference vβ 0 . The components of the energy-momentum tensor with complex arguments on the right-hand side are defined by and T orders x j increasingly from the right to the left. The identity in (3.20) comes from swapping the two circles that play a symmetric role in the Euclidean version of the theory. It still holds for the Luttinger model if the finite-volume theory corresponds to antiperiodic boundary conditions for the fermionic fields ψ ± (x), as we assumed, even if the corresponding CFT does not have full modular invariance. When L → ∞, the right-hand side of (3.20) tends to the vacuum expectations providing the dual representation of the equilibrium expectations in the thermodynamic limit: Besides, such vacuum expectations are universal because they receive contribution only from the tensor product of the two vacuum highest-weight representations of the Virasoro algebra. They factorize according to In the theory on the circle with circumference vβ 0 , (3. 26) We infer that the identity in (3.17) holds in the thermodynamic limit for infinite-volume profiles β(x) > 0 with arbitrary asymptotic values β L and β R and the function f defined by (1.10). By scaling, the right-hand side is then independent of the choice of β 0 . In particular, (3.18) together with (3.24) show that in the limit L → ∞, (3.27) This is equivalent to (1.7) with F given by the right-hand side of (1.8) multiplied by the central charge c.
The above result allows one to easily extract the value of the thermal Drude weight D th , see (1.16), which may be obtained [18,48] from where the nonequilibrium expectation is calculated for the inverse-temperature profile β(x) interpolating between the asymptotic values β L and β R with β = β L −β R . The factor −β 2 0 is there to relate D th to the response to temperature rather than inverse-temperature change.
The space integral receives the contribution from the region of length ≈ 2vt between the two ballistically separating heat waves where the heat current takes the long-time value which is proportional to the temperature 1/β 0 . The result may be also obtained using the partitioning protocol that leads to the same steady state. It agrees with the calculation of the thermal conductivity by the Green-Kubo formula, where we used (3.25) and (3.26) to express the infinite-volume 2-point correlation function of the heat current. Note that the regular part of Re κ th (ω) vanishes, confirming the absence of diffusive heat transport in nonequilibrium CFT.

Other Correlators
The finite-volume relation in (3.15) between the nonequilibrium and equilibrium expectations may be rendered explicit also for observables where j (x − j , x + j ) are primary fields whose transformation laws under the Diff + (S 1 ) symmetry take the form for f given by (1.10) and β 0 by (3.12). Note that these are simpler relations than for the energymomentum tensor components T ± since those fail to be primary fields with conformal weights (2, 0) and (0, 2) due to the Schwarzian derivative term in (3.7) reflecting the conformal anomaly. In a similar way as for observables built from the operators T ± , the relations (4.3) hold also in the thermodynamic limit which may be controlled like before. Also, as before, β 0 may be taken arbitrary in the infinite volume. For example, the Luttinger model has a conserved U(1) current with the light-cone components J ± (x ∓ ) = √ Kρ ± (x ∓ ), see (3.5), with conformal weights (1, 0) and (0, 1), respectively, and (renormalized) fermionic fields ψ ± (x − , x + ) with conformal weights (4.4) accompanied by their Hermitian conjugates ψ † ± (x − , x + ) with the same conformal weights. 9 Their infinite-volume equilibrium 2-point correlation functions have the form and Thus, it follows from (4.3) that the corresponding nonequilibrium correlation functions are and (4.8) We note that the latter agrees with Eq. (19) in [33] to first order in the expansion parameter in (2.1) and exactly reproduces Eq. (10) in [33] in the long-time limit.

Temperature and Velocity Profiles
It is straightforward to generalize the argument of Sect. 3.2 to nonequilibrium states as in (1.3) with G given by (1.12), which is more conveniently written as Our argument above goes through as it stands but with Choosing them as one straightens out both profiles 1/β ± (x), replacing (3.15) by the identity (4.12) define the expectations in a simple example of a generalized Gibbs state with different temperatures for the right and the left movers. In particular, (3.17) becomes (4.13) The thermodynamic limit of the expectations in (4.12) of the observables in (3.16) may still be conveniently studied by going to the dual picture which, upon setting β 0,± = β 0 (1±ν 0 /v), takes the form 10 where on the right-hand side the expectation is in the theory on the circle with circumference vβ 0 . We infer that in the thermodynamic limit where on the right-hand side the first (second) vacuum expectation is in the theory on the circle with circumference vβ 0,+ (vβ 0,− ), and the last equality follows by the factorization in (3.23) of the vacuum expectations and their rescaling. In a similar way as in Sect. 3.3, this shows that the identity in (4.13) holds in the thermodynamic limit with the infinite-volume profiles β ± (x) > 0 with arbitrary positive asymptotic values β ±,L and β ±,R and f ± (x) given by the first of the equations in (4.10) with arbitrary β 0,± > 0. In particular, (4.16) implying that with the functions There is a conceptual gain from the consideration of the nonequilibrium states with different temperature profiles for the right and left movers: unlike the states with equal profiles, such states are preserved by the Schrödinger-picture dynamics. Indeed, for G given by (4.9), (4.19) so that under the Schrödinger-picture evolution the profiles β ± (x) move ballistically to the right and to the left, respectively. This still holds in the thermodynamic limit and makes it clear why for long times such states converge to the generalized Gibbs state in (4.12) with inverse temperatures β 0,+ = β +,L and β 0,− = β −,R . Conversely, we may view the nonequilibrium states with G given by (4.9) as the generalized Gibbs state with local profiles β ± (x) whose time evolution under the Schrödinger-picture dynamics reduces to the time evolution of the profiles governed by the equations This is reminiscent of the time evolution of the generalized Gibbs states with local profiles in the generalized hydrodynamic picture of integrable models out of equilibrium [14][15][16][17][18][20][21][22][23]. In CFT, however, no hydrodynamic-scale closure is needed to obtain the hydrodynamic evolution equation (4.20) [49].

Temperature and Chemical-Potential Profiles
Suppose that our CFT possesses a U(1) symmetry generated by a conserved current with components ρ and j satisfying ∂ t ρ + ∂ x j = 0. (4.21) For instance, ρ = ρ + + ρ − and j = K v(ρ + − ρ − ) for the local Luttinger model, see, e.g., Appendix B in [36] (we recall that K is the Luttinger parameter). One may then consider states with both temperature and chemical-potential 11 profiles by taking (4.22) where ρ(x) is the zero time density and μ(x) is a periodic chemical-potential profile. Suppose that the light-cone components J ± = (1/2)(ρ ± v −1 j) of the current depend only on x ∓ , respectively, and satisfy the U(1) current algebra: In terms of the Fourier modes, (4.25) and the commutation relations take the form (4.26) For concreteness, we shall restrict our discussion to the case of the local Luttinger model in which case J ± = √ Kρ ± and κ = K , see (3.5), or to the case when J ± is one of the Cartan subalgebra components of the two current algebras of the level k WZW theory based on a compact Lie group [40] (e.g., U(1) or SU(N )) in which case κ = k/2. It follows from [47,50] that in these cases there exist, besides the two projective representations U ± of Diff + (S 1 ) considered above for which (4.27) also two commuting projective representations V ± of the additive gauge group of periodic smooth maps h(x) = h(x + L) on R generated infinitesimally by J ± , for h(x) = εξ(x). Under the adjoint action of V ± (h), (4.30) Recall from Sect. 3.2 that the conjugation with operators U ( f ) = U + ( f )U − ( f ) for f ∈ Diff(S 1 ) satisfying (1.10) with β 0 given by (3.12) was previously used to straighten out the periodic inverse-temperature profile β(x) to a constant one given by β 0 . We shall keep the function f as before and choose Then conjugating G in (4.22) first with V (h) = V + (h)V − (h) and then with U ( f ) straightens out the chemical-potential and inverse-temperature profiles to μ 0 and β 0 , respectively, leading to the identity (4.33) The thermodynamic limit can be controlled as before using the dual representation of the equilibrium expectations (that involves now the theory on a circle of circumference vβ 0 with twisted boundary conditions). In the infinite volume, one can again treat profiles with arbitrary asymptotic values. Using the fact that in the theories under consideration, in the thermodynamic limit, one obtains the identities in (1.7) with as well as the formulas For the matrix D of the Drude weights obtained from the nonequilibrium expectations with respect to states with profiles with small kinks of heights β = β L −β R and μ = μ L −μ R around the constant values β 0 and μ 0 of the inverse temperature and chemical potential, the formula in (3.28) generalizes to (4.39) Note that the thermal Drude weight D 11 is independent of μ 0 and the density Drude weight D 22 is independent of the temperature. In particular, for the Luttinger model, The lack of symmetry of D is due to the asymmetric way in which the temperature and the chemical potential enter into the Gibbs state. 12 As before, one may also consider nonequilibrium states with different local profiles β ± (x) and μ ± (x) for right and left movers, which are defined by (1.3) with (4.40) This leads to the replacement of the functions F and G in (4.35) and (4.37) with and respectively. In this case, the expectation values are given by (4.17) for E(x, t) ∞ neq and J (x, t) ∞ neq as well as The nonequilibrium states in (4.40) form the family of generalized Gibbs states with local profiles that correspond to the commuting conserved charges H ± of (4.12) and Such a family of states is again preserved by the Schrödinger evolution that displaces the local profiles β ± (x) and μ ± (x) ballistically. Their time evolutions are governed by (4.20) together with ∂ t μ ± ± v∂ x μ ± = 0. (4.44) In the long-time limit, one obtains a genuine (i.e., with constant profiles) generalized Gibbs state which is the thermodynamic limit of the state with In particular, generalizing (1.15). As a specific example, consider the Luttinger model with the initial state given by (4.40) with β ± (x) = β 0 and with μ ± (x) interpolating when L → ∞ between the asymptotic values μ ±,L and μ ±,R . The charges Q J ± = √ K ρ ± (x)dx are then the number operators of right and left moving plasmons. One can also consider in that case the number operators of electrons (and holes) Q e ± = ρ ± (x)dx, which are different from Q J ± , although the total charges are equal, i.e., Q e + + Q e − = Q J + + Q J − . Indeed, we infer from (3.5) that (4.50) Unlike for Q J ± , the spectra of Q e ± are composed of integers. In terms of Q e ± , the generalized Gibbs state appearing in the long-time limit of the evolution will correspond to From (4.49) we infer that the value of the permanent current in the limiting nonequilibrium steady state is It was argued in [51][52][53] that μ e ± that couple in the steady state to the electron charges correspond to the chemical potentials of free electrons of wide leads connected to a Luttinger wire, at least if μ + (x) = μ − (x). As was discussed in [36], the second equality of (4.52) would then provide an explanation for the experimental measurements [54] of conductance in quantum wires that gave results close to the universal constant e 2 / h, equal to 1/2π in the units h = e = 1 that we are using. This universal value is different from the conductance K e 2 / h predicted in [55] which, instead, is consistent with the first equality of (4.52) that uses the asymptotic values of the imposed chemical-potential profiles that couple in the steady state to the plasmon charges. The above extends the derivation of the universal result obtained in [36] to states with constant temperature and with chemical-potential profiles possibly different for the right and left movers. 13 The Luttinger model possesses also a conserved axial current with ρ A = ρ + − ρ − and j A = (v/K )(ρ + + ρ − ) satisfying ∂ t ρ A + ∂ x j A = 0. Note that (4.48) implies that the permanent axial current in the limiting nonequilibrium steady state considered above takes the value with the universal coefficient both when expressed in terms of the asymptotic values of the profiles μ ± (x) and in terms of μ e ± .

Equilibrium Dynamics and Relation to Euclidian CFT
In this section, we discuss the relation between the articles [29][30][31][32] and the present paper. In [29] it was argued that the kernel of the 1-particle density matrix in the ground state of a nonrelativistic high density Fermi gas in a trap may be described on mesoscopic scales by the 2-point function of the fermionic massless free field whose Fermi velocity varies in space. These results were generalized in [30][31][32] to certain nonrelativistic systems of interacting 1d bosons in traps. Despite similarities, there are several differences with the approach of the present paper. First, the arguments in [29] were based on the analysis of the ground-state Euclidiantime correlators in the presence of a trap and these were shown to correspond to Euclidian CFT correlators in an appropriate curved metric (in [32] also the coupling to a gauge field appeared implicitly). In this paper, we consider positive temperatures, but the correspondence of [29] generalizes to low-temperature states leading to the compactification of the Euclidian time direction in CFT, just as for homogeneous equilibria. Hence, for specific CFTs, the states in (1.3) may, indeed, be viewed as describing on mesoscopic scales the 1d nonrelativistic lowtemperature matter in traps, with β(x) having the interpretation of the position-dependent Fermi velocity in appropriate units. Second, the argument of [29] was done for the equilibrium dynamics (although some nonequilibrium situations were also considered), whereas in the bulk of the present paper we study the dynamics generated by the homogeneous Hamiltonian that does not preserve the states in (1.3). Our considerations may, however, be generalized to dynamics induced by inhomogeneous Hamiltonians [56], in particular to the ones that preserve the states in (1.3). Third, we use the Minkowski version of CFT, whereas the papers [29][30][31][32] employed the Euclidian CFT. That is usually considered as an innocent distinction handled by the Wick rotation. Indeed, for the primary fields, the correlators in the Euclidian theory in the metric considered in [29][30][31][32] and with compactified time do agree, up to the Wick rotation, with the corresponding correlators in the states of (1.3) for which the time dependence is generated by the equilibrium dynamics. As shown below, however, that does not hold directly for the correlators of the energy-momentum tensor components which are of main interest in this paper. This points to the need of caution when one applies Euclidian techniques in the study of systems that are inhomogeneous in space or/and time.
To be more concrete, let us briefly discuss the equilibrium dynamics for the states in (1.

3) that is generated by the inhomogeneous HamiltonianH
we immediately obtain from (3.15) and (3.11) the relations Note the difference of the right-hand sides with those of (3.17) and (4.3) corresponding to the dynamics generated by the homogeneous Hamiltonian H of (3.8) which results in the time dependence in the arguments of function f . After the Wick rotation t j = −iτ j , the correlators in (5.3) become the correlation functions of the same primary fields in the Euclidian CFT on the torus S 1 × S 1 parameterized by (x mod L , τ mod β 0 ) and equipped with the Riemannian metric where h 0 = dzdz for the complex coordinate z = f (x) + ivτ on the torus with f (x) given by (1.10) and β 0 by (3.12), and where σ (z,z) = −2 ln f (x). In other words, where on the left-hand side is the Wick rotated (5.3) and on the right-hand side the Euclidian correlation functions on the torus S 1 × S 1 with the Riemannian metric h of (5.4). Such a relation is well known for f (x) = x. Its generalization to general f follows from the identity j j (z j ,z j ) The ground-state version of the relations in (5.5) provided the basis for the use of Euclidian CFT in the description of the trapped 1d fermions or bosons on mesoscopic scales in [29][30][31][32], with the interpretation of vβ(x)/β 0 as the position-dependent Fermi velocity clearly reflected in the form of the metric in (5.4). Let us pass to the discussion of the energy-momentum correlators. For f (x) = x, the Wick-rotated correlators in (5.2) are represented by the Euclidian correlation functions of the energy-momentum components T + = T zz (dz) 2 and T − = Tzz (dz) 2 on the torus S 1 × S 1 with metric h 0 . However, for general f (x) one has [57], in the notation z + = z, z − =z, j T r j (z j ,z j ) and the right-hand side does not represent correctly the Schwarzian-derivative terms of the Wick-rotated (5.2). In the Euclidian domain, the Schwarzian derivative appears in the transformation law of the energy-momentum components when one deals with holomorphic transformations z → f (z) [40,57], but this is not the case here. A closer examination shows that where R(z,z) = −β (x)/β(x) is the scalar curvature of the metric e σ h 0 . This substantiates the comment made above about the need of caution.

Conclusions
We elaborated on the formula of [33] giving the full time evolution of the energy density and heat current from a nonequilibrium state with a preimposed temperature profile in the Luttinger model with local interactions. The formula was obtained in [33] by expanding the nonequilibrium state around the equilibrium to all orders. More details on the perturbative computation involving the exact calculation of complicated integrals, that may be interesting in its own right, were given. The main part of the paper was devoted to showing how the formula of [33], a result of the resummation of the perturbative series, may be obtained using Minkowskian conformal symmetries of the local Luttinger model. The idea was to use conformal transformations to map spatially inhomogeneous situations to homogeneous ones, straightening out a nonuniform temperature profile to a constant one. This led to a direct relation between nonequilibrium and equilibrium states, yielding the remarkable formula of [33] as a corollary. The CFT argument holds for a general class of unitary CFTs and could be applied to a wider class of nonequilibrium states that are preserved by the Schrödingerpicture evolution. The states in this class may be viewed as particular examples of simple generalized Gibbs states with local profiles, and they tend to ordinary generalized Gibbs states at long times, somewhat similarly as in the scenario recently advocated for integrable models where the evolution at certain length and time scales could be described by generalized hydrodynamics [14][15][16][17][18][20][21][22][23]. We obtained similar results also for CFTs with a U(1) current algebra (including the local Luttinger model itself) where we treated nonequilibrium states with temperature and chemical-potential profiles. Moreover, our results permit a more detailed analysis within CFT, compared to using the partitioning protocol studied before [27], of how a system starting in a state that looks like two different equilibria joint together evolves in time towards a nonequilibrium steady state described by a generalized Gibbs state.
As was discussed in Sect. 5, at least some families of the CFT nonequilibrium states that we studied in the present paper could be interpreted as providing a mesoscopic-scale description of dense nonrelativistic 1d matter in macroscopic traps [29][30][31][32]. The dynamical correlators in such CFT states should similarly describe the corresponding nonrelativistic correlators at mesoscopic time scales both for after-quench and for equilibrium dynamics. The other way of arriving at the family of nonequilibrium CFT states that we considered is by reversing the logic of this paper. In the periodic-volume Minkowski CFT, the conformal symmetries (together with the gauge symmetries if a U(1) current algebra is present) are broken in the usual equilibria that are not preserved by the symmetries. Instead, the application of the symmetry transformations to the equilibrium states generates the family of nonequilibrium states that were studied here.
In the infinite volume, the states with kink-like profiles give access to the full counting statistics of the energy or charge transfers through the kinks, similarly to the states arising in the partitioning protocol [27]. Although such statistics in both approaches differ at finite times, they have the same long-time large deviations. This will be discussed elsewhere as it requires using different boundary conditions for finite volumes that allow one to avoid the duplication of kinks in the profiles. Finally, another interesting exercise, which was abundantly discussed in the similar context of quantum quenches [25], concerns the evolution of the entanglement entropy or negativity starting from states with profiles of the type consider here. By the replica trick, the latter may be extracted from nonequilibrium correlators of the twist primary fields in the replicated theory, to which our approach gives direct access. The analysis of the corresponding formulas is left for future research.
To compute these integrals we insert the Taylor series f ν ( p + Q j ) = f ν ( p) + f (1) j+1 (A.2) into the integrand and obtain In the second equality we wrote (iν) n−m 0 +1 = (iν − vp + vp) n−m 0 +1 to expand into a binomial series, and in the third we summed the bosonic Matsubara frequencies ν using the Mittag-Leffler series of coth(vβ p/2) and changed variables to s = vβ p. We note that, for M ≥ 1, the integrand of the last integral is singular, but the singularity is removable, i.e., one can replace coth(s/2) by coth(s/2) − 2/s without changing the result. To further compute these integrals we use partial integrations. We find that the integrals in (A.6) are zero for M ≥ 3, which implies R n = 0. The remaining integrals I