## Dynamics of a homogeneous active dumbbell system

Antonio Suma, 1 Giuseppe Gonnella, 2 Gianluca Laghezza, 2,3 Antonio Lamura, 4 Alessandro Mossa, 2 and Leticia F. Cugliandolo 5

1 SISSA - Scuola Internazionale Superiore di Studi Avanzati,

Via Bonomea 265, 34136 Trieste Italy ∗

2 Dipartimento di Fisica, Universit` a di Bari and

INFN, Sezione di Bari, via Amendola 173, Bari, I-70126, Italy

†

3 Current address: Department of Physics, Theoretical Physics, University of Oxford,

1 Keble Road, Oxford OX1 3NP, United Kingdom ‡

4 Istituto Applicazioni Calcolo, CNR, via Amendola 122/D,

Bari, I-70126, Italy §

5 Sorbonne Universit´ es, Universit´ e Pierre et Marie Curie - Paris 6,

Laboratoire de Physique Th´ eorique et Hautes Energies,

4, Place Jussieu, Tour 13, 5` eme etage, 75252 Paris Cedex 05, France ´ ¶

We analyse the dynamics of a two dimensional system of interacting active dumbbells. We characterise the mean-square displacement, linear response function and deviation from the equilibrium fluctuation-dissipation theorem as a function of activity strength, packing fraction and temperature for parameters such that the system is in its homogeneous phase. While the diffusion constant in the last diffusive regime naturally increases with activity and decreases with packing fraction, we exhibit an intriguing non-monotonic dependence on the activity of the ratio between the finite density and the single particle diffusion constants. At fixed packing fraction, the time-integrated linear response function depends non-monotonically on activity strength. The effective temperature extracted from the ratio between the integrated linear response and the mean-square displacement in the last diffusive regime is always higher than the ambient temperature, increases with increasing activity and, for small active force it monotonically increases with density while for sufficiently high activity it first increases to next decrease with the packing fraction. We ascribe this peculiar effect to the existence of finite-size clusters for sufficiently high activity and density at the fixed (low) temperatures at which we worked. The crossover occurs at lower activity or density the lower the external temperature. The finite density effective temperature is higher (lower) than the single dumbbell one below (above) a cross-over value of the P´ eclet number.

PACS numbers: 05.70.Ln, 47.63.Gd, 66.10.C-

## I. INTRODUCTION

Active matter is constituted by self-propelled units that extract energy from internal sources or its surroundings, and are also in contact with an environment that allows for dissipation and provides thermal fluctuations. The locally gained energy is partially converted into work and partially dissipated into the bath. The units can interact via potential forces or through disturbances in the medium. This new class of (soft) matter is the focus of intense experimental, theoretical and numerical studies for practical as well as fundamental reasons. Several review articles summarise the current understanding of active systems [1-9].

Due to the consumed energy, detailed balance is broken in active matter, and these systems are inherently out of equilibrium. Natural examples are bird flocks, schools of fish, and bacterial colonies. Artificial self-propelled particles have also been realised in the laboratory by using, for instance, granular materials [10, 11] or colloidal particles with specific surface treatments [12-14] and are especially suited for experimental tests. Different models have been proposed to mimic

∗ antonio.suma@sissa.it

† gonnella@ba.infn.it

‡ gianluca.laghezza@physics.ox.ac.uk

§ a.lamura@ba.iac.cnr.it

¶ leticia@lpthe.jussieu.fr

these systems. For instance, run-and-tumble motion is used to model Escherichia coli bacteria [5] while active Brownian particles are used to model Janus colloidal particles [15].

Self-propulsion is responsible for many interesting, and also sometimes surprising, collective phenomena. Some of these are: the existence of orientationally ordered states in two spatial dimensions [16, 17], spatial phase separation into an aggregated phase and gas-like regions for sufficiently large packing fractions in the absence of attractive interactions [1827], giant density fluctuations [19, 22, 28, 29] and accumulation at boundaries [30-32], spontaneous collective motion [33], glassy features [34-37], unexpected rheological properties [38] and non-trivial behavior under shear [39-43].

Active systems, being essentially out of equilibrium, also pose many fundamental physics questions such as whether thermodynamic concepts could apply to them in their original setting or with simple modifications.

The effective temperature notion was proposed to describe some macroscopic aspects of slowly relaxing macroscopic physical systems, such as glassy systems and gently sheared super-cooled liquids [44, 45]. This (potentially) thermodynamic intensive parameter is defined as the parameter replacing ambient temperature in the (multi) time dependent fluctuation-dissipation relations between induced and spontaneous out of equilibrium fluctuations of the system. In systems with multiple time-scales special care has to be taken in the choice of the time-regime in which a thermodynamiclike parameter could be extracted. More precisely, experi-

ence shows that one may identify it in the time-regime of (large) structural relaxation, while at short-time scales the microscopic dynamics imposes the system's fluctuations (be them quantum, active or thermal). To retain a thermodynamic sense, the effective temperature should also be measurable with suitable choices of thermometers such as well-chosen tracer particles and it should be the same for all observables evolving in the same time regime.

The effective temperature idea has been explored, to a certain extent, in the context of active matter. The effective temperature of a bacterial bath was estimated from the StokesEinstein relation of a tracer particle in [46]. The effective temperature notion was used to characterise crystallization effects known to occur under large active forces [34, 47, 48] and the emergence of collective motion [33]. The deviations from the equilibrium fluctuation-dissipation theorem in equilibrium were used to reveal the active process in hair bundles [49] and model cells [50]. In biological systems the nature of the microscopic active elements is difficult to study directly. The fluctuation-dissipation relations could be useful to characterise the active forces in active matter in general, and in living cells in particular. With this idea in mind, BenIsaac et al. analysed the perturbed and spontaneous dynamics of blood cells in the lab, and compared the outcomes to the ones of a single particle Langevin model with a special choice of the statistics of the active forces analytically [51]. Similar ideas were used in [52, 53] to characterise motor activity in living cells and actin-myosin networks.

We stressed the fact that the effective temperature is not a static parameter in the sense that it cannot be read from a system's snapshot. It should be determined from dynamic measurements in which the separation of time-scales has to be very carefully taken into account to obtain sensible results [44, 45]. Having said this, the effective temperature has been found to play a role similar to ambient temperature in the celebrated experiment of Perrin now performed with active particles. Indeed, the sedimentation of an active colloidal suspension of Janus particles under a gravity field exhibits the same exponential density distribution as a standard dilute thermal colloidal system. The only difference is that the parameter that replaces the thermal system's temperature in the active case is equal to the effective temperature inferred from an independent measurement of the long-time diffusive motion of an active colloidal particle [14]. This problem was studied analytically with a run-and-tumble model [54] and a Langevin process for a tagged active point-like particle with a suitable choice of activity [55].

Deviations from the equilibrium fluctuation-dissipation theorem, as well as other ways of measuring the effective temperature by using tracers, were analysed numerically by Loi et al. in (relatively loose) systems of active point-like particles [56, 57] and long molecules [57, 58]. All these measurements yielded consistent results. In this paper, we will follow this kind of analysis in a model of active matter that we now discuss.

In their simplest realisation, active units are taken to be point-like. However, the importance of the shape and polarity of the self-propelled particles for their collective be- haviour has been stressed in the literature [59-63]. Active units, whether synthetic or natural, are typically rodlike or elongated. This is the case of most bacteria, chemically propelled nanorods, and actin filaments walking on molecular motor carpets. The length-scale of these units is of the order of several micrometers.

A simple way to model a shortly elongated swimmer is to use a dumbbell, consisting of two colloids linked by a Hookean spring (see [64] and references therein). Such passive dumbbell models have been used to mimic the viscolelastic behavior of linear flexible polymers suspended in Newtonian fluids [65, 66]. Hydrodynamic interactions between dumbbell swimmers have been considered in [67-69]. In this study we add activity to Brownian dumbbells in the form of a constant propulsive force acting on the direction connecting the two beads. We include potential interactions between the dumbbells but we do not impose any alignment rule. This model was used to describe bacterial systems [70]. Compared with self-propelled spherical particle models, it phase separates at smaller densities [24, 25]. Clustering and phase separation are here due to the out of equilibrium drive exerted by the persistent local energy input that breaks detailed balance. Moreover, together with spontaneous aggregation, dumbbells break chiral invariance, and rotate displaying nematic order with spiral patterns.

We focus on the dynamic behaviour in the homogeneous phase of the two dimensional system. We study its dynamic properties at various fixed temperatures and dumbbell parameters as will be introduced below, but for a broad range of values of the surface density and strength of the active force. We analyse the behaviour of the single passive and active molecule analytically and we later use molecular dynamics simulations to study the many-body system. More precisely, we compute the mean-square displacement, linear response function, their relation and its implications on the effective temperature ideas [45] that we discuss in detail in the main text.

The body of the paper is organised as follows. In Sec. II we introduce the model. In Section III we give some details on the numerical algorithm that we use to study the problem numerically. Section IV is devoted to the study of the passive dumbbell system, both in its single molecule limit and manybody case. In Sec. V we present the analysis of a single active dumbbell molecule and in Sec. VI the one of a system of interacting and active dumbbell molecules. Finally, in Sec. VII we summarise our results and present our conclusions.

## II. THE MODEL

A dumbbell is a diatomic molecule made of two spherical colloids with diameter σ d connected by a spring of elastic constant k that one can mimic, in its simplest form, with Hooke's law

$$V _ { \text{H} } ( r ) = \frac { 1 } { 2 } k r ^ { 2 }$$

with r the distance between their centers of mass. An additional repulsive force is added, derived from just the repul-

sive part of a Lennard-Jones potential, that ensures that the two colloids cannot overlap. This potential is called WeeksChandler-Anderson (WCA) and it is given by [71]

$$V _ { \text{wca} } ( r ) = \begin{cases} V _ { \text{LJ} } ( r ) - V _ { \text{LJ} } ( r _ { c } ) & \quad r < r _ { c } \\ 0 & \quad r > r _ { c } \end{cases}$$

with

$$V _ { \text{LJ} } ( r ) = 4 \epsilon \left [ \left ( \frac { \sigma _ { \text{d} } } { r } \right ) ^ { 1 2 } - \left ( \frac { \sigma _ { \text{d} } } { r } \right ) ^ { 6 } \right ] \,. \quad \quad ( 2 )$$

glyph[epsilon1] is an energy scale and σ d is, once again, the diameter of the spheres in the dumbbell. r c is the minimum of the LennardJones potential, r c = 2 1 / 6 σ d . We neglect hydrodynamic interactions.

The equation of motion for a single dumbbell immersed in a liquid is the Langevin equation

$$m _ { \text{d} } \ddot { \mathbf r } _ { i } & = - \gamma \dot { \mathbf r } _ { i } - k ( \mathbf r _ { i } - \mathbf r _ { j } ) \\ & \quad - \frac { \partial V _ { \text{wca} } ( r _ { i j } ) } { \partial r _ { i j } } \frac { ( \mathbf r _ { i } - \mathbf r _ { j } ) } { r _ { i j } } + \eta _ { i } \quad \quad ( 3 )$$

with i, j = 1 2 , labelling the two spheres in the molecule, r i the position of the i -th monomer with respect to the origin of a Cartesian system of coordinates fixed to the laboratory, r ij = r i -r j and r ij = | r ij | . m d is the mass of each sphere in the dumbbell and γ is the friction coefficient. The Gaussian noise has zero mean and it is delta-correlated

$$\langle \eta _ { i a } ( t ) \rangle = 0 \, \quad \ \ \ \ ( 4 ) \quad \ a e p _ { \quad \ \. \. }$$

$$\langle \eta _ { i a } ( t ) \eta _ { j b } ( t ^ { \prime } ) \rangle = 2 \gamma k _ { B } T \delta _ { i j } \delta _ { a b } \delta ( t - t ^ { \prime } ) \, \quad ( 5 ) \quad \text{aom}$$

with k B the Boltzmann constant and T the temperature of an equilibrium environment in which the dumbbells move. a, b label the coordinates in d dimensional space. Note that an effective rotational motion is generated by the random torque due to the white noise acting independently on the two beads.

We add now the active force to Eq. (3). It acts in the direction of the spring linking the two colloids, i.e. in the direction ˆ n of the straight line passing by the two centers of mass, and it is constant in modulus. It reads

$$\mathbf F _ { \text{act} } = F _ { \text{act} } \, \hat { \mathbf n } \,.$$

Having established the single molecule stochastic model, we extend it to consider a system of N such bi-atomic molecules in interaction immersed in a bidimensional space with surface S . The molecule number density is n d = N/S , and the surface fraction is

$$\phi = N \, \frac { S _ { \text{d} } } { S }$$

with S d the area occupied by an individual dumbbell. The spring is supposed to be massless and void of surface. Therefore, in d = 2 we have S d = πσ / 2 d 2 .

In order to model the many-body system we introduce an inter-molecular potential and we slightly modify the interaction between the spheres in the same molecule. The interaction between dumbbells is purely repulsive and avoids the superposition of different molecules, i.e. ensures the excluded volume condition in the WCA way. We use then a LennardJones potential still truncated to have only the repulsive part as in Eq. (2). The elastic potential between spheres belonging to the same dumbbell is modified to be of the finite extensible non-linear elastic (FENE) kind to avoid the unlimited separation of the colloids belonging to the same molecule:

$$\mathbf F _ { \text{fene} } = \frac { k \mathbf r } { 1 - ( r ^ { 2 } / r _ { 0 } ^ { 2 } ) } \. \quad \quad ( 8 )$$

The denominator ensures that the spheres cannot go beyond the distance r 0 .

The dynamic equations for one dumbbell in the system are glyph[negationslash]

$$m _ { d } \ddot { \mathbf r } _ { i } ( t ) & = - \gamma \dot { \mathbf r } _ { i } ( t ) - \mathbf F _ { \text{fene} } ( \mathbf r _ { i, i + 1 } ) + \eta _ { i } \\ & \quad - \sum _ { \substack { j = 0, j \neq i \\ j = 0, j \neq i } } ^ { 2 N } \frac { \partial V _ { \text{wca} } ^ { i j } } { \partial r _ { i j } } \frac { \mathbf r _ { i j } } { r _ { i j } } + \mathbf F _ { \text{act} } _ { i } \,, \\ \therefore \quad \, \cdot \quad \, \cdot \quad \, \cdot \quad \, \cdot \,.$$

glyph[negationslash]

$$m _ { d } \ddot { \mathbf r } _ { i + 1 } ( t ) & = - \gamma \dot { \mathbf r } _ { i + 1 } ( t ) + \mathbf F _ { \text{fence} } ( \mathbf r _ { i, i + 1 } ) + \mathbf \eta _ { i + 1 } \\ & \quad - \sum _ { j = 0, j \neq i + 1 } ^ { 2 N } \frac { \partial V _ { \text{wca} } ^ { i + 1, j } } { \partial r _ { i + 1, j } } \frac { \mathbf r _ { i + 1, j } } { r _ { i + 1, j } } + \mathbf F _ { \text{act} } \,, ( 1 0 )$$

with i = 1 3 , , ... 2 N -1 and V ij wca ≡ V wca ( r ij ) with V wca defined in Eq. (2). The statistics of the noise η is Gaussian with average and correlation given by Eqs. (4) and (5), respectively. As the active force's direction lies along the molecular axis it depends on the diatomic molecule but is the same for the two atoms. This is the reason why we label it i in the two equations above. Note that once the active force is attached to a molecule a sense of back and forth atoms is attributed to them (see Fig. 1). The active forces are applied to all molecules in the sample during all their dynamic evolution. F act is a timedependent vector since, although its modulus is constant, it does change direction together with the molecule's rotation. For each dumbbell F act is directed from the colloid i (tail) to the colloid i + 1 (head). Note the difference between this kind of activity and the random one used in [57, 58] for the numerical study of active polymers.

Figure 1. (Color online.) A sketch of an active dumbbell molecule.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-afa4a18c-7ae0-423d-9e77-354b7d0835d1.jpg)

The P´ eclet number, Pe , is a dimensionless ratio between the advective transport rate and the diffusive transport rate. For particle flow one defines it as Pe = Lv/D , with L a typical length, v a typical velocity, and D a typical diffusion constant in the problem. We choose L → σ d , v → F act /γ

and D → D pd cm = k B T/ (2 γ ) of the passive dumbbell [see Eq. (21)]; then,

$$P e = \frac { 2 \sigma _ { d } F _ { \text{act} } } { k _ { B } T } \,. \quad \quad \quad ( 1 1 ) \quad \underset { \text{cif} } { \text{du} }$$

Another important parameter is the active Reynolds number

$$R e _ { \text{act} } = \frac { m F _ { \text{act} } } { \sigma _ { \text{d} } \gamma ^ { 2 } }$$

defined in analogy with the usual hydrodynamic Reynolds number Re = Lv/ν , where ν is the kinematic viscosity of a given fluid, representing the ratio between inertial and viscous forces. Here we set L → σ d , v → F act /γ and ν → γσ /m 2 d d .

## III. THE MOLECULAR DYNAMICS ALGORITHM

We solved the stochastic Langevin equation with an algorithm due to Vanden-Eijinden and Ciccotti [72] that is exact to order (∆ ) t 2 with ∆ t the time-step. We used a square bidimensional box with periodic boundary conditions.

The units of mass, length and energy are m d , σ d , and glyph[epsilon1] , respectively. The adimensional elastic constant is defined as k ∗ = kσ /glyph[epsilon1] 2 d and we take a large value k ∗ = 30 to avoid the excessive extension of the dumbbell, and also prevent vibrations. The length r 0 in the FENE potential is rendered adimensional as r ∗ 0 = r /σ 0 d and we used r ∗ 0 = 1 5 . . The thermal energy is normalized by the energy scale in the LennardJones potential, k B T ∗ = k B T/glyph[epsilon1] . The friction constant has dimensions of mass/time and we can make it adimensional as γ ∗ d = γ/ √ glyph[epsilon1]m d /σ 2 d . Physical realizations are usually in the over-damped regime which is ensured by choosing a value of γ ∗ d such that Re act glyph[lessmuch] 1 . We preferred to consider a model with inertial terms in order to have access to the crossover between ballistic and diffusive regimes. Concretely, we used γ ∗ d = 10 for which the molecular oscillations are strongly inhibited. Finally, the adimensional temperatures k B T ∗ are in the range 0 001 . -1 . Once all parameters are expressed in terms of reduced units we can effectively set m d = σ d = k B = glyph[epsilon1] = 1 .

The optimal choice of time-step is delicate. We first identified some relevant time-scales in the problem and we later chose ∆ t . The time-scale for the oscillations of the free dumbbell is τ osc = 2 π √ m /k d = 2 π/ √ k ∗ √ m σ /glyph[epsilon1] d 2 d (and equals 1.14 for our choice of parameters). The inertial time scale is t I = m /γ d = 1 /γ ∗ d √ m σ /glyph[epsilon1] d 2 d (and equals 0.1 for our choice of parameters). We will show below that there is another time-scale, typically longer, associated to the angular diffusion, t a = γR 2 ∞ / (2 k B T ) ≈ γ / ∗ d (2 k B T ∗ ) √ m σ /glyph[epsilon1] d 2 d where we used R ∞ ≈ σ d the typical length of the dumbbell molecule. We chose to work with ∆ = t t I / 100 to see details of the ballistic regime and ∆ = t τ osc / 100 to enter the later diffusive and active regimes. Times are measured in units of √ ( σ m 2 d d ) /glyph[epsilon1] .

## IV. PASSIVE SYSTEM

In this Section we review the behaviour of the passive dumbbell single molecule and interacting system, paying special attention to the dependence of the dynamic regimes on the surface fraction φ .

## A. Passive single dumbbell

One can simply show that the equation of motion for the position of the center-of-mass, r cm = ( r 1 + r 2 ) / 2 , of a single dumbbell governed by Eq. (3) under the same force f acting on each bead is

$$2 m _ { \text{d} } \ddot { \mathbf r } _ { \text{cm} } ( t ) = - 2 \gamma \dot { \mathbf r } _ { \text{cm} } ( t ) + 2 \mathbf f ( t ) + \xi ( t ) \quad ( 1 3 )$$

with the new noise ξ ( ) t ≡ η 1 ( ) + t η 2 ( ) t with vanishing average, 〈 ξ a ( ) t 〉 = 0 , and correlation

$$\langle \xi _ { a } ( t ) \xi _ { b } ( t ^ { \prime } ) \rangle = 4 \gamma k _ { B } T \, \delta _ { a b } \delta ( t - t ^ { \prime } ) \. \quad \ ( 1 4 )$$

This is the Langevin equation of a point-like particle with mass 2 m d , under a force 2 f in contact with a bath with friction coefficient 2 γ at temperature T . Equivalently, one can divide this equation by two and obtain a Langevin process for a point-like particle with mass m d under a force f in contact with a bath with friction coefficient γ at temperature T/ 2 . From these analogies one recovers several results on the statistics of the center-of-mass position and velocity. Under no external force, f = 0 , the center-of-mass velocity is distributed according to the Maxwell distribution for a particle with mass 2 m d at equilibrium at temperature T (or mass m d at temperature T/ 2 ). The center-of-mass mean-square displacement between two times t ′ and t after preparation

$$\Delta ^ { 2 } ( t, t ^ { \prime } ) \equiv \langle [ { \mathbf r } _ { \mathrm c m } ( t ) - { \mathbf r } _ { \mathrm c m } ( t ^ { \prime } ) ] ^ { 2 } \rangle \quad \quad ( 1 5 )$$

can be calculated as

$$\text{can be calculated as} \\ \Delta ^ { 2 } ( t, t ^ { \prime } ) = d & \Big { [ } \Big { ( } v _ { 0 } ^ { 2 } - \frac { k _ { B } T } { 2 m _ { d } } \Big { ) } \frac { ( e ^ { - \frac { \gamma } { m _ { d } } t ^ { \prime } } - e ^ { - \frac { \gamma } { m _ { d } } t } ) ^ { 2 } } { ( \gamma / m _ { d } ) ^ { 2 } } \\ & \quad + \frac { k _ { B } T } { \gamma } ( t - t ^ { \prime } ) \\ & \quad - \frac { m _ { d } k _ { B } T } { \gamma ^ { 2 } } ( 1 - e ^ { - \frac { \gamma } { m _ { d } } ( t - t ^ { \prime } ) } ) \Big { ] }. \quad \quad ( 1 6 ) \\ \text{where} \, v _ { n } \text{ is the velocity of the particle at the initial time } t = 0.$$

where v 0 is the velocity of the particle at the initial time t = 0 . Given the inertial time

$$t _ { I } = m _ { \text{d} } / \gamma \,,$$

one obtains different limits in relation to different values of t and t ′ . At short times 0 ≤ t ′ ≤ t glyph[lessmuch] t I , by expanding all exponentials at small arguments, we obtain ballistic behaviour,

$$\Delta ^ { 2 } ( t, t ^ { \prime } ) = d v _ { 0 } ^ { 2 } \, ( t - t ^ { \prime } ) ^ { 2 } \. \quad \quad \ \ ( 1 8 )$$

At long total times t ≥ t ′ glyph[greatermuch] t I and short time-delay ( t -t ′ ) glyph[lessmuch] t I , one also obtains ballistic behaviour,

$$\Delta ^ { 2 } ( t, t ^ { \prime } ) = d \, \frac { k _ { B } T } { 2 m _ { d } } \, ( t - t ^ { \prime } ) ^ { 2 } \, \quad \quad ( 1 9 )$$

with the initial velocity v 2 0 replaced by its average in equilibrium 〈 v 2 〉 = k B T/ (2 m d ) (for a particle with mass 2 m d ). In both cases the time-delay dependence crosses over to diffusive motion

$$\Delta ^ { 2 } ( t, t ^ { \prime } ) = 2 d D _ { \text{cm} } ^ { \text{pm} } \left ( t - t ^ { \prime } \right ) \,, \quad \quad ( 2 0 ) \quad \text{and}$$

for long time delay ( t -t ′ ) glyph[greatermuch] t I , with a diffusion constant taking the form

$$D _ { c m } ^ { \text{pd} } = k _ { B } T / ( 2 \gamma ) \. \quad \ \ ( 2 1 ) \quad \ \ r e _ { F } ^ { \text{cl} }$$

This is the diffusion constant used in the definition of the P´ eclet number in Eq. (11).

The length of the dumbbell molecule, R ( ) = t r 1 -r 2 , with r 1 and r 2 the position of the centers of the two spheres, is also a fluctuating quantity. For the parameters used one shows that 〈 R t ( ) 〉 approaches R ∞ = 0 96 . ≈ σ d . The angular degrees of freedom can also be simply analyzed, especially under the assumption that R t ( ) is constant which is rather accurate since R t ( ) does not fluctuate more than 3% around its mean value. In this approximation one finds that the angle diffuses with an angular diffusion constant equal to

$$D _ { a } = 2 k _ { B } T / ( \gamma R _ { \infty } ^ { 2 } ) \. \text{ \quad \ \ } ( 2 2 )$$

The linear instantaneous response, R , quantifies the effect of a small impulsive perturbation, say h ( t ′′ ) applied at time t ′′ , on the dynamics of the system of interest. For a dumbbell moving in a plane, the a -th component of the center-of-mass position perturbed by a force acting on its b -th component, from time t ′ until time t , is

$$\langle r _ { \text{cm} } ^ { a } ( t ) \rangle _ { \text{h} } = \langle r _ { \text{cm} } ^ { a } ( t ^ { \prime } ) \rangle + \int _ { t ^ { \prime } } ^ { t } d t ^ { \prime \prime } \, R _ { a b } ( t, t ^ { \prime \prime } ) \, h _ { b } ( t ^ { \prime \prime } ) \,. \quad ( 2 3 )$$

Equivalently,

$$R _ { a b } ( t, t ^ { \prime } ) \equiv \frac { \delta \langle r _ { c m } ^ { a } ( t ) \rangle _ { \mathbf h } } { \delta h _ { b } ( t ^ { \prime } ) } \Big | _ { \mathbf h = 0 } \. \quad \quad ( 2 4 )$$

Motivated by the interpretation of Eq. (13) that keeps the temperature of the noise unaltered and equal to T , we take the perturbation to the center-of-mass to be h ( t ′′ ) = 2 ( f t ′′ ) . In order to focus on the long time limit of interest, t glyph[greatermuch] t I , and to simplify the expressions, we drop the inertia term from Eq. (13). In this limit, the center-of-mass position is given by

$$\mathbf r _ { \text{cm} } ( t ) = \mathbf r _ { \text{cm} } ( t ^ { \prime } ) + \frac { 1 } { 2 \gamma } \int _ { t ^ { \prime } } ^ { t } d t ^ { \prime \prime } \left [ 2 \mathbf f ( t ^ { \prime \prime } ) + \xi ( t ^ { \prime \prime } ) \right ]. \quad ( 2 5 ) \quad \text{dumb} \right ]$$

The total linear response to a step-like perturbation that is applied from t ′ to t , with independent components 2 f a on each Cartesian spatial direction such that δ (2 f a ( t ′′ )) /δ (2 f b ( t ′ )) = δ ab , is then

$$\chi ( t, t ^ { \prime } ) \equiv \int _ { t ^ { \prime } } ^ { t } d t ^ { \prime \prime } \, R _ { a a } ( t, t ^ { \prime \prime } ) = \frac { d } { 2 \gamma } \, ( t - t ^ { \prime } ) \,. \quad ( 2 6 ) \quad \text{In}$$

Summation over repeated indices was used to go from the second to third member of this equation. Comparing now to the mean-square displacement, ∆ ( 2 t, t ′ ) , one finds that in the diffusive regime

$$2 k _ { B } T \chi ( t, t ^ { \prime } ) = \Delta ^ { 2 } ( t, t ^ { \prime } ) \, \quad \quad ( 2 7 )$$

and both functions depend on the two times only through their difference, t -t ′ . This relation is the same as the one for a point-like Brownian particle in contact with a bath at temperature T . We will use it as a reference to define the effective temperature from deviations in which the bath temperature is replaced by another parameter. (The comparison between χ and the time-derivative of the correlation function between the position measured at different times, which in equilibrium is proportional to the inverse temperature of the bath, yields a non-trivial relation that violates the equilibrium FDT for an unconfined Brownian particle [73].)

Anumerically convenient way to extract the linear response from Eq. (23), especially useful for the interacting systems studied in the following, consists in taking the perturbing forces f to be uncorrelated with the center-of-mass position and random, with zero mean, [ f a ] = 0 , and correlation, [ f f a b ] = f 2 δ ab . Multiplying Eq. (23) by 2 f c and taking the average over their distribution one finds

$$[ \langle r _ { \text{cm} } ^ { a } ( t ) \rangle 2 f _ { c } ] = ( 2 f ) ^ { 2 } \int _ { t ^ { \prime } } ^ { t } d t ^ { \prime \prime } R _ { a c } ( t, t ^ { \prime \prime } ) \. \quad ( 2 8 )$$

Setting a = c and summing over components one obtains

$$\chi ( t, t ^ { \prime } ) = \frac { [ \langle r _ { c m } ^ { a } ( t ) \rangle 2 f _ { a } ] } { ( 2 f ) ^ { 2 } } \. \quad \quad ( 2 9 )$$

Then, one derives from Eq. (25)

$$\left [ \langle r _ { \text{cm} } ^ { a } ( t ) \rangle 2 f _ { a } \right ] = \frac { 1 } { 2 \gamma } \, \delta _ { a a } \, ( 2 f ) ^ { 2 } \, ( t - t ^ { \prime } ) \, \quad ( 3 0 )$$

that yields

$$\chi ( t, t ^ { \prime } ) = \frac { d } { 2 \gamma } ( t - t ^ { \prime } ) \, \text{ \quad \ \ } ( 3 1 )$$

consistently with Eq. (26).

The way in which we probe the linear response of the interacting system is by applying a force of modulus 2 f in a random direction to the center-of-mass, and not on the rotational and vibrational degrees of freedom, which gives a linear response proportional to 1 / (2 γ ) in the case of a single passive dumbbell. Averaging over different angles, one recovers a result proportional to 1 / (2 γ ) . To obtain the linear response to independent perturbations applied on the spatial d directions, that is to say χ , we simply multiply the result by d .

## B. Passive many-body dumbbell system

In order to study the many-body system we performed numerical simulations in the form described in Sec. III. We averaged over 100 different realizations of a system with linear size l = 100 and adimensional temperature T = 0 001 . (to

make the notation lighter here and in what follows we avoid the asterisks to denote adimensional parameters). Measurements were done after an equilibration time of order t eq = 10 2 , starting from an initial random configuration of our system, in both the position of the center-of-mass and the angular direction of the dumbbells, and we show data gathered until a time equal to 10 6 .

Figure 2. (Color online.) Mean-square displacement as a function of the delay time hereafter called t measured after t ′ = 100 , in a system of interacting passive dumbbells with the surface fractions given in the key at T = 0 001 . (double logarithmic scale).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-bb71e66d-d417-4708-9c09-9a7b26ccfdc5.jpg)

First, we studied a very loose system, with φ = 0 01 . (data not shown). From the analysis of the system's global meansquare displacement we found that the ballistic regime ends at the inertial time t I = m /γ d glyph[similarequal] 0 1 . when it crosses over to a diffusive regime. The ballistic regime is characterized by ∆ ( ) 2 t glyph[similarequal] d v 〈 2 〉 t 2 with 〈 v 2 〉 = k B T/ (2 m d ) the thermal velocity of a Brownian particle with mass 2 m d , which for the parameters used takes the value 〈 v 2 〉 glyph[similarequal] 5 10 · -4 . The diffusion constant in the free-diffusive regime is very close to the one of the single passive dumbbell, D cm glyph[similarequal] D pd cm glyph[similarequal] 5 · 10 -5 .

Next, we studied systems with five higher surface fractions: φ = 0 1 . , 0 2 . , 0 3 . , 0 4 . , 0 5 . , see Fig. 2. We found that, as in the single molecule case, the ballistic regime crosses over at t I = m /γ d = 0 1 . to a first free diffusive regime, approximatively in the interval 1 ≤ t glyph[lessorsimilar] 10 2 , that is now followed by a second diffusive regime that feels the effect of the dumbbell concentration [74]. Quite naturally, the dynamics slows down for denser systems as shown by the spread of curves at the late stages in the plot.

In Fig. 3 data are represented in such a way that the last diffusive regime appears as a plateau at 1 /t → 0 with height equal to the diffusive constant. The diffusive constant in the interacting regime, normalized by the one of the center-ofmass of the free passive dumbbell, D pd cm , is plotted in Fig. 4 against the surface fraction. For comparison we also similarly display normalised data obtained for a system of colloids coupled to the same thermal bath and interacting via the WCA potential (2). Within numerical accuracy the colloidal and dumbbell systems have the same behaviour until φ glyph[lessorsimilar] 0 3 . while the dumbbell data are slightly above the colloidal ones

Figure 3. (Color online.) Rescaled mean-square displacement in systems of interacting passive dumbbells with the surface fractions given in the key. The data in Fig. 2 are divided by 2 dt and presented as a function of 1 /t to identify the diffusive regime as a plateau with height given by the diffusion constant. The single passive dumbbell diffusive constant takes the value D pd cm = 5 10 -5 for these parameters, consistently with the trend of the numerical results for the many-particle system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-c9923653-eb14-43a2-be03-024685e723bb.jpg)

for larger values of φ . We find very good agreement with the prediction of Tokuyama and Oppenheim,

$$\frac { D _ { \text{cm} } ( \phi ) } { D _ { \text{cm} } ( 0 ) } = \frac { 1 } { 1 + H ( \phi ) } \, \text{ \quad \quad \ } ( 3 2 )$$

with and

$$b ( \phi ) = \sqrt { 9 \phi / 8 } \quad c ( \phi ) = 1 1 \phi / 1 6 \,, \quad \ \ ( 3 4 )$$

derived from a perturbative calculation for low-concentrated colloidal systems [75]. This expression has no free parameters. The agreement is very good for all surface fractions for colloids, and until φ glyph[similarequal] 0 3 . for dumbbells (continuous line in the figure) if we use the three terms above in the expansion. For denser systems the shape of the molecules starts playing a role. We stress that the comparison is made between normalized curves. The third and fourth sets of data in this figure correspond to systems of active dumbbells and will be discussed in Sec. VI.

## V. ACTIVE SINGLE MOLECULE

We turn now to the dynamics of a single active molecule as a reference case for the interacting system to be discussed in the next section.

$$H ( \phi ) & = \frac { 2 [ b ( \phi ) ] ^ { 2 } } { 1 - b ( \phi ) } - \frac { c ( \phi ) } { 1 + 2 c ( \phi ) } \\ & \quad - \frac { b ( \phi ) \, c ( \phi ) \, [ 2 + c ( \phi ) ] } { [ 1 + c ( \phi ) ] \, [ 1 - b ( \phi ) + c ( \phi ) ] } \quad \text{(33)}$$

Figure 4. (Color online.) Upper panel: in linear-log scale, the diffusion constant of an interacting system, as a function of surface fraction φ , in four cases: the passive dumbbell system displayed with (green) crosses × ; the passive colloidal system shown with (red) pluses + ; the active dumbbell system with F act = 0 001 . and F act = 0 1 . drawn with (pink) triangles glyph[triangle] and (blue) squares glyph[square] , respectively. In all cases the systems are in contact with the same thermal bath at temperature T = 0 001 . and interact with the same WCApotential given in Eq. (2). Lower panel: in double linear scale, the same data relative to its single particle limit, as a function of surface fraction φ . The solid line represents the analytic prediction by Tokuyama-Oppenheim (TO) for a colloidal system with the same parameters [75].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-a0f2c1dc-51bf-4c77-9b93-d90e5a687a33.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-76e4b128-297c-42da-8301-72bc70b89d56.jpg)

Take a single active molecule constrained to move in two dimensional space. The length of the molecule is not altered by the active force and it still approaches R ∞ ≈ σ d . From the Langevin equation for the center-of-mass position that acquires a forcing term one calculates the average square velocity of the center-of-mass. For t glyph[greatermuch] t I the stationary value for its a -th component is given by

$$\langle v _ { c m a } ^ { 2 } ( t ) \rangle & \to \frac { 1 } { 2 m _ { d } } \left ( k _ { B } T + \frac { F _ { \text{act} } ^ { 2 } } { \gamma ( \frac { 1 } { t _ { I } } + \frac { 1 } { t _ { a } } ) } \right ) \quad \text{$\max_{\text{live}$}$} } { \quad } \quad \text{than} \\ & \equiv \frac { k _ { B } T _ { \text{kin} } } { 2 m _ { d } }$$

with the inertial and angular time-scales

$$t _ { I } = \frac { m _ { \text{d} } } { \gamma } \, \quad \ \ t _ { a } = \frac { 1 } { D _ { a } } = \frac { \gamma \sigma _ { \text{d} } ^ { 2 } } { 2 k _ { B } T } \, \quad ( 3 6 )$$

that will find a clear meaning when studying the mean-square displacement below.

The formula (35) gives an expression to what is called the kinetic or granular temperature [76-78] of the system. Quite naturally, T kin = T for F act = 0 . Note that T kin depends on the two time scales t I and t a . It approaches

$$\bigcap _ { \alpha \, e } \quad k _ { B } T _ { \text{kin} } \simeq \begin{cases} \ k _ { B } T + t _ { I } F _ { \text{act} } ^ { 2 } / \gamma \quad & \text{for} \quad & t _ { a } \gg t _ { I } \,, \\ \ k _ { B } T + t _ { a } F _ { \text{act} } ^ { 2 } / \gamma \quad & \text{for} \quad & t _ { I } \gg t _ { a } \,. \end{cases}$$

In our problem t a glyph[greatermuch] t I and the relevant limit is

$$k _ { B } T _ { \text{kin} } = k _ { B } T + \frac { m _ { \text{d} } } { \gamma ^ { 2 } } F _ { \text{act} } ^ { 2 } \. \quad \quad ( 3 7 )$$

Moreover, for the parameter values later used in the simulation, m /γ d 2 = 0 01 . .

Figure 5. (Color online.) Mean-square displacement of active dumbbells for various values of F act given in the key at temperature T = 0 001 . and linear size of the system l = 100 . The behaviour of a single dumbbell is shown with thinner lines while the one of a system of interacting dumbbells with φ = 0 1 . is shown with thicker lines. The average has been taken over 200 thermal histories. The inertial and angular time-scales are signaled with vertical black lines at t I = 0 1 . and t a = 5000 , respectively. The floating time-scale t ∗ is proportional to the inverse square active force, F -2 act . See the text for a discussion of its effect.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-fa5a64ce-3314-46d5-9ba9-e3583db7320e.jpg)

After some simple calculations one finds that the meansquare displacement of a single dumbbell molecule moving on a plane under the longitudinal active force has an initial ballistic regime for times shorter than the inertial time-scale t I . Beyond this time-scale, the over-damped limit takes over and inertia can be neglected. One can then quite easily derive the time-behaviour of the self-diffusion. For times shorter than the angular time scale t a one finds

$$\Delta ^ { 2 } ( t ) \simeq 2 d D _ { \text{cm} } ^ { \text{pd} } \, t + \left ( \frac { F _ { \text{act} } } { \gamma } \right ) ^ { 2 } t ^ { 2 } \, \quad \ \ ( 3 8 )$$

with D pd cm = k B T/ (2 γ ) the diffusion constant of the passive diatomic molecule (see also [79] where a similar calculation for an active ellipsoid has been performed). Within this regime one can still identify two sub-regimes. For times shorter than t ∗ = 2 dD pd cm γ /F 2 2 act the first term dominates and the dynamics is diffusive as in the absence of the active force. If, instead, times are longer than t ∗ the dynamics becomes ballistic again and it is controlled by the active force. For even longer times, going beyond t a , a new diffusive regime establishes,

$$\Delta ^ { 2 } ( t ) \simeq 2 d D _ { A } \, t \, \quad \ \ \ ( 3 9 ) \quad \text{in}$$

with a diffusion constant that depends upon the active force and in d = 2 is equal to

$$D _ { A } = D _ { \text{cm} } ^ { \text{pd} } + \frac { 1 } { 2 } \left ( \frac { F _ { \text{act} } } { \gamma } \right ) ^ { 2 } \frac { 1 } { D _ { a } }. \quad \ \ ( 4 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

In order to make the parameter dependence explicit we call this diffusion constant

$$D _ { A } ( F _ { \text{act} }, \phi = 0 ) = \frac { k _ { B } T } { 2 \gamma } \left [ 1 + \frac { 1 } { 2 } \left ( \frac { F _ { \text{act} } \sigma _ { d } } { k _ { B } T } \right ) ^ { 2 } \right ] \quad ( 4 1 ) \quad \text{a sm} \\ \quad \text{cle(s)}$$

with R ∞ replaced by σ d . This expression shows a dependence on the square of the P´ eclet number that will later appear in the effective temperature as well. The reason why the dynamics slow down in this last regime although there is still the active force acting on the dumbbell is that the angular motion goes against the translational one. All these regimes are summarized in

↦

$$\text{ballistic} & \mapsto \text{diffusive} \mapsto \text{ballistic} & \mapsto \text{diffusive} \\ t _ { I } & \quad t ^ { * } \quad t _ { a }$$

↦

↦

with the time-scales t I and t a introduced in Eq. (36). Note that the 'intermediate' time-scale t ∗ is inversely proportional to the square of the active force, F -2 act , and can go below the inertial time-scale t I for sufficiently strong non-equilibrium forcing or beyond t a for sufficiently weak one.

The generic features described in the previous paragraph can be seen in Fig. 5, especially in the curves for F act = 0 1 . , 0 01 . , where t ∗ &lt; t a and the regimes separated by t ∗ are distinct. t ∗ goes below t I for the largest active force F act = 1 used in Fig. 5 for which the data shown have already reached the second ballistic regime and cross over, at t a , to the last diffusive regime. At the other extreme, for the smallest active force F act = 0 001 . the time-scale t ∗ is very large going beyond the time-scale t a . The data show the first ballistic regime, the cross-over to the first free-diffusive regime and a very smooth cross-over to the last diffusive regime with a different diffusion constant. There is no second ballistic regime in this case. For the intermediate active forces, F act = 0 1 . , 0 01 . , the four regimes can be observed in the data curves. The numerical results yield values for the crossover times, diffusion constants and velocities in the ballistic regimes that are indistinguishable - within numerical accuracy - from the analytic predictions.

A similar crossover in the mean-square displacement was observed in systems and models of spherical particles. On the theoretical side, reference to such a crossover is made in [19, 26, 79]. In [80] the dynamics of artificial swimmers made of polystyrene microspheres coated with platinum on one size, while keeping the second half as the non-conducting polystyrene, was studied. In this system, the platinum catalyzes the reduction of a 'fuel' of hydrogen peroxide to oxygen and water propelling the particles in a preferred direction. Particle tracking was used to characterise the particle motion as a function of hydrogen peroxide concentration. At short times, the particles move predominantly in a directed way, with a velocity that depends on the concentration of the fuel molecules, while at long times the motion randomises and becomes diffusive, with a diffusion constant that also depends on the fuel concentration. The analytic results for the single dumbbell mean-square displacement are consistent with the experimental results. A similar crossover was also observed in the motion of beads propelled by adsorbed bacteria [81] and in bacterial baths [46].

Another, independent probe of the dynamics of an out of equilibrium system is the displacement induced by applying a small constant external force to one (or more) tagged particle(s) in the system:

$$\text{$\quad $\quad$\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ \ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad \ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ \quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{\$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{ $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ \ \ $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\	 $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ $ $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ :: $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ = \quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\ =\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$.\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \text{$\quad $\ } \$$

This equation generalises (29) to the many-body system.

The perturbing force f i = glyph[epsilon1] i f is applied to every monomer (in the same way as the active force) at time t 0 and kept constant until the measuring time t . f is its modulus and glyph[epsilon1] i its direction which is uniformly distributed and is the same for the two monomers of a given dumbbell. In the case of a single active dumbbell N = 1 and one finds again

$$\chi ( t ) = d ( 2 \gamma ) ^ { - 1 } \, t$$

a result that is independent of the activation force F act . We will write it as

$$\chi ( t ) = d \mu ( F _ { \text{act} }, \phi = 0 ) \ t \text{ \quad \ \ } \text{(44)}$$

with µ F ( act , φ = 0) = µ F ( act = 0 , φ = 0) . Using the relation between mean-square displacement, ∆ 2 , and induced displacement, χ , in equilibrium, 2 k B Tχ t ( ) = ∆ ( ) 2 t for all t , to define a possibly time-dependent effective temperature out of equilibrium

$$2 k _ { B } T _ { \text{eff} } ( t ) \chi ( t ) = \Delta ^ { 2 } ( t ) \, \quad \quad \ \ ( 4 5 )$$

one finds, in the late diffusive regime t glyph[greatermuch] t a ,

$$k _ { B } T _ { \text{eff} } ( t ) & \to k _ { B } T _ { \text{eff} } ( F _ { \text{act} }, \phi = 0 ) \\ & = \frac { D _ { A } ( F _ { \text{act} }, \phi = 0 ) } { \mu ( F _ { \text{act} }, \phi = 0 ) }$$

that becomes

$$\text{tan becomes} \\ T _ { \text{eff} } ( F _ { \text{act} }, \phi = 0 ) & = T \left [ 1 + \frac { 1 } { 2 } \left ( \frac { F _ { \text{act} } \sigma _ { \text{d} } } { k _ { B } T } \right ) ^ { 2 } \right ] \quad \text{Th} \\ & = T \left ( 1 + \frac { 1 } { 8 } \text{Pe} ^ { 2 } \right ) \,. \quad \text{(47)} \quad \text{is ve} \\ \text{mean} \\ \text{We can error} \, \text{if the same} \, \text{if in} \, \text{$T$} \, from the $E$} \, \text{in} \, \text{$Later}$$

We can arrive at the same result by defining T eff from the Einstein relation between diffusion constant and mobility. We stress the fact that T eff takes a constant value in this late diffusive regime. Note the different bath-temperature dependence in T eff and T kin . We reckon also that T eff ≥ T kin and that the two expressions coincide in the (unphysical) case t I glyph[greatermuch] t a . For the parameter values later used in the simulations k B T eff = k B T +0 5 . F 2 act / k ( B T ) and the factor in the F 2 act term in k B T eff equals 10 at T = 0 05 . (while the one in k B T kin is only 0.01).

The response-displacement relation in a stochastic model for an active particle was studied in [51]. In this paper, an overdamped Langevin equation for a randomly kicked pointlike particle with a variable number of kicking motors producing pulses of force ± f 0 acting at Poisson distributed times during intervals of duration ∆ τ , and in contact with a thermal environment, was analyzed. In the low frequency limit (long time-delays) the effective temperature approaches a constant T eff → T + N f m 2 0 ∆ τ 2 / τ ( +∆ ) τ with N m the number of motors and τ the averaged waiting-time between the pulses. The dependence on f 2 0 is similar to the dependence on Pe 2 that we find for the single active dumbbell.

Again from the T eff perspective, Szamel [55] studied a different model in which the self-propulsion of a single pointlike active particle is modelled as a fluctuating force evolving according to the Ornstein-Uhlenbeck process, independently of the state of the particle. The particle moves in a viscous medium that is assumed to force overdamped motion. The free diffusion properties of the particle lead to an effective temperature defined from an extension of the Einstein relation linking the diffusion coefficient of the free particle to the variance of the fluctuating term in the OrnsteinUhlenbeck process for the active force and the friction coefficient of the medium. The mobility and self-diffusion of an active Janus particle were monitored with micro-tracking by Palacci et al. [14] to infer from them the effective temperature. In both cases, as for our active single dumbbell, the effective temperature increases with the activation parameter as Pe 2 .

Palacci et al. [14], Tailleur and Cates [18, 54] and Szamel [55] measured and calculated the stationary probability distribution of the active particles' positions in a linear external potential (mimicking gravity) in a regime such that the sedimentation velocity is small with respect to the swimming velocity. They found an exponential form that corresponds to a Boltzmann distribution under gravity, with the parameter associated to the sedimentation length given by the equilibrium one with the temperature being replaced by the effective one of the free active particle. Szamel also studied whether the ambient temperature is simply replaced by the single particle effective one in the position probability distribution function of the Ornstein-Uhlenbeck active particle under a different external potential and found that this is not the case [55].

## VI. ACTIVE MANY-BODY DUMBBELL SYSTEM

The collective behavior of an ensemble of dumbbell molecules in interaction and under the effect of active forces is very rich [24, 25]. In this section we first describe some general features of this behavior in the homogeneous phase. Later, we will focus on the evaluation of the diffusion constant, linear response and effective temperature.

## A. Homogeneous phase

Figure 6. (Color online.) Schematic phase diagram of the dumbbell system. Inside the curve, at high P´ eclet numbers, the system undergoes phase separation into two phases characterized by two different densities. Here, typically, large and stable clusters are observed, as shown in the snapshot on the right. The dumbbells are freezed and point preferentially towards the center of the cluster. The enlargement shows the border of such a cluster, with the green part, corresponding to the head of the dumbbells, pointing inside, and the tails (red part) pointing outside. For small P´clet numbers the system e does not show the formation of such large clusters and a typical configuration is presented. The location of the transition line is based on the results of simulations with F act = 1 (see [24, 25] for further details); the right snapshot is taken at T = 0 01 . , the one on the left at T = 0 05 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-314e12e2-74b0-4abf-acff-d017055aae2e.jpg)

Active dumbbell systems exhibit a transition between a phase at small P´ eclet number, that includes the passive limit, and a phase with stable aggregates of dumbbells, not existing without activity, at high P´ eclet number. A schematic representation of the phase diagram is given in Fig. 6. We will not discuss the details of this phase diagram here but we simply state that we will work with sufficiently low P´ eclet number so as to stay in the region that we call the homogeneous phase although clustering phenomena, that will be important for the effective temperature, are already present. Typical configurations for the two phases are also shown in Fig. 6. In the aggregated phase the clusters reach the size of the system, while in the homogeneous phase they remain of finite size and are not stable. The phase separating kinetics in the phase with

large-scale aggregation was studied in [24]. One can also observe that, while in the clusters of the homogeneous phase the dumbbells are not stucked and can move, they are frozen inside the aggregates in the high P´ eclet phase.

Figure 7 shows six typical snapshots of the system, all at the same density φ = 0 4 . and temperature T = 0 05 . but for different active forces, F act = 0 01 . , 0 1 . , 0 3 . , 0 5 . , 0 7 . , 1 (in reading order). The corresponding P´ eclet numbers vary in the interval [0 4 . , 40] which is well inside the homogenous phase in the phase diagram of Fig. 6. On the four last snapshots, i.e. beyond F act glyph[similarequal] 0 3 . we start seeing clustering effects. Their presence becomes less relevant moving away from the critical surface in parameter space. Configurations of the system for different densities at fixed active force strength are shown in Fig. 8.

In Fig. 9 we show the static structure factor

$$S ( q ) = \frac { 1 } { N } \sum _ { i = 1 } ^ { N } \langle | e ^ { i { \mathbf q } \cdot { \mathbf r } _ { i } } | ^ { 2 } \rangle \quad \quad ( 4 8 )$$

of a system with φ = 0 1 . (upper panel) and φ = 0 4 . (lower panel) for different active forces, F act = 0 01 . , 0 1 . , 0 3 . , 0 5 . , 0 7 . , 1 following the line code given in the keys. The positions r i here are taken to be the ones of the two beads in the diatomic molecule. This function char-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-e3158a0b-de7c-466e-a279-6ceaa10671eb.jpg)

Figure 7. (Color online.) Clustering effects in the active manybody system. Density is φ = 0 4 . , temperature is T = 0 05 . and from left to right and top to bottom the active forces are F act = 0 01 . , 0 1 . , 0 3 . , 0 5 . , 0 7 . , 1 . The side of the box is l = 100 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-9b4c1a37-e792-4318-9420-9e40ec9c8c14.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-542cf012-2bd4-4d48-865f-b522c8eab02d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-cf4b2240-41ca-4d44-a958-3544a1cfc61e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-330d6234-b0a9-4259-b6e9-ecaf098a8c9f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-f2bdf314-02f0-4aa2-980d-71941d4d6503.jpg)

Figure 8. (Color online.) Clustering effects in the active many-body system. In all snapshots the active force strength is F act = 0 5 . and the temperature is T = 0 05 . , while the densities are φ = 0 1 . , 0 2 . , 0 3 . , 0 4 . (in reading order). The side of the box is l = 100 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-6f0d3c86-d4aa-4c09-855e-517194c6202a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-8b932dea-8efd-46b1-aeb6-83975a09567a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-0241be31-608c-412c-b1af-69a0d01b69a7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-58176dc1-b8e5-4570-8d7e-c0c974c52191.jpg)

acterises the strength of density fluctuations at a length scale of the order of 2 π/q . Let us first discuss the data for φ = 0 4 . (lower panel). For small active force strength, F act = 0 01 . , the system is very close to a simple fluid made of diatomic molecules. Consequently, two peaks are visible in the curve. One represents the molecule elongation, i.e. the distance between the two beads in the dumbbell, and it is located at q glyph[similarequal] 6 21 . that corresponds to glyph[lscript] glyph[similarequal] 1 . The other one signals the typical distance between beads belonging to different dumbbells and is located at q glyph[similarequal] 3 5 . , that is to say glyph[lscript] glyph[similarequal] 1 79 . . For increasing values of the active force, clustering is favoured, and molecules in them tend to be closer to each other. Therefore, the first peak (at longer distances in real space) progressively disappears while the second one (for distances of the order of glyph[lscript] glyph[similarequal] 1 ) increases its weight. However, the most important new feature in the curves is the appearance and growth of the structure factor close to vanishing q . We note that the curves intercept at q glyph[similarequal] 2 26 . and that the curves move upwards with increasing F act for q smaller than this value. This increase quantifies the presence and growth of the clusters with increasing activity. The larger values of the structure factor at small q observed for F act ≥ 0 3 . correspond to the clustering effects observed in this range of F act in Fig. 7.

Similar features in the structure factor were observed experimentally in [82] and numerically in [19, 26, 79]. In the dilute polymer melt active sample studied in [57] no such important increase of the lowq structure factor was observed and the conclusion was that active forces were making the sample more compact but uniformly, with no clustering effects. The special q values for the dumbbell system with packing fraction φ = 0 1 . are: the molecular elongation peak is here located at q glyph[similarequal] 6 67 . , the curves intercept at q glyph[similarequal] 1 48 . and the second peak at F act = 0 01 . is not really visible.

Although we have not shown it analytically, the probabil-

ity distribution functions (pdf) of the velocity components is well represented by a Gaussian pdf with a kinetic temperature given by Eq. (35) for small values of φ , as long as phase separation does not occur. These pdfs are shown with continuous and dashed lines in Fig. 10. Increasing the surface concentration, jamming effects slightly reduce the overall T kin .

In essentially out of equilibrium systems with non-potential forces that drive the dynamics, such as driven granular matter [76-78], vortex lattices [83] and the active matter we here study, the kinetic temperature should be higher than the ambient temperature. In our case we measured T kin glyph[similarequal] 0 05 . for F act = 0 01 . and F act = 0 1 . , so that T kin glyph[similarequal] T in these two cases. For the strongest F act used, F act = 1 , we measured T kin = 0 0596 . &gt; T . We can compare this value to the one expected for the single active dumbbell under the same conditions, according to Eq. (35): T single kin glyph[similarequal] 0 06 . which is very close to the value measured for the ensemble.

Figure 9. (Color online.) The structure factor of the active sample with φ = 0 1 . (upper panel) and φ = 0 4 . (lower panel) for different active force strengths given in the key. See the text for a discussion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-f8055089-2e7b-4ad1-9928-8a8b9f31f22b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-362d75df-7a82-4836-bd0d-4a2ba092bba1.jpg)

Figure 10. (Color online.) Pdf of the horizontal component of the center-of-mass dumbbell velocity, v CMx , in a system with φ = 0 1 . at T = 0 05 . , for three values of the active force, F act = 0 01 . , 0 1 . , 1 . Data are shown in linear-log scale. The values of the kinetic temperature, T kin = 0 0500 . , 0 0501 . , 0 0596 . , respectively, used in the Gaussian fit shown with continuous and dashed lines are very close to the result for a single molecule, see Eq. (35). At φ = 0 3 . , for the same set of active forces, we have checked that the pdf is also well represented by a Gaussian distribution with T kin given by Eq. (35).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-7a5a75de-12bf-46d0-aebe-121af64996b1.jpg)

## B. Mean-square displacement and asymptotic diffusion

## 1. Very low temperature

In order to appreciate the effects of the dumbbell interaction, in Fig. 5 we showed with thicker lines the meansquare displacement in the active dumbbell system with φ = 0 1 . and the same four values of the active force F act = 0 001 . , 0 01 . , 0 1 . , 1 used in the single dumbbell case. Data for single and collective systems under the same active force are shown with the same color and line style. For the two smaller applied forces, F act = 0 001 . , 0 01 . , we do not see any difference in ∆ 2 between the single and interacting case in this scale. For the two larger forces, F act = 0 1 . , 1 , the interaction between the molecules slows down the dynamics in the sense that the mean-square displacement of the interacting system, at the same time-lag, is smaller than the one for the single molecule. For the forces F act = 0 001 . , 0 01 . , 0 1 . we still see the angular time scale t a , a feature that does not exist in models of point-like active particles as the ones studied in [26, 56, 57]. Note that for F act = 1 and such a low temperature, T = 0 001 . , the system is in the phase-separated phase, the velocity-component pdf is no longer Gaussian (not shown) and the mean-square displacement is much slower than in the single molecule case.

A first presentation of the diffusion constant, D A ( F act , φ ) obtained by studying the system at the different surface fractions at a very low temperature, T = 0 001 . , and for the active forces, F act = 0 1 0 001 . , . , was given in Fig. 4. In this plot D A ( F act , φ ) was normalized by the center-of-mass diffusion

Figure 11. (Color online.) The diffusion constant in the active system at T = 0 05 . , D A ( F act , φ ) , as a function of F act for five values of the system's surface fraction given in the key. Normal scale is used. The continuous (orange) line is the analytic expression for φ = 0 , see Eq. (41). The dashed (red) line is a fit of the data to the form D A (0 0 2) + , . aF 2 act with D A (0 , 0 2) = 0 00153 . . obtained from the passive data at this temperature and a = 0 273 . (the F 2 act dependence found for the single molecule) for φ = 0 2 . (green crosses × ). The dashed (blue) line is a fit to the form D A (0 0 2) + , . aF α act with D A (0 , 0 2) = 0 00153 . . , a = 0 259 . and α glyph[similarequal] 1 56 . (the fit used in [56, 57], though α glyph[similarequal] 2 3 . was found in that case for the low density used). The corrections to this fit will be studied in the next figures and discussed in the text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-c475d990-ebe5-423a-a0ef-2209eae82df8.jpg)

constant of a single active molecule under the same conditions, D A ( F act , 0) . At F act = 0 1 . , the data fall-off from 1 at φ = 0 very quickly as soon as φ &gt; 0 . At this very low temperature and relatively high active force small clusters appear and the diffusion properties are altered by them. On the other hand, at F act = 0 001 . the data are close to those for the passive system but they lie above them. They are not well represented by the TO expression and, moreover, they give a first indication of non-monotonic dependence of the normalised data with the active force, that we will discuss in more detail below when using a higher bath temperature.

## 2. Working temperature

We will show in the following results for a higher temperature, T = 0 05 . , such that, for all the active forces considered, the system is in the homogeneous phase, even though, as showed before, non-trivial fluctuation effects are clearly observable at sufficiently high values of F act . We also studied the system at other higher temperatures obtaining results similar to those found at T = 0 05 . . (The values of the effective temperature found at T = 0 1 . will be reported in Table 1.)

The diffusion constant in the last diffusive regime of the active system is shown, as a function of active force, in Fig. 11 for five densities φ = 0 1 . , 0 2 . , 0 3 . , 0 4 . , 0 5 . . D A ( F act , φ ) was extracted from the analysis of the mean-square displacement in the late t glyph[greatermuch] t a regime. The figure demonstrates that

D A ( F act , φ ) decreases with increasing φ and increases with increasing F act . The analytic expression for the single particle limit, φ = 0 , given in Eq. (41) is added with a continuous (orange) line. Clearly, the finite density induces a reduction of the diffusion constant for all F act . A quadratic fit of the data for φ = 0 2 . , which conserves the F 2 act dependence of the single particle limit, is shown with a dashed (red) line. More precisely, the form used is D A (0 , 0 2) + . aF 2 act with D A (0 0 2) = 0 00153 , . . obtained from the passive data at this temperature and a = 0 273 . for φ = 0 2 . (green crosses × ). Had we left the F act = 0 intercept as a free parameter in this fit we would have obtained 0 0098 . which is relatively close to the actual passive limit, as it should. This fit is acceptable at small F act but a deviation at large values of the active force is clear. Another fit, with an additional fit parameter as the power in the F act dependence, is also shown with a dashed (blue) line. This fit was used in [56, 58] to describe the diffusion constant of a rather dilute ensemble of interacting active point-like particles. We find here that a power law with α glyph[similarequal] 1 56 . fits the dumbbell data at φ = 0 2 . rather correctly for all F act shown. The analysis of the φ -dependence will be done next, where a more convenient way of describing these data will be proposed.

glyph[negationslash]

In Fig. 12 we show the same diffusion constant, D A ( F act , φ ) , as a function of the surface fraction φ . The values of the active force F act are given in the key. We have already stressed, when showing the very low temperature data, that as soon as F act = 0 the rather complex TO packing fraction dependence breaks down. We therefore seek for a different (and simpler) φ -dependence of the diffusion constant of the active interacting sample. In the linear-ln presentation the data (up to φ glyph[similarequal] 0 4 . ) are rather well fitted by a straight line, suggesting

$$D _ { A } ( F _ { \text{act} }, \phi ) \simeq D _ { A } ( F _ { \text{act} }, 0 ) \, e ^ { - b ( F _ { \text{act} } ) \phi } \quad \ \ ( 4 9 )$$

with D A ( F act , 0) the single active molecule diffusion constant. (This is confirmed by a double linear presentation of data.) The exponential dependence on φ is much simpler than the TO expression and it cannot be taken as a formal proposal for the behavior of D A as we do not have an analytic justification for it. We simply stress here that it provides an acceptable description of data at this temperature.

In the second panel in Fig. 12 the finite φ diffusion constant is normalised by the single molecule one, in the form used in Fig. 4. The plot shows an unexpected non-monotonic behavior as a function of F act suggesting that the effective b F ( act ) should be non-monotonic. It is interesting to observe that the diffusion constant ratio increases at small values of F act starting to decrease when clustering effects become relevant for F act ≥ 0 3 . .

In the upper panel in Fig. 13 we display the exponential factor e -b F ( act ) φ . We use open symbols for the result of the ratio D A ( F act , φ ) /D A ( F act , 0) between numerical data, and filled symbols for the the fit of D A ( F act , φ ) /D A ( F act , 0) , as a function of F act . Three values of φ were used and are given in the key. Extracting b F ( act ) from here one confirms the nonmonotonic dependence with F act with b varying in the interval [1 , 2 7] . , circa. One observes that the fit proposed in Eq. (49)

Figure 12. (Color online.) The diffusion constant, as extracted from the mean-square displacement of the dumbbells center-of-mass, as a function of φ , in a system of active molecules at T = 0 05 . . Upper panel: A linear-ln scale is used due to the large variation in the absolute values of D A ( F act , φ ) . The values of F act are given in the key. The values at φ = 0 , shown with larger symbols, correspond to the theoretical value of D A for a single dumbbell given by Eq. (41). The data for relatively small φ are rather well fitted by a straight line, suggesting ln D A ( F act , φ ) = ln D A ( F act , 0) -b F ( act ) φ with D A ( F act , 0) the single active molecule diffusion constant. In the lower panel the data are normalised by the diffusion constant of the single molecule and they are presented in linear scale. Nonmonotonic behaviour with respect to F act is observed suggesting that the fitting function b F ( act ) should be non-monotonic.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-690f6d04-ad0c-499f-b2b9-7d2ddb02e5ea.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-1f9c5a56-6c70-43d2-b165-a403dfe7b036.jpg)

works reasonably well for all the values considered for F act . The approximate exponential fit of data is compared to the TO prediction for F act = 0 in the lower panel. The two are very close for φ glyph[lessorsimilar] 0 3 . while they deviate considerably for higher density.

To the best of our knowledge this non-monotonic behaviour has not been observed in the literature yet.

Figure 13. (Color online.) Upper panel: the exponential e -b F ( act ) φ for three densities given in the key. Original data are represented by open symbols while the filled symbols represent the values obtained from the fit in Eq. (49) of D A ( F act ,φ ) D A ( F act , 0) performed in Fig. 12. The lines are guides-to-the-eye obtained with a spline of the filled points. Lower panel: the correct TO description of data, D A (0 ,φ ) D A (0 0) , = [1 + H φ ( )] -1 for F act = 0 (with red solid line), confronted to the exponential approximation in Eq. (49), that in this case yields b F ( act = 0) = 2 4 . (with green dashed line). The two curves are very close to each other up to φ glyph[lessorsimilar] 0 3 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-ff774cc7-476d-4f23-bc45-9639b37028ff.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-a72e7c52-c802-4126-8a8e-ea8df5c78b29.jpg)

## C. The linear response function

We now turn to the study of the linear response function, integrated over a period of time, as defined in Eq. (42).

We started with an analysis of the amplitude of the applied perturbation, so as to determine the optimal value to be used for each set of parameters. Indeed, one has to use a small enough applied force to keep the response within the linear regime, but large enough to reduce the fluctuations. We found that f = 0 01 . yields the best results.

Figure 14 shows the time integrated linear response function χ t ( ) as a function of time t for various active many-body systems at T = 0 05 . and F act = 0 001 . , 0 01 . , 0 1 . , 1 . In the two panels ( φ = 0 1 . and φ = 0 3 . ) double linear scales are

Figure 14. (Color online.) The induced displacement or linear susceptibility for a system of active dumbbells at T = 0 05 . and surface fractions φ = 0 1 . (above) and φ = 0 3 . (below). The different values of the active force used are given in the key. The magnitude of the applied field is f = 0 01 . . One sees a non-monotonic dependence of µ F ( act , φ ) on F act , similarly to what we observed in D A ( F act , φ ) /D A ( F act , 0) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-6a47562e-504b-4fc3-a209-ddf374d75a83.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-b6bb1baf-6586-494b-8ae0-7d2235cce927.jpg)

used and the dependence of the linear response integrated over time on F act is now made visible. The induced displacement at a given time t decreases with the dumbbells' concentration. While the F act = 0 001 . and F act = 0 01 . curves lie on top of each other (within numerical accuracy), we notice the nonmonotonic dependence on F act for higher activities. Indeed, in both panels the curve for F act = 0 1 . lies above the ones for F act = 0 001 . , 0 01 . and F act = 1 . Moreover, the data for F act = 0 01 . and F act = 1 appear in different order for φ = 0 1 . and φ = 0 3 . . We have not included data for intermediate active forces 0 1 . &lt; F act &lt; 1 in this plot to ease the visualization but we have analysed them to extract the asymptotic slope (see the data presented in Fig. 15 below).

The non-monotonic dependence of χ F ( act , φ ) on F act for sufficiently large φ is reminiscent of - though not the same as - a negative resistivity, that is to say a decreasing dependence of a current on the applied field that drives it, observed,

Figure 15. (Color online.) Upper panel: The asymptotic slope of the integrated linear response µ F ( act , φ ) as a function of φ for various active forces given in the key. The two curves are exponential fits to the data for F act = 0 01 . , 0 1 . ; the pre factor is 1 (2 / γ ) = 0 05 . . Lower panel: the ratio µ F ( act , φ ) /µ F ( act , 0) is plotted as a function of F act for the density values given in the key. The scale on the right side gives the values of µ F ( act , φ ) . With open symbols, the ratios between numerical data; with filled symbols the result of the exponential fit of data. The curves are non-monotonic and the relative motility is enhanced for active forces of the order of 0 1 . glyph[lessorsimilar] F act glyph[lessorsimilar] 0 4 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-e996f950-6a2f-4e32-aa90-6221fb3cc47a.jpg)

for instance, in sufficiently dense kinetically constrained systems [84].

We extract the slope of these time-dependent curves for times such that t glyph[greatermuch] t a and we call it dµ F ( act , φ ) . We already know that at fixed F act , it is a decreasing function of the dumbbell concentration φ , and that at fixed density it is a non-monotonic function of F act .

The study of the dependence of µ on these two parameters is performed in Fig. 15. In the upper panel we plot µ F ( act , φ ) as a function of φ for several active forces given in the key. We confirm the monotonic decay with increasing density. The two dashed curves are exponential fits as function of φ for two choices of the active force, F act = 0 01 . , 0 1 . .

However, the dependence on F act is less straightforward. The curves for different F act cross. For instance, there is an

inverting value of φ , say φ ∗ , so that µ F ( act = 1 , φ &lt; φ ∗ ) is smaller that µ F ( act = 0 1 . , φ &lt; φ ∗ ) while one observes µ F ( act = 1 , φ &gt; φ ∗ ) &gt; µ F ( act = 0 1 . , φ &gt; φ ∗ ) . At fixed density, the relative motility is enhanced for active forces of the order of 0 1 . glyph[lessorsimilar] F act glyph[lessorsimilar] 0 4 . but the trend is reversed for higher active forces. The maximal effect is seen for φ glyph[similarequal] 0 2 . .

In the second panel we show the F act -dependence of the ratio µ F ( act , φ ) /µ F ( act , 0) for three packing fractions given in the key. The factor µ F ( act , 0) = 1 / (2 γ ) is just a constant independent of F act and φ . We use open data points to represent the ratio between the numerical data, and filled symbols with joining lines to represent the outcome of the fit to an exponential:

$$\mu ( F _ { \text{act} }, \phi ) = \mu ( F _ { \text{act} }, 0 ) \, e ^ { - c ( F _ { \text{act} } ) \, \phi } \,. \quad \ \ ( 5 0 )$$

Note that if c F ( act ) were given by b F ( act ) as for the diffusion constant, we would have a φ -independent T eff , and the same as for φ = 0 . We obtain a non-monotonic dependence on F act for φ glyph[lessorsimilar] 0 4 . while the data for the highest density, φ = 0 5 . , may approach a plateau at the value reached at F act = 0 1 . .

## D. Fluctuation-dissipation relation

In Fig. 16 we display the parametric plot ∆ ( 2 χ ) for the three choices of activity, F act = 0 01 . , 0 1 . , 1 and two surface fractions φ = 0 1 . , 0 3 . all at T = 0 05 . . We see a very weak dependence on φ and a strong dependence on F act . All data points at F act = 0 (not shown) fall on top of each other (independently of φ ) as the passive system is in equilibrium with the bath. The data for F act = 0 01 . are very close to these and the data fit gives values of the effective temperature (the exponential of the vertical axis off-set of the slope of the curves) that are near the ambient temperature as well (see the numerical estimates in Table I).

For larger values of the activity, the parametric construction displays the familiar shoulder separating short from long timedelay behaviour [44, 45]. The short time-delay fluctuations are much affected by the microscopic dynamics controlled, in this case, by the ambient temperature (the same for all sets of curves) and the activity (very different for the three sets of curves shown with different colour and line style).

After the shoulder, at long time-delays, the structural dynamics where interactions between dumbbells are important sets in. The effective temperature will be extracted from the slope of the parametric curves (in linear scale) in this timedelay regime, or from the vertical axis off-set of the slope in the late-time regime in double logarithmic scale, as explained above. That is to say, the asymptotic relation between linear integrated response and mean-square displacement allows one to extract the effective temperature as

$$k _ { B } T _ { \text{eff} } ( F _ { \text{act} }, \phi ) = \frac { D _ { A } ( F _ { \text{act} }, \phi ) } { \mu ( F _ { \text{act} }, \phi ) } \. \quad \ \ ( 5 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The effective temperature values are larger than the environmental temperature for F act &gt; 0 . The numerical dependence of T eff on F act is reported in Fig. 17, where T eff against F act is shown. Note that T eff is consistently larger than T and increases with F act . The solid line in the plot represents the theoretical result for a single active dumbbell, T eff glyph[similarequal] T + cF 2 act /T . The quadratic dependence on F 2 act is similar to what was found for interacting active point-like particles [56, 57] and interacting active polymers [57, 58] at low density, comparable to φ glyph[similarequal] 0 1 . in our case where no aggregation effects exist. The square power-law dependence on F act also applies to the dumbbell data for φ = 0 1 . . Although in this double logarithmic scale the data seem to be very close to this form for all φ , as we will see, a closer look at them shows a non-trivial φ dependence.

Figure 16. (Color online.) Parametric plot ∆ ( 2 χ ) for three values of the activity, F act = 0 01 . , 0 1 . , 1 , two values of the surface concentration φ = 0 1 . , 0 3 . , and T = 0 05 . . Double logarithmic representation. The effective temperature can be read from the off-set in the y direction of the projection of the straight line in the long-time regime after the shoulder, ln T eff = ln ∆ 2 -ln 2 χ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-ecc2a220-6bdd-45a6-98fe-96ba818552fc.jpg)

The results for T eff for two working temperatures, T = 0 05 . and T = 0 1 . , are summarised in Table I. At very small active force we find T eff glyph[similarequal] T at all densities, as expected as the system is near equilibrium (see the first rows in the two sets of data shown in Table I). For intermediate active forces, e.g. F act = 0 1 . (third rows in the tables), we see that T eff is significantly larger than T and that it weakly increases with increasing φ . For still larger active forces the dependence on φ changes, as T eff decreases with φ for sufficiently large values of φ ( φ &gt; ∼ 0 1 . for T = 0 05 . and φ &gt; ∼ 0 3 . for T = 0 1 . , when F act = 1 ). Consistently with what discussed in the previous paragraph, for fixed density, T eff increases with F act for the two working temperatures. Moreover, at fixed density and active force F act &gt; 0 5 . we see that T eff is lower at higher temperature.

In Fig. 18 we display the ratio between the finite density effective temperature and the single molecule limit one as a function of the strength of the active force, T eff ( F act , φ ) /T eff ( F act , 0) vs. F act . If we could separate the diffusion coefficient of a single dumbbell from the φ dependence as proposed in Eq. (49), and if we had the same den-

Figure 17. (Color online.) Fitted values of T eff from the parametric plot in Fig. 16 for five values of φ given in the key, all at the temperature T = 0 05 . . In this representation, all data seem to be in good agreement with the T eff formula for a single dumbbell, Eq. (47), that is represented with the solid line. However, the actual values of T eff are given in Table I where one sees that there is, though, a weak difference in the data for different φ : for F act = 0 1 . , T eff tends to increase with increasing φ , while for F act = 1 , T eff tends to decrease with increasing φ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-c7c98c00-ae0b-4105-ab39-978fa2fd533a.jpg)

Table I. Values of T eff from the analysis at different active forces and system's densities, at T = 0 05 . and T = 0 1 . .

| T = 0 . 05   | T = 0 . 05   | T = 0 . 05   | T = 0 . 05   | T = 0 . 05   |
|--------------|--------------|--------------|--------------|--------------|
| F act        | φ = 0        | φ = 0 . 1    | φ = 0 . 3    | φ = 0 . 5    |
| 0.001        | 0.0500       | 0.0499       | 0.0488       | 0.0502       |
| 0.01         | 0.0509       | 0.0498       | 0.0502       | 0.0522       |
| 0.1          | 0.142        | 0.152        | 0.179        | 0.223        |
| 0.5          | 2.35         | 2.51         | 2.59         | 2.44         |
| 1            | 9.27         | 9.68         | 7.83         | 6.56         |
| T = 0 . 1    | T = 0 . 1    | T = 0 . 1    | T = 0 . 1    | T = 0 . 1    |
| F act        | φ = 0        | φ = 0 . 1    | φ = 0 . 3    | φ = 0 . 5    |
| 0.001        | 0.100        | 0.102        | 0.100        | 0.097        |
| 0.01         | 0.100        | 0.100        | 0.100        | 0.101        |
| 0.1          | 0.146        | 0.158        | 0.167        | 0.201        |
| 0.5          | 1.25         | 1.34         | 1.56         | 1.77         |
| 1            | 4.71         | 5.03         | 5.10         | 4.92         |

sity dependence of the µ F ( act , φ ) as in the remaining factor in D A ( F act , φ ) , meaning b F ( act ) = c F ( act ) , then T eff ( F act , φ ) should be independent of φ and just equal to T eff ( F act , 0) . The data show that this holds for φ = 0 1 . where the data points are (within numerical accuracy) constant but it does not for higher densities. For φ ≥ 0 1 . one observes a non-monotonic dependence of T eff ( F act , φ ) /T eff ( F act , 0) on F act , reflecting in some way the non-monotonic behavior first shown for the diffusion and response ratios D A ( F act , φ ) /D A ( F act , 0) and µ F ( act , φ ) /µ F ( act , 0) . The finite density effective temperature is higher than the single molecule one for active forces in the interval [0 , 0 5] . while it is lower than the single molecule one for active forces in the interval [0 5 . , 1] . A maximum is reached at around F act glyph[similarequal] 0 1 . . We ascribe the change in be- haviour to the presence of finite size clusters in the sample for F act &gt; ∼ 0 5 . and sufficiently large density. These clusters (see the last panels in Fig. 7) would respond differently from the homogeneous bulk and their dynamics will be different. We also note that these curves seem to cross at the active force strength value F act = 0 5 . . A similar analysis of this ratio at a different temperature, T = 0 1 . , shows that the crossing occurs at a different value of the active force strength, F act glyph[similarequal] 1 , that corresponds to the same P´ eclet number Pe = 20 (within numerical accuracy).

Figure 18. (Color online.) The ratio between the finite density effective temperature and the single molecule limit one as a function of the strength of the active force, T eff ( F act , φ ) /T eff ( F act , 0) , where T eff ( F act , 0) = T [1 + Pe 2 / 8] . The temperature is T = 0 05 . . See the text for a discussion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0638v2-figure-c4358931-2a5b-4c03-a862-49eb05119dd3.jpg)

Finally, we note that T eff &gt; T kin as soon as the active force is intense enough to see a considerable deviation of both from the ambient temperature. Confront, for instance, the analytic values T eff = 0 142 . and T kin glyph[similarequal] 0 050 . at φ = 0 and F act = 0 1 . , and the numeric values T eff glyph[similarequal] 9 68 . and T kin glyph[similarequal] 0 06 . at φ = 0 1 . and F act = 1 .

## VII. CONCLUSIONS

In this paper we studied the dynamic properties of isolated and interacting, passive and active, dumbbell systems.

The model studied here differs from other models of active matter previously studied in the literature with numerical techniques from a similar viewpoint. For instance, we do not impose polar alignment mechanisms as in [36, 85]. The active forces used here are different from the ones used in the analysis of active macromolecules presented in [57, 58]: while in the dumbbell system the active forces act along the main axis of the molecule, in the polymer system the active forces act only on the center monomer and in random directions. The Langevin dynamics used in this paper is different from the Monte Carlo rule chosen in [26] and from the run-and-tumble bacteria models in [18, 54]. Moreover, we did not distinguish

translation and rotational noise, as done in [19], but we simply added independent Gaussian white noise (related to the friction dissipative term in the usual way) to the Langevin equations for the positions of the two atoms in the dumbbell.

For fixed temperature and active force strength, this very simple model has a phase transition between a homogeneous and a phase separated phase [24, 25]. In this paper we focused on the low density phase in which the system is homogeneous on average with, possibly, giant density fluctuations for sufficiently high density [25]. We studied the dynamics for three different temperatures and a wide range of densities and active forces in this phase.

We first presented a detailed study of the mean-square displacement of the interacting active sample. We analysed the single passive and active dumbbell dynamics and we investigated how the finite density affects the various dynamic regimes and, in particular, the diffusive properties in the late time-delay limit that goes beyond the angular diffusion time of the single molecule. As one had expected, we found that the diffusion constant decreases with increasing density and increases with increasing active force. The ratio between the finite-density and the single particle diffusion constant exhibits an intriguing non-monotonic dependence on the active force. We also analysed how the Tokuyama-Oppenheim density-dependence expression for this ratio [75] is modified by self-propulsion, finding that a simpler exponential decay (with a non-monotonic active force dependent factor in the exponential) fits the data reasonably well.

We then studied the linear response function of the dumbbell displacement to infinitesimal perturbations that push them in random directions. As usual, to minimise the numerical error, we focused on the linear response integrated over time (instead of the instantaneous response that fluctuates much more). We found that, at fixed time delay, the integrated linear response decreases monotonically with increasing density. The dependence on the active force is again non-trivial, with a maximum response reached for active forces with strength in 0 1 . glyph[lessorsimilar] F act glyph[lessorsimilar] 0 3 . for densities in 0 1 . glyph[lessorsimilar] φ glyph[lessorsimilar] 0 3 . .

glyph[negationslash]

The kinetic temperature is extracted from the asymptotic value of a one-time observable, the kinetic energy. As already discussed in detail in the context of glassy systems and granular matter [44, 45], although one proves that a system is not in equilibrium with the thermal bath whenever T kin = T , the reverse is not true. In glassy systems in relaxation T kin = T while these systems evolve out of equilibrium. The reason is that the kinetic temperature carries information about a fast observable, the kinetic energy, that can be able to quickly equilibrate with the environment, while other observables, being slower, can still be far from equilibrium. Indeed, the kinetic temperature does not characterise the large scale structural relaxation of glassy or driven systems.

The effective temperature notion, as obtained from the deviation from the equilibrium fluctuation-dissipation theorem linking spontaneous and induced fluctuations, has proven to be very useful to understand the dynamics of slowly relaxing passive systems, such as glasses and gently sheared supercooled liquids [45]. This concept has been explored in some active systems as well as we explained in the introduction. In this work we studied the fluctuation-dissipation ratio for the low-density active dumbbell system and we characterised it as a function of density and activity. The effective temperature is always higher than the ambient temperature, it increases with increasing activity and, for small active force it monotonically increases with density while for sufficiently high activity it first increases to next decrease with the packing fraction. This effect should be due to the existence of finite-size clusters for sufficiently high activity and density at the fixed (low) temperatures at which we worked. The crossover occurs at lower activity or density the lower the external temperature. The finite density effective temperature is higher (lower) than the single dumbbell one below (above) a cross-over value of the P´ eclet number.

In the active dumbbell system we measured T kin values that are very close to T for F act = 0 01 . , 0 1 . , 1 . The existence of the kinetic temperature characterising the short-time delay behaviour of fluctuations, and the effective temperature extracted from the fluctuation-dissipation relations at long timedelays, does not invalidate the possible thermodynamic interpretation of the latter. One simply has to focus on the dynamics of the systems in one or another dynamic regime [44, 45] and test its thermodynamic properties within in it.

We want to stress once again that the effective temperature is not a parameter that characterises the statistical properties of the activity but an intensive parameter that tells us about the dynamic properties of the interacting and active many-body system.

Fily and Marchetti [19] argue that the effective tempeature notion cannot apply to active matter systems in their dense phase. They base their claim on the comparison of instantaneous snapshots of typical clustered configurations and equivalent thermal equilibrium ones with similar overlap between particles at the same packing fraction. This argument cannot be used to refute the effective temperature ideas in this context (nor in the glassy context either) as its very definition is dynamic and correlation and linear response functions at different times need to be calculated to derive T eff . Having said this, the analysis of T eff in a dense active phase has not been performed yet and we cannot claim that the concept will still be valid in this clustered phase.

The dumbbell system is particularly interesting as it allows one to study translational and, simultaneously, also rotational degrees of freedom. A recent analysis of the effective temperature ideas in a passive dumbbell system has shown that the relation between the two can depend non-trivially on the elongation of the molecule [86]. It would be very interesting to explore the effect of this parameter in the results that we showed and to investigate the fluctuation-dissipation relations of rotational degrees of freedom in the active dumbbell system as well.

The values of the effective temperatures found from the fluctuation-dissipation relation of the active system should be confirmed with alternative measurements to test their thermodynamic meaning. For instance, one could use spherical passive particles as tracers and measure the effective temperature from their diffusion properties (as in [46]) and kinetic energy fluctuations. We will report on the use of tracers in this con-

text in a separate publication.

Finally, it would be interesting to analyse the effect that external potentials may have on the dynamics of this active dumbbell system. Tailleur and Cates derived a diffusive approximation for active run-and-tumble particles [18, 54] and found that, for weak external potentials that modify only slightly the velocity of the particles, the system should follow equilibrium dynamics with the equilibrium bath temperature replaced by T eff . It will be interesting to check whether this result holds for the active dumbbell system under appropriate external forces and how it is modified by the finite density and

- [1] J. Toner, Y. Tu, and S. Ramaswamy, Ann. of Phys. 318 , 170 (2005).
- [2] D. A. Fletcher and P. L. Geissler, Ann. Rev. Phys. Chem. 60 , 469 (2009).
- [3] G. I. Menon, Rheology of Complex Fluids p. 193 (2010).
- [4] S. Ramaswamy, Ann. Rev. Cond. Matt. Phys. 1 , 323 (2010).
- [5] M. E. Cates, Rep. Prog. Phys. 75 , 042601 (2012).
- [6] T. Vicsek and A. Zafeiris, Phys. Rep. 517 , 71 (2012).
- [7] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Rev. Mod. Phys. 85 , 1143 (2013).
- [8] U. S. Schwarz and S. A. Safran, Rev. Mod. Phys. 85 , 1327 (2013).
- [9] G. de Magistris and D. Marenduzzo, to appear in Physica A (2014).
- [10] V. Narayan, S. Ramaswamy, and N. Menon, Science 317 , 105 (2007).
- [11] J. Deseigne, O. Dauchot, and H. Chat´ e, Phys. Rev. Lett. 105 , 098001 (2010).
- [12] W. F. Paxton, K. C. Kistler, C. C. Olmeda, A. Sen, S. K. S. Angelo, Y. Cao, T. E. Mallouk, P. E. Lammert, and V. H. Crespi, J. Am. Chem. Soc. 126 , 13424 (2004).
- [13] L. Hong, S. Jiang, and S. Granick, Langmuir 22 , 9495 (2006).
- [14] J. Palacci, C. Cottin-Bizonne, C. Ybert, and L. Bocquet, Phys. Rev. Lett. 105 , 088304 (2010).
- [15] P. Romanczuk, M. B¨ ar, W. Ebeling, B. Lindner, and L. Schimansky-Geier, Eur. Phys. J. Special topics 202 , 1 (2012).
- [16] A. Czir´ ok, H. E. Stanley, and T. Vicsek, J. Phys. A 30 , 1375 (1997).
- [17] T. Vicsek, A. Czir´k, E. Ben-Jacob, I. Cohen, and O. Shochet, o Phys. Rev. Lett. 75 , 1226 (1995).
- [18] J. Tailleur and M. E. Cates, Phys. Rev. Lett. 100 , 218103 (2008).
- [19] Y. Fily and M. C. Marchetti, Phys. Rev. Lett. 108 , 235702 (2012).
- [20] Y. Fily, S. Henkes, and M. C. Marchetti, Soft Matter 10 , 2132 (2014).
- [21] G. S. Redner, M. F. Hagan, and A. Baskaran, Phys. Rev. Lett. 110 , 055701 (2013).
- [22] I. Buttinoni, J. Bialk´ e, F. K¨ ummel, H. L¨ owen, C. Bechinger, and T. Speck, Phys. Rev. Lett. 110 , 238301 (2013).
- [23] J. Stenhammar, A. Tiribocchi, R. Allen, D. Marenduzzo, and M. Cates, Phys. Rev. Lett. 111 , 145702 (2013).
- [24] G. Gonnella, A. Lamura, and A. Suma, Int. J. Mod. Phys. C (2013).
- [25] A. Suma, D. Marenduzzo, G. Gonnella, and E. Orlandini, submitted to EPL (2014).
- [26] D. Levis and L. Berthier, Phys. Rev. E 89 , 062301 (2014).

interactions within the sample, and strong external potential forces.

Acknowledgments: We thank D. Loi and S. Mossa for very useful discussions. We warmly thank J. Tailleur for very useful comments on the manuscript. LFC is a member of Institut Universitaire de France. GG acknowledges the support of MIUR (project PRIN 2012NNRKAF).

- [27] R. Wittkowski, A. Tiribocchi, J. Stenhammar, R. Allen, D. Marenduzzo, and M. Cates, Nat. Comm. 5 , 4351 (2014).
- [28] S. Ramaswamy, R. A. Simha, and J. Toner, Europh. Lett. 62 , 196 (2003).
- [29] F. Ginelli, F. Peruani, M. B¨ ar, and H. Chat´ e, Phys. Rev. Lett. 104 , 184502 (2010).
- [30] J. Elgeti and G. Gompper, EPL 85 , 38002 (2009).
- [31] J. Elgeti and G. Gompper, EPL 101 , 48003 (2013).
- [32] A. Bricard, J.-B. Caussin, N. Desreumaux, O. Dauchot, and D. Bartolo, Nature 503 , 95 (2013).
- [33] S. Wang and P. G. Wolynes, Proc. Nac. Acad. Sc. USA 108 , 15184 (2011).
- [34] T. Shen and P. G. Wolynes, Proc. Nac. Acad. Sc. USA 101 , 8547 (2004).
- [35] T. E. Angelini, E. Hannezo, X. Trepat, M. Marquez, J. J. Fredberg, and D. A. Weitz, Proc. Nat. Acad. Sc. USA 108 , 4714 (2011).
- [36] S. Henkes, Y. Fily, and M. C. Marchetti, Phys. Rev. E 84 , 040301 (R) (2011).
- [37] L. Berthier and J. Kurchan, Nature Phys. 9 , 310 (2013).
- [38] H. Chat´ e, F. Ginelli, and R. Montagne, Phys. Rev. Lett. 96 , 180602 (2006).
- [39] M. E. Cates, S. M. Fielding, D. Marenduzzo, E. Orlandini, and J. M. Yeomans, Phys. Rev. Lett. 101 , 068102 (2008).
- [40] G. P. Saracco, G. Gonnella, D. Marenduzzo, and E. Orlandini, Phys. Rev. E 80 , 051126 (2009).
- [41] M. A. Bees and O. A. Croze, Proc. Roy. Soc. A: Math. Phys. and Eng. Sci. 466 , 2057 (2010).
- [42] R. N. Bearon, M. A. Bees, and O. A. Croze, Phys. Fluids 24 , 121902 (2012).
- [43] G. P. Saracco, G. Gonnella, D. Marenduzzo, and E. Orlandini, Central European Journal of Physics 10 , 1109 (2012).
- [44] L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. E 55 , 3898 (1997).
- [45] L. F. Cugliandolo, J. Phys. A: Math. and Theor. 44 , 483001 (2011).
- [46] X.-L. Wu and A. Libchaber, Phys. Rev. Lett. 84 , 3017 (2000).
- [47] T. Shen and P. G. Wolynes, Phys. Rev. E 72 , 041927 (2005).
- [48] S. Wang and P. G. Wolynes, J. Chem. Phys. 135 , 051101 (2011).
- [49] P. Martin, A. J. Hudspeth, and F. J¨licher, Proc. Nac. Acad. Sc. u USA 98 , 14380 (2001).
- [50] D. Mizuno, C. Tardin, C. F. Schmidt, and F. C. MacKintosh, Science 315 , 370 (2007).
- [51] E. Ben-Isaac, Y.-K. Park, G. Popescu, F. L. Brown, N. S. Gov, and Y. Shokef, Phys. Rev. Lett. 106 , 238103 (2011).
- [52] P. Bohec, F. Gallet, C. Maes, S. Safaverdi, P. Visco, and F. van Wijland, EPL 102 , 50005 (2013).
- [53] E. Fodor, M. Guo, N. S. Gov, P. Visco, D. A. Weitz, and F. van

- Wijland (2014).
- [54] J. Tailleur and M. E. Cates, EPL 86 , 60002 (2009).
- [55] G. Szamel, Phys. Rev. E 90 , 012111 (2014).
- [56] D. Loi, S. Mossa, and L. F. Cugliandolo, Phys. Rev. E 77 , 051111 (2008).
- [57] D. Loi, S. Mossa, and L. F. Cugliandolo, Soft Matter 7 , 10193 (2011).
- [58] D. Loi, S. Mossa, and L. F. Cugliandolo, Soft Matter 7 , 3726 (2011).
- [59] F. Peruani, A. Deutsch, and M. B¨ ar, Phys. Rev. E 74 , 030904(R) (2006).
- [60] Y. Yang, V. Marceau, and G. Gompper, Phys. Rev. E 82 , 031904 (2010).
- [61] F. Peruani, J. Starruss, V. Jakovljevic, L. Sogaard-Andersen, A. Deutsch, and M. B¨ ar, Phys. Rev. Lett. 108 , 098102 (2012).
- [62] S. R. McCandlish, A. Baskaran, and M. F. Hagan, Soft Matter 8 , 2527 (2012).
- [63] H. H. Wensink, J. Dunkel, S. Heindereich, K. Drescher, R. Goldstein, H. Lowen, and J. Yeomans, Proc. Nat. Acad. Sc. USA 109 , 14308 (2012).
- [64] H. H. Wensink, H. Lowen, M. Marechal, A. Hartel, R. Wittkowski, U. Zimmermann, A. Kaiser, and A. M. Menzel, Eur. Phys. J. - Sp. Topics 222 , 3023 (2013).
- [65] A. V´ asquez-Quesada, M. Ellero, and P. Espa˜ol, Phys. Rev. E n 79 , 056707 (2009).
- [66] B. Kowalik and R. Winkler, J. Chem. Phys. 138 , 104903 (2013).
- [67] G. Alexander and J. Yeomans, Europh. Lett. 83 , 34006 (2008).
- [68] V. B. Putz, J. Dunkel, and J. Yeomans, Chemical Physics 375 , 557 (2010).
- [69] V. Gyrya, I. S. Aranson, L. Berlyand, and D. Karpeev, Bull. Math Biol. Chem. Phys. 72 , 148 (2010).
- [70] C. Valeriani, M. Li, J. Novosel, J. Arlt, and D. Marenduzzo, Soft Matter 7 , 5228 (2011).
- [71] J. D. Weeks, D. Chandler, and H. C. Andersen, J. Chem. Phys. 54 , 5237 (1971).
- [72] E. Vanden-Eijnden and G. Ciccotti, Chem. Phys. Lett. 429 , 310 (2006).
- [73] L. F. Cugliandolo, J. Kurchan, and G. Parisi, J. Phys. I (France) 4 , 1641 (1994).
- [74] J. K. Dhont, An introduction to dynamics of colloids , vol. 2 (Elsevier Science, 1996).
- [75] M. Tokuyama and I. Oppenheim, Phys. Rev. E 50 , 16 (1994).
- [76] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. Mod. Phys. 68 , 1259 (1996).
- [77] I. S. Aranson and L. S. Tsimring, Rev. Mod. Phys. 78 , 641 (2006).
- [78] O. Pouliquen and Y. Forterre, Annu. Rev. Fluid Mech. 40 , 1 (2008).
- [79] B. ten Hagen, S. van Teeffelen, and H. L¨ owen, J. Phys.: Condens. Matter 23 , 194119 (2011).
- [80] R. Howse, R. A. L. Jones, A. J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, Phys. Rev. Lett. 99 , 048102 (2007).
- [81] N. Darnton, L. Turner, K. Breuer, and H. C. Berg, Biophys. J. 86 , 1863 (2004).
- [82] I. Theurkau, C. Cottin-Bizonne, and J. Palacci, Phys. Rev. Lett. 108 , 268303 (2012).
- [83] A. B. Kolton, R. Exartier, L. F. Cugliandolo, D. Dom´ ınguez, and N. Grønbech-Jensen, Phys. Rev. Lett. 89 , 227001 (2002).
- [84] M. Sellitto, Phys. Rev. Lett. 101 , 048301 (2008).
- [85] B. Szab´ o, G. J. Sz¨ll¨si, B. G¨nci, Z. Jur´nyi, D. Selmeczi, and o o o a T. Vicsek, Phys. Rev. E 74 , 061908 (2006).
- [86] A. S. Ninarello, N. Gnan, and F. Sciortino (2014), arXiv:1406.6226.
Antonio Suma Giuseppe Gonnella Gianluca Laghezza Antonio Lamura Alessandro Mossa Leticia F. Cugliandolo
Dynamics of a homogeneous active dumbbell system
We analyse the dynamics of a two dimensional system of interacting active
dumbbells. We characterise the mean-square displacement, linear response
function and deviation from the equilibrium fluctuation-dissipation theorem as
a function of activity strength, packing fraction and temperature for
parameters such that the system is in its homogeneous phase. While the
diffusion constant in the last diffusive regime naturally increases with
activity and decreases with packing fraction, we exhibit an intriguing
non-monotonic dependence on the activity of the ratio between the finite
density and the single particle diffusion constants. At fixed packing fraction,
the time-integrated linear response function depends non-monotonically on
activity strength. The effective temperature extracted from the ratio between
the integrated linear response and the mean-square displacement in the last
diffusive regime is always higher than the ambient temperature, increases with
increasing activity and, for small active force it monotonically increases with
density while for sufficiently high activity it first increases to next
decrease with the packing fraction. We ascribe this peculiar effect to the
existence of finite-size clusters for sufficiently high activity and density at
the fixed (low) temperatures at which we worked. The crossover occurs at lower
activity or density the lower the external temperature. The finite density
effective temperature is higher (lower) than the single dumbbell one below
(above) a cross-over value of the Peclet number.