Active Fluids Within the Unified Coloured Noise Approximation

Active matter is made of active particles which are able to convert energy from the environment into directed persistent motion. They can be modelled by stochastic differential equations subject to persistent noise. Run and tumble and active Brownian particle (ABP) models have been first proposed and still are considered closer to experimental observations but do not allow for much analytical progress. The Gaussian coloured noise (GCN) model, introduced as a time coarse-grained version of the ABP can be tuned to have the same variance of the active force as the ABP, which leads to a simpler analytical treatment. Finally, the UCNA can be considered as a Markovian reduction of the GCN. We give a simple derivation of the governing equation and analyse some of its recent applications ranging from the study of the swim pressure, its relation to the mobility, to the state induced by a moving object.

Run and tumble [4,5] and active Brownian particle (ABP) models [6] have been initially proposed to interpret experiments conducted on bacterial suspensions. More recently, the Gaussian coloured noise (GCN) model has gained a lot of attention. It was introduced with the idea of capturing the peculiar aspect of run and tumble and ABP models (i.e. the persistence of the trajectories of the active particles) and of reducing their mathematical complexity. In the GCN the components of the active force have a Gaussian distribution and are exponentially correlated in time with a characteristic time, τ . By applying to the GCN model an adiabatic elimination of the fast degrees of freedom one obtains the UCNA [7,8].
The UCNA [9] has the special property that its configurational steady state distribution is known, and that many stationary properties can be estimated. Employing this approximation, we present a description of a model of N mutually interacting active particles in the presence of external fields and characterise its steady state behaviour. Within the UCNA, we show that it is possible to develop a statistical mechanical approach similar to the one employed in the study of equilibrium liquids and to obtain the explicit form of the many-particle distribution function by means of the multidimensional unified coloured noise approximation. Such a distribution plays a role analogous to the Gibbs distribution in equilibrium statistical mechanics and provides a complete information about the microscopic steady state of the system. From here we develop a method to determine the oneand two-particle distribution functions in the spirit of the Born-Green-Yvon (BGY) equations of equilibrium statistical mechanics [10]. The resulting equations which contain extra-correlations induced by the activity allow determining the stationary density profiles in the presence of external fields, the pair correlations and the pressure of active fluids. In the low-density regime we obtain the effective pair potential φ acting between two isolated particles separated by a distance, r, showing the existence of an effective attraction. We apply the equations to different problems ranging from the study of the swim pressure, its relation to the mobility, to the investigation of the stationary state induced by a moving object in a "bath" of active particles.
Before closing this short introduction, we mention the fact that the UCNA method has been applied to the study of the effect of self-propulsion on a meanfield order-disorder transition [11]. Starting from a φ 4 scalar field theory subject to an exponentially correlated noise, the UCNA allows us to map the non-equilibrium active dynamics onto an effective equilibrium one. One can study the evolution of the second-order critical point as a function of the noise parameters: the correlation time, τ , and the noise strength, D. Our results suggest that the universality class of the φ 4 equilibrium model remains unchanged.

The Genesis of the UCNA Model of Active Particles
In order to understand the physical motivations of the model we shall discuss, it is necessary to give a brief historical account. In modern times, H.C. Berg was the first to introduce a model to describe the motion of bacteria in a viscous medium at small Reynolds number, the so-called run and tumble model, where the bacteria swim with constant velocity until a random tumble event suddenly decorrelates the orientation [4]. The active Brownian particles (ABP) model introduced to make analytical progress describes particles swimming at fixed speed u that rotates by slow angular diffusion. The two models have been shown to possess the same coarse-grained fluctuating hydrodynamics by Cates and Tailleur [12]. An advantage of the ABP is the possibility of taking into account external fields acting on the bacteria such as obstacles or gravity and interactions among them. For the n-th particle one hasṙ n (t) = v 0 e n (t) − ∇ n U γ , (8.1) where U represents the total potential energy of an N particle systems, whereas γ v 0 e i is the so-called active force, whose modulus is fixed, but whose direction e n (t) changes in time by rotational diffusion according to the laẇ e n (t) = D r η n × e n (t), (8.2) where η n (t) are Gaussian distributed with zero mean and have time correlations In spite of the great progress achieved using the ABP, the so-called active Ornstein-Uhlenbeck (AOU) or Gaussian coloured model has gained a great popularity because it has a simpler mathematical structure and lends itself to some analytical treatments due to the Gaussian character of the fluctuations of the active force [13]. The governing equations of such model are very similar to Eq. (8.1) with the difference that the active force γ v 0 e n (t) is replaced by γ u n (t), where the components of u n (t) vary between −∞ and ∞.
Within the AOUP we can obtain a series of useful results and in some cases we can solve exactly the equations, As, for instance, in the case of harmonic potentials where the equilibrium distribution is known. In the free-particle case U = 0 the free mean squared displacement is (r(t)−r(0)) 2 = 2Dτ [t +τ (1−e −t/τ )]. Thus, a free particle moves ballistically with typical speed v = √ D/τ at short times (t τ ) and diffusively with diffusion constant D at long times (t τ ). The typical distance travelled by a particle during a ballistic flight is the persistence length L = √ Dτ . If one observes the system on scales larger than L its properties will be almost indistinguishable from those of a system subject to standard thermal noise with an effective temperature T = Dγ .

The Unified Coloured Noise Approximation (UCNA)
In the following, we consider the evolution equation relative to the GCN and from this we shall derive the UCNA equation. For the sake of simplicity, we introduce the vector x of components x i of index i ≡ (α, n), where α is the Cartesian component associated with the coordinate of the n-th particle. We first differentiate w.r.t. time Eq. (8.3) and eliminate the active force γ u n (t) using Eq. (8.4). The resulting equation has the form of an underdamped Langevin equation with space dependent friction matrix: Neglecting the acceleration term in Eq. (8.6) we shall obtain the so-called unified coloured noise approximation, which is analogous to the Kramers to Smoluchowski reduction and is exact in the limits τ → 0 and τ → ∞. One can derive the UCNA equation by the original method of Hänggi and Jung [7]: on a new time scale s = tτ −1/2 one can recast the Langevin equation into the form If det ij is positive definite the damping is large for both small and large correlation times τ and in both cases one can set d 2 x i ds 2 = 0 and obtain a Markovian approximation of the coloured noise process of the form (8.9) which is to be interpreted in the Stratonovich sense.

Kinetic Approach
It is straightforward to write the equation of evolution for the N -particle probability distribution density of positions x associated with the overdamped limit d 2 x i dt 2 = 0. It reads It is, however, instructive to derive such an equation from a kinetic argument. We consider Eq. (8.6), define the velocity variable v i =ẋ i and write the following (stochastically equivalent to Eq. (8.6)) Kramers equation describing evolution of phase-space distribution of N particles, N (x i , v i , t): This kind of equation occurs in the study of colloidal solutions and is treated by multiple time scale methods. In general, we cannot solve Eq. (8.11) which involves both the velocity and the position variables. However, we can attack the problem by assuming that the velocity degrees of freedom evolve much faster than the positional degrees of freedom. This type of assumption is done when one reduces the Kramers phase-space equation to the Smoluchowski configurational equation. In fact, we make the ansatz that the phase-space distribution factorises in a spatial part and a velocity part. We construct a time-independent trial phase-space distribution having a factorised form: where is the conditional velocity distribution when the particles positions are fixed at x. P N (x) corresponds to the distribution of the particles and is the marginalised distribution giving the distribution of positions of the particles regardless their velocities: . (8.13) In order to determine P N we integrate Eq. (8.11) with respect to all velocities and obtain the continuity equation relating the probability density P N and the probability current J i where the current J i (x, t) is the dN-dimensional vector: After multiplying Eq. (8.11) by v i and integrating over the dN velocities, we obtain the momentum balance equations where p ik (x, t) ≡ dvv i v j N Eqs. (8.14) and (8.16) are an Nd + 1 system which is not closed because one does not know the explicit form of the tensor p ik (x, t) in terms of P N and J i . A simplifying ansatz is to assume that the velocities have a local distribution similar to the one they would have in an equilibrium system; this is the following multivariate Gaussian distribution: It is important to notice that the variance depends on the positions of the particles in contrast with equilibrium systems and is consistent with the fact that the friction is position dependent. Within the Gaussian ansatz we can rewrite the balance Eq. (8.16) as Finally, we assume that the time derivative of the current vanishes on a faster time scale than the time derivative of the density so that dropping the time derivative in Eq. (8.18) and expressing J i in terms of P N in Eq. (8.14) we obtain Eq. (8.10) The advantage of such a derivation is that we have obtained not only the distribution function of positions, but also the approximate form of the distribution of the velocities of the particles. The latter is peculiar because it depends on the positions of all the particles at variance with the equilibrium case.

Stationary Solution in the Absence of Current
Let us consider the configurational distribution function P N (x) in the steady state associated with Eq. (8.10). In order to realise the steady state there are two possibilities, namely when the divergence of the probability flux vanishes, i ∂ i J i = 0, or when the flow, J, itself vanishes. Since only the configurational space is considered in such a reduced description and the positional variables, x i , are even under time-reversal transformation, the condition J i = 0 for arbitrary i is equivalent to the detailed balance condition [14]. In detail, if the matrix −1 is non-singular, which can be rewritten as: The detailed balance implies a stronger condition than the one represented by having a stationary distribution, since it implies that there is no net flow of probability around any closed cycle of states. Such a situation is no longer true when we consider the phase-space probability density of the original GCN problem, as discussed in Appendix 2.
From the above equation, one can find an exact expression for the probability density, which reads 22) where Z N is a normalisation constant. Such a formula, in principle, fully describes within the unified colour approximation the steady state distribution of a system of interacting particles subject to coloured noise. In the white-noise limit τ → 0 the formula reduces to the Boltzmann distribution corresponding to the potential U . For finite values of τ , instead, the distribution maintains a Boltzmann-like distribution but with the effective potential given by Eq. (8.22). The presence of the additional terms ∂U(x) ∂x k 2 and ln det has repercussions in the form of the steady state configuration. Such a form of distribution is at a first glance surprising since in equilibrium systems, energy is exchanged reversibly with the environment and the form of P N (x) is determined by the potential and the temperature of the environment. On the contrary, in nonequilibrium systems, energy is exchanged irreversibly with the environment and in general there is no one-to-one correspondence between potential and P N (x). The vanishing of all components of the probability current (see Eq. (8.20)) is tantamount of the existence of the detailed balance condition, i.e. of the microscopic reversibility in the dynamics of the active system. This is reflected in the Boltzmannlike form of the distribution function. One may ask whether this is an artefact of the UCNA treatment of the dynamics or is a genuine property of the system. As we shall discuss below, by considering an elementary case, the detailed balance condition is violated by the original GCN dynamics by terms proportional to the persistence time, τ .
For a total potential, U , consisting of the sum of purely repulsive pair potentials the overall result is to create a sort of effective attractive potential among the particles. The origin of such an attraction can be understood as follows: the drag force on each particle is determined by the bare friction with the solvent medium plus an additional contribution stemming from the interactions. The non-equilibrium force is an attraction between self-propelled particles causing them to cluster. In the case of J i = 0 and i ∂ i J i = 0, it is not possible in general to obtain explicit solutions apart from some special cases which will be discussed later.

Fox Approximation
The approximate treatment obtained by applying the UCNA method is not unique. An alternative method has been put forward by Fox [15] who employed functional calculus in order to derive the effective equation for the distribution function P (x, t) corresponding to the GCN model. The resulting equation of evolution is valid in the small τ regime and has been applied to active fluids by Farage et al. [16]. It reads Interestingly, the Fox and the UCNA approaches in the case of a single coloured noise yield the same steady state distribution function, whereas the approach to such a solution is different in the two cases. In the case where the particles are subject to different types of noises, each characterised by its own relaxation time, the UCNA approximation does not give the correct equation of motion even in the small τ limit, whereas the Fox method correctly reproduces such a limit. Therefore, in order to describe mixtures of active particles or of passive and active particles it is convenient to apply Fox's approach in spite of the fact that it only describes the small τ regime [17,18].

Entropy Production in UCNA
The detailed balance requires that the probability of making a transition forward in time equals the probability of making the reverse transition, backward in time, when the system is in the steady state. It is easy to verify that within the UCNA approximation the condition of detailed balance holds if the probability current vanishes, J = 0 in the steady state. The vanishing of J implies the existence of an effective potential U eff which fully determines the distribution. We shall see in Sect. 8.7 that this is not the case when J = 0. This is the reason why the UCNA steady state distribution has a form similar to a Boltzmann distribution, although with an effective potential which depends on the persistence time.
A measure of the distance from thermodynamic equilibrium is provided by the entropy production, so that it is interesting to study such a quantity in the steady state of the UCNA evolution equation. To this purpose, let us consider the rate of change of the Shannon entropy (for the sake of simplicity we study the case with N = 1 of Eq. (8.10)).
where we obtained the last equality by partial integration. We decomposeṠ into two contributions:Ṡ whereṠ s is the entropy production due to irreversible processes occurring inside the system andṠ m is the entropy flux from the environment to the system. We shall show thatṠ s is positive definite, whereas S m can have either sign. In the steady state the rate of change of the entropy vanishes so thatṠ m = −Ṡ s . From Eq. (8.19) for (8.25) with the following probability current: where We can eliminate ∇P from Eq. (8.24) and obtain the following expression in terms of the current: We now identify the first term in Eq. (8.27) as an entropy production rate always non-negative, and the second terṁ with the entropy flux due to heat exchanges between the system and the surroundings and the temperature T = Dγ . We identify the heat flux with the average change of effective potential energy, U eff , of the system in the unit time evaluated as follows: After an integration by parts we obtain and by comparing Eqs. (8.31) and (8.29), we find the following relation: Finally, we haveṠ Notice that now the temperature entering the formula connectingṠ m and Q is uniform and given by T . Let us remark that due to the detailed balance condition bothṠ s andṠ m vanish in the steady state UCNA, showing that the UCNA method maps the underlying GCN non-equilibrium description into an equilibrium one. At variance with the UCNA, in the GCN bothṠ s andṠ m are non-vanishing in the steady state.

H-Theorem
The following calculation proves the approach to the stationary distribution in terms of the entropy functional. One sees immediately that the entropy fluẋ The quantityṠ s is nothing else but the rate of change of the Kullback-Leibler entropy, S KL ≡ − dxP (x, t) ln(P (x, t)/P steady (x)), which is positive due to the sign ofṠ s and vanishes at equilibriuṁ Thus the Kullback-Leibler entropy of the UCNA process is an ever increasing function and satisfies an H-theorem. The relative entropy S KL (t) is a functional of the non-equilibrium probability distribution and generalises the ordinary thermodynamic entropy which is defined for equilibrium states.

Born-Green-Yvon Hierarchy in the Steady State
We go back now to the multidimensional case and adopt indices to specify components and particles. We focus attention on the steady state properties of the system as described by the UCNA. We must remark that formula, Eq. (8.22) refers to N particles and therefore is not of practical use when the particles are mutually interacting. We need to derive from it expressions for the one-body and two-body distribution functions. The procedure is similar to the one employed in equilibrium statistical mechanics. We shall use the steady condition, Eq. (8.20), to derive a set of equations similar to the BGY hierarchy for distribution functions in equilibrium systems. The hierarchy becomes of practical utility in conjunction with a suitable truncation scheme in order to eliminate the dependence on the higher order correlations. In the following, the Cartesian components (from 1 to d) are identified by the indices α and β, and the particles are identified by Latin indices. The total potential is assumed to be the sum of the mutual pairwise interactions w(r − r ) between the particles and of the potential exerted by the external field u(r) : U(r 1 , . . . , r N N P N (r 1 , r 2 , . . . , r N ).
The equation for P (2) N has a structure similar to that of a standard equilibrium gas but for the term containing −1 α1,βn and unless we introduce approximations it is of little practical use.
We write 39) and we remark that in the limit of small (τ/γ ) the matrix −1 α1,βn can be approximated as [19]: where u αβ ≡ ∂ 2 u(r) ∂r α ∂r β and w αβ ≡ ∂ 2 w(r) ∂r α ∂r β . We substitute this approximation and recast Eq. (8.38) in terms of the n-th order density distributions ρ (n) (r 1 , r 2 , . . . , r n ) = N ! (N −n)! P (n) N (r 1 , r 2 , . . . , r n ) and find which represents the BGY equation for the pair density distribution ρ (2) . By integrating also over the coordinate 2 we find the BGY equation for the one-body density: (8.42) that in the limit of τ → 0 is just the BGY equation for the single-particle distribution function.
The r.h.s. of Eq. (8.42) contains the coupling to the external field and the socalled direct interaction among the particles, whereas the l.h.s. besides the ideal gas term contains a term proportional to the activity parameter.

Active Pressure
A natural way to define the pressure in a system of active spherical particles driven by coloured noise is by using the virial theorem which relates the virial of the external forces confining the particles in a given volume to the pressure exerted on the walls by the particles themselves. The forces exerted by the bounding walls of the container are macroscopically described as external pressure [20,21]. Each oriented area element dS exerts a force −p(r)dS so that wherep is the average pressure over the boundary surface, r is the position vector of the surface element and the last equality follows from the divergence theorem (∇ · r = d). Now, in order to evaluate the external force virial, Eq. (8.43), we multiply Eq. (8.38) by r α1 , integrate over r 1 and r 2 and sum over indices. After an integration by parts we obtain the following equation: where the forces are separated in two parts: wall and interparticle forces, F ext αi = − ∂u(r 1 ) ∂r α1 and F int αi = − k ∂w(r 1 −r k ) ∂r α1 , respectively, and the symbol · stands for an average over the stationary distribution P N .
In the case where the confining vessel has constant curvature one finds:p = p(r). In general, when the linear size of the vessel is much larger than the persistence length, the standard virial definition of pressure based on the assumption of the constancy of the pressure on the boundary of the system is correct [20,21]. In order to obtain a closed expression for the pressure, we write the term stemming from the internal forces as: and approximate the average of the trace of −1 as: where in the second equality we have used Eq. (8.40). We now write The first term in the r.h.s. of Eq. (8.48) represents an ideal gas-like contribution to the pressure, T N/V , also referred to as the swim pressure, due to the rotational degrees of freedom. The second and third term in Eq. (8.48) represent indirect interaction contributions, and take into account the slowing down of active fluids near a boundary and in regions of high density, respectively. The indirect interaction pressure involves the interplay between the rotational degrees of freedom and the interparticle forces and is a non-equilibrium effect. In fact, in the limit of τ → 0 the quantity Eq. (8.48) reduces to T t N/V , the ideal gas contribution to the pressure. Besides the virial method, for the UCNA there exist two other approaches to evaluate not only the pressure but also the surface tension. In the first of approach, these can be identified with the volume and area derivatives, respectively, of the partition function associated with the stationary non-equilibrium distribution. The second alternative method is a mechanical approach and is related to the work necessary to deform the system. The pressure is obtained by comparing the expression of the work in terms of local stress and strain with the corresponding expression in terms of microscopic distribution. This work is determined from the force balance encoded in the Born-Green-Yvon equation and can be used to obtain a formula for the local pressure tensor and the surface tension even in inhomogeneous situations. Nicely, the three procedures lead to the same values of the pressure, and give support to the idea that the UCNA partition function is more than a formal property of the system, but determines the stationary non-equilibrium thermodynamics of the model. For further details the reader may consult ref. [9].

Velocity Correlations
The kinetic derivation of the UCNA has shown that a system of active particles displays velocity correlations. Within the present treatment these correlations have been approximated by means of a Gaussian multivariate distribution whose variance depends on the potential.
We consider N interacting particles in 1d. We perform numerical simulations of systems with N = 1000 composed by GCN-driven particles interacting via the potential φ(x) = i>j (x i − x j ) −12 for several values of the density ρ = N/L, of D and τ .
The velocity variance depends on the configuration of the particles, so that by averaging it over the positions we obtain the overall velocity variance of a GCNdriven system: The v 2 computed numerically via Eq. (8.49) is plotted in Fig. 8.1b (full lines) as a function of the 1d density ρ = N/L of the system and for several values of D (at fixed τ ). In all these simulations we compute the variance ẋ 2 and report the results in Fig. 8.1 as connected symbols.
To test the validity of the Gaussian ansatz for the velocity distribution given by Eq. (8.17), we compute the average over positions in Eq. (8.49) directly from the coordinates obtained numerically, instead of using the theoretical P N of Eq. (8.22). This is plotted in Fig. 8.1 as dashed lines and follows well the numerical curves, although some expected deviation [7] is observed upon increasing D to very high values. If we assume a uniform density and long-ranged interactions (mean-field approximation) the velocity distribution Eq. (8.17) simplifies substantially since all the out-of-diagonal term of ∂ α ∂ β w(r i − r j ) are of order one and can be neglected with respect to the terms on the main diagonal that are of order N [19]. This yields the density-dependent variance: is the mean potential curvature integrated from the diameter σ and L = √ Dτ is the characteristic length of the active motion. In the last equality of Eq. (8.50) we have used the fact that, for a generic repulsive potential, σ corresponds roughly to the distance where the interaction force balances the selfpropulsion force (i.e. |w (σ )| ≈ γ v = γ √ D/τ ) and w 2 = w (∞) − w (σ ) ≈ γ √ D/τ . This is plotted in Fig. 8.1 as a thick line for the largest D and follows well the data when L is large. To first order in τ we obtain the results plotted as dasheddotted lines in Fig. 8.1 and the theory compares well with the numerical simulations. However, by fixing τ and increasing D this approximation deviates strongly from the simulations. We can now derive an expression for pressure by a kinetic argument by identifying it with

Active Elastic Dumbbells
Let us consider N mutually noninteracting elastic dumbbells, i.e. two point particles bound together by an elastic spring of constant α 2 , moving in a vessel represented by a harmonic weak confining potential, of spring constant ω 2 [23,24]. Such a model, similar to the harmonic trap model [25][26][27], was proposed long ago by Riddell and Uhlenbeck. It contains the minimal ingredients to observe the competition between internal forces and confining potential and can be solved without introducing further approximations. The potential energy reads [28]: with w(r) = 1 2 α 2 r 2 . By setting u(r) = k 2 r 2 L 2 , one introduces a volume dependence in the spring constant associated with the confining potential and for simplicity of notation we shall use ω 2 = k L 2 . The virial pressure is obtained by applying the general formula Eq. (8.47) By simple algebraic manipulations we find . (8.53) In the limit of τ → 0 and α → 0 the pressure reduces to the expected ideal gas pressure of a system of 2N noninteracting particles in a vessel of volume L d . On the other hand, one can see that the pressure decreases with increasing values of τ , i.e. if the persistence length L = √ Dτ exceeds the typical size of the vessel the particles do not explore the whole space of the vessel, but remain localised at the bottom.

Pressure of N Noninteracting Active Particles Surrounded by Harshly Repulsive Walls
As a second example we consider an assembly of N noninteracting active particles constrained in a region of space near the origin by a spherically symmetric external potential in three dimensions. Using Eq. (8.42) one can derive the following exact formula expressing the mechanical balance condition: where the components of the pressure tensor normal (N) and tangential (T) to the walls are p N (r) = T s ρ (1) , (8.55) and p T (r) = T s ρ (1) The density profile according to Eq. (8.22) can be written explicitly as: (8.57) so that we can fully determine the components of the pressure tensor.

Active Particles in a Time-Dependent Potential
The results presented in the previous sections concern static cases, where the external potential is constant in time. In this section we address the interesting issue of a time-dependent external potential [29]. In particular, we shall consider a shifting potential U(x, t) = U(x − ct) in one dimension, moving at constant speed c and inducing a stationary current in the system. The potential barrier interacts with a fluid of active particle, see the sketch depicted in Fig. 8.2. The effect of a moving potential on a particle fluid is a general problem in modelling the motion of a driven obstacle in a medium, in several different fields, such as in the active microrheology of colloidal systems or in the translocation dynamics of polymer chains through nanopores. Fig. 8.2 Sketch of the system: a potential barrier moves at constant velocity c in a channel with a noninteracting active particle fluid, producing a density profile which is non-uniform along the x direction [29] In the case of the GCN model, we will show that the coupling of self-propulsion (namely a finite persistence time τ ) with the stationary current gives place to an effective dynamical potential, which vanishes in both the limits of c → 0, and τ → 0 (passive particles). The main physical effect we observe in this model, which is accounted for by a generalised UCNA scheme, consists in a much enhanced accumulation of active particles at the interface fluid/obstacle, with respect to the static case or with respect to the behaviour of a passive particle fluid.

Effective Potential
As shown in Sect. 8.3, the first step in deriving UCNA equations is to take the time derivative in Eq. (8.3). For the sake of simplicity, we consider a one-dimensional system and a shifting time-dependent potential of the form U(x, t) = U(x − ct). We haveẋ By comparing with Eq. (8.6) we note that a new term appears in Eqs. (8.61): an effective force F * (x), which reduces to −dU/dx when c = 0. As we shall show, this additional contribution in the force term due to the finite velocity of the obstacle c > 0 is responsible for new dynamical effects.

Dynamical UCNA and Particle Density Profile
In order to show how these effects can be described within a generalised UCNA scheme, it is useful to consider the associated Fokker-Planck equation. It is time saving to adopt non-dimensional variables for positions, velocities and time, and rescale forces accordingly. We define v T = √ D/τ , measure lengths using the characteristic length, , of the potential, and introduce the following nondimensional variables: where ζ plays the role of a non-dimensional friction. To lighten the notation we shall drop the bar over the non-dimensional variables without incurring in ambiguities.
For the probability distribution of position and velocity (y, v) we thus obtain 63) where we have introduced the shifted variable y = x − ct. We, now, look for an approximate solution to this equation. We start by eliminating the v dependence of the phase-space distribution (y, v), by multiplying by powers of v and integrating w.r.t. v. Thus, one obtains a set of coupled first order ordinary differential equations, the so-called Brinkman hierarchy, whose first two members are the continuity equation and the momentum balance equation, respectively: Here we have introduced the density ρ(y), the current J (y) and the momentum current (y), defined as: whereρ is a constant such that the solution is periodic at ρ(L) = ρ(−L), where 2L is the system size. As we shall see later, for large systems L l,ρ ≈ ρ(±L) and the current is almost vanishing at the boundaries.
From the analysis of the case of static potentials, discussed in the previous sections, we know that the solution of Eq. (8.63) in regions where F * (y) = 0 and (y) = 1 can be written as (see Eq. (8.17)): is a Hermite function of zero order. Indeed, by substituting the form Eq. (8.70) in Eq. (8.63) (with F * = 0), we obtain a solution provided ρ(y) satisfies the following condition: Next, we insist in looking for a solution of Eq. (8.63) even in the region where F * (y) = 0 of the form: where we have introduced the following (non-uniform) Hermite function, which is position dependent through the trial function β(y): Substituting now the trial distribution Eq. (8.73) into Eq. (8.63), we get If we continue the projection procedure beyond the first order in (v − c) there will be an error in the equation for the second moment, which becomes inconsistent with the value of the second moment imposed by the trial distribution (which, in fact, is already fixed by the trial form and therefore does not contain enough parameters to satisfy the extra conditions.) The ansatz for the phase-space distribution gives the following expression for the momentum flux: .
We empirically set β(y) = (y) in the regions where (y) ≥ 0, and β(y) = 0 otherwise. Then, the expression (8.84) can be evaluated numerically and in Fig. 8.3 we compare the analytical prediction with numerical simulations, in the case of the following external potential: characterised by the steepness 1/ξ .

Average Drag Force
Our analytical approach allows us to obtain an estimate for the average drag force exerted by the active fluid on the moving wall, defined as: The comparison of the analytical prediction with numerical simulations for the drag force is shown in Fig. 8.4. Note the non-monotonic behaviour of the force-velocity relation is characterised by a maximum value of the force F max for a particular velocity c * . In order to highlight the new physical effects arising due to the coupling of selfpropulsion and a stationary current, it is useful to compare the behaviour observed in the active particle model with the one obtained in the case of a moving potential in a (passive) thermal bath. In the latter situation, the noise term acting in the stochastic equation for the particle velocity is a delta-correlated noise of amplitude 2/ζ . In the right panel of Fig. 8.4 we show the maximum value of the drag force F max as a function of 1/(ξ ζ ) in both models. The qualitative difference between the two behaviours relies on the observation that in the active case the average drag force can increase indefinitely by reducing the parameter ξ , which characterises the steepness of the travelling potential.

Conclusions
In this chapter, we have reviewed the recent developments of the theory of active particles driven by coloured noise within an approximate scheme, the UCNA. Such a method has the great advantage of providing predictions and equations much simpler with respect to other methods. The reason is that the adiabatic approximation at the basis of the method eliminates the faster degrees of freedom, the velocity in the case of the GCN and of the ABP.
Let us define the temperature T = D/τ and the local temperature θ(x) = T / (x). Now, the total time derivative of the entropy Eq. (8.88) can be written as the sum, s(t) =ṡ s (t) +ṡ m (t) (see Sect. 2.4). Explicitly, we find after integrating by parts Eq. (8.88) the following expressions: and the entropy fluẋ The dimensional form of the total energy (x, v) and of the heat flux are, respectively: and It is suggestive to rewriteq with a local density of heat flux defined aṡ where n(x, t) = dvp(x, v, t) and n(x, t) v 2 x = dvv 2 p(x, v, t) , with v 2 x the mean squared velocity at given position. Expression Eq. (8.97) represents an interesting connection between the local entropy production of the medium (or entropy flux) and the local heat flux divided by the same local temperature θ(x) = T / (x) featuring in the approximate detailed balance solution, Eq. (8.103). Define the dissipative components of the current as: It is clear that in the GCN in the steady state while the time derivative of the entropy vanishesṡ = 0, its two contributions are not necessarily zero.
vĨ v (x, v, t) (8.100) s s is an entropy production rate which is always non-negative, whileṡ m is the entropy flux due to heat exchanges between the system and the surroundings and can have either signs.

Appendix 2: Absence of Detailed Balance Condition in the GCN
As in Appendix 1, let us consider the elementary case of a single active particle in one dimension driven by Gaussian coloured noise. The probability current, I(x, v), is the two dimensional vector, comprising both reversible and irreversible contribution given by Eq. (8.90 p s must be the product of a function of the position, π(x), times a "local" Maxwellian whose velocity variance is also position-dependent: and conclude that a function π(x) satisfying Eq. (8.104) only exists when (x) is a constant 0 . This condition occurs for τ → 0, which is the equilibrium limit of the model, or when the potential is a linear or parabolic function of x. In conclusion, apart from the special case 0 , the Kramers equation (8.11) does not satisfy the detailed balance condition. However, we can determine an approximate steady solution consistent with the UCNA under the form: In spite of the fact that such trial solution is not divergence-free, i.e. divI = 0, the first three velocity moments, obtained by multiplying the divergence by (1, v, v 2 ), respectively, and integrating w.r.t. v, vanish, so that in this subspace the zero divergence condition holds. Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence and indicate if changes were made.
The images or other third party material in this chapter are included in the chapter's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.