Preprint
Article

Biased Random Process of Randomly Moving Particles with Fixed Mean and Group Velocities

Altmetrics

Downloads

279

Views

187

Comments

1

Tao Guo  *

This version is not peer-reviewed

Submitted:

14 October 2023

Posted:

16 October 2023

You are already at the latest version

Alerts
Abstract
In a particle swarm with a fixed kinetic energy, the constituent particles move randomly and their velocity magnitudes follow a Maxwell distribution characterized by specific parameters. Within a certain time frame, particles in a sub-swarm may exhibit a bias in their movement direction. Thus, the particles may exhibit a special characteristic (biased random motion), where the group velocity in a particular direction remains constant u. For this biased particle swarm, composed of randomly moving particles with a fixed average speed c, particles have a higher probability of moving in a specific direction. Conversely, their likelihood of moving in all other directions is equally reduced. This study starts from the perspective of biased random walks and using Mathematica software proves that the diffusion rate of particles (in the reference frame Ru observed from R0) in all directions is slower than that of an unbiased particle swarm with the same average speed c. Specifically, the degree of reduction is determined by the Lorentz-like factor. Lastly, we present the Ito equation for this biased random motion and provide associated verification examples. This study aims to serve as a reference for comprehending the underlying principles of the special relativity effect.
Keywords: 
Subject: Physical Sciences  -   Mathematical Physics

Introduction

Most random phenomena in nature exhibit characteristics of bias 1 . This includes random motion, as such phenomena tend not to occur in a strictly neutral or uniform manner. In a group composed of randomly moving particles, the motions of particles in its sub-groups tend to be biased. Among these biased random motions, several distinct motion processes warrant particular attention. These special motions include particles with a higher probability of movement in a specific direction, while maintaining an equivalent yet lower probability of movement in other directions 2 ; particles with greater probabilities of movement toward a point, while maintaining an equivalent yet lower probability of movement in other directions 3 ; and biased random processes of particles capable of generating rotation 4 . The study of such phenomena typically begins with the simplest biased random processes.
Biased random processes, akin to the first scenario previously described, have been the subject of extensive research[5]. The nature of the bias can shape the form of the expression for a random process. For example:
d x ( t ) = u S t d t + σ S t d w ( t ) ,
where w ( t ) denotes a Brownian motion process. In Eq. 1, when the drift (u) and diffusion ( σ ) terms grow synchronously according to S t , geometric Brownian motion will be produced. There will be fractional Brownian motion when particles follow the motion law in Eq. 2.
d x ( t ) = θ [ u x ( t ) ] d t + σ d w ( t ) ,
where θ is the regression speed. Additionally, there is a simpler scenario where particles with the same jump speed have a higher probability of moving in a specific direction and equal probabilities of moving in other directions. This is a common biased phenomenon in physics. In a system with constant kinetic energy, it is generally challenging to alter the distribution of particle speeds, while the directions of motion can change more readily. When describing such systems using Eq. 3
d x ( t ) = u d t + σ d w ( t ) ,
the system’s kinetic energy will increase as the drift term u increases. It is obvious that such an equation does not reflect real scenarios. In contrast, when the drift term u increases, it may be a better solution to reduce the diffusion term σ accordingly. Unfortunately, there is currently no clear quantitative relationship available. This article will delve into this issue in detail and derive a quantitative relationship of Lorentz-like factor between the drift (u) and diffusion ( σ ) terms. This kind of problem is similar to the problem that the particle speed or time slows down in a moving frame in special relativity. Some researchers 6 - - 9 demonstrated the diffusion process (random walk) with a special relativistic effect. However, it is the diffusion or random walk of particles when their speed is affected by the special relativity effect. The case involved in this study is that the particle velocity is determined, and the diffusion in the biased random particle swarm formed by them is restricted by the Lorentz-like factor. They are not the same kind of problem.
In my previous research 2 , I demonstrated the special relativistic-like effect of biased random motion when particles move at a uniform speed of c. This article extends that work by considering particles whose speeds follow a Maxwell distribution with equal average speed c. I establish that such biased random motion retains its special relativistic-like effect. Additionally, I present the corresponding Itô equations and simulate the motion of related examples. My findings offer insights into the intricate relationship between biased random processes and the special relativity effect.

Results and Discussion

Biased Random Walk

In order to elucidate the relationship between the drift and diffusion terms, this study commences with an examination of random walks. Among various biased random walk scenarios, we specifically focus on a situation where the probability of a particle moving in one of the six possible directions in a 3-dimensional space is denoted as p, while the probability in the remaining five directions is given by 1 p 5 , subject to the condition 1 6 < p < 1 . Although the case when 0 < p < 1 6 presents an intriguing scenario, it falls outside the purview of this study. Given a particle’s jump speed as c, the collective or group velocity of all particles can be described by the equation:
u = 6 p 1 5 c .
This can be expressed equivalently in the form X 1 ( t ) , X 2 ( t ) , X 3 ( t ) , t T , where X 1 ( t ) , X 2 ( t ) and X 3 ( t ) represent 1-dimensional random walks with a jump speed of c 3 . We introduce a 3-dimensional rectangular coordinate system, aligning its (1, 1, 1) direction parallel to the group velocity u . Consequently, the 1-dimensional random walk model along the x i -axis, which is one of the three equivalent coordinate axes, is given by:
P x , bia = 1 + u / c 2 , n > 0 , 1 u / c 2 , n < 0 , 0 , otherwise ,
where n N . The corresponding unbiased 1-dimensional random walk model on the x i -axis among the three equivalent coordinate axes is
P x , unbia = 1 2 , n > 0 , 1 2 , n < 0 , 0 , otherwise ,
where n N . The relative standard deviation of the two 1-dimensional random walk processes along the two equivalent x i -axes is
σ bia σ unbia = 12 25 ( p 1 ) ( 3 p + 2 ) .
Consider the 3-dimensional random walk vectors, representing velocities, formed by X 1 ( t ) , X 2 ( t ) and X 3 ( t ) over a unit time interval. Their magnitudes follow the discrete Maxwell distribution, and the average magnitude is proportional to the standard deviation of any component vector along the equivalent axes 2 . Consequently, the average velocity magnitude of the biased 3-dimensional random walk is c 2 u 2 c times that of its unbiased counterpart, with u defined in Eq. 4 (see Part 1 of the Supplementary Information for the detailed Mathematica code). Evidently, the biased case represents a decelerative process. This also illuminates the relationship between the drift term and the attenuation rate of diffusion in a continuous process with constant energy from an alternative perspective. For a more rigorous depiction of this relationship in continuous time, further proofs are warranted. Yet, directly proving this for a continuous-time stochastic process poses challenges. We can transpose to the following strategy.

Continuous Biased Random Process

When particles in swarm A possess speeds governed by a Maxwell distribution with a mean of c, the associated scale parameter for this distribution is 1 2 π 2 c . Particles in A can be categorized into subgroups A i , based on their speeds, where each subgroup’s velocity vectors terminate on a series of spheres of radius r i ( i = 1 , 2 , 3 , ). The endpoints of these vectors for particles in A i are uniformly distributed over the sphere with radius r i . The number of particles in A distributed across these spheres follows a Maxwell distribution with respect to r, with the scale parameter 1 2 π 2 c and mean radius r ¯ = c . Hence, the scenario where particle speeds are uniformly c can be viewed as a specialized case ( r i = c ). The collective motion characteristics of particles whose velocity vectors terminate on spheres of radius r i embody the motion traits of particles with speeds dictated by the Maxwell distribution. The particle counts in both A and A i are substantial, with a uniform distribution. When a subgroup A of the parent particle swarm—where the magnitudes of the particle velocities follow a Maxwell distribution—moves along the z-axis (with a reference frame based on the parent particle swarm) at an average speed of u, it is equivalent to each sub-particle group A i ( i = 1 , 2 , 3 , ) moving along the z-axis with an average speed of r i u c .
Let R follow a Maxwell distribution with a scale parameter λ . Then a R will follow a Maxwell distribution with a scale parameter a λ . Specifically, along the z-axis, all the component speeds u i in particle swarm A i , given by u i = r i u c (where i = 1 , 2 , 3 , ), follow a Maxwell distribution with a mean value of u. Moreover, the standard deviations σ i = r i 1 u 2 / c 2 3 (where i = 1 , 2 , 3 , ) of the mixed distributions (we can regard the particles in A i as a particle swarm possessing a mixed distribution with weight w = c + u 2 c ) across all spherical layers also conform to a Maxwell distribution. For the case along the x- or y-axis, the particles in A i can also be viewed as a particle swarm possessing a mixed distribution with weight w = c + u 2 c . Adopting a method analogous to that in my previous work 2 , it can be demonstrated that the ratio of the standard deviation of the mixed distribution along the x- or y-axis to that of the unbiased moving particle on the same axis satisfies the relationship c 2 u 2 c , and the ratio of the standard deviation along the z-axis also satisfies this relationship (see Part 2 of the Supplementary Information for the detailed Mathematica code).

Itô Equation of Biased Random Processes

The above results indicate that for particles possessing a group velocity of u and with speeds following the Maxwell distribution, their speeds within R u exhibit a deceleration characterized by a Lorentz-like factor. Despite this change in speed, their distributions still follow the Maxwell distribution within R u . This behavior permits their motion to be described by the following Itô equation when viewed from R 0 . It should be noted that the definitions of R 0 and R u align with those presented in my previous study 2 .
d x 1 ( t ) = u 1 d t + σ c 2 u 2 c d w 1 ( t ) , d x 2 ( t ) = u 2 d t + σ c 2 u 2 c d w 2 ( t ) , d x 3 ( t ) = u 3 d t + σ c 2 u 2 c d w 3 ( t ) ,
where u 1 , u 2 and u 3 are the drift speeds and u 1 2 + u 2 2 + u 3 2 = u ; c is the average speed of particles; w 1 ( t ) , w 2 ( t ) and w 3 ( t ) are Brownian motion processes; and σ is their standard deviation. Then, the slice distribution of Eq. 8 is
c 3 e c 2 ( x i t u i ) 2 2 σ 2 t ( c 2 u 2 ) 2 2 π 3 2 σ 3 t 3 2 ( c 2 u 2 ) 3 2 ,
where i takes the values 1, 2 and 3, in accordance with the Einstein summation convention. This means that ( x i t u i ) 2 is equivalent to 1 3 ( x i t u i ) 2 . Moreover, the forward equation of Eq. 8 is
p ( x , t ) t = σ 2 ( c 2 u 2 ) 2 c 2 2 p ( x , t ) u i p ( x , t ) x i ,
where p ( x , t ) is the probability of a particle at position x ( x 1 , x 2 , x 3 ) and time t. This is also a drift-diffusion (or advection-diffusion) equation. Only the diffusion coefficient of this equation is limited by the advection term.
This result (3-dimensional biased Brownian motion) can be simulated by three 1-dimensional biased Brownian motions X 1 ( t ) , X 2 ( t ) , X 3 ( t ) , t T . Fig. 1 depicts the trajectories of 300 particles simulated to move randomly over 10 s with a time step of 0.01 s. The randomly moving particles are subject to scenarios with group speeds of 6 (Figure 1 a ) and 0 (Figure 1 b ), respectively. Analyzing the diffusion effects, there appears to be no significant difference when compared to an unbiased scenario unaffected by the relativistic-like effect. Further, when considering diffusion velocity, scenarios with a notable group speed show a reduced rate of diffusion.

Conclusions

For a particle swarm composed of particles with a fixed energy, assumed to be in reference frame R 0 , the magnitudes of the particle velocities follow a Maxwell distribution with parameter 1 2 π 2 c , yielding an average speed of c. When a sub-swarm of particles (assumed to be in reference frame R u ), exhibits a higher probability of motion in a specific direction and equal yet lower probabilities in other directions (with the group velocity of this sub-swarm denoted as u), as observed from R 0 , the velocities of particles within this sub-swarm decelerate, adhering to a pattern determined by the Lorentz-like factor c 2 u 2 c . Consequently, the biased motion of such particles can be described by a set of Itô equations (Eq. 8). Both the time slice distribution and the forward equation for these biased particles differ from those in the reference frame R 0 .

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org Mathematica code for necessary calculation process and graphics. The following file is available free of charge.

Data Availability Statement

All data generated or analysed during this study are included in this published article and its supplementary information files

Acknowledgments

The author is thankful for Wolfram Inc

Conflicts of Interest

There are no conflicts to declare

References

  1. Tapiero, C.S.; Vallois, P. Run length statistics and the Hurst exponent in random and birth-death random walks. Chaos, Solitons & Fractals 1996, 7, 1333–1341. [Google Scholar]
  2. Guo, T. Study on the average speed of particles from a particle swarm derived from a stationary particle swarm. Scientific Reports 2021, 11, 1–4. [Google Scholar] [CrossRef] [PubMed]
  3. Guo, T. Dynamics of stochastic-constrained particles. Scientific Reports 2023, 13, 2759. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, T.; Guo, T. The angular speed distribution of randomly moving-particle group. AIP Advances 2022, 12, 045005. [Google Scholar] [CrossRef]
  5. Codling, E.A.; Plank, M.J.; Benhamou, S. Random walk models in biology. Journal of the Royal society interface 2008, 5, 813–834. [Google Scholar] [CrossRef] [PubMed]
  6. Debbasch, F.; Chevalier, C. Relativistic stochastic processes. AIP Conference Proceedings. American Institute of Physics, 2007, Vol. 913, pp. 42–48. [CrossRef]
  7. Dunkel, J.; Hänggi, P. Relativistic brownian motion. Physics Reports 2009, 471, 1–73. [Google Scholar] [CrossRef]
  8. Hakim, R. Relativistic stochastic processes. Journal of Mathematical Physics 1968, 9, 1805–1818. [Google Scholar] [CrossRef]
  9. Dunkel, J.; Talkner, P.; Hänggi, P. Relativistic diffusion processes and random walk models. Physical Review D 2007, 75, 043001. [Google Scholar] [CrossRef]
Figure 1. Simulation results of the trajectories of 300 particles randomly moving for 10 s ( c = 10 , u = 6 and σ = 2 ). a , group speed of u. b , group speed of 0.
Figure 1. Simulation results of the trajectories of 300 particles randomly moving for 10 s ( c = 10 , u = 6 and σ = 2 ). a , group speed of u. b , group speed of 0.
Preprints 87804 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated