PaperlayerView exact revision on arXiv

Source-linked AI summary

Autophoresis of a Janus particle near a planar wall: a lubrication limit

Tachin Ruangkriengsin, Günther Turk, Howard A. Stone

arXiv:2603.00791v2cond-mat.softphysics.flu-dyn

TL;DR

Near-wall Janus-particle dynamics are difficult to resolve numerically because confinement produces steep solute gradients and complex interactions. This paper uses lubrication-limit asymptotics for axisymmetric and slightly tilted configurations, finding that cap size controls whether tilt is restored or amplified near contact.

  • Problem

    Numerical studies have difficulty resolving Janus-particle flow and transport near walls, where confinement and steep solute gradients complicate the dynamics.

  • Method

    The paper develops an asymptotic lubrication analysis for a Janus particle near an impermeable wall, treating mixed boundary conditions and slightly tilted orientations.

  • Results

    The angular velocity changes sign at Φ ≈4.60, producing restoring rotation for 𝜙/𝜖1/2 < 4.60 and destabilizing rotation for 𝜙/𝜖1/2 > 4.60.

  • Takeaways & Limitations

    Cap size relative to the lubrication region determines near-contact rotational stability, while even a small active or inert region can strongly influence wall-induced motion.

  • Takeaways & Limitations

    The analysis is restricted to small tilt angles, Ψ ≪1, and therefore cannot capture finite-angle steady skating states.

Abstract

from arXiv · show

We study the self-diffusiophoresis of a spherical chemically active particle near a planar, impermeable wall, with a focus on the influence of particle orientation on propulsion. We analyze a Janus particle with asymmetric surface chemical activity, consisting of a small inert region within a catalytically active cap. While numerical simulations have been used to study such particles, they encounter difficulties resolving the flow and transport in the near-wall regime due to geometric confinement and steep solute concentration gradients. We address this limitation through an asymptotic analysis in the lubrication limit, where the gap between the particle and the wall is narrow. In particular, we consider the distinguished limit in which the inert region is asymptotically comparable in size to the lubrication region. We analyze an axisymmetric configuration in which the inert face is oriented parallel to the wall and extend the analysis to slightly tilted orientations. We find that the cap size determines whether a tilted particle rotates back toward the axisymmetric state or continues to reorient, thereby characterizing its rotational stability in the near-contact regime.

I. INTRODUCTION

The paper studies self-diffusiophoresis of chemically patterned Janus particles near walls, where confinement alters motion and existing numerical methods struggle at very small gaps. It uses lubrication-limit asymptotics to analyze axisymmetric and slightly tilted configurations.

  • Background: Autophoresis arises when surface reactions generate solute-concentration gradients that induce diffusiophoretic slip and self-propulsion.Active regions impose fixed solute flux, whereas inert regions impose no flux.
  • Motivation: Near-wall confinement can produce skating, hovering, or reflection, while chemically patterned boundaries can guide and sort Janus particles.
  • Existing approaches: Existing boundary-element, multipole, bispherical-coordinate, and Galerkin methods face limitations in resolving steep gradients at extreme near-wall separations.Some methods require separations above 1.1 particle radii or artificial repulsive potentials, while bispherical expansions are restricted to simple geometries.
  • Contribution: The paper applies asymptotic analysis to a Janus particle in the near-wall regime, where transport is governed by the narrow lubrication gap.For 1–5-micron colloids, the regime corresponds to separations of a few hundred nanometers or less.
  • Problem setup: The analysis begins with an inert face parallel to the wall and extends to slightly tilted orientations in the near-contact regime.The gap distance is 𝜖a, with 𝜖≪1.
  • Problem setup: The model considers a large active cap with a small inert face, while also allowing the complementary small-active-cap configuration.The inert-face angular extent is assumed small, 𝜙≪1.

A. Dimensionless formulation

The formulation nondimensionalizes solute transport, phoretic slip, and Stokes flow, then introduces stretched lubrication coordinates in which the active and inert regions meet at R=Φ. The leading concentration is determined by regional solutions joined through transition matching.

  • A. Dimensionless formulation: Lengths are scaled by particle radius a, while concentration, velocity, pressure, and force use scales set by α, D, β, and μ.The characteristic scales are aα/D, αβ/D, μαβ/(aD), and μaαβ/D, respectively.
  • A. Dimensionless formulation: The particle boundary is rewritten in dimensionless variables before the analysis turns to lubrication-scale coordinates.
  • A. Dimensionless formulation: For the prescribed-flux model, both solute transport and flow are dominated by the narrow gap, unlike kinetic models requiring matched particle-scale and gap-scale regions.
  • B. Scalings in the lubrication limit: In the lubrication region, radial distance scales as 𝜖^1/2, gap height as 𝜖, and pressure as 𝜖^-2.
  • B. Scalings in the lubrication limit: The axisymmetric force acts only in the wall-normal direction, and the gap pressure determines the corresponding vertical force.
  • B. Scalings in the lubrication limit: The active–inert boundary lies at R=Φ, separating the active region R>Φ from the inert region R<Φ.The two regions are solved separately and matched near R=Φ.
  • A. Solute concentration: In the axisymmetric gap, the concentration and flow are independent of θ, and the solute flux conditions distinguish the active and inert regions.
  • A. Solute concentration: The concentration is expanded separately as C+ and C−, with leading-order fields independent of Z and regional regularity conditions imposed at the origin and far field.

B. Flow analysis

The lubrication analysis expands velocity and pressure in the small gap, using mass conservation and pressure matching to determine the leading-order flow and force. In the small-inert-face limit, the vertical hydrodynamic force approaches the fully active-particle result.

  • Velocity and pressure are expanded in powers of the gap parameter 𝜖.
  • The leading-order Stokes equations constitute the lubrication approximation, with pressure independent of the axial coordinate 𝑍.
  • Radial mass conservation in the fixed axisymmetric gap relates the pressure gradient to the solute concentration gradient.
  • Matching the inner gap pressure to the outer pressure imposes decay at large radial distance and permits explicit integration of the leading-order pressure.
  • In the limit Φ ≪1, the leading-order vertical hydrodynamic force approaches 3𝜋, consistent with a fully active particle near a wall.The gap dynamics are dominated by the active face when the inert portion is small relative to the lubrication region.

C. A complementary configuration

The complementary configuration places a small active cap toward the wall and applies the same lubrication framework. When the active region is large relative to the gap, the vertical force again approaches 3𝜋, while hovering lies outside the asymptotic regime considered.

  • The complementary configuration has a small active cap and large inert face, with the active cap oriented toward the wall.
  • The leading-order excess concentrations are defined separately inside and outside the active-cap boundary, with matching imposed through continuity of the radial concentration gradient.
  • The leading-order vertical hydrodynamic force is obtained by integrating the pressure over the gap-scale region.
  • In the limit Φ ≫1, the leading-order vertical hydrodynamic force approaches 3𝜋 because the active portion is large relative to the lubrication region.Even a small active cap can dominate near-wall gap dynamics when its size exceeds the lubrication scale.
  • The plotted vertical velocity is normalized by 𝛼𝛽/𝐷 and shown as a function of Φ for both original and complementary configurations.
  • Hovering lies outside the asymptotic regime 𝜙≪1, 𝜖≪1, with Φ = 𝜙/𝜖1/2 = 𝑂(1), and does not occur in the force-free setup.The text associates near-field hovering with inert faces that are not asymptotically small or with larger separations.

IV. SLIGHTLY TILTED JANUS PARTICLE NEAR A PLANAR WALL

The analysis extends the parallel inert-face configuration to a particle tilted by a small angle about the 𝑦-axis. Tilt is scaled like the inert-face size, and the general gap problem is treated perturbatively when the scaled tilt is small.

  • The inert surface is tilted about the 𝑦-axis by a small angle 𝜓≪1 relative to the wall.
  • When 𝜓 = 0, the tilted-configuration expression reduces to the axisymmetric result.
  • Because the active–inert boundary combines 𝜙 and 𝜓, the tilt is assigned the same scaling used to compare 𝜙 with 𝜖1/2.
  • The leading-order gap boundary is R = Φ + Ψ cos 𝜃, but analytical progress requires the perturbative limit Ψ ≪Φ.This domain perturbation expands the boundary about R = Φ.

A. Solute concentration

For the tilted particle, the leading-order gap concentration retains angular dependence and is determined through a perturbed active–inert boundary. A local transition analysis supplies the flux matching needed to solve the perturbative concentration problem.

  • Unlike the axisymmetric formulation, the tilted problem retains angular dependence in both flow and solute concentration.
  • The active–inert boundary in the gap is shifted to R = Φ + Ψ cos 𝜃, separating the C+ and C− concentration domains.
  • The leading-order concentration is independent of the axial coordinate 𝑍 and follows from the O(𝜖) problem and the particle-surface flux condition.
  • The governing equations alone do not determine angular variations, so a local analysis near the active–inert boundary is required.
  • For Ψ ≪1, domain perturbation expands the concentration in Ψ, with the leading term equal to the axisymmetric solution.
  • The perturbed concentration uses angular ansätze that convert boundary conditions into radial problems solved separately inside and outside R = Φ.The interior solution remains regular at the origin, while the exterior solution is regular in the far field.
  • The resulting interior and exterior solutions have hypergeometric-function forms, with constants fixed by the perturbed boundary conditions.

B. Flow analysis

The flow analysis extends the lubrication framework to first order in tilt, coupling solute-driven slip to velocity and pressure fields in the narrow gap. It derives a Reynolds pressure equation and determines pressure matching across the inert–active transition.

  • B. Flow analysis: The velocity and pressure fields are expanded at leading order in ε and first order in the tilt parameter Ψ.This is denoted the (01)-order in the preceding perturbation framework.
  • B. Flow analysis: The tilted velocity and pressure fields use angular dependences matched to the solute concentration, reducing the lubrication Stokes equations to radial–axial differential equations.The assumed forms are U_θ(R,Z) sin θ for velocity and P(R,Z) cos θ for pressure.
  • B. Flow analysis: Slip boundary conditions at the particle surface and no-slip conditions at the wall determine the velocity profile, while gap-integrated mass conservation yields a differential equation for pressure.The pressure is independent of Z in the tilted correction.
  • B. Flow analysis: The Reynolds equation’s right-hand side vanishes away from the inert–active edge R = Φ, so pressure is solved separately inside and outside that radius.The inner and outer pressures P− and P+ satisfy regularity at the origin and decay in the far field, respectively.
  • B. Flow analysis: Pressure remains continuous at R = Φ, whereas its radial gradient has a jump set by the discontinuity in the solute concentration gradient.Weak integration across R = Φ supplies the jump condition and determines the pressure constants Q− and Q+.

C. Phoretic motion

The phoretic-motion analysis converts the tilted near-wall flow and pressure fields into force-free and torque-free particle velocities. It finds positive translation in x for all Φ, while angular velocity reverses at Φ ≈ 4.60, changing rotational stability.

  • C. Phoretic motion: Force and torque in the x- and y-directions vanish by symmetry for the axisymmetric configuration but arise at leading order from tilt.The force and torque are normalized by μaαβ/(Dε) and μa^2αβ/(Dε), respectively.
  • C. Phoretic motion: The particle velocities V_x and Ω_y are obtained by combining numerically integrated force and torque with mobility coefficients under force-free and torque-free conditions.The calculation solves for concentration and pressure constants before evaluating velocity, force, torque, and mobility relations.
  • C. Phoretic motion: For all Φ, the slightly tilted particle translates in the positive x- and z-directions, while its rotation direction depends on Φ.Figure 5 summarizes the resulting motion directions.
  • C. Phoretic motion: Φ ≈ 4.60 marks the angular-velocity sign change: Φ < 4.60 restores the axisymmetric orientation, whereas Φ > 4.60 destabilizes it.For fixed inert-face angle ϕ, the ratio Φ = ϕ/ε^1/2 increases as the particle approaches the wall.
  • C. Phoretic motion: The predicted motion directions qualitatively agree with a bispherical-coordinate analysis, although detailed comparison is limited to the slightly tilted regime.The related analysis also reports a reversal of rotational direction with separation.
  • C. Phoretic motion: In the complementary configuration with a small active cap facing the wall, V_x is negative and rotation changes from counterclockwise at moderate closeness to clockwise at extreme near-wall separation.The vertical velocity V_z retains the original direction but has a different magnitude.

V. CONCLUSION

The lubrication analysis identifies how cap size and particle tilt govern near-wall translation and rotational stability, while also clarifying the model’s scope and physical limitations.

  • Axisymmetric configuration: The parameter Φ = 𝜙/𝜖1/2 compares cap size with gap distance and captures how strongly the wall effectively sees the active or inert region.Even a small inert face or active cap can substantially influence motion at sufficiently small gaps.
  • Framework and extensions: The mixed-boundary lubrication framework extends to patterned activity and establishes continuity of solute concentration and its first derivative through catalytic edges.The authors identify spatially varying mobility or activity and substrate-driven flows as analytically tractable extensions.
  • Slightly tilted configuration: For slight tilts, the wall-parallel translational velocity remains positive for all Φ, while rotational velocity changes sign at a finite Φ.Thus, the axisymmetric state can be rotationally stable or unstable depending on particle–wall separation.
  • Limitations: The analysis is limited to small tilt angles, so it describes local stability near an axisymmetric state rather than arbitrary orientations or finite-angle skating.A full Ψ = O(1) solution would be needed to determine whether steady skating occurs in the three-dimensional lubrication limit.
  • Limitations: Brownian motion is neglected; rotational diffusion may still perturb tilted particles and reduce near-wall retention times, despite suppressed wall-normal fluctuations at small gaps.This caveat is especially relevant for experimentally sized particles and near-wall rotational behavior.
  • Limitations: The zeroth-order kinetic model predicts mechanical power scaling as O(𝜖−1), suggesting that finite-rate reaction kinetics may be required in sufficiently narrow gaps.The divergence arises even though particle velocity remains O(1), because slip-driven and motion-induced lubrication forces each scale as O(𝜖−1).

Appendix A: Solute concentration: transition-region analysis

The transition-region analysis matches inner and gap-scale solute solutions across the active–inert boundary, determining leading-order regularity and continuity conditions.

  • Transition-region setup: Stretched coordinates resolve an O(𝜖) transition region centered at the active-cap–inert-face boundary, where radial and vertical scales are comparable.This scaling leads to a leading-order Laplace problem in the transition region.
  • Transition-region setup: The concentration is expanded in powers of 𝜖1/2, and matching conditions connect the transition-region solution to the gap-scale fields on both sides of R = Φ.The expansion introduces C(0), C(1), and higher-order terms used to enforce matching.
  • Leading-order solution: At leading order, C(0) satisfies Laplace’s equation with no-flux conditions at the particle and wall, and the bounded solution is constant.The far-field constants therefore coincide with C(0).
  • Matching and regularity: The transition analysis establishes continuity of C(0) and dC(0)/dR across R = Φ, while d2C(0)/dR2 need not be continuous there.This regularity result follows from matching rather than directly patching the outer solutions.

Appendix B: Solute concentration: a connection formula across a general catalytic edge

The general catalytic-edge analysis extends the transition-region matching procedure to angularly dependent boundaries and derives continuity conditions for the leading-order concentration and its normal derivative.

  • General edge geometry: The inert-face region Ω is bounded by a smooth catalytic edge Γ, with the axisymmetric and tilted cases represented by R < Φ and R < Φ + Ψ cos 𝜃.This formulation allows the edge geometry to vary in the lateral plane.
  • Matching result: The local transition analysis yields continuity of the leading-order solute concentration and its normal derivative across Γ.These conditions provide the closure needed for the angularly dependent tilted problem.
  • Local transition region: Near Γ, arclength and signed-distance coordinates define stretched variables for resolving the local transition region.The leading-order gap height is approximated by its value on the edge.
  • Local transition region: The concentration is expanded in powers of 𝜖1/2, and the leading-order transition problem satisfies Laplace’s equation over a strip with homogeneous Neumann conditions.The strip spans the stretched normal coordinate and the local gap height.
  • Matching result: The leading-order problem admits only a constant bounded solution, yielding continuity of the leading-order solute concentration across Γ.Matching to the gap solutions supplies the far-field conditions for the transition problem.
Loading 2603.00791v2…