Controlled motion of Janus particles in periodically phase-separating binary fluids

We numerically investigate the propelled motions of a Janus particle in a periodically phase-separating binary fluid mixture. In this study, the surface of the particle tail prefers one of the binary fluid components and the particle head is neutral in the wettability. During the demixing period, the more wettable phase is selectively adsorbed to the particle tail. Growths of the adsorbed domains induce the hydrodynamic flow in the vicinity of the particle tail, and this asymmetric pumping flow drives the particle toward the particle head. During the mixing period, the particle motion almost ceases because the mixing primarily occurs via diffusion and the resulting hydrodynamic flow is negligibly small. Repeating this cycle unboundedly moves the Janus particle toward the head. The dependencies of the composition and the repeat frequency on the particle motion are discussed.


I. INTRODUCTION
Self-propelled motions of micro-and nano-particles have attracted much interest from a wide range of viewpoints.They will provide us with important applications, such as nanomachines and drug delivery [1][2][3][4][5][6][7][8].Recent, focus has been on their collective dynamics because they are very fascinating in a growing field of non-equilibrium physics, i.e., active matter [9][10][11][12][13].Self-propelled particles use some energy or nutrients to generate the selfpropulsion force.For example, biological molecules, such as ATPase and myosin, convert chemical energy to mechanical motion through chemomechanical coupling [14].In non-biological systems, Marangoni effect can induce spontaneous motions of liquid droplets [15,16].
Janus particles, which have heterogeneous surface properties, are often employed as artificial self-propelled systems [17,18].For example, the self-propelled motions are modelled by asymmetric nanoparticles partially coated with platinum.The catalytic decomposition of hydrogen peroxide, which occurs selectively on the Ptsurface, drives the nanoparticles [1,2,5,19,20].Interfacial phoretic effects [21,22] are another possible mechanism of micro-swimmers [17].Jiang et al. demonstrated that a Janus particle can create an asymmetric temperature gradient around it in a defocused laser beam.The induced gradients lead to spontaneous drift motions of the Janus particles [7].Also, local heating by illumination light induces active motions of Janus particles in a binary mixture of lower critical solution temperature [23,24].
Besides the thermo-and diffusio-phoresis motions [7,23,24], Janus particles can move in phase-separating binary mixtures because of the coupling between the wetting and phase separation [25,26].The phase separation of binary fluid mixtures has been well studied [27][28][29].During the later stage of the phase separation, the domain patterns grow with time.The particles in the phase-separating mixtures are trapped in one of the phases or at the interfaces.Even when Brownian motions and external forces are absent, the particles move with the resulting coarsening of the domain patterns [30][31][32][33].
Because the Janus particles have asymmetric wettability, we expect that their motions also become asymmetric.The direction of the asymmetric motion will more or less depend on the particle direction.Janus particles with two distinct wettabilities are occasionally used as surfactants to stabilize the phase-separated domains [34][35][36].
In the final stage of the phase separation, the particle motions will be frozen.
In this article, we demonstrate a possible mechanism of spontaneous motions of a Janus particle in periodically phase-separating mixtures.By continually varying the temperature or pressure slightly above and below the transition point, one can cause periodic processes of phase separation and mixing [37][38][39][40].By resetting the binary mixtures to the one-phase state, we expected that we could continuously propel the particle.Here, we examine this expectation by means of numerical simulations.The dependencies of the particle motion on the average composition and the duration of the cycle are discussed.
In Sec.II, we explain our numerical model, which is based on fluid particle dynamics method [41,42].Numerical results are shown and are discussed in Sec.III.We summarize our study and discuss some remarks in Sec.IV.

A. Free energy functional
We consider the case, in which a spherical Janus particle is suspended in a binary mixture.The surface of the particle is heterogeneous in wettability.For numerical simulations, dealing with the particle in a continuous manner is convenient.We express it with a smooth shape function as [41,42], Here, r is the coordinate in a lattice space and R is the position of the particle in an off-lattice space.a is the radius of the particle and d represents the width of the smooth interface.In the limit of d → 0, ψ is unity and zero in the interior and exterior of the particle, respectively.We also define the surface distribution as The free energy function comprises two parts as [32,43,44], ( The first part F mix is the mixing free energy given by where φ(r) is the local concentration of the first component of the mixture.T is the temperature and we set the Boltzmann constant to unity.v 0 is the molecular volume.
The coefficient of the gradient term C is related to the interface tension and is of the order of v 2/3 0 [29].f BW (φ) is the Bragg-Williams type of the mixing free energy as where χ is the interaction parameter between two fluid components.Under the mean field approximation, χ = 2 and φ = 0.5 give the critical point.The third term in the integrand of Eq. ( 3) is introduced to prevent the solvent from penetrating into the particle [32].χ p and φ 0 are the control parameters.
F sur is the surface free energy, which is given by where n is the unit vector toward the orientation of the particle.W represents the heterogeneity of the particle wettability.As illustrated in Fig. 1(a), we use where W 0 and W 1 are the material constants for the wetting.If W < 0, the component of larger φ tends to wet the surface [25,26].In the early stages of phase separation, the surface of W = 0 largely influences the pattern formation (see below).

B. Time development equations
The hydrodynamic equation for the flow field v is given by FIG. 1.(a) A sketch of a spherical particle with a heterogeneous surface.The orientation is described by a unit vector n.We set W0 = −1 and W1 = 1 in Eq. ( 6), so that we have W = 0 and W = −2 at the two poles.We refer to the pole with W = 0 as "head" and the other as "tail".The particle head is neutral in wettability and the tail prefers component A. (b) A typical phase diagram of a binary fluid mixture.We change the χ parameter by using a square wave function of t.Phase separation occurs during the period with χ = χ d .However, when χ = χm, the phase-separated domains are mixed.
Here ρ is the material density.In this work, we assume that all the materials have the same density.Σ is the viscous stress tensor, which is given by In the spirit of fluid particle dynamics (FPD), we assume that the viscosity η depends on the particle distribution as [41,42], where η 0 is the viscosity of the binary fluid and ∆η is the viscosity difference between the solvent and the particles.
In the limit of ∆η → ∞, the particles will behave as solid particles [42].p is a part of pressure, which imposes the incompressibility condition: f is the force field stemming from the particle interactions and is given by, where Ω is the effective volume defined as Ω = drψ R (r).It is approximated as Ω ≈ 4πa 3 /3 in threedimensional systems.g originates from the torque acting on the particle, and is given by The particle motions are caused by the hydrodynamic flow and their kinetics are described as Furthermore, the time development equation of the concentration field is where L(ψ) is the kinetic coefficient, and in which we set L(ψ) = L 0 (1 − ψ) to eliminate the flux inside the particles.L 0 is the kinetic coefficient of the bulk mixture.ζ represents the thermal fluctuation satisfying the fluctuation-dissipation relation.In this model, the diffusion flux does not contribute to the particle motion directly.In each state, the total free energy, including the kinetic energy, should decrease with time.Its temporal change is described in Appendix A. Our model can be applied to many-particle systems, where Janus particles behave as surfactants [34][35][36].We hope that we will report our studies on them elsewhere in the near future.

C. Numerical simulations
We numerically solve the above equations using the Maker and Cell method with staggered grid [45].We discretize the space by d and set v 0 = d 3 and C = d 2 .Also, we set the particle radius to a = 6d.The simulation box is a three-dimensional system (64 3 ) with periodic boundary conditions.The time increment is 0.005 t 0 , where t 0 is a typical diffusion time defined by t 0 = d 2 T /L 0 .For the wettability, we set W 0 = −1 and W 1 = 1 in this study.Hence, the head of the particle is neutral in the wettability and the tail prefers the component of large φ (see Fig. 1(a)).Hereafter, we describe the more and less wettable components as A and B, respectively.The viscosity parameters in Eq. ( 8) are η 0 = ρL 0 /T and ∆η = 49η 0 .In Eq. (3), we set φ 0 to be equal to the average composition, φ , and χ p = 20.Because Reynold numbers in colloidal systems are very small, we iterate to integrate Eq. ( 7) without updating R, n, and φ until |ρ(∂/∂t + v • ∇)v| becomes less than 10 −3 η 0 L 0 /(T 0 d 3 ).The intensity of the thermal fluctuation is given by ζ i (r, t)ζ j (r , t ) = 0.05T L(ψ)δ(r − r )δ(t − t )δ ij , where i and j stand for x, y and z.

D. Periodic phase separation
To induce periodic phase separation, we change the χ parameter uniformly in space using square wave functions of time (see Fig. 1(b)).t m and t d denote the durations for mixing and demixing, respectively.We also define t t = t m + t d .For simplicity, because we consider mixtures near the phase-separation point, we assume that the other physical parameters are constant, independent of χ.In durations of mt t ≤ t < mt t + t m , χ is fixed at χ m below the critical point.Here, m is an integer.The phase-separated domains will be mixed during these times.In mt t + t m ≤ t < (m + 1)t t , we retain χ at χ = χ d , above the critical point, so that phase separation proceeds in the bulk.Throughout the simulations presented in this article, we fix the χ parameters to χ d = 2.7 and χ m = 1.3.In the demixing periods with χ d , the coexistence concentrations of A-and B-rich phases are φ ∼ = 0.893 and 0.107, respectively.The correlation length, also referred to as the interface thickness, is ξ d = C/(χ − χ c ) ∼ = 1.20d.However, during the mixing periods with χ m , the correlation length is given by ξ m = C/{2(χ c − χ)} ∼ = 0.845d.

A. Domain patterns
Figure 2 shows snapshots of the domain patterns in the case of t m = t d = 100 t 0 .The snapshots were obtained at t = 200 t 0 and t = 1000 t 0 , which correspond to the times when the first and fifth demixing periods have just finished, respectively.The average concentrations, φ , of component A are (a) 0.3, (b) 0.5, and (c) 0.7.In the demixing periods with χ d , the volume fractions of the A-rich phase are (a) 24.6%,(b) 50% and (c) 75.4%, respectively.The light blue surfaces represent the isosurfaces of φ = 0.5 and the dark blue sphere represents the Janus particle.At t = 0, the particle is directed toward the z-axis, i.e., n(t = 0) = (0, 0, 1).In each case, we observed that the particle follows the background flows, which are caused by the interface tensions during the phase separations.The mechanisms of the particle motions are discussed below.

B. Particle trajectories
Figure 3 shows typical trajectories of the Janus particle from t = 0 to t = 10 4 t 0 in the binary mixtures of different concentrations.The time intervals are t m = t d = 500 t 0 .The trajectories show that the particle tends to move directionally toward the particle head.Because the particle surface has a heterogeneous affinity to the components, the phase separation proceeds asymmetrically around the particle.This asymmetry of the phase separation dynamics may cause the directional motion of the Janus particle.Interestingly, the particle moves toward its head in all the mixtures.Figure 3 also indicates that the trajectories are not completely straight.The degree of the directionality and the particle speed depend on parameters such as the average concentrations and the time intervals.
Figures 4(a) and (b) show the temporal changes in the trajectory length and velocity toward the head.Since the particle orientation, n, changes with time, the trajectory length d and the velocity V toward the particle orientation are calculated as Positive and negative values of V represnt the forward and backward motions of the particle position, respectively.In Fig. 4, the time intervals are set to t m = t d = 10 3 t 0 .Here, each curve was obtained from one simulation run.We simulated nine average concentrations from φ = 0.1 to 0.9.In the mixtures of φ ≤ 0.1 and φ ≥ 0.8, we did not observe any drastic motion of the particle; hence, their curves have not been included in In Fig. 4(a), the trajectory lengths indicate stepwise motions.The particle is almost fixed in the mixing periods.However, during the demixing periods, the particle shows forward displacements.By repeating these cyclic motions, the particle continuously propels in a periodically phase-separating binary mixture.The displacement in each cycle is of the order of the particle diameter.Although the onsets of the motion in the demixing periods are not clearly seen in Fig. 4, the detailed analyses indicate that the particle does not start moving simultaneously with the quenching into the demixing states.It moves most largely after a certain incubation time t i , which is discussed later.
The cyclic behaviors are also clearly displayed in the particle velocity.After the initial incubation time in each demixing period, the particle velocity shows large positive values.After this transient deterministic motion, the velocity decreases gradually with some fluctuations.In particular, it can have negative values in more symmetric mixtures, with φ = 0.5.A similar stepwise motion is observed in a system where a Janus particle with metallic surfaces is trapped at a liquid-air interface [20].The particle's stepwise motion is due to spontaneous cyclic bursts of bubbles.In our system, the stepwise motion is due to the controlled changes of the interaction parameter.
Figures 5(a) and (b) show the temporal changes of the average-concentration differences and the velocity intensities.They are calculated as where Ω t is the system volume.In the early stages of the demixing periods, the phase separation starts via diffusion of the components.As shown in Fig. 5, the hydrodynamic flow develops simultaneously with the phase separation, and thus, it is small during the early stages [46,47].Because the particle is transported by the background flow, the incubation time of the particle motion t i in Fig. 4 corresponds to the duration of the early stage of the phase separation, t e .Here, t e depends on the average concentration.This is because the growth rate of the concentration field depends on the average concentration.
Inside the spinodal regime, the phase separation proceeds via spinodal decomposition.As the average concentration approaches the spinodal points ( φ = 0.5 ± 0.255 for χ d = 2.7), the growth rate is decreased to zero.However, in the binodal regime, the phase separation occurs via the nucleation of the droplets of the minority phase.
The nucleation rate is also decreased to zero as φ approaches the equilibrium concentration.In both the processes, more symmetric mixtures are more unstable and the durations of the early stages are shortened.This may suggest that the particle moves faster in the symmetric mixtures because the incubation time, t i , during which the particle is at rest, is reduced.We define the duration of the early stage t e as ∆φ 2 (mt t + t m + t e ) = ∆φ 2 eq /2.Here, ∆φ eq is the concentration difference in the equilibrium demixing state.From the simulation results, we obtain t e = 34.6 t 0 for φ = 0.5, t e ≈ 74.0 t 0 for φ = 0.3 and 0.7, and t e ≈ 292 t 0 for φ = 0.2.These are indicated by the arrows in Fig. 6(a).
Figs. 6(a) and (b) show plots of the averaged speeds toward the particle head, V , and perpendicular to it, V ⊥ , as a function of the time interval.Here, we set t m = t d for simplicity.The parallel and perpendicular velocity in the m-th cycle are defined as From these, the averaged velocities V and V ⊥ are obtained with averaging M = 10 cycles as where X refers to X = and ⊥.The error bars in Figs.6(a) and (b) represent the standard deviations of V m and V ⊥m .The ratio of V ⊥ to V is plotted in Fig. 6(c).In Fig. 6(a), each curve of V is non-monotonic with maxima.The maxima peaks suggest that we can choose efficient time intervals for propelling the particle.As discussed above, the hydrodynamic flow has not developed yet in the early stage of the demixing periods.In the cases of small time intervals, the phase-separating times are too short for the hydrodynamic flows to develop sufficiently; hence, the particle is not dragged, implying that we have to maintain the system in the phase-separated state till the conclusion of the early stage t e .However, for large time intervals, the average speed along the orientation becomes low.In the mixtures of φ = 0.2, 0.3, and 0.7, the parallel speed becomes largest at approximately t d ∼ = t e .Thus, larger time intervals are not needed for propelling the particle with high speeds.However, in the symmetric mixture ( φ = 0.5), the parallel speed becomes largest around t d ≈ 200 t 0 , which is larger than the duration of an early stage t e .This is attributed to large background flows emerging from other domains near the particle, as discussed below.
As has been shown in Fig. 3, we observe the particle fluctuations, which are characterized by the perpendicular motion V ⊥ in Fig. 6(b).The particle changes its orientation and the resulting direction of the particle motion.Thus, similar to the V ⊥ changes, the temporal changes of the orientation are also considered as a measure of the particle fluctuation.Fig. 7(a) demonstrates the autocorrelation of the orientation vector.It is calculated as where we set t max = 5 × 10 3 t 0 .In Fig. 7(a), we plot the autocorrelation for φ = 0.3 as a typical example.They decrease with time indicating that the memory of the orientation is gradually lost.Figure 7(a) shows that the decay rate is increased with increasing t d .The autocorrelations for different φ behave in the same way.In Fig. 7(b), we plot the change rate of the orientation ( ṅ2 ) 1/2 , which is calculated for M = 10 as As shown in Figs.6(b) and 7(b), both V ⊥ and ( ṅ2 ) 1/2 are large for large values of t d , in contrast to those for small t d .In the demixing periods, the phase separation proceeds throughout the bulk.The fluctuations of the particle motions and orientation stem from the hydrodynamic flows accompanied by the spontaneous growth of the domains surrounding the particle.These background flows are independent of the particle; thus, they disturb the particle motion and change its orientation.Therefore, if we use large time intervals, the particle motion is likely to deviate from the straight line along the initial orientation n(t = 0).As suggested in Figs. 6 and  7, the ratios V ⊥ / V and ( ṅ 2 ) 1/2 increased slightly with t d .Also, in this sense, the large time interval is not preferred for the controlled propulsions.
Figure 6(a) shows that the propelling speed is higher in the symmetric mixtures ( φ = 0.5) than in the asymmetric mixtures.However, the fluctuations of the particle motion are also large as shown in Figs.6(b) and 7(b).Therefore, the symmetric mixtures are not suitable for keeping the straight line motions.The particle moves largely; however, its motion easily loses the directionality with time.Furthermore, in the asymmetric mixtures, the fluctuations of the particle motion and the orientation are relatively small for the preferred t d ; thus, the asymmetric mixtures are more suitable to control the particle motions.We consider the mechanisms of the propelled motions in the two types of asymmetric mixtures separately.

C. The wettable component-rich mixture
Figure 8 shows the patterns of the evolutions of the concentration and flow field in the mixture of φ = 0.7.As shown in Fig. 3, the particle moves rather straightforward.In the early stage of the phase separation during the demixing periods, the A-rich phase wets the halfportion of the particle and a wetting layer is formed on it.Because of the directional fluxes of the A-component toward the surface, the A-component is depleted near the outside of this first wetting layer.This accumulated layer structure is similar to the oscillating profiles of the concentration field near a flat wall or a particle with a homogeneous surface [47,48].To compensate for the depleted region, non-spherical droplets of the B-rich phase are formed as shown in Fig. 8(b).The droplets grow near the surface via coagulation and coalescence.Because the processes of the coagulation and coalescence occur asymmetrically near the particle tail, the associated hydrodynamic flow likely pushes the particle toward the particle head.This process is shown in Figs.8(c) and  (d).During the phase separation, many small droplets of the B-rich phase are also formed in the bulk.They grow to their typical sizes with time via coalescence and coagulation, or evaporation and condensation processes [27,29].Around coalescing droplets, other hydrodynamic flows are induced, which disturb the particle motion.The disturbing flows are relatively weaker; thus, the particle motion remains straightforward, in contrast to that in the symmetric mixtures.
The mixture of φ = 0.8 should be phase-separated at equilibrium, when χ = χ d (= 2.7).However, because the nucleation rate is very small, the time intervals we employed (t d ≤ 10 3 t 0 ) were not enough to induce the phase separation.Thus, we did not observe any motion of the Janus particle in the mixture of φ ≥ 0.8.

D. The wettable component-deficient mixture
The motions in the B-rich mixtures are highly directed toward the particle orientation as shown in Figs. 3 and  6. Figure 9 shows the pattern evolutions of the concentration and flow field in the mixture of φ = 0.3.In the early stages of the phase separation, a cap-shaped domain of the A-rich phase emerges from the wettable portion of the particle surface.It covers a large amount of the interfacial area, although its volume is not so large.Then, it tends to change its shape to a sphere for reducing the interfacial energy after the early stage of the phase separation.This process induces a pumping hydrodynamic flow around the tail of the Janus particle and the resulting pumping flow pushes the particle toward the head as shown in Figs.9(c) and (d).The decay time of the localized pumping flow, t h , is estimated as t h ≈ ηa/σ, where σ is the interface tension.
The mixtures of φ = 0.3 and 0.7 have the same stability for the phase separation in the bulk.Because the tail of the particle prefers component A, this asymmetry leads to the difference in the particle motions between the φ = 0.3 and the φ = 0.7 cases.Figure 6(a) shows that the highest parallel speed in in the φ = 0.7 mixture is approximately twice of that in the φ = 0.3 case.However, Figure 6(c) indicates that the ratio V ⊥ / V in the φ = 0.7 case is also approximately twice of that in the φ = 0.3 case at the maximum parallel speed.Also, Figure 6(a) shows that the standard deviations of V ⊥m for the φ = 0.3 mixture are smaller than those for the φ = 0.7 case.The particle is suggested to move more steadily and smoothly in the φ = 0.3 mixture.Thus, we conclude that the B-rich mixtures are more preferred to induce more straight motions.
The stability for the bulk phase-separation in the mixture of φ = 0.2 is the same as that in the φ = 0.8 case.However, we have not observed any motions in the φ = 0.8 mixture, although the particle moves straightforward in the φ = 0.2 case.This difference suggests that the domain formation in the mixture of φ = 0.2 is attributed to the heterogeneous nucleation at the particle surface [49].As the average concentration approaches the binodal line, the thermal nucleation rate in the bulk is strongly decreased [29].However, the rate of heterogeneous nucleation on the wetting surface is large enough to induce it during our demixing periods.Because the number of the surrounding droplets is decreased, particle motion becomes more straightforward in more asymmetrically B-rich mixtures (see the case of φ = 0.2 in Fig. 6).
In both cases, the hydrodynamic flow around the particle tail pushes the particle toward the head.In this sense, our particle motion may be categorized to a pusher in the active matter field [11].However, more detailed analyses on the flow pattern are required before this conclusion can be derived.
Figures 4(a) and (b) indicate that the particle does not show any large changes of the position and orientation during the mixing periods.This is because diffusion dominates the mixing but it does not contribute to the particle motion.However, here, we should note that the mixing periods are very important for resetting the binary mixtures.As noted above, the hydrodynamic flow caused by the interface tension is large in the demixing periods, whereas it is small in the mixing periods.This difference in the hydrodynamic flows present between the mixing and demixing periods induces the continuous propelled motion of the Janus particle.
In the above simulations, we set t m = t d for simplicity.However, considering whether this mixing interval is sufficient to reset the binary mixtures is important.The characteristic length of the phase-separated domain increases with time algebraically.In the symmetric mixture, a bicontinuous pattern is formed and the domain grows obeying [28] (t) ≈ cσt/η.
However, in the asymmetric mixtures, the minority phase forms droplets.The droplets grow with time via coalescence and coagulation as [27,29] (t) ≈ c (T t/η) 1/3 .
Here, c and c are non-dimensional numbers, which depend on the volume fractions.Then, the characteristic length at the ends of the demixing periods would be given by (t d ).If the volume fraction is quite small, the droplets grow via nucleation and growth mechanism, and the above scaling relation (Eq.( 27)) is replaced by (t) ∝ (L 0 σt) 1/3 .Because its growth exponent is the same as that of Eq. ( 27), we consider only cases described by Eq. ( 27) above.
In the mixing periods, these phase-separated domains should be dissolved into a homogeneous state via diffusion.The diffusion time is estimated as t dif ∼ (t d ) 2 /D m , where D m = 2L 0 (χ sp −χ m )/T is the diffusion constant at χ = χ m and χ sp is the interaction parameter at the spinodal point.When the time interval of the mixing period is longer than the diffusion time, the system can be reset for the next demixing period.However, if t dif t m , the mixtures show non-steady states [37][38][39][40] and the particle would not move directionally.The condition t dif t m is rewritten as t m (T /η)

IV. SUMMARY AND REMARKS
We proposed a possible mechanism of propelled motions of Janus particles in periodically phase-separating binary mixtures using numerical simulations.Because the particle has an asymmetric surface in wettability, phase separation proceeds heterogeneously around the particle.The resulting asymmetric hydrodynamic flow leads to directional motions of the particle.After a usual one-step quench, the directional motion slows down and will stop eventually.By changing the interaction parameter, the phase-separated system recovers to a one-phase mixing state.Therefore, under periodic changes of the interaction parameter, we can induce a continuous motion along the particle orientation.
We found that the particle propels more directionally in asymmetric binary mixtures.In symmetric mixtures, the hydrodynamic flow emerging from other domains surrounding the particle is so large that it disturbs the directional motions of the Janus particle.We also found that the propelling speed and the directionality depend on the frequency of the change in the interaction parameter.In the cases of short time intervals, the hydrodynamic flow does not develop well and it cannot drive the particle.However, in the cases of long time intervals, the hydrodynamic flows from the surrounding domains disturb the directional motion.We can efficiently move the particle in the intermediate time intervals, which are comparable to the characteristic time of the spinodal decomposition or nucleations in the early stage of the phase separation.The pumping hydrodynamic flow that is localized around the particle decays with the relaxation time t h .t h is of the order of ηa/σ.The interval of the demixing period should be larger than t i + t h .However, the large time interval for the demixing periods leads to the loss of the directionality as discussed above.Thus, the most efficient time interval would be t d ≈ t i + t h .The displacement of the particle in each cycle is of the order of the particle diameter.Then, the maximum particle speed could be increased to V ∼ a/(t i + t h ).
Unfortunately, the simulations performed in this study are limited owing to the numerical costs.We need to deepen our understanding of the propelled motions of such a Janus particle, and we hope that we will present a report on them in the near future.We make some critical remarks to improve our study as follows.
(1) In this article, we show only the simulations with the particle diameter 2a = 12d.Here, d is comparable to the correlation length; hence, our particle is rather small.In actual phase-separating mixtures, such small particles would show drastic Brownian motions and the directional motions we discovered might be smeared out.Our preliminary simulations with larger particles indicated that the displacement of the particle in each cycle is of the order of the particle size.In other words, they suggest that the particle speed can be increased linearly with its size by employing appropriate temporal changes of the interaction parameter.
(2) The wettability of our particle changes smoothly on the surface (see Eq. ( 6)).However, an actual Janus particle usually has two distinct surfaces and the wettability changes abruptly at the equator.Instead of Eq. ( 6), we performed some simulations with an alternate surface function given by Here, d h is introduced to avoid the singularity at the equator of the Janus particles.A particle described by Eq. ( 28) with small d h would behave more realistically.The preliminary simulations using Eq. ( 28) with d h = 0.05 demonstrated essentially similar results as those using Eq. ( 6).This similarity can be derived from the fact that our particle motions are induced by the secondary effect of the phase separation, i.e., the asymmetric growth of the hydrodynamic flow around the particle.Thus, we consider that our findings are robust for the heterogeneous surface pattern.To improve the efficiency of the particle motions, simulations with a variety of the surface structures would be interesting.We consider that the propelled motions in this article are robust also for the details of the free energy function.We numerically confirmed that the particle moves in the same way in binary mixtures described by the Ginzburg-Landau free energy instead of Eq. ( 4).
(3) The easiest method to realize our findings in actual systems would be a direct observation of them with an optical microscope equipped with a temperature control hot stage [40].However, changing the temperature with high frequencies might be experimentally difficult because the thermal diffusion constant is finite.To in-duce a high speed propulsion, large differences of the temperature from the binodal point are preferred for both the mixing and demixing periods.However, for deep quenches in the demixing periods, the incubation times for the phase separation would be reduced up to microscopic timescales.
The pressure control [50] is considered to be an alternative method to induce the continuous propulsions.The combination of the temperature control and illuminationinduced phase separation would also be able to induce periodic phase separation with high speeds [23,24].In this study, we employ only the square waves of the interaction parameter to induce periodic phase separation.We should perform more simulations with other types of wave functions to find more efficient propulsion schemes.
(4) In our model, particle motion is caused by the hydrodynamic flow.However, studying the behaviors of a Janus particle in solid mixtures where v = 0 would also be interesting.In such mixtures, the particle motion is caused by the thermodynamic forces (∝ ∇δF/δφ).Because the coarsening behavior of the phase separation domain pattern depends on the fluidity of the mixtures (see Eqs.( 26) and ( 27)), whether our findings are applicable to the solid mixture is not trivial.
(5) Because the intensity of the hydrodynamic flow is proportional to 1/η 0 (see Eq. ( 7)), we expected that the particle speed is also proportional to 1/η 0 .However, our simulations with different solvent viscosities indicated that V is approximately proportional to 1/η 0 in symmetric mixtures, whereas it is approximately independent of η 0 in asymmetric mixtures (not shown here).This is because the rate-limiting process of the phase separation in asymmetric mixtures is the nucleation and growth of the minority phase; their rates are dominated by the diffusion constant, not the solvent viscosity .Thus, the results reported in this article are not quantitatively universal.They will depend on the system parameters such as the solvent viscosity and interface tension.Although we consider that the propelled motion in periodic phase separation can be qualitatively realized in actual systems, we have to perform more simulations with wide ranges of the system parameters.
For simplicity, we also assumed a constant viscosity η 0 in the surrounding fluids.However, the solvent viscosity is generally a function of the local composition.As noted above, the particle speed depends on the solvent viscosity in symmetric mixtures; thus, studying the influences of the viscosity difference would also be interesting.
where we neglect the contribution of the system boundary.For the force f and torque g fields, we obtained the following equations in our scheme.
Then, we finally obtain where we assumed ∇ • v = 0. Thus, we confirmed that G decreases with time only via diffusion and viscous dissipation.

Fig. 4 .
Fig.4.Also, we did not plot the curves for the mixtures of φ = 0.4 and φ = 0.6 because they essentially demonstrate behaviors similar to those with the symmetric mixture φ = 0.5.

FIG. 6 .
FIG. 6. Plots of the averaged velocity (a) toward the particle head V and (b) perpendicular to the orientation V ⊥ , with respect to the time intervals.The arrows in (a) indicate the durations of the early stage te.The average concentrations are φ = 0.2, 0.3, 0.5 and 0.7.The error bars represent the standard deviation of the particle motions.(c) The ratio of the perpendicular speed V ⊥ to the parallel speed V .

FIG. 7 .
FIG. 7. (a) The autocorrelation function of the orientation N (t).The average concentration was constant at φ = 0.3, whereas the time intervals changed.(b) The dependence of the change rate of the particle orientation ( ṅ2 ) 1/2 on the time interval.

FIG. 9 .
FIG. 9. Snapshots of typical pattern evolutions around the particle in the B-rich mixtures ( φ = 0.3).The arrows show the flow field.The red and blue domains are the A-and B-rich phases, respectively.The time intervals are t d = tm = 300 t0.
/D m for the bicontinuous patterns.If we set t d (T /η) 2 D −3 m , we have t dif t d in the asymmetric mixtures.Therefore, long annealing times for the mixing are not required.The total time interval t t = t m + t d can be reduced to t t ∼ = t d .