← Stochastic Differential Equations

Brownian motion

Brownian motion is the basic random process used to model continuously evolving stochastic fluctuations in many stochastic differential equations.

Core idea. Brownian motion \(W(t)\) changes continuously through time, but its future movement is random. Over a short interval, it can move up or down, and the size of that random movement is governed by the length of the interval.

Start with the simplest question

Suppose we know the current value \(W(t)\). What happens during the next short time interval \(\Delta t\)?

The change is

\[\boxed{\Delta W=W(t+\Delta t)-W(t)}.\]

For Brownian motion,

\[\boxed{\Delta W\sim N(0,\Delta t)}.\]

This one statement contains the central rule.

PartMeaning
mean \(0\)the increment has no preferred upward or downward direction
variance \(\Delta t\)the spread of the random increment grows with the time interval
standard deviation \(\sqrt{\Delta t}\)the typical size of a Brownian increment is of order \(\sqrt{\Delta t}\)

Why do we write \(\Delta W=\sqrt{\Delta t}\,Z\)?

If

\[Z\sim N(0,1),\]

then multiplying by \(\sqrt{\Delta t}\) gives

\[\boxed{\Delta W=\sqrt{\Delta t}\,Z}.\]

Its mean remains zero:

\[E[\Delta W]=0,\]

and its variance becomes

\[\operatorname{Var}(\Delta W)=\Delta t.\]
This is why the square root appears. Variance scales with the square of the multiplier. To obtain variance \(\Delta t\), the normal random number must be multiplied by \(\sqrt{\Delta t}\).

A numerical example

Suppose

\[\Delta t=0.04.\]

Then

\[\sqrt{\Delta t}=0.2.\]

If a standard normal draw gives

\[Z=1.3,\]

then

\[\Delta W=0.2(1.3)=0.26.\]

If instead \(Z=-0.8\), then

\[\Delta W=0.2(-0.8)=-0.16.\]

So the process may move either upward or downward.

What does a Brownian path look like?

Three numerical Brownian paths generated using independent increments \(\Delta W=\sqrt{\Delta t}Z\) with a small \(\Delta t\). They start from the same point but quickly separate because their random increments differ.

Every path starts at

\[W(0)=0.\]

After that, each path follows its own random sequence of increments.

The plotted curves are numerical approximations. A computer can only simulate Brownian motion at finitely many time points. The mathematical Brownian path is defined continuously for every \(t\ge0\).

Brownian motion is continuous, not discrete

A simulation may use points such as

\[0,\Delta t,2\Delta t,3\Delta t,\ldots\]

to approximate the process. That does not make Brownian motion a discrete-time process.

The true mathematical process is indexed by continuous time:

\[\boxed{W(t),\qquad t\ge0}.\]

The discrete grid is only a computational approximation used to draw or simulate the continuous-time process.

The formal properties

Standard Brownian motion satisfies four key properties.

PropertyMeaning
\(W(0)=0\)the process starts at zero
independent incrementschanges over non-overlapping time intervals are independent
normal increments\(W(t)-W(s)\sim N(0,t-s)\)
continuous pathsthe path has no jumps

Independent increments

Consider two non-overlapping intervals:

\[[t_1,t_2]\qquad\text{and}\qquad[t_3,t_4],\qquad t_2\le t_3.\]

The increments

\[W(t_2)-W(t_1)\]

and

\[W(t_4)-W(t_3)\]

are independent.

Knowing what Brownian motion did in one interval does not tell us the random increment in a separate future interval.

The distribution depends only on interval length

For \(0\le s\[\boxed{W(t)-W(s)\sim N(0,t-s)}.\]

The distribution depends on the interval length \(t-s\), not on the absolute starting time \(s\).

Example. The Brownian increment from time 1 to time 2 has the same distribution as the increment from time 10 to time 11: both are \(N(0,1)\).

Brownian motion at a fixed time

Because \(W(0)=0\), taking \(s=0\) gives

\[\boxed{W(t)\sim N(0,t)}.\]

Therefore

\[E[W(t)]=0,\qquad\operatorname{Var}(W(t))=t.\]

Its standard deviation is

\[\sqrt{t}.\]

So the possible Brownian values spread out as time increases.

The distribution at different times

Because \(W(t)\sim N(0,t)\), the distribution remains centred at zero while its spread increases with time.

At small \(t\), the distribution is narrow. At larger \(t\), it is wider because the variance equals \(t\).

Why Brownian motion does not simply average back to zero

The expected value is zero:

\[E[W(t)]=0.\]

But this does not mean that one Brownian path stays near zero. It means that if we average many independent paths at the same time \(t\), the positive and negative values balance on average.

Expected value and individual path are different ideas. A particular trajectory can wander far from zero even though the ensemble mean remains zero.

The size of a Brownian increment

For a small interval \(\Delta t\),

\[\operatorname{SD}(\Delta W)=\sqrt{\Delta t}.\]

For example:

\(\Delta t\)standard deviation of \(\Delta W\)
11
0.010.1
0.00010.01

As the time interval shrinks, Brownian increments shrink—but only like the square root of time.

Why Brownian paths are so rough

An ordinary differentiable function changes approximately like

\[\Delta x\propto\Delta t.\]

Brownian motion changes on the larger scale

\[\Delta W\propto\sqrt{\Delta t}.\]

Therefore

\[\frac{\Delta W}{\Delta t}\sim\frac{1}{\sqrt{\Delta t}},\]

whose typical magnitude grows without bound as \(\Delta t\to0\).

This explains the roughness. Brownian motion is continuous, but almost surely nowhere differentiable. It has no ordinary instantaneous velocity \(dW/dt\).

Continuous but nowhere differentiable

These two statements are not contradictory:

ContinuousNot differentiable
there are no jumps in the paththe path is too irregular to have a finite ordinary slope
nearby times have nearby valueszooming in does not make the path look smooth

This unusual combination is exactly why ordinary calculus must be modified when Brownian motion appears in differential equations.

Brownian motion and ordinary time behave differently

In ordinary calculus, a small time increment is \(dt\). Brownian motion has a random increment \(dW_t\).

Their characteristic sizes are different:

\[dt\sim dt,\qquad dW_t\sim\sqrt{dt}.\]

This means Brownian fluctuations are much larger than \(dt\) on very small time scales.

Why does \((dW_t)^2\) behave like \(dt\)?

Since a Brownian increment has typical size \(\sqrt{dt}\), squaring gives

\[(dW_t)^2\sim dt.\]

In Itô calculus this becomes the formal rule

\[\boxed{(dW_t)^2=dt}.\]

This does not mean that ordinary algebra has suddenly made a random number exactly equal to time. It is shorthand for a precise limiting property called quadratic variation.

Do not treat \(dW_t\) as an ordinary differential. Rules involving \(dW_t\) belong to stochastic calculus.

Quadratic variation intuition

Divide \([0,T]\) into many short intervals. Ordinary smooth increments satisfy squared changes that vanish when summed over a very fine partition. Brownian increments behave differently:

\[\boxed{\sum_k(\Delta W_k)^2\longrightarrow T}.\]

This non-zero quadratic variation is the mathematical reason that second-order terms survive in Itô's formula.

Brownian motion as accumulated random shocks

Brownian motion can be viewed informally as the continuous-time limit of many tiny independent random shocks.

many tiny independent fluctuations→accumulate through time→continuous random path→Brownian motion

This viewpoint explains why Brownian motion is useful for representing unresolved environmental or demographic fluctuations when a continuous approximation is appropriate.

Connection to stochastic differential equations

A stochastic differential equation often has the form

\[\boxed{dX_t=f(X_t,t)\,dt+g(X_t,t)\,dW_t}.\]
TermRole
\(f(X_t,t)dt\)deterministic drift
\(g(X_t,t)dW_t\)random fluctuation
\(W_t\)Brownian motion driving the noise

The Brownian term does not mean that the biological variable itself is Brownian motion. It means Brownian increments are being used to drive random fluctuations in that variable.

A simple SDE step

For

\[dX_t=\mu\,dt+\sigma\,dW_t,\]

over a small interval \(\Delta t\),

\[\Delta X\approx\mu\Delta t+\sigma\Delta W.\]

Using

\[\Delta W=\sqrt{\Delta t}Z,\qquad Z\sim N(0,1),\]

gives

\[\boxed{\Delta X\approx\mu\Delta t+\sigma\sqrt{\Delta t}Z}.\]

This is the basic structure behind the Euler–Maruyama method developed later.

Brownian motion versus a CTMC path

CTMCBrownian motion
state often changes by discrete jumpspath is continuous
event times are individually identifiablefluctuation is represented continuously
useful for event-level population modelsuseful for diffusion and SDE approximations

A CTMC may be more natural when individual infections or recoveries matter explicitly. Brownian-driven SDEs can be useful when populations are large enough that a continuous random approximation is appropriate.

What Brownian motion does not mean

Brownian motion is a mathematical model of continuous stochastic fluctuation. It does not claim that every biological source of randomness is literally Gaussian, independent or memoryless. Whether Brownian noise is appropriate depends on the biological mechanism and the modelling scale.

Key idea. Brownian motion is a continuous-time stochastic process with \(W(0)=0\), independent Gaussian increments and \(W(t)-W(s)\sim N(0,t-s)\). Its increments scale like \(\sqrt{\Delta t}\), making paths continuous but extremely rough. This unusual scaling is the foundation of the \(dW_t\) term in stochastic differential equations and of Itô calculus.