Matthias Krüger

Alexandre Solon

Vincent Démery

Christian M Rohwer

David S Dean

Matthias Kr Üger

Vincent D Émery

Stresses in non-equilibrium fluids: Exact formulation and coarse-grained theory

ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

I. INTRODUCTION

The physics of out-of-equilibrium fluids is fascinating in its complexity and the great variety of phenomena they exhibit. This is in part due to the many ways a system can be forced out of equilibrium, e.g., by a continuous external driving such as an applied shear, [START_REF] Dhont | An Introduction to Dynamics of Colloids[END_REF][START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] or an imposed temperature gradient, [START_REF] Ortiz De Zárate | Hydrodynamic Fluctuations in Fluids and Fluid Mixtures[END_REF] or by a (sudden) change of a control parameter as in a temperature quench producing a supercooled liquid (see, e.g., Ref. 4).

In general, describing the static properties of fluids in thermal equilibrium is already a formidable task due to the inherent many-body interactions. [START_REF] Hansen | Theory of Simple Liquids[END_REF] In this respect, density functional theory (DFT) [START_REF] Evans | [END_REF]7 has been very successful in the study of static density profiles and correlation functions. Concerning the dynamics (or non-equilibrium situations), insights have been obtained in the dilute limit from exactly solvable two-body models, 1,8-10 e.g., from the Boltzmann or the Smoluchowski equation. For multi-body dynamics, approximate treatments include mode coupling theory [START_REF] Fuchs | [END_REF]12 or dynamical density functional theory (DDFT), [13][14][15] including power functional theory. 16 Non-equilibrium molecular fluids have also been studied extensively by computer simulations. 17-20 a) Current address: Institute for Theoretical Physics, Georg-August-Universität Göttingen, 37073 Göttingen, Germany.

While the above-mentioned approaches start from the microscopic details of the fluid (at the particle level), a very different approach starts from effective field theories (such as Landau-Ginzburg theory [21][START_REF] Onuki | Phase Transition Dynamics[END_REF][START_REF] Kardar | Statistical Physics of Fields[END_REF], where microscopic details are neglected in order to study only the large-scale phenomena, notably in near-critical fluids. Non-equilibrium scenarios also yield insights and challenges in this context: Several computations have been concerned with critical Casimir forces away from equilibrium, i.e., for temperature quenches, [START_REF] Gambassi | [END_REF][25][26] moving objects, [27][28][29] or shear in confinement. 30 Non-equilibrium fluctuations arising from conservation laws have also been demonstrated to lead to Casimir forces in various setups. [31][32][33][34][35][36][37][38][39] Moreover, the pressure and stresses exerted by active systems have attracted growing attention for their unusual properties. 40,41 In all these studies, the computation of forces and stresses in out-of-equilibrium situations is non-trivial, and the form of the applicable stress tensor has been discussed (partly controversially). 26,42,43 In this manuscript, we derive several expressions for the stress tensor of a liquid. In particular, starting from the microscopic dynamics of interacting Brownian particles, we compute the exact stress tensor for any given (snapshot) density realization. Importantly, the stress tensor's average agrees with the form obtained from the Smoluchowski equation. [44][START_REF] Kreuzer | Nonequilibrium Thermodynamics and its Statistical Foundations[END_REF][START_REF] Aerov | [END_REF] Before averaging, the fluctuating form is naturally well suited 1 Stresses in non-equilibrium fluids: Exact formulation and coarse-grained theory for use in Green-Kubo fluctuation relations. 47,48 As a second step, we derive a non-equilibrium version of the stress tensor for a Landau-Ginzburg Hamiltonian and identify the various contributions, which we relate to those found in the exact microscopic version. Using the (non-equilibrium) stress tensor to compute the force on an embedded object, the results are in agreement with the existing literature. 26 We then use the insights obtained from the analysis of Brownian particles and Landau-Ginzburg theory and derive a form of the stress tensor which is valid for a general energy functional, in particular reproducing the two aforementioned cases. As an application, we use the expression for the instantaneous stress tensor to compute temporal and spatial stress correlations; this can be done exactly to leading order in the interaction potential V. We show that, after integration over time, the stress-stress correlations decay as power laws in space. These correlations have been of recent interest 49 (see also Ref. 50) and are expected to influence the flow in confinement 51 or for inhomogeneous flow velocities. 52 Finally, we use the Green-Kubo formula to deduce the viscosity of the liquid from the stress correlations.

The manuscript is organized as follows: Sec. II introduces the system under consideration and defines several observables of interest. Section III gives the exact instantaneous stress tensor for the density operator, while Sec. IV analyzes the stresses in a standard local field theory (with only zeroth and first derivative terms). Section V derives an expression for the stress tensor for a general energy functional. Section VI provides an application: We compute the correlations of the stress tensor.

II. THE SYSTEM AND OBSERVABLES

A. System

To investigate non-equilibrium liquids, we choose the well-studied and experimentally relevant model system of over-damped spherical Brownian particles. This system has the advantage that, even if driven far from equilibrium, the solvent stays equilibrated in any situation and acts as a bath at the given temperature. This way, a well-characterized out-of-equilibrium state is obtained.

Using Brownian dynamics directly implies a canonical or grand-canonical description, where the solvent acts as a bath at the given temperature. We generally consider systems for which canonical and grand-canonical descriptions are equivalent due to a very large (infinite) particle number (for example, a semi-infinite system bound by a planar surface).

The Brownian particles with positions at x µ are subject to a potential Ψ({x µ }) that includes pairwise interactions (denoted by V ) as well as an external potential (denoted U),

Ψ({x µ }) = µ<ν V (x µ -x ν ) + µ U(x µ ). (1)
Indices µ and ν run over all particles. The thermal energy scale is denoted by k B T ≡ β 1 , with Boltzmann's constant k B and the (solvent-imposed) temperature T. The bare diffusivity (in the absence of interactions) of the Brownian particles is denoted by D. Each particle thus obeys the over-damped Langevin equation

dx µ dt = D βF µ + √ 2Dξ µ , (2)
where ξ µ is a Gaussian white noise with zero mean and with correlations ξ µ,i (t)ξ ν ,j (t ) = δ µν δ ij δ(t t ), and F µ = -∇ x µ Ψ is the force acting on particle µ due to the potential Ψ. (Throughout, i and j label spatial components, while Greek indices label particles.) If Ψ = 0, each particle performs isotropic Brownian motion. Equation ( 2) is a standard theoretical model, but it neglects hydrodynamic interactions in comparison with the experimental situation of Brownian suspensions so that it is sometimes referred to as the free draining limit.

B. Observables-Mean and fluctuating-In and out of equilibrium

We summarize the important observables for this manuscript in Table I. The basic quantity is the density operator, [START_REF] Hansen | Theory of Simple Liquids[END_REF] ρ(x, t)

= µ δ(x -x µ (t)). (3)
It is the starting point for all considerations that follow. Averaged over the equilibrium distribution, one obtains the mean equilibrium density ρ, defined as

ρ(x) ≡ ρ(x) eq = µ δ(x -x µ ) eq . (4)
Here we have introduced the equilibrium average . . . eq , which, for the over-damped system, is exactly given by . . . eq = ∫ dΓ . . . e -βΨ(Γ) ∫ dΓe -βΨ(Γ) , (5) where Γ ≡ {x µ }. As noted above, for large systems the grand canonical average agrees with the canonical one given here. We introduce the density fluctuation field φ(x, t), which quantifies the deviation of the density operator from its equilibrium mean φ(x, t) = ρ(x, t)ρ(x).

Fluctuations can be characterized by their two-point correlation function

C(x, x , t, t ) = φ(x, t)φ(x , t ) , (7)
Density operator: ρ(x, t) = µ δ(x x µ (t)) ρ(x) = ρ(x) eq
Mean density in equilibrium φ(x, t)

Fluctuation of density about its equilibrium value, φ(x, t) = ρ(x, t)ρ(x) φ(x, t)φ(x , t ) Time-dependent correlations of density fluctuations ρ(x, t)

Mean density ρ [START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] (x, x ) Averaged two-body density:

ρ (2) (x, x ) = µ ν δ(x -x µ )δ(x -x ν ) g(r)
Pair correlation function in bulk:

g(r) = 1 N ρ ν µ δ(r -x ν + x µ )
where we introduced the average . . . over noise realizations given a (possibly non-equilibrium) ensemble of initial conditions. The correlation function in Eq. ( 7) is thus well-defined in or out of equilibrium. In stationary state, C is a function of t t only [START_REF] Risken | The Fokker-Planck Equation[END_REF] and depends only on the relative coordinate x x in homogeneous systems. Its spatial Fourier transform, C is the intermediate scattering function. [START_REF] Dhont | An Introduction to Dynamics of Colloids[END_REF] Another important quantity (related to C) is the two-body density, defined by 5

ρ (2) (x, x , t) = ν µ δ(x -x µ (t))δ(x -x ν (t)) . ( 8
)
For bulk systems, ρ [START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] depends only on one coordinate r (see Ref. 5) and can be expressed via the pair correlation function g,

g(r) = 1 ρ 2 ρ (2) (x + r, x) = 1 N ρ µ ν δ(r + x µ -x ν ) . (9)
Again, the above average for g can be evaluated in or out of equilibrium. One example for a non-equilibrium pair correlation function is found in systems under shear, where g is distorted compared with that of the reference equilibrium case. 1

III. EXACT MICROSCOPIC STRESS TENSOR

In this section, we derive and discuss the exact microscopic stress tensor (including its fluctuations) for the system of Brownian particles.

A. Microscopic theory

The stress tensor is related to forces in the system. These forces can be read off directly on the equation of motion, which for the density operator ρ is given by 54

∂ ρ ∂t (x, t) = ∇ • Dρ(x, t)∇ δ βE δ ρ(x, t) + 2D ρ(x, t)η(x, t) . ( 10
)
Equation (10), interpreted with the Itō convention, is an exact reformulation of Eq. ( 2). The term η is a vectorial Gaussian white noise field with zero mean and correlations

η i (x, t)η j (x , t ) = δ ij δ(x -x )δ(t -t ). ( 11
)
E is the (free) energy functional, which contains an ideal gas part, a contribution from interactions via the inter-particle potential V, and the external potential U,

E[ρ(x)] = k B T dx ρ(x) ln(ρ(x)) + 1 2 dxdx ρ(x)V (x -x )ρ(x ) + dx ρ(x)U(x). (12)
The reader should note that E is not the free energy functional of DFT. [START_REF] Evans | [END_REF] Indeed, one important difference of this work with respect to DDFT [13][14][15] is the presence of the noise in Eq. (10).

The stress tensor for DDFT was discussed in Ref. 46. Equation ( 10) may now be rewritten for identification of the stress tensor σ: The divergence of the stress tensor appears directly, [START_REF] Kreuzer | Nonequilibrium Thermodynamics and its Statistical Foundations[END_REF][START_REF] Aerov | [END_REF] ∂ ρ ∂t

(x, t) = -βD∇ • [∇ • σ(x, t) -ρ(x, t)∇U(x)] + ∇ • 2D ρ(x, t)η(x, t) . ( 13
)
Equation ( 13) is a force balance between external and interparticle forces, where the latter are expressed via the stress tensor. The divergence of the stress tensor is thus identified by comparing Eqs. ( 13) and (10). Using Eq. ( 12), we obtain

∇ • σ(x) = -k B T ∇ρ(x) -ρ(x)∇ dx V (x -x )ρ(x ). ( 14
)
The divergence of the stress tensor at position x thus has a local entropic or osmotic contribution involving the density at x, and an interaction term which involves the potential V.

We emphasize that Eq. ( 14) gives the instantaneous stress tensor, which is valid for any given (snapshot) configuration of particles.

The noise-averaged form of Eq. ( 14) can be rewritten using the two-body density of Eq. ( 8), which can also be expressed via

ρ(x)ρ(x ) = ρ(x) δ(x -x ) + ρ (2) (x, x ). ( 15
)
Noting that the first term on the rhs of Eq. ( 15) does not contribute to Eq. ( 14) (reflecting the fact that a particle cannot exert a force on itself), we obtain the familiar form for the divergence of σ, 45

∇ • σ(x) = -k B T ∇ ρ(x) -dx [∇V (x -x )]ρ (2) (x, x ). ( 16
)
This expression may, for example, be found starting from the Smoluchowski equation. [START_REF] Aerov | [END_REF] The expression for the stress tensor itself, both instantaneous and averaged, can be obtained from Eqs. ( 14) and ( 16), respectively. Indeed, we show in Appendix A that for a spherical potential V (r) = V (r) with r = |r|, the following expression of σ leads to the correct force balance,

σ(x) = -k B T ρ(x)I + 1 2 1 0 dλ dr × rr r V (r)ρ(x + (1 -λ)r)ρ(x -λr). (17)
Averaging Eq. ( 17), we find for the mean stress tensor

σ(x) = -k B T ρ(x) I + 1 2 1 0 dλ dr × rr r V (r)ρ (2) (x + (1 -λ)r, x -λr) , ( 18
)
which is the celebrated Irving-Kirkwood formula 44 for the stress tensor; Eq. ( 17) extends it to individual microscopic configurations. Note that adding a divergence-free term to the stress tensor does not change the force balance in Eq. ( 13) so that different expressions of the stress tensor are acceptable. However, the expression of Eq. ( 18) can be argued to possess the most physical symmetries. 55

B. Stress tensor for the field

For further use, we also give the form of the stress tensor in terms of the fluctuating field φ, such that ρ = ρ + φ. From Eq. ( 17), this directly gives

σ(x) = -k B T ρ(x) + φ(x) I + 1 2 1 0 dλ dr × rr r V (r) ρ + φ (x + (1 -λ)r) ρ + φ (x -λr). (19)
If the correlations of the field φ are known (e.g., assuming Gaussian fluctuations 56 ), Eq. ( 19) can then be used to compute the correlations of the stress tensor and the viscosity via Green-Kubo relations. Equation ( 19) can be simplified further to a form which displays clearly the off-diagonal components, and which will be useful for computing shear viscosity. Up to a divergencefree term, one gets for bulk systems, where ρ is not a function of x, that [this may be directly inferred from Eq. ( 14)]

σ(x) = -k B T φ(x) + ρ dx V (x -x )φ(x ) I + 1 2 1 0 dλ dr rr r V (r)φ(x + (1 -λ)r)φ(x -λr). (20)
The above equation shows that only the terms quadratic in φ contribute to the off-diagonal components of the stress tensor in bulk (and thus, as we will see later, to the viscosity). While Eq. ( 20) is useful for computation of off-diagonal elements, care should be taken for the diagonal term in this representation, because the integral over the potential may become arbitrarily large [Eq. ( 20) omits diagonal terms due to ρ, which are necessary to regularize the pair correlation for large V ].

C. The different terms in the force balance

The terms in Eq. ( 13) are interpreted physically as force densities. There is the force (density) acting on the external potential U (for example, the force acting on a wall which bounds the fluid),

f (U) (x) = ρ(x)∇U(x). (21)
In equilibrium, there is no net particle current in the system so that the external force balances on average the divergence of σ f (U) (x) eq = ∇ • σ(x) eq . ( 22) Equation ( 22) reflects the well-known fact that, in equilibrium, the stress tensor is directly related to the force acting on walls or embedded objects. Out of equilibrium, the mismatch of f (U ) and ∇ • σ gives rise to particle currents.

Taking the time derivative of the density operator in Eq. ( 3), one obtains

∂ ρ ∂t (x, t) = -∇ •        µ v µ (t)δ(x -x µ (t))        , ( 23
)
where v µ = dx µ /dt. The instantaneous current is thus identified to be j

(x, t) = µ v µ (t)δ(x -x µ (t)). ( 24
)
Comparing Eqs. ( 23) and ( 13), we obtain for the mean current

k B T D j ≡ f (j) = ∇ • σ -f (U) . (25)
Equation (25) shows that ∇ • σ and f (U ) balance the frictional force f (j) of that current. 1 Equation ( 25) has yet another important consequence: Out of equilibrium, the stress tensor σ may not generally be used to compute forces on walls or embedded objects because one needs to account for the contribution of the currents.

IV. STRESS TENSOR IN LANDAU-GINZBURG THEORY

In this section, we consider the stress tensor for field theories obeying the Landau-Ginzburg Hamiltonian. In this phenomenological description, no fundamental relations for forces exist, and, indeed, the form of the stress tensor in this context has been a subject of recent debate. 26,42,43 The form for the stress derived in this section is in close analogy to the microscopic stress tensor of Eq. ( 17), and the found equation of motion is analogous to Eq. ( 25). This derivation, in direct comparison to the exact one, thus sheds light on the aforementioned debate and demonstrates how non-equilibrium forces in Landau-Ginzburg theory can be found unambiguously. It also demonstrates the main differences between computations of forces in or out of equilibrium.

Phenomenological field theories are particularly well suited to investigate large-scale generic phenomena. One may, for example, study the universal aspects of systems near critical points [21][START_REF] Onuki | Phase Transition Dynamics[END_REF][START_REF] Kardar | Statistical Physics of Fields[END_REF] or the properties of long-ranged correlations which are present due to out-of-equilibrium initial conditions as, e.g., in Ref. 38. In such scenarios, one seeks expressions independent of the microscopic details (such as the interaction potential V ). We thus investigate in this section the possibility of expressing the stress tensor directly at the coarse-grained level of the field theory, based purely on the Landau-Ginzburg Hamiltonian. Let us consider the Hamiltonian for a scalar field Φ in d dimensions,

H[Φ] = dx κ 2 (∇Φ) 2 + U(x) ≡ dx H(x). (26)
Although higher orders in ∇Φ can be included based on symmetry arguments, [START_REF] Kardar | Statistical Physics of Fields[END_REF] we do not consider this case here. U can be a general polynomial of Φ, but the simplest example of the above field theory is the Gaussian case where

U(x) = m(x)Φ(x) 2 /2.
Here, m(x) can be a function of position so that it may include contributions from external potentials. For bulk, with m constant in space, the correlation length is then set by √ κ/m. [START_REF] Kardar | Statistical Physics of Fields[END_REF] As mentioned, at the level of effective field theories such as Eq. ( 26), the definition of mechanical quantities is not obvious. Indeed, one needs to define the nature of the field Φ and specify whether it corresponds to a matter field (for instance, a particle density, a spin density, or local charge density) or a potential field (for instance, the local electrostatic potential or chemical potential). If the field Φ is in an arbitrary non-equilibrium configuration, in the mechanical sense, then there are locally unbalanced body forces in the system. We proceed by applying a fictitious external field so that the system is in local mechanical equilibrium. The force exerted by this fictitious field on a volume therefore cancels out exactly the local body forces generated by the internal interactions in the system. Concretely, we apply an external field h which shifts the overall energy (Hamiltonian) of Φ to

H h [Φ] = H[Φ] + dx h(x)Φ(x). ( 27
)
Next, for a given configuration of Φ, the external field is chosen to ensure local mechanical equilibrium of the given configuration when the field is applied, i.e., δH h δΦ(x) = 0. The required field h is thus

h(x) = - δH[Φ] δΦ(x) , (28)
from which we can extract the force density f h acting on the sub-volume V due to the imposed field, which follows from its fundamental definition from the Hamiltonian H h (see e.g., Ref. 26),

f h (x) = -Φ(x)∇h(x) = Φ(x)∇ δH[Φ] δΦ(x) . ( 29
)
If Φ is a density field, the force density may be viewed intuitively as the product of the density and the gradient of the chemical potential associated with

H[Φ].
Because the force due to h balances the force density that is not due to h, we obtain for the (body) force density f (j) when h = 0:

f (j) (x) = -Φ(x)∇ δH[Φ] δΦ(x) . ( 30
)
We added the superscript j in order to emphasize the similar nature of the forces in Eqs. ( 30) and ( 25), which will become more apparent below. We note that Eq. ( 30) is exactly the force density in Eq. (10), when one makes the identification Φ ≡ ρ and H = E. It is important to note that if the field theory is written down in terms of the fluctuations of the field Φ about its average value Φ, as Φ = Φ + φ, and the Hamiltonian for the fluctuations is H f [φ], the local body force is given by

f (j) (x) = -Φ(x) + φ(x) ∇ δH f δφ(x) . ( 31
)
In the limit of small fluctuations the above can be used to derive an approximate model B dynamics for interacting particle systems, [57][58][59][60] where only linear terms in φ are kept, yielding

f (j) (x) = -Φ(x)∇        dx δ 2 H f δφ(x)δφ(x ) φ=0 φ(x )        . ( 32
)
Next, we identify the force on the external potential, as in Eq. ( 21). Here, the force density f (U ) is identified by the change of H under a small displacement of the potential or, equivalently, a small displacement of the object giving rise to the potential U. 26 Defining the vector X from the origin to a randomly chosen point on the object, we have

f (U) = -∇ X H. ( 33
)
Manipulating the partial and functional derivatives as detailed in Appendix B, we obtain the relation between the force density in Eqs. ( 30) and (33),

f (j) = ∇ • T -f (U) , ( 34
)
where we have introduced the stress tensor T,

T ij (x) = δ ij H(x) -Φ(x) δH δΦ(x) - ∂H ∂∇ j Φ(x) ∇ i Φ(x). (35)
Equation ( 35) is the main result of this section. (In this section, we use T instead of σ, following the usual field theory notation.) The structure of Eq. ( 34) is identical to that in Eq. ( 25) so that the field theory stress tensor from Eq. ( 35) is on the same footing as the microscopic one. As mentioned below Eq. ( 25), the computation of forces on external objects, f (U ) , is non-trivial out of equilibrium because of the term involving the current in Eq. ( 34). This will be relevant, for example, when computing Casimir forces between walls or objects. As a crosscheck, we show in Appendix B that Eq. ( 34) gives the same value for f (U ) as the expression previously derived in Eq. ( 18) of Ref. 26. A particular case worth mentioning is that of no-flux boundary conditions, where the surface normal component of f (j) is forced to vanish. For the geometry of two parallel plates, the force (or pressure) acting on the plates can then be computed by evaluating the stress tensor in Eq. ( 35) at the given surface. This was used for computation of non-equilibrium Casimir forces in Ref. 38 and tested quantitatively with simulations of interacting Brownian particles in Ref. 39.

In global equilibrium, one can show 26 that f (j) vanishes on average [see Eq. (B6)], as was the case in the microscopic theory in Eq. (25). Our formulation (and the form of T) then agrees with the commonly accepted equilibrium definition of the stress tensor T eq , i.e., ∇ • T eq = ∇ • T eq eq = f (U) eq (36) and T eq expressed as [using Eq. (B6)] 61

T eq ij = δ ij H - ∂H ∂∇ j Φ ∇ i Φ. (37)

V. GENERAL FORMULA FOR THE STRESS TENSOR

Inspired by the similarities of the formulations in Secs. III and IV, we now proceed to derive a form for the stress tensor for a more general case, which encompasses energy functionals such as Eqs. ( 12) and (26). To this end, we introduce a Hamiltonian for a field Φ [which, for the case of Eq. ( 12), takes the role of the density ρ],

H[Φ] = 1 2 Φ(x)V (x -x )Φ(x )dxdx + U(Φ(x), x)dx.
(38) Here, the first term incorporates non-local interactions, while the second part describes local interactions, as well as external potentials. U(Φ, x) can be any function of Φ(x) and x for the following derivation to be valid. To make connections to the previous sections, we note that the local part for Brownian particles of [Eq. ( 12)] is recovered with

U(Φ, x) = k B T Φ log(Φ) + UΦ, ( 39
)
while the non-local part of Eq. ( 38) already takes the form of Eq. ( 12). On the other hand, the non-local part of the Landau-Ginzburg case [Eq. ( 26)] is recovered for

V (x) = -κ∇ 2 δ(x), (40)
which, used in Eq. ( 38), yields κ 2 (∇Φ) [START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] . (Note that by "nonlocal" we mean any contributions giving rise to correlations of the field Φ with a non-zero correlation length in equilibrium.) This shows that the cases of Secs. III and IV are contained in Eq. (38).

In order to derive the stress tensor from Eq. ( 38), we again rely on the body force density of Eq. ( 30), which is equally valid for the Hamiltonian of Eq. ( 38). This yields

f (j) (x) = -Φ(x)∇ δH[Φ] δΦ(x) , = -Φ(x)∇ dx V (x -x )Φ(x ) -Φ(x)∇ ∂U ∂Φ . ( 41
)
We now use Eq. (B1) for U,

∇ i U(Φ(x), x) = (∇ i Φ) ∂U ∂Φ - ∂U ∂X i = (∇ i Φ) ∂U ∂Φ + f (U) i , (42)
where, as before, -∂U ∂X i is the force density in direction i acting on the object described by the coordinate X. We now find the divergence of the stress tensor,

∇ • σ(x) = f (j) (x) + f (U) (x) (43) = -Φ(x)∇ dx V (x -x )Φ(x ) +∇ U -Φ(x) ∂U ∂Φ(x) , (44)
where we could equivalently have used the symbol T for the stress tensor, as in Sec. IV. Equation (44) shows that the local part of H yields a diagonal contribution to the stress tensor, while its non-local part further yields off-diagonal components. This also demonstrates that stresses are only propagated in the presence of a correlation length in the system. In order to invert the divergence of the non-local part, we use the steps of Appendix A. We arrive at

σ ij (x) = δ ij U(Φ(x), x) -Φ(x) ∂U ∂Φ(x) + 1 2 dr r i r j r V (r) 1 0 dλΦ(x -λr)Φ(x + [1 -λ]r). (45)
This is the main result of this section. To repeat, this form of the stress tensor can be used for energy functionals which contain a quadratic non-local part, and an arbitrary local part, as given in Eq. (38). As mentioned before, this result reproduces, in particular, the stress tensor for the exact description of Brownian particles, Eq. ( 12), as well as for the phenomenological Landau-Ginzburg theory. The two theories are thus solidly linked, and the claimed validity of Eq. ( 35) is underpinned. This is hopefully enlightening regarding the debate of stresses in phenomenological theories. In order to demonstrate the above-mentioned recovery of Eqs. ( 17) and ( 35), we start with the Brownian case. Using Eq. ( 39) in Eq. ( 45) yields

U(Φ(x), x) -Φ(x) ∂U ∂Φ(x) = -k B T Φ(x), (46)
which is the correct result-see Eq. ( 17).

To reproduce Eq. ( 35), we use V (x) = κ∇ 2 δ(x) in Eq. (45). After some computation steps, we find

1 2 dr r i r j r V (r) 1 0 dλΦ(x -λr)Φ(x + [1 -λ]r) = κ δ ij 6 2Φ∇ 2 Φ -(∇Φ) 2 + 2 3 Φ∇ i ∇ j Φ - 1 3 ∇ i Φ∇ j Φ . ( 47
)
On the other hand, evaluating Eq. ( 35) for H = κ 2 (∇Φ) [START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] , one obtains

κ δ ij 1 2 [∇Φ(x)] 2 + Φ(x)∇ 2 Φ(x) -∇ i Φ(x)∇ j Φ(x) . ( 48
)
While Eqs. ( 47) and ( 48) are not equal, their difference

2κ 3 ∇ i (Φ∇ j Φ) -δ ij ∇ k (Φ∇ k Φ) (49)
is divergence-free. This demonstrates that Eq. ( 45) indeed reproduces Eqs. ( 17) and (35). The observed difference in divergence-free terms is acceptable and has no physical consequences as regards stresses and forces. Indeed, these differences are introduced when inverting the divergence of Eq. ( 44), which is a not unique process.

VI. APPLICATION: STRESS CORRELATIONS OF BROWNIAN SUSPENSIONS

In this section, we make use of the expressions derived for the microscopic stress tensor Eq. ( 19) to compute the two-point correlations of the stress at different positions and times for a bulk equilibrium system. We provide the limit of high temperatures, which corresponds to the leading-order term in the external potential V. We also compare with the results obtained from the Gaussian field theory, using Eq. ( 35) to define the stress tensor.

Note that this two-point correlation function can be used, via Green-Kubo relations, [START_REF] Kubo | Statistical Physics II: Nonequilibrium Statistical Mechanics[END_REF] to compute the viscosity (in linear response to a flow field). Because of this, and because it incorporates the fluctuations around equilibrium, the correlation function is a quantity whose computation requires a form for the stress tensor which is valid out of equilibrium. [For example, the form of Eq. ( 37) should not be used.]

A. Diagonal components-Pressure fluctuations

We start with the diagonal part of the stress tensor in Eq. ( 19), denoting δσ = σ σ eq . The leading term at small V (or high T ) results from the ideal contribution in Eq. (19). It reads (without summing repeated indices) lim βV →0 δσ ii (x, t)δσ kk (0, 0) eq = (k B T ) 2 φ(x, t)φ(0, 0) eq . ( 50) In order to determine the leading order in V, the correlation function in Eq. ( 50) should be evaluated for V = 0, which corresponds to the ideal gas. The correlations for the ideal gas are denoted • id and are computed in Appendix C. The two-point correlation reads

Φ(x, t)Φ(0, 0) id = ρ (4πDt) 3/2 e -|x| 2 4Dt . ( 51
)
We obtain lim βV →0

δσ ii (x, t)δσ kk (0, 0

) eq = (k B T ) 2 ρ (4πDt) 3/2 e -|x| 2 4Dt . ( 52
)
We observe from Eq. ( 52) that the temporal correlations are due to a diffusion process. The long-ranged character of the correlations in space becomes more apparent after integrating Eq. ( 52) over time

∞ 0 dt lim βV →0 δσ ii (x, t)δσ kk (0, 0) eq = (k B T ) 2 ρ 4πD|x| . (53)
This quantity decays slowly, as the inverse of the distance.

Because the time integral of such equilibrium correlations appears typically in linear response formulae, we expect the long-ranged form of Eq. ( 53) to be relevant for non-equilibrium perturbations.

As a final note, the equal-time fluctuations of pressure from Eq. ( 52) are

dx lim βV →0 δσ ii (x, 0)δσ kk (0, 0) eq = (k B T ) 2 ρ. ( 54
)
These are the pressure fluctuations of an ideal gas, which have a partly controversial history, starting with Gibbs. [START_REF] Gibbs | Elementary Principles in Statistical Mechanics[END_REF] We note that the stress tensor of Landau-Ginzburg theory, Eq. ( 35), together with a Gaussian Hamiltonian, Eq. ( 26), yields a qualitatively different result for spatial pressure correlations. The reason is that Eq. ( 35) yields a stress tensor which is purely quadratic in the fluctuating field Φ, resulting from the symmetry assumptions underlying the field theoretic Hamiltonian in Eq. ( 26). [START_REF] Kardar | Statistical Physics of Fields[END_REF] In contrast to that, it is the term ∼ρ ln ρ in Eq. ( 12) which yields the correct ideal gas contribution.

B. Diagonal-off-diagonal correlations

As visible in Eq. ( 20), the off-diagonal component of the stress involves the potential V. The leading term of the correlation between diagonal and off-diagonal components of σ is thus linear in V. This term reads (for i j)

lim βV →0 δσ ij (x, t)δσ kk (0, 0) eq = k B T 1 0 dλ dr × r i r j r V (r) φ(x + (1 -λ)r), t)φ(x -λr, t)φ(0, 0) id , (55)
where, again, the correlation of φ is evaluated for the ideal gas. Since the density of the ideal gas has Poissonian statistics, [START_REF] Velenich | [END_REF] the three-point correlator in Eq. ( 55) is per se non-zero. Using Eq. (C7), we find the correlation in terms of Fourier transforms,

φ(k, t) φ(k , t) φ(q, 0) id = (2π) 3 ρδ(k + k + q)e -Dq 2 t , (56)
where, in this paper, Fourier transforms are defined by

f (k) = dxf (x)e -ik•x . ( 57
)
While there is indeed a non-zero three-point correlation in Eq. ( 56), its contribution to the correlation of stresses in Eq. ( 55) is zero (Appendix C). The reason is that Eq. ( 56) results from a single Brownian particle and its self-correlations (as is apparent from the linearity in ρ), and this particle cannot exert a force upon itself. Some relations useful to derive this statement mathematically are given in Appendix C. We thus note that there is no correlation between diagonal and off-diagonal stress components at the given order in V.

C. Off-diagonal components-Shear stress fluctuations

Here we start by computing the exact leading term in V in the microscopic theory. We then also derive the shear stress fluctuations in the Gaussian Landau-Ginzburg theory.

Exact leading order in V

To leading order, the off-diagonal component of the stress correlator is proportional to V 2 . For ease of notation, we abbreviate,

lim βV →0 δσ ij (x, t)δσ kl (0, 0) = Σ ijkl (x, t). ( 58
)
For i j and k l, this expression reads

Σ ijkl (x, t) = 1 4 drdr 1 0 dλdλ r i r j r V (r) r k r l r V (r ) × φ(x + r -λr, t)φ(x -λr, t) ×φ(r -λ r , 0)φ(-λ r , 0) id . (59)
In Fourier space, the correlation reads (again, see Appendix C)

φ(k, t) φ(k , t) φ(q, 0) φ(q , 0) id = ρ(2π) 3 δ(k + k + q + q )C(k + k , t) + (2π) 6 ρ 2 δ(k + k )δ(q + q ) + δ(k + q)δ(k + q )C(k, t)C(k , t) + δ(k + q )δ(k + q)C(k, t)C(k , t) , (60)
with C(k, t) = e -Dk 2 t . The four-point correlation function in Eq. ( 60) contains a term linear in ρ, which, as in Eq. ( 56), shows the non-Gaussian behavior of ideal particles. Again, this term does not contribute to the stress correlator. The last three terms in Eq. ( 60) are equivalent to a Gaussian decoupling, where all pairs of functions φ appear as in Wick's theorem. The first of these three is time independent and does not contribute to Eq. ( 59) by symmetry. We then find for i j and k l (61) where for i j

Σ ijkl (x, t) = dk (2π) 3 e ik•x 1 0 dλdλ ρ 2 2 dp (2π) 3 × Ãij (λk -p) Ãkl (-λ k + p)e -D[p 2 +(k-p) 2 ]t ,
Ãij (k) = k i k j k Ṽ (k). ( 62
)
We continue by analyzing Eq. ( 61) in real space, which is found easily from the Gaussian decoupling of the four-point function in Eq. ( 59). We thus have, using Eq. ( 51) (again i j),

Σ ijij (x, t) = 1 2 drdr 1 0 dλdλ r i r j r V (r) r i r j r V (r ) × ρ2 (4πDt) 3 e -|x+r-λr-r +λ r | 2 -|x-λr+λ r | 2 4Dt . ( 63
)
We now address the large distance behavior of the stress correlations. For a short-ranged potential, the values of |r| and |r | that contribute to the above integral are constrained to that range. For values of |x| much larger than the interaction range, we may expand the exponential in Eq. ( 63) for |r| |x| and |r | |x|. Because of the symmetry of the integrands in Eq. ( 63), we seek the leading terms of order r i r j r i r j of the exponential in Eq. (63). The integrals over λ and λ can then be performed to obtain

Σ ijij (x, t) = G[V ] ρ2 10D 2 t 2 -Dt(x 2 i + x 2 j ) + x 2 i x 2 j 9216(Dt) 7 π 3 e -|x| 2 2Dt , ( 64
)
where we have defined the functional

G[V ] = 8π 2 225 ∞ 0 drr 5 V (r) 2 . ( 65
)
Here, too, we perform the time integral, which is the important quantity in Green-Kubo relations, and obtain

∞ 0 dtΣ ijij (x, t) = G[V ] ρ2 5 + 40 x 2 i x 2 j |x | 4 -4 x 2 i +x 2 j |x | 2 48Dπ 3 |x| 8 . ( 66
)
Shear stress fluctuations in Eq. ( 66) and pressure fluctuations in Eq. ( 53) thus decay as power laws in space. These observations should have implications for the rheology of suspensions in confined systems. 30,50,51 The long-ranged property of the shear stress correlations has been noted in Ref. 49 via a Zwanzig-Mori approach. In particular, a Fourier inversion of Eq. ( 6) from Ref. 49 with s = 0 yields the analogue of ∫ ∞ 0 dtΣ ijij (x, t), giving contributions ∼ |x| 3 . These show a different power than those in Eq. ( 66), ∼ |x| 8 , which could be due to the absence of momentum conservation in our formulation. This must be investigated in the future. While the framework of DDFT has been extended beyond Brownian dynamics (e.g., including inertia 15,65 ), the approach of this work may be further developed in analogy to "model H." 21

From Landau-Ginzburg theory

Next we consider the off-diagonal stress correlations from Landau-Ginzburg theory, using the stress tensor T in Eq. (35). This, together with a Gaussian decoupling of the four-point correlation (again, i j), gives

T ij (x, t)T ij (0, 0) ≈ κ [∇ i ∇ j C(x, t)] 2 + [∇ i ∇ i C(x, t)][∇ j ∇ j C(x, t)] . ( 67
)
We are interested in the case where |x| is much larger than the correlation length. In this limit, the correlation function reads in model B dynamics [START_REF] Kardar | Statistical Physics of Fields[END_REF]56 (see Appendix D for the leading order in correlation length)

C(x, t) = k B T m 1 (4π Dt) 3/2 e -|x| 2 4 Dt , ( 68
)
where m is the coefficient appearing in the Gaussian Hamiltonian in Eq. ( 26) of the field theory. Explicitly, m may also be expressed in terms of the small wave-vector limit of the direct correlation function c (2)5,56 (see Appendix D),

m = k B T 1 ρ -c (2) 0 , (69)
which may also be related to the isothermal compressibility. [START_REF] Hansen | Theory of Simple Liquids[END_REF] We have also introduced D ≡ D/S 0 , with S 0 = (1ρc (2) 0 ) -1 = k B T / ρm being the small wave-vector limit of the structure factor; D is an approximation for the diffusion coefficient, [START_REF] Dhont | An Introduction to Dynamics of Colloids[END_REF] which emerges when taking the small wave-vector limit of Eq. (D4). Equation (67) yields

T ij (x, t)T ij (0, 0) = κ 2 (k B T ) 2 m 2 2 D2 t 2 -Dt(x 2 i + x 2 j ) + x 2 i x 2 j 512( Dt) 7 π 3 e -|x| 2 2Dt . ( 70
)
Apart from prefactors, this expression is qualitatively equal to Eq. ( 63): It can be shown that the dependencies on x and t of Eqs. ( 70) and ( 63) differ only by terms which arise from divergence-free parts of the stress tensor.

D. Viscosity

We finish with the evaluation of the shear viscosity η to leading order in the interaction potential. Using the Green-Kubo relation, [START_REF] Kubo | Statistical Physics II: Nonequilibrium Statistical Mechanics[END_REF] the viscosity is given in terms of the stress correlator (i j),

η = β V ∞ 0 dt dxdx σ ij (x, t)σ ij (x , 0) eq , ( 71
)
where V is the volume of the system. Performing one spatial integral, one obtains

η = β dx ∞ 0 dt σ ij (x, t)σ ij (0, 0) eq . (72)
Importantly, one should not integrate the result of Eq. ( 66) over x to obtain the result of Eq. (72), because Eq. ( 66) is only valid for large x. Instead, the integration over x is straightforward from Eq. ( 61), since it selects the mode k = 0 in the integral, leading to

η = β ρ 2 4D dp (2π) 3 p 2 i p 2 j p 4 Ṽ (p) 2 = β ρ 2 120π 2 D ∞ 0 dpp 2 Ṽ (p) 2 .
(73) This is the exact result for the viscosity of Brownian suspensions to leading order in the potential V, which is naturally quadratic. This expression is also in agreement with perturbation computations, where the perturbed pair distribution is found under shear. While such results have been obtained for hard spheres, 9 we have not been able to find the result of Eq. (73) in the literature for a direct check.

As pointed out, the stress tensor is not uniquely defined, but only up to divergence-free terms. The viscosity in Eq. ( 73) is, however, a measurable material property, which should be computable unambiguously. We therefore demonstrate that divergence-free terms do not contribute to Eq. ( 72). First we note that for any vector (or tensor) field A(x) which decays quicker than 1/|x| 3 for large |x|, we have

dxA i (x) = dx(∇ j x i )A j (x) = -dxx i ∇ j A j (x). (74)
We have used that ∇ j x i = δ ij . In the last step, partial integration was performed using the decay properties of A(x) at large |x| so that boundary terms vanish. This yields the relation

dxA(x) = -dxx∇ • A. (75)
Proceeding by applying this to the auto-correlator of stress in Eq. (72), we denote A(x) = ∫ ∞ 0 dt δσ(x, t)δσ ij (0, 0) eq . If this expression decays quickly enough, any divergence-free parts do not contribute to the viscosity in Eq. (72). Using the form of σ from Eq. ( 20) with i j, we see that, to the considered order, the stress correlation decays quickly enough, and that therefore the viscosity in Eq. ( 73) is unambiguous.

VII. SUMMARY

The exact stochastic equation for the evolution of the density operator [START_REF] Dean | [END_REF] allows one to identify the (stochastic) form of the stress tensor directly. The mean of this quantity naturally agrees with the one found from the Smoluchowski equation. It is interesting to note that non-equilibrium states with a finite particle current give rise to an additional force contribution so that the stress tensor has to be applied with care when computing forces on boundaries or objects. The non-equilibrium body force in a general field theory, where the field Φ is associated with a density operator, can be derived using force balance arguments. The result of this agrees with that obtained from the evolution of the density operator for the field theory given by Eq. ( 12) with the identification Φ = ρ. This allows a systematic derivation of the field theoretic stress tensor, valid out of equilibrium, and the associated procedure for computing non-equilibrium forces on immersed objects/potentials.

We further computed the spatial and temporal correlations of the stress tensor exactly to leading order in interaction potential V. These correlations (or more generally forces) are more long ranged than the equilibrium correlation length, a fact which we attribute to the conservation of particles in the system considered here. In contrast to Ref. 49, our system conserves particle density, but not momentum. Offdiagonal stress correlations from the field theory were shown to agree qualitatively with those from the microscopic stress tensor.

Future work can thus involve computing these correlations for the case of particle density and momentum conservation. The impact of such correlations on flow in confinement poses further interesting questions. 30,51 We have added a term in the Hamiltonian, which depends explicitly on a position X, the position of the object giving rise to the external potential U. Another explicit form is

δH δΦ(x) = ∂H ∂Φ -∇ j ∂H ∂∇ j Φ . (B2)
Putting this together yields for the force in Eq. ( 30),

f (j) i = -∇ i Φ(x) δH δΦ(x) + ∇ i H + ∂H ∂X i -∇ j (∇ i Φ) ∂H ∂∇ j Φ . ( B3
)
We can also identify -∂H ∂X i with the force density in direction i acting on the object described by the coordinate X. This gives Eqs. (34) and (35) in the main text.

In order to demonstrate the agreement between Eq. ( 34) and Ref. 26, we use

∇ i Φ(x) δH δΦ(x) = Φ(x)∇ i δH δΦ(x) + δH δΦ(x) ∇ i Φ(x) (B4)
to arrive at

f (U) i = ∇ j δ ij H -∇ i Φ ∂H ∂∇ j Φ - δH δΦ(x) ∇ i Φ(x), (B5)
which is identical to Eq. ( 18) of Ref. 26.

In global equilibrium, one can prove 26 that the body force vanishes exactly and, furthermore, that

∇ i Φ(x) δH δΦ(x) eq = 0 = Φ(x)∇ i δH δΦ(x) eq . ( B6
)
where we introduced C(k, t) = e -Dk 2 t . With this, Eq. ( 60) in the main text is found, and Eq. ( 61) (for i j and k l) via σij (k, t) σkl (q, 0) = 1 4 1 0 dλdλ dpdr (2π) 6 Ãij (λkp) Ãkl (λ qr) φ(p, t) φ(kp, t) φ(r, 0) φ(qr, 0) . (C11)

APPENDIX D: EFFECTIVE HAMILTONIAN OF REF. 56

The effective Gaussian Hamiltonian for the system of Brownian particles proposed in Ref. .

(D3)

The right-hand side of Eq. (D3) is to be understood in the sense of inverse operators. The non-trivial part in the Hamiltonian in Eq. (D1) is the direct correlation function c [START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] , which is an important and well-studied object in the theory of liquids. [START_REF] Hansen | Theory of Simple Liquids[END_REF] The corresponding Model B equation of motion reads 56

∂φ ∂t = D k B T ∇ • ρ ∇ δH δφ + ∇ • 2D ρ η , = R∆φ + ∇ • 2D ρ η , ( D4
)
with noise correlations given in Eq. (11). With this, Eq. (68) of the main text is found for small wave-vectors.

  56 (see Ref.67 for a partly related approach) isH = 1 2 dxdyφ(x)∆(x, y)φ(y). (D1)This form contains a quadratic potential which plays the role of an effective interaction potential between densities,∆(x, y) ≡ k B T 1 ρ(x) δ(xy)c (2) (x, y) . (D2)As discussed in Ref.56, the Hamiltonian in Eq. (D1) yields the correct result for the first two moments of the fluctuating field φ in equilibrium, i.e., φ eq = 0 and φ(x)φ(y) eq = 1 ρ(x) δ(xy)c[START_REF] Larson | The Structure and Rheology of Complex Fluids[END_REF] (x, y)-1

TABLE I

 I

	. Observables relevant for this manuscript. Note that the mean den-
	sity, the two-body density, and the pair correlation can be evaluated both in
	and out of equilibrium.
	Symbol	Meaning
	ρ(x, t)

ACKNOWLEDGMENTS

This research was supported by the DFG Grant No. KR 3844/2-1 and the ANR project FISICS. A.S. acknowledges funding through a PLS fellowship of the Gordon and Betty Moore foundation. C.M.R. gratefully acknowledges S. Dietrich for financial support. V.D. thanks M. Schindler for interesting discussions.

APPENDIX A: DERIVATION OF THE STRESS TENSOR FOR ANY REALIZATION OF THE DENSITY FIELD [EQ. (17)]

Here we show that the divergence of Eq. ( 17) equals Eq. ( 14) under minimal assumptions and without any thermodynamic averaging. An alternative derivation can be found in Ref. 66. Making the trivial decomposition into an ideal part and the interacting term, we can write

where σ i is the ideal gas contribution,

The remaining non-ideal term due to inter-particle interactions can be written as

where we used a translation of the integration variable, and in the second line we have symmetrized the integral over the variables r and r. At this point it is important to choose interaction potentials obeying the reflection symmetry ∇V (r) = ∇V(r)

[this is clearly satisfied by isotropic potentials with V (r) = V (r)], for which one may write

This leads to

we can write

where the gradient should be taken with respect to x. Putting this together for isotropic potentials then yields

The equality of the divergence of Eqs. ( 17) and ( 14) then follows.

APPENDIX B: RELATIONS USED FOR THE DERIVATION OF THE STRESS TENSOR IN THE FIELD THEORY

The following relation is useful for derivation of Eq. ( 34) (noting that the Hamiltonian depends on the field and its gradients):

APPENDIX C: CORRELATIONS OF AN IDEAL GAS

For the computation of the stress tensor correlations to leading order in potential V, we require the time-dependent correlation functions of the ideal gas, which can be computed exactly. [START_REF] Velenich | [END_REF] We write the Fourier transform of the density field for N ideal particles as

where x µ0 denotes the position of particle µ at time t = 0 and X µ (t) is its subsequent displacement at time t. In this approximation, all the X µ (t) are independent Brownian motions with diffusion constant D. Two averages are now taken. Firstly, the average over the initial coordinate x i0 takes the form

where V is the volume of the system. The second average is over the Brownian motions X µ (t). The first two correlation functions are ρ(k, t) id = ρ(2π) 3 δ(k) (C3) and ρ(k, t) ρ(q, 0) id = (2π) 3 ρδ(k + q)e -Dq 2 t + (2π) 6 ρ 2 δ(k)δ(q). (C4)

From these correlations, we can get the correlations of the fluctuations. First, we note that the decomposition ρ(x, t) = ρ + φ(x, t) reads in Fourier space ρ(k, t) = (2π) 3 ρδ(k) + φ(k, t). We deduce that φ(k, t) id = 0, (C5)

Equation ( C6) is used to get Eq. ( 51) from Eq. ( 50) in the main text.

The three-point function is found similarly. After a few computation steps, we get

The three-point correlation of the fluctuations reduces to Eq. ( 56) in the main text. A useful relation is

with inverse Fourier transform

We also note that Ãij (k) = Ãij (k). The following relation reflects the fact that a particle cannot exert forces on itself:

With these relations, one may demonstrate that Eq. ( 56) gives no contribution to Eq. ( 55). The four-point correlation function needed for evaluation of Eq. ( 59) is found to be ρ(k, t) ρ(k , t) ρ(q, 0) ρ(q , 0) id = ρ(2π) 3 δ(k + k + q + q )C(k + k , t) + (2π) 6 ρ 2 δ(k + k )δ(q + q ) + δ(k + q)δ(k + q )C(k, t)C(k , t) + δ(k + q )δ(k + q)C(k, t)C(k , t) + (2π) 6 ρ 2 δ(k + k + q)δ(q )C(k + k , t) + δ(k + k + q )δ(q)C(k + k , t) + δ(k + q + q )δ(k )C(k, t) + δ(k + q + q)δ(k)C(k , t) + (2π) 9 ρ 3 δ(k + k )δ(q)δ(q ) + δ(k + q)δ(k )δ(q )C(k, t) + δ(k + q )δ(k )δ(q)C(k, t) + δ(k + q)δ(k)δ(q )C(k , t)

+ δ(k + q )δ(k)δ(q)C(k , t) + δ(q + q )δ(k )δ(k) + (2π) 12 ρ 4 δ(k)δ(k )δ(q)δ(q ) ,