Force renormalization for probes immersed in an active bath

Langevin equations or generalized Langevin equations (GLEs) are popular models for describing the motion of a particle in a fluid medium in an effective manner. Here we examine particles immersed in an inherently nonequilibrium fluid, i.e., an active bath, which are subject to an external force. Specifically, we consider two types of forces that are highly relevant for microrheological studies: A harmonic, trapping force and a constant,"drag"force. We study such systems by molecular simulations and use the simulation data to derive an effective GLE description. We find that, in an active bath, the external force in the GLE is not equal to the physical external force, but rather a renormalized external force, which can be significantly smaller. The effect cannot be attributed to the mere temperature renormalization, which is also observed.


I. INTRODUCTION
Microrheology focuses on linking the microscopic dynamics of a system to its macroscopic response.The goal of microrheological studies is then to investigate the behavior of an immersed and experimentally detectable probe particle to gain knowledge about the surrounding medium which is made up of smaller particles [1][2][3].Microrheological techniques have been used to study a variety of different mediums ranging from complex fluids [4][5][6][7] to granular matter [8][9][10][11][12].
In general, microrheological studies can be classified in two categories: passive and active [13,14].Passive microrheology monitors the diffusive motion of an immersed probe to extract information about fluctuations and transport properties of the medium.Active microrheology relies on subjecting a probe to an external force to determine the properties of the immersive medium.This is true for both experimental and computational microrheological studies.The external forces often take the form of either constant, drag forces [15][16][17] or harmonic, trapping forces [18].Microrheological studies employing these forms of external forces have also recently been applied to probe the properties of active systems [19][20][21][22][23][24].
Since microrheology focusses on investigating a small number of immersed probes, usually described by their positions and velocities, the extracted information will correspond to a strongly coarse-grained description, in which the surrounding medium is replaced by effective friction coefficients and fluctuations.In previous studies, for example, the system of a passive probe immersed in an active bath and subject to a harmonic potential has been mapped onto an overdamped Langevin equation [25][26][27].However, in other studies, it has been shown that, under certain conditions, the dynamics of a probe immersed in an active bath cannot be assumed to be Markovian [28,29].In this case, the Langevin equation is no longer sufficient and, instead, a generalized Langevin equation (GLE) is required to capture the system dynamics in a coarse-grained model [30,31].
A simple and popular Ansatz for such a GLE is to assume the form [17,[31][32][33][34][35][36] M V(t) = f (R, t) − t T ds K(t − s)V(s) + Γ(t), (1) where M is the particle mass, V(t) its velocity, R(t) its position, f (R, t) describes the drift induced by the external forces (in passive microrheology, f ≡ 0), K(t − s) is a memory kernel, and Γ(t) a colored stochastic noise with zero mean.Vroylandt and Monmarché [35,37] have recently shown that it is possible to derive this form of the GLE -possibly with a space-dependent memory kernel -from first principles by projection operator techniques [30,38] for Hamiltonian systems at equilibrium.In the absence of external forces (at f (R, t) = 0), the noise term in the GLE as obtained from the projection operator obeys a so-called second fluctuation dissipation theorem (2FDT) ⟨Γ(t 0 )Γ(t)⟩ eq = k B T K(t − t 0 ).The effect of including external forces in a GLE description [36,[39][40][41], non-stationary dynamics [42,43], and stochastic microscopic dynamics [44,45] has also been discussed using first principle theory.Important conclusions from these studies are that the 2FDT is no longer generally valid in such cases: namely, in the presence of non-linear or stochastic interactions [45,46], external forces, or general non-equilibrium settings [47].Furthermore, additional terms in this relation can emerge [35,36,45,48].
In this work we will take a complementary approach and use the GLE (1) as a starting point to extract effective equations of motion for the coarse-grained variables from microscopic simulations.Thus, our goal is not to derive a GLE model (1) from a microscopic Hamiltonian (which may not be rigorously possible), but to map an ensemble of microscopic trajectories onto the model (1) such that important static properties, such as the distribution of particle positions, and important dynamic properties, such as the average drift and velocity correlation functions, are reproduced by the dynamical system defined by the GLE.All averages, ⟨...⟩, will therefore refer to the average over this, possibly non-equilibrium, ensemble.One should note that one has some freedom as to how to distribute the total force between f , K and Γ [17,30,49].At equilibrium, a popular Ansatz [33,48,50,51] is to assume that f subsumes all reversible interactions (which also include medium-mediated interactions) and therefore correspond to the force emerging from static coarse-graining.The other two terms then characterize the thermal coupling with the medium; i.e., they can be seen as a thermostat [52] with colored noise.
What happens in active systems far from equilibrium?This question has been intensely discussed in the literature [19,22,23,47,49,[53][54][55][56][57], often focusing on the question of whether the 2FDT is still fulfilled.However, as mentioned above, one has some freedom as to how to distribute the net bath forces between the different terms of the GLE.In the absence of drift forces, it has been shown [17,58] that many memory reconstruction methods [18,[58][59][60][61][62] will yield GLEs that automatically satisfy a generalized 2FDT, (in tensorial notation).Putting it differently: One can impose an orthogonality condition between the coarsegrained variables and the stochastic noise ⟨V (0)Γ(t)⟩ = 0, where t = 0 refers to some arbitrary "initial time", as an additional condition on the memory kernel of the force-free GLE.Mathematically, this implies that the memory kernel and the velocity autocorrelation function depend on each other via a so-called Volterra equation [31,58].This imposed constraint is reminiscent of the typical projector used in projection operator techniques, which would lead to similar orthogonality conditions [30,36,37,45].We will impose this condition throughout the manuscript for all our analytical and numerical analysis.However, we would like to emphasize that we do not apply any projection operator techniques, but rather solely apply this orthogonality condition.
In the present paper, we then address the following question: How does the active bath affect the drift force f on the probe in the presence of external drag forces?As prototype test systems, we consider a passive probe in a homogeneous bath of active Langevin particles [63][64][65] subject to two typical types of external forces F ext : A constant drag force and a harmonic trap force.
Our main result can be summarized as follows: We can obtain a consistent GLE mapping of this system, but only if we renormalize the force: The drift force differs from the applied force by a factor α which is much smaller than one.The factor does not depend on whether the potential is linear or harmonic.Physically, this means that the active bath particles help the probe particle to access high potential regions not only by increasing its kinetic energy (i.e., effectively increasing its temperature), but also by actively pushing it into these regions.
By a mathematical analysis, we identify a generalized equipartition theorem which characterizes the relation between the renormalized force and the kinetic energy of the particle.The analysis also shows that the generalized 2FDT in our model may acquire an extra term in nonequilibrium systems, consistent with Refs.[36,54,55], which however vanishes within the error in the systems considered here.
Our results are of high importance for the interpretation of active microrheology in active systems, as they show how to extract information about an active bath by observing transport of immersed tracer particles.Our results are particularly pertinent to systems in which the constituents are of non-negligible mass, i.e.where the constituents are meso-or macroscopic.Examples of such systems include synthetic active robots [66][67][68] and granular active matter [69][70][71][72].Whereas there exist previous studies on the microrheological properties of active granular matter [8][9][10][11][12], synthetic active robots is an emerging field and we are not aware of any studies of microrheological studies of such systems.

II. MODEL AND SIMULATION DETAILS
We consider a three dimensional system of a passive probe immersed in a bath of active Langevin particles (ALPs) of mass m and radius R, which propel themselves with a constant force F 0 subject to rotational diffusion with a diffusion constant D R [63][64][65][73][74][75][76][77][78].An analogous model in the overdamped limit (a probe immersed in a bath of active Brownian particles with mass m → 0) has often been focused on in previous studies [25,[79][80][81][82][83][84].These ABP fluids have a number of similar properties to ALP fluids, including the quadratic scaling of the effective temperature (or effective diffusion constant in the case of an ABP) with the bath activity and the velocity correlations [65].
The ALPs are coupled to a thermal bath with temperature k B T via a Markovian, Langevin thermostat, and they interact with each other and with the immersed probe by repulsive hard core interactions of the Weeks-Chandler-Anderson (WCA) type [85].Specifically, we follow Refs.[63][64][65] and set the rotational inertia of ALPs to zero for simplicity.The equations of motion for an ALP, n, in the bath, are thus given by where F 0 is the propulsion force of the active particle, e(t) is the orientation of the active particle, and γ = 6πηR is the damping constant for an ALP radius R in a thermal bath with viscosity η.The term ξ n (t) represents a stochastic force that mimics implicit collisions of the ALPs with thermal bath particles.These collisions are modelled as Gaussian white noise with mean zero and variance given by a fluctuation-dissipation relation The resulting translational diffusion coefficient of isolated ALPs is given by D T = k B T /γ.Finally, the terms −∇U WCA (r n − R) and − n̸ =m ∇U WCA (r n − r m ) in Eq.
(3) describe the WCA interactions with the probe and with all other ALPs, respectively.The time evolution of the orientation of the ALP, e(t), is governed by rotational diffusion, where N(t) is again Gaussian white noise with a mean of zero and a variance (another fluctuation-dissipation relation) Here D R is the rotational diffusion constant, which is given by D R = 3D T /4R 2 for a particle of radius R. The immersed probe only experiences the externally applied force and forces from interactions with the surrounding ALPs.Unlike the ALPs, it is not coupled to the thermal bath.Ref. [65] has shown that the behavior of a probe immersed in an ALP fluid remains qualitatively the same regardless of whether it is coupled to the thermal bath.In particular, the probe acquires an enhanced kinetic temperature which scales quadratically with the activity of the ALP fluid even in the presence of a thermostat.All simulations are performed using LAMMPS [86].The length, energy, and mass scales in the simulation system are defined by the Lennard-Jones (LJ) diameter σ, energy ϵ, and ALP mass m, respectively, which defines the LJ time scale t LJ = σ m/ϵ.We use truncated and shifted LJ potentials with the energy scale ϵ which are cut off at r c = 2 1 6 σ, resulting in purely repulsive WCA interactions as described above.The cubic simulation box has periodic boundary conditions in all three dimensions and a side length based on the desired density of the bath.The probe has a mass of M = 100m and is defined as a rigid body with a radius R p = 3σ.The body of the probe is constructed so that its surface is smooth, resulting in full slip boundary conditions for the LJ fluid.The active bath consists of ALPs with a mass of m ALP = 1m and a radius of R = 0.5σ.The number of ALPs in the bath is determined by the desired density of the bath.The parameters of the thermal bath are chosen such that k B T = 1 ϵ and η = 1 m/(σt LJ ), resulting in γ = 3πm/t LJ , D T = (3π) −1 σ 2 /t LJ , and D R = π −1 /t LJ .For the active bath, we choose an active force of F 0 = 50 mσ 2 /t lJ because Ref. [65] showed that this force is sufficiently large to induce non-equilibrium probe dynamics.but not large enough to cause motility induced phase separation in the active bath.
In the following, for consistency with previous work [65], we will often also use dimensionless quantities t = t γ/m, v = v m/k B T , F0 = F 0 1 γ m/k B T , and DR = D R m/γ.Dimensionless quantities will denoted with a hat.Density values and the spring constant of the harmonic potential, and all distances are given in LJ units.In these cases, the LJ units are explicitly written with the value.

III. RESULTS
The behavior of free particles in the ALP bath has been reported in previous work [65] for different bath densities and active forces F 0 .One important result of this study was that the bath transmits a kinetic temperature to the probe (k B T eff = M ⟨V 2 ⟩/3 in 3D), which, somewhat unexpectedly, can be significantly higher than the effective kinetic temperature of the ALP bath particles.In the present article, we will analyze the probe dynamics in the presence of external forces as described below.

A. Harmonic potential
We trap the immersed probe in a harmonic potential such that its equation of motion is: (7) where M is the mass of the probe, V(t) is its velocity, and U WCA (R − r n ) is the WCA potential due to an ALP particle, n, at position 2 is the harmonic, trapping potential with trapping constant k and trap center R 0 .We test three different values of k: k = 1, 5, and 10 in LJ units of ϵσ −2 .We map our system dynamics onto the GLE [18,31]: (8) where M is the particle mass, V(t) is its velocity, K(t−s) is the memory kernel, Γ(t) is the stochastic force, and we allow for the possibility that U eff Harm = 1 2 kR 2 may differ from U Harm , as is known also for non-linear systems [30,87].Here and throughout, we will use the stationary GLE [58], i.e., the initial time T in Eq. ( 1) is taken to be in the infinite past.Following a strategy outlined in Ref. [58] (see also appendix C), we use Eq.(8) to construct a Volterra equation of the second kind, which relates the memory kernel with the force autocorrelation function C F (t) = ⟨F(t)•F(0)⟩ (with F = M V), the velocity autocorrelation function (VACF) C V (t) = ⟨V(t)• V(0)⟩, and the cross correlation Eq. ( 9) can then be used to determine K(t) following a modification of an algorithm presented in Ref. [58].Details are shown in appendix C.
The corresponding 2FDT for GLE (8) in the presence of a harmonic potential is derived in Appendix A: with I(t) = ∞ t ds K(s).It differs from Eq. (2) by an additional term which is proportional to the instantaneous correlations between the position and the velocity.This term is zero at equilibrium due to time reversal symmetry, but could, in principle, be relevant in stationary active systems.In the systems considered in the present work, however, ⟨VR⟩ was always found to be zero within the statistical error (data not shown).We infer that the additional term can be neglected and the coarse-grained GLE for our systems satisfies the usual generalized 2FDT, Eq. (2).

Position probability distribution
In equilibrium, a harmonically trapped particle has a Gaussian position distribution centered at the trap center with variance ⟨(R i − R 0,i ) 2 ⟩ = k B T /k for each axis i.As expected, this matches what we find for a probe immersed in a passive bath, as we can see in Fig. 1a), which shows the position distribution of the probe in the x-direction relative to the trap center.Although we only show the position distribution in the x-direction, since the probe is confined along all axes, we find the same form of position distribution in both the y-and z-directions as well.
We see in Fig. 1b), that the position distribution of the probe remains Gaussian and centered around the trap center even in an active bath, as in equilibrium.However, the variance of the probe position distribution is higher than that of a probe immersed in a passive bath.These results are consistent with previous studies of a harmonically confined probe immersed in an active bath [88][89][90].We note that, for higher values of k when the range of confinement is comparable to the persistence length of the active particles, the distribution has been shown to become non-Gaussian [90].However, all our values of k remain within this limit.
These previous studies have attributed the enhanced variance of the position distribution to an enhanced temperature.However, when we use the definition of the enhanced kinetic temperature as quoted earlier from Ref. [65], we find that k B T eff /⟨x 2 ⟩ ̸ = k, as can be seen in Fig. 2 where we plot k B T eff /⟨x 2 ⟩ as a function of k.Furthermore, we see in Fig. 2 that the value of k B T eff /⟨x 2 ⟩ for a given constraint k is not density independent in the active case, as we would expect in equilibrium.
Nevertheless, it does appear that k B T eff /⟨x 2 ⟩ ∝ k.We therefore posit that, due to the activity of the bath particles, the external force on the probe needs to be renormalized in our coarse-grained GLE model.Then, k B T eff /⟨x 2 ⟩ = k for the active bath, where k = αk is the renormalized spring constant of the harmonic potential with renormalization factor α.
To determine the renormalization factor α, we fit our data in Fig. 2 to a line with zero offset.As expected, in the passive case, α ≈ 1 for both bath densities, i.e., the external force does not need to be renormalized and is density independent.In the active case, we find that for a bath with density ρ = 0.4σ −3 , α 0.4 = 0.190 ± 0.001 and for a bath with density ρ = 0.8σ −3 , α 0.8 = 0.079 ± 0.002.In both cases, α < 1, meaning that the activity of the bath effectively decreases the spring constant of the harmonic potential and, thereby, decreases the trapping force exerted on the probe.In Appendix B, we show how α varies as a function of F0 .We find that α decreases from the equilibrium value α = 1 as a function of F0 .Furthermore, we find that a renormalization factor is necessary even for very small bath activities ( F0 = 1).

Velocity autocorrelation function and memory kernel
We now examine the velocity autocorrelation function (VACF) and memory kernel of the trapped probe.In calculating the memory kernel, we use the renormalized spring constant k obtained from Section III A 1. Interestingly, we find that the reconstructed memory kernel fails to go to zero in the long time limit if we insert the bare spring constant or omit the spring potential altogether.
In Appendix A, we analyze this problem more generally.We show that, for any confining renormalized external potential U eff , inverting the Volterra equation to extract the corresponding memory kernel is bound to fail unless the following (tensorial) condition is fulfilled: Eq. ( 11) can be seen as a generalized version of the equipartition theorem in nonequilibrium systems.If condition ( 11) is not met, the reconstructed memory kernel will not converge to a zero value in the limit t → ∞.Equation ( 11) is another way to characterize the renormalized potential.From Figs. 3a) and c), we see that, when immersed in a passive bath, the VACF/memory kernel of a harmonically constrained probe only differs slightly from that of an unconstrained probe.Namely, the oscillation of the VACF for a harmonically constrained probe is slightly deeper than that for an unconstrained probe.The increased depth of this oscillation relates to the restoring force on the probe; therefore, we expect that its depth would increase with an increased value of k, which agrees with what we see in Fig. 3a).The memory kernel of a free vs. trapped probe in a passive bath exhibits even fewer differences.In fact, the differences between the memory kernel of a free vs. trapped particle all fall within our error bounds.
When immersed in an active bath, the VACF of the harmonically constrained probe depends more significantly on the confinement strength than in the passive case.In the inlay of Fig. 3b), it can be observed that the trapped probe has a decreased local kinetic temperature, C V (0), in comparison with that of a free probe.Furthermore, the local kinetic temperature decreases as the trapping strength k increases.Within the range of k values we studied, however, this effect is relatively small (two orders of magnitude smaller than that of the kinetic temperature) and cannot be readily seen in the probe velocity distributions (see Appendix D).Therefore, this change in the kinetic temperature cannot account for the necessary rescaling of k.The memory also slightly depends on k at small times, but the differences are again smaller than the error at times t > 0.1.

Stochastic force distribution
Once we have determined the memory kernel of the probe, it is straightforward to determine the stochastic force Γ(t): (12) In Fig. 4, we see that in both the case of a passive and an active bath, the stochastic force distribution calculated from simulation data (solid lines) approximately matches a zero-centered Gaussian with the same standard deviation (dashed lines).The only deviations from the zero-centered Gaussian are an enhanced peak at the center of the distribution, which we infer results from the low density of the bath and which was also seen for an unconstrained probe in Ref. [65].

Characterization of active bath
We now examine the properties of the ALP fluid in the vicinity of the probe.To this end, we analyze the angle dependent density distribution ρ(r) of ALPs around the probe in a comoving and corotating frame with origin at the probe center, R(t), and z axis always aligned in the direction of the instantantenous probe velocity V(t).To quantify the angular dependence, we expand the density distribution ρ(r) in spherical harmonics,  15)) and b) dipole moment (Ω 0 1 (r), see Eq. ( 14)), surrounding a free probe immersed in a bath with activities F0 = 0.0 (indigo) and F0 = 5.3 (orange).The probe has a radius Rp = 3.0σ and is immersed in a bath with global density ρ = 0.4σ −3 .The black dashed lines mark the global density in a) and the zero baseline in b).
The coefficients Ω m l (r) can be determined from simulations according to where the sum n∈δr runs over all bath particles in a spherical shell δr around the particle (i.e., their distance from the probe center lies within the interval [r − δr/2 : r + δr/2]), and V is the volume of the shell, V(r) = 4/3π((r + δr/2) 3 − (r − δr/2) 3 ).Specifically, the radial average of ρ(r) can be obtained from We first characterize the bath in the case that the probe is free.In comparing the radially averaged density curves for a passive and an active bath surrounding a free probe in Fig. 5a), we see that the first peak of the density curve is higher in an active bath.Hence adding activity to the bath leads to a higher density of bath particles in the vicinity of the probe.We also see that the first shell of active particles is shifted slightly closer to the probe center in comparison with that of the passive bath.Furthermore, the spacing between the first and second shells of the active bath particles are closer together.This -in combination with the forward shifted first shell -leads to the coincidence of a peak in the density profile of the active bath with a trough in the density profile of the active bath at r ≈ 4.9σ.These closer shells indicate that ALPs are able to move closer to the probe, and each other, than passive bath particles.We infer that this ability stems from the higher kinetic energy of ALPs due to their active force, which allows them to overcome more of the repulsive potential from WCA interactions with the probe.
In a passive bath, all higher order moments of the bath are zero.However, in an active bath, the active fluid acquires a negative dipole moment surrounding the probe, which is sustained to large values of r, as can be seen in Fig. 5b).This negative dipole moment indicates that ALPs collect behind the probe relative to its instantaneous velocity, even at large distances from the probe.The phenomenon of active particles gathering behind an immersed passive probe has previously been seen for a probe dragged through an ABP bath in Ref. [20], where it is framed as a difference between the forces behind and in front of the probe.We find no significant structure in higher order spherical harmonics for either the active or passive bath.
Returning to the system of a probe immersed in an active bath and subject to a harmonic constraint, we find that in spite of the additional harmonic force on the probe, the density profile of the ALP bath (as calculated from Eq. ( 15)) remains identical for a trapped and a free probe (see Fig. 6a)).Additionally, its dipole moment remains the same whether the probe is trapped or free (see Fig. 6b)).This is most likely due to the fact that the forces exerted by the harmonic trap are relatively small, i.e., within the linear regime.We will see in section III B that the situation changes if we subject the probe to very large forces.

B. Constant drag force
Next we consider a passive probe which is pulled through a bath of interacting ALPs with a constant 'drag' force along the x-axis: F = F ext e x (see Fig. 7).Thus, the equation of motion for the probe is: where M is the mass of the probe, V(t) is its velocity, and U WCA (R − r n ) is the WCA potential due to an ALP particle, n, at position r n .Aside from the constant external force, the immersed probe only experiences forces from interactions with the surrounding ALPs.The ALP bath has a global density of 0.4σ −3 .In the steady mode, the constant drag force induces a steady drift of the probe in x direction with average velocity ⟨V⟩ neq ≡ V.
In our coarse-graining approach we map this system onto a stationary GLE of the form where f is the effective drift force.We use a general tensorial memory kernel here as K may become anisotropic for large drag forces [17].Following Ref. [17], we split this equation into two parts: Here, V(t) = U(t) + V such that V is the drift velocity of the probe.Eq. ( 19) has exactly the same form as the GLE for an unconstrained probe.The only difference is that the velocity U(t) is now the velocity of the probe in a co-moving frame with velocity V. Since the equations are of the same form, Eq. ( 19) can also be transformed into a Volterra equation, which can then be inverted to find the memory kernel of the probe.The generalized 2FDT is fulfilled by construction [17].Using such a mapping, the effective drift force on the probe can then be calculated from Eq. 18, as f = At equilibrium, the behavior of the probe at small drag velocities can be described by linear response theory.One purpose of the present study is to examine whether a similar behavior can be found in an active bath, if we renormalize the force according to the prescription found in the previous section, III A. In linear response theory, the drift velocity of the probe in response to a constant drag force F ext is given by the linear relation ⟨V⟩ = µ eq F ext with the mobility µ eq = ∞ 0 K eq (s)ds −1 .Here K eq is the equilibrium memory kernel evaluated at F ext = 0.The analogous behavior in an active bath would be ⟨V⟩ = µ 0 F eff ext , where µ 0 is determined from the integrated memory kernel of the free probe and F eff ext = αF ext is renormalized according to Section III A. We should note that that linear response theory is not generally valid in our system [65].Our goal here is to probe the existence of a linear regime where ⟨V⟩ ∝ F ext , and test the relation ⟨V ⟩/F ext = αµ 0 in that regime.
In Ref. [20], a probe dragged through a 2D bath of ALPs has been studied for drag forces well beyond the linear regime.Here, we focus on small to intermediate forces and possible transitions between linear and nonlinear regimes.

Static mobility
In the following, we use dimensionless quantities Fext = 1 γ m kBT F ext and μ = µ 1 γ .Furthermore, to unify the notation for active and passive bath particles, we define a dimensionless renormalized force Fext , which is given by Fext = Fext in a passive bath and Fext = α Fext in an active bath.We will consider three values of the mobility (in dimensionless form): (i) the drag mobility μDrag defined by where ⟨ Vx ⟩ = V x is the dimensionless average velocity along the x-axis; (ii) the linear GLE mobility μ0 defined by where K0 is the (isotropic) dimensionless memory kernel of a free probe (at Fext = 0); and (iii) the tensorial GLE mobility μK defined by where K is the memory kernel in the presence of the external force.To properly compare the ranges of the linear regime in the passive and active bath case, we introduce a dimensionless Peclet number for the probe: Pe = ⟨ Vx ⟩/ Vdiff , where the diffusion velocity scale Vdiff is given by Vdiff = D eff m/k B T /R p and depends on the probe radius R p = 3.0σ and the effective diffusion coefficient of an isolated ALP, as calculated in Ref. [65].In terms of Fext , the Peclet number is given by Pe = μDrag Fext / Vdiff (see Appendix E).
In Fig. 8a), we graph ⟨ Vx ⟩ as a function of Fext to determine the drag mobility μDrag within the linear regime using Eq.(20).We note that, for low values of F ext in an active bath, the stochastic motion of the active bath particles dominates the applied external force.Hence, the deviation of ⟨ Vx ⟩ is significantly larger for low values of F ext , resulting in larger errors in our data.
We extract a mobility μDrag = 0.02057 ± 0.00006 for a probe immersed in a passive bath ( F0 = 0.0) with a density ρ = 0.4σ −3 .For a probe in an active bath ( F0 = 5.3) of density ρ = 0.4σ −3 , we extract a mobility of μDrag = 0.154 ± 0.006.Comparing these values, we see that increasing the activity of the bath also increases the mobility of the probe particle.
On the other hand, using the memory kernels of a free probe and Eq. ( 21), we obtain the linear GLE mobility μ0 = 0.0226±0.0003for a probe in a passive bath of density ρ = 0.4σ −3 and μ0 = 0.147 ± 0.002 for a probe in an active bath of density ρ = 0.4σ −3 with F0 = 5.3.These values for the mobility agree well with those calculated from Eq. (20).The slight discrepancy between values is most likely due to the reconstruction of the memory, which becomes less accurate at longer times.We emphasize that renormalizing the drift force in Eq. ( 20) is essential for reaching this agreement.
For high forces, the drag mobility is no longer a fixed value.In Fig. 8b), we show the force dependent drag mobility (μ Drag ( Fext ) = ⟨ Vx ⟩/ Fext ) as a function of the Peclet number.We expect μDrag ( Fext ) to be constant within the linear response regime for a probe immersed in a passive bath.The figure shows that this is indeed the case for Peclet number up to Pe ∼ 1.For probes immersed in an active bath, we observe a very similar behavior, i.e., a crossover from roughly constant drag mobility to force dependent drag mobility around Pe ∼ 1. Identifying this crossover is more difficult in this case because the measured values for μDrag fluctuate strongly for low values of Pe.These fluctuations reflect the fact that the activity of the bath induces larger variations in the probe velocity.
We should note that, due to the enhanced diffusion of the active bath particles in comparison to the passive bath particles, equivalent Peclet numbers for a probe in an active and a passive bath correspond to different external drag forces (see also Appendix E).Much larger drag forces are required to achieve the same Peclet number for an active bath in comparison with a passive bath.
Figs. 8c) and d) show that the behavior beyond the linear response regime is qualitatively different for a probe immersed in an active vs. a passive bath.A probe dragged through a passive bath exhibits thinning be-  havior (increased mobility), whereas a probe dragged through a bath of ALPs exhibits thickening behavior (reduced mobility).Due to these opposite behaviors, the values of μDrag ( Fext ) for a passive and an active bath at the same Peclet number approach each other beyond the linear regime.We will discuss and analyze this further at the end of Section III B 5.

Velocity autocorrelation function and memory kernel
Having looked at the active microrheological properties of the active bath, we turn to the dynamic properties of the dragged probe: namely, the VACF and the memory kernel.We map the dynamics of our system onto Eq. (17).We examine the VACF in the co-moving frame because, for a probe dragged through a passive bath, it should be roughly identical for all drag forces within the linear response regime [17].We would like to evaluate the universality of the co-moving VACF for a probe dragged through an active bath as well.
As we would expect, within the linear regime (Pe < 1), the VACF and memory kernel are isotropic and independent of Fext for a probe dragged through a passive bath.We note that, in Figs.10a) and b), the memory kernel for a probe dragged with Pe ≈ 1.5 through a passive bath already deviates slightly from the memory kernel of the undragged probe.This suggests that the transition from the linear to the nonlinear regime in a passive bath is a relatively sharp transition at Pe = 1.Because these deviations are primarily visible in the memory kernel as opposed to the VACF, this suggests that the transition to the nonlinear regime is primarily due to the stochastic force distribution.In fact, for a probe dragged through a passive bath at Pe ≈ 1.5, we can already see a slight asymmetry in the probe's stochastic force distribution (see Appendix F).This asymmetry is a signature of the nonlinear regime, as will be discussed in the next section, 11.
For the active bath, we also see that the VACF and memory kernel are isotropic and independent of Fext within the linear regime.Furthermore, in this case, the linear regime appears to hold for Peclet numbers of up to at least Pe ≈ 1.5 (Figs. 9 and 10), suggesting that the tranistion from the linear to nonlinear regime is less sharp in an active bath.
Beyond the linear regime, however, the VACF and memory kernel are neither isotropic (see Figs. 9e,f) and 10e,f)) nor independent of Fext for both passive and active systems.In both cases, a large drag force causes the VACF to decay more rapidly in both its parallel and perpendicular components.The memory kernel is primarily affected at short times, where it is greater than the memory kernel of a free probe.This is, again, true for both a probe dragged through a passive and an active bath.However, the other aspects of the VACF/memory kernel form differ for passive and active baths.In a passive bath, we see that the dragged probe acquires an increased local kinetic temperature in both its parallel and perpendicular components.In an active bath, on the other hand, the dragged probe only acquires an increased local kinetic temperature in its parallel component.In its perpendicular component, the probe dragged through an active bath actually experiences a decrease in local kinetic temperature.
As a consequence of the changes in the memory ker- nel, the GLE mobility becomes anisotropic at high Peclet numbers, and both its parallel and perpendicular components increase (data not shown).In the linear regime, however, the tensorial and linear GLE mobility are identical: μK = 1μ 0 for Pe < 1, where 1 is the unit tensor.

Stochastic force distribution
Once we have determined the memory kernel of the probe, we can determine the stochastic force through a trivial re-ordering of Eq. ( 19) [58]: where F(t) is the instantaneous force acting on the particle, which is calculated during the simulated trajectory.The probability distribution of the stochastic force is shown in Fig. 11.In both the case of a passive and an active bath, for low Peclet numbers (within the linear regime), the stochastic force distribution calculated from simulation data (solid lines) approximately matches a zero-centered Gaussian with the same standard deviation (dashed lines).The only deviations from the zerocentered Gaussian are, again, an enhanced peak in the distribution, particularly in the case of the passive bath, which we infer results from the low density of the bath.However, at higher drag forces, the distribution calculated from simulation data develops an asymmetry in both components parallel and perpendicular to the drag force.This asymmetry has already been seen for a pas- sive bath in Ref. [17].Our results here show that it also occurs for an active bath, although the magnitude of the asymmetry appears to be lower in the active case.We infer that the lower asymmetry in the active case stems from the already enhanced fluctuations due to the bath activity.
The asymmetry takes the form of a long tail in the negative forces of the distribution.This long tail occurs because, outside of the linear regime, bath particles in the vicinity of the probe do not necessarily travel at the same relative velocity as the probe [91].This leads some bath particles to 'crash,' with relatively high forces, into the probe, promoting the long tail.Because the average of the stochastic force distribution remains zero, the distribution also exhibits a slightly enhanced probability for positive forces.The distribution outside the linear regime can be described by a split normal Gaussian [17].15)) of the immersed probe in baths with activities F0 = 0.0 (indigo) and F0 = 5.3 (orange).Open circles with dotted lines show ρ(r) for a probe which is not dragged, whereas solid lines with stars show ρ(r) for a probe dragged with Pe ≈ 7.2.The probe has a radius Rp = 3.0σ.The black dashed line marks the global density, ρ = 0.4σ −3 .

Characterization of the active bath
We hypothesize that the behavior of the probe mobility beyond the linear response regime results primarily from the different number of contacts experienced by a dragged probe.Fig. 12 shows that both for a passive and an active bath, the density profiles (calculated from Eq. ( 15)) for a free vs. a dragged probe particle differ from each other.However, the effect of the drag force on the density profile is opposite in the two cases: Whereas a probe dragged through a passive bath has a higher number of contacts in its immediate vicinity, a probe dragged through an active bath has a lower number of contacts.
It has previously been shown that a probe dragged with sufficient force through a passive bath experiences an increased density of particles in front of it [17,92].By examining the dipole moment in the vicinity of the probe (see Fig. 13, calculated from Eq. ( 14)), we verify this buildup of particles in front of a probe dragged through a passive bath with Pe > 1.A similar frontal buildup is exhibited by a probe dragged through an active bath with Pe ≈ 7.2.However, for lower Peclet numbers (Pe ≈ 1.5 and Pe ≈ 4.5), it is no longer observed.In spite of the drag force, we still see an accumulation of ALPs behind the probe for Pe < ∼ 1.5.We reason this qualitative difference in the bath behavior at the probe interface following the arguments of Ref. [20].For small Peclet numbers, the speed of the ALPs is greater than that of the dragged probe.Thus, the ALPs are able to accumulate behind the probe and push, as was the case for a probe without an external drag force.Consequently, we still see a negative dipole moment in the active bath for Pe ≈ 1.5 in Fig. 13a).For intermediate Peclet numbers, the dragged probe and the ALPs move at similar speeds.Therefore, the dipole moment of the active bath approaches zero, which we see happening for Pe ≈ 4.5 in Fig. 13b).Although for this value of Pe the dipole moment is not entirely zero, we can see that its magnitude and range are smaller.For large Peclet numbers, ALPs can no longer catch up to the probe and, consequently, the active bath acquires a positive dipole moment in the vicinity of the probe, as seen for Pe ≈ 7.2 in Fig. 13c).

Effective drift force at high drag forces
To further characterize the transition from the linear to the nonlinear regime, we now evaluate the dimensionless effective drift force in our coarse-grained GLE, Here, K∥ (t) is the parallel component of the dragged probe's memory kernel shown in Figs.10a) and b) and we have exploited the fact that K and hence μ are diagonal.In the linear regime, we have already established μDrag = μK,∥ = μ0 in Secs III B 1 and III B 2, hence fx = Fext .Our goal in this section is to test this relation in the nonlinear regime.To this end, we plot C fx versus Pe = C Fext in Fig. 14.Here we rescale by C = μDrag / Vdiff (see Appendix E) in order to be able to compare the nonlinear effects in the passive and active bath case.
The figure shows that the relation fx = Fext holds in the linear regime as expected, but breaks down at high Peclet numbers.In the case of the a passive bath, fx is smaller than Fext = Fext at large drag forces.In the case of an active bath, fx is larger than Fext = α Fext .As a result, the values of fx in the active and passive bath at given physical drag force Fext approach each other.This is shown in the inset of Fig. 14, where the ratio fx / Fext is plotted as a function of the high Peclet number.We conclude that the change of the drag mobility μDrag = ⟨ Vx ⟩/F ext in the nonlinear regime, which has been reported in Fig. 8, can be decomposed into two factors: First the GLE mobility μK decreases with increasing Peclet number, both in a passive and in an active bath (cf.Fig. 10).Second, the force renormalization factor fx / Fext increases in an active bath, and decreases in a passive bath.In the nonlinear regime, F ext thus must be renormalized even in the passive system.In a passive bath, both factors result in thinning behavior.In an active bath, they compete with each other, but the force renormalization effect dominates; hence one effectively observes thickening.The latter can be associated with the redistribution of bath particles around the probe, from back end to front end, shown in Fig. 13.
In the literature, both 'thinning' (increased mobility) [20,24,27,[92][93][94][95][96] and 'thickening' (decreased mobility) [17,92,95,97,98] behavior have been observed for different systems beyond the linear regime.In general, systems which are not dominated by hydrodynamics exhibit thinning behavior [95].This aligns with our finding that a probe dragged through a passive bath exhibits thinning behavior.However, in spite of the fact that our system neglects hydrodynamics, we find that a probe dragged through a bath of ALPs exhibits thickening behavior.This thickening behavior is in contrast to what has been previously observed in an ABP bath, where particle inertia is also neglected [20,24,27].Thickening has, however, been observed in microrheological studies of granular active systems, where inertia plays an impor-tant role [12,98].

IV. DISCUSSION AND OUTLOOK
In sum, we have analyzed the dynamics of a probe particle immersed in a bath of active Langevin particles in the presence of external forces by molecular dynamics simulations.We found that it is possible to map this system onto a GLE that satisfies a 2FDT.However, in order to achieve full consistency, we must (i) rescale the temperature such that k B T eff = M ⟨V 2 ⟩/d, with d being the dimensionality of the system.More generally, we must replace the equilibrium 2FDT by a generalized 2FDT, Eq. ( 2).
(ii) renormalize the external force, F eff ext = αF ext .It must be chosen such that the distribution of probe positions in a confining renormalized potential U eff (R) satisfies a generalized equipartition theorem, Eq. (11).We have identified a linear regime of small to moderate forces where the renormalization factor α does not depend on the amplitude of the force.
If these two conditions are fulfilled, the positions of the probe in a confining potential are Boltzmann distributed according to a Boltzmann-like distribution, P (R) ∝ exp(−U eff /k B T eff ), the mobility of the particle can be determined consistently as semi-infinite integral over the memory kernel, and the 2FDT ensures that the mobility µ and long-time diffusion constant D of a free particle are connected by a generalized Einstein relation, D = µ k B T eff (see Eq. (A9) in the appendix).
The ability to map this non-equilibrium system onto a GLE that satisfies the generalized 2FDT enables the use of many recently developed coarse-grained methods which rely on its fulfillment [61,99,100].
We emphasize that both temperature and force renormalization are necessary to achieve full consistency.This is because the effective temperature and the force renormalization depend on each other via the generalized equipartition theorem.On less formal grounds, we can also argue as follows: If we only rescale the temperature, as is commonly done, the probe particle will have the correct velocity distribution in the GLE model, but the position distribution will deviate from the true distribution.If we only renormalize the potential, the probe will not have the correct velocity distribution.
It is instructive to discuss our findings in relation to another popular type of coarse-grained model for particles immersed in an active fluid, i.e., the overdamped Markovian Langevin equation (OLE) [27].An OLE for such a system should reproduce three target quantities: (i) The diffusion constant of a free probe, (ii) the drift velocity in the presence of a slowly varying drag force, (iii) the position distribution in the presence of a confining potential.Since the OLE can be derived from a GLE in certain limits, we can use our results to propose a suitable OLE: Here Ξ(t) is a Gaussian distributed white noise satisfying ⟨Ξ(t)Ξ(t ′ )⟩ = 2q1δ(t − t ′ ) with q = k B T eff µ K = k B T μ and we have defined the OLE mobility μ = µ K α and the OLE temperature k B T = k B T eff /α.This analysis shows that, in the overdamped limit, it is not possible to distinguish between a renormalization of the mobility and a renormalization of the force, as long as the renormalized temperature is adjusted accordingly.Physically, this means that the separate effect of force renormalization -as opposed to just temperature renormalization [101] -only becomes relevant if inertial effects are important.As such, our results are most experimentally relevant for systems in which the constituents are of the meso-or macroscopic scale, such as active granular matter systems [8][9][10][11][12][69][70][71][72] or systems of synthetic active robots [66][67][68].However, in principle force renormalization can be distinguished from temperature renormalization for any system in which you can independently measure the kinetic temperature.
The coarse-grained GLE for a single probe particle in a renormalized potential is equivalent to a GLE for an equilibrium system, therefore we cannot expect to find interesting nonequilibrium effects in such a model.However, the renormalization of temperature and forces has interesting consequences for systems of several interacting particles.We found that the kinetic temperature of the probe particle not only depends on the parameters of the bath such as density and activity level, but also on the properties of the probe, such as the size [102].A similar behavior can also be expected for the force renormalization factor α. If probes of different size interact with each other, we therefore expect the forces to be rescaled differently, implying that the effective forces are non-reciprocal.This is a common signature of nonequilibrium [103][104][105][106][107][108][109] whose consequences for our particular type of system yet remain to be explored.
Furthermore, we have seen that the GLE mapping with constant renormalization parameter α breaks down for very large external forces.If α depends on the force, the renormalization may turn a conservative force into a non-conservative drift force, i.e., it may no longer be possible to rewrite the renormalized force as a gradient of a renormalized potential.Developing coarse-grained methods for external forces beyond the linear response regime will be an important area for future research.Moreover, theoretical arguments (appendix A) suggest that it might sometimes be necessary to modify the 2FDT by adding another force-dependent term, e.g., as in Eq. (10).In the case studies considered in the present work, the additional term was found to vanish; in general, it might be nonzero and generate true nonequilibrium signatures even in a single-particle model.
It has been shown that a harmonically confined active particle exhibits different behavior from a harmon-ically confined passive particle [110][111][112].Namely, the radial probability distribution of the active particle exhibits two 'phases' which depend on the particle activity and the trap strength: a 'passive' phase described by a Boltzmann-like distribution around the trap center, and an 'active' phase described by a non-Boltzmann distribution which is peaked away from the trap center [112].Given the shared behaviors of a probe immersed in an active bath with an active particle itself [82,113,114], it would be interesting to determine whether there exists an 'active' phase probability distribution for a probe immersed in an active bath for other values of activity and spring potentials, and whether this can be described by a a GLE with a modified 2FDT and/or a force-dependent α-parameter.
where we have used standard relations for the Fourier transform and C V (0) = ⟨VV⟩.
(II) We first discuss the case of a free particle with U eff ≡ 0. In that case, the Volterra equation in the limit ω → 0 reads M ⟨VV⟩ = K(ω = 0) C V (ω = 0) and hence we obtain a generalized Einstein relation between the mobility and diffusion tensors, In isotropic systems, one has M ⟨VV⟩ = 1 k B T eff , D = 1 D, and µ K = 1 µ K , and the generalized Einstein relation simplifies to the more common form D = k B T eff µ K .From Eq. (A8), we can deduce that K(0) is not bounded if CV (0) = 0.In this case, both the diffusion constant and the mobility are zero.
(III) Next we consider a confined particle subject to a confining potential U (R), such that ⟨R 2 ⟩ < ∞, e.g., a particle in an external potential with infinite barrier height.It does not diffuse, hence C V (ω = 0) = 0. Since K(ω = 0) is finite (one of our basic assumptions) and C V (ω = 0) = 0, the memory term in Eq. (A8) vanishes, and so does the term −iωM C V on the left hand side of the equation.Thus the Volterra equation only has a solution if M ⟨VV⟩ = − lim ω→0 C fV (ω), which implies which is listed as Eq. ( 11) in the main text.We note that this condition is equivalent to requesting θ R (0) = 0.This can be seen by multiplying Eq. (A1) with R(t), inserting f = −∇U eff , taking the trajectory ensemble average, and using 0 = d dt ⟨RV⟩ = ⟨R V⟩ + ⟨VV⟩.Applying the general result (11) to the harmonic oscillator potential, U eff (R) = 1 2 k R 2 , we find that the effective spring constant must fulfill k = M ⟨V 2 ⟩/⟨R 2 ⟩ = k B T eff /⟨x 2 ⟩, which corresponds to the renormalization of the spring constant in the main text.
(IV) Finally, we analyze the implications of (A6) for the fluctuations of the stochastic noise Γ(t).We have analyzed the force free case in Ref. [17] and shown that the Volterra equation then enforces a generalized 2FDT (Eq. ( 2) in the main text).Here we use Eqs.(A1) and (A5) as a starting point and follow the same steps as Ref. [17] (Appendix A) to derive a corresponding equation in the presence of a drift force.First, We write Eq. (A1) in terms of a time variable t 0 and then multiply this equation by Γ(t).We then take the ensemble average to obtain that Here, f (t) stands for f (t) = f (R(t)), and we use the Einstein summation convention.For the remainder of the derivation, we will indicate the insertion of previous results with brackets ... .With the definition of θ V (Eq.(A4)), we can express the first term on the right hand side (r.h.s.) of (A12) as We can additionally re-express the last term on the r.h.s. of Eq. (A12) using the definition of θ V given in Eq. (A4): we rewrite the second term on the r.h.s. of (A14) as In the last step, we have inserted Eq. (A5) and used the symmetry relation for C V once more.Combining (A12)-(A15), we obtain the following general tensorial expression for the noise fluctuations: The last three terms vanish because θ V ≡ 0. However, the generalized 2FDT, Eq. (A11), is only recovered if one has, in addition, θ f (t − t 0 ) ≡ 0 (noise orthogonality with respect to the force).As discussed in Ref. [17], this condition can be enforced if the force is constant or only depends explicitly on time, as in the example of the constant drag force.On the other hand, it may be broken if the force depends on the position of the particle [43,46].In such cases, the 2FDT is modified according to We now specifically consider the case of the (possibly renormalized) harmonic potential, where the drift force can be expressed as f = − kR.Then we have θ f (t − t 0 ) harm.= − kθ R (t − t 0 ).θ R (t − t 0 ) can be calculated explicitly as follows: We first take the derivative of θ R (t − t 0 ) with respect to t 0 , and obtain where we have used θ V (t − t 0 ) ≡ 0. Integrating over t 0 and using θ R (t − t 0 ) → 0 at (t − t 0 ) → ∞ gives an expression for θ R , which is generally valid: Inserting this result into Eq.(A17), we finally obtain the following expression for the modified 2FDT in the presence of a harmonic potential: At equilibrium, we have ⟨VR⟩ = 0 due to time reversal symmetry.Therefore, the last term vanishes and one recovers the force-free generalized 2FDT, Eq. ( 2) in the main text.In active systems, this is not necessarily the case.For example, active particles have a tendency to accumulate at boundaries in a state where the active velocity points towards the boundary, which creates a clear correlation between position and velocity.If such a correlation were present for a harmonically trapped particle, omitting the last term in Eq. (A20) in the corresponding coarse-grained GLE would effectively map the system onto an equilibrium system, where ⟨VR⟩ = 0 by construction.In such cases, we can only recover ⟨VR⟩ ̸ = 0 by using the full modified 2FDT.We note that it may be possible to define projection operators which formally lead to other definitions of the memory kernel and thus other FDTs in which the ⟨VR⟩ term may or may not be present.Here, however, we do not apply projection operator techniques, but rather analyze the GLE under the sole assumption that the memory is given by the Volterra equation.In this case, provided that the system is not at equilibrium, the ⟨VR⟩ term cannot generally be eliminated.
Numerically, however, we found the additional term in the 2FDT to be negligible in the systems considered in the present work.The value of ⟨VR⟩ was always zero within the error.We note that, due to Eq. (A19), this also implies θ R (t − t 0 ) ≡ 0 (noise orthogonality with respect to the position).On the other hand, the force renormalization in an active bath was substantial, as reported in the main paper, and moreover, the correct implementation of the force renormalization was crucial for the memory reconstruction, see appendix C. Setting t 0 = 0, this gives the Volterra equation of the second kind quoted in the main text, Eq. (C2) is inverted using [58]: with the initial condition K(0) = 1 M C F (0) C V (0) − k.In Section III A 2, we used the the renormalized harmonic constraint k in Eq.C3 to solve for the probe memory kernel based on our findings in Section III A 1. Additionally, we tried solving the memory kernel without a harmonic constraint ( KNH ( t)) as well as with the imposed harmonic constraint k ( Kk ( t)).The results are shown in Figs.16a) and b) respectively.Although Figs. 16a) and b) both look qualitatively the same as Fig. 3d), comparing the Figs.16a ′ ) and b ′ ) with Fig. 3d ′ ) reveals that they are not the same.In Fig. 3d ′ ), where we define K( t) with k, the difference between the memory kernel of the constrained probe and that of a free probe goes to zero.However, this difference converges to a finite value at late times in the case that we either neglect the harmonic constraint (Fig. 16a ′ )) or we define K( t) with k (Fig. 16b ′ )).This finding is similar to the renormalization/linearization within the Mori-Zwanzig formalism in the presence of non-linear external potentials [30,87].We know that the memory kernel of the free probe goes to zero in the long time limit [65].Therefore, the fact that KNH ( t) and ∆ Kk ( t) go to finite values, indicates that they do not approach zero in the long time limit as they should.Comparing Fig. 16a ′ ) to that of Fig. 16b ′ ), we see that KNH ( t) approaches a smaller (in terms of magnitude) finite value than ∆ Kk ( t).This can be understood given that, in our case, the value of k is closer to zero than to k.We conclude that consistent force renormalization

FIG. 1 .
FIG. 1. Position distribution in x-direction (P (x)) relative to the trap center of a harmonically trapped probe immersed in a) a passive ( F0 = 0.0) and b) an active bath ( F0 = 5.3).Solid lines show simulation data, whereas dashed lines show a Gaussian distribution with mean 0 and the variance a) expected for a probe in a harmonic potential with each respective spring constant k and b) of the simulation data.Each bath has average density ρ = 0.4σ −3 .Values of k are shown in the legend.

3 FIG. 2 .
FIG.2.kBT eff /⟨x 2 ⟩ as a function of k for a probe immersed in passive baths ( F0 = 0.0, indigo) and active baths ( F0 = 5.3, orange) with density ρ = 0.4σ −3 (triangles) and ρ = 0.8σ −3 (x-symbol).Error bars on the data points are plotted, but are too small to be visible.Solid colored lines show the line of best fit for baths with density ρ = 0.4σ −3 , whereas dashed colored lines show the line of best fit for baths with density ρ = 0.8σ −3 .The plot also contains a solid black line showing kBT eff /⟨x 2 ⟩ = k.

FIG. 3 .
FIG. 3. Velocity autocorrelation function (VACF) and memory kernel of an immersed probe trapped in a harmonic potential.The larger graphs in the left column show a) the VACF and c) the memory kernel for a probe immersed in a passive bath ( F0 = 0.0), whereas the larger graphs in the right column show the b) VACF and d) memory kernel for a probe immersed in an active bath ( F0 = 5.3).Light shading indicates error bars.Each bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.Smaller graphs correspond to the larger graph directly above them and show the difference between the VACF/memory kernel of a free and trapped probe.

FIG. 4 .
FIG. 4. Stochastic force distribution of the harmonically trapped probe in a bath of density ρ = 0.4σ −3 .The solid lines show simulation data, whereas the dotted lines show zero-centered Gaussian distributions with the same standard deviation.a) The stochastic force distribution for a passive bath.b) The stochastic force distribution for an active bath with activity F0 = 5.3.In b), the Gaussian distributions overlap the simulation data.

2 FIG. 6 .
FIG.6.a) Radially averaged density (Eq.(15))) and b) dipole moment (Ω 0 1 (r), see Eq. (14)), surrounding a harmonically confined probe immersed in a bath with global density ρ = 0.4σ −3 and activity level F0 = 0.0 (indigo) and F0 = 5.3 (orange).The probe has a radius Rp = 3.0σ.Different symbols show different harmonic trapping strengths k = 1ϵσ −3 (plus symbol), k = 5ϵσ −3 (star symbol), k = 10ϵσ −3 (x symbol).The results for a free probe are shown in a) by open circles and in b) by bars.For ease of visualization, in a) we add different vertical shifts to the density profiles for different trapping strengths and mark the global density with black dashed lines.

t−∞
FIG. 7. Cartoon of a passive probe dragged by a constant, external force through a bath of active Langevin particles.

FIG. 8 .
FIG. 8. Drag mobility of the probe as a function Peclet number for probes in an active (orange) and passive (indigo) bath.The circles show the simulation data.a) Average velocity along the x-axis, ⟨ Vx⟩, as a function of the renormalized external drag force, Fext.The solid lines show the best linear fit for these data, the slope of which is μDrag in the linear regime.b) Force dependent mobility as a function of Peclet number.c), d) Zoom in on the mobility as a function of Peclet number in the c) active case and d) passive case (in c), the first data point is not shown).The dashed lines mark the value of μK.Light shading indicates error bars.Each bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.

FIG. 9 .
FIG. 9. Velocity autocorrelation function (VACF) of the dragged, immersed probe.Each bath has global density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.a,b) The top row shows the component of the VACF which is parallel to the drag force, whereas c,d) the middle row shows the perpendicular component.e,f ) The bottom row shows the difference between parallel and perpendicular components of the memory kernel.a,c,e) The left column is for a passive bath, whereas b,d,f ) the right column is for an active bath ( F0 = 5.3).The black curves in plots a-d) show the VACF for a free probe and insets show the difference between the VACF of a free and a dragged probe.Light shading indicates error bars.

FIG. 10 .
FIG. 10.Memory kernel of the dragged, immersed probe.Each bath has global density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.a,b) The top row shows the component of the memory kernel which is parallel to the drag force, whereas c,d) the middle row shows the perpendicular component.e,f ) The bottom row shows the difference between parallel and perpendicular components of the memory kernel.a,c,e) The left column is for a passive bath, whereas b,d,f ) the right column is for an active bath ( F0 = 5.3).The black curves in plots a-d) show the memory kernel for a probe which is not dragged and insets show the difference between the memory kernel of a free and a dragged probe.Light shading indicates error bars.

2 FIG. 11 .
FIG. 11.Stochastic force distribution of the dragged, immersed probe in a bath of density ρ = 0.4σ −3 .a,b) The top row shows the component of the stochastic force which is parallel to the drag force, whereas c,d) the bottom row shows one of the perpendicular components (the other perpendicular component is identical).a,c) The left column is for a passive bath, whereas b,d) the right column is for an active bath ( F0 = 5.3).The solid lines show simulation data, whereas the dotted lines show zero-centered Gaussian distributions with the same standard deviation.

2 FIG. 13 .
FIG.13.Dipole moment (Ω 0 1 (r), see Eq. (14)) of the bath as a function of distance from the center of the immersed probe (r) for a probe dragged with a Peclet number a) Pe ≈ 1.5, b) Pe ≈ 4.5, and c) Pe ≈ 7.2 for both a passive ( F0 = 0.0, indigo) and an active ( F0 = 5.3, orange) bath.Each bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.The black dashed line marks the zero baseline.

FIG. 14 .
FIG. 14. Effective drift force ( fx) as a function of the Peclet number, rescaled by C = µDrag/ Vdiff (see text for explanation) .Pe = C Fext (see Appendix E).The dashed, black line shows expected behavior fx = Fext in the linear regime.The inset shows fx/Fext as a function of the Peclet number, where the dashed lines mark fx/Fext = 1 (indigo) and fx/Fext = α0.4(orange).Each bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.

FIG. 15 .
FIG. 15.The renormalization constant α as a function of the bath activity F0.The dashed black line shows α = 1, what we would expect for a passive bath.The bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.

FIG. 16 .
FIG. 16.Memory kernel of a probe immersed in an active bath ( F0 = 5.3) and trapped by a harmonic potential with strength k, where the memory kernel is calculated a) without including the harmonic potential, b) using the imposed value of the harmonic constraint k rather than the renormalized value k.Light shading indicates error bars.Smaller graphs (a ′ ) and b ′ )) show the difference between the memory kernel of a free probe and a trapped probe, for the memory kernels shown in the graph immediately above.The bath has average density ρ = 0.4σ −3 and the probe has a radius Rp = 3.0σ.

5 FIG. 19 .
FIG. 19.Stochastic force distribution of the dragged, immersed probe in a bath of density ρ = 0.4σ −3 .a,b) The top row shows the component of the stochastic force which is parallel to the drag force, whereas c,d) the bottom row shows one of the perpendicular components (the other perpendicular component is identical).a,c) The left column is for a passive bath, whereas b,d) the right column is for an active bath ( F0 = 5.3).The solid lines show simulation data, whereas the dotted lines show zero-centered Gaussian distributions with the same standard deviation.