Confined Brownian Motion Tracked With Motion Blur: Estimating Diffusion Coefficient and Size of Confining Space

Mesoscopic environments and particles diffusing in them are often studied by tracking such particles individually while their Brownian motion explores their environment. Environments may be, e.g., a domain in a cell membrane, an interior compartment of a cell, or an engineered nanopit. Particle trajectories are typically determined from time-lapse recorded movies. These are recorded with sufficient exposure time per frame to be able to detect and localize particles in each frame. Since particles move during this exposure time, particles image with motion blur. This motion blur can compromise estimates of diffusion coefficients and the size of the confining domain if not accounted for correctly. We do that here. We give explicit and exact expressions for the variance of measured positions and the mean-squared displacement of a Brownian particle confined in, respectively, a 1D box, a 2D box, a 2D circular disc, and a 3D sphere. Our expressions are valid for all exposure times, irrespective of the size of the confining space and the value of the diffusion coefficient. They apply also in the common case where the exposure time is smaller than the time-lapse due, e.g., to “dead time” caused by the readout process in the camera. These expressions permit determination of diffusion coefficients and domain sizes for given movies for the simple geometries we consider. More important, the trends observed in our exact results when parameter values are varied are valid also for more complex geometries for which no exact analytical solutions exist. Wherever the underlying physics is the same, the exact quantitative description of its consequences provided here is portable as a qualitative and semi-quantitative understanding of its consequences in general. The results may also be useful for other types of reflected Brownian motion than those occurring in single-particle tracking, e.g., in nuclear magnetic resonance imaging techniques. For use in that particular context, we briefly discuss the effects of confinement on anisotropic Brownian motion imaged with motion blur.


INTRODUCTION
In many branches of science, individual particles or molecules are labeled fluorescently in order to track and characterize their motion. To this end, they must be sufficiently isolated from each other in time and/or space, e.g., by sparse labeling [1] or super-resolution microscopies [2]. Only then may their centers be found with a precision that increases with the number of photons recorded [3]. When it does, nanometer precision is obtained routinely. In this manner, single-particle tracking is used extensively to uncover, e.g., the organization and dynamics of biology at its shortest length scale [4] and its interactions with nanoparticles [5]. One key advantage of single-particle tracking is that, with sufficiently long trajectories, it yields single-particle results. With sample-averaging unnecessary, it can detect heterogeneity in a population. Moreover, if such heterogeneity is absent, it can detect anomalous behavior in individual particles. Thus spotted before averaging, such outliers can be excluded from samples before sample-averages are calculated. This can improve accuracy and precision of estimates substantially.
Since the molecules and particles are studied under ambient conditions, they undergo Brownian motion. Various methods of various qualities are used to estimate diffusion coefficients from recorded particle trajectories describing free Brownian motion [6,7]. Here, we address how to do similar estimates for Brownian motion in confined spaces: for example, a domain in a cell membrane, a compartment inside a cell, or an engineered nanopit [4,[8][9][10]. Such motion in confined domains and the properties of the domains themselves are important. Indeed, nano-domains in cell membranes, such as lipid nano-domains and nano-domains created by cellular filaments, are involved in multiple biological mechanisms, such as signal processing, membrane trafficking, and various diseases [11][12][13]. Moreover, single-particle tracking of nano-particles is often used to characterize lab-on-a-chip devices and thus the importance of single-particle tracking for characterization of particles and their confining environments.
The implementation of single-particle tracking may, however, itself confound results obtained with it. In order to track a fluorescently labeled particle or molecule, one records a timelapse movie of micrographs of the particle. Typically, the dim signal from its fluorescent label is recorded with a highly sensitive camera. Each frame of the time-lapse movie must have recorded enough photons from the label to permit detection and localization of the particle in that frame [14]. So, a sufficiently long exposure time is used. The particle, however, moves during exposure. Since typical exposure times last tens of milliseconds, they are much longer than the characteristic time-scale of the particle's inertial motion. The particle thus undergoes Brownian motion while being imaged. This creates motion blur, and the observed position essentially is the time-average over positions visited by the particle during exposure. This has consequences for the statistics of motion, such as variances of measured positions of the particle and its mean-squared displacement (MSD).
While the effects of motion blur on the statistics of free Brownian motion have been studied in detail [6,[14][15][16], less effort has been devoted to its effects in the case of confined Brownian motion. Clearly, if the exposure time is very long, the particle visits all of a confining domain many times during a single exposure time. Thus, if the domain is smaller than the diffraction limit, super-localization techniques will locate the particle at the centroid of the confining domain if at all. Such consequences of motion blur have been observed [17], and approximate formulas that correct for these effects have been developed for specific geometries [18]. Their quality, unfortunately, depends on the values of experimental parameters. Therefore, we here derive exact expressions for the variances of positions and the MSD for Brownian particles in confinement recorded with motion blur. Specifically, we treat the cases of 1D Brownian motion confined to a box, 2D Brownian motion confined, respectively, to a box and to a circular disc, and 3D Brownian motion confined to a sphere. Our formulas are valid for all frame rates, irrespective of the value of the diffusion coefficient and the size of the confining domain. We also model the common case of "dead time" making the exposure time shorter than the inverse frame rate. While our exact results are derived for uniform and isotropic Brownian motion, we also discuss the implications of confinement and motion blur on anisotropic Brownian motion confined to a 2D disc. We do this using simulated trajectories. We compare the statistics obtained from those trajectories to the theory for the isotropic case.
We note that confined (sometimes called reflected) Brownian motion also plays an important role in other research fields, e.g., in probability theory and nuclear magnetic resonance (NMR) imaging. Time-averaging also may influence results there, e.g., in pulsed field-gradient spin-echo NMR, where the diffusion of the spins during the finite duration of the pulse causes that the pores in materials appear smaller than they are (for a review, see [19, Sec. II.B]).

True Positions Versus Measured Positions
In order to track a fluorescently labeled particle in a fluorescence microscopy experiment, one records a movie of it. Typically, one uses a constant time-lapse, Δt, between frames. 1 If a particle's true position was x 0 at time t 0 , we let x(t|x 0 , t 0 ) denote its true position at time t. If positions could be recorded instantaneously, we could sample positions x j x(t j |x 0 , t 0 ) at equidistant discrete points in time, t j jΔt. That is, however, not how positions are recorded in camera-based single-particle tracking. Instead, Δt is typically divided into two complementary parts, Δt open and Δt closed . Thus, Δt Δt open + Δt closed . The camera shutter is open during a finite exposure time, Δt open , in order to record a sufficient number of photons to determine the particle's position with the desired precision [3]. The time interval Δt closed may be a necessary "dead time" between frames due to finite duration of readout processes in the camera. In general, Δt open ≤ Δt, and equality is achieved only if the camera shutter is left open during the entire time-lapse.
Finite exposure time can be accounted for by introducing the shutter function s j (t), a non-negative function that integrates to unity over the time interval Δt [15]. For the present case, the shutter function is a step function with s j (t) 1/Δt open for t ∈ [t j − Δt open , t j ] and 0 otherwise. More advanced shutter functions can be considered [15], but they are less relevant in practice.
Since the particle moves during the recording of its position, its measured position is affected by motion blur. That position, x msr (t j |x j−1 , t j−1 ), is found as Δt open tj tj−Δtopen x t|x j−1 , t j−1 dt, (1) provided that it started at position x j−1 at t j−1 and that the estimated position is the temporal average over the positions in the trajectory of the particle, x(t|x j−1 , t j−1 ), during the exposure time Δt open . The difference between the true positions, x j , and the measured positions, x msr (t j ), is illustrated in Figure 1A for free diffusion.
Below, we explore the consequences of this time-averaging on the variance and the mean-squared displacement (MSD) for Brownian motion in confining spaces. For free Brownian motion, these effects are discussed at length in, e.g., [6,14,15].
In real measurements, the measured position of a particle also contains a localization error due to photon shot noise. In practice, localization errors may be important in the analysis of experimental data, but we present our formulas without localization errors to keep them simple. In Section 4.5, we briefly describe how to account also for localization errors.

Free Brownian Motion
We consider a particle embedded in a fluid at rest at finite temperature, subject only to the fluctuating thermal forces from the fluid surrounding it. We assume that the time-lapse Δt, with which the particle is observed, is long relatively to the inertial time-scale of the particle's thermal motion. Then, Einstein's original description of free Brownian motion applies. In this case, the time-dependence of the probability density P to find the Brownian particle at position x at time t, given that it was at position x 0 at time t 0 , is governed by Fick's second law of diffusion (1855), the Fokker-Planck equation [20], with the initial condition P(x, t 0 |x 0 , t 0 ) δ(x − x 0 ). Here, D is the diffusion constant.
With the boundary conditions P(x, t|x 0 , t 0 ) → 0 for x → ± ∞, for all t > t 0 ( closed boundaries at x ± ∞), the solution to Eq. 2 is That is, the conditional probability density is normal with mean x 0 and variance σ 2 2D(t − t 0 ).

Statistics of Displacements for Instantaneously Sampled Positions
As a consequence, consecutive sampled positions of a Brownian particle are related by displacements Δx j ≡ x j − x j−1 that are normally distributed with mean and covariance, An example of time-lapse sampled positions of a Brownian particle is shown in Figure 1A.
Equation 2 also gives the mean-squared displacement, with its signature proportionality to the time lag nΔt with constant of proportionality 2D. This property is often used to estimate the diffusion coefficient from trajectory data. To this end, 2D is determined as the slope in a fit of a straight line to the MSD estimated from data. This and other methods for determination of D are still discussed in the literature [7,21], but for trajectories with signal-to-noise ration SNR DΔt √ /σ loc > 1, D is certainly determined optimally by comparing the sample variance of the particle's positions with Eq. 5 [6,7].

Motion Blur Reduces Variances of Displacements and Correlates Them
Since motion blur effectively acts as a low-pass filter, the statistical properties of displacements Δx msr, j x msr, j − x msr, j−1 for measured positions are different from those of the true positions Δx j x j − x j−1 . For the shutter function s j (t) introduced above, the displacements of the measured positions have the properties [6,7,15] 〈 Δx msr, 〈Δx msr, j Δx msr, ℓ 〉 0 for j − ℓ > 1. (10) So, the motion blur reduces the variance of the displacements and introduces a nearest-neighbor correlation between them. Both effects are proportional to Δt open . The variance decreases from its maximal value, 2DΔt, at Δt open 0 to its minimum, 4DΔt/3, at Δt open Δt, see Figure 1B. Consequently, the value for D is underestimated even in the absence of localization errors if motion blur is present in measurements due to finite exposure time but is ignored, and the sample variance of the displacements of the measured positions is used as a naive estimate of the diffusion coefficient D.
The MSD of the measured positions for n ≥ 1, however, is still linear in the time lag nΔt with slope 2D, though offset by a constant due to the motion blur, Consequently, a fit of a straight line to the MSD estimated from positions subject to motion blur will return the true value of D on average. However, unless data are very rich, superior alternatives for the determination of D exist also for this case [6,7].

Brownian Motion in 1D Confinement
To account for the statistical effects of confinement, we first consider the simplest possible case, a particle diffusing in 1D in a box of length L. Let x denote the coordinate of the particle with axis chosen such that the box has x ∈ [0, L]. Equation 2 holds also in this box, but with closed boundary conditions, i.e., For initial condition P(x, t 0 |x 0 , t 0 ) δ(x − x 0 ) with 0 < x 0 < L, P(x, t|x 0 , t 0 ) is found by separation of variables [22] and is Here, we have introduced the characteristic time-scale τ ≡ (L/2) 2 /D. Below, we use Eq. 13 to calculate an exact expression for the variance and the MSD for positions measured using finite exposure time.
Brownian motion confined to a 2D square or rectangular box is described by the independent motion of each of two Cartesian coordinates in their respective 1D boxes.

Variance of Positions Recorded With Motion Blur for a Brownian Particle in a 1D Box
In the limiting case of instantaneous recording of exact positions, there is uniform probability density 1/L to observe the particle anywhere in the 1D box ( Figures 1D, 2A). Consequently, the variance of the particle's position is var(x msr ) L 2 /12. However, for Δt open /τ non-negligible, the Frontiers in Physics | www.frontiersin.org January 2021 | Volume 8 | Article 583202 distribution of measured particle positions depends on the value of Δt open /τ ( Figures 2B,C). This is because motion-blurred positions are averages over consecutive true positions, and it is unlikely that all contributions to such an average are close to a boundary when Δt open /τ is non-negligible. As a result, the confining domain appears smaller than it is [17]. On the other hand, if Δt open /τ ≫ 1, the particle visits most of the confining domain during each exposure and all measured positions thus appear in the center of the box with a vanishing positional variance, var(x msr )x0, as consequence.
In practice, Δt open /τ is non-negligible. In order to obtain expressions for, respectively, the variance and the MSD for a particle in confinement including the effect of motion blur, valid for all values of Δt open /τ, we follow Wong and Halvorsen [23], who used a similar calculation for a particle in a harmonic potential. In the process, we also account for the possibility that Δt open < Δt.
The variance of the measured positions is where the average 〈 · 〉 is over all possible trajectories. Here, 〈x msr 〉 L/2, due to symmetry. Assuming a uniform distribution of initial positions, i.e., p(x 0 ) 1/L for 0 < x 0 < L and zero otherwise, the variance is where 〈 · 〉 0 means averaging over all trajectories starting at the position x 0 . We then insert x msr (t|x 0 ) from Eq. 1 and write The time-ordered auto-correlation function can be expressed via the conditional probability density in Eq. 13, The calculation is straightforward, but lengthy (Supplementary Material). Using Eqs. 15, 16, 17, we get that the variance for a particle in 1D confinement is This expression is exact and agrees with simulated data for Brownian motion in 1D confinement, as shown in Figure 2D. The result for the variance was also derived in [24] in the context of the apparent pore sizes in materials when studied with pulsed gradient-field spin-echo NMR.
We note that the sum in Eq. 18 converges rapidly. Terms with p > 0 are only important for Δt open /τ → 0 and even for, e.g., Δt open /τ 0.1, the relative difference between the exact expression and a truncation of it, which includes the first term only, is less than 0.5% (see Supplementary Figure S1 In the opposite limit, Δt open /τ ≪ 1, we Taylor expand the terms in the curly brackets in Eq. 18 and find Note that for Δt open /τ → 0, the variance is identical to the variance L 2 /12 for a particle in a box with side length L and uniform probability distribution, as it should be. These limiting behaviors are all indicated in Figure 2D. Since both parameters of interest, D and L, appear in the expressions for the variance, it is not possible to determine both of them independently based on the variance alone. Access to the full distribution of the measured positions would, in principle, allow this (Figure 2A-C). The MSD, however, also grants access to this information, as we shall see below.

MSD of Positions Recorded With Motion Blur, for a Brownian Trajectory in a Box
Obviously, the MSD must be bounded for confined Brownian motion. Here, we consider this MSD for positions recorded with motion blur: since 〈(x msr, n ) 2 〉 〈(x msr, 0 ) 2 〉 var(x msr ). Here, var(x msr ) is given by Eq. 18 and the remaining average in Eq. 21 is again over all trajectories. The expression for the time-ordered autocorrelation function is calculated similarly to the variance above, and the result for the MSD is The details of the calculation may be found in the Supplementary Material. This expression is exact and agrees with simulated data for a Brownian particle confined to a 1D box (Figure 3). For nΔt/τ ≪ 1, an expansion to lowest order in Δt/τ and Δt open /τ shows that the expression for the MSD in Eq. 22 in this limit reduces to the MSD for free diffusion in Eq. 11. That is, the influence of the confining walls is negligible in this case.
For nπ 2 Δt/4τ ≫ 1, we immediately get from Eq. 22 that MSD(n)x2 var(x msr ), where var(x msr ) is found in Eq. 18. Additionally, still in this limit, using Eq. 20 valid for Δt open /τ ≪ 1, we find that Thus, motion blur, described to lowest order in Δt open /τ, reduces the asymptotic value of MSD(n) by a factor (1 − Δt open /τ) relatively to its value in the absence of motion blur (Figure 3). Figure 3 also illustrates that both D and L, in principle, can be determined from data, provided that Δt < τ and n is sufficiently large. If the first condition is not satisfied, information about D is lost, and the MSD essentially assumes its plateau value already for n 1. If, on the other hand, n is limited and Δt < τ, e.g., due to poor sampling of a particle's trajectory, information about the confining domain is lost, since the particle, typically, does not encounter the boundaries in this case.
In the absence of motion blur and localization errors-i.e., if positions are recorded instantaneously with infinite precision-the MSD can be calculated exactly for all values of n directly from Eq. 13 ( Figure 3) [22], Results for the MSD similar to our Eq. 22 were previously found by Destainville and Salomé [18] in order to improve analysis of data from [17]. Considering a fully open shutter (Δt Δt open ), they elegantly arrived at an expression for the MSD by considering the lowest-order expansion of the correlation between the true positions by adjusting the characteristic time-scales to ensure that their expression obeys the correct scaling properties for short and longtime-scales. In our notation, their Eq. 8 reads 2 This is similar, but not identical to Eq. 22, see Figure 3, likely because their parameters were chosen to produce the correct scaling properties from a single term, i.e., similar to our term with p 0 in Eq. 22. Moreover, for n → ∞, the MSD reaches a plateau value equal to twice the measured variance, and keeping only the first term in the sum for the variance in Eq. 22 gives a more accurate expression for the plateau value than Eq. 25 when Δt open /τ > 0.08 (see Supplementary Figure S1). For analysis of data, they suggest to fit the measured MSD to a simpler expression (their Eq. 2) followed by an adjustment of the fitted parameters to account for the finite exposure time, rather than to fit directly with the MSD of Eq. 25. It is, however, simpler and more accurate to estimate D and L by fitting the measured MSDs with our exact expression, Eq. 22, using D and L as fitting parameters.

2D Brownian Motion Confined to a Disc
Another relevant case to consider is that of a Brownian particle confined to a 2D circular disc with radius a L/2 and reflecting walls at the boundary. Without loss of generality, we let the disc be centered at the origin. In this case, the problem does not factorize into a component for each Cartesian coordinate. Instead, we write the particle's position in polar coordinates as r r(r, ϕ) with 0 ≤ r < a and 0 ≤ ϕ < 2π. If particles are distributed uniformly over the disc, the probability density is in polar coordinates p(r, ϕ) r/(πa 2 ). The mean position is at the center, 〈r〉 (0, 0), in Cartesian coordinates, and the variance is For the case of Brownian motion confined to a 2D disc, the solution to the diffusion Fokker-Planck equation for the conditional probability density for the particle's position r at time t, given it was r 0 at time t 0, is, with closed boundary conditions at r a [22], where τ a 2 /D and α ℓm > 0 is the mth positive root, J ′ ℓ (α ℓm ) 0, of the derivative of the Bessel function of the first kind of order ℓ, J ℓ . The zeros are arranged in ascending order of magnitude: 0 < α ℓ1 < α ℓ2 < . . ..

Variance of Positions Recorded With Motion Blur, for a Brownian Trajectory on a 2D Disc
For the 2D case, we define the measured position similarly to Eq. 1. The variance of the measured positions is calculated in a manner similar to the 1D case, and the result is For details of the calculation, see Supplementary Material. This exact expression agrees with simulated data (Figure 4A-C) of particles undergoing Brownian motion confined to a 2D disc and subject to motion blur ( Figure 4D).
In the limit where the exposure time, Δt open , is negligible compared to the characteristic time, τ, Δt open τ , (29) which is equal to the variance of a particle's positions under instantaneous (and noise-free) recording (see Supplementary  Material).
On the other hand, for Δt open ≫ τ, we get var(r msr )x 7 24 To arrive at this expression, we evaluated the relevant sum following [25], see Supplementary Material.
These limiting behaviors are also indicated in Figure 4D. Qualitatively, the implications of the variance of particle positions confined to a 2D circular domain and with positions recorded using finite exposure time are similar to the case of a particle confined to a box.

MSD of Positions Recorded With Motion Blur, for a Brownian Trajectory on a 2D Disc
This expression is exact for all values of Δt open /τ and n ≥ 1, and it agrees with simulation results for particles undergoing 2D Brownian motion confined to a disc ( Figure 4E-G). To our knowledge, the literature contains no similar explicit expression for the MSD under these circumstances, let alone an exact one. When fitting Eq. 28 to a mean-squared displacement calculated from data, the sum converges rapidly. So, we give the first five values of α 1m for completeness: α 11 1.8412, α 12 5.3314, α 13 8.5363, α 14 11.7060, and α 15 14.8636.
In the case of instantaneous and noise-free recording of positions, it has been shown that the MSD can be calculated from Eq. 27 [22,25]. The result is That is, Eq. 32 is the 2D version of Eq. 24, valid for instantaneous (and noise-free) sampling. We note that this expression is recovered from our result, Equation (31), in the limit of Δt open /τ → 0 (Supplementary Material). This expression also agrees with simulated data of a particle undergoing 2D Brownian motion confined to a disc, with instantaneous and noise-free recording of simulated positions ( Figure 4E-G).

3D Brownian Motion Confined to a 3D Sphere
We now consider a Brownian particle confined to a 3D sphere with radius a L/2. Similar to the 2D case confined to a circle, solutions to the diffusion equation factorize when expressed in spherical coordinates with the center of the sphere chosen as origin. Then, r r(r, θ, ϕ) with 0 ≤ r < a, 0 ≤ θ < π, and 0 ≤ ϕ < 2π. For particles distributed uniformly within the sphere, the probability density distribution in Cartesian coordinates is p(x, y, z) 3/(4πa 3 ) for x 2 + y 2 + z 2 < a 2 . This becomes p(r, θ, ϕ) 3r 2 sinθ/(4πa 3 ) in spherical coordinates. So, the variance is var(r) 〈r 2 〉 The solution to the diffusion Fokker-Planck equation for the conditional probability density for the particle's position r at time t, given it was r 0 at time t 0, with closed boundary conditions at r a, is [22] P(r, r ′ , t) Here, τ a 2 /D and j ℓ (x) π/(2x) J ℓ+1/2 (x), where j ℓ (x) is a spherical Bessel function of the first kind and ℓ is a positive integer. The functions Y ℓm (θ, ϕ) are spherical harmonics, and β ℓn are the (non-zero) zeros of the derivatives of the spherical Bessel functions, j ′ ℓ (β ℓn ) 0, which are arranged in ascending order of magnitude 0 < β ℓ1 < β ℓ2 < . . ..

Variance of Positions Recorded With Motion Blur, for a Brownian Trajectory in a 3D Sphere
For the 3D case, we define once again the measured position similarly to Eq. 1. The calculation of the variance of the measured positions is similar to the 1D case leading to the result (for details, see Supplementary Material), Δt open τ .
Details of the calculation are presented in Supplementary Material.

MSD of Positions Recorded With Motion Blur, for a Brownian Trajectory in a 3D Sphere
A calculation similar to the one done for the 1D case gives the MSD for a 3D Brownian trajectory confined within a 3D sphere. The result is This expression is exact for all values of Δt open /τ and n ≥ 1. Figure 5E-G show excellent agreement between Eq. 37 and the MSD calculated from simulated Brownian trajectories confined within a 3D sphere and imaged with motion blur. This constitutes another result that, to our knowledge, does not exist in the literature. As for the 1D and 2D cases, the sum converges rapidly, which is useful knowledge when applying the expression in a fit to experimental data. So, for practical completeness and the reader's convenience, we again give the first five values of β 1n : β 11 2.08158, β 12 5.94037, β 13 9.20584, β 14 12.4044, and β 15 15.5792.
In the absence of motion-blur-i.e., in the case of instantaneous recording of exact positions-the MSD can be calculated directly from the propagator in Eq. 34 with the result [22] This result is also indicated in Figure 5E-G.

2D Anisotropic Brownian Motion
Confined to a 2D Disc particle with D y > D x relatively to those of a particle with D y D x ( Figure 6E). This effect flips as Δt/τ increases and motion blur in conjunction with D y > D x causes the observed MSDs to fall below the theory for isotropic Brownian motion ( Figure 6F,G). This illustrates that for anisotropic Brownian motion, the coupling between confinement and motion blur is non-trivial and care must be taken to rule out anisotropic diffusion before our exact results for isotropic Brownian motion are used to interpret experimental results.

DISCUSSION
In this work, we presented exact formulas for the variance and MSD of positions of a particle undergoing Brownian motion in confining spaces. These formulas remain valid when random Gaussian localization errors make true and recorded positions differ, provided that these localization errors are independently and identically distributed, which often is the case or may be assumed. To account for the localization error, our formulas should be offset merely by 2σ 2 loc for each coordinate contributing to the quantity of interest, the variance, or the MSD (Section 4.5). This localization error may be found using established methods [3,6].
We assumed normal Brownian motion throughout this article, and our formulas are guaranteed to apply only for such motion. On the other hand, we used only second moments of distributions, so results may in principle generalize to other cases with finite second moment if they exist. In some cases, particles may exhibit more complicated, types of motion, such as sub-diffusion due to interactions with other molecules in a membrane or "hop-diffusion": particles may be almost confined to nano-domains but "hop" between neighboring domains [17,27].
Other complication we ignored. Particles experience increased viscous drag near surfaces, as described by Faxén [28], and may experience attractive depletion forces near surfaces if many smaller particles-e.g., macromolecules-are present. Results of simulations and experiments, in which these and other effects are, respectively, simulated or suspected, can now be compared to the exact statistics given above for cases of steric confinement alone.
Additionally, while the approach described here delivers results for individual molecules, precision on estimates for parameters requires long trajectories [4]. In applications, the number of data points in a trajectory depends, e.g., on the fluorescent label and may be much lower than that used in our simulations. In such cases, one may compromise as regards true single-molecule results and average over a sample of several or many molecules to improve statistics of results obtained with our formulas [4].
Our formulas should facilitate that as much information as possible can be extracted about the mobility of the particles and the confining domains. To do this, one should account for the fact that the estimated MSD-values are correlated, since they are obtained from the same trajectory. This makes optimal fitting to MSDs complicated [16]. Additionally, in the present case of confined Brownian motion, such correlations most likely also depend on the position of the particle relative to the boundaries. Our formulas could, in principle, also facilitate a more accurate use of the MSD as a basis for classifying the type of motion [27].
The covariance-based estimator (CVE) introduced and demonstrated in [6,7,29] is an attractive alternative; however, away from surfaces, it estimates the bulk value of D, does that even from short segments of the trajectory, and does this optimally, provided its SNR DΔt √ /σ loc > 1 [6,7]. So, our exact MSD is probably best used with experimental MSDvalues as a fairly simple way to gauge deviation from bulk behavior as function of time-scale of measurement.
In conclusion, there is still work to do in the field of analysis of single-molecule and single-particle trajectories. Trajectories of self-propelled particles, e.g., are better understood as persistent random motion. Tools for analysis of such trajectories with accounting for motion blur and/or effects of confinement would amount to quite an oeuvre even for the simplest model, the Ornstein-Uhlenbeck process, we predict, based on [30], though methods developed in [31,32] may help. Here, we have exhausted a small fraction of this field with exact formulas for the variances and MSDs for simple Brownian motion, valid when particle positions are recorded using finite exposure time, which is the typical case by far. These formulas should help researchers to get the most out of their data for motion of Brownian particles and the domains that confine them, e.g., for particles diffusing in cell membranes [4,10,17].

Simulations
Equation (3) describes the propagator P(x, t|x 0 , t 0 ) of the diffusion equation. It shows how to propagate free diffusion between fixed points in time, t j jΔt (j integer): given that the particle is at position x j at time t j , simply draw a random number from a Gaussian distribution with zero mean and variance σ 2 2DΔt and add this random number to x j . The sum is x j+1 . Thus, we simulate time-lapse sampled positions on a trajectory of Brownian motion by iterating where (η j ) j 1,2,... is the normalized Gaussian white noise, i.e., for all i, j, 〈η j 〉 0 and 〈η i η j 〉 δ ij . The displacements Δx j x j − x j−1 are consequently by construction independent and normally distributed, with zero mean and variance σ 2 2DΔt. Since independent, they are uncorrelated, 〈Δx i Δx j 〉 0 for i ≠ j; compare Eqs. 4-6. Dotted lines through colored symbols: Exact results for isotropic Brownian motion, D y /D x 1, excluding (Eq. 32) and including (Eq. 31) the effect of motion blur. The horizontal black dashed line marks a 2 , i.e., twice the variance for a uniform distribution on a circle.