Active Matter, Microreversibility, and Thermodynamics

Active Matter, Microreversibility, and Thermodynamics

PDF

Pierre Gaspard¹, Raymond Kapral²

Research. Vol 2020 Article ID 9739231

Less

Research. Vol 2020 Article ID 9739231

• Research Article •

Active Matter, Microreversibility, and Thermodynamics

Full

Pierre Gaspard¹, Raymond Kapral²

Affiliations

¹Center for Nonlinear Phenomena and Complex Systems, Université Libre de Bruxelles (U.L.B.), Code Postal 231, Campus Plaine, B-1050 Brussels, Belgium

²Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada, M5S 3H6

Published: 2020-05-20 doi: 10.34133/2020/9739231

Outline

Abstract

Less

Active matter, comprising many active agents interacting and moving in fluids or more complex environments, is a commonly occurring state of matter in biological and physical systems. By its very nature, active matter systems exist in nonequilibrium states. In this paper, the active agents are small Janus colloidal particles that use chemical energy provided by chemical reactions occurring on their surfaces for propulsion through a diffusiophoretic mechanism. As a result of interactions among these colloids, either directly or through fluid velocity and concentration fields, they may act collectively to form structures such as dynamic clusters. A general nonequilibrium thermodynamics framework for the description of such systems is presented that accounts for both self-diffusiophoresis and diffusiophoresis due to external concentration gradients, and is consistent with microreversibility. It predicts the existence of a reciprocal effect of diffusiophoresis back onto the reaction rate for the entire collection of colloids in the system, as well as the existence of a clustering instability that leads to nonequilibrium inhomogeneous system states.

Cite this Article

Pierre Gaspard, Raymond Kapral. Active Matter, Microreversibility, and Thermodynamics[J]. Research, 2020 , 2020 (5) : 9739231 - . DOI: 10.34133/2020/9739231

Full Text

Less

1. 　Introduction

Less

Active matter is composed of motile entities or agents interacting with each other either directly or through the velocity and concentration fields of the medium in which they move. Such interactions lead to collective dynamics giving rise to states of matter that may differ from those in equilibrium systems. The study of such collective behavior presents challenges and is currently a topic of considerable scientific interest. Systems with many complex agents can be investigated in different ways. One way is to describe collective dynamics at the macroscale in terms of fields representing the distribution of the agents across the system. These fields are ruled by partial differential equations that are established using general symmetries and experimental observations. Another approach is to model active matter as being composed of active particles moving in space according to specific rules that are postulated on the basis of empirical considerations.

Both of these approaches have been used to explore the origins and types of collective dynamics that can be found in active matter systems, and research on this topic ranges from studies of simple active particle models, often satisfying minimal rules, to suspensions of more complex active synthetic or biological agents [1–11]. The collective behavior in systems where the active agents are chemically propelled colloids, the subject of this paper, has also been the topic of experimental and theoretical research [12–25].

Systems containing colloidal particles are governed by physicochemical laws, so that their time evolution can be understood from first principles using statistical-mechanical methods. This approach was pioneered by Einstein [26] and Smoluchowski [27–29] at the beginning of the twentieth century and systematically developed since then for passive colloidal particles [30–33]. In active matter, the colloidal particles are propelled with energy supplied by the surrounding solution, so that the description should be extended to include the molecular concentrations of fuel and product powering their motion, in addition to the velocity field of the fluid. Through such an approach, active matter can be described from the scale of a single colloidal motor moving in the surrounding fluid, up to the macroscale where many colloidal motors generate collective motion by interaction. At the macroscale, collective dynamics is described in terms of the distribution function giving the orientation as well as the position of the colloidal motors. This statistical-mechanical approach has the advantage that the parameters characterizing active matter at the macroscale can be deduced from the microscopic level of description. The knowledge of these parameters in terms of the properties of materials composing the colloidal motors and the surrounding solution is fundamental for engineering active systems.

The present paper contributes to the statistical-mechanical and nonequilibrium thermodynamic approaches for active matter systems [34–42], and considers systems whose active agents are Janus colloids with catalytic and noncatalytic faces moving by diffusiophoresis generated by chemical reactions taking place on their catalytic faces or caps [43, 44, 40]. Because of diffusiophoresis, the velocity and concentration fields are coupled together in the fluid around the Janus particle [45]. We start from the calculation of the diffusiophoretic force and torque on a single Janus particle moving in a fluid in the presence of molecular species corresponding to the fuel and the product of the reaction taking place on its catalytic surface. The concentrations of these molecular species may develop gradients on large scales under nonequilibrium conditions, and these gradients should be included in the calculation of the force and torque. The resulting diffusiophoretic force and torque enter the coupled Langevin equations ruling the displacement, rotation, and overall reaction of a single active particle.

Next, the evolution equation is established for the distribution function of the ensemble of active particles in a dilute colloidal solution. In order to be consistent with microreversibility, the principles of nonequilibrium thermodynamics are used to relate the thermodynamic forces or affinities to the current densities with linear response coefficients satisfying Onsager’s reciprocal relations [46–53]. This method allows us to obtain all the possible couplings compatible with microreversibility, including a priori unexpected reciprocal effects. Moreover, this method provides an expression for the entropy production rate density for active matter in agreement with the second law of thermodynamics and including the contribution of the reaction powering activity. Through this procedure, macroscopic evolution equations are obtained that govern the collective dynamics of colloidal motors coupled to the molecular concentrations of fuel and product. These equations can be shown to generate the reciprocal effect of diffusiophoresis back onto the reaction rate that has been obtained previously for a single particle [39, 40], but now at the macroscale. Furthermore, pattern formation due to a clustering instability manifests itself under nonequilibrium conditions induced by a bulk reaction replenishing the solution with fuel.

The paper is organized as follows. Section 2 is devoted to the dynamics of a single colloidal motor. The force and torque due to diffusiophoresis are deduced by solving the diffusion equations for the molecular concentrations coupled to the Navier-Stokes equations for the fluid velocity, including the contributions of concentration gradients at large distances from the particle. These contributions were neglected previously [39, 40] and are calculated in detail here. In Section 3, the diffusiophoretic force and torque obtained in Section 2 are incorporated into the evolution equation for the distribution function describing the ensemble of colloidal motors, and the entropy production rate density is explicitly obtained. Two implications of these results are presented in Sections 4 and 5. First, the reciprocal effect due to the diffusiophoretic coupling of an external force and torque back onto the reaction rate is recovered, now at the level of the collective dynamics. Second, a clustering instability leading to pattern formation is shown to manifest itself. The conclusions of the research are given in Section 6. The appendices provide additional details of the calculations.

2. 　Diffusiophoresis and Colloidal Motors

Less

This section describes the motion of a single spherical Janus colloidal motor of radius

R

that is propelled by self-diffusiophoresis generated by a reversible reaction,

(1) A + C ⇄ κ − κ + B + C,

with rate constants

κ ±

taking place on its catalytic surface, as depicted in Figure 1. In this reaction,

A

is the fuel and

B

the product, which are present in the solution surrounding the particle. Moreover, the concentrations of the

A

and

B

molecular species are assumed to have gradients

g k

with

k = A, B

at large distances from the particle that also contribute to motion by diffusiophoresis; thus, the motion of the particle is determined by processes in the fluid surrounding the particle.

2.1. 　Chemohydrodynamics around a Colloidal Motor

In order to determine the force and the torque due to diffusiophoresis, as well as the overall reaction rate, the velocity of the fluid and the concentrations of fuel

A

and product

B

should be obtained by solving the Navier-Stokes equations for the fluid velocity

v = v fluid

coupled to the advection-diffusion equations for the molecular concentrations

c k

with

k = A, B

(2) ρ ∂ t v + v · ∇ v = − ∇ p + η ∇ 2 v, (3) ∇ · v = 0, (4) ∂ t c k + v · ∇ c k = D k ∇ 2 c k,

where

ρ

is the constant mass density (the fluid being assumed to be incompressible),

p

the hydrostatic pressure,

η

the shear viscosity, and

D k

the molecular diffusivity of species

k

The coupling between the velocity and concentration fields is established with the boundary conditions [40, 54]

(5) n · v − v solid R = 0, (6) 1 ⊥ · v − v solid R = 1 ⊥ · b ∇ v S − ∑ k b k ∇ c k R, (7) D k n · ∇ c k R = − ν k κ + c A − κ − c B R,

where

n

is the unit vector normal to the solid surface,

1 ⊥ ≡ 1 − n n

b

is the slip length,

∇ v S = ∇ v + ∇ v T

T

denotes the transpose,

b k

is the diffusiophoretic coefficient of species

k

coupling the velocity field to the corresponding concentration field because of different interactions between the solid surface with the molecules of different species. The velocity field inside the solid particle is given by

v solid = V + Ω × r − R

in terms of the translational and angular velocities of the particle, respectively, denoted by

V

and

Ω

, and position

R

of the center of mass of the particle. Equations (5) and (6) are the boundary conditions on the components of the velocity field that are, respectively, normal and tangential to the interface

Σ t

, which is located on the sphere

∥ r − R ∥ = R

. The last equation, i.e., equation (7) is the boundary condition for the two reacting species

k = A, B

, where

ν k

is the stoichiometric coefficient of species

k

in the reaction (

ν A = − 1

and

ν B = + 1

), and

κ ±

are the forward and reverse surface rate constants per unit area.

The velocity field is assumed to vanish at large distances from the particle, so that the fluid is at rest except in the vicinity of the colloid. With the aim of obtaining mean-field equations for a dilute suspension of active particles, we also assume that the concentration fields can have nonvanishing gradients on large spatial scales. Accordingly, the concentration gradients

∇ c k ∞ = g k

are taken to exist at large distances from the colloidal particle.

We suppose that the diffusiophoretic coefficients take the values

b k c

and

b k n

on the catalytic and noncatalytic hemispheres, respectively, while the surface rate constants per unit area take positive values

κ ± c

on the catalytic hemisphere and vanish on the noncatalytic hemisphere,

κ ± n = 0

. Using spherical coordinates

θ, φ

with polar angle

θ

defined with respect to the axis of the cylindrical symmetry of the Janus particle, we have

(8) b k θ, φ = ∑ h = c, n b k h H h θ, κ ± θ, φ = ∑ h = c, n κ ± h H h θ,

where

H h θ

denotes the Heaviside function such that

H h θ = 1

on hemisphere

h

and is zero otherwise. The catalytic hemisphere is taken as

0 ≤ θ ≤ π / 2

, and the noncatalytic hemisphere as

π / 2 < θ ≤ π

Solving equations (2)–(4) with the boundary conditions (5)–(7), the velocity and concentration fields can be obtained in the vicinity of the colloidal motor [45, 40, 41]. Accordingly, the force and the torque exerted by the fluid on the motor, as well as the overall reaction rate at its catalytic surface, are given by the following surface integrals at the fluid-colloid interface

Σ t

(9) F = − ∫ ∑ t P · n d Σ, T = − ∫ ∑ t r − R × P · n d Σ, W = ∫ ∑ t κ + c A − κ − c B d Σ,

where

P = p 1 − η ∇ v S

is the pressure tensor of the fluid. The fluctuating contributions from thermal noise can also be included [40, 41].

2.2. 　Coupled Langevin Equations for the Motor

The orientation of the Janus particle is described by the unit vector

u

attached to the axis of the cylindrical symmetry of the Janus particle and pointing towards the catalytic hemisphere. Accordingly, the displacement and rotation of the particle are ruled by

(10) d R d t = V, d u d t = Ω × u,

in terms of the translational and rotational velocities. These velocities, as well as the number

N

of reactive events taking place on the particle, are governed by the following coupled Langevin equations [39, 40, 41]:

(11) M d V d t = − γ t V + F d + F ext + F f l t, (12) I · d Ω d t = − γ r Ω + T d + T ext + T f l t, (13) d N d t = W r x n + W d + W f l t,

where

M

and

I

denote the mass and inertia tensor of the motor,

γ t = 6 π η R 1 + 2 b / R / 1 + 3 b / R

is the translational friction coefficient,

γ r = 8 π η R 3 / 1 + 3 b / R

the rotational friction coefficient,

F d

and

T d

the diffusiophoretic force and torque,

F ext

and

T ext

the external force and torque exerted on the particle, while

F f l t

and

T f l t

are the contributions to the force and torque due to thermal fluctuations. The overall net reaction rate is

W r x n

W d

is the reciprocal contribution of diffusiophoresis back onto the reaction rate, and

W f l t

is the fluctuating reaction rate. If the Janus particle has a magnetic dipole

μ

and is subjected to an external magnetic field

B

, then the external torque would be given by

T ext = μ u × B

. In the overdamped regime, the coupled Langevin equations are obtained by neglecting the inertial terms in equations (11) and (12).

Solving the Navier-Stokes equations (2) and (3) coupled to equation (4) for molecular concentrations with the boundary conditions (5)–(7), the force and the torque exerted on a spherical particle of radius

R

in a fluid with shear viscosity

η

and the overall net reaction rate are given by [40]

(14) F d = 6 π η R 1 + 3 b / R ∑ k b k 1 ⊥ · ∇ c k ¯ s, (15) T d = 12 π η R 1 + 3 b / R ∑ k b k r × ∇ c k ¯ s, (16) W r x n = 4 π R 2 κ + c A − κ − c B ¯ s,

expressed in terms of the surface average

(17) · ¯ s = 1 4 π ∫ · r = R d cos θ d φ .

The expressions (14) and (15) find their origin in the generalization of Faxén’s theorem to a sphere moving in a time-dependent velocity field [55, 56]. When writing these equations, we have taken into account the possibility that the diffusiophoretic coefficients

b k

and the surface rate constants

κ ±

may be nonuniform on the particle surface.

2.3. 　Motion in Molecular Concentration Gradients

If molecular diffusion is fast enough so that the concentration fields adopt stationary profiles around the catalytic particle in the concentration gradients

g k

, the diffusiophoretic translational and rotational velocities can be written as follows (see Appendix A):

(18) V d = F d γ t = V s d u + ∑ k ξ k 1 + ε k Q u · g k, (19) Ω d = T d γ r = ∑ k λ k u × g k,

where the parameters

ξ k

ε k

, and

λ k

are given in equations (A.35)–(A.40) in terms of the diffusiophoretic coefficients

b k h

, the rate constants per unit area

κ ± c

, the slip length

b

, the molecular diffusivities

D k

, and the geometry of the Janus particle. The

3 × 3

identity matrix is

1

, while

(20) Q u ≡ u u − 131 .

The self-diffusiophoretic velocity, expressed in terms of the molecular concentrations

c ¯ k

extrapolated to the center of the particle, is

(21) V s d = ∑ k ζ k c ¯ k = ς κ + c c ¯ A − κ − c c ¯ B,

since the parameters

ζ k

may be written in the forms

ζ A = ς κ + c

and

ζ B = − ς κ − c

(see Appendix A, equations (A.33)-(A.34).

In the absence of a reaction, we recover the diffusiophoretic velocities given in Refs. [57, 58]:

(22) V d = ∑ k ξ k 0 g k with ξ k 0 = b k c + b k n 2 1 + 2 b / R, Ω d = ∑ k λ k 0 u × g k with λ k 0 = 9 16 R b k c − b k n .

Moreover, if the diffusiophoretic coefficients are the same on both hemispheres

b k c = b k n

, the angular velocity is equal to zero,

Ω d = 0

In the presence of a reaction, but without gradients (

g k = 0

), we have

κ + c c ¯ A ≠ κ − c c ¯ B

and the linear velocity reduces to the contribution of self-diffusiophoresis,

V d = V s d u

, characterizing the activity of the Janus particle.

The overall reaction rate can be written as follows:

(23) W r x n = k + c ¯ A − k − c ¯ B + ϖ k + g A − k − g B · u,

in terms of rate constants

k ± = Γ κ ± c

and a parameter

ϖ = O R

given in equation (A.45). In the absence of the concentration gradients, we recover the expression obtained in Ref. [40]. In the presence of the concentration gradients

g k

, there is an extra contribution depending on the direction

u

of the Janus particle. However, this last term is normally negligible because we typically have

R ∥ g k ∥ ≪ c ¯ k

for micrometric particles and macroscopic gradients of molecular concentrations.

We note that both the self-diffusiophoretic velocity (21) and the leading term of the reaction rate (23) are proportional to each other. Their ratio defines the self-diffusiophoretic parameter

χ

which was introduced in Refs. [39, 40],

(24) χ ≡ V s d k + c ¯ A − k − c ¯ B = ς Γ,

where the last equality was obtained using

k ± = Γ κ ± c

3. 　Active Suspension of Colloidal Motors

Less

3.1. 　Onsager’s Reciprocal Relations

We now show that Onsager’s principle of nonequilibrium thermodynamics [46–53] can be used to establish coupled diffusion-reaction equations of motion for active matter that are consistent with microreversibility. According to Onsager’s principle, currents are related to thermodynamic forces (or affinities) by

(25) J α = ∑ β L α β A β,

where the linear response coefficients satisfy the Onsager reciprocal relations,

(26) L α β = L β α,

if the affinities are even under time reversal. The thermodynamic entropy production rate density is given by

(27) σ s = k B ∑ α J α A α = k B ∑ α β L α β A α A β ≥ 0,

where

k B

is Boltzmann’s constant.

3.2. 　Mean-Field Equations for the Active Suspension

The system we consider is a dilute solution containing the reactive molecular

A

and

B

species together with colloidal motors

C

in an inert solvent

S

. The motors are spherical Janus particles and, as described in Section 2, have hemispherical catalytic surfaces where the reaction

A ⇄ B

takes place. Moreover, we suppose that the solution is globally at rest, so that the velocity field is equal to zero on scales larger than the size of colloids. The solution is described at the macroscale in terms of the molecular densities

n A r, t

and

n B r, t

, as well as the distribution function of the colloidal motors,

f r, u, t

, where

r = x, y, z

is the position and

u = sin θ cos φ, sin θ sin φ, cos θ

is the unit vector giving the orientation of the Janus particles (expressed in spherical coordinates in the laboratory frame). The distribution function is defined as

(28) f r, u, t ≡ ∑ i = 1 N C δ 3 r − r i t δ 2 u − u i t,

where

r i, u i i = 1 N C

are the positions and orientational unit vectors of the colloidal motors. For a dilute suspension, the evolution equation of this distribution function can be deduced from the Fokker-Planck equation for the probability that a single colloidal motor is located at the position

r

with the orientation

u

[39, 40, 41]. Once, this distribution function is known, we can obtain the successive moments of

u

(29) n C r, t ≡ ∫ f r, u, t d 2 u, (30) p r, t ≡ ∫ u f r, u, t d 2 u, (31) q r, t ≡ ∫ Q u f r, u, t d 2 u, ⋮

where

d 2 u = d cos θ d φ

n C

is the density or concentration of colloidal motors,

p

is the polar order parameter or polarization of the colloidal motors, and

q

is the traceless order parameter analogous to that for apolar nematic liquid crystals expressed in terms of the tensor (20) and, thus, satisfies

t r q = 0

At the macroscale, the reaction is

(32) A + C ⇄ k − k + B + C,

with the rate constants

k ±

. For the colloidal suspension treated here, the reaction should be described by a reaction rate density

w

that is proportional to the distribution function of colloidal motors and determined by the surface reaction taken into account with the boundary conditions (7) in Section 2.

The mean concentrations of molecular species are defined by

n k = 1 − ϕ c ¯ k

, where

ϕ = 4 π R 3 n C / 3

is the volume fraction of the suspension. Their corresponding gradients are related to those considered in Section 2 by

∇ n k = 1 − ϕ g k

for a dilute enough suspension. The coupled diffusion-reaction equations for the different species take the following forms:

(33) ∂ t n k + ∇ · j k = ν k w k = A, B, (34) ∂ t f + ∇ · V f − D t ∇ f = D r L r f,

where

j k

are the molecular current densities,

V

is the total drift velocity obtained by adding the drift velocity due to the external force

V ext = F ext / γ t

to the diffusiophoretic velocity (18) giving

(35) V = V s d u + ∑ k ξ k 1 + ε k Q u · ∇ n k + β D t F ext,

with the self-diffusiophoretic velocity (equation (21))

(36) V s d = ∑ k ζ k n k = ς κ + c n A − κ − c n B,

now expressed in terms of the mean concentrations

n k

, and the inverse temperature

β = k B T − 1

. In equation (34),

D t

is an effective translational diffusion coefficient related to the effective translational friction coefficient by Einstein’s formula

D t = k B T / γ t

and

D r

is an effective rotational diffusion coefficient related to the effective rotational friction coefficient by

D r = k B T / γ r

. Since the shear viscosity increases as

η ≃ η 0 1 + 2.5 ϕ

with the volume fraction

ϕ

of the suspension [26, 31], both friction coefficients

γ t

and

γ r

also increase, and the diffusion coefficients decrease. In particular, it is known that

D t ≃ D t 0 1 − 2.1 ϕ

[31]. A similar dependence on the volume fraction

ϕ

is expected for the parameters

ς

ξ k

ε k

, and

λ k

given in Appendix A, since these parameters are proportional to the diffusiophoretic coefficients

b k h

that are known to be inversely proportional to shear viscosity,

b k h ∝ η − 1

[54, 57, 58]. The effects of this dependence would manifest themselves if the colloidal suspension became dense enough. Here, such effects are assumed to play a negligible role.

The Janus particles have a spherical shape so that their random rotational and translational motions are decoupled. In this case, the rotational diffusion operator is given by

(37) L r f = 1 sin θ ∂ θ sin θ e − β U r ∂ θ e β U r f + 1 sin 2 θ ∂ φ e − β U r ∂ φ e β U r f,

expressed in terms of the rotational energy associated with the torque exerted by an external magnetic field

B

on some magnetic dipole

μ

of the particle [52] and that due to the diffusiophoretic effect, we have the following:

(38) U r = − μ B · u − γ r ∑ k λ k ∇ n k · u .

3.3. 　Translational and Rotational Current Densities

Before proceeding with nonequilibrium thermodynamics, we need to identify in equation (34) the current densities associated with the translational and rotational movements of the colloidal motors. The distribution function

f r, u

for colloidal Janus particles is defined in the five-dimensional space

x, y, z, θ, φ

, where

x, y, z

are the Cartesian coordinates for the position

r

and

θ, φ

the spherical coordinates for the orientation

u

. Vector calculus is used in these coordinates to obtain the corresponding gradients and divergences [59].

For the rotational degrees of freedom we have

(39) d u 2 = d θ 2 + sin 2 θ d φ 2 = g i j d q i d q j with g i j = 10 0 sin 2 θ .

The scalar product between a pair of rotational vectors

a r, b r ∈ ℝ 2

is given by

a r · b r = ∑ i, j = θ, φ g i j a r i b r j

, and the scalar product of such a vector with itself is denoted

a r 2 = a r · a r

. In spherical coordinates, the rotational gradient and divergence are given, respectively, by [59]

(40) grad r X = ∂ θ X 1 sin 2 θ ∂ φ X, (41) div r X r = 1 sin θ ∂ θ X r θ sin θ + ∂ φ X r φ .

In the five-dimensional space, the gradient is given by

(42) grad X = ∇ X grad r X with ∇ X = ∂ x X ∂ y X ∂ z X,

and the divergence of a five-dimensional vector,

X = X t, X r T

, is

(43) div X = ∇ · X t + div r X r .

Using these notations, equation (34) can be written in the form of a local conservation law involving the five-dimensional current density,

J C = j t, j r T

, as

(44) ∂ t f + div J C = 0, or ∂ t f + ∇ · j t + div r j r = 0,

with translational current density

(45) j t = V f − D t ∇ f = f V s d u + f ∑ k ξ k 1 + ε k Q u · ∇ n k − D t ∇ f + f β ∇ U t,

where

U t r = − F ext · r

is the translational potential energy due to the external force

F ext

, rotational current density

(46) j r = − D r e − β U r grad r e β U r f = f ∑ k λ k ∇ n k · grad r u − D r grad r f − f β μ B · grad r u,

and their translational and rotational divergences,

∇ · j t

and

div r j r = − D r L r f

, where

L r

is the operator (37).

3.4. 　Nonequilibrium Thermodynamics of the Active Suspension

Local thermodynamic equilibrium is assumed on scales larger than the size of the colloidal motors where the description by the mean-field equations (33) and (34) is valid and the fluid is at rest. According to this assumption, thermodynamic quantities can be locally expressed in the active suspension in terms of the molecular densities

n k r, t

and the distribution function

f r, u, t

. Furthermore, we suppose that the system is isothermal and isobaric and the solution is dilute in the species

A

B

, and

C

. The appropriate thermodynamic potential is thus Gibbs’ free energy given by the following volume integral of the corresponding density:

(47) G = ∫ d 3 r n S ψ S + ∑ k = A, B n k ψ k + n k k B T ln n k e n S + ∫ d 2 u f ψ C + f k B T ln f 4 π e n S + f U t r − f μ B · u,

where the first term is the contribution from the solvent

S

of density

n S

, the next terms from fuel

k = A

and product

k = B

dilutely dispersed in the solvent [52], and the last terms from all the orientations

u

of the colloidal motors moving in the mechanical potential energies due to the external force

F ext

and the external torque exerted by the magnetic field

B

on the magnetic dipoles of the colloids. We can thus deduce the following chemical potentials:

(48) μ S = δ G δ n S = ψ S − k B T n S n A + n B + n C, (49) μ k = δ G δ n k = ψ k + k B T ln n k n S k = A, B, (50) μ C = δ G δ f = ψ C + k B T ln f 4 π n S + U t r − μ B · u .

Here,

ψ k = μ k 0 + k B T ln n S / n 0

, where

μ k 0

is the standard chemical potential of species

k

and

n 0 = 1 mole / liter

is the standard concentration. Since the solution is dilute, we have taken the solvent density

n S

to be essentially uniform in space and constant in time.

Next, we use the principles of nonequilibrium thermodynamics in order to express the current densities in terms of the affinities or thermodynamic forces given in Table 1. For the reaction (32), the affinity is given by

(51) A r x n = 1 k B T μ A − μ B = ln k + n A k − n B,

and the corresponding current density is the rate density

w

introduced in equation (33). At chemical equilibrium, we have

A r x n = 0

w = 0

, and

k + n A, eq = k − n B, eq

. In the linear regime close to equilibrium where

δ n k = n k − n k, eq

, the chemical affinity (51) can be approximated as

(52) A r x n = δ n A n A, eq − δ n B n B, eq = 1 D r x n k + δ n A − k − δ n B,

where

(53) D r x n ≡ 12 k + n A + k − n B eq,

is the diffusivity of the reaction taking place on the colloidal motors [39, 40]. For the diffusion processes of species

k

, the affinity associated with the current density

j k

is given by

A k = − grad μ k / k B T

in terms of the chemical potential

μ k

. For molecular species, the gradient is tridimensional in Euclidian space, so that

A k = − ∇ μ k / k B T = − n k − 1 ∇ n k

. For the colloid with chemical potential (50), the affinity is given by the five-dimensional gradient (42) as

(54) A C = − grad μ C k B T = − 1 f ∇ f + f β ∇ U t grad r f − f β μ B · grad r u,

if the magnetic field

B

is uniform. In this five-dimensional space, the associated current density

J C = j t, j r T

given by equations (45) and (46) can thus be written in the following form:

(55) J C = f V s d u 0 + f ∑ k = A, B ξ k 1 + ε k Q u λ k grad r u · ∇ n k − f D t 1 0 0 D r 1 r · grad μ C k B T,

where

1 r

is the

2 × 2

identity matrix. In this form, we see that the first term is related to the reaction affinity since the self-diffusiophoretic velocity can be written as

V s d = χ D r x n A r x n

. The next two terms can be related to the affinities of molecular species, and the last term to the affinity of the colloidal species.

According to the Curie principle, there is no coupling between processes with different tensorial characters. However, the Janus particles have a director given by the unit vector

u

and we have adopted a description in terms of the distribution function

f r, u, t

for the Janus particles. Consequently, it is possible that a vectorial process such as diffusion may be coupled to a scalar process such as reaction if it is polarized by the unit vector

u

. If we introduce the densities

N C = f Δ 2 u

for Janus particles having their orientation

u

in cells with a size of

Δ 2 u

, along with the associated current densities,

(56) J C = J C Δ 2 u,

we may write a general coupling (25) of the following form:

(57) w j A j B J C = L r r L r A L r B L r C L A r L A A L A B L A C L B r L B A L B B L B C L C r L C A L C B L C C · A r x n − ∇ μ A k B T − ∇ μ B k B T − grad μ C k B T,

up to possible nonlinear contributions that may be required in order for the reaction rate to obey the mass-action law. In equation (57), we have that

L r r

1 × 1

L r k

1 × 3

L r C

1 × 5

L k r

3 × 1

L k l

3 × 3

L k C

3 × 5

L C r

5 × 1

L C k

5 × 3

, and

L C C

5 × 5

(for

k, l = A, B

According to Onsager’s reciprocal relations (26), the linear response coefficients should obey

(58) L r k = L k r T, L r C = L C r T, L k l = L l k T, L C C = L C C T, L k C = L C k T,

for

k = A, B

and where

T

again denotes the transpose.

We assume that the molecular species

A

and

B

undergo Fickian diffusion without cross-diffusion, so that

(59) L k l = D k n k δ k l 1,

and that the reaction rate does not depend on the gradients

∇ n A

∇ n B

, whereupon

(60) L r A = L r B = 0 .

This last assumption consists in neglecting the terms with the coefficient

ϖ

in equation (23), which is usually justified as mentioned in Section 2.

The scalar coefficient associated with the reaction can be identified as

(61) L r r = D r x n n C,

and the linear response coefficients

L C r

L C k

, and

L C C

in equation (57) can be determined using the current density (55), as described in Appendix B. As a consequence of Onsager’s reciprocal relations, we can conclude that the reaction rate and the current densities should be given by

(62) w = D r x n n C A r x n − χ D r x n ∫ u · ∇ f + f β ∇ U t d 2 u, (63) j k = − D k ∇ n k + n k ∫ ξ k 1 + ε k Q u · ∇ f + f β ∇ U t + λ k grad r u · grad r f − f β μ B · grad r u d 2 u .

In equation (62), the second term describes the reciprocal effects of diffusiophoresis back onto reaction. The second term in equation (63) is due to cross-diffusion between the molecular and colloidal species due to diffusiophoresis. We see that the linear response coefficients depend on the unit vector

u

in a manner similar to that already shown in Refs. [39, 40].

With respect to standard expressions, the terms involving the integral

∫ d 2 u

in equation (63) are required in order to satisfy Onsager’s reciprocal relations and for these quantities to be compatible with microreversibility. However, these extra terms can be shown to be negligible, although the reciprocal terms are not negligible in equations (45) and (46). In order to show that the extra terms are negligible, we suppose that the self-diffusiophoretic and diffusiophoretic velocities take the typical value

V s d ~ V d ~ 10 μ m / s

[60]. According to Ref. [58], the molecular gradients used in experiments of diffusiophoresis are of the order of

∇ n k ~ 105 mol / m 4

, so that diffusiophoretic parameters have the value

ξ k, ε k ~ 10 − 10 m 5 s − 1 mol − 1

. Moreover, we have

λ k ~ ξ k / R

, but since

grad r f ~ R ∇ f

, the effect of the coefficients

λ k

is again of the same order of magnitude as

ξ k

and

ε k

. Molecular diffusivities typically have the value

D k ~ 1 0 − 9 m 2 / s

, while the translational diffusion coefficient of a micrometric colloidal particle is of the order of

D t ~ 1 0 − 13 m 2 / s

. The molecular concentrations used in experiments on self-diffusiophoresis are about

n k ~ 103 mol / m 3

, while the density of micrometric colloidal particles is approximately

n C ~ 1018 m − 3 ~ 1 0 − 6 mol / m 3

, or lower. If we assume that the molecular and colloidal gradients take comparable values

∇ n k / n k ~ ∇ f / f

, the ratio between the extra term and the standard molecular diffusion term in equation (63) is given by

(64) n k ξ k ∇ f D k ∥ ∇ n k ∥ ~ ξ k f D k ~ 1 0 − 8,

which shows that the second term in equation (63) is negligible. Accordingly, the standard Fickian expressions

j k ≃ − D k ∇ n k

are very well justified for the molecular current densities. In the presence of colloidal motors, the expressions compatible with microreversibility are nevertheless given by equations (62) and (63). In contrast, the terms associated with the diffusiophoretic parameters in the colloidal current density (55) have effects that are not negligible.

The conclusion from these considerations is that active matter can be described as generalized diffusion-reaction processes in complete compatibility with microreversibility and Onsager’s reciprocal relations. In this way, the program of nonequilibrium thermodynamics is complete and application of equation (27) gives the following expression for the thermodynamic entropy production rate density:

(65) k B − 1 σ s = D r x n n C A r x n 2 + ∑ k = A, B D k ∇ n k 2 n k − 2 χ D r x n A r x n ∫ u · ∇ f + f β ∇ U t d 2 u − 2 ∫ ∑ k = A, B ∇ n k · ξ k 1 + ε k Q u · ∇ f + f β ∇ U t d 2 u − 2 ∫ ∑ k = A, B λ k ∇ n k · grad r u · grad r f − f β μ B · grad r u d 2 u + D t ∫ 1 f ∇ f + f β ∇ U t 2 d 2 u + D r ∫ 1 f grad r f − f β μ B · grad r u 2 d 2 u ≥ 0 .

The second law is satisfied if

D t ≫ χ 2 D r x n > 0

D k D t ≫ n C n k ξ k 2 > 0

D k D t ≫ n C n k ε k 2 > 0

, and

D k D r ≫ n C n k λ k 2 > 0

, which is expected.

The results derived in this section provide the basis for the analysis of collective effects in suspensions of active Janus particles. In Sections 4 and 5, we describe two collective phenomena that emerge from this theoretical framework: the effect of an external force and torque on the reaction rate, and a clustering instability.

4. 　Effect of External Force and Torque

Less

Using a thermodynamic formulation that is consistent with microreversibility, we showed earlier [39–41, 61] how the application of an external force and torque on a single colloidal motor can change the reaction rate on its surface and even lead to a net production of fuel rather than product. Now we show how these considerations can be extended to a suspension of such motors.

4.1. 　Local Evolution Equations

We suppose that the colloidal motors are subjected to an external force

F ext = F ext 1 z

and an external torque induced by an external magnetic field

B = B 1 z

exerted on the magnetic moment

μ

of the colloidal particles, both oriented in the

z

-direction. If

β μ B

is large enough, the distribution function is given by

(66) f r, u, t = n r, t β μ B 4 π sinh β μ B exp β μ B cos θ,

so that

p r, t = 1 z u z n C r, t

with

u z = coth β μ B − β μ B − 1

. Moreover, the terms with the coefficients

ξ k

ε k

, and

λ k

are assumed to be negligible in equation (34). If the concentrations are uniform in the

x

- and

y

-directions, the process is ruled by

(67) ∂ t n C + ∂ z χ u z k + n A − k − n B n C − D t ∂ z n C − β F ext n C = 0,

obtained by integrating equation (34) over the orientation

u

. This equation for

n C

is coupled to equation (33) with the Fickian molecular current densities

j k ≃ − D k ∇ n k

and the local reaction rate

(68) w = k + n A − k − n B n C − χ u z D r x n ∂ z n C − β F ext n C,

given by equation (62), as predicted by Onsager’s reciprocal relations.

4.2. 　Global Evolution Equations

Defining the mean value of the

z

-coordinate for the colloidal motors as

(69) z ≡ ∫ z n C d 3 r ∫ n C d 3 r,

and using equation (67), we obtain the evolution equation,

(70) d z d t = χ u z w r x n + β D t F ext,

with mean reaction rate,

(71) w r x n ≡ k + ∫ n A n C d 3 r ∫ n C d 3 r − k − ∫ n B n C d 3 r ∫ n C d 3 r .

Furthermore, integrating equation (33) with

k = B

over the position

r

with the local rate (68), we get the total reaction rate

(72) d N r x n d t = d N B d t = − d N A d t = N C w r x n + χ u z D r x n β F ext,

where

N C = ∫ n C d 3 r

is the total number of colloidal motors in the suspension. Equations (70) and (72) have precisely the same structure as for a single colloidal motor. However, one should note that the mean reaction rate in equation (71) contains spatial correlations between the solute and colloid concentration fields. Given the structure of the equations, the results obtained in Refs. [39, 40] also apply here. In particular, there exists a regime where the entire ensemble of colloidal motors is propelled and carries out work against the external force by consuming fuel. In addition, there is also a regime where fuel is synthesized if the external force that opposes motion is sufficiently large to reverse the reaction

A ⟶ B

. The efficiencies of these processes are given by the same expressions as in Refs. [39, 40].

5. 　Clustering Instability and Pattern Formation

Less

The equations of motion developed in Section 3 that describe a dilute suspension of colloidal motors moving in a dilute solution of fuel

A

and product

B

molecular species will be shown in this section to lead to a clustering instability. This instability can be described by the mean-field equations obtained above for the concentrations of the molecular species and the distribution function of the colloidal motors in the absence of an external force and torque (

F ext = 0

and

B = 0

). A number of other mean-field descriptions that predict instabilities and the formation of various clustering states of collections of diffusiophoretic colloidal particles have appeared in literature [62, 18, 6, 19, 20], and use techniques involving coupled moment equations similar to those adopted in this section.

5.1. 　Molecular Diffusion and Reaction

The equation for the colloidal motors is coupled to the reaction-diffusion equations for the molecular species

A

and

B

, accounting for the fact that the reaction

A ⇄ B

occurs both at the surface of the catalytic hemisphere of the colloids and in the bulk:

(73) ∂ t n k = D k ∇ 2 n k + ν k w tot,

where the total reaction rate density is given by

(74) w tot = k + n A − k − n B n C + k + 2 n A − k − 2 n B .

The system is driven out of equilibrium if

k + / k − ≠ k + 2 / k − 2

[61].

5.2. 　Colloidal Density and Polarization

If the second moment (31) as well as higher moments are assumed to be negligible, the evolution equations for the density of colloidal motors (29) and the polarization (30) are given by

(75) ∂ t n C + ∇ · n C ∑ k ξ k ∇ n k + V s d p = D t ∇ 2 n C, (76) ∂ t p + ∇ · p ∑ k ξ k ∇ n k + 13 ∇ V s d n C + 15 Δ : ∇ ∑ k ε k ∇ n k p = D t ∇ 2 p − 2 D r p + 23 n C ∑ k λ k ∇ n k,

in terms of the fourth-order tensor

Δ

with the following components:

Δ i j m n = δ i j δ m n + δ i m δ j n − 2 / 3 δ i n δ j m

D r

is large enough so that

2 D r p

dominates the other terms involving

p

in equation (76), we can neglect these other terms and this equation can be inverted to obtain

(77) p ≃ 1 6 D r − V s d ∇ n C + n C ∑ k 2 λ k − ζ k ∇ n k,

under which circumstances the field

p

is driven by the gradients of the colloid and species densities. Substituting this result into equation (75) for the density

n C

of colloidal particles, we find

(78) ∂ t n C + ∇ · n C ∑ k ξ k + V s d 6 D r 2 λ k − ζ k ∇ n k − D t eff ∇ n C = 0,

with the effective diffusion coefficient

(79) D t eff ≡ D t + V s d 2 6 D r,

expressing the enhancement of diffusivity due to the self-diffusiophoretic activity [17].

5.3. 　Coupled Colloidal and Molecular Diffusion-Reaction Equations

In the following, we suppose that the diffusion coefficient is the same for both molecular species:

D ≡ D A = D B

. Consequently,

n 0 = n A + n B

remains uniform during the time evolution if initially so. Therefore,

n B = n 0 − n A

is known and only

n A

needs to be determined. Introducing the notations

(80) a ≡ n A, c ≡ n C,

we have the following coupled equations describing the system:

(81) ∂ t a = D ∇ 2 a − W tot, (82) W tot = c K a − K 0 + K 2 a − K 20, (83) ∂ t c = ∇ · D t + τ r V s d 2 ∇ c − ξ + σ V s d c ∇ a, (84) V s d = ζ a − V 0,

with

(85) K ≡ k + + k −, K 0 ≡ k − n 0, K 2 ≡ k + 2 + k − 2, K 20 ≡ k − 2 n 0, τ r ≡ 6 D r − 1, ξ ≡ ξ A − ξ B, λ ≡ λ A − λ B, ζ ≡ ζ A − ζ B = ς κ + c + κ − c, V 0 ≡ − ζ B n 0 = ς κ − c n 0, σ ≡ τ r 2 λ − ζ .

Moreover, consistency with the existence of equilibrium requires that

ζ / K = V 0 / K 0 = ς / Γ = χ

is equal to the diffusiophoretic parameter (24) that is the ratio between the self-diffusiophoretic velocity (21) and the leading term of the overall reaction rate (23).

For this system, there exists a uniform nonequilibrium steady state, where

c

keeps its initial uniform value

c 0

and the molecular density is also uniform at the value

(86) a 0 = c 0 K 0 + K 20 c 0 K + K 2,

in order to satisfy the stationary condition

W tot = 0

. For this molecular concentration

a = a 0

, we notice that the rate

K a − K 0

of the catalytic reaction on the colloids is not vanishing under the nonequilibrium condition

k + / k − ≠ k + 2 / k − 2

5.4. 　Pattern Formation

To analyze the stability of this homogeneous steady state, for simplicity we consider a one-dimensional system where the fields

a

and

c

only depend on the variable

z

. Accordingly, the gradients

∇

can be replaced by partial derivatives

∂ z

in equations (81) and (83). The set of equations (81)–(84) is then numerically integrated by spatial discretization over the grid

z = i Δ z

with

i = 1, 2, ⋯, M

with

Δ z = 0.1

and

M = 500

. The integration is performed with a Runge-Kutta algorithm of varying order 4-5 over a long enough time interval to reach a steady state. Figure 2 shows numerical results for the parameter values:

(87) c 0 = 1, n 0 = 10, D = 1, D t = 1, τ r = 1, ξ = − 3, σ = − 2, V 0 = 0.5, ζ = 0.1, K = 0.2, K 0 = 1, K 2 = 0.3,

and increasing values of

K 20

. We observe the formation of clusters of colloidal motors in regions where the fuel

A

is depleted, as expected.

5.5. 　Linear Analysis of the Clustering Instability

The threshold of this clustering instability can be found from a linear stability analysis. Linearizing the equations around the uniform steady state, we find that the perturbations obey

(88) ∂ t δ a δ c = D ∂ z 2 − K ~ − w − ρ ∂ z 2 D t eff ∂ z 2 δ a δ c,

with

(89) K ~ ≡ c 0 K + K 2, w ≡ K a 0 − K 0, ρ ≡ c 0 ξ + σ V s d, D t eff ≡ D t + τ r V s d 2, V s d ≡ ζ a 0 − V 0 .

Supposing that the perturbations behave as

δ a, δ c ~ exp i q z + s t

, we obtain the dispersion relations

(90) s ± q = − 12 K ~ + D + D t eff q 2 ± 12 K ~ + D − D t eff q 2 2 − 4 ρ w q 2 .

These dispersion relations are depicted in Figures 3(a)–3(c), respectively, below, at, and beyond the threshold. The leading dispersion relation is associated with the conserved unstable mode of the colloidal motors because

s + 0 = 0

. The subleading dispersion relation is associated with the reactive mode of the molecular species because

s − 0 = − K ~

. We notice that, since

K ~ > 0

, there is no possibility for a Hopf bifurcation to uniform oscillatory behavior, which would be the case if the eigenvalues

s ± 0

were complex. We also note that there is no wavelength selection at the level of linear stability analysis in this clustering instability.

Therefore, instability manifests itself if

(91) ρ w + K ~ D t eff < 0,

and the threshold is given by the condition

(92) ρ w + K ~ D t eff = 0,

which leads to the value

K 20 ≃ 1.89817

for the parameter values (87).

The dispersion relations can also be obtained from the evolution equation (34) for the distribution function. Supposing that

f = f z, θ

and

a = a z

, we have

(93) ∂ t f + ∂ z V d z f − D t ∂ z f = D r sin θ ∂ θ sin θ ∂ θ f + 2 λ cos θ f ∂ z a,

where

(94) V d z = V s d cos θ + ξ + ε cos 2 θ − 13 ∂ z a,

with

ε ≡ ε A − ε B

. Equation (93) is coupled to the diffusion-reaction equation (81) for the concentration field

a

with

c = ∫ f d 2 u

The linear stability analysis can be carried out for equation (93) coupled to equation (81) with the rate (82) in a similar manner to that for equations (81)–(84). This analysis is presented in Appendix C. The dispersion relations can be computed by truncating the infinite matrix (C.5) in order to obtain the eigenvalues as a function of the wave number

q

. The result converges to the dispersion relations shown in Figures 3(d)–3(f) below, at, and beyond the threshold, for the parameter values (87) and

ε = 1

. The convergence occurs faster for the leading dispersion relation than for the subleading ones. For the chosen parameter values, we can see that the approximation where we suppose that the vector field

p

is driven by the gradients (which corresponds to truncating to a

2 × 2

matrix) constitutes a good approximation to describe the instability. Indeed, the leading dispersion relation of Figure 3(b) is already very close to that in Figure 3(e).

The conclusion is that equations (81)–(84) provide a robust description of the clustering instability and of the emerging patterns.

6. 　Conclusion

Less

Autonomous motion is not possible at equilibrium and active matter relies on the presence of nonequilibrium constraints to drive the system out of equilibrium. As a result the theoretical formulations provided by nonequilibrium thermodynamics and statistical mechanics are a natural starting point for the description of such systems.

Many of the active matter systems currently under study involve active agents such as molecular machines or self-propelled colloidal particles with linear dimensions ranging from tens of nanometers to micrometers. The transition from microscopic to macroscopic description for fluids containing active agents of such sizes takes place in the upper range of this scale. Suspensions of active colloidal particles are interesting in this connection since, as described earlier in this paper, the colloidal particles are large compared to the molecules of the medium in which they reside. The dynamics of the suspension can then be described by considering the equations for the positions, velocities, and orientations of the colloidal particles in the medium, or through field equations that describe the densities of these particles.

Nonequilibrium thermodynamics provides a set of principles that these systems must obey. Most important among these is microreversibility that stems from the basic time reversal character of the microscopic dynamics. On the macroscale, this principle manifests itself in Onsager’s reciprocal relations that govern what dynamical processes are coupled and how they are described. For example, for single Janus particles propelled by a self-diffusiophoretic mechanism, microreversibility implies the existence of reciprocal effect where the reaction rate depends on an applied external force [39–41, 61].

This paper extended the nonequilibrium thermodynamics formulation to the collective dynamics of ensembles of diffusiophoretic Janus colloids. In particular, we considered Janus colloids driven by both self-diffusiophoresis arising from reactions on the motor catalytic surface as well as motion arising from an external concentration gradient. This latter contribution is essential for the extension of the theory to collective motor dynamics. The resulting formulation is consistent with microreversibility and an expression for the entropy production is provided. From this general formulation of collective dynamics, one can show that if an external force and torque are applied to the system, the overall reaction rate depends on the applied force. In addition, a stability analysis of the equations governing the collective behavior predicts the existence of a clustering instability seen in many experiments of Janus colloids. Such considerations can be extended to ensembles of thermophoretic Janus colloids [63].

Funding

Less

Fonds de la Recherche Scientifique - FNRS(T.0094.16)
Université Libre de Bruxelles
Compute Canada
Natural Sciences and Engineering Research Council of Canada

References

Less

Vicsek

, and Zafeiris

, “Collective motion,” Physics Reports, vol. 517, no. 3-4, pp. 71–140, 2012

Romanczuk

, Bär

, Ebeling

, Lindner

, and Schimansky-Geier

, “Active Brownian particles,” The European Physical Journal Special Topics, vol. 202, no. 1, pp. 1–162, 2012

Aronson

I. S.

, “Collective behavior in out-of-equilibrium colloidal suspensions,” Comptes Rendus Physique, vol. 14, no. 6, pp. 518–527, 2013

Cates

M. E.

, and Tailleur

, “When are active Brownian particles and run-and-tumble particles equivalent? Consequences for motility-induced phase separation,” EPL, vol. 101, no. 2, article 20010, 2013

Elgeti

, Winkler

R. G.

, and Gompper

, “Physics of microswimmers—single particle motion and collective behavior: a review,” Reports on Progress in Physics, vol. 78, no. 5, article 056601, 2015

Speck

, Menzel

A. M.

, Bialké

, and Löwen

, “Dynamical mean-field theory and weakly non-linear analysis for the phase separation of active Brownian particles,” The Journal of Chemical Physics, vol. 142, no. 22, article 224109, 2015

Bechinger

, Di Leonardo

, Löwen

, Reichhardt

, Volpe

, and Volpe

, “Active particles in complex and crowded environments,” Reviews of Modern Physics, vol. 88, no. 4, article 045006, 2016

Marchetti

M. C.

, Fily

, Henkes

, Patch

, and Yllanes

, “Minimal model of active colloids highlights the role of mechanical interactions in controlling the emergent behavior of active matter,” Current Opinion in Colloid & Interface Science, vol. 21, pp. 34–43, 2016

Zöttl

, and Stark

, “Emergent behavior in active colloids,” Journal of Physics: Condensed Matter, vol. 28, no. 25, article 253001, 2016

Ramaswamy

, “Active matter,” Journal of Statistical Mechanics: Theory and Experiment, vol. 2017, no. 5, article 054002, 2017

Gompper

, Winkler

R. G.

, Speck

, Solon

, Nardini

, Peruani

, Löwen

, Golestanian

, Kaupp

U. B.

, Alvarez

, Kiørboe

, Lauga

, Poon

W. C. K.

, DeSimone

, Muiños-Landin

, Fischer

, Söker

N. A.

, Cichos

, Kapral

, Gaspard

, Ripoll

, Sagues

, Doostmohammadi

, Yeomans

J. M.

, Aranson

I. S.

, Bechinger

, Stark

, Hemelrijk

C. K.

, Nedelec

F. J.

, Sarkar

, Aryaksama

, Lacroix

, Duclos

, Yashunsky

, Silberzan

, Arroyo

, and Kale

, “The 2020 motile active matter roadmap,” Journal of Physics: Condensed Matter, vol. 32, no. 19, article 193001, 2020

Theurkauff

, Cottin-Bizonne

, Palacci

, Ybert

, and Bocquet

, “Dynamic clustering in active colloidal suspensions with chemical signaling,” Physical Review Letters, vol. 108, no. 26, article 268303, 2012

Gao

, Pei

, Dong

, and Wang

, “Catalytic iridium-based Janus micromotors powered by ultralow levels of chemical fuels,” Journal of the American Chemical Society, vol. 136, no. 6, pp. 2276–2279, 2014

Wang

, Duan

, Ahmed

, Sen

, and Mallouk

T. E.

, “From one to many: dynamic assembly and collective behavior of self-propelled colloidal motors,” Accounts of Chemical Research, vol. 48, no. 7, pp. 1938–1946, 2015

Ginot

, Theurkauff

, Detcheverry

, Ybert

, and Cottin-Bizonne

, “Aggregation-fragmentation and individual dynamics of active clusters,” Nature Communications, vol. 9, no. 1, article 696, 2018

Thakur

, and Kapral

, “Collective dynamics of self-propelled sphere-dimer motors,” Physical Review E, vol. 85, no. 2, article 026121, 2012

Kapral

, “Perspective: nanomotors without moving parts that propel themselves in solution,” The Journal of Chemical Physics, vol. 138, no. 2, article 020901, 2013

Saha

, Golestanian

, and Ramaswamy

, “Clusters, asters, and collective oscillations in chemotactic colloids,” Physical Review E, vol. 89, no. 6, article 062316, 2014

Pohl

, and Stark

, “Self-phoretic active particles interacting by diffusiophoresis: a numerical study of the collapsed state and dynamic clustering,” The European Physical Journal E, vol. 38, no. 8, article 93, 2015

Liebchen

, Marenduzzo

, and Cates

M. E.

, “Phoretic interactions generically induce dynamic clusters and wave patterns in active colloids,” Physical Review Letters, vol. 118, no. 26, article 268001, 2017

Huang

M.-J.

, Schofield

, and Kapral

, “Chemotactic and hydrodynamic effects on collective dynamics of self-diffusiophoretic Janus motors,” New Journal of Physics, vol. 19, no. 12, article 125003, 2017

Colberg

P. H.

, and Kapral

, “Many-body dynamics of chemically-propelled nanomotors,” The Journal of Chemical Physics, vol. 147, no. 6, article 064910, 2017

Robertson

, Huang

M.-J.

, Chen

J.-X.

, and Kapral

, “Synthetic nanomotors: working together through chemistry,” Accounts of Chemical Research, vol. 51, no. 10, pp. 2355–2364, 2018

Huang

M.-J.

, Schofield

, Gaspard

, and Kapral

, “From single particle motion to collective dynamics in Janus motor systems,” The Journal of Chemical Physics, vol. 150, no. 12, article 124110, 2019

Stark

, “Artificial chemotaxis of self-phoretic active colloids: collective behavior,” Accounts of Chemical Research, vol. 51, no. 11, pp. 2681–2688, 2018

Einstein

Investigations on the Theory of the Brownian Movement, Dover, New York, 1956

Smoluchowski

M. V.

, “Über Brownsche Molekularbewegung unter Einwirkung äußerer Kräfte und deren Zusammenhang mit der verallgemeinerten Diffusionsgleichung,” Annalen der Physik, vol. 353, no. 24, pp. 1103–1112, 1916

von Smoluchowski

, “Drei Vorträge über Diffusion, Brownsche Molekularbewegung und Koagulation von Kolloidteilchen,” Physikalische Zeitschrift, vol. 17, pp. 557–571, 1916

Smoluchowski

M. v.

, “Versuch einer mathematischen Theorie der Koagulationskinetik Kolloider Lösungen,” Zeitschrift für Physikalische Chemie, vol. 92U, no. 1, pp. 129–168, 1918

Kapral

, “Kinetic theory of chemical reactions in liquids,” Advances in Chemical Physics, Wiley, vol. 48, pp. 71–179, 1981

Happel

, and Brenner

Low Reynolds Number Hydrodynamics, Martinus Nijhoff Publishers, The Hague, 1983

Dhont

J. K. G.

An Introduction to Dynamics of Colloids, Elsevier, Amsterdam, 1996

Nägele

, “Colloidal hydrodynamics,” Physics of Complex Colloids, Proc. of the Int. School of Physics Enrico Fermi, Bechinger

, Sciortino

, and Ziherl

, Eds., IOS Press, Amsterdam, pp. 507–601, 2013

Gerritsma

, and Gaspard

, “Chemomechanical coupling and stochastic thermodynamics of the F1-ATPase molecular motor with an applied external torque,” Biophysical Reviews and Letters, vol. 5, no. 4, pp. 163–208, 2011

Seifert

, “Stochastic thermodynamics, fluctuation theorems and molecular machines,” Reports on Progress in Physics, vol. 75, no. 12, article 126001, 2012

Speck

, “Stochastic thermodynamics for active matter,” EPL, vol. 114, no. 3, article 30006, 2016

Jülicher

, Grill

S. W.

, and Salbreux

, “Hydrodynamic theory of active matter,” Reports on Progress in Physics, vol. 81, no. 7, article 076601, 2018

Fodor

É.

, and Cristina Marchetti

, “The statistical physics of active matter: from self-catalytic colloids to living cells,” Physica A: Statistical Mechanics and its Applications, vol. 504, pp. 106–120, 2018

Gaspard

, and Kapral

, “Communication: mechanochemical fluctuation theorem and thermodynamics of self-phoretic motors,” The Journal of Chemical Physics, vol. 147, no. 21, article 211101, 2017

Gaspard

, and Kapral

, “Fluctuating chemohydrodynamics and the stochastic motion of self-diffusiophoretic particles,” The Journal of Chemical Physics, vol. 148, no. 13, article 134104, 2018

Gaspard

, and Kapral

, “Thermodynamics and statistical mechanics of chemically powered synthetic nanomotors,” Advances in Physics: X, vol. 4, no. 1, article 1602480, 2019

Speck

, “Thermodynamic approach to the self-diffusiophoresis of colloidal Janus particles,” Physical Review E, vol. 99, no. 6, article 060602, 2019

Golestanian

, Liverpool

T. B.

, and Ajdari

, “Designing phoretic micro- and nano-swimmers,” New Journal of Physics, vol. 9, no. 5, article 126, 2007

Oshanin

, Popescu

M. N.

, and Dietrich

, “Active colloids in the context of chemical kinetics,” Journal of Physics A: Mathematical and Theoretical, vol. 50, no. 13, article 134001, 2017

Reigh

S. Y.

, Huang

M.-J.

, Schofield

, and Kapral

, “Microscopic and continuum descriptions of Janus motor fluid flow fields,” Philosophical Transactions of the Royal Society A, vol. 374, no. 2080, article 20160140, 2016

Onsager

, “Reciprocal relations in irreversible processes I,” Physical Review, vol. 37, no. 4, pp. 405–426, 1931

Onsager

, “Reciprocal relations in irreversible processes II,” Physical Review, vol. 38, no. 12, pp. 2265–2279, 1931

Prigogine

Introduction to Thermodynamics of Irreversible Processes, Wiley, New York, 1967

Haase

Thermodynamics of Irreversible Processes, Dover, New York, 1969

de Groot

S. R.

, and Mazur

Nonequilibrium Thermodynamics, Dover, New York, 1984

Nicolis

, “Irreversible thermodynamics,” Reports on Progress in Physics, vol. 42, no. 2, pp. 225–268, 1979

Landau

L. D.

, and Lifshitz

E. M.

Statistical Physics, Part 1, Pergamon Press, 3rd, Oxford, UK, 1980

Callen

H. B.

Thermodynamics and an Introduction to Thermostatistics, Wiley, 2nd, New York, 1985

Gaspard

, and Kapral

, “Nonequilibrium thermodynamics and boundary conditions for reaction and transport in heterogeneous media,” The Journal of Chemical Physics, vol. 148, no. 19, article 194114, 2018

Mazur

, and Bedeaux

, “A generalization of Faxén's theorem to nonsteady motion of a sphere through an incompressible fluid in arbitrary flow,” Physica, vol. 76, no. 2, pp. 235–246, 1974

Albano

A. M.

, Bedeaux

, and Mazur

, “On the motion of a sphere with arbitrary slip in a viscous incompressible fluid,” Physica A: Statistical Mechanics and its Applications, vol. 80, no. 1, pp. 89–97, 1975

Anderson

J. L.

, “Colloid transport by interfacial forces,” Annual Review of Fluid Mechanics, vol. 21, no. 1, pp. 61–99, 1989

Anderson

J. L.

, and Prieve

D. C.

, “Diffusiophoresis caused by gradients of strongly adsorbing solutes,” Langmuir, vol. 7, no. 2, pp. 403–406, 1991

Brenner

, “Coupling between the translational and rotational brownian motions of rigid particles of arbitrary shape I. Helicoidally isotropic particles,” Journal of Colloid Science, vol. 20, no. 2, pp. 104–122, 1965

Wang

Nanomachines: Fundamentals and Applications, Wiley-VCH, Weinheim, Germany, 2013

Huang

M.-J.

, Schofield

, Gaspard

, and Kapral

, “Dynamics of Janus motors with microscopically reversible kinetics,” The Journal of Chemical Physics, vol. 149, no. 2, article 024904, 2018

Speck

, Bialké

, Menzel

A. M.

, and Löwen

, “Effective Cahn-Hilliard equation for the phase separation of active Brownian particles,” Physical Review Letters, vol. 112, no. 21, article 218304, 2014

Gaspard

, and Kapral

, “The stochastic motion of self-thermophoretic Janus particles,” Journal of Statistical Mechanics: Theory and Experiment, vol. 2019, no. 7, article 074001, 2019

Appendix

Less

Year 2020 volume 2020 Issue 5

PDF

Cite this Article

BibTeX

Article Info

doi: 10.34133/2020/9739231

Receive Date：2020-03-15
Online Date：2026-10-27
Published：2020-05-20

Article Data

Affiliations

History

Received：2020-03-15
Accepted：2020-04-19

Funding

Fonds de la Recherche Scientifique - FNRS(T.0094.16)

Université Libre de Bruxelles

Compute Canada

Natural Sciences and Engineering Research Council of Canada

Affiliations

¹Center for Nonlinear Phenomena and Complex Systems, Université Libre de Bruxelles (U.L.B.), Code Postal 231, Campus Plaine, B-1050 Brussels, Belgium

²Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada, M5S 3H6

References

Share

https://castjournals.cast.org.cn/joweb/research/EN/10.34133/2020/9739231

Share to

Scan QR to access full text

Cite this article

BibTeX

Citations

表12种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

Table 1: Current densities and corresponding affinities or thermodynamic forces in the active suspension: $w$ is the reaction rate density introduced in equation (33) corresponding to the affinity (51), $j A$ and $j B$ are the molecular current densities of fuel $A$ and product $B$ , $j t$ is the translational current density of colloids given by equation (45), $j r$ is the rotational current density of colloids given by equation (46). The translational and rotational current densities of colloids form the five-dimensional current density (55) according to $J C = j t, j r T$ . Similarly, the translational and rotational affinities of colloids form the five-dimensional affinity (54). $β = k B T − 1$ denotes the inverse temperature, $μ k$ the chemical potentials, $∇ = ∂ x, ∂ y, ∂ z$ the gradient in Euclidean space, and $grad r$ the rotational gradient (40). We note that colloidal motors with the given orientation $u$ are considered as so many independent species in the free energy (47), which is expressed by equation (56).

Process	Current	Affinity	Dimension
Reaction	$w$	$β μ A − μ B$	$1$
Molecular diffusion of fuel	$j A$	$− ∇ β μ A$	$3$
Molecular diffusion of product	$j B$	$− ∇ β μ B$	$3$
Translational diffusion of colloids	$j t$	$− ∇ β μ C$	$3$
Rotational diffusion of colloids	$j r$	$− grad r β μ C$	$2$

Figure 1: Schematic representation of a Janus particle with its catalytic (C) and noncatalytic (N) hemispheres where the surface reaction (1) takes place between fuel $A$ and product $B$ supplied by the solution surrounding the particle. The particle is also subjected to some external force $F ext$ and torque $T ext$ . The position of its center of mass is $R$ , and $u$ is the unit vector giving its orientation and pointing in the direction of the catalytic hemisphere.

Figure 2: Nonequilibrium steady state of the one-dimensional system for the parameter values (87): (a) with $K 20 = 2$ ; (b) with $K 20 = 2.5$ .

Figure 3: Dispersion relations of linear stability analysis for $K 20 = 1$ , $K 20 = 1.89817$ , and $K 20 = 3$ , respectively, below, at, and beyond the threshold of clustering instability. The dispersion relations are obtained in (a)–(c) with the approximation (90), and in (d)–(f) by truncating equation (C.5) into a $5 × 5$ matrix. The other parameter values are given in equation (87) and $ε = 1$ .

Articles: Latest Articles; Most Read; Collections

Updates: Events; News; Multimedia

About: About Us

Contact

No. 86 Xueyuan South Road, Haidian District, Beijing

100081

010-62662699(Operation and Technical Support)

qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House