Active Brownian particle in homogeneous media of diﬀerent viscosities: numerical simulations

Self-propelled colloids, active polymers and membranes, driven (vibrated) granular layers and hybrid synthetic-biological systems are striking examples of systems containing synthetic active Brownian particles. Such particles autonomously convert the available energy of the environment into their own directed mechanical motion. In most studies the self-propelled Brownian particles move in overdamped media. Recently, experiments with Janus particles in a low-pressure plasma have appeared. A distinctive feature of such a medium is an extremely low viscosity at which the inertial eﬀects play a significant role, resulting in underdamped Brownian motion. At present, there is a lack of statistical theory describing the underdamped Brownian motion of self-propelled particles at all time scales. This paper presents the numerical simulation results of active Brownian motion in homogeneous media of diﬀerent viscosities. The calculations are performed using a mathematical model of a self-propelled Brownian sphere with translational and rotational inertia. The time-dependent mean square displacement and mean linear displacement (the noise-averaged trajectory) of the particle are investigated as a function of medium viscosity, self-propulsion velocity and moment of inertia. Our simulation reveals that the dynamics of a self-propelled spherical particle significantly depends on two independent dimensionless parameters of the particle: the ratio of the self-propulsion velocity to the characteristic thermal velocity and the ratio of the friction coeﬃcient to the rotational diﬀusion coeﬃcient. The obtained statistical characteristics of active Brownian motion are compared with the known theoretical models in a wide range of medium viscosities. We propose simple corrections to the basic theory of overdamped active Brownian motion, which allow one to calculate the eﬀective diﬀusion coeﬃcient and the persistence length of a self-propelled Brownian particle in a medium with any dynamic viscosity. The results obtained are discussed in relation to active particles in a colloidal plasma and superfluid of self-propelled particles. This paper presents the numerical simulation results of active Brownian motion in media of diﬀerent viscosities. The calculations were performed using a mathematical model of a self-propelled Brownian sphere with translational and rotational inertia. The obtained statistical characteristics of active Brownian motion are compared with the known theoretical models in a wide range of medium viscosities. Simple semiempirical corrections are proposed to the basic theory of overdamped active Brownian motion, which allow one to calculate the eﬀective diﬀusion coeﬃcient and the persistence length of a self-propelled Brownian particle in media with any dynamic viscosities.


Introduction
Active matter is a state of condensed matter, which represents an ensemble of active Brownian nano-or microparticles in a dissipative medium. Such particles autonomously convert the available (chemical, electromagnetic, thermal, etc.) energy of the environment into their own directed mechanical motion. [1][2][3][4][5][6][7][8][9][10][11] Many examples of active matter are of biological origin; however, a significant part of the current experimental and theoretical research is devoted to artificial systems: self-propelled colloids, active polymers and membranes, driven (vibrated) granular layers and hybrid synthetic-biological systems. Synthetic active matter is of great fundamental interest in statistical mechanics, opens up new areas in materials science, and has many promising applications: from purely technical to biomedical. 2,3 A striking example of a synthetic active particle is the Janus particle. 12 It is a bifunctional particle the two parts of which differ in chemical composition and as a result have different surface and/or volume properties. The interaction of the Janus particle with the environment generates a driving force that causes the particle to move itself.
Depending on the parameters of the medium and the particle (and the associated characteristic relaxation times), the active Brownian motion can be underdamped and overdamped. 13 In the latter case, inertia does not noticeably effect on the motion of the particle. Most experiments deal with overdamped motion of Janus particles. 2,3 Their selfpropulsion can be caused, for example, by diffusiophoresis, 4,5 thermophoresis, 6,7 electrophoresis (electroosmosis) [8][9][10] or bubble thrust, 11 depending on the specific mechanism of the ''particlemedium'' interaction. However, more recently, experimental studies of Janus particles in a low-pressure plasma have appeared. 14,15 A distinctive feature of such a medium is an extremely low viscosity at which the inertial effects play a significant role, resulting in underdamped Brownian motion. The driving force of Janus particles in these experiments is the photophoretic force. 14, [16][17][18] In most theoretical research, the active Brownian particles are considered to be moving in an overdamped medium. 1 However, there is a lack of theory describing underdamped Brownian motion of self-propelled particles. This paper presents the numerical simulation results of active Brownian motion in media of different viscosities. The calculations were performed using a mathematical model of a self-propelled Brownian sphere with translational and rotational inertia. The obtained statistical characteristics of active Brownian motion are compared with the known theoretical models in a wide range of medium viscosities. Simple semiempirical corrections are proposed to the basic theory of overdamped active Brownian motion, which allow one to calculate the effective diffusion coefficient and the persistence length of a self-propelled Brownian particle in media with any dynamic viscosities.

Basic equations
To demonstrate basic features of the problem, we restrict our analysis here to two spatial dimensions. In this case, a single particle trajectory is described by the particle position vector r(t) = [x, y] and the particle orientation vector n(t) = [cos y, sin y], where y(t) is the angle between the x-axis and n, see Fig. 1. Let us assume that due to some driving force acting along n, the particle moves in a medium with a velocity V S = V S n. The direction of n is determined by the random rotation of the particle, the characteristic frequency of which o is equal to the rotational diffusion coefficient. Neglecting the drag effect of the medium, for purely thermal noise and a spherically symmetric particle of radius R in bulk situations we have where is the translational short-time diffusion coefficient. Here T is the temperature of the medium (expressed in energy units), n is the friction coefficient of the particle, and M is the particle mass.
In the absence of deterministic (other internal and external) forces and neglecting the drag effects of the medium, the twodimensional translational motion of a self-propelled spherical Brownian particle is described by the following Langevin equations: where x x (t) and x y (t) are independent delta-correlated Gaussian processes with zero means and autocorrelation functions representing white noise, i.e., hx The random rotation of the particle caused by fluctuations of the tangential component of the momentum transmitted from the environmental particles is described using the following stochastic differential equation: 19,20 € where m = T/Io is the rotational damping rate of the particle with the moment of inertia I; and x y (t) is a delta-correlated Gaussian process with hx y i = 0 and hx y (0)x y (t)i = 2Tmd(t)/I. In the classical description of Brownian motion, the temperatures corresponding to the translational (3) and rotational (4) noise are assumed to be equal. We stick to this reasonable assumption here. However, note that in some special cases, effective temperatures of the translational and rotational motion may differ, for example: for hot Brownian motion of a heated colloidal particle in a solvent, 21 or for a dust grain absorbing the ambient gas. 22

Overdamped Brownian motion of a self-propelled particle
In most theoretical research, active Brownian particles are considered to be moving in an overdamped medium. In this case, it is usually assumed that the left-hand sides of eqn (3) and (4) responsible for inertia are equal to zero, i.e., ẍ(t) = 0, ÿ(t) = 0 and € y(t) = 0. The inertialess Langevin equations provide a quantitative description of the overdamped particle dynamics at the long-time scale (t c n À1 ). In particular, the mean square displacement (MSD) of an overdamped self-propelled Brownian particle in an arbitrary direction is described by the following equation: 4,23 Further, we assume x(0) 0. Note that eqn (5) is valid for times significantly longer than the momentum relaxation time t p = n À1 . Eqn (5) leads to the following effective (long-time) diffusion coefficient: x axis is directed along the vector of the particle's self-propulsion velocity at the moment t = 0. The x* axis is directed along the vector of the initial total velocity of the particle.
The first term of eqn (6) is the translational short-time diffusion coefficient (2) associated with passive Brownian motion. The second term represents the contribution from the driving force due to the coupling between translational and rotational motions: it is a well-known solution of the so-called telegraph equation, used to describe a persistent random walk of living cells 24 in the absence of velocity fluctuations, i.e., for | : r(t)| V S . Coupling of rotational diffusion with directed active motion is quantified by the persistence length where the inverse rotational diffusion coefficient o À1 is a characteristic time scale (the persistence time) during which the particle is more likely to move in a similar direction to its initial orientation. If we direct the x axis along the vector of the particle's self-propulsion velocity at t = 0, as shown in Fig. 1, then the mean linear displacement (MLD) of an active colloidal particle is described by the following equation: 3,13,23 This noise-averaged function is also called the 'average trajectory' 3 or 'mean position' 23 of the particle.

Underdamped Brownian motion of a self-propelled particle with inertia
For a correct description of the overdamped Brownian motion at very short time scales (t { t p ), the translational inertia must be taken into account. In this limit the well-known Ornstein-Uhlenbeck formula 25 describes the MSD of a passive particle (with V S 0): Mn 2 e Ànt À 1 ð Þ : Note that, in contrast to eqn (5), formula (9) is valid at all times of observation. However, there is currently no analytical equation describing the time-dependent MSD for a self-propelled Brownian particle with inertia. In order to correctly describe the underdamped active Brownian motion (even at long-time scales), the rotational inertia must be taken into account. 13 Solving the system of eqn (3) and (4) leads to the following generalized form of the effective (long-time) diffusion coefficient: 13,26,27 where F is the dimensionless function that depends on the ratio Io 2 /T and corrects eqn (6). The asymptotic dependences of F(Io 2 /T) can be derived analytically for two limits: 26,27 for Io 2 /T { 1, and for Io 2 /T c 1.
In contrast to the diffusion, the effect of inertia on the mean linear displacement and the persistence length of a selfpropelled Brownian particle has not been studied at all.

Simulation parameters
Hereinafter, for the moment of inertia of a spherical particle relative to the axis of rotation passing through its center of mass the following parameter is used: a o 1, which depends on the mass distribution over the volume of the particle (here we assume that the distribution is symmetric with respect to the axis of rotation). We also introduce two independent dimensionless parameters responsible for the particle dynamics: the inertial delay number 13 which is the ratio of the friction coefficient to the rotational diffusion coefficient and indicates the inertial delay between the self-propulsion direction and the particle velocity, originating from the finite relaxation times in the system; and the dimensionless self-propulsion velocity which is the ratio of the self-propulsion velocity to the characteristic thermal velocity. Note that in colloids, the Péclet number is commonly used as a measure of selfpropulsion: 3,28,29 The dynamics of an active Brownian particle was simulated by numerically integrating the equations of motion (3) and (4). Discrete random changes in the velocities of translational and rotational motion were set by where Dt is an integration step and z i is a Gaussian random number with zero mean and unit variance. The integration step Dt was less than 0.05 min[m À1 ,o À1 ]. The calculations were performed for a particle with diameter 2R = 1 mm and mass M = 7.9 Â 10 À13 g (corresponding to the averaged density r = 1.5 g cm À3 ) in a viscous medium with a temperature T = 0.025 eV. The mass distribution parameter a was varied from 0.01 to 1. The friction coefficient n was varied in the range from 3.9 to 3.9 Â 10 6 s À1 , which, for example, corresponds to a change in pressure from B10 À1 to B 10 5 Pa in the argon atmosphere. The Péclet number was in the range from 0.1 to 1000 and the ratio of characteristic frequencies b was in the range from 10 À6 to 10 6 . At the initial moment, each particle was located at the origin x(0) = 0, y(0) = 0, and the angle y(0) represented a uniformly distributed random variable in the range from 0 to 2p. Fig. 2 shows exemplary trajectories of simulated particles for a fixed Péclet number and different inertial delay numbers. In order to accumulate good statistics, we got up to 10 6 independent particle trajectories for each set of simulation parameters.

Mean-square displacement.
To illustrate the numerically calculated MSD of a self-propelled Brownian sphere, Fig. 3 represents the mass transfer functions hx 2 (t)i/2t for fixed a and various sets of b and g (or Pe). Note that in comparison with MSD, the mass transfer function has several advantages: 30,31 (I) more pronounced transitions between ballistic and diffusion regimes; (II) equality to the effective diffusion coefficient at long observation times. Along with the simulation results, Fig. 3 shows the mass transfer functions for eqn (5) and (9).
With a relatively high friction (b c 1, see Fig. 3a) at times t { o À1 , an active Brownian particle moves in the same way as a passive one: motion of both particles is ballistic, i.e., hx 2 (t)i p t 2 , at the shortest time scale t { n À1 ; a crossover to a short-time diffusive regime with hx 2 (t)i p t occurs at t 1 B n À1 . However, at longer time scales, the behavior of active and passive particles differs: there is a second crossover from a short-time diffusive to a propulsive ballistic regime at t 2 B g À2 n À1 , and, finally, diffusive motion occurs after the third crossover at t 3 B o À1 . Note that similar crossovers have been obtained earlier in theoretical studies of a persistent random walk. 23,32 Some crossover corresponding to point t 3 was also observed in the experiment with the Janus colloid particle. 4 For b c 1 and for all values of Pe (used in the simulations), the numerically calculated MSD of the self-propelled Brownian sphere at all observation times is well described by the following equation: which is a combination of eqn (5) and (9). Fig. 3b and c show that when the friction coefficient is comparable or less than the rotational diffusion coefficient, i.e., b t 1, eqn (17) is inapplicable for describing the MSD of an active Brownian particle. In this case, a random walk of an active Brownian particle corresponds to a passive one with some enhanced diffusion coefficient, i.e., D 4 D T . Fig. 4 shows the dimensionless correction F for the effective (long-time) diffusion coefficient (10) in dependence on the ratio b/a. For numerical data, function F was calculated as F = 2Pe À2 (D/D T À 1), where D = hx 2 (t)i/2t was measured at times t \ 100max[n À1 ,o À1 ]. For comparison, the plots of  eqn (11) and (12) are also presented in Fig. 4 taking into account that Io 2 /T 3a/4b. It is easy to see that eqn (11) is consistent with the simulation results at b/a c 1, and eqn (12) is in satisfactory (but not good) agreement with our simulations at b/a { 1. To approximate the numerical data for all values of a and b, we propose the following expression for the dimensionless correction of the effective diffusion coefficient (see symbols in Fig. 4): The difference between eqn (18) and the simulation results shown in Fig. 4 does not exceed 1.5% which is comparable with the statistical measurement error in the numerical experiment. Substituting eqn (18) into eqn (10) in the limit b/a -N, we obtain eqn (6). For b/a { 1, eqn (18) provides a better agreement with the simulation results than eqn (12). In the latter case (in contrast to eqn (6)), the contribution from the driving force (self-propulsion) to the effective (long-time) diffusion coefficient is independent on the rotational diffusion coefficient: Up to this point, we have calculated the MSD for the x-projection. Since the medium is considered homogeneous, then hx 2 i = hy 2 i, and for the particle position vector, we simply get hr 2 (t)i = 2hx 2 (t)i = 4Dt.
3.2.2 Persistent random walk in the direction of the initial driving force. Fig. 5 represents the mean linear displacement (MLD) of a self-propelled Brownian sphere along the direction of its initial orientation (see x in Fig. 1) which was calculated numerically for fixed a and various sets of (b, g), and normalized to L S V S /o. In addition to the simulation results, Fig. 5 also shows the plot of eqn (8) derived for an overdamped active Brownian particle. 13,23 This comparison reveals that eqn (8) is valid only for b/a c 1. The persistence length L of a selfpropelled Brownian particle was calculated as an average value of MLD(t) at times t \ 100max[n À1 ,o À1 ]. The normalized persistence length L/L S is shown in Fig. 6 in dependence on the inertial delay number b for various sets of a and g. The numerically obtained values of L/L S turned out to be independent of the dimensionless self-propulsion velocity g. To fit the numerical data for all values of a and b, we propose the following dimensionless correction for the persistence length of a self-propelled Brownian sphere: where The difference between eqn (21) and the simulation results shown in Fig. 6 is within 5% at b B 1 and does not exceed 2% for large and small values of b. In the limit b -N, we get L E L S . And it turns out that for small b, the persistence length   of a self-propelled Brownian sphere is independent on the rotational diffusion coefficient: 3.2.3 Persistent random walk in the direction of the initial movement. In the above section, as in previous studies, 3,13,23 the mean linear displacement of an active particle was determined along its initial self-propulsion velocity, i.e., along the orientation vector n(0). However, in real experiments with selfpropelled spherical particles, it could be difficult to distinguish between the thermal and self-propulsion components of the total particle velocity. In this case, the mean linear displacement could be determined along the direction of the initial total velocity of the particle, i.e., along : r(0), as shown in Fig. 1. Here, the mean linear displacement of the simulated particle in the direction of its initial motion (MLD*) was calculated as follows: Fig. 7 represents the numerically calculated MLD* of the selfpropelled Brownian sphere for fixed a and various sets of b and g. In addition to the simulation results, Fig. 7 also shows the plots of eqn (8) and theoretical MLD* for a passive Brownian particle, which is written as: where L T = hV T in À1 is the effective persistence length of a passive particle with an initial velocity equal to the arithmetic mean velocity of its Brownian motion in two-dimensional Fig. 7 reveals that eqn (8) describes the MLD* only for b c 1 and g c 1.
At long observation times (t c max[n À1 ,o À1 ]) the MLD* function tends to the effective persistence length L* for all values of b and g. The following simple expression for L* can be used to approximate the numerical data for g t 1 and all values of b, as well as for b c 1 and all values of g: The numerically calculated L* normalized to {L T + L S /(1 + 1/g)} is plotted in Fig. 8 as a function of b and clearly shows the validity range of eqn (25). For numerical data, L* was obtained as an average value of MLD*(t) at times t \ 100max[n À1 ,o À1 ]. The effective persistence length of an active particle with low activity (g { 1) tends to the persistence length of a passive particle (i.e., L* E L T ) only at relatively low friction (b { 1), and in the opposite case (b c 1) it depends on the Péclet number: For strong activity (g c 1) and high friction (b c 1), we get L* E L S . However, we have not been able to find a simple and physically sound approximation for the effective persistence length of an active particle for b { 1 and   (8), and the thin dotted line plots eqn (24) for the passive particle (g = 0). Here the time is normalized to the momentum relaxation time t p = n À1 . g c 1. In this case, additional information on the correlation of the particle velocities caused by two different random processes (x x and x y ) is required. This is the subject of a separate analytical study and is beyond the scope of this work. Fig. 3 and 5 show that the MSD and MLD functions obtained for the same Péclet number and different inertial delay numbers differ markedly. Therefore, the Péclet number cannot be a universal comparative characteristic (measure) of selfpropulsion. Our simulation reveals that the MSD and MLD functions of a self-propelled spherical particle with a given mass distribution over its volume (a) significantly depend on two independent dimensionless parameters: the inertial delay number (b) and the dimensionless self-propulsion velocity (g). Moreover, the ratio of long-time to short-time diffusion coefficients is completely determined by (a, b, g): from eqn (2), (10) and (18)

Discussion
For both underdamped and overdamped classical Brownian motion, the inertialess (overdamped) Langevin equation provides a quantitative description of the dynamics of a passive Brownian particle at long time scales (t c n À1 ) regardless of the damping rate in the medium. As for active Brownian motion, comparison of the numerically calculated MSD and MLD functions with the existing analytical solutions (5) and (8) clearly shows that the overdamped Langevin equations fail to describe the long-time behavior of an active Brownian particle when its friction coefficient is not so high compared to the rotational diffusion coefficient, i.e., b t 1. For b c 1 the translational inertia of particle's active Brownian motion is important only at the shortest time scale (t t n À1 ), while the rotational inertia can be completely neglected. In the opposite case (b t 1), the rotational inertia of an active particle plays an essential role at all time scales.
For overdamped active Brownian motion (b c 1) with a small dimensionless self-propulsion velocity (g { 1), both the MSD and MLD functions exhibit four consecutive regimes of motion separated by three prominent crossovers (see Fig. 3a and 7a). It should be noted here that in experiments with colloidal particles, 2-4,28 the first two crossovers at short and intermediate time scales are beyond the temporal and spatial resolutions of typical observation equipment. As g increases, the intermediate crossover merges with the others, and for g c 1, only one crossover is found at t B o À1 . In this limit, the numerically calculated MSD agrees with a persistent random walk of a living cell and a non-living active particle in the absence of velocity fluctuations. 24,33 Based on our simulations, we can conclude that a selfpropelled particle with a large self-propulsion velocity (g c 1) and a large inertial delay number (b c 1) behaves like a passive Brownian particle in a medium with significantly less effective friction n eff B n/b and high temperature T eff E g 2 T/2.
As the particle's rotational diffusion increases, motion becomes more random and less directed. For underdamped motion with b/a t 1 and any self-propulsion velocity, we obtain that both the MSD and MLD functions have single crossovers at t 1 B n À1 (see Fig. 3b and 7b). Thus, in the limit b { 1, our simulation does not support the theoretical prediction 32 that a single crossover of the MSD function should be at point t 1 B o À1 . Our finding is that a random walk of an active Brownian particle in the underdamped limit corresponds to a passive particle in a medium with the same viscosity but with a higher temperature T eff % T Â 1 þ 0:25g 2 ffiffiffiffiffiffiffi ffi 5ab p À Á . In the classical description, the translational diffusion coefficient of a spherical Brownian particle is written as where Z is the dynamic viscosity of a medium. Given eqn (1), (2) and (26), it is easy to get For a typical micron-sized colloidal particle in a liquid medium with a dynamic viscosity of 10 À3 C10 À1 Pa s at room temperature, we obtain b B 10 7 C10 11 . In experiments with laboratory complex plasmas, argon at pressures P B 1C10 Pa is usually used as a buffer gas. 34 In this case, the coefficient of dynamic viscosity can be calculated using the Epstein formula: 35 where d = 0.5 + p/8 for specular reflection of neutrals from the particle surface and d = 1 + p/8 for complete accommodation, 36,37 and u n is the average thermal velocity of neutrals. Taking into account eqn (28), we find that at low gas pressures, b is proportional to P 2 and R 3 . Thus, in highly rarefied media, the inertial delay number can vary over a wide range. For example, b B 0.1 for a particle with R = 10 mm at a pressure of 1 Pa, and b B 10 for R = 10 mm and P = 10 Pa. Here, we assume that T is equal to room temperature. However, it is worth noting that the kinetic temperature of a single Brownian particle levitating in a gas discharge at room temperature can be higher (usually up to two times). [38][39][40] With such combinations of experimental parameters, the existing statistical theories developed for the underdamped and overdamped limits of active Brownian motion are inapplicable, see Fig. 4. These simple estimates show that in experiments with active Brownian particles in rarefied gas media, including weakly ionized plasmas, both overdamped and underdamped motions are possible. The transition from overdamped to underdamped motion is easily accomplished by changing the buffer gas pressure. In addition to Janus particles in a gas-discharge plasma, 14,15 suitable conditions can be realized for self-propelled particle pairs 41,42 formed by the effective nonreciprocal forces of the wake-mediated interaction between microparticles. 43 Another unusual example of underdamped active Brownian motion is the so-called vibrobots, which convert vibrational energy into directed motion using their tilted elastic legs. 27,44 A non-obvious question to ask is what kind of motion would one expect for a self-propelled particle in 'nonviscous' fluids. 3 A promising example would be the so-called superfluid microswimmer: a Janus particle suspended in superfluid helium self-propels due to the transformation of the superfluid to a normal fluid on the heated side of the particle. 45 This situation can arise both in bulk 46 and under a free surface of superfluid helium. 47 At a temperature of 1.6 K the viscosity of liquid helium is 1.3 Â 10 À6 Pa s. 48 Of course, for a realistic description of particle motion under such conditions, the 'particle-quantum vortex' interaction must be taken into account. 49 However, using our simple (non-quantum) formalism, we can estimate that the momentum relaxation time of such a particle with R = 10 mm can be t p = (2/9)R 2 r/Z B 20 ms and the inertial delay number b can be B10 5 . Although this time scale is available for a typical visualization system, instead of both short-time ballistic and diffusion regimes, an enhanced propulsive (ballistic) regime will be observed due to the huge self-propulsion velocity of the Janus particle (g c 1). 45 This directional movement, in turn, can be interrupted by periodic releases of quantum vortices and be limited by the finite dimensions of the experimental setup.
The two-dimensional model for active Brownian motion considered here can be generalized to three spatial dimensions. Taking into account the third dimension in the individual dynamics of a particle, does not bring any new physical effects but can change (by a factor of t 3) the values of the calculated constants. 29 In any case, this may become a topic for the next study. And in contrast to the individual dynamics, the collective dynamics can differ significantly for two and three dimensions. 1,3,29,34

Conclusions
In summary, at present, there is a lack of statistical theory describing underdamped Brownian motion of a self-propelled particle at all time scales. We carried out the numerical simulation of active Brownian motion in media of different viscosities.
For this purpose, we used a mathematical model a self-propelled Brownian sphere with translational and rotational inertia. It was shown that in contrast to underdamped passive Brownian motion, the inertialess Langevin equations fail to describe the diffusive regime of underdamped active Brownian motion. The time-dependent mean square displacement and mean linear displacement (the noise-averaged trajectory) of the particle were investigated as a function of medium viscosity, self-propulsion velocity and moment of inertia. It was shown that the Péclet number, which is often used in active colloids, cannot be applied as a comparative characteristic (measure) of self-propulsion of active Brownian particles in media with very different viscosities. Our simulation revealed that the dynamics of a self-propelled Brownian spherical particle with a given mass distribution over its volume significantly depends on two independent dimensionless parameters of the particle: the ratio of the self-propulsion velocity to the characteristic thermal velocity and the ratio of the friction coefficient to the rotational diffusion coefficient.
We propose simple corrections to the basic theory of overdamped active Brownian motion, which allow one to calculate the effective diffusion coefficient and the persistence length of a self-propelled Brownian particle in media with any dynamic viscosities. The obtained corrections perfectly describe the results of numerical simulation over the entire viscosity range. Moreover, in the underdamped limit, our diffusion correction provides a better agreement with the simulation than the correction proposed earlier in ref. 13, 26 and 27. We also proposed to calculate the mean linear displacement of an active particle not along its initial self-propulsion velocity, but along the initial orientation of the total particle velocity (including its thermal component). This approach is justified, because in real experiments with self-propelled Brownian particles of a spherical shape, it could be difficult to distinguish between the thermal and self-propulsion velocity components. Based on the results of numerical simulation, a simple semiempirical correction was also proposed, which makes it possible to calculate the effective persistence length of a selfpropelled Brownian particle in this reference frame.
Our results can be useful for describing motion of selfpropelled particles on a vibrating plate or suspended in gaseous media and nonviscous fluids. In particular, laboratory experiments with self-propelled particles in a weakly ionized gas discharge plasma have great potential for studying active Brownian motion since they enable to control the damping rate of the medium over a wide range.

Conflicts of interest
There are no conflicts to declare.