Diffusion velocity modulus of self-propelled spherical and circular particles in the generalized Langevin approach

This paper presents a generalized Langevin framework to model the averaged velocity modulus of self-propelled spherical and disk-shaped Brownian particles in a harmonic potential, revealing that while an internal Ornstein-Uhlenbeck mechanism induces spontaneous velocity fluctuations, these effects diminish over time as the system evolves.

The Gist

Here is an explanation of the paper using simple language and creative analogies.

The Big Picture: A Self-Driving Robot in a Stormy Ocean

Imagine a tiny, self-driving robot (a particle) floating in a giant, warm ocean (a thermal fluid). This ocean isn't empty; it's filled with billions of tiny, invisible waves and currents (thermal noise) that constantly bump into the robot, pushing it around randomly. This is what physicists call Brownian motion.

Usually, if you want this robot to move in a specific direction, you have to push it with an external force (like a fan blowing on it). But in this research, the robot has a built-in engine. It can generate its own speed from the inside. The author, Pedro Colmenares, wanted to figure out exactly how fast this "self-driving" robot moves on average when it's also being pushed by the ocean waves and trapped in a gentle, elastic cage (a harmonic potential).

The Two-Step Dance

The author realized that describing this movement is too complicated to do all at once. So, he split the robot's life into two distinct "dances":

1. The Internal Engine (The "Self-Propulsion"):
   Imagine the robot has a nervous system made of three independent, jittery springs (mathematically called Ornstein-Uhlenbeck processes). These springs randomly push the robot's legs, giving it an initial burst of speed.

   • The Analogy: Think of a person trying to walk while their legs are being randomly kicked by invisible elves. The person's speed isn't constant; it jitters up and down. This is the "inner mechanism" that gives the particle its initial push.

2. The Ocean Drift (The "Diffusion"):
   Once the robot has that initial push, it enters the ocean. Here, it gets hit by the water molecules (thermal bath) and is also pulled back by a giant, invisible rubber band (the external potential) that tries to keep it in one spot.

   • The Analogy: Imagine you are running on a treadmill that is also on a boat in a storm. You are running (self-propulsion), but the boat is rocking (thermal noise) and a rope is pulling you back (the external field).

The paper combines these two steps into a single, complex equation (a Generalized Langevin Equation) to predict the robot's Velocity Modulus—which is just a fancy way of saying "how fast it is going on average, regardless of which way it's pointing."

The Shape Matters: Ball vs. Coin

The author didn't just look at one shape; he looked at two:

• The Sphere (3D Ball): A fully round particle that can roll in any direction (up, down, left, right, forward, backward).
• The Disk (2D Coin): A flat particle that can only slide around on a flat surface (like a coin on a table).

The Discovery:
When the author simulated these shapes, he found something interesting.

• The Jittery Start: At the very beginning, the robot's speed fluctuates wildly. It speeds up, slows down, and changes direction because of its internal "jittery springs."
• The Calm Down: As time goes on, the internal engine settles down. The wild fluctuations fade away, and the robot reaches a steady, predictable average speed.
• The Difference: The 3D ball showed much more complex and "bumpy" speed patterns than the 2D coin. The coin's speed was a smooth, steady climb to a finish line. The ball, having more directions to wiggle in, had a more chaotic journey before settling down.

Why This Matters

This research is like creating a better map for nanotechnology.

• Real-world application: Scientists are building tiny machines (nanomotors) that can swim inside our bodies to deliver medicine. These machines are often self-propelled.
• The Problem: Current models often assume these machines move at a constant speed or ignore how they interact with the "temperature" of their environment.
• The Solution: This paper provides a more accurate "physics engine" for these simulations. It accounts for the fact that these tiny machines start with a random internal kick and then interact with the fluid around them.

The "Takeaway" Metaphor

Think of the particle as a drunk sailor on a swaying ship (the thermal bath) who is also trying to walk a tightrope (the external potential).

• Old models assumed the sailor was either perfectly sober or just randomly stumbling.
• This paper says: "Wait, the sailor has a caffeine pill (the internal engine) that makes him jittery at first, but he eventually finds his rhythm."

The author calculated exactly how fast this sailor moves on average, considering both the caffeine jitter and the swaying ship. The result is a more realistic way to predict how tiny, self-driving particles behave in the real world, which is crucial for designing future medical nanobots and understanding how bacteria move.

In Summary

• The Goal: Predict the speed of a self-driving particle in a fluid.
• The Method: Split the problem into "internal jitter" (engine) and "external drift" (fluid).
• The Result: The particle starts with wild speed fluctuations but eventually settles into a steady pace. The 3D shape (sphere) behaves more chaotically than the 2D shape (disk).
• The Impact: Better tools for designing microscopic machines and understanding biological movement.

Technical Summary

Here is a detailed technical summary of the paper "Diffusion velocity modulus of self-propelled spherical and circular particles in the generalized Langevin approach" by Pedro J. Colmenares.

1. Problem Statement

The paper addresses the theoretical modeling of accelerated self-propelled diffusive particles (ASPDP). While Active Brownian Particles (ABPs) are commonly modeled with constant propulsion velocities, this research focuses on particles that gradually accelerate due to an internal mechanism before reaching a diffusive stationary state.

The specific challenges addressed are:

• Deriving the Diffusion Velocity Modulus (VM) for particles with internal propulsion mechanisms described by stochastic processes, rather than constant velocity.
• Accounting for the particle's interaction with a thermal bath of harmonic oscillators and an external harmonic potential (parabolic field).
• Distinguishing between the internal propulsion dynamics and the external diffusive dynamics, specifically for 3D spherical and 2D circular (disk) geometries.
• Critiquing existing "energy depot" models for Active Matter, suggesting they require re-examination for proper Hamiltonian descriptions.

2. Methodology

The author employs a Generalized Langevin Equation (GLE) framework, splitting the system's dynamics into two distinct stochastic processes:

A. The Two-Process Framework

1. Internal Propulsion (Process 1):

   • The internal self-velocity is generated by a set of independent Ornstein-Uhlenbeck Processes (OUP).
   • This mechanism is assumed to be insensitive to external perturbations.
   • The OUPs provide the initial time-dependent velocity vector, $\mathbf{v}_p(t)$, which acts as the initial condition for the diffusion process.
   • The dynamics are governed by: $dv_j(t) = -\kappa_j v_j(t)dt + \epsilon_j dW_j(t)$, where $\kappa$ is hydrodynamic drag and $\epsilon$ is noise intensity.

2. External Diffusion (Process 2):

   • The particle diffuses in a thermal bath under an external parabolic potential.
   • The dynamics are described by a modified Generalized Langevin Equation (GLE) derived in the author's previous work [25].
   • Crucially, the reservoir Hamiltonian includes a term for the interaction between bath particles and the external field, ensuring thermodynamic consistency.
   • The initial conditions for this GLE are the velocity and position vectors generated by the OUPs (Process 1).

B. Mathematical Derivation

• Velocity Correlation: The total velocity $\mathbf{v}_d(t)$ is expressed as a sum of the modulated initial propulsion velocity and a colored noise term arising from the thermal bath.
• Velocity Modulus (VM): The author calculates the root-mean-square (RMS) diffusion velocity, $s_d(t) = \sqrt{\langle \mathbf{v}_d(t) \cdot \mathbf{v}_d(t) \rangle}$.
• Coordinate Transformations:
  • For the Sphere (3D): The OUP system is transformed into spherical coordinates $(s_p, \theta, \phi)$. The author derives a new set of stochastic differential equations (SDEs) for the velocity modulus and angular phases using Itô calculus and complex variable transformations (adapting Gardiner's polar method).
  • For the Disk (2D): The 3D equations are reduced to polar coordinates, recovering known results for planar systems.

3. Key Contributions

• New Analytical Model: The paper provides the first analytical derivation of the diffusion velocity modulus for a self-propelled particle where the propulsion is an accelerating OUP process rather than a constant velocity.
• 3D Spherical Derivation: It presents a novel derivation of the OUP in spherical coordinates, resulting in a coupled system of SDEs for the velocity modulus and angular phases. This derivation had not appeared in the literature prior to this work.
• Critique of Energy Depot Models: The author argues that standard "energy depot" models for Active Matter (which assume a specific balance of energy conversion) may be flawed when viewed through a strict Hamiltonian perspective, as the functional derivative of the local squared velocity can be multi-valued.
• Thermodynamic Consistency: The use of the modified GLE with a field-interacting bath ensures the system satisfies the fluctuation-dissipation theorem and yields physically consistent thermodynamic properties.

4. Results

• General Velocity Modulus Equation:
  The derived RMS velocity modulus for the ASPDP model is:
  $$s_d(t) = \left[ \chi^2(t) \sum_{j} \frac{\epsilon_j^2}{2\kappa_j}(1 - e^{-2\kappa_j t}) + \frac{3k_B T}{M}(1 - \chi^2(t)) \right]^{1/2}$$
  Where $\chi(t)$ is the system susceptibility to the external field.
• Spontaneous Fluctuations:
  • The system exhibits spontaneous fluctuations in the diffusion velocity modulus at short times due to the internal OUP mechanism.
  • These fluctuations are caused by changes in the direction of the propulsion speed.
  • As time increases ($t \to \infty$), these momentary fluctuations fade out, and the system reaches a constant VM, consistent with passive Brownian motion behavior at long times.
• Geometry Dependence (Sphere vs. Disk):
  • Sphere (3D): Numerical simulations of the derived SDEs show rich, non-monotonic behaviors (e.g., initial bumps, varying saturation times) depending on the anisotropy of drag ($\kappa$) and noise ($\epsilon$) parameters in the $x, y, z$ directions.
  • Disk (2D): The velocity modulus for the disk is a monotonic saturation curve regardless of parameter combinations. The "richness" of shapes seen in the 3D sphere is lost in the 2D case due to reduced parameter space and symmetry.
• Convergence: Simulations indicate that while the velocity modulus converges to a stable average after sufficient realizations (e.g., >2000), the angular phases ($\theta, \phi$) do not converge to a single value, reflecting the stochastic nature of the orientation.

5. Significance

• Theoretical Advancement: This work bridges the gap between passive Brownian motion and active matter by introducing a physically consistent model for accelerating self-propelled particles. It moves beyond the simplification of constant propulsion velocity.
• Applicability: The model is proposed as a framework for describing nano-motors and other micro-swimmers that do not instantly reach a terminal velocity but accelerate through an internal mechanism.
• Methodological Rigor: By explicitly including the interaction between the thermal bath and the external field in the Hamiltonian, the study ensures that the derived dynamics are thermodynamically sound, addressing limitations in previous overdamped Langevin approaches.
• Geometric Insight: The comparison between spherical and circular particles highlights how dimensionality and internal parameter anisotropy significantly alter the transient dynamics of self-propelled systems, a factor often overlooked in simplified 2D models.

In summary, Colmenares provides a rigorous stochastic framework for analyzing the transient and steady-state velocity of self-propelled particles, demonstrating that internal acceleration mechanisms lead to complex, geometry-dependent fluctuations before settling into a diffusive steady state.