A self-propelled thermophoretic microgear †

An asymmetric microgear will spontaneously and unidirectionally rotate if it is heated in a cool surrounding solvent. The resulting temperature gradient along the edges of the gear teeth translates in a directed thermophoretic force, which will exert a net torque on the gear. By means of computer simulations, the validity of this scenario is proved. The rotational direction and speed are dependent on gear – solvent interactions, and can be analytically related to system parameters like the thermal di ﬀ usion factor, the solvent viscosity, or the temperature di ﬀ erence. This microgear provides a simple way to extract net work from non-isothermal solutions, and can become a valuable tool in micro ﬂ uids.


I. Introduction
Molecular motors usually refer to biological systems that operate in an environment where thermal uctuations are signicant, and that have an internal mechanism to convert energy into directed motion or mechanical work.These motors are ubiquitous in nature and play a crucial role in the transport within biological organisms.Examples are motor proteins moving along laments 1 or bacteria swimming in low Reynolds number environments. 2 Such biological active matter has also been used to design nanomachines, as is the case of microgears pushed by a bacterial bath. 3,4Inspired by nature, purely synthetic micromotors display self-propelled motion without biological components.These motors are receiving a rapidly increasing attention as a very promising avenue to realize purposeful functions on the micro and nanoscales.5][16][17][18][19][20] Phoresis refers to the directed dri motion that suspended particles experience in inhomogeneous conditions.Important examples of such inhomogeneities are gradients of temperature (thermophoresis), concentration (diffusiophoresis), and electric potential (electrophoresis).In cases where the gradients are locally produced by the particles themselves, self-propulsion can occur.This is the case of a thermophoretic swimmer realized by laser heating a colloidal sphere half metalcoated. 18Due to their simplicity and controllability, a benchmark investigation of the properties of active colloids has been experimentally performed with phoretic swimmers. 21,22Phoretic micromotors have been designed until now by considering heterogeneous surface properties, as is the case of Janus particles and heterodimers to make swimmers, of twin and tethered Janus particles to build a rotor, 18 or of partially coated gears. 14o nd alternative designs of phoretic motors in general, and with homogeneous surfaces in particular, is challenging from a fundamental viewpoint, and has a great technological interest.
In this paper, we show that an asymmetric microgear with homogeneous surface properties rotates when heated in a cool surrounding solvent.The speed and direction of the microgear rotation are determined by its geometry, the interactions with the solvent, and the applied temperature differences.This can be experimentally realized by heating an asymmetric microgear with larger thermal conductivity than the solvent.Our results provide a novel route to design phoretic micromotors with homogeneous surfaces, which can be fueled by local heating.

II. Mesoscopic model
A hybrid simulation scheme in two dimensions is employed to model both the microgear and the surrounding solvent, whose typical time and length scales are separated by orders of magnitude.The nature of the problem imposes the necessity of reproducing stable temperature gradients, the effect of hydrodynamic interactions, and the conservation of angular momentum. 23A particle-based mesoscopic simulation technique known as multiparticle collision dynamics (MPC) [24][25][26][27][28] is employed for the solvent, while the microgear and its interactions with the solvent are simulated by standard molecular dynamics (MD).The MPC dynamics consists of alternating streaming and collision steps.In the streaming step, the solvent particles of mass m move ballistically for a time h.In the collision step, particles are sorted into a square lattice with lattice size a, and interchange momentum relative to the center-ofmass velocity of each collision cell.In our simulations the stochastic rotation collision rule with variable collision angle a introduced by Ryder and Yeomans 29,30 is employed.This collision rule locally conserves mass, linear momentum, energy and angular momentum.It can therefore be proved that the algorithm properly captures hydrodynamic interactions, thermal uctuations, and the sustainability of temperature inhomogeneities.Simulation units are reduced by setting a ¼ 1, m ¼ 1, and k B T ¼ 1 with k B the Boltzmann constant and T a reference temperature.We employ standard MPC parameters h ¼ 0.1 and the mean number of solvent particles per cell r ¼ 10.From the kinetic theory of the MPC with a xed rotation angle a ¼ 90 , 31 we can approximately calculate the Schmidt number as Sc ¼ 12, which corresponds to a liquid-like behavior.In our simulations, the kinematic viscosity is obtained to be n ¼ 0.46, which is 40% smaller than the previous estimation.
The considered microgear is a solid structure where the surface is a sequence of sawteeth in a closed circular shape (Fig. 1).In our simulations, a gear with 8 teeth is used, with an internal radius R 1 ¼ 19a and an external radius R 2 ¼ 25a.The short edge of each sawtooth is in the radial direction such that the tooth has angles q 1 ¼ 40 and q 2 ¼ 90 .The microgear is surrounded by MPC solvent which is conned inside a circular wall with radius R w ¼ 45a.To obtain the solid gear structure two components are considered.One is a rigid gear with the sawteeth prole, with a momentum of inertia I ¼ 10 6 ma 2 .The rigid gear is free to rotate around its center xed at the center of the simulation setup.Then a single-layer of monomer beads is mounted along the edges of the rigid gear, where the separation between neighboring beads is a.Each bead is attached to the rigid gear by a harmonic spring of constant k ¼ 600k B T/a 2 .There are no further interactions between different beads.The external wall is similarly constructed by xing beads with springs along an external circle.The coupling of the microgear and the solvent takes place through the MD bead-solvent particle interactions.The employed interaction is a Lennard-Jones (LJ) type potential 32,33 Here r is the distance between the bead center and the solvent particle, 3 refers to the potential intensity, s to the bead radius, and n to a positive integer describing the potential stiffness.The attractive or repulsive LJ potentials are obtained respectively by taking c ¼ 0 or c ¼ 3 with the corresponding cutoff r c .The bead radius is taken as s ¼ 1.25a, and 3 ¼ k B T. For efficiently exchanging energy with the surrounding solvent, the considered bead mass is M ¼ m.A hard repulsive potential (n ¼ 24, c ¼ 3) is chosen for the external wallsolvent interactions, while both repulsive and attractive potentials are considered for the microgear-solvent interactions.Note that given the large overlap between neighboring beads (the separation between beads is 0.4 times their diameter) the solvent particles (not shown in Fig. 1) remain conned between the microgear and the circular wall.The equations of motion are integrated with a velocity-Verlet algorithm and a time step The simulated microgear temperature T g is uniformly imposed by independently thermostatting every bead in the gear edges every ten MD steps with a Maxwellian velocity distribution of temperature T g , which is similar to the Andersen thermostat. 34The thermosttating operation violates the conservation of the microgear angular momentum, which is then restored by adding or subtracting the corresponding small overall angular momentum.This compensation does not affect the microgear rotation, since the angular momentum variation in the thermostat operation slightly uctuates around zero.Moreover, energy is drained from the system by thermostatting the wall beads with temperature T w .By imposing the gear temperature increment DT ¼ T g À T w , a steady-state temperature distribution is quickly established (Fig. 2a).The environment of the solvent particles close to the summit and the cle of each gear tooth is quite different (different size of the heating areas), such that the solvent temperature is different in both positions and varies along the edges.8][39] This temperature discontinuity could enhance the geometry-induced temperature gradient along the edges.We refer to VT l and VT s as the temperature gradients along the long and the short edges of each gear tooth (Fig. 2b).To the leading order, the gradients are expected to be proportional to the gear temperature increment, e.g.|VT l | ¼ l 1 |DT|, with l 1 a positive coefficient determined by the solid-solvent coupling and the gear geometry.In the example shown in Fig. 2c, l 1 x (120a) À1 .In the radial direction, the temperature varies logarithmically, as shown in Fig. 2d, which is a consequence of the conservation of energy.
The simulations performed here enforce the microgear constant temperature.Experimentally this corresponds to a microgear fabricated with a material of thermal conductivity much higher than that of the solvent, as it would be the case of a metal or a metal-coated microgear in water solution.However, the heat transport within the microgear is disregarded in our simulations given that the temperature is imposed by the use of a local thermostat.This is not relevant for our purpose, since the way in which the microgear constant temperature is imposed does not affect the solvent temperature distribution, nor the solvent-gear interactions, and hence nor the gear motion.Alternatively, thermophoretic microgears can be constructed with materials of low or moderate thermal conductivity, and simulated with bead-bead interactions.Such microgears will not display a homogeneous temperature distribution, but a central temperature higher than that at the gear edges.The temperature at each summit will be lower than the temperature at the cles.As a result, the temperature gradient of the solvent along the gear edge is still present, which is the crucial point for the motion of the hot microgear.Depending on the material properties, this temperature gradient can in principle be larger or smaller than in the case of the gear with constant temperature, which corresponds to a different value of l 1 and therefore to a different gear rotation speed.

III. Results and discussion
A qualitative discussion of the microgear motion is rst presented, followed by the simulation results, and a quantitative discussion.1][42] This force is generated by the inhomogeneous interactions of the surrounding solvent with the particle. 43When a solvent with a temperature gradient is in contact with a planar wall, only the tangential temperature gradient induces a thermophoretic force, which is then parallel to the wall.Thus, for our microgear VT l and VT s will respectively result in the thermophoretic forces on the long edges f T,l and the short edges f T,s , parallel to the edges.Depending on the gear-solvent interactions, 32 the thermophoretic forces can be along (thermophilic) or against (thermophobic) the temperature gradient.Due to the gear geometry, the thermophoretic forces exert a non-vanishing torque ) on the gear which results in its unidirectional rotation.Here, R l and R s respectively refer to the coordinates of the center of force on the long and short edges, and the summation accounts for multiple teeth.Besides the thermophoretic forces, the solvent exerts the standard pressure forces normal to the gear edges.
The pressure forces on the short and long edges of the gear teeth produce torques in opposite directions that cancel each other.Note that in our simulations, given the very large overlap between neighboring beads, the edges of the microgear can be regarded as at walls.Otherwise, for rough walls a normal temperature gradient will generate a thermophoretic force perpendicular to the wall.However, similar to the pressure forces, these normal forces would not result in any contribution to the net gear torque.Furthermore, in most of our simulations, the short tooth edge (hence R s ) and long tooth edge are, respectively, parallel and perpendicular to the radial direction.This means that only the forces applied in the long edges contribute to the total torque, and |T | x 8R 1 |f T,l |.A hot gear (DT > 0) built of a thermophilic material will then rotate with the long teeth edges forward, which in the geometry of Fig. 1 is clockwise.Rotation in the opposite direction is then expected in the case of a cold microgear (DT < 0), or when the material is thermophobic (see Fig. 3 and ref. 44).For a gear in equilibrium (DT ¼ 0) only pressure forces will be present, such that no net rotation is predicted.The microgear rotation in the simulations is characterized by measuring rotation angle 4 as illustrated in Fig. 1, where a positive 4 corresponds to a clockwise motion.Simulation results of a thermophilic gear show in Fig. 3a an example of forward rotation for a hot gear (DT > 0), backward rotation for a cold gear (DT < 0), and no rotation in the case of a non-heated gear (DT ¼ 0).The averaged quantities consider a minimum of 8 independent runs.Fig. 3c shows an instantaneous gear trajectory where the unidirectional rotation can be observed to be simultaneously accompanied by thermal uctuations (see also ref. 44).The thermophilic microgear is simulated by considering repulsive bead-solvent interactions (n ¼ 3, c ¼ 3). 32Hot microgears with thermophobic behavior show in Fig. 3b the expected anticlockwise rotation.These thermophobic gears are simulated by attractive interactions (c ¼ 0) of two different kinds of soness (n ¼ 6, and n ¼ 10).In all cases the self-induced rotation is due to the breakdown of the spatial symmetry produced by the asymmetric geometry of the heated microgear.We perform additional simulations for a microgear with symmetric teeth as displayed in the inset of Fig. 3e.The thermophoretic forces along both sides of each tooth are then symmetric with respect to the microgear radial direction, which results in a zero torque and vanishing net rotation (Fig. 3b).
In order to provide an expression for the rotation of the self-propelled microgear in terms of the material properties, the thermophoretic force on the long edge of the microgear needs to be explicitly calculated.For an isolated large suspended particle, f T is well-accepted to be proportional to the temperature gradient VT with the so-called thermodiffusion factor a T , [41][42][43] By denition a T > 0 corresponds to a thermophobic particle, and a T < 0 to a thermophilic particle.A bead embedded on the microgear interacts with the solvent only partially, such that its thermodiffusion factor a T,g can be related to that of the isolated bead a T,g ¼ l 2 a T , with the dimensionless correction factor l 2 (0 < l 2 < 1).Independent simulations with a single bead are performed to quantify a T for the different gear-solvent interactions used in Fig. 3e.This has been implemented by directly measuring the thermophoretic force on one isolated bead xed in a solvent with an externally imposed temperature gradient. 32,45The thermophoretic force on the long edge of the microgear then reads f T,l ¼ ÀN l l 2 a T k B VT l with N l ¼ 20 the number of beads on each long edge.Therefore, the effective thermodiffusion factor of the long edge is a ef T ¼ N l l 2 a T .This is a convenient concept when the constituent surface beads cannot be clearly identied, as in the case of most experimentally available systems.The torque exerted on the microgear is then with ẑ the unit vector towards the plane.Here a T is assumed to be temperature independent.The resulting angular velocity u ¼ m r T , is proportional to the microgear rotational mobility m r .
Regarding the microgear as a disk with hydrodynamic radius R H , the mobility can be identied by m r ¼ 1/(4phR H 2 ). 46The angular velocity is then The linear dependence of the rotation angle 4 with time shown in Fig. 3a-c allows us to quantify the angular velocity u of the gear in our simulations.The data in Fig. 3d and e are nicely consistent with linear dependence predicted by eqn (3) on DT and on a T .A quantitative comparison of our simulation results with eqn (3) is non-trivial, since we do not really have a reliable measurement of parameters l 2 and R H in eqn (3).In the case of the repulsive gear with DT ¼ 1.0, we have measured l 1 , and the thermal diffusion factor a T ¼ À1.0.The hydrodynamic radius can be considered to be the external gear radius R H x R 2 , which together with the measured solvent viscosity, and u, determines the factor l 2 x 0.1.This value is consistent with the fact that only 13% of the area of the microgear beads is in contact with the solvent.On the other hand it is important to note that, the essential mechanisms of this self-propelled microgear are the thermophoretic effect and the geometry-induced temperature gradient along the microgear edges, which are rather general and universal.Therefore, extensions of the model, like the consideration of the surface beads with internal degrees of freedom, and/or the gear with a temperature gradient inside (moderate heat conductivity), would cause only quantitative changes, leaving the essence of the device unchanged.
In order to emphasize the experimental feasibility and potential of the thermophoretic microgear, it is interesting to discuss a possible estimation of the orders of magnitude of the gear rotation speed u.The rotational mobility of a 3 dimensional microgear is calculated as m r ¼ 1/(4phR H 2 h H ), 46 with h H the gear thickness.The gear hydrodynamic radius R H , and the internal radius R 1 in eqn (2) will be of the same order of magnitude, such that the size dependence can be summarized as u $ a ef T /(R H h H ). The thermodiffusion factor a T is well-known to be strongly dependent on particle size in general, 40 and in particular for diluted spherical colloids, 47,48 such that a signicant dependence is also expected for at surfaces.A polystyrene particle with 1 mm diameter in water has been characterized with a T $ 5000. 49Although there are no available experimental data to determine the relationship between a T and a ef T , we can for example consider a gear with the radius R H ¼ 50 mm, thickness h H ¼ 1 mm, and then in a similar spirit to our simulations we assume l 2 ¼ 0.1 and N l ¼ 100.This would correspond to a linear increase of a ef T with the length of the long edge, which would also be about R H . Considering now the water viscosity h $ 0.001 kg ms À1 and a temperature gradient VT ¼ 0.1 K mm À1 , it could be concluded that the microgear rotates $1 round per second, which can be easily observed in experiments.
][52][53][54] This radiometer is driven by thermal creep and works therefore for rareed gases with typical sizes of millimeter.In contrast, the microgear presented here is driven by the thermophoretic effect in liquids, which is expected to work in microscales.The Crookes radiometer is built upon vanes with sides of different heat absorption, and therefore different temperatures.The rotation only happens on the cool side of the blade in the front.Meanwhile, the described thermophoretic microgear can rotate in both directions.Both the rotational direction and speed will change not only with the applied temperature increment but also with many other factors, related to the nature of the thermodiffusion factor.This factor is determined by composition of the gear and the uid contained between the walls, and it will be affected by additional substances diluted in the uid, or external conditions like pressure or average temperature.Moreover, the competition of thermophoresis with other effects like thermoelectricity 55,56 has interestingly shown the existence of materials whose properties vary from thermophobic to thermophilic.All these effects can provide a large versatility to this device.Correspondingly, the thermophoretic gear can become a very valuable tool to investigate the thermophoretic properties of a wide class of systems.Until now a requirement to determine the thermodiffusion factor, or equivalently the Soret coefficient, has been that the investigated system should be a solution.Therefore, materials systems that would for example precipitate in solution like gold in water could be investigated now by means of this new device.

IV. Conclusions
In conclusion, we have proved that an asymmetric microgear with homogeneous surface properties can spontaneously and unidirectionally rotate in solution, when the gear temperature is clearly distinguished from that of the environment.Our proof has been performed by means of computer simulations, but the effect is not restricted to the peculiarities of our model.Similar microrotor can also be realized experimentally, for example by heating a metallic microgear with a laser.When coupling the microrotor to an external device, net work could be extracted from nonisothermal solutions.As an example of other practical applications, our microrotor could be used as a stirring device, which could be locally controlled.Furthermore, when keeping the microgear xed, the reaction of the thermophoretic force on the tooth edges can result in the motion of the surrounding uid. 57This can be employed to construct a thermophoretic pump, whose uid motion is perpendicular to the applied temperature gradient. 58Our ndings provide an alternative strategy to design synthetic micromotors, which have become a promising tool in the eld of microuidics, and also in the investigation of thermal diffusion.

Fig. 1
Fig. 1 Simulation setup of the eight-teeth microgear within a circular bead wall.Parameters are described in the main text.

Fig. 2
Fig. 2 (a) Steady state distribution of the solvent temperature with gear and wall temperatures fixed as T g ¼ 1.5 T and T w ¼ 0.5 T, and gearsolvent soft repulsive interactions (n ¼ 3, c ¼ 3).(b) Zoom-in in the neighborhood of one gear tooth, indicating the temperature gradients VT l and VT s in the long and short edges.(c) Temperature along one of the long tooth edges, where r l is the distance from one summit to the corresponding cleft.(d) Temperature as a function of r w , the distance from the external wall to one of the tooth summits.Symbols in (c and d) correspond to simulation values and lines to data fits, in (c) the fit is linear and in (d) logarithmic.

Fig. 3
Fig. 3 (a-c) Rotational angle of the microgear as a function of time.Positive 4 refers to rotation with the long edge forward.Symbols are simulation results and lines linear fits.(a) Averaged angle for repulsive gear repul (n ¼ 3, c ¼ 3) and various temperature increments.(b) Averaged angle for DT ¼ 1.0 and various interactions, att-a (full squares) refers to soft attractive gear (n ¼ 6, c ¼ 0), att-b (open squares) to hard attractive gear (n ¼ 10, c ¼ 0), and sym (circles) to a gear with symmetry geometry as depicted in the inset of (e).(c) Single run with DT ¼ 1.0 and repul.(d and e) Angular velocity as a function of DT, and the thermodiffusion factor a T of a single surface bead in solution.Symbols correspond to the same parameters as in (a) and (b).