Diffusion Coefficients as Vector Fields
In a multidimensional stochastic differential equation, each independent Brownian noise pushes the state in a particular direction, and that direction can vary with the state. The diffusion coefficient for one noise is therefore naturally viewed as a vector field. When the direction changes from point to point, the Stratonovich-to-Itô conversion adds a drift: each field is differentiated in its own direction. This additional term is called the noise-induced drift.
Work on a filtered probability space \((\Omega,\mathcal{F},(\mathcal{F}_t)_{t\ge0},\mathbb{P})\) carrying a standard \(m\)-dimensional Brownian motion \(\mathbf{B}_t=(B_t^1,\ldots,B_t^m)\) with independent components. Let the state \(X_t\) take values in \(\mathbb{R}^d\), and let \(V_k:\mathbb{R}^d\to\mathbb{R}^d\) be the diffusion vector field corresponding to \(B^k\). The notation \(DV_k(x)\) denotes the \(d\times d\) derivative matrix of \(V_k\) at \(x\).
For clarity, the \(i\)-th coordinate of the correction from one field is \((DV_k(x)V_k(x))_i=\sum_{j=1}^d \partial_j V_k^i(x)V_k^j(x)\). It is the directional derivative of the \(i\)-th component of \(V_k\) along \(V_k\) itself. With independent Brownian components, there is one such term for each field, and the terms are summed. There are no cross terms between distinct fields in this conversion.
All stochastic integrals below are Itô integrals in the class \(\mathcal{L}^2\): their integrands are progressively measurable and satisfy \(\mathbb{E}\int_0^T |H_s|^2\,ds<\infty\) for each finite \(T\). Stratonovich integrals are understood through the Itô integral plus one half of the relevant covariation. For locally bounded coefficients, the same conversion can be applied up to exit times from bounded regions; the displayed global examples have linear, globally Lipschitz coefficients. The standard existence-and-uniqueness result for globally Lipschitz coefficients of at most linear growth gives their solutions on finite time intervals.
Why the Directional Derivative Appears
Proof. By the definition of the Stratonovich integral, its \(i\)-th coordinate for field \(k\) is \(\int_0^t V_k^i(X_s)\,dB_s^k+\tfrac12[V_k^i(X),B^k]_t\). We compute the covariation. Write the Itô diffusion coefficient of \(X\) in coordinate \(j\) against \(B^\ell\) as \(V_\ell^j(X_s)\); a finite-variation drift has zero covariation with Brownian motion. Apply Itô's formula to \(g(x)=V_k^i(x)\). Its first derivatives are \(\partial_jg=\partial_jV_k^i\), and its second derivatives are \(\partial_{j r}g=\partial_{j r}V_k^i\). The formula includes the correction explicitly:
The second-derivative term and the finite-variation part of \(dX_t^j\) have zero covariation with \(B^k\). Since \(d[B^\ell,B^k]_t\) is \(dt\) when \(\ell=k\) and zero otherwise, the remaining term gives \(d[V_k^i(X),B^k]_t=\sum_j\partial_jV_k^i(X_t)V_k^j(X_t)\,dt\). Consequently the correction in coordinate \(i\) is half the sum over \(k\) of these directional derivatives. This is exactly \(\tfrac12\sum_k DV_k(X_t)V_k(X_t)\), proving the asserted drift. \(\square\)
This argument uses the Itô multiplication rules \(dB_t^k\,dB_t^\ell=\mathbf{1}_{\{k=\ell\}}dt\), \(dB_t^k\,dt=0\), and \(dt\,dt=0\). In particular, the correction is not found by differentiating the drift \(b\), nor by differentiating one field along a different independent noise field.
Worked Two-Dimensional Systems
Worked Example: A Triangular Field Produces a Constant Drift
Consider one Brownian motion and the field \(V(x,y)=(y,1)\), with zero Stratonovich drift and initial state \((0,0)\). Find the Itô equation and verify its solution.
The derivative matrix and its action on the field are \(DV(x,y)=\begin{pmatrix}0&1\\0&0\end{pmatrix}\) and \(DV(x,y)V(x,y)=(1,0)\). Therefore the Itô correction is \((\tfrac12,0)\), and the equations become \(dX_t=\tfrac12dt+Y_t\,dB_t\), \(dY_t=dB_t\).
The integrands \(Y_s=B_s\) are progressively measurable and satisfy \(\mathbb{E}\int_0^T B_s^2ds=\int_0^T s\,ds=T^2/2<\infty\). The second equation gives \(Y_t=B_t\). The first then gives \(X_t=\tfrac12t+\int_0^t B_s\,dB_s\). To evaluate the integral, use Itô's formula with \(f(x)=x^2\), \(f'(x)=2x\), \(f''(x)=2\), and the explicit correction \(\tfrac12 f''(B_t)dt=dt\): \(d(B_t^2)=2B_t\,dB_t+dt\). Since \(B_0=0\), this yields \(X_t=\tfrac12B_t^2\).
Check: The candidate \(X_t=\tfrac12B_t^2\) has differential \(dX_t=B_t\,dB_t+\tfrac12dt\), which agrees with the converted equation because \(Y_t=B_t\). Also, the original Stratonovich relation is \(dX_t=Y_t\circ dB_t\); with \(Y_t=B_t\), the Stratonovich chain rule gives \(X_t=\tfrac12B_t^2\), confirming the same result.
Worked Example: Rotational Noise and Radial Drift
Let \(V(x,y)=(-y/2,x/2)\), with zero Stratonovich drift and initial state \((1,0)\). Convert the system and check what happens to the radius.
Here \(DV=\begin{pmatrix}0&-1/2\\1/2&0\end{pmatrix}\), so \(DV(x,y)V(x,y)=(-x/4,-y/4)\). The Itô correction is \((-x/8,-y/8)\), giving \(dX_t=-Y_t\,dB_t/2-X_t\,dt/8\) and \(dY_t=X_t\,dB_t/2-Y_t\,dt/8\).
An explicit candidate solution is \(X_t=\cos(B_t/2)\), \(Y_t=\sin(B_t/2)\). For \(f(u)=\cos(u/2)\), \(f'(u)=-\tfrac12\sin(u/2)\) and \(f''(u)=-\tfrac14\cos(u/2)\); Itô's formula gives \(dX_t=-\tfrac12\sin(B_t/2)dB_t-\tfrac18\cos(B_t/2)dt\). For \(g(u)=\sin(u/2)\), \(g'(u)=\tfrac12\cos(u/2)\) and \(g''(u)=-\tfrac14\sin(u/2)\); it gives \(dY_t=\tfrac12\cos(B_t/2)dB_t-\tfrac18\sin(B_t/2)dt\). These are precisely the converted equations, and the initial values are \((1,0)\). The integrands are bounded, so they belong to \(\mathcal{L}^2\) on every finite interval.
Check: Directly, \(X_t^2+Y_t^2=1\). In the Itô equation the drift contribution to \(d(X_t^2+Y_t^2)\) is \(-\tfrac14(X_t^2+Y_t^2)dt+\tfrac14(X_t^2+Y_t^2)dt=0\); the first term comes from twice the radial drift and the second from the squared diffusion coefficients. The Brownian terms cancel. Thus the Itô correction balances the quadratic-variation contribution, as required for pure rotation in Stratonovich form.
Worked Example: Two Independent Coordinate Fields
Let \(B^1,B^2\) be independent Brownian motions and take \(V_1(x,y)=(x,0)\), \(V_2(x,y)=(0,y)\), with zero Stratonovich drift and initial state \((1,2)\). Compute the conversion and verify an explicit solution.
The derivative matrices are \(DV_1=\begin{pmatrix}1&0\\0&0\end{pmatrix}\) and \(DV_2=\begin{pmatrix}0&0\\0&1\end{pmatrix}\). Thus \(DV_1V_1=(x,0)\) and \(DV_2V_2=(0,y)\). The Itô equations are \(dX_t=\tfrac12X_tdt+X_t\,dB_t^1\) and \(dY_t=\tfrac12Y_tdt+Y_t\,dB_t^2\).
The candidate solution is \(X_t=e^{B_t^1}\), \(Y_t=2e^{B_t^2}\). For \(f(u)=e^u\), \(f'(u)=e^u\) and \(f''(u)=e^u\), so Itô's formula, including its \(\tfrac12 f''dt\) term, gives \(d(e^{B_t^1})=e^{B_t^1}dB_t^1+\tfrac12e^{B_t^1}dt\). Applying the same calculation to \(2e^{B_t^2}\) gives the second equation. The initial values are \(1\) and \(2\). Each integrand has finite expected squared integral: for example, \(\mathbb{E}\int_0^T e^{2B_s^1}ds=\int_0^T e^{2s}ds=(e^{2T}-1)/2<\infty\), and the corresponding expectation for \(4e^{2B_s^2}\) is \(4\int_0^T e^{2s}ds<\infty\).
Check: The two correction terms affect separate coordinates, because each field points along its own coordinate axis. In particular, there is no term involving a derivative of \(V_1\) along \(V_2\) or of \(V_2\) along \(V_1\). The verified solutions also have the expected initial values \(e^0=1\) and \(2e^0=2\).
The Generator View
The same correction can be recognized by comparing the differential operators associated with the two descriptions. For a twice continuously differentiable test function \(f:\mathbb{R}^d\to\mathbb{R}\), define \(V_k f(x)=\nabla f(x)\cdot V_k(x)\). In Stratonovich form the generator is \(b\cdot\nabla+\tfrac12\sum_k V_k(V_k f)\), where the field acts as a directional derivative. The corresponding Itô generator uses the corrected drift.
Proof. For each \(k\), differentiate \(V_k f=\sum_i V_k^i\partial_i f\) along \(V_k\): \(V_k(V_k f)=\sum_{j,i}V_k^j(\partial_jV_k^i)\partial_i f+\sum_{j,i}V_k^jV_k^i\partial_{ji}f\). The first sum is \((DV_kV_k)\cdot\nabla f\), and the second is \(V_k^\mathsf{T}(D^2f)V_k\). Summing these identities, multiplying by one half, and adding \(b\cdot\nabla f\) gives the displayed equality. \(\square\)
This equality explains why the correction matters beyond rewriting an equation: it ensures that the Stratonovich and Itô descriptions give the same infinitesimal action on smooth test functions. It also provides a useful independent check of a conversion: the first-order drift term must contain exactly the directional derivatives of the diffusion fields.
A Common Mistake
A frequent error is to differentiate each coordinate of a field with respect to the coordinate bearing the same index, rather than calculate the full directional derivative. For \(V(x,y)=(y,1)\), that shortcut might incorrectly report \(DV\,V=(0,0)\), and hence no drift. But the derivative matrix is \(\begin{pmatrix}0&1\\0&0\end{pmatrix}\), and multiplying it by \((y,1)\) gives \((1,0)\). The correct correction is therefore \((1/2,0)\). The worked solution \(X_t=B_t^2/2\) verifies the correction: its Itô differential contains \(\tfrac12dt\), which the incorrect conversion would omit.
The distinction is between differentiating a field and differentiating it in the direction of that field. In coordinates, every partial derivative that contributes must be multiplied by the corresponding component of the same diffusion field. When there are several independent Brownian drivers, calculate this product field by field and then sum.
| Stratonovich diffusion field | Directional derivative | Itô drift correction |
|---|---|---|
| \(V(x,y)=(y,1)\) | \(DV\,V=(1,0)\) | \((1/2,0)\) |
| \(V(x,y)=(-y/2,x/2)\) | \(DV\,V=(-x/4,-y/4)\) | \((-x/8,-y/8)\) |
| \(V_1=(x,0),\ V_2=(0,y)\) | \(DV_1V_1=(x,0),\ DV_2V_2=(0,y)\) | \((x/2,y/2)\) |
Check Your Understanding
For each question, identify the diffusion field or fields before calculating the correction.
- For \(V(x,y)=(x+y,x)\), compute \(DV(x,y)V(x,y)\) and state the Itô drift correction.
- For \(V_1(x,y)=(y,0)\) and \(V_2(x,y)=(0,x)\) driven by independent Brownian motions, calculate each field's own directional derivative and their summed correction.
- Explain why \(DV_1V_2\) is not a term in the conversion formula for independent Brownian drivers.
- For the rotational field \(V(x,y)=(-y/2,x/2)\), verify algebraically that the Itô drift contribution to \(d(X_t^2+Y_t^2)\) cancels the quadratic-variation contribution.
- For \(V(x,y)=(y,1)\), use Itô's formula for \(B_t^2\) to check the solution \(X_t=B_t^2/2\) against the converted equation.