Mean squared displacement explained

In statistical mechanics, the mean squared displacement (MSD, also mean square displacement, average squared displacement, or mean square fluctuation) is a measure of the deviation of the position of a particle with respect to a reference position over time. It is the most common measure of the spatial extent of random motion, and can be thought of as measuring the portion of the system "explored" by the random walker. In the realm of biophysics and environmental engineering, the Mean Squared Displacement is measured over time to determine if a particle is spreading slowly due to diffusion, or if an advective force is also contributing.^[1] Another relevant concept, the variance-related diameter (VRD, which is twice the square root of MSD), is also used in studying the transportation and mixing phenomena in the realm of environmental engineering.^[2] It prominently appears in the Debye–Waller factor (describing vibrations within the solid state) and in the Langevin equation (describing diffusion of a Brownian particle).

The MSD at time

t is defined as an ensemble average:

$$\backslash text\; \backslash equiv\; \backslash left\backslash langle\; \backslash left|\backslash mathbf(t)-\backslash mathbf\backslash right|^2\backslash right\backslash rangle=\; \backslash frac\backslash sum\_^N\; \backslash left|\backslash mathbf(t)\; -\; \backslash mathbf(0)\backslash right|^2$$

where N is the number of particles to be averaged, vector

x ^{(i)} ( 0 ) = x ^{(i)} _{0} is the reference position of the i -th particle, and vector x ^{(i)} (t) is the position of the i -th particle at time t.^[3]

Derivation of the MSD for a Brownian particle in 1D

The probability density function (PDF) for a particle in one dimension is found by solving the one-dimensional diffusion equation. (This equation states that the position probability density diffuses out over time - this is the method used by Einstein to describe a Brownian particle. Another method to describe the motion of a Brownian particle was described by Langevin, now known for its namesake as the Langevin equation.)$$\backslash frac\; =\; D\backslash frac,$$given the initial condition

p(x,t =0 \mid x _{0)} = \delta(x - x _{0)} ; where x(t) is the position of the particle at some given time, x _{0} is the tagged particle's initial position, and D is the diffusion constant with the S.I. units m ^{2s} ^{-1} (an indirect measure of the particle's speed). The bar in the argument of the instantaneous probability refers to the conditional probability. The diffusion equation states that the speed at which the probability for finding the particle at x(t) is position dependent.

The differential equation above takes the form of 1D heat equation. The one-dimensional PDF below is the Green's function of heat equation (also known as Heat kernel in mathematics):$$P(x,t)\; =\; \backslash frac\; \backslash exp\backslash left(-\backslash frac\backslash right).$$This states that the probability of finding the particle at

x(t) is Gaussian, and the width of the Gaussian is time dependent. More specifically the full width at half maximum (FWHM)(technically/pedantically, this is actually the Full duration at half maximum as the independent variable is time) scales like$$\backslash text\backslash sim\backslash sqrt.$$Using the PDF one is able to derive the average of a given function, L , at time t :$$\backslash langle\; L(t)\; \backslash rangle\; \backslash equiv\; \backslash int^\backslash infty\_\; L(x,t)\; P(x,t)\; \backslash ,\; dx,$$where the average is taken over all space (or any applicable variable).

The Mean squared displacement is defined as$$\backslash text\; \backslash equiv\; \backslash left\backslash langle\; \backslash left(x(t)\; -\; x\_0\; \backslash right)^2\; \backslash right\backslash rangle,$$expanding out the ensemble average$$\backslash left\backslash langle\; \backslash left(x\; -\; x\_0\backslash right)^2\; \backslash right\backslash rangle\; =\; \backslash left\backslash langle\; x^2\; \backslash right\backslash rangle\; +\; x\_0^2\; -\; 2\; x\_0\; \backslash langle\; x\; \backslash rangle,$$dropping the explicit time dependence notation for clarity. To find the MSD, one can take one of two paths: one can explicitly calculate

\langle x ^{2\rangle} and \langle x\rangle , then plug the result back into the definition of the MSD; or one could find the moment-generating function, an extremely useful, and general function when dealing with probability densities. The moment-generating function describes the moment of the PDF. The first moment of the displacement PDF shown above is simply the mean: \langle x\rangle . The second moment is given as \langle x ^{2\rangle} .

So then, to find the moment-generating function it is convenient to introduce the characteristic function:$$G(k)\; =\; \backslash langle\; e^\backslash rangle\backslash equiv\; \backslash int\_I\; e^P(x,t\backslash mid\; x\_0)\; \backslash ,\; dx,$$one can expand out the exponential in the above equation to give$$G(k)\; =\; \backslash sum^\backslash infty\_\backslash frac\backslash mu\_m.$$By taking the natural log of the characteristic function, a new function is produced, the cumulant generating function,$$\backslash ln(G(k))\; =\; \backslash sum^\backslash infty\_\backslash frac\backslash kappa\_m,$$where

\kappa _{m} is the cumulant of x . The first two cumulants are related to the first two moments, \mu , via \kappa _{1} = \mu _{1;} and \kappa _{2} = \mu _{2-\mu} ^{2 ,} ^{1} where the second cumulant is the so-called variance, \sigma ^{2} . With these definitions accounted for one can investigate the moments of the Brownian particle PDF,$$G(k)\; =\; \backslash frac\backslash int\_I\; \backslash exp\backslash left(i\; k\; x\; -\; \backslash frac\backslash right)\; \backslash ,\; dx;$$by completing the square and knowing the total area under a Gaussian one arrives at$$G(k)\; =\; \backslash exp(i\; k\; x\_0\; -\; k^2\; D\; t).$$Taking the natural log, and comparing powers of ik to the cumulant generating function, the first cumulant is$$\backslash kappa\_1\; =\; x\_0,$$which is as expected, namely that the mean position is the Gaussian centre. The second cumulant is$$\backslash kappa\_2\; =\; 2Dt,\; \backslash ,$$the factor 2 comes from the factorial factor in the denominator of the cumulant generating function. From this, the second moment is calculated,$$\backslash mu\_2\; =\; \backslash kappa\_2\; +\; \backslash mu\_1^2\; =\; 2Dt\; +\; x\_0^2.$$Plugging the results for the first and second moments back, one finds the MSD,$$\backslash left\backslash langle\; \backslash left(x(t)\; -\; x\_0\backslash right)^2\; \backslash right\backslash rangle\; =\; 2Dt.$$

Derivation for n dimensions

For a Brownian particle in higher-dimension Euclidean space, its position is represented by a vector

x = (x _{1,x} _{2,\ldots,x} _{n)} , where the Cartesian coordinates x _{1,x} _{2,\ldots,x} _{n} are statistically independent.

The n-variable probability distribution function is the product of the fundamental solutions in each variable; i.e.,

$$P(\backslash mathbf,t)\; =\; P(x\_1,t)P(x\_2,t)\backslash dots\; P(x\_n,t)=\backslash frac\backslash exp\; \backslash left\; (-\backslash frac\; \backslash right).$$

The Mean squared displacement is defined as

$$\backslash mathrm\; \backslash equiv\; \backslash left\backslash langle\; |\backslash mathbf-\backslash mathbf|^2\backslash right\backslash rangle=\; \backslash left\backslash langle\; \backslash left(x\_1(t)-x\_1(0)\backslash right)^2\; +\; \backslash left(x\_2(t)-x\_2(0)\backslash right)^2\; +\; \backslash dots\; +\; \backslash left(x\_n(t)-x\_n(0)\backslash right)^2\; \backslash right\backslash rangle$$

Since all the coordinates are independent, their deviation from the reference position is also independent. Therefore,

$$\backslash text=\; \backslash left\backslash langle\; \backslash left(x\_1(t)-x\_1(0)\backslash right)^2\; \backslash right\backslash rangle\; +\; \backslash left\backslash langle\; \backslash left(x\_2(t)-x\_2(0)\backslash right)^2\; \backslash right\backslash rangle\; +\; \backslash dots\; +\; \backslash left\backslash langle\; \backslash left(x\_n(t)-x\_n(0)\backslash right)^2\; \backslash right\backslash rangle$$

For each coordinate, following the same derivation as in 1D scenario above, one obtains the MSD in that dimension as

2 Dt . Hence, the final result of mean squared displacement in n-dimensional Brownian motion is:

$$\backslash text\; =\; 2nDt.$$

Definition of MSD for time lags

In the measurements of single particle tracking (SPT), displacements can be defined for different time intervals between positions (also called time lags or lag times). SPT yields the trajectory

\vec r(t) = [x(t),y(t)] , representing a particle undergoing two-dimensional diffusion.

Assuming that the trajectory of a single particle measured at time points

1 \Delta t, 2 \Delta t,\ldots,N\Delta t , where \Delta t is any fixed number, then there are N(N -1 )/ 2 non-trivial forward displacements \vec d _{ij} = \vec r _{j} - \vec r _{i} (1 \leqslant i < j \leqslant N , the cases when i = j are not considered) which correspond to time intervals (or time lags) \Delta t _{ij} = (j - i)\Delta t . Hence, there are many distinct displacements for small time lags, and very few for large time lags, {\rm MSD} can be defined as an average quantity over time lags:^{[4] [5]}

$$\backslash overline\; =\backslash frac\; 1\; \backslash sum\_^)^2\; \backslash qquad\; n=1,\backslash ldots,N-1.$$

Similarly, for continuous time series :

$$\backslash overline\; =\; \backslash frac\; 1\; \backslash int\_0^\; [r(t\; +\; \backslash Delta)\; -\; r(t)]^2\; \backslash ,\; dt$$

It's clear that choosing large

T and \Delta \ll T can improve statistical performance. This technique allow us estimate the behavior of the whole ensembles by just measuring a single trajectory, but note that it's only valid for the systems with ergodicity, like classical Brownian motion (BM), fractional Brownian motion (fBM), and continuous-time random walk (CTRW) with limited distribution of waiting times, in these cases, \overline {\delta ^{2(\Delta)}} = \left\langle [r(t) - r( 0 )] ^{2} \right\rangle (defined above), here \left\langle ⋅ \right\rangle denotes ensembles average. However, for non-ergodic systems, like the CTRW with unlimited waiting time, waiting time can go to infinity at some time, in this case, \overline {\delta ^{2(\Delta} )} strongly depends on T , \overline{\delta ^{2(\Delta)}} and \left\langle [r(t) - r( 0 )] ^{2} \right\rangle don't equal each other anymore, in order to get better asymptotics, introduce the averaged time MSD:

$$\backslash left\backslash langle\; \backslash right\backslash rangle\; =\; \backslash frac\; \backslash sum\; \backslash overline$$

Here

\left\langle ⋅ \right\rangle denotes averaging over N ensembles.

Also, one can easily derive the autocorrelation function from the MSD:

$$\backslash left\backslash langle\; \backslash right\backslash rangle\; =\; \backslash left\backslash langle\; r^2(t)\; \backslash right\backslash rangle\; +\; \backslash left\backslash langle\; r^2(0)\; \backslash right\backslash rangle\; -\; 2\backslash left\backslash langle\; r(t)r(0)\; \backslash right\backslash rangle,$$where

\left\langle r(t)r( 0 ) \right\rangle is so-called autocorrelation function for position of particles.

MSD in experiments

Experimental methods to determine MSDs include neutron scattering and photon correlation spectroscopy.

The linear relationship between the MSD and time t allows for graphical methods to determine the diffusivity constant D. This is especially useful for rough calculations of the diffusivity in environmental systems. In some atmospheric dispersion models, the relationship between MSD and time t is not linear. Instead, a series of power laws empirically representing the variation of the square root of MSD versus downwind distance are commonly used in studying the dispersion phenomenon.^[6]

See also

• Root-mean-square deviation of atomic positions: the average is taken over a group of particles at a single time, where the MSD is taken for a single particle over an interval of time
• Mean squared error

Notes and References

1. Tarantino. Nadine. Tinevez. Jean-Yves. Crowell. Elizabeth Faris. Boisson. Bertrand. Henriques. Ricardo. Mhlanga. Musa. Agou. Fabrice. Israël. Alain. Laplantine. Emmanuel. 2014-01-20. TNF and IL-1 exhibit distinct ubiquitin requirements for inducing NEMO–IKK supramolecular structures. J Cell Biol. en. 204. 2. 231–245. 10.1083/jcb.201307172. 0021-9525. 3897181. 24446482.
2. Book: B., Fischer, Hugo. Mixing in inland and coastal waters. 1979-01-01. Academic Press. 9780080511771. 983391285.
3. Frenkel, Daan & Smit, Berend. Understanding molecular simulation: From algorithms to applications. Academic Press, 196 (2nd Ed.), p. 97.
4. Michalet . Xavier . Mean square displacement analysis of single-particle trajectories with localization error: Brownian motion in an isotropic medium . Physical Review E . 20 October 2010 . 82 . 4 . 041914 . 10.1103/PhysRevE.82.041914 . 21230320 . 3055791 . 2010PhRvE..82d1914M .
5. Qian . H. . Sheetz . M. P. . Elson . E. L. . Single particle tracking. Analysis of diffusion and flow in two-dimensional systems . Biophysical Journal . 1 October 1991 . 60 . 4 . 910–921 . 10.1016/S0006-3495(91)82125-7 . 1742458 . 1260142 . 1991BpJ....60..910Q . en . 0006-3495.
6. Davidson. G. A.. 1990-08-01. A Modified Power Law Representation of the Pasquill-Gifford Dispersion Coefficients. Journal of the Air & Waste Management Association. 40. 8. 1146–1147. 10.1080/10473289.1990.10466761. 1047-3289.