Hydrodynamic effective field theory and the analyticity of hydrostatic correlators

We study one-loop corrections to retarded and symmetric hydrostatic correlation functions within the Schwinger-Keldysh effective field theory framework for relativistic hydrodynamics, focusing on charge diffusion. We first consider the simplified setup with only diffusive charge density fluctuations, and then augment it with momentum fluctuations in a model where the sound modes can be ignored. We show that the loop corrections, which generically induce non-analyticities and long-range effects at finite frequency, non-trivially preserve analyticity of retarded correlation functions in spatial momentum due to the KMS constraint, as a manifestation of thermal screening. For the purposes of this analysis, we develop an interacting field theory for diffusive hydrodynamics, seen as a limit of relativistic hydrodynamics in the absence of temperature and longitudinal velocity fluctuations.


Introduction
| 3 from a hydrostatic generating functional, W = T −1 0 d d x √ −g F(g µν , A µ ), expressed in terms of the background metric g µν and gauge field sources A µ . Here T 0 is the constant temperature of the global thermal state and d is the number of spatial dimensions. To wit, T µν eq δg µν + J µ eq δA µ . (1.1) The associated conservation equations ∇ µ T µν = F νλ J λ and ∇ µ J µ = 0 follow directly from the diffeomorphism and gauge invariance of W. An important ingredient in implementing the ensuing constraints on the hydrodynamic constitutive relations is the finiteness of the static correlation length. This feature is expected for generic thermal systems that are not at second order critical points and do not have spontaneously broken symmetries, and implies that retarded hydrostatic (time-independent) correlation functions should fall-off exponentially at large distances. Equivalently, the Fourier space zero-frequency retarded correlation functions must be analytic in an expansion in small spatial momentum. The generating functional thus admits a derivative expansion, and one can obtain the hydrodynamic constitutive relations in the static limit order-by-order in this expansion.
Moving away from thermal equilibrium, the full EFT formalism of hydrodynamics should capture the physics of small fluctuations about the hydrostatic equilibrium state. Thus, hydrodynamic correlation functions should reduce to those derived from the equilibrium generating functional on integrating out all the fluctuation modes. Of particular interest is the requirement of analyticity at small spatial momenta noted above, that follows directly from thermal screening and the finite static correlation length. In the presence of stochastic interactions, it is known that hydrodynamic correlation functions are generically non-analytic out of equilibrium due to the presence of long-time tails [12]. The requirement that these long-distance or infrared non-analyticities must drop out at all loop orders as we specialize to hydrostatic configurations, is quite non-trivial and poses a rigid consistency check for the EFT framework. This question is particularly important because the tree-level propagators for the fluctuation degrees of freedom in these EFTs lack any intrinsic scale that could define the finite spatial correlation length. In fact, the lack of such a scale is precisely what leads to the physics of long-time tails out of equilibrium.
A crucial ingredient in the EFT framework is the discrete Kubo-Martin-Schwinger (KMS) symmetry [3]. It is motivated from its namesake KMS condition in thermal field theory [13,14], and ensures that the non-linear fluctuation-dissipation theorem [15] is satisfied by the EFT framework. While the KMS condition captures the fact that the system in question is fluctuating around the thermal state, by itself it does not guarantee the finiteness of the spatial correlation length. However, combined with the structure of hydrodynamic EFTs, we find that the KMS symmetry does, in fact, conspire in a nontrivial way to ensure the analyticity of hydrostatic correlators. The goal of this paper is to illustrate how this analyticity emerges explicitly at one-loop order in relativistic hydrodynamics for density-density two-point functions.
For simplicity, we will work in a kind of "incompressible limit", where sound modes are frozen and the theory only has longitudinal charge diffusion and transverse shear diffusion modes in the low-energy spectrum. Since sound modes are typically "higher energy" compared to diffusive modes, we expect that the simplified model faithfully captures the low-energy behaviour of the EFT. Existing notions of "incompressibility" in the literature, however, are only consistent with linearised hydrodynamics. As we include interactions, due care has to be taken while defining this limit so as to not violate the underlying structure of the EFT, in particular the KMS condition. To reaffirm the role of KMS symmetry in ensuring the analyticity of hydrostatic correlators, for the majority of this paper we will lift the KMS symmetry by untying the dissipative and stochastic fluctuation parameters in the EFT. The simplified model that we develop in this paper, which we dub diffusive hydrodynamics, will also be helpful more generally to probe stochastic signatures in hydrodynamics in the presence of interacting momentum modes. In the remainder of this introductory section, we will elaborate further on the analyticity features of interest, and also summarize our results.

Thermal correlators and summary of results
Recall that for a generic set of bosonic Hermitian operators O I , the retarded and symmetric correlators are defined respectively as The thermal expectation values are taken in the grand canonical ensemble with the Hamiltonian where Q is a conserved charge. These correlators are not independent. For example, working in Fourier space and specializing to the hydrodynamic regime ω/T 0 ≪ 1, one can derive the relation This is an incarnation of the fluctuation-dissipation theorem, which ties together statistical fluctuations in a thermal system to dissipation, and follows from the KMS periodicity condition in Euclidean time for equilibrium states.
For generic systems in thermal equilibrium, at least away from second order critical points and in the absence of any spontaneously broken symmetries, screening implies a finite static or spatial correlation length. It follows that retarded (or Euclidean) correlation functions fall off exponentially at large distances, or equivalently the Fourier space retarded correlators at zero-frequency are analytic at low spatial momentum k, i.e.
We refer to this condition as spatial analyticity. Importantly, it allows for a well-defined derivative expansion, used previously to develop a hydrostatic generating functional for hydrodynamics [10,11]. In the hydrodynamic regime, where ω/T 0 ≪ 1, |k|/T 0 ≪ 1 and we consider perturbations about an equilibrium state, the combination of eqs. (1.3) and (1.4) implies a nontrivial structure for the retarded correlator.
In this paper, we will focus on the correlations of a conserved charge n = J 0 , which is the time-component of a conserved U (1) current J µ . We can investigate the above features by looking at the symmetric and retarded correlators for n. At the tree level in the hydrodynamic EFT (i.e. in the linear response theory), they are known to take the form [16] G R,tree In terms of the chemical potential µ and conductivity σ, the coefficients appearing here are the static charge susceptibility χ = ∂n/∂µ| T and the diffusion constant D = σ/χ. The coefficient σ controls the strength of stochastic interactions. In thermal equilibrium, the KMS condition setsσ = σ, leading to these tree-level correlators obeying the fluctuation-dissipation theorem in eq. (1.3). We note that the symmetric correlator at ω = 0, which is nothing but the treelevel static propagator of n, has a pole at k = 0. This is consistent with our discussion above because symmetric correlators are not restricted by any analyticity requirement. Such poles are, of course, a generic feature of gapless degrees of freedom. The relevant point here is that despite the gapless nature of these tree-level correlations, the retarded correlator is indeed analytic at zero frequency in accordance with eq. (1.4).
While the KMS condition does not play any role at tree-level in ensuring the analyticity of the retarded correlator, it becomes crucial as we start to include loop corrections. In the hydrodynamic EFT, upon accounting for density and transverse-momentum fluctuations within diffusive hydrodynamics, and specializing to d = 3 spatial dimensions, we find that the correlators in eq. (1.5) admit one-loop corrections that behave at ω = 0 as Here Z χ,Λ are UV-sensitive renormalization factors with polynomial dependence on the momentumcutoff of the EFT. These are expected in the low-energy effective field theory due to the irrelevant nature of hydrodynamic interactions. We shall not compute these factors in detail in this paper; a thorough discussion restricted to the scalar diffusion model can be found in [9].
The remaining coefficients in eq. (1.7) are interesting physically, because they quantify the infrared and non-analytic corrections to the hydrodynamic correlation functions. These can be explicitly computed in terms of three hydrodynamic parameters: pressure p(T, µ), shear viscosity η(T, µ), and conductivity σ(T, µ), and two stochastic parametersη(T, µ) andσ(T, µ). The enthalpy density w(T, µ) and charge density n(T, µ) are defined using the usual thermodynamic relations dp = (w − µn)/T dT + ndµ. The bulk viscosity ζ(T, µ) and the associated stochastic parameterζ(T, µ) do not appear in the diffusive limit of hydrodynamics. Using the charge susceptibility χ = ∂n/∂µ| T and the diffusion constant D = σ/χ from above, and further defining the shear diffusion constant γ η = η/w and the cross-susceptibility χ ǫ = ∂w/∂µ| T , we find All the coefficients here are understood to be evaluated on the equilibrium configuration. The explicit functional form of these coefficients is not very important. The important point is that both the retarded and symmetric correlators are generically non-analytic in the static limit. However, upon imposing the KMS condition, which in this case requires settingσ = σ and η = η, we find We see that in the presence of the KMS condition, which ensures that the correlators satisfy the fluctuation-dissipation theorem, the non-analytic pieces in the retarded correlator drop out in the hydrostatic limit, reaffirming the consistency of hydrostatic equilibrium in the EFT framework.
On the other hand, the symmetric correlator in the hydrostatic limit remains non-analytic. Full analysis of the finite-ω behaviour of these correlation functions and the ensuing long-time tails will appear in a companion paper [17].
The discussion for arbitrary spatial dimensions follows in a similar manner. In d = 1, 2 spatial dimensions, 1 the corrections inside the brackets in eq. (1.7) typically behave as |k| d−2 , |k| d , . . ., while in d = 2 spatial dimensions, they behave as log(k 2 ), |k| 2 log(k 2 ), . . .. While we will not derive the explicit form of the associated coefficients in general, we will show in the course of this paper that the one-loop non-analytic corrections to the retarded function drop out in arbitrary dimensions due to the KMS symmetry, maintaining consistency with the existence of thermal equilibrium.
The paper is organized as follows. To sketch out the road-map and develop a qualitative understanding of the results, we will first ignore the momentum modes altogether in section 2, and start with a simplified EFT model of a single diffusive conserved charge as proposed in [8,9]. However, unlike the original references, we will lift the KMS symmetry from the model and compute the one-loop corrections to the symmetric and retarded correlation functions directly using the EFT generating functional. Note that we cannot deduce the two correlators using one-another in the absence of KMS symmetry, because the fluctuation-dissipation theorem in eq. (1.3) does not apply. We will then set up the theory of diffusive hydrodynamics in section 3, and repeat the computation of one-loop correlation functions in the presence of transverse momentum fluctuations. We will finish with a discussion in section 4. In appendix A, we provide further details of the derivation of the EFT for diffusive hydrodynamics and the explicit one-loop computations.

| Density correlation functions in simple diffusion
In this section, we investigate the issue of spatial analyticity of hydrostatic correlators in a simple EFT model with a single diffusive charge, which has all the relevant qualitative features of the full problem. The generalization to full hydrodynamics will be presented in the next section.

Simple stochastic diffusion
Thermal fluctuations are conventionally introduced into the hydrodynamic setup by including stochastic noise sources in the constitutive relations with short-range Gaussian correlations. One way to realize this effective theory is to add a fluctuation field for each conserved hydrodynamic variable in the theory and require that the ensuing effective action, when inserted inside a path integral, reproduces the full set of real-time correlation functions in thermal equilibrium [19][20][21] (see [16] for a review). Focusing on a single conserved charge n(µ) with the chemical potential µ, the respective classical conservation equation at one derivative order is given by where J i is the charge flux expressed in terms of µ and its derivatives. In the following we will often use the "covariant" current J µ with J 0 = n.
is the electric field associated with the background gauge field A µ coupled to J µ . The charge conductivity σ(µ) is an arbitrary function of µ, required to be non-negative by the second law of thermodynamics. The equilibrium state is given by µ = µ 0 being a constant.
In the EFT language, the chemical potential µ is related to the fundamental "chemical shift" field Λ β as µ = T 0 Λ β + A rt , where A rµ is one of the Schwinger-Keldysh background gauge fields along with A aµ . 2 To account for stochastic fluctuations, we add an auxiliary field ϕ a to partner the hydrodynamic field Λ β . The effective action S of the theory is fixed so that it reproduces all the two-point response functions at tree-level as predicted by classical linearized hydrodynamics. Introducing the effective Lagrangian density through S = d d+1 x L, we have [3,8] Here T 0 is the constant temperature of the thermal equilibrium state. We have introduced an arbitrary coefficientσ(µ) that characterizes the strength of stochastic fluctuations. This is the most general Lagrangian consistent with the Schwinger-Keldysh EFT framework for a single diffusive charge at one-derivative order. However, for this effective Lagrangian to correctly reproduce the hydrodynamic response functions, we need to impose the constraint from the fluctuation-dissipation theorem. This equates the strength of stochastic fluctuationsσ to the dissipative transport coefficient σ (see [16] for more details), Assuming the underlying microscopic theory to be CPT-invariant, in the Schwinger-Keldysh (SK) effective field theory framework of [3], this follows as a consequence of the dynamical KMS symmetry of the effective action In addition, it implies the charge conjugation properties for the coefficients n(−µ) = −n(µ), In the following, we find that the KMS condition is critical for ensuring the analyticity of retarded correlators in spatial momentum. Hence, to keep track of the non-analytic behaviour explicitly, we will perform the majority of forthcoming manipulations assumingσ and σ to be independent.

Linear expansion and interactions
To be able to use this theory in a loop expansion, we need to expand the Lagrangian (2.2) orderby-order in fluctuations about the equilibrium state Λ β = µ 0 /T 0 , ϕ a = 0; see [8] for details. It is practically easier to work with the fluctuations in the density δn = n(T 0 Λ β ) − n(µ 0 ) rather than Λ β . Ignoring the background fields at first, and ignoring certain total-derivative terms, we find the "free" Gaussian part of the effective Lagrangian to be where D = σ/χ is the diffusion constant, with χ = n ′ (µ) being the charge susceptibility, and We have dropped the subscript "0" for clarity and the coefficients are understood to be evaluated at equilibrium, i.e. µ = µ 0 . This leads to the tree-level propagators Here p = (ω, k) collectively denotes frequency and momentum. We have denoted δn by solid and ϕ a by wavy lines. Note that the δnϕ a propagator is purely retarded, ϕ a δn propagator is purely advanced, while δnδn propagator can be expressed as a sum of retarded and advanced pieces. This observation generalizes to full hydrodynamics, and will be important in the explicit loop calculations in the following discussion.
The free diffusive action (2.5) is well-known in the literature (see e.g. [16]), however, it does not account for interactions. To this end, we can expand eq. (2.2) to the next order in δn and ϕ a to obtain 3 where the two coupling constants λ andλ are defined via We can work-out the momentum-space Feynman rules for the associated vertices The same procedure can be iterated to arbitrarily high orders in fluctuations, depending on the sensitivity required. This form of the effective action, albeit for energy diffusion instead of charge diffusion, was derived recently in [9].
To be able to compute the correlations functions using the non-equilibrium generating functional [3], we also need the structure of the background field couplings in the EFT. Since, in this work, we are only interested in density correlators, we only turn on the A rt , A at components of the gauge fields. We directly borrow the results from [8], leading to where χ ′ = ∂χ/∂µ. Denoting A rt by a dotted and iA at by a dashed line, we can represent these diagrammatically as The first two of these are the usual linear couplings between operators and sources in quantum field theories, while the remaining are non-linear couplings. Frequency and momenta going away from the vertex are taken to be positive.
Let us briefly comment on the dimension counting scheme that we implement in the following discussion. For diffusive processes, it is typically argued that ∂ t ∼ ∂ i ∂ i . Following e.g. [24], we can demand that [∂ i ] = 1 and require that the couplings in spatial kinetic terms in eq. (2.5) are dimensionless. If we take [∂ t ] = 2 to ensure consistent scaling, and thus [L] = d + 2, we find [ϕ a ] = [δn] = d/2, and it follows that [λ] = [λ] = −d/2 and thus both interactions in eq. (2.7) are irrelevant. We therefore expect that loops will induce power-like sensitivity to the UV cutoff. To count loops (L), we can use the relation L = 1 + (N/2 − 1)V − E/2 for N th order vertices in diagrams with E external lines. In the case of interest below, E = 2, and so at one-loop order, V = 2/(N − 2). Thus, one-loop corrections require diagrams with 2 cubic vertices. In principle, diagrams with a single quartic vertex can also contribute at this order, however the contributions from such diagrams to the two-point functions are purely UV-cutoff dependent and only renormalize the classical hydrodynamic parameters [8,9].

One-loop corrections from density fluctuations
Defining the Schwinger-Keldysh generating functional as [3] we can compute the retarded and symmetric correlators of n via G R nn (p) = −iδ 2 W δA at (p)δA rt (−p) = χ − ωχ δn(p)ϕ a (−p) + (source coupling diagrams), Using the propagators in eq. (2.6), it is trivial to check that we reproduce the tree-level results for the two correlation functions given in eq. (1.5).
We will now compute the stochastic loop corrections to the classical tree-level propagators in eq. (1.5). This amounts to integrating out the stochastic noise field ϕ a at one-loop order using the interaction vertices defined above in eqs. (2.7) and (2.10). Let us start with the retarded function. At one-loop order, we have two diagrams contributing to the δnϕ a propagator (see [8] for details) In addition, we have two diagrams involving non-linear background source couplings (2.14b) We can compute these diagrams to find the one-loop correction to the retarded density correlator; see appendix B for calculational details. We find the one-loop correction to the retarded correlator to be where p = (ω, k) and I 1 (p) is the result of the loop momentum integral The value of this integral depends on the UV-cutoff Λ and the number of spatial dimensions d. Above, we have explicitly computed it for d = 2, 3 for reference. We see that, generically, the oneloop corrections at finite frequency exhibit IR-sensitive non-analytic dependence on ω, reflecting the physics of long-time tails, which invalidates the hydrodynamic derivative expansion at a given order [25]. However, if we restrict to the hydrostatic (ω = 0) limit, we have argued that thermal screening should render the correlator analytic. Setting ω = 0, the one-loop correction reduces to which vanishes trivially upon imposing the KMS conditionσ = σ, irrespective of the dimensionality. Hence, we see that the KMS condition is crucial to reproduce the analyticity of retarded correlators in the hydrostatic limit.
We can repeat the same procedure for the symmetric correlator. In this case, we do not have any background coupling diagrams at one-loop order. We only need to consider the seven diagrams correcting the δnδn propagator Relegating the explicit computations to the appendix, we find that It can be explicitly checked that, on imposing the KMS condition, these corrections satisfy the fluctuation-dissipation theorem (1.3) along with the retarded correlator corrections in eq. (2.15). Finally, in the ω = 0 limit, the symmetric correlator behaves as which survives the KMS limitσ = σ. This indicates that IR-singular non-analytic terms are generated in the hydrostatic limit, as a manifestation of the lack of any intrinsic mass scale in the tree-level propagators (2.6). This may appear surprising, given the earlier discussion of thermal screening. However, the symmetric correlator is not determined explicitly by the hydrostatic generating functional and does not have to meet the analyticity requirement. The leading IR-singularity, proportional to 1/k 2 , to the symmetric correlator is still given by the tree-level propagator.

| Density correlation functions in diffusive hydrodynamics
The discussion of charge diffusion above can be extended to the full theory of relativistic hydrodynamics. We will again focus on the transport of a single conserved charge, but consistently incorporate additional hydrodynamic momentum modes in the thermal bath. We will develop an effective field theory for diffusive hydrodynamics, where the sound modes have been frozen and the theory only contains diffusive charge and shear modes. This allows us to compute the hydrodynamic correlators in a slightly simpler setting, while still incorporating the non-trivial effects of momentum modes. Our goal is not to provide a comprehensive analysis, but to test the spatial analyticity of retarded charge density correlators in a system with additional gapless modes, specifically those associated with momentum fluctuations.

Stochastic diffusive hydrodynamics
The classical dynamical equations for relativistic hydrodynamics are given by charge and energymomentum conservation, along with the constitutive relations expressing J µ and T µν in terms of the hydrodynamic variables: temperature T , chemical potential µ, fluid four-velocity u µ , and their derivatives. For instance, in the Landau frame up to first order in derivatives we have 5 We have introduced the background gauge field A µ coupled to J µ , and the associated field strength F µν = 2∂ [µ A ν] , but have avoided introducing the metric source for T µν that we do not require for our purposes. Here ∆ µν = η µν + u µ u ν is the projector transverse to the fluid velocity. The thermodynamic pressure p, energy density ǫ, and charge density n are functions of T and µ, and related to each other via the thermodynamic relations dp = sdT + ndµ and ǫ + p = T s + µn for some entropy density s. The charge conductivity σ, shear viscosity η, and bulk viscosity ζ are non-negative functions of T and µ. We shall be interested in fluctuations about an equilibrium state with µ = µ 0 at a constant global temperature T = T 0 in the rest frame u µ = (1, 0).
The EFT for hydrodynamics is formulated in terms of a thermal vector β µ and "chemical shift" field Λ β , along with their stochastic noise partners X µ a and ϕ a . These are related to the aforementioned hydrodynamic fields as u µ /T = β µ and µ/T = Λ β + β µ A rµ . Similar to the diffusive EFT, we introduce the Schwinger-Keldysh double copies of background gauge fields A rµ , A aµ . With all the ingredients in place, the full non-linear effective action for relativistic hydrodynamics takes the form [1] where £ Xa denotes a Lie derivative along X µ a . We have introduced the stochastic coefficients σ,η, andζ, which are also arbitrary non-negative functions of T and µ. We can recover the previous diffusive model in eq. (2.2) by simply setting X µ a = 0 and u µ /T = δ µ t /T 0 . This is the most general effective action for relativistic hydrodynamics, in Landau frame, compatible with the Schwinger-Keldysh framework of [3], truncated at one-derivative order. We utilize the usual derivative counting scheme for relativistic hydrodynamics: in the "r" sector we take u µ , T, µ, A rµ ∼ O(∂ 0 ), while in the "a" sector we have X µ a , ϕ a ∼ O(∂ 0 ), A aµ ∼ O(∂ 1 ). The dynamical KMS symmetry in eq. (2.4) generalizes to a covariant version, where £ β denotes a Lie-derivative along β µ . In addition to the charge conjugation properties of various coefficients, the KMS symmetry relates the stochastic fluctuation coefficients to the dissipative transport coefficients via KMS :σ = σ,η = η,ζ = ζ. However, we shall not implement these constraints for now to investigate the role of KMS in ensuring the analyticity of hydrostatic retarded correlators.
Similar to our discussion in section 2, we can expand the Lagrangian order-by-order in the hydrodynamic fluctuations δβ µ = β µ − δ µ t /T 0 , δΛ β = Λ β − µ 0 /T 0 , and associated noise X µ a , ϕ a . The full theory turns out to be quite complicated in practice. To make our lives simpler, we will focus on a sub-sector, which we call diffusive hydrodynamics, where the temperature and the longitudinal velocity fluctuations have been frozen and so can be ignored. This is similar to the well-known "incompressible limit" of non-relativistic hydrodynamics and amounts to ignoring the sound modes and only focusing on the diffusive modes. We take The last two conditions are required for consistency with the KMS condition (3.5). This removes the energy conservation equation and the trace of the spatial stress tensor from the effective action. In terms of the original hydrodynamic variables, the first two conditions imply 6 where u 2 = u i u i . We have suppressed terms that are cubic or higher order in δβ i . To impose these conditions in the EFT, we include additional constraint terms in the effective action for arbitrary Lagrange multipliers Ψ r , Ψ a . For consistency with the KMS condition (3.5), these multipliers must transform as

Linear expansion and interactions
It turns out to be more convenient to work with density fluctuations δn = J t | Ar=0 − n(µ 0 ) + 2iTσu i ∂ i ϕ a and momentum fluctuations π i = T ti + 2iTη(u k ∂ k X i a + u k ∂ i X k a ). 7 The free part of the effective Lagrangian (3.4), along with the Lagrange multiplier terms, takes the form The explicit calculational details can be found in appendix A. As before, D = σ/χ is the charge diffusion constant with χ = ∂n/∂µ| T being the charge susceptibility, while γ η = η/w is the shear diffusion constant with w = ǫ + p being the enthalpy density. We still have the same tree-level propagators in the charge sector, eq. (2.6), but we also find new momentum dependent propagators given by (3.12) 6 If we were to decompose the relativistic fluid velocity according to u µ = (1, v)/ √ 1 − v 2 , these conditions will imply ∂iv i = 0 exactly, making contact with the incompressible limit of hydrodynamics. 7 These definitions of δn and π i , including the imaginary parts, are the true charge and momentum densities that are conserved in the presence of stochastic fluctuations, unlike J 0 and T 0i that are only conserved classically; see appendix A for details.

| 14
We have denoted π i with a bold and X i a with a coiled line. All the cross-propagators between the momentum and the charge sectors are zero. This is a feature of the simplified diffusive hydrodynamic theory eq. (3.7), and will not be the case in full hydrodynamics in the presence of sound modes. These propagators are made transverse to momentum by the presence of the Lagrange multipliers. The propagators for the multipliers themselves can be computed, but they do not have any physical poles and drop out for our case of interest.
At the next order in fluctuations, we have the three-point interaction Lagrangian which generalizes eq. (2.7). Various coupling constants are defined as (3.14) The vertices arising from the first two interactions in eq. (3.13) are the same as given in eq. (2.9). In addition, we get six more vertices There is another interaction vertex involving the Lagrange multiplier Ψ a . However, all the loop diagrams involving this vertex include a purely retarded/advanced loop integral and hence do not contribute to the problem at hand.
Finally, we need the coupling structure to the background fields. Turning on only the A rt , A at components of the gauge fields, we find See the appendix A for details of the derivation. Most of these are the same as the ones found in the diffusive case in eq. (2.10), but we do find four additional background interaction vertices, given by

One-loop corrections from momentum fluctuations
We can use this theory to work out the corrections to the density correlation functions arising from momentum fluctuations. The variational formulas in eq. (2.13) are still valid, but with the modified generating functional (see footnote 4) This still leads to the same tree-level two-point correlators as reported in eq. (2.6).

Retarded correlation function
For the retarded function at one-loop order in diffusive hydrodynamics, we still get contributions from the previous diagrams in eq. (2.14). In addition, we have the following three one-loop diagrams involving momentum fluctuations contributing to the δnϕ a propagator: With the incorporation of momentum fluctuations, we also have four new background coupling diagrams: We have ignored the possible diagrams which involve a purely retarded loop and hence do not contribute to the final result. The first diagram in eq. (3.19a) and the first two diagrams in eq. (3.19b) involve leading derivative non-dissipative couplings and dominate in the IR. The remaining four diagrams here and the ones in eq. (2.14) lead to subleading corrections in k 2 . The explicit computation of these diagrams can be found in the appendix B. Combining the result from eqs. (2.14) and (3.19), we find that the full one-loop correction to the retarded Hence, we see that the non-analyticities drop out of the retarded correlator in the hydrostatic limit, even after the incorporation of momentum fluctuations, provided the KMS condition is imposed.

Symmetric correlation function
We can repeat the same procedure for the symmetric correlation function as well. In addition to the 7 diagrams in eq. (2.18), we find 10 more in the presence of momentum fluctuations: The first diagrams in eqs. (3.24a) to (3.24c) contribute to leading order corrections in k 2 , while the remaining diagrams, along with those in eq. (2.18), contribute at the subleading order. We have computed these diagrams in Appendix B, which yield a complicated expression given by When the KMS symmetry is respected, the result simplifies and the remaining contributions match the relevant terms in the retarded correlation function, in accordance with the fluctuationdissipation theorem (1.3). We have computed these terms in the ω = 0 limit and verified that the answer is well-defined and still non-analytic, even in the presence of KMS. For d = 3, the results have been explicitly presented in the introduction.

| Discussion
In this paper, we explored the spatial analyticity of hydrodynamic correlation functions at loop-level, using the Schwinger-Keldysh EFT formalism that incorporates nonlinear stochastic fluctuations consistent with the KMS condition. We computed the complete one-loop correction to the retarded and symmetric two-point functions of charge density, due to diffusive density and transverse momentum fluctuations. The symmetric correlation function is found to be generically non-analytic, owing to the absence of a characteristic mass scale in the EFT. Despite this, we find that the retarded correlation function non-trivially respects the spatial analyticity requirement expected in the hydrostatic limit due to thermal screening. Notably, the analyticity behaviour crucially depends on the dynamical KMS symmetry in the EFT, that leads to a series of non-trivial cancellations of non-analytic terms in the hydrostatic limit. Our final results are summarized in section 1.1.
We performed the majority of loop calculations in this paper in generic spatial dimensions and away from the hydrostatic limit (ω = 0). This analysis can be recycled to study the phenomenon of long-time tails at finite-ω finite-k in relativistic hydrodynamics in the presence of momentum fluctuations, closely following a similar analysis for scalar diffusion model in [9]. These results will appear in a companion paper [17].
During the course of our analysis, we developed a novel EFT model for "incompressible" diffusive hydrodynamics that is consistent with the KMS condition beyond the linearised fluctuations, accounting for the presence of arbitrary non-linear interactions. This model can be understood as systematically integrating out the "higher-energy" sound modes from the hydrodynamic setup and allows us to study the effects of momentum fluctuations in a controlled setting, which is expected to have far-reaching applications beyond the context of the present work.
Although not the focus of this paper, in computing the symmetric correlation functions, we have exhibited specific non-analyticities that are computable, independent of the UV cutoff, and determined by thermodynamic parameters and hydrodynamic transport coefficients. As such, symmetric correlators are interesting quantities in themselves. One approach to computing transport coefficients within the microscopic theory utilizes Euclidean correlators on a finite size lattice. Although the physical content of Euclidean correlators most directly maps to Minkowskian retarded correlators, considering the static limit of the symmetric correlator may be helpful for extracting transport coefficients via finite size scaling.
In terms of future directions, its important to note that we have limited our considerations in this paper to generic thermal systems that are assumed to have a finite screened spatial correlation length. There are of course physical systems in which long-range correlations survive and may lead to physical non-analyticities in retarded correlators. Examples include conformal fixed points, and superfluids, which are known to require additional degrees of freedom within the hydrostatic generating functional. It would be interesting to explore how spatial non-analyticities re-emerge in such theories. The analysis can also be extended in the context of non-relativistic hydrodynamics using recently developed field theory methods [6].
Finally, we have argued that the KMS condition, associated with the existence and properties of the thermal equilibrium state, is the critical ingredient in ensuring the spatial analyticity of the retarded correlators in the hydrostatic limit. While this is, of course, consistent with the assumptions underlying the derivative expansion for the hydrostatic generating functional, it is natural to ask whether it can be imposed in a more direct and local manner within the EFT itself. Currently, the KMS condition is formulated as a non-local discrete symmetry condition to be imposed on the Schwinger-Keldysh EFT, and it would be interesting to understand if a more covariant formalism would allow for it to be imposed directly in writing down the EFT Lagrangian. In particular, we have only showed the analyticity of retarded correlators at oneloop level. It will be interesting to explore whether such a result can be derived at the full non-perturbative level in the EFT formalism of hydrodynamics.

Acknowledgments
AS would like to acknowledge the hospitality of the Centre for Particle Theory at Durham University, and the Department of Physics at McGill University, where parts of this work were carried out. This work is supported in part by NSERC, Canada. Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Colleges and Universities.

A | The effective action for diffusive hydrodynamics
In this appendix, we provide further details on the computation of the perturbative effective action for relativistic diffusive hydrodynamics presented in section 3.2, starting from the Schwinger-Keldysh effective field theory.
Let us start with the effective action for full hydrodynamics in eq. (3.4), with the constraints eq. (3.7). We need to expand the effective action order-by-order in the dynamical fields δβ µ , δΛ β , X µ a , ϕ a , and the background fields A rt , A at . Segregating space and time derivatives in the effective Lagrangian eq. (3.4), we can express it as where we have only kept terms up to cubic order in dynamical and background fields, and defined Here w = ǫ+p is the enthalpy density. The subscript "0" denotes that the corresponding quantity is evaluated at equilibrium. These are the true non-classical conserved charge and momentum densities. We can further split these into with δn and π i being the respective versions in the absence of background fields,

B | Loop calculations
In this appendix, we provide explicit details of the calculations involved in computing the loop integrals presented in the main text.
One-loop corrections to the retarded and symmetric correlation functions of charge density can be computed using the variational formulae in eq. (2.13). Starting with the retarded function, we first need to compute the one-loop diagrams correcting the δnϕ a propagator given in eqs. (2.14a) and (3.19a). Let us denote the amputated version, i.e. diagrams ignoring the external legs, by iΓ ar (p). Next, we have diagrams involving a non-linear background field coupling from eqs. (2.14b) and (3.19b), whose amputated version is denoted by iΓ ′ ar (p). The full one-loop retarded correlation function is given in terms of these as For completeness, we note that the one-loop correction to the advanced correlation function will be given by its complex conjugate, On the other hand, for the symmetric function, we need to compute the diagrams in eqs. (2.18) and (3.24). Note that the diagrams in eqs. (2.18a) and (3.24a) are the same as those for the retarded function in eqs. (2.14a) and (3.19a), except with different external legs. Hence, their amputated version is still given by iΓ ar (p). Furthermore, the diagrams in eqs. (2.18b) and (3.24b) are merely their complex conjugates, and hence their amputated version is given by −iΓ * ar (p). Therefore, the only diagrams that we need to compute for the symmetric function are given in eqs. (2.18c), (3.24c) and (3.24d). We denote their amputated version by Γ aa . In total, the one-loop correction to the symmetric function is given by G S,1-loop nn (p) = −Γ aa (p) ω 2 + D 2 k 4 − 2Tσk 2 Γ ar (p) (iω − Dk 2 )(ω 2 + D 2 k 4 ) + 2Tσk 2 Γ * ar (p) (iω + iDk 2 )(ω 2 + D 2 k 4 ) .