![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
Colloidal motors with micrometer dimensions and no moving parts can be propelled by self-diffusiophoresis. Coupling between molecular concentration gradients generated by asymmetric surface chemical reactions and the velocity slip between colloidal particle and the surrounding fluid solution is responsible for propulsion. The interfacial properties involved in this propulsion mechanism can be described by nonequilibrium thermodynamics and statistical mechanics, disclosing the fundamental role of microreversibility in the coupling between motion and reaction. Among other phenomena, the approach predicts that propulsion by fuel consumption has the reciprocal effect of fuel synthesis by mechanical action.
GRAPHICAL ABSTRACT
![](/cms/asset/cbde911d-8560-4aa1-9c63-ad0e99ab5917/tapx_a_1602480_uf0001_oc.jpg)
1. Introduction
Artificial or natural engines, machines, and motors with dimensions ranging from molecular to macroscopic scales are able to convert free energy into motion by various mechanisms that depend on their size. At the smallest scale, we find synthetic molecular machines made of catenanes and rotaxanes. These nanometric machines are driven by successive changes in the surrounding concentrations and/or by photochemical means [Citation1–Citation4]. On the same scale, organometallic complexes carrying out homogeneous catalysis also undergo directional motion, in particular during polymerization, or enhanced diffusion. These machines are so small that the substrate molecules remain in the periphery of their catalytic centers.
Biomolecular motors are larger molecular assemblies with enzymatic activity. Examples are rotary motors such as the ATP synthase and bacterial flagellar rotary motors, and the actin-myosin, kinesin-microtubule, and dynein-microtubule linear motors [Citation5,Citation6]. Since they have sizes of the order of 10 nm or more, like other enzymes, they can accommodate internal catalytic sites or channels. They are powered by ATP hydrolysis or transmembrane ionic currents.
Synthetic motors with sizes ranging from tens of nanometers to micrometers have been made and studied [Citation7]. They are solid particles composed of several materials, and some portions of the surfaces of these motors are catalytically active. The processes that drive propulsion take place at the interface between the solid particle and the surrounding fluid, and the properties of the interface thus play a key role in the mechanism. Since the thickness of the interface is much smaller than the particle diameter, a continuous-medium description may be used. Indeed, it is known that continuous-medium equations such as the Navier-Stokes equations are valid for fluid flows down to the nanoscale [Citation8]. At this level of description, the velocity or concentration fields should satisfy boundary conditions incorporating the relevant interfacial properties, which are themselves determined by the molecular structure and molecular dynamics of the interface.
The colloidal motors that are the focus of this review are paradigmatic examples of synthetic motors because their propulsion is autonomous once they are immersed in a suitable solution; in addition, the collective behavior of many motors has interesting properties. In this regard, a comparison with biological systems is instructive. In muscles, the actin-myosin linear motors compose sarcomeres, which are the intracellular structures performing work during muscle contraction. Although the contraction velocity of a single myosin head running along an actin filament is of the order of a micrometer per second, the muscle can reach a velocity multiplied by the number of sarcomeres composing the muscular fibers. This example shows what collective motion of micrometric machines can achieve by cooperative action.
The paper is organized as follows. The phenomenology of synthetic chemically powered motors is described in Section 2. Fluctuating chemohydrodynamics is developed in Section 3 and solved for spherical Janus particles in Section 4, in order to obtain the coupled overdamped Langevin equations ruling their dynamics in Section 5. Fluctuation theorems and microreversibility are discussed in Section 6. A comparison with biomolecular motors is given in Section 7. Section 8 is devoted to the implications for collective behavior. Concluding remarks are given in Section 9.
2. Synthetic chemically powered motors
Some of the first micron-scale synthetic motors were bimetallic rods using hydrogen peroxide fuel [Citation9, Citation10]. Subsequently, motors with other shapes made from other materials were constructed and studied [Citation7,Citation11–Citation14]. An important class of colloidal particles is made of silica spheres with a platinum cap. Heterogeneous catalysis will take place on this metallic surface once immersed in a solution containing some reactant, such as hydrogen peroxide. These colloidal motors can be described by simplified models assuming the particle is spherical with a catalytic portion or dimers composed of linked catalytic and non-catalytic spheres [Citation15].
It is worth noting that liquid droplets can also exhibit self-propulsion [Citation16–Citation19]. Here, the propulsion mechanism is the Marangoni effect due to the dependence of the surface tension between the droplet and the surrounding solution on tensio-active species. The surface tension is an important interfacial property that determines the boundary conditions on the velocity fields inside and outside the droplet in the framework of hydrodynamics. Although the surface tension as well as other interfacial properties are determined by the molecular structure of the interface, the interface is assumed to be arbitrarily thin in the continuous-medium description provided by hydrodynamics for the colloidal particles immersed in a fluid.
An important difference between machines with dimensions that are above and below micrometers is that the thermal and molecular fluctuations become negligible at macroscopic scales although their effects are essential on micrometric or submicrometric scales, as it is well known for Brownian motion. Therefore, small machines and motors should be described in terms of stochastic processes and statistical mechanics. Indeed, these theories provide a way to describe random motions in terms of probabilities, mean values, variances, and other statistical moments or cumulants. The effects of thermal fluctuations on continuous velocity or concentration fields can also be described, for instance, using fluctuating hydrodynamics [Citation20,Citation21]. Instead, on macroscopic scales, the amplitudes of the thermal fluctuations become smaller than the size of the machine, so that their effects may be neglected.
A distinction should be made between motors without moving parts that are propelled by the reactions taking place at their interfaces, and those that are propelled by cyclic changes in their shapes [Citation22]. For motors such as bacteria propelled by flagellae, amoebae moving by the dynamics of their cytoskeleton, unicellular organisms swimming with ciliae, as well as animals and macroscopic engines, it is the time dependence of the shape that performs work and leads to propulsion. For motors without moving parts, propulsion has its origin in slippage of the velocity at the interface between the solid particle and the surrounding fluid due to phoretic effects that include thermophoresis involving temperature gradients along the interface, diffusiophoresis involving concentration gradients, and electrodiffusiophoresis in electrolyte solutions. More specifically, these phoretic effects are induced by the molecular interaction forces between the solid particle and the solute species in the solution. For these colloidal motors, surface reactions generate inhomogeneous concentration fields of the reactant and product species in the solution around the solid particle giving rise to concentration gradients along the surface of the particle. Through diffusiophoresis, these concentration gradients induce a slippage between the fluid and solid velocities, giving rise to the propulsion of the particle in the fluid.
A key issue is to understand the mechanism by which the transduction of the available free energy into motion takes place. Every motor uses some source of free energy, which can be heat, chemical reaction, or electron transfer provided by the voltage of a battery. Biomolecular motors are often powered by the hydrolysis of adenosine triphosphate (ATP) into adenosine diphosphate (ADP) and an inorganic phosphate (Pi): , taking place inside their catalytic sites and generating conformational changes of the motor protein complex. Colloidal motors made of metal catalyst often use hydrogen peroxide as fuel. The free energy is provided by aqueous reaction,
, taking place on their metal catalytic surfaces. With hydrogen peroxide fuel, micrometric colloidal particles can reach speeds of the order of 10
m/s [Citation23,Citation24]. Higher speeds of about 20
m/s have been observed [Citation25], using hydrazine as fuel and iridium as catalyst.
These surface reactions, as well as diffusiophoresis, surface tension, and hydrophobicity, are interfacial properties that can be tuned by changing the materials composing the interface. Nonequilibrium thermodynamics provides a way to identify the material and interfacial properties that are involved in the mechanisms of propulsion. The next section is devoted to this key issue.
3. Fluctuating chemohydrodynamics
3.1 Catalytic particle moving in a reactive solution
The goal is to find equations of motion for an active colloidal particle and its surrounding fluid environment that account for thermal fluctuations as well as the processes that give rise to propulsion by diffusiophoresis. For inactive particles, this task reduces to the well-known problem of determining the stochastic evolution equations for Brownian motion. A similar approach can be adapted to account for active motion, and we shall illustrate it below by considering Janus motors whose surfaces have catalytic and non-catalytic hemispheres.
The starting point of the calculation is the computation of the force and torque exerted by the fluid on the Janus particle using methods developed in the seventies for Brownian motion [Citation26,Citation27]. The fluctuating Navier-Stokes equations are linearized, which is justified since the flow is laminar around micrometric particles, and they are solved using the method of the induced force densities that efficiently takes the boundary conditions into account. We now outline how this method can be extended to include coupling to concentration fields so that it is able to describe diffusiophoretic effects (or thermophoretic effects from temperature gradients) on the motion of the Janus particle.
The force exerted on the Janus particle by the fluid is given by the surface integral of the fluctuating pressure tensor at the interface
between the fluid and the Janus particle and, if present, an external force
, e.g. caused by a gravitational field. As a consequence, Newton’s equation for the Janus particle takes the form,
where is the mass of the Janus particle and
its mass density.
In a similar manner, a torque is exerted by the fluid on the Janus particle so that the angular velocity obeys the equation,
where the inertia tensor of the Janus particle has the components
with
, where
is the position of the center of mass of the Janus particle, and
is an external torque [Citation28–Citation31]. The external torque may be due to an external magnetic field
exerted on a magnetic dipole
attached to the particle [Citation32], in which case
, or due to a gravitational field acting on the nonuniform mass density of the Janus particle [Citation33].
In order to compute the force and torque on the Janus particle, we must evaluate the surface integrals involving the fluid pressure tensor at the catalytic surface. This computation requires a knowledge of the fluctuating fluid equations, along with boundary conditions that prescribe how the fluid fields are coupled to the concentration fields on the particle surface. These equations must be consistent with the basic conditions for microscopic time reversibility and other symmetries of the problem. Nonequilibrium thermodynamics combined with the theory of stochastic processes provides a method that can be used to construct these equations.
3.2. Thermodynamics and fluctuations
In nonequilibrium thermodynamics, the densities, , of any quantities such as the concentrations
of particles of species
, the mass density
, the linear momentum density
where
is the fluid velocity, the energy density
including the internal energy density
and the kinetic energy density
, as well as the entropy density
, obey balance equations with the general form
where are the corresponding current densities defined with respect to the barycentric motion and
are source densities [Citation34–Citation36].
Thermal fluctuations are incorporated in such continuum descriptions of the dynamics by adding noise terms to the mean currents
for the different irreversible processes
(the index
denoting the type and components of the scalar, vector or tensorial quantities):
[Citation20,Citation21,Citation37,Citation38] The mean currents are assumed to satisfy linear phenomenological laws
in terms of linear response coefficients and the thermodynamic forces also called affinities
corresponding to the currents. The linear dependence (4) remains valid as long as the nonequilibrium driving processes take place over scales larger than the molecular mean free paths. According to microreversibility, the linear response coefficients obey the Onsager–Casimir reciprocal relations
where
when
is even or odd under time reversal [Citation35, Citation39–Citation44]. Only coefficients
that couple processes with the same parity under time reversal contribute to the entropy production rate,
where is the Boltzmann constant. The fluctuating currents
are assumed to be Gaussian white noise processes characterized by
on time scales longer than molecular correlation times [Citation20,Citation21,Citation37,Citation38]. The second equations in (6) are called fluctuation–dissipation relations [Citation45]. Notice that they vanish if processes with opposite parities under time reversal are coupled together by coefficients such that . In this case, there is no associated noise to consider.
This general stochastic formulation may be applied to a Janus motor propelled by a diffusiophoretic mechanism. We suppose that the motor is suspended in a multi-component fluid containing solute species, labeled by the index , that interact with the motor through short-range intermolecular potentials
.
We further assume that the reversible reactions occur on the catalytic hemisphere of the motor, and call species
the fuel and
the product. (Generalizations to other reaction schemes such as those mentioned in Sec. 2 can be carried out.) These chemical reactions produce inhomogeneous
and
concentration fields in the motor vicinity that lead to a body force on the motor. If no external forces act on the system and momentum is conserved, fluid flows arise in the surrounding medium and are responsible for motor propulsion. The forms that the stochastic equations for the fluid velocity and concentration fields take in the solution and on the surface are discussed below.
3.3. Stochastic equations in the bulk phases
The bulk fluid phase equations are well known. If the exothermicity of the reaction is negligible, we may suppose that the system remains isothermal with an invariant and uniform temperature . The hydrodynamic and diffusive processes in the solution surrounding the catalytic Janus particle are described by the coupled Navier-Stokes and diffusion equations. The fluid is assumed to be incompressible,
, so that the mass density remains uniform. The fluctuating Navier-Stokes equations for the velocity field
are given by
where is the hydrostatic pressure, and
is the viscous part of the pressure tensor,
, where
is the shear viscosity, which is related to the corresponding Onsager coefficient by
. The Gaussian white noise fields
are characterized by
and
.
The fluctuating diffusion equations for the concentration fields of the different solute species
have the form
where the current densities can be expressed as , in terms of the molecular diffusivity
of species
and Gaussian white noise fields
satisfying
and
, with
denoting the
identity matrix [Citation21,Citation38]. The correlation functions are written to account for the fact that the affinities associated with the diffusive part of current densities are given by the gradients of chemical potentials
so that the Onsager coefficients are given by
[Citation35]. Moreover, the noise terms on the tensorial pressure and vectorial current densities are uncorrelated as a consequence of the Curie symmetry principle [Citation21,Citation35]:
.
Inside a solid Janus particle of radius , the velocity and concentration fields take the values,
, for
.
3.4. Stochastic equations at the interface
Any density can be decomposed as
, in terms of the densities
in the two bulk phases on both sides of the interface and the excess surface density
. Here
are Heaviside functions for the two bulk phases and the surface Dirac distribution
restricts the quantity it multiplies to the interface. Given this decomposition, stochastic equations that are consistent with microscopic reversibility can be written for interfacial quantities.
Several irreversible processes take place at the interface between the fluid and the Janus particle. First, there is the reaction with stoichiometric coefficients
and
. This reaction takes place on the catalytic hemisphere of the Janus particle with local rate,
where are the rate constants that are positive on the chemically active surface and zero elsewhere, and
is the interfacial noise associated with the surface reaction that satisfies
and
, where
and
are restricted to the interface by the delta distribution
. We note that the rate constants
have the units of m/s. Equivalently, we can use the rate constants
along with the diffusive rate constants
, which have the units of m
/s.
Second, there is a coupling between the frictional force along the interface associated with the partial slip of the velocity field between the fluid and the solid particle, and the interfacial diffusive transport of the excess surface densities . Diffusiophoresis results from the coupling between these two processes [Citation42,Citation43,Citation46,Citation47]. We have
where is a unit vector normal to the interface,
, the velocity slip between the fluid and the solid particle is given by
, for
, while
and
are Gaussian white noise terms. Equation (10) provides a boundary condition on the
component of the velocity field at the interface with the Janus particle. In addition to this boundary condition, the normal component of the velocity field obeys
also at the interface [Citation48].
Equations (10) and (11) are simplifications of more general results [Citation49], in which we assumed that the boundary layer and neighboring bulk phase are in local equilibrium so that the surface chemical potentials take their bulk values, . In this case, the surface Onsager coefficients are related to the coefficient of sliding friction,
, and the surface diffusion coefficients,
, while the surface Onsager coefficients relating the slip velocity and the surface current density are given by
where the
are diffusiophoretic constants and are given by
Here the slip length is and
, where
is the interfacial thickness,
is the interaction potential between species
and the colloid surface, and
is the inverse temperature[Citation50–Citation52].
In terms of these Onsager coefficients the Gaussian white noise terms and
satisfy the conditions
,
and
,
, where
and
are here also restricted to the interface. Since the interfacial slippage has an affinity that is odd under time reversal while the gradients of interfacial chemical potentials are even, the Onsager–Casimir reciprocal relations imply that the diffusiophoretic coupling terms do not appear in the expression for the excess surface entropy production, and the corresponding noise terms are uncorrelated,
. (We may make an alternative choice taking
as the current and
that is even under time reversal as the associated affinity. With this choice, the linear response coupling coefficients obey the standard symmetric Onsager reciprocal relations, and the diffusiophoretic coupling terms do appear in the expression for the excess surface entropy production.)
The boundary conditions on the concentration fields are determined by the transport of species in the direction normal to the interface, the surface reaction rate (9), as well as the reciprocal effect of diffusiophoresis. For species , the boundary condition can be expressed as
also for . Here
is a sink into the boundary pool of adsorbed species with
the surface velocity, and
the surface current density (11), which is determined, in particular, by the reciprocal effect of diffusiophoresis.
3.5. Langevin equations
Using the above results for the fluctuating fluid fields along with the boundary conditions for these fields on the surface of the particle, computation of the surface averages yields generalized Langevin equations for the force and torque.
For simplicity, we consider here time scales that are longer than the characteristic time of sound, , so that the fluid is effectively incompressible. In addition, on time scales longer than the hydrodynamic time,
, friction may be supposed to be time independent and given by the Stokes friction coefficient. This is the case for micrometric particles. Under these conditions, the Langevin equations take the form [Citation53],
where is the translational friction coefficient for arbitrary slip, and the diffusiophoretic force is
expressed in terms of the surface average . In this equation, we kept the deterministic diffusiophoretic force, so that the fluctuating force
also includes any possible noisy contributions from the concentration field in addition to those from fluid fluctuations.
Similarly, from the computation of the torque exerted by the fluid on the particle, we find the Langevin equation for the angular velocity,
where is the time-independent rotational friction coefficient,
is the deterministic diffusiophoretic torque and is a random torque. The diffusiophoretic constants
are included in the surface averages since they may be non-uniform on the particle surface in general.
4. Spherical Janus particle
4.1. Diffusion and reaction
For a micron-size Janus particle moving at the velocity in an aqueous solution where the solute molecular diffusion coefficients are of order
, the Péclet numbers,
, take values smaller than unity. In this regime, the diffusion Equation (8) reduces to
because the term
is relatively proportional to the Péclet number while the noise term goes as its square root. Discarding the noise and sink term
, the boundary conditions (13) on the concentration fields become
with where the Heaviside function
takes the values
on the catalytic hemisphere and
on the chemically inactive hemisphere, and
are the uniform concentrations in the solution far from the particle. We see that there is no coupling to the velocity field so that the concentration fields can be determined at this level of approximation independently of the velocity field. This approximation is justified because diffusion is fast enough around a micrometric particle for the concentration profiles to adjust themselves with respect to the orientation of the particle. Indeed, the time scales for rotational or translational diffusion,
and
, respectively, are longer than that for molecular diffusion,
.
The concentration fields of species around a spherical Janus particle of radius
are given by
where the function satisfies the diffusion equation
, which can be solved by expressing it in a series of Legendre polynomials,
, with the boundary conditions
and
in spherical coordinates
aligned parallel to the particle axis. The boundary condition at the particle surface
involves the dimensionless Damköhler number,
. In the reaction-limited regime
and
, while in the diffusion-controlled regime
and
.
Examples of concentration profiles for the fuel A are shown in along the axis of the Janus particle. We see the depletion of fuel near the catalytic hemisphere as the Damköhler number increases from the reaction- to the diffusion-limited regime. In this latter regime, we have that
near the catalytic hemisphere. Accordingly, the relative fuel concentration behaves as
and reaches the value
in the limit
. In the same limit, the relative product concentration
also reaches the value
. We note that these concentration fields are obtained by averaging over the fluctuations.
Figure 1. Relative concentration profiles of the fuel A along the axis of the Janus particle for three values of the Damköhler number. C and N denote the catalytic and non-catalytic hemispheres, respectively. The profiles are obtained for
.
![Figure 1. Relative concentration profiles of the fuel A along the axis z of the Janus particle for three values of the Damköhler number. C and N denote the catalytic and non-catalytic hemispheres, respectively. The profiles are obtained for cˉB=0.](/cms/asset/0425543b-0267-44da-83bd-8515f49acb4a/tapx_a_1602480_f0001_b.gif)
The reaction rate is given by in terms of the fluctuating local rate (9), where the surface integral is carried out over the catalytic hemisphere of the Janus particle because the rate constants
vanish on the non-catalytic hemisphere. Using the above results, its mean value may be written as
with where
is the first coefficient of the expansion of
, which behaves approximately as
with the Damköhler number.
At thermodynamic equilibrium, the concentrations should satisfy the Guldberg-Waage condition,
, where
is the standard free energy of the reaction
. The free energy of the reaction is related to the concentrations by
.
The reaction is driven out of equilibrium if the concentrations are not in their equilibrium ratio. In this respect, the nonequilibrium control parameter of the reaction is defined in general as the dimensionless chemical affinity,
which is positive (resp. negative) for the reaction running in the direction (resp.
), and vanishes at equilibrium since
.
In the following, we consider the reaction in the linear regime close to equilibrium where the deviations of the concentrations from their equilibrium values, are small:
. In this regime, the chemical affinity can be approximated as
, up to terms of second order in the concentration deviations
. Introducing the reaction diffusivity
associated with the reaction rate (20), the chemical affinity may also be written close to equilibrium as
, up to terms with higher powers in the reaction rate.
4.2. Fluid flow
Once the concentration fields are known, the velocity field can be calculated to get the force and the torque on the particle. For micrometric particles moving at self-diffusiophoretic velocity m/s in an aqueous solution of viscosity
, the fluid flow is laminar. Indeed, the Reynolds number is much smaller than unity,
, taking the kinematic viscosity
. In this regime, the Navier-Stokes Equations (7) reduce to the Stokes equations
because the left-hand side of (7) is relatively proportional to the Reynolds number, while the noise term behaves as the square root of the Reynolds number.
Using the methods of [Citation54] and [Citation55], for stick boundary conditions (), and uniform diffusiophoretic constants
and
, the velocity field is given by the solution of Stokes equations as [Citation56–Citation57],
where ,
is the unit vector in the direction of the Janus particle axis,
are the associated Legendre polynomials of degree
and order
,
, and
are again the coefficients of the expansion of
. An example of the velocity field (22) is depicted in [Citation56].
Figure 2. Streamlines of the velocity field (22) around a Janus particle in the laboratory frame where the particle moves with the velocity for
,
, and
from [Citation56]. The vertical axis
corresponds with
of the Janus particle. The upper hemisphere is catalytic. The motor here behaves as a pusher with
. The blue arrow denotes the direction of particle motion and black arrows the fluid flow.
![Figure 2. Streamlines of the velocity field (22) around a Janus particle in the laboratory frame where the particle moves with the velocity Vsd=2Υa1/3 for Re≃0.013, Pe≃0.17, and Da≃6.2 from [Citation56]. The vertical axis z corresponds with u of the Janus particle. The upper hemisphere is catalytic. The motor here behaves as a pusher with Υa2<0. The blue arrow denotes the direction of particle motion and black arrows the fluid flow.](/cms/asset/9fb57ee6-f535-4f9e-b593-07e5543b7816/tapx_a_1602480_f0002_oc.jpg)
At large distances, the velocity field (22) behaves as the Stokes dipole
which is a puller if and a pusher if
. We notice that this velocity field goes as
and thus contributes to the pressure tensor. In contrast, if there are concentration gradients
but no reaction, the velocity field is exactly given by
where is the diffusiophoretic velocity [Citation50,Citation58]. In this passive case, the velocity field (24) decreases as
faster than in the active case (23).
We recall that, if an external force is exerted on the colloidal particle, the velocity field is given by
where is the particle velocity. Here, the velocity field decreases as
. The behavior as
for the particle driven by active self-diffusiophoresis is thus intermediate between its driving by an external force and by passive diffusiophoresis. We note that
in all the cases.
4.3. Self-diffusiophoretic force and torque
Since the concentration fields generated by the reaction are given by Equation (19), we may write an explicit formula for the deterministic part of the force (15), giving the following self-diffusiophoretic force for arbitrary slip if the diffusiophoretic constants are uniform on the spherical surface of the Janus particle:
Accordingly, this force is oriented in the direction of the Janus particle axis . It is proportional to the diffusiophoretic constants (12), and also to the reaction rate (20), which expresses the fact that this force is active. An important aspect is that this force can be reversed depending on the concentrations of fuel and product molecules in the vicinity of chemical equilibrium. This feature is shared with molecular motors [Citation59].
For the same particle the diffusiophoretic torque (17) vanishes by cylindrical symmetry, so that and only the frictional torque due to viscosity remains in Equation (16).
We note that these results can be extended to Janus particles with diffusiophoretic constants taking different values on both hemispheres [Citation60].
5. Stochastic dynamics of Janus particle self-propulsion
5.1. Coupled overdamped Langevin equations
For micrometer-sized particles inertia is unimportant and henceforth we consider the overdamped limits of the Langevin equations. Letting , the Langevin Equation (14) for a spherical Janus particle reduces to
in terms of the self-diffusiophoretic velocity , the translational diffusion coefficient of the particle
, and the fluctuating velocity
satisfying
and
.
Using the expression (26) for the self-diffusiophoretic force and noting that the self-diffusiophoretic velocity can be written as in terms of the unit vector
, for a spherical Janus particle
is given by
In the last equality, we expressed the velocity in terms of the mean reaction rate (20) and the diffusiophoretic parameter, . The self-diffusiophoretic velocity (28) and the parameter
remain finite in the limits of perfect stick (
) and perfect slip (
) boundary conditions.
Since there is no self-diffusiophoretic torque for a spherical Janus particle, the overdamped limit of the rotational Langevin Equation (16) is given by
The rotational diffusion coefficient is related to the rotational friction coefficient by . Therefore, because the particle is spherical, this equation does not depend on the particle position or the reactive state (if the external torque is spatially uniform). Consequently, this stochastic equation is autonomous; thus, it drives the direction
independently of what happens for translation and reaction.
In addition to these two equations, we must also consider the stochastic equation for reaction since the concentration fields depend on reactive processes on the surface of the particle. The number of reactive events during the time interval
since the beginning of observation satisfies the stochastic differential equation:
where, in addition to the mean rate , we have a contribution
from self-diffusiophoresis to be determined. The fluctuating rate
satisfies
and
. The rate (30) can be written as
in terms of the numbers of molecules A and B in the solution.
Equations (27) and (30) are coupled equations for the currents of the variables
that are associated with the corresponding mechanical and chemical affinities
. We can use the Onsager symmetry principle in order to determine
. Since the variables
and
are even under time reversal, the Onsager coefficients must satisfy
. In view of this property, the coupled equations must take the form,
where the vector of the noise terms is denoted by . Onsager symmetry dictates that
. From the fluctuation–dissipation relations (6), we also have the property,
, for the noise correlations.
Examples of random trajectories obtained from simulations of these coupled overdamped Langevin equations are shown in . In (left), we see the self-propelled motion of the Janus particle by consumption of fuel in the absence of an external force. The particle is oriented by an external magnetic field in the -direction. In (right), an external force is exerted that is opposed to the direction of self-propulsion and large enough in magnitude so that the reaction is reversed and the reactant is synthesized from the product, instead of being consumed.
Figure 3. Janus particle subjected to an external force and magnetic field oriented in the -direction: The random rescaled displacement
, the number of reactive events
, and the particle orientation
, versus the rescaled time
for parameter values
,
, and
and (left) a zero external force, and (right) a rescaled external force equal to
.
![Figure 3. Janus particle subjected to an external force and magnetic field oriented in the z-direction: The random rescaled displacement z∗=zDr/Dt, the number of reactive events n∗=nDr/Drxn, and the particle orientation uz, versus the rescaled time t∗=Drt for parameter values βμB=2, Wrxn/DrxnDr=0.8, and χDrxn/Dt=0.8 and (left) a zero external force, and (right) a rescaled external force equal to f=βFDt/Dr=−2.5.](/cms/asset/282d98a8-511a-47dc-abb1-490b5e863624/tapx_a_1602480_f0003_oc.jpg)
If the external force and the magnetic field are oriented in the -direction,
and
, the particle is oriented on average in the same direction:
, while the averages
and
. Depending on the values of the mechanical and chemical affinities, the mean velocity
and rate
can take positive, vanishing, or negative values. By definition, the mean velocity vanishes at the stall force
, which is proportional to the self-diffusiophoretic force
. Also, the mean reaction rate is equal to zero at the force
. These two conditions are depicted in that shows the mean values of the velocities and rate as a function of the rescaled external force
for a positive value of the reaction rate
[Citation61]. In this case, the propulsion driven by the reaction exerts a mechanical work if the force is in the range
, corresponding to the domain I in . If the force is sufficiently opposed to propulsion to satisfy
, the mean reaction rate can become negative in the domain II of fuel synthesis from the product.
Figure 4. Janus particle subjected to an external force and magnetic field oriented in the -direction [Citation61]: The mean values of the fluctuating rescaled velocities
and rate
versus the rescaled magnitude of the external force
for the parameter values
,
, and
. The dots show the results of a numerical simulation with a statistics of
trajectories integrated over the time interval
.
denotes the rescaled stall force and
the threshold between fuel synthesis and consumption. The Janus particle is propelled against the external force in the interval I. Fuel synthesis happens in the interval II.
![Figure 4. Janus particle subjected to an external force and magnetic field oriented in the z-direction [Citation61]: The mean values of the fluctuating rescaled velocities r˙∗=r˙/DtDr and rate n˙∗=n˙/DrxnDr versus the rescaled magnitude of the external force f=βFDt/Dr for the parameter values βμB=2, Wrxn/DrxnDr=0.8, and χDrxn/Dt=0.8. The dots show the results of a numerical simulation with a statistics of 105 trajectories integrated over the time interval t∗=10. fstall denotes the rescaled stall force and f0 the threshold between fuel synthesis and consumption. The Janus particle is propelled against the external force in the interval I. Fuel synthesis happens in the interval II.](/cms/asset/50d0b6c6-7a36-4392-8813-e94828fa8e62/tapx_a_1602480_f0004_oc.jpg)
5.2. Comparison with microscopically reversible particle-based simulations
These theoretical predictions were compared with particle-based simulations of Janus particle dynamics [Citation62]. The Janus particle is constructed from catalytic and non-catalytic beads that interact with the fluid species by repulsive intermolecular potentials while the motions of the fluid particles were described by multiparticle collision dynamics [Citation63–Citation65]. The reactive and nonreactive dynamics of the system satisfied microscopic reversibility, and the nonequilibrium conditions were established by reservoirs with constant concentrations far from the particle. The specification of the intermolecular potential and multiparticle collision parameters determined the transport properties of the system. As in the Langevin simulations, an external force along
was applied and the magnetic field in the
-direction controlled the orientation of the Janus particle.
The results of these simulations for and
are shown in as a function of the external force. The fits to these graphs indicated by the lines are in accord with the theoretical predictions. In particular, we note that the reaction rate varies with the external force as predicted by the theory,
. For the given parameters,
while the fit yields
. The results of another simulation, with parameters chosen so that the system is at equilibrium in the absence of an external force, are also presented in . Since the diffusiophoretic parameter
is nonzero, an external force can change the reaction rate. This is seen in the figure, and we also see that for negative values of the external force the Janus particle converts the product to fuel as predicted by the theory.
Figure 5. Simulations with microscopically reversible kinetics of the Janus motor subjected to an external force and a magnetic field both oriented in the
-direction: Plots of the dependence on the rescaled force
of the rescaled average motor velocity in the
-direction,
(left), and of the rescaled reaction rate,
(right) for systems with (a)
and
(circles) and (b)
and
(squares). The results for (a) and (b) systems were obtained from averages over 200 and 100 realizations of the dynamics, respectively. The fits to the data are indicated by (a) upper and (b) lower lines. See Huang et al. [Citation62], for additional information.
![Figure 5. Simulations with microscopically reversible kinetics of the Janus motor subjected to an external force Fext and a magnetic field both oriented in the z-direction: Plots of the dependence on the rescaled force f of the rescaled average motor velocity in the z-direction, ⟨z˙∗⟩ (left), and of the rescaled reaction rate, ⟨n˙∗⟩ (right) for systems with (a) cˉA=10 and cˉB=9 (circles) and (b) cˉA=10 and cˉB=10 (squares). The results for (a) and (b) systems were obtained from averages over 200 and 100 realizations of the dynamics, respectively. The fits to the data are indicated by (a) upper and (b) lower lines. See Huang et al. [Citation62], for additional information.](/cms/asset/514f8a0c-77b2-4ece-8bef-6f277bbe8606/tapx_a_1602480_f0005_oc.jpg)
We note that nonlinear dependencies due to corrections in powers of the Péclet number manifest themselves for larger values of the external force magnitude.
5.3. Thermodynamic efficiency
Similar to molecular motors [Citation59], the efficiency of the mechanical power of the Janus motor can be characterized by , and the efficiency of the chemical process of synthesis by
. Since the thermodynamic entropy production rate of the coupled processes
is non-negative according to the second law of thermodynamics, the mechanical and chemical efficiencies are bounded by and
in their respective domains of application.
The stall force where mean motor velocity vanishes and the force
where the mean reaction rate vanishes are depicted in the
-
plane of the mechanical and chemical affinities in . Domain I in this figure corresponds to affinity values where self-diffusiophoretic propulsion occurs as a result of fuel consumption, while in domain II fuel is synthesized from the product.
Figure 6. Schematic representation of the different regimes of the active particle in the plane of the mechanical and chemical affinities for [Citation53]. In domain I, self-diffusiophoretic mechanical work is powered by the reaction. In domain II, an external force of sufficient magnitude acting in a direction opposite to that of the Janus particle velocity can yield the synthesis of fuel from the product. For
, the slopes of the lines
and
are instead positive.
![Figure 6. Schematic representation of the different regimes of the active particle in the plane of the mechanical and chemical affinities for χ>0 [Citation53]. In domain I, self-diffusiophoretic mechanical work is powered by the reaction. In domain II, an external force of sufficient magnitude acting in a direction opposite to that of the Janus particle velocity can yield the synthesis of fuel from the product. For χ<0, the slopes of the lines ⟨z˙⟩=0 and ⟨n˙⟩=0 are instead positive.](/cms/asset/134372d6-270d-487f-918f-654cd5fe58d5/tapx_a_1602480_f0006_oc.jpg)
For a given chemical affinity, the maximal value of the mechanical efficiency is given by . Accordingly, the efficiency of self-propulsion increases as the square of the diffusiophoretic coupling
. A similar expression holds for the chemical efficiency [Citation61]. The locations where the efficiencies reach their maximal values are depicted as dashed lines in . The efficiency of self-propelled motors has been studied in [Citation66], and in [Citation67], by evaluating in detail the dissipation due to viscous friction in the surrounding solution.
5.4. The Fokker–Planck equation and its consequences
The Fokker–Planck equation that follows from the coupled overdamped Langevin Equations (29) and (31) and governs the time evolution of the joint probability density can be written as
where the rotational diffusion operator is
with a rotational energy associated with the external torque exerted on the particle, for instance, due to an external magnetic field.
If we suppose that there is no external force and integrate over the position
and the orientation
of the Janus particle, we obtain the equation,
, for the time evolution of the probability
that
reactive events have happened during the time interval
. Therefore, the probability starting from
at
is Gaussian
describing a random walk with a drift at the reaction rate and the diffusivity
.
A similar result holds if the external force is non-vanishing. Indeed, integrating over time the stochastic differential equation for in Equation (30) with an external force
in the
-direction and calculating the mean value and variance of
, we find the following effective rate and diffusivity
in terms of the mean value of the orientation along the -direction and the reorientation time
, with
. By the central limit theorem, the probability
is given after a long enough time by the Gaussian (35) with the effective quantities instead of the bare ones.
Similarly, if we suppose that there is an external force in the
-direction while the reaction rate vanishes
and integrate over the variables
,
,
, and
, we find that the probability density to observe the particle with the position
at time
is also Gaussian
here describing a random drift of mean velocity and diffusivity
. Now, if the reaction rate is not equal to zero, we can integrate the stochastic differential Equation (27) over time and calculate the mean value and variance of
to get the following effective drift velocity and diffusivity
also in terms of . Again, by the central limit theorem, the probability density
becomes after a long enough time the Gaussian (37) with the effective quantities instead of the bare ones.
In the absence of an external magnetic field , there is no preferential orientation
, so that the rotational motion remains diffusive and controlled by the rotational diffusion time
. In this case, we recover the known result [Citation22] that the effective diffusion coefficient is given by
which is quadratic in the reaction rate and thus also in the self-diffusiophoretic velocity
. The larger the magnitude of the self-diffusiophoretic velocity, the more enhanced the diffusive random walk of the Janus particle. Furthermore, the velocity autocorrelation function is given in this case by
which justifies the use of colored noise models for active particles [Citation68–Citation70].
6. Fluctuation theorems and microreversibility
The fluctuation theorems are recent results of nonequilibrium statistical mechanics, showing that the ratio of the probabilities of opposite fluctuations in the particle or energy transfers across an open system are related by
to the affinities or thermodynamic forces driving the system out of equilibrium [Citation71–Citation76]. These relations hold independently of the details of the process and arbitrarily far from thermodynamic equilibrium beyond the regimes of linear response. At equilibrium, the affinities are equal to zero
and the ratio of the probabilities goes to one, so that we recover the principle of detailed balance. The fluctuation theorems are established on the basis of the time-reversal symmetry of the underlying microscopic dynamics. Remarkably, they imply that the entropy production
is always non-negative in agreement with the second law of thermodynamics.
6.1. Chemical fluctuation theorem
A fluctuation theorem can be obtained for the Gaussian probability distribution (35). Indeed, we have that
with the chemical affinity , in the case without an external force when the mechanical affinity is equal to zero
.
As shown in , the probabilities and
to find the number
of reactive events at time
shift away from the origin
in opposite directions with the rates
. Consequently, their overlap rapidly decreases as time
increases. However, their ratio continues to vary as
in the long-time limit. This remarkable property explains that the probability of observing an opposite fluctuation quickly becomes very small, so that the process is observed to be irreversible, although the opposite fluctuations are still present.
6.2. Mechanical fluctuation theorem
Using the Gaussian probability distribution (37), we have a similar fluctuation relation for the mechanical displacement ,
with the mechanical affinity , in the case where
.
If , by using the central limit theorem with the effective quantities (38), we recover the effective fluctuation relation [Citation77], for the displacement along the
-direction
which is expressed in terms of an effective force resulting from the external and self-diffusiophoretic forces, and the effective temperature
. The relation (45) reduces to Equation (44) in the absence of chemical reaction. Here, we see that fluctuation theorems can be used to justify the introduction of effective temperatures on the basis of nonequilibrium driving mechanisms.
6.3. Mechanochemical fluctuation theorem
The previous results can be generalized to the joint probability of a displacement
of the motor after
reactive events have occurred during the time interval
, which obeys the mechanochemical fluctuation theorem
This result has been proved in [Citation61], using the method of the cumulant generating function by modifying the Fokker–Planck equation to include the counting variables. The fluctuation theorem (46) implies the non-negativity of the entropy production (32). A fundamental result is that it also implies the Onsager symmetry of Equation (31) in consistency with microreversibility.
shows that the mechanochemical fluctuation theorem is satisfied. This theorem is also valid in the nonlinear regime away from equilibrium for reactions more complicated than those considered here, where it predicts generalizations of the Onsager reciprocal relations to the higher cumulants and their responses [Citation72,Citation76].
Figure 8. Janus particle subjected to an external force and magnetic field oriented in the -direction [Citation61]: (a) Probability density
with
versus the rescaled displacement
at the rescaled time
for the parameter values
,
,
, and
. (b) Verification of the mechanochemical fluctuation relation (46) in the same conditions. The probability ratio is calculated if
and
are larger than
. The dots are the results of a numerical simulation with an ensemble of
trajectories and an integration with the time step
. The lines depict the theoretical expectations.
![Figure 8. Janus particle subjected to an external force and magnetic field oriented in the z-direction [Citation61]: (a) Probability density P(z∗,n∗;t∗) with n∗=nDr/Drxn=2,1,0,−1,−2 versus the rescaled displacement z∗=zDr/Dt at the rescaled time t∗=Drt=1 for the parameter values βμB=1, βFDt/Dr=−1, Wrxn/DrxnDr=0.4, and χDrxn/Dt=0.4. (b) Verification of the mechanochemical fluctuation relation (46) in the same conditions. The probability ratio is calculated if P(z∗,n∗;t∗) and P(−z∗,−n∗;t∗) are larger than 10−4. The dots are the results of a numerical simulation with an ensemble of 107 trajectories and an integration with the time step dt∗=10−3. The lines depict the theoretical expectations.](/cms/asset/9f5b9e7e-3911-43e0-9bba-282e635a23e4/tapx_a_1602480_f0008_oc.jpg)
6.4. Finite-time fluctuation theorem
In general, the aforementioned fluctuation theorems hold in the long-time limit. They may also hold at finite times for the Gaussian processes of Equations (35) or (37). Remarkably, this is also the case for more complicated processes involving linear reactions such as possibly coupled to diffusion [Citation78–Citation80]. Finite-time fluctuation theorems can be established for diffusion-influenced surface reactions on immobile catalysts. These processes are ruled by the fluctuating diffusion Equations (8) with the surface reaction rate (9), but a fluid at rest (
). In such systems, we have that
with a time-dependent affinity
that converges towards
in the long-time limit. Analytic expressions can be obtained for
by solving a linear diffusion-reaction problem with special boundary conditions.
shows the time-dependent affinity for the diffusion-influenced surface reaction on an immobile Janus particle for three different values of the Damköhler number. We see that the convergence towards the asymptotic value is slower in the diffusion-limited than in the reaction-limited regime.
Figure 9. Time-dependent affinity of the finite-time fluctuation theorem for an immobile Janus particle of radius
with surface reaction influenced by the diffusion of fuel and product from a spherical reservoir at the distance
from the particle center. The asymptotic value of the affinity is equal to
(dotted line). The diffusion coefficients take the value
and the rate constants
. The Damköhler number is given by
.
![Figure 9. Time-dependent affinity At of the finite-time fluctuation theorem for an immobile Janus particle of radius R=1 with surface reaction influenced by the diffusion of fuel and product from a spherical reservoir at the distance L=10 from the particle center. The asymptotic value of the affinity is equal to Arxn=1 (dotted line). The diffusion coefficients take the value D=DA=DB=1 and the rate constants κ=κ+=κ−. The Damköhler number is given by Da=2κR/D.](/cms/asset/013d94d1-7ddb-4487-a1fd-3e2012f6d087/tapx_a_1602480_f0009_b.gif)
7. Comparison with biomolecular motors and enzymes
7.1. Efficiencies and fluctuation theorems for biomolecular motors
Energy transduction and mechanochemical coupling are general features common to colloidal and biomolecular motors [Citation81]. Indeed, that shows the efficiencies of a self-diffusiophoretic motor in the linear regime near equilibrium is comparable to similar diagrams obtained for molecular motors in the linear or nonlinear regimes [Citation59,Citation82–Citation84]. In these latter cases, the different domains in the plane of chemical and mechanical affinities are separated by curved instead of straight lines due to the kinetic nonlinearities. Otherwise, for molecular motors, there also exists a domain where propulsion is powered by the reaction and another domain where fuel is synthesized by the external force, as in for the self-diffusiophoretic motor.
Furthermore, mechanochemical fluctuation theorems hold for biomolecular motors[Citation82,Citation84,Citation85]. Accordingly, the ratio of probabilities of opposite fluctuations in the movement and reaction of the motors behaves as
in terms of the mechanical and chemical affinities, and
, which have similar definitions as for the self-diffusiophoretic motor. These theorems are the expression of microreversibility, and they have consequences for the Onsager reciprocal relations of the linear response coefficients, as well as their generalizations involving the nonlinear response coefficients [Citation72,Citation73,Citation76].
These results show the analogy between self-diffusiophoretic active particles and molecular motors. Similar results also hold for enzymes.
7.2. Enhanced diffusion linear in the reaction rate
According to Equation (39), the enhancement of diffusion by the self-diffusiophoretic mechanism is quadratic in the reaction rate. We may wonder if there is a mechanism that is instead linear in the reaction rate. Here, we consider propulsion by the release in a funnel of several product molecules from the fuel molecule. As these molecules move away from the reaction center, the phase-space domain widens in the funnel as the translational, rotational, and vibrational degrees of freedom of the molecules are freed. In this regard, the available phase-space volume increases as with the number
of degrees of freedom and the distance
from the reaction center in the funnel. Therefore, the free-energy potential has an entropic contribution going as
and the force exerted on the motor is
[Citation86–Citation88]. During the release, the motor is thus propelled according to
, giving the displacement
over the time interval
, where
is the diffusion coefficient corresponding to Stokes’ friction coefficient
. If the reaction proceeds at the rate
, this displacement is repeated every time step
on average, generating the propulsion velocity
. Because of rotational diffusion at the rate
, the effective diffusion coefficient is thus given by
since the ratio of translational to rotational diffusion coefficients is equal to in terms of Stokes’ hydrodynamic radius
of the motor. Therefore, this mechanism leads to a linear dependence on the reaction rate instead of the quadratic dependence of Equation (39).
8. Collective motion
Some of the most fascinating properties of self-propelled particles are associated with their collective behavior. These include active self-assembly processes and a variety of dynamical cluster states, and there is a growing literature on this topic for systems of active particles [Citation32,Citation89–Citation93]. We provide a sketch of how the collective behavior of diffusiophoretic Janus motors can be described.
Systems of self-propelled colloidal particles can experience different kinds of interactions: (1) Direct interparticle interactions already manifest themselves at equilibrium where they determine the phase diagram of colloidal systems [Citation94]. They play no role in dilute suspensions, but they become increasingly important as the colloidal density increases. (2) Chemotactic interactions arise from nonequilibrium gradients of species concentration fields. It is important to distinguish between diffusiophoretic effects due to concentration gradients in the absence of reaction and self-diffusiophoretic effects in the presence of reaction. (3) Hydrodynamic interactions between particles are mediated by the velocity field. They have a nonequilibrium origin since they involve the gradients of the velocity field; also hydrodynamic flows accompany the force-free dynamics of self-diffusiophoretic particles. It has been shown that the nature of the collective dynamics of diffusiophoretic motors depends strongly on chemotactic interactions, and hydrodynamic interactions also play an important role [Citation95]. Simple motility-induced-phase-separation models, [Citation69,Citation96–Citation98], do not incorporate these effects and new mechanisms involving concentration gradients and hydrodynamic coupling operate to determine the forms that the collective dynamics take. In keeping with the theme of this review, we confine our attention to the collective behavior of systems containing colloidal motors that move by diffusiophoretic mechanisms.
Here, we consider an isothermal dilute suspension of colloidal motors moving in a dilute solution of fuel A and product B molecular species. At lowest order, the motors are isolated from each other and they move in a dilute solution of the molecular species. As in the Chapman-Enskog expansion for the solution of Boltzmann’s equation, we may suppose that there exist local concentration gradients of fuel and product species and local velocity gradients of the fluid. Consequently, we should first solve the problem where every colloidal particle moves in concentration and velocity fields that are not uniform far from the particle, but manifest concentration and velocity gradients. To simplify the discussion, we assume that hydrodynamic effects play a negligible role and focus on the nature of the concentration fields.
8.1. A colloidal particle moving in concentration gradients
First, we have to solve the same problem as in Subsec. 4.1, but with the boundary condition at a large distance from the particle, instead of
. The solution may again be expressed as in Equation (19), but with a different function
now depending on the concentration gradients
. The force and the torque exerted on the Janus particle in such concentration fields can be calculated by Equations (15) and (17). In this way, we obtain the diffusiophoretic velocity and angular velocity,
The concentration gradients have been written as
because the concentration fields
in the presence of many colloidal particles may differ on large scale from the local concentration profiles
around every particle and; moreover, these gradients are no longer uniform on large scales as supposed with the notation
, but are given instead by spatial derivatives of non-uniform fields
.
In order to have non-vanishing torque coefficients, , we suppose that the diffusiophoretic constants
are not uniform on the surface of the particle, so that Equation (17) now gives the non-zero angular velocity (50). Moreover, this non-uniformity also contributes to the other coefficients
,
, and
.
If the diffusiophoretic constants are uniform on the particle surface, we have
and no angular velocity,
. Then, in a uniform background of concentrations, we have
and the velocity reduces to the active self-diffusiophoretic velocity
. Instead, at chemical equilibrium where
, the self-diffusiophoretic contributions to the force vanish,
and
, and there remains the passive diffusiophoretic velocity,
with
. Out of chemical equilibrium, the coefficients
are non-vanishing and there are active contributions of the reaction to the other coefficients
and
as well.
8.2. Ensemble of colloidal motors
A dilute suspension of colloidal motors moving in a dilute solution of fuel A and product B molecular species may be described with the distribution function of the colloidal motors, , expressed in terms of 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 (33) for the probability that a single motor is located at the position
with the orientation
by using the diffusiophoretic velocity (49) and angular velocity (50) to get
in terms of effective translational and rotational diffusion coefficients and
, and the operator (34) describing rotational diffusion of the colloidal motor in the effective rotational energy
, including the contribution of the diffusiophoretic torque in addition to the contribution from an external torque due to a magnetic field
acting on a magnetic dipole
of the Janus particle.
Once this equation is established, we can obtain the equations for the successive moments of , including the density of colloidal particles
and the polarizability or polar order parameter of the colloidal motors
. Higher moments can be neglected if only the smallest isotropies are considered. These equations are coupled to the effective reaction-diffusion equations for the molecular species.
Such sets of coupled equations already manifest clustering instabilities, even without the coupling to hydrodynamics [Citation99,Citation100]. In order to study hydrodynamic interactions, we need to couple the previous diffusion-reaction equations to the Navier-Stokes equations for the velocity field and understand how these latter are modified by the presence of the colloidal particles. It is already known that the viscosity coefficient is modified according to Einstein’s formula and corrections [Citation101], but we also need to determine how the activity of the motors affects the velocity field. The methods developed since Einstein’s 1906 paper, [Citation102–Citation107], can be extended to determine not only the dissipative but also the active contributions to the pressure tensor. Several of them have already been obtained [Citation108], but methods are available to deduce them systematically. Collective dynamics of such diffusiophoretic motors has been studied in several types of system [Citation95,Citation109–Citation111].
9. Conclusion and perspectives
Thermodynamics and statistical mechanics provide fundamental approaches to study the propulsion mechanisms of self-diffusiophoretic motors. The coupling between motion and reaction required for energy transduction should be consistent with the underlying microreversibility of the processes considered. This led to phenomena such as the prediction that the reaction rate for diffusiophoretic motors depends on the external force, and formed the basis for the fluctuation theorems for these systems. The phenomena described here should be observable in experiments.
The same general considerations enter the description of collective behavior. As noted above, the equations for a single self-diffusiophoretic motor in an external concentration gradient can be constructed to be consistent with microscopic reversibility, guaranteeing that all coupling contributions appear in the evolution equations. These equations may then be used as input for the construction of consistent mean-field equations for the collective dynamics. In addition, the way in which the nonequilibrium state is established by coupling to reservoirs is an essential part of a dynamically consistent formulation of collective dynamics. The present methods provide a way to systematically deduce the equations of motion of these assemblies of colloidal motors, allowing their precise engineering and aid in the development of a full theory that goes beyond a mean-field formulation.
The study also shows that the design of colloidal motors requires engineering material interfacial properties, besides the shape of the colloids. These properties, which are determined by the coating of the colloid, can be systematically identified with equilibrium and nonequilibrium interfacial thermodynamics. They include the diffusiophoretic constants, together with the surface reaction rates, the slip length characterizing hydrophobicity, or the surface tension.
The approach can be extended to non-spherical particles, nonlinear surface reactions, electrodiffusiophoresis, or thermophoresis. Among the other perspectives, one can envisage situations where the medium surrounding the active particle is a rarefied or dilute gas, instead of a liquid solution, as described by the Boltzmann equation, instead of the Navier-Stokes and diffusion equations. The effects of long-time tails on diffusiophoresis could be investigated in such situations [Citation112]. The study of these effects concerns transport and reaction of atmospheric aerosols [Citation113] or particles in interstellar clouds. The medium can also be a plasma, a viscoelastic fluid, or a liquid crystal, having different interfacial properties involved in the propulsion mechanism.
Acknowledgments
The Authors thank Patrick Grosfils and Mu-Jie Huang for fruitful discussions. Financial support from the International Solvay Institutes for Physics and Chemistry, the Université libre de Bruxelles (ULB), the Fonds de la Recherche Scientifique - FNRS under the Grant PDR T.0094.16 for the project “SYMSTATPHYS”, and the Natural Sciences and Engineering Research Council of Canada is acknowledged.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Kay ER, Leigh DA, Zerbetto F. Synthetic molecular motors and mechanical machines. Angew. Chem. Int. Ed. 2007;46:72–191.
- Sauvage JP, Gaspard P, editors. From non-covalent assemblies to molecular machines. Weinheim: Wiley-VCH; 2011.
- van Leeuwen PWNM. Homogeneous catalysis. Dordrecht: Kluwer Academic Publishers.; 2004.
- Dey KK, Pong FY, Breffke J, et al. Dynamic coupling at the Ångstrom scale. Angew Chem Int Ed. 2016;55:361–397.
- Alberts B, Bray D, Johnson A, et al. Essential cell biology: an introduction to the molecular biology of the cell. New York: Garland Science; 1998.
- Jones RAL. Soft machines: nanotechnology and life. Oxford: Oxford University Press; 2004.
- Wang J. Nanomachines: fundamentals and applications. Weinheim, Germany: Wiley-VCH; 2013.
- Mareschal M, Malek Mansour M, Puhl A, et al. Molecular dynamics versus hydrodynamics in a two-dimensional Rayleigh-Bénard system. Phys Rev Lett. 1988;61:2550–2553.
- Paxton WF, Kistler KC, Olmeda CC, et al. Catalytic motors: autonomous movement of striped nanorods. J Am Chem Soc. 2004;126:13424–13431.
- Fournier-Bidoz S, Arsenault AC, Manners I, et al. Synthetic self-propelled nanorotors. Chem Commun. 2005;441–443.
- Wang W, Duan W, Ahmed S, et al. Small power: autonomous nano- and micromotors propelled by self-generated gradients. Nano Today. 2013;8:531–554.
- Sánchez S, Soler L, Katuri J. Chemically powered micro- and motors. Angew Chem Int Ed. 2014;54:1414–1444.
- Colberg PH, Reigh SY, Robertson B, et al. Chemistry in motion: tiny synthetic motors. Acc Chem Res. 2014;47:3504–3511.
- Abdelmohsen L, Peng F, Tu Y, et al. Micro- and nano-motors for biomedical applications. J Mater Chem B. 2014;2:2395.
- Rückner G, Kapral R. Chemically powered nanodimers. Phys Rev Lett. 2007;98:150603.
- Sumino Y, Magome N, Hamada T, et al. Self-running droplet: emergence of regular motion from nonequilibrium noise. Phys Rev Lett. 2005;94:068301.
- Chen YJ, Nagamine Y, Yoshikawa K. Self-propelled motion of a droplet induced by Marangoni-driven spreading. Phys Rev E. 2009;80:016303.
- Thutupalli S, Seemann R, Herminghaus S. Swarming behavior of simple model squirmers. New J Phys. 2011;13:073021.
- Izri Z, van der Linden MN, Michelin S, et al. Self-propulsion of pure water droplets by spontaneous Marangoni-stress-driven motion. Phys Rev Lett. 2014;113:248302.
- Landau LD, Lifshitz EM. Statistical physics, part 2. Oxford: Pergamon Press; 1980.
- Ortiz de Zárate JM, Sengers JV. Hydrodynamic fluctuations in fluids and fluid mixtures. Amsterdam: Elsevier; 2006.
- Kapral R. Perspective: motors without moving parts that propel themselves in solution. J Chem Phys. 2013;138:020901.
- Valadares LF, Tao YG, Zacharia NS, et al. Catalytic motors: self-propelled sphere dimers. Small. 2010;6:565–572.
- Ke H, Ye S, Carroll RL, et al. Motion analysis of self-propelled Pt-silica particles in hydrogen peroxide solutions. J Phys Chem A. 2010;114:5462–5467.
- Gao W, Pei A, Dong R, et al. Catalytic iridium-based Janus micromotors powered by ultralow levels of chemical fuels. J Am Chem Soc. 2014;136:2276–2279.
- Mazur P, Bedeaux D. A generalization of Faxén’s theorem to nonsteady motion of a sphere through an incompressible fluid in arbitrary flow. Phys A. 1974;76:235–246.
- Bedeaux D, Mazur P. Brownian motion and fluctuating hydrodynamics. Phys A. 1974;76:247–258.
- Hills BP. A generalized Langevin equation for the angular velocity of a spherical Brownian particle from fluctuating hydrodynamics. Phys A. 1975;80:360–368.
- Felderhof BU. Force density induced on a sphere in linear hydrodynamics: I. Fixed sphere, stick boundary conditions. Phys A. 1976;84:557–568.
- Felderhof BU. Force density induced on a sphere in linear hydrodynamics: II. Moving sphere, mixed boundary conditions. Phys. A 1976;84:569–576.
- Bedeaux D, Albano AM, Mazur P. Brownian motion and fluctuating hydrodynamics II; A fluctuation-dissipation theorem for the slip coefficient. Phys A. 1977;88:574–582.
- Bechinger C, Leonardo RD, Löwen H, et al. Active particles in complex and crowded environments. Rev Mod Phys. 2016;88:045006.
- Campbell AI, Ebbens SJ. Gravitaxis in spherical Janus swimming devices. Langmuir. 2013;295:14066–14073.
- Prigogine I. Introduction to thermodynamics of irreversible processes. New York: Wiley; 1967.
- de Groot SR, Mazur P. Nonequilibrium thermodynamics. New York: Dover; 1984.
- Nicolis G. Irreversible thermodynamics. Rep Prog Phys. 1979;42:225–268.
- Landau LD, Lifshitz EM. Statistical physics, part 1. 3rd ed. Oxford: Pergamon Press; 1980.
- Gardiner CW. Handbook of stochastic methods for physics, chemistry and the natural sciences. 3rd ed. Berlin: Springer; 2004.
- Onsager L. Reciprocal relations in irreversible processes I. Phys Rev. 1931;37:405–426.
- Onsager L. Reciprocal relations in irreversible processes II. Phys Rev. 1931;38:2265–2279.
- Casimir HBG. On Onsager’s principle of microscopic reversibility. Rev Mod Phys. 1945;17:343–350.
- Waldmann L. Non-equilibrium thermodynamics of boundary conditions. Z Naturforschg A. 1967;22:1269–1280.
- Bedeaux D, Albano AM, Mazur P. Boundary conditions and non-equilibrium thermodynamics. Phys A. 1976;82:438–462.
- Haase R. Thermodynamics of irreversible processes. New York: Dover; 1969.
- Callen HB, Welton TA. Irreversibility and generalized noise. Phys Rev. 1951;83:34–40.
- Kovac J. Non-equilibrium thermodynamics of interfacial systems. Phys A. 1977;86:1–24.
- Bedeaux D. Nonequilibrium thermodynamics and statistical physics of surfaces. Adv Chem Phys. 1986;64:47–109.
- Albano AM, Bedeaux D, Mazur P. On the motion of a sphere with arbitrary slip in a viscous incompressible fluid. Phys A. 1975;80:89–97.
- Gaspard P, Kapral R. Nonequilibrium thermodynamics and boundary conditions for reaction and transport in heterogeneous media. J Chem Phys. 2018;148:194114.
- Anderson JL. Colloid transport by interfacial forces. Ann Rev Fluid Mech. 1989;21:61–99.
- Anderson JL, Prieve DC. Diffusiphoresis caused by gradients of strongly absorbing solutes. Langmuir. 1991;7:403–406.
- Ajdari A, Bocquet L. Giant amplification of interfacially driven transport by hydrodynamic slip: diffusio-osmosis and beyond. Phys Rev Lett. 2006;96:186102.
- Gaspard P, Kapral R. Fluctuating chemohydrodynamics and the stochastic motion of self-diffusiophoretic particles. J Chem Phys. 2018;148:134104.
- Lighthill MJ. On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers. Commun Pure Appl Math. 1952;5:109–118.
- Blake JR. A spherical envelope approach to ciliary propulsion. J Fluid Mech. 1971;46:199–208.
- Reigh SY, Huang MJ, Schofield J, et al. Microscopic and continuum descriptions of Janus motor fluid flow fields. Phil Trans R Soc A. 2016;374:20160140.
- Campbell AI, Ebbens SJ, Illien P, et al. Experimental observation of flow fields around Janus spheres. arXiv:1802.04600. 2018.
- Anderson JL. Transport mechanisms of biological colloids. Ann N Y Acad Sci. 1986;469:166–177.
- Jülicher F, Ajdari A, Prost J. Modeling molecular motors. Rev Mod Phys. 1997;69:1269.
- Huang MJ, Schofield J, Gaspard P, et al. From single particle motion to collective dynamics in Janus motor systems. J Chem Phys. 2019;150:124110.
- Gaspard P, Kapral R. Mechanochemical fluctuation theorem and thermodynamics of self-phoretic motors. J Chem Phys. 2017;147:211101.
- Huang MJ, Schofield J, Gaspard P, et al. Dynamics of Janus motors with microscopically reversible kinetics. J Chem Phys. 2018;149:024904.
- Malevanets A, Kapral R. Mesoscopic model for solvent dynamics. J Chem Phys. 1999;110:8605–8613.
- Kapral R. Multiparticle collision dynamics: simulation of complex systems on mesoscales. Adv Chem Phys. 2008;140:89–146.
- Gompper G, Ihle T, Kroll DM, et al. Multi-particle collision dynamics: a particle- based mesoscale simulation approach to the hydrodynamics of complex fluids. Adv Polym Sci. 2009;221:1–87.
- Tao YG, Kapral R. Dynamics of chemically powered nanodimer motors subject to an external force. J Chem Phys. 2009;131:024113.
- Sabass B, Seifert U. Dynamics and efficiency of a self-propelled, diffusiophoretic swimmer. J Chem Phys. 2012;136:064508.
- Farage TFF, Krinninger P, Brader JM. Effective interactions in active Brownian suspensions. Phys Rev E. 2015;91:042310.
- Marconi UMB, Gnan N, Paoluzzi M, et al. Velocity distribution in active particles systems. Sci Rep. 2016;6:23297.
- Wittmann R, Maggi C, Sharma A, et al. Effective equilibrium states in the colored-noise model for active matter I. Pairwise forces in the Fox and unified colored noise approximations. J Stat Mech Exp. 2017;11:113207.
- Lebowitz JL, Spohn H. A Gallavotti-Cohen-type symmetry in the large deviation functional for stochastic dynamics. J Stat Phys. 1999;95:333–365.
- Andrieux D, Gaspard P. A fluctuation theorem for currents and non-linear response coefficients. J Stat Mech Theory Exp. 2007;2007:P02006.
- Andrieux D, Gaspard P. Fluctuation theorem and Onsager reciprocity relations. J Chem Phys. 2004;121:6167–6174.
- Jarzynski C. Equalities and inequalities: irreversibility and the second law of thermodynamics at the nanoscale. Annu Rev Condens Matter Phys. 2011;2:329–351.
- Seifert U. Stochastic thermodynamics, fluctuation theorems and molecular machines. Rep Prog Phys. 2012;75:126001.
- Gaspard P. Multivariate fluctuation relations for currents. New J Phys. 2013;15:115014.
- Falasco G, Pfaller R, Bregulla AP, et al. Exact symmetries in the velocity fluctuations of a hot Brownian swimmer. Phys Rev E. 2016;94:030602(R).
- Andrieux D, Gaspard P. Temporal disorder and fluctuation theorem in chemical reactions. Phys Rev E. 2008;77:031137.
- Gaspard P, Kapral R. Finite-time fluctuation theorem for diffusion-infuenced surface reactions. J Stat Mech Exp. 2018;2018:083206.
- Gaspard P, Grosfils P, Huang MJ, et al. Finite-time fluctuation theorem for diffusion-influenced surface reactions on spherical and Janus catalytic particles. J Stat Mech Exp. 2018;2018:123206.
- Hill TL. Free energy transduction and biochemical cycle kinetics. New York: Dover; 2005.
- Lacoste D, Lau AWC, Mallick K. Fluctuation theorem and large deviation function for a solvable model of a molecular motors. Phys Rev E. 2008;78:011915.
- Astumian RD. Thermodynamics and kinetics of molecular motors. Biophys J. 2010;98:2401–2409.
- Gerritsma E, Gaspard P. Chemomechanical coupling and stochastic thermodynamics of the F1-ATPase molecular motor with an applied external torque. Biophys Rev Lett. 2010;5:163–208.
- Andrieux D, Gaspard P. Fluctuation theorems and the nonequilibrium thermodynamics of molecular motors. Phys Rev E. 2006;74:011906.
- Zwanzig R. Diffusion past an entropy barrier. J Phys Chem. 1992;96:3926–3930.
- Reguera D, Rubí JM. Kinetic equations for diffusion in the presence of entropic barriers. Phys Rev E. 2001;64:061106.
- Rubí JM, Lervik A, Bedeaux D, et al. Entropy facilitated active transport. J Chem Phys. 2017;146:185101.
- Theurkauff I, Cottin-Bizonne C, Palacci J, et al. Dynamic clustering in active colloidal suspensions with chemical signaling. Phys Rev Lett. 2012;108:268303.
- Elgeti J, Winkler RG, Gompper G. Physics of microswimmers–single particle motion and collective behavior: a review. Rep Prog Phys. 2015;78:056601.
- Wang W, Duan W, Ahmed S, et al. From one to many: dynamic assembly and collective behavior of self-propelled colloidal motors. Acc Chem Res. 2015;48:1938–1946.
- Zöttl A, Stark H. Emergent behavior in active colloids. J Phys Condens Matter. 2016;28:253001.
- Liebchen B, Marenduzzo D, Pagonabarraga I, et al. Clustering and pattern formation in chemorepulsive active colloids. Phys Rev Lett. 2015;115:258301.
- Jones RAL. Soft condensed matter. Oxford UK: Oxford University Press; 2002.
- Huang MJ, Schofield J, Kapral R. Chemotactic and hydrodynamic effects on collective dynamics of self-diffusiophoretic Janus motors. New J Phys. 2017;19:125003.
- Cates ME, Tailleur J. When are active Brownian particles and run-and-tumble particles equivalent? consequences for motility-induced phase separation. EPL. 2013;101:20010.
- Bialké J, Löwen H, Speck T. Microscopic theory for the phase separation of self-propelled repulsive disks. EPL. 2013;103:30008.
- Speck T, Menzel AM, Bialké J, et al. Dynamical mean-field theory and weakly non-linear analysis for the phase separation of active Brownian particles. J Chem Phys. 2015;142:224109.
- Saha S, Golestanian R, Ramaswamy S. Clusters, asters, and collective oscillations in chemotactic colloids. Phys Rev E. 2014;89:062316.
- Pohl O, Stark H. Self-phoretic active particles interacting by diffusiophoresis: a numerical study of the collapsed state and dynamic clustering. Eur Phys J E. 2015;38:93.
- Bedeaux D, Kapral R, Mazur P. The effective shear viscosity of a uniform suspension of spheres. Phys A. 1977;88:88–121.
- Einstein A. Investigations on the theory of the Brownian movement. New York: Dover; 1956.
- Landau LD, Lifshitz EM. Fluid mechanics. 2nd ed. Oxford: Pergamon Press; 1987.
- Lebenhaft JR, Kapral R. Diffusion-controlled processes among partially absorbing stationary sinks. J Stat Phys. 1979;20:25–56.
- Happel J, Brenner H. Low Reynolds number hydrodynamics. The Hague: Martinus Nijhoff Publishers; 1983.
- Dhont JKG. An introduction to dynamics of colloids. Amsterdam: Elsevier; 1996.
- Nägele G. Colloidal hydrodynamics. In: Bechinger C, Sciortino F, Ziherl P, editors. Physics of complex colloids, proc. of the int. school of physics “Enrico Fermi”. Amsterdam: IOS Press; 2013. p. 507–601.
- Ramaswamy S. The mechanics and statistics of active matter. Annu Rev Condens Matter Phys. 2010;1:323–345.
- Thakur S, Kapral R. Collective dynamics of self-propelled sphere-dimer motors. Phys Rev E. 2012;85:026121.
- Colberg P, Kapral R. Many-body dynamics of chemically-propelled motors. J Chem Phys. 2017;147:064910.
- Robertson B, Huang MJ, Chen JX, et al. Synthetic motors: working together through chemistry. Acc Chem Res. 2018;51:2355–2364.
- Li T, Raizen MG. Brownian motion at short time scales. Ann Phys. 2013;525:281–295.
- Li W, Davis EJ. The effects of gas and particle properties on thermophoresis. J Aerosol Sci. 1995;26:1085–1099.