2-dimensional random walk of a silveradatom on an Ag(111) surface[1]Simulation of the Brownian motion of a large particle, analogous to a dust particle, that collides with a large set of smaller particles, analogous to molecules of a gas, which move with different velocities in different random directions
Brownian motion is the random motion ofparticles suspended in a medium (aliquid or agas).[2] The traditional mathematical formulation of Brownian motion is that of theWiener process, which is often called Brownian motion, even in mathematical sources.
This motion pattern typically consists ofrandom fluctuations in a particle's position inside a fluid sub-domain, followed by a relocation to another sub-domain. Each relocation is followed by more fluctuations within the new closed volume. This pattern describes a fluid atthermal equilibrium, defined by a giventemperature. Within such a fluid, there exists no preferential direction of flow (as intransport phenomena). More specifically, the fluid's overalllinear andangular momenta remain null over time. Thekinetic energies of the molecular Brownian motions, together with those of molecular rotations and vibrations, sum up to the caloric component of a fluid'sinternal energy (theequipartition theorem).[3]
This motion is named after the Scottish botanistRobert Brown, who first described the phenomenon in 1827, while looking through a microscope atpollen of the plantClarkia pulchella immersed in water. In 1900, the French mathematicianLouis Bachelier modeled the stochastic process now called Brownian motion in his doctoral thesis, The Theory of Speculation (Théorie de la spéculation), prepared under the supervision ofHenri Poincaré. Then, in 1905, theoretical physicistAlbert Einstein publisheda paper in which he modelled the motion of the pollen particles as being moved by individual watermolecules, making one of his first major scientific contributions.[4]
The direction of the force of atomic bombardment is constantly changing, and at different times the particle is hit more on one side than another, leading to the seemingly random nature of the motion. This explanation of Brownian motion served as convincing evidence thatatoms and molecules exist and was further verified experimentally byJean Perrin in 1908. Perrin was awarded theNobel Prize in Physics in 1926 "for his work on the discontinuous structure of matter".[5]
Themany-body interactions that yield the Brownian pattern cannot be solved by a model accounting for every involved molecule. Consequently, only probabilistic models applied tomolecular populations can be employed to describe it.[6] Two such models of thestatistical mechanics, due to Einstein and Smoluchowski, are presented below. Another, pure probabilistic class of models is the class of thestochastic process models. There exist sequences of both simpler and more complicated stochastic processes which converge (in thelimit) to Brownian motion (seerandom walk andDonsker's theorem).[7][8]
Reproduced from the book ofJean Baptiste Perrin,Les Atomes, three tracings of the motion of colloidal particles of radius 0.53 μm, as seen under the microscope, are displayed. Successive positions every 30 seconds are joined by straight line segments (the mesh size is 3.2 μm).[9]
Observe what happens when sunbeams are admitted into a building and shed light on its shadowy places. You will see a multitude of tiny particles mingling in a multitude of ways... their dancing is an actual indication of underlying movements of matter that are hidden from our sight... It originates with the atoms which move of themselves [i.e., spontaneously]. Then those small compound bodies that are least removed from the impetus of the atoms are set in motion by the impact of their invisible blows and in turn cannon against slightly larger bodies. So the movement mounts up from the atoms and gradually emerges to the level of our senses so that those bodies are in motion that we see in sunbeams, moved by blows that remain invisible.
Although the mingling, tumbling motion of dust particles is caused largely by air currents, the glittering, jiggling motion of small dust particles is caused chiefly by trueBrownian dynamics; Lucretius "perfectly describes and explains the Brownian movement by a wrong example".[10]
The discovery of this phenomenon is credited to the botanistRobert Brown in 1827. Brown was studying plant reproduction when he observedpollen grains of the plantClarkia pulchella in water under a microscope. These grains contain minute particles on the order of 1/4000th of an inch (6 microns) in size. He observed these particles executing a jittery motion. By repeating the experiment with particles of inorganic matter he was able to rule out that the motion was life-related, although its origin was yet to be explained.[11][12]
The mathematics of much of stochastic analysis including the mathematics of Brownian motion was introduced byLouis Bachelier in 1900 in his PhD thesis "The theory of speculation", in which he presented an analysis of the stock and option markets. However this work was largely unknown until the 1950s.[13][14]: 33
Albert Einstein (in one of his1905 papers) provided an explanation of Brownian motion in terms of atoms and molecules at a time when their existence was still debated. Einstein proved the relation between the probability distribution of a Brownian particle and thediffusion equation.[14]: 33 These equations describing Brownian motion were subsequently verified by the experimental work ofJean Baptiste Perrin in 1908, leading to his Nobel prize.[15]Norbert Wiener gave the first complete and rigorous mathematical analysis in 1923, leading to the underlying mathematical concept being called aWiener process.[14]
The instantaneous velocity of the Brownian motion can be defined asv = Δx/Δt, whenΔt <<τ, whereτ is the momentum relaxation time. In 2010, the instantaneous velocity of a Brownian particle (a glass microsphere trapped in air withoptical tweezers) was measured successfully. The velocity data verified theMaxwell–Boltzmann velocity distribution, and the equipartition theorem for a Brownian particle.[16]
There are two parts to Einstein's theory: the first part consists in the formulation of a diffusion equation for Brownian particles, in which the diffusion coefficient is related to themean squared displacement of a Brownian particle, while the second part consists in relating the diffusion coefficient to measurable physical quantities.[17] In this way Einstein was able to determine the size of atoms, and how many atoms there are in a mole, or themolecular weight in grams, of a gas.[18] In accordance toAvogadro's law, this volume is the same for all ideal gases, which is 22.414 liters at standard temperature and pressure. The number of atoms contained in this volume is referred to as theAvogadro number, and the determination of this number is tantamount to the knowledge of the mass of an atom, since the latter is obtained by dividing themolar mass of the gas by theAvogadro constant.
The characteristic bell-shaped curves of the diffusion of Brownian particles. The distribution begins as aDirac delta function, indicating that all the particles are located at the origin at timet = 0. Ast increases, the distribution flattens (though remains bell-shaped), and ultimately becomes uniform in the limit that time goes to infinity.
The first part of Einstein's argument was to determine how far a Brownian particle travels in a given time interval.[4] Classical mechanics is unable to determine this distance because of the enormous number of bombardments a Brownian particle will undergo, roughly of the order of 1014 collisions per second.[2]
He regarded the increment of particle positions in time in a one-dimensional (x) space (with the coordinates chosen so that the origin lies at the initial position of the particle) as arandom variable () with someprobability density function (i.e., is the probability density for a jump of magnitude, i.e., the probability density of the particle incrementing its position from to in the time interval). Further, assuming conservation of particle number, he expanded thenumber density (number of particles per unit volume around) at time in aTaylor series,where the second equality is by definition of. Theintegral in the first term is equal to one by the definition of probability, and the second and other even terms (i.e. first and other oddmoments) vanish because of space symmetry. What is left gives rise to the following relation:Where the coefficient after theLaplacian, the second moment of probability of displacement, is interpreted asmass diffusivityD:Then the density of Brownian particlesρ at pointx at timet satisfies thediffusion equation:
Assuming thatN particles start from the origin at the initial timet = 0, the diffusion equation has the solution:[19]This expression (which is anormal distribution with the mean and variance usually called Brownian motion) allowed Einstein to calculate themoments directly. The first moment is seen to vanish, meaning that the Brownian particle is equally likely to move to the left as it is to move to the right. The second moment is, however, non-vanishing, being given byThis equation expresses the mean squared displacement in terms of the time elapsed and the diffusivity. From this expression Einstein argued that the displacement of a Brownian particle is not proportional to the elapsed time, but rather to its square root.[17] His argument is based on a conceptual switch from the "ensemble" of Brownian particles to the "single" Brownian particle: we can speak of the relative number of particles at a single instant just as well as of the time it takes a Brownian particle to reach a given point.[20]
The second part of Einstein's theory relates the diffusion constant to physically measurable quantities, such as the mean squared displacement of a particle in a given time interval. This result enables the experimental determination of the Avogadro number and therefore the size of molecules. Einstein analyzed a dynamic equilibrium being established between opposing forces. The beauty of his argument is that the final result does not depend upon which forces are involved in setting up the dynamic equilibrium.
In his original treatment, Einstein considered anosmotic pressure experiment, but the same conclusion can be reached in other ways.
Consider, for instance, particles suspended in a viscous fluid in a gravitational field. Gravity tends to make the particles settle, whereas diffusion acts to homogenize them, driving them into regions of smaller concentration. Under the action of gravity, a particle acquires a downward speed ofv =μmg, wherem is the mass of the particle,g is the acceleration due to gravity, andμ is the particle'smobility in the fluid.George Stokes had shown that the mobility for a spherical particle with radiusr is, whereη is thedynamic viscosity of the fluid. In a state of dynamic equilibrium, and under the hypothesis of isothermal fluid, the particles are distributed according to thebarometric distributionwhereρ −ρo is the difference in density of particles separated by a height difference, of,kB is theBoltzmann constant (the ratio of theuniversal gas constant,R, to theAvogadro constant,NA), andT is theabsolute temperature.
Perrin examined the equilibrium (barometric distribution) of granules (0.6microns) ofgamboge, a viscous substance, under the microscope. The granules move against gravity to regions of lower concentration. The relative change in density observed in 10 microns of suspension is equivalent to that occurring in 6 km of air.
Dynamic equilibrium is established because the more that particles are pulled down bygravity, the greater the tendency for the particles to migrate to regions of lower concentration. The flux is given byFick's law,whereJ =ρv. Introducing the formula forρ, we find that
In a state of dynamical equilibrium, this speed must also be equal tov =μmg. Both expressions forv are proportional tomg, reflecting that the derivation is independent of the type of forces considered. Similarly, one can derive an equivalent formula for identicalcharged particles of chargeq in a uniformelectric field of magnitudeE, wheremg is replaced with theelectrostatic forceqE. Equating these two expressions yields theEinstein relation for the diffusivity, independent ofmg orqE or other such forces:Here the first equality follows from the first part of Einstein's theory, the third equality follows from the definition of theBoltzmann constant askB =R /NA, and the fourth equality follows from Stokes's formula for the mobility. By measuring the mean squared displacement over a time interval along with the universal gas constantR, the temperatureT, the viscosityη, and the particle radiusr, the Avogadro constantNA can be determined.
The type of dynamical equilibrium proposed by Einstein was not new. It had been pointed out previously byJ. J. Thomson[21] in his series of lectures at Yale University in May 1903 that the dynamic equilibrium between the velocity generated by aconcentration gradient given by Fick's law and the velocity due to the variation of the partial pressure caused when ions are set in motion "gives us a method of determining Avogadro's constant which is independent of any hypothesis as to the shape or size of molecules, or of the way in which they act upon each other".[21]
An identical expression to Einstein's formula for the diffusion coefficient was also found byWalther Nernst in 1888[22] in which he expressed the diffusion coefficient as the ratio of theosmotic pressure to the ratio of thefrictional force and the velocity to which it gives rise. The former was equated to the law of van 't Hoff while the latter was given byStokes's law. He writes for the diffusion coefficientk′, where is the osmotic pressure andk is the ratio of the frictional force to the molecular viscosity which he assumes is given by Stokes's formula for the viscosity. Introducing theideal gas law per unit volume for the osmotic pressure, the formula becomes identical to that of Einstein's.[23] The use of Stokes's law in Nernst's case, as well as in Einstein and Smoluchowski, is not strictly applicable since it does not apply to the case where the radius of the sphere is small in comparison with themean free path.[24]
Confirming Einstein's formula experimentally proved difficult.Initial attempts byTheodor Svedberg in 1906 and 1907 were critiqued by Einstein and by Perrin as not measuring a quantity directly comparable to the formula.Victor Henri in 1908 took cinematographic shots through a microscope and found quantitative disagreement with the formula but again the analysis was uncertain.[25] Einstein's predictions were finally confirmed in a series of experiments carried out by Chaudesaigues in 1908 and Perrin in 1909.[26][11] The confirmation of Einstein's theory constituted empirical progress for thekinetic theory of heat. In essence, Einstein showed that the motion can be predicted directly from the kinetic model ofthermal equilibrium. The importance of the theory lay in the fact that it confirmed the kinetic theory's account of thesecond law of thermodynamics as being an essentially statistical law.[27]
Brownian motion model of the trajectory of a particle of dye in water
Smoluchowski's theory of Brownian motion[28] starts from the same premise as that of Einstein and derives the same probability distributionρ(x,t) for the displacement of a Brownian particle along thex in timet. He therefore gets the same expression for the mean squared displacement:. However, when he relates it to a particle of massm moving at a velocityu which is the result of a frictional force governed by Stokes's law, he findswhereμ is the viscosity coefficient, anda is the radius of the particle. Associating the kinetic energy with the thermal energyRT/N, the expression for the mean squared displacement is64/27 times that found by Einstein. The fraction 27/64 was commented on byArnold Sommerfeld in his necrology on Smoluchowski: "The numerical coefficient of Einstein, which differs from Smoluchowski by 27/64 can only be put in doubt."[29]
Smoluchowski[30] attempts to answer the question of why a Brownian particle should be displaced by bombardments of smaller particles when the probabilities for striking it in the forward and rear directions are equal.If the probability ofm gains andn −m losses follows abinomial distribution,with equala priori probabilities of 1/2, the mean total gain is
Ifn is large enough so that Stirling's approximation can be used in the formthen the expected total gain will be[citation needed]showing that it increases as the square root of the total population.
Suppose that a Brownian particle of massM is surrounded by lighter particles of massm which are traveling at a speedu. Then, reasons Smoluchowski, in any collision between a surrounding and Brownian particles, the velocity transmitted to the latter will bemu/M. This ratio is of the order of10−7 cm/s. But we also have to take into consideration that in a gas there will be more than 1016 collisions in a second, and even greater in a liquid where we expect that there will be 1020 collision in one second. Some of these collisions will tend to accelerate the Brownian particle; others will tend to decelerate it. If there is a mean excess of one kind of collision or the other to be of the order of 108 to 1010 collisions in one second, then velocity of the Brownian particle may be anywhere between10–1000 cm/s. Thus, even though there are equal probabilities for forward and backward collisions there will be a net tendency to keep the Brownian particle in motion, just as the ballot theorem predicts.
These orders of magnitude are not exact because they do not take into consideration the velocity of the Brownian particle,U, which depends on the collisions that tend to accelerate and decelerate it. The largerU is, the greater will be the collisions that will retard it so that the velocity of a Brownian particle can never increase without limit. Could such a process occur, it would be tantamount to a perpetual motion of the second type. And since equipartition of energy applies, the kinetic energy of the Brownian particle,, will be equal, on the average, to the kinetic energy of the surrounding fluid particle,.
In 1906, Smoluchowski published a one-dimensional model to describe a particle undergoing Brownian motion.[31] The model assumes collisions withM ≫m whereM is the test particle's mass andm the mass of one of the individual particles composing the fluid. It is assumed that the particle collisions are confined to one dimension and that it is equally probable for the test particle to be hit from the left as from the right. It is also assumed that every collision always imparts the same magnitude ofΔV. IfNR is the number of collisions from the right andNL the number of collisions from the left then afterN collisions the particle's velocity will have changed byΔV(2NR −N). Themultiplicity is then simply given by:and the total number of possible states is given by2N. Therefore, the probability of the particle being hit from the rightNR times is:
As a result of its simplicity, Smoluchowski's 1D model can only qualitatively describe Brownian motion. For a realistic particle undergoing Brownian motion in a fluid, many of the assumptions don't apply. For example, the assumption that on average there are an equal number of collisions from the right as from the left falls apart once the particle is in motion. Also, there would be a distribution of different possibleΔVs instead of always just one in a realistic situation.
Thediffusion equation yields an approximation of the time evolution of theprobability density function associated with the position of the particle going under a Brownian movement under the physical definition. The approximation becomes valid on timescales much larger than the timescale of individual atomic collisions, since it does not include a term to describe the acceleration of particles during collision. The time evolution of the position of the Brownian particle over all time scales described using theLangevin equation, an equation that involves a random force field representing the effect of thethermal fluctuations of the solvent on the particle.[16] At longer times scales, where acceleration is negligible, individual particle dynamics can be approximated usingBrownian dynamics in place ofLangevin dynamics.
Instellar dynamics, a massive body (star,black hole, etc.) can experience Brownian motion as it responds togravitational forces from surrounding stars.[32] The rms velocityV of the massive object, of massM, is related to the rms velocity of the background stars bywhere is the mass of the background stars. The gravitational force from the massive object causes nearby stars to move faster than they otherwise would, increasing both andV.[32] The Brownian velocity ofSgr A*, thesupermassive black hole at the center of theMilky Way galaxy, is predicted from this formula to be less than 1 km s−1.[33]
An animated example of a Brownian motion-likerandom walk on a 2D surface with periodic boundary conditions. In thescaling limit, random walk approaches the Wiener process according toDonsker's theorem.
An alternative characterisation of the Wiener process is the so-calledLévy characterisation that says that the Wiener process is an almost surely continuousmartingale withW0 = 0 andquadratic variation.
A third characterisation is that the Wiener process has a spectral representation as a sine series whose coefficients are independent random variables. This representation can be obtained using theKosambi–Karhunen–Loève theorem.
The Wiener process can be constructed as thescaling limit of arandom walk, or other discrete-time stochastic processes with stationary independent increments. This is known asDonsker's theorem. Like the random walk, the Wiener process is recurrent in one or two dimensions (meaning that it returns almost surely to any fixedneighborhood of the origin infinitely often) whereas it is not recurrent in dimensions three and higher. Unlike the random walk, it isscale invariant.A d-dimensionalGaussian free field has been described as "a d-dimensional-time analog of Brownian motion."[35]
The French mathematicianPaul Lévy proved the following theorem, which gives a necessary and sufficient condition for a continuousRn-valued stochastic processX to actually ben-dimensional Brownian motion. Hence, Lévy's condition can actually be used as an alternative definition of Brownian motion.
LetX = (X1, ...,Xn) be a continuous stochastic process on aprobability space(Ω, Σ,P) taking values inRn. Then the following are equivalent:
X is a Brownian motion with respect toP, i.e., the law ofX with respect toP is the same as the law of ann-dimensional Brownian motion, i.e., thepush-forward measureX∗(P) isclassical Wiener measure onC0([0, ∞);Rn).
The spectral content of a stochastic process can be found from thepower spectral density, formally defined aswhere stands for theexpected value. The power spectral density of Brownian motion is found to be[38]whereD is thediffusion coefficient ofXt. For naturally occurring signals, the spectral content can be found from the power spectral density of a single realization, with finite available time, i.e.,which for an individual realization of a Brownian motion trajectory,[39] it is found to have expected valueandvariance[39]
For sufficiently long realization times, the expected value of the power spectrum of a single trajectory converges to the formally defined power spectral density, but its coefficient of variation tends to. This implies the distribution of is broad even in the infinite time limit.
Brownian motion is usually considered to take place inEuclidean space. It is natural to consider how such motion generalizes to more complex shapes, such assurfaces or higher dimensionalmanifolds. The formalization requires the space to possess some form of aderivative, as well as ametric, so that aLaplacian can be defined. Both of these are available onRiemannian manifolds.
Riemannian manifolds have the property thatgeodesics can be described inpolar coordinates; that is, displacements are always in a radial direction, at some given angle. Uniform random motion is then described by Gaussians along the radial direction, independent of the angle, the same as in Euclidean space.
Thenarrow escape problem is a ubiquitous problem in biology, biophysics and cellular biology which has the following formulation: a Brownian particle (ion,molecule, orprotein) is confined to a bounded domain (a compartment or a cell) by a reflecting boundary, except for a small window through which it can escape. The narrow escape problem is that of calculating the mean escape time. This time diverges as the window shrinks, thus rendering the calculation asingular perturbation problem.
^Meyburg, Jan Philipp; Diesing, Detlef (2017). "Teaching the Growth, Ripening, and Agglomeration of Nanostructures in Computer Experiments".Journal of Chemical Education.94 (9):1225–1231.Bibcode:2017JChEd..94.1225M.doi:10.1021/acs.jchemed.6b01008.
^abFeynman, Richard (1964)."The Brownian Movement".The Feynman Lectures of Physics, Volume I. p. 41.
^Pathria, RK (1972). Statistical Mechanics. Pergamon Press. pp. 43–48, 73–74. ISBN 0-08-016747-0.
^Tsekov, Roumen (1995). "Brownian motion of molecules: the classical theory".Ann. Univ. Sofia.88: 57.arXiv:1005.1490.Bibcode:1995AUSFC..88...57T.the behavior of a Brownian particle is quite irregular and can be described only in the frames of a statistical approach.
^Davis, Mark H. A.; Bachelier, Louis; Etheridge, Alison (2011).Louis Bachelier's Theory of Speculation: The Origins of Modern Finance. Princeton: Princeton University Press.ISBN978-1-4008-2930-9.
^See p. 535 inSommerfeld, A. (1917). "Zum Andenken an Marian von Smoluchowski" [In Memory of Marian von Smoluchowski].Physikalische Zeitschrift (in German).18 (22):533–539.
Chaudesaigues, M. (1908). "Le mouvement brownien et la formule d'Einstein" [Brownian motion and Einstein's formula].Comptes Rendus (in French).147:1044–6.
Henri, V. (1908). "Études cinématographique du mouvement brownien" [Cinematographic studies of Brownian motion].Comptes Rendus (in French) (146):1024–6.
Perrin, J. (1909). "Mouvement brownien et réalité moléculaire" [Brownian movement and molecular reality].Annales de chimie et de physique. 8th series.18:5–114.
Danish version: "Om Anvendelse af mindste Kvadraters Methode i nogle Tilfælde, hvor en Komplikation af visse Slags uensartede tilfældige Fejlkilder giver Fejlene en 'systematisk' Karakter".
French version: "Sur la compensation de quelques erreurs quasi-systématiques par la méthodes de moindre carrés" published simultaneously inVidensk. Selsk. Skr. 5. Rk., naturvid. og mat. Afd., 12:381–408, 1880.