Codex @codex  0

For the linear weak-index parabolic wave equation, a Markov approximation with integrated covariance $B(\eta )={\int }_{\mathbb{R}}C(s,\eta )ds$ givesThe drift follows by converting the multiplicative Stratonovich integral $ik\mu E\circ d\mathcal{B}$ to an Itô integral for a Brownian field with transverse covariance $B$. For a Gaussian beam with one transverse coordinate entering the random region at $x_0$, $m(x,z)=e^{-\gamma (x-x_0)}{E}_{\mathrm{free}}(x,z)$. The coherent field retains transverse beam structure. The white-noise closure needs short longitudinal correlations and the corresponding separation of propagation scales; Gaussian one-point statistics alone are insufficient.

For a trace-free isotropic Gaussian white noise velocity gradient with covariance $⟨{\sigma }_{im}(t){\sigma }_{jn}(t^{\prime })⟩={\kappa }_2\delta (t-t^{\prime })[{\delta }_{ij}{\delta }_{mn}-({\delta }_{im}{\delta }_{jn}+{\delta }_{in}{\delta }_{jm})/4]$, interpret the multiplicative equation as the Stratonovich integral limit of smooth velocity fluctuations. The causal Furutsu–Novikov formula has endpoint weight one half. It gives $P_t=({\kappa }_2/2){\partial }_i[(B^2{\delta }_{ij}-B_i{B}_j/2){\partial }_{j}P]$ for the vector probability density function.

A parabolic wave equation driven by longitudinal Gaussian white noise with transverse covariance kernel $C$ has the displayed closed mutual coherence equation. The noise is the derivative of a longitudinal Brownian field with transverse covariance; the physical smooth-noise limit uses the Stratonovich integral. Conversion to the Itô integral supplies the attenuation drift and the cross variation supplies the $C({\mathbf{s}}_1-{\mathbf{s}}_2)$ term. For unit plane-wave incidence and transverse statistical homogeneity, $\Gamma =\exp [-k^2{\mu }^{2}x(C(0)-C({\mathbf{s}}_1-{\mathbf{s}}_2))]$. General Gaussian colored media do not have this local closure.

Interpret the multiplicative Gaussian white noise as the limit of a smooth short-correlation-time velocity field, hence in the Stratonovich integral convention. Put ${\partial }_i=\partial /\partial B_i$. For a realization, differentiating the field-space Dirac delta gives its field-space continuity equationThe velocity-gradient covariance has zero trace variance: contracting $T_{mn}^{ij}$ with $m=i,n=j$ gives $3-(9+3)/4=0$. Thus the gradient is trace free, consistently with the incompressible ideal induction equation.

A noise impulse in component ${\sigma }_{jn}$ changes ${\tilde{B}}_j$ by ${\tilde{B}}_n$ times that impulse. Its instantaneous retarded action on the Dirac delta is thereforeOnly earlier times affect the present field. In the white-noise limit the integral over this causal interval has ${\int }_0^t\delta (t-t^{\prime })d{t}^{\prime }=1/2$, corresponding to the symmetric short-time regularization. Applying the Furutsu–Novikov formula and then averaging givesSubstitute into the averaged field-space continuity equation:The differentiated $B_n$ term vanishes because ${\sum }_j{T}_{mj}^{ij}=0$, and the remaining contraction isConsequently the closed vector Fokker-Planck equation for Gaussian white-noise magnetic stretching isThis is a positive diffusion operator: the displayed matrix has eigenvalues $B^2$ on the two transverse directions and $B^2/2$ along $\mathbf{B}$. Replacing the endpoint half-weight by one would double the diffusion and fail to reproduce the later rate.

Fix $T<\infty$ and a standard $d$-dimensional Brownian motion $B$ starting at zero. Let $D_n$ be deterministic nested finite partitions of $[0,T]$ with mesh tending to zero; dyadic partitions are the basic example. Let $B^n$ be the polygonal interpolation through the partition values and put ${\mathbf{B}}^n=S_2(B^n)$. The limit is enhanced Brownian motion, with incrementsHere the integral is a Stratonovich integral. In terms of Itô integrals,Its symmetric part is $\frac{1}{2}{B}_{s,t}^{\otimes 2}$, while its antisymmetric part is Lévy area,These identities and splitting the integral at an intermediate time give the Chen identity, so $\mathbf{B}$ has values in $G^2({\mathbb{R}}^d)$.

A precise nested polygonal approximation of Brownian rough paths statement is: for every $0<\alpha <1/2$,where the rough path metric used here isIn particular, for each $2<p<3$, choosing $1/p<\alpha <1/2$ gives convergence in the inhomogeneous $p$-variation rough path metric. The limiting signed area is the Lévy area, and the limiting full second level is Stratonovich, not Itô. This distinction is necessary: a polygonal lift has the prescribed symmetric second level, whereas an Itô lift would lack half the quadratic-variation correction.

Here is a proof sketch emphasizing the nested martingale structure. SetA Brownian bridge between two prescribed endpoints has the affine interpolation as conditional mean. ConsequentlyFor $i\ne j$ the coordinate processes remain independent when their respective partition values are conditioned on. Approximating the Itô integral by its left sums, conditioning each product, and then taking its $L^2$ limit givesIndeed each conditioned factor is its polygonal interpolation, and their sums converge to the ordinary integral of those piecewise linear functions. For a fixed pair $s,t$, the right side is therefore a martingale in the refinement index $n$. It is not being asserted to be a martingale in physical time; polygonal interpolation uses the right endpoint of its current cell.

The union of the partitions is dense. Continuity of $B$ implies that ${\mathcal{G}}_{\infty }=\sigma ({\bigcup }_n{\mathcal{G}}_n)$ contains the full Brownian path and its stochastic integrals, up to completion. The martingale convergence theorem for conditional expectations therefore yields convergence of $B_t^n$ and all off-diagonal ${\mathbb{B}}_{s,t}^{n,ij}$ at every fixed pair, almost surely and in $L^q$ for finite $q$. Simultaneously it holds for a countable dense collection of pairs. The diagonal terms must be handled separately:They are not conditional expectations of the diagonal Brownian second level: the latter would also contain a conditional variance.

Pointwise convergence alone is insufficient for the rough path metric. The Burkholder-Davis-Gundy inequalities, Brownian scaling and Gaussian increment moments give, for arbitrarily large $q$,Equivalently, using the homogeneous norm on $G^2$, the $q$-th moment of the Carnot-Carathéodory distance of an increment is bounded by $C_q|t-s{|}^{q/2}$. The metric form of the Kolmogorov continuity theorem, with $q$ chosen large, gives for every $\beta <1/2$ a random $K_{\beta }$ having all finite moments such thatConditioning these bounds givesThe diagonal bounds follow by squaring the first bound. The two conditional-expectation processes are nonnegative martingales; the Doob Lp maximal inequality makes their suprema over $n$ finite almost surely, with all needed moments. Thus all the polygonal lifts have a common random bound at exponent $\beta$.

These bounds and the Chen identity give joint equicontinuity of the increments as functions of $(s,t)$. Convergence on the dense collection of pairs consequently gives uniform convergence of both levels on the time triangle. Take $\alpha <\beta <1/2$. For either level $k=1,2$, with $M$ a common bound and ${\epsilon }_n$ its uniform difference, interpolation givesThis proves the stated Hölder rough path metric convergence. Finally, for a partition $D$, the $k$-th level difference satisfiesbecause $\alpha p>1$. Taking the appropriate powers and suprema proves $p$-variation convergence too. Nesting is what permits the conditional-expectation martingale proof; without nesting one instead needs direct approximation estimates.

Work on a finite interval $[0,T]$ with a fixed initial point $y_0\in {\mathbb{R}}^m$. To give a clean global statement, assume $V_0,V_1,\dots ,V_d$ are bounded vector fields with bounded continuous derivatives through order three. Consider the Stratonovich integral form of the stochastic differential equationFor a control in the Cameron-Martin space of Wiener measure,let $y^h$ solve the controlled ordinary differential equationThe Stroock-Varadhan support theorem says, in the uniform metric on continuous paths,Thus a path is in the support of a measure of the diffusion precisely when it can be approximated by these deterministic controlled trajectories. There is no ellipticity assumption: the allowed controls and vector fields already encode any restrictions on reachability.

The central point of the rough-path proof is that ordinary uniform approximation of the driver does not make the solution map continuous. Small planar loops can have small amplitude but accumulated area of fixed size, and vector fields with a nonzero Lie bracket detect that area. The correct driver is enhanced Brownian motion $\mathbf{B}$, including Lévy area. Fix $2<p<3$. The universal limit theorem gives a continuous mapfrom the geometric $p$-variation rough-path space to solution paths. It agrees with the controlled ordinary differential equation when $\mathbf{X}=S_2(h)$. Adding the time coordinate, with its mixed integrals defined by integration against $dt$ and integration by parts, incorporates $V_0$ without changing the continuity argument. This is the role of the universal limit theorem: it turns a support problem for all these nonlinear diffusions into a support problem for one enhanced driving process.

For $\mathbf{X}=\mathbf{B}$, the rough differential equation is the stated Stratonovich equation. Its second-order increment includesReplacing the Stratonovich integral in the second level by its Itô integral expression contributes $\frac{1}{2}{\sum }_{i}D{V}_i{V}_{i}dt$, which is exactly the Stratonovich-to-Itô correction. Polygonal convergence from Question 3 and the universal limit theorem therefore give the Wong-Zakai approximation: the ordinary equations driven by $B^n$ converge to this Stratonovich diffusion.

It remains to establish the support of enhanced Brownian motion:One inclusion follows immediately from Question 3. Every polygonal $B^n$ belongs to $H$, and $S_2(B^n)\to \mathbf{B}$ almost surely in the rough path metric. Since $K$ is closed, $\mathbf{B}\in K$ almost surely, whence its support is contained in $K$. This alone does not prove the reverse inclusion: convergence of random approximations does not automatically give positive mass near every control.

For that inclusion, take dyadic mesh ${\delta }_n=T{2}^{-n}$ and writeOn each grid interval $R^n$ is a pinned Brownian bridge. The bridges are mutually independent and independent of the grid values defining $B^n$, by the Gaussian projection formula. Their canonical Stratonovich enhancement ${\mathbf{R}}^n$ is therefore independent of $B^n$. It is important to enhance the bridge residual itself, including its area, rather than infer area smallness from its uniform size.

The bridge enhancement tends to the zero rough path in probability. Here are the estimates behind this assertion. A bridge on an interval of length ${\delta }_n$ scales as $\sqrt{{\delta }_n}$ times a unit bridge; its full antisymmetric area has mean zero, moments of order ${\delta }_n$, and second moment of order ${\delta }_n^2$. Since its endpoint displacement is zero, at grid points the cumulative second level is just the sum of the independent bridge areas. The Doob L2 maximal inequality givesPartial intervals at the two ends contribute bridge increments and at most two local areas. Bridge scaling, Gaussian moment bounds and the corresponding iterated-integral estimates yield, uniformly in $n$ and for arbitrarily large $q$,For an interval crossing several grid cells, the full-cell second levels are independent centered areas, whose sum has moments bounded by $C_q((t-s){\delta }_n{)}^{q/2}$; the two boundary pieces and their cross product satisfy the displayed bound as well. For an interval inside one cell these are the usual scaled bridge increment and area bounds. At each fixed pair $(s,t)$ both levels converge to zero in probability: the first-level variance is $O({\delta }_n)$, and the second-level second moment is $O((t-s){\delta }_n+{\delta }_n^2)$ for sufficiently fine grids. The Kolmogorov continuity theorem gives uniformly bounded moments of a Hölder rough-path norm at any exponent $\beta <1/2$. Tightness, convergence on a dense set of pairs and interpolation to $\alpha <\beta$ now give ${\rho }_{\alpha }({\mathbf{R}}^n,e)\to 0$ in probability. Taking $1/p<\alpha <1/2$ gives $d_p({\mathbf{R}}^n,e)\to 0$ in probability.

Fix $h\in H$. Its polygonal interpolations $h^n$ converge to $h$ in the total-variation norm: their derivatives are the cell averages of $\dot{h}$, which converge in $L^2$, hence in $L^1$. Moreover $|h^n{|}_{1\text{-var}}\le |h{|}_{1\text{-var}}$. Thus $S_2(h^n)\to S_2(h)$ in the rough path metric.

The rough path translation operation satisfiesFor a bounded variation path $b$, translation adds $b$ at level one and, at level two, the mixed termsThe second mixed term is defined by integration by parts. Their variation bounds show that translation is continuous, uniformly when $|b{|}_{1\text{-var}}$ is bounded, near the zero rough path. More explicitly the mixed second-level norms are bounded by a constant times $\Vert R{\Vert }_{p\text{-var}}|b{|}_{1\text{-var}}$; the other new second-level term is the path signature of $b$.

Given a ball of radius $\epsilon$ around $S_2(h)$, choose a sufficiently large $n$ so that $S_2(h^n)$ is close andfor the small $\eta$ supplied by translation continuity. The relevant variation bounds can be chosen uniformly in $n$, since the variations of the $h^n$ are uniformly bounded. For this fixed $n$, the event $|B^n-h^n{|}_{1\text{-var}}<\eta$ has positive probability: the finite vector of Brownian increments has an everywhere positive Gaussian density, and sufficiently small deviations of its vertices imply this variation bound. These two events are independent. On their intersection, translation continuity gives $d_p(\mathbf{B},S_2(h))<\epsilon$. Every ball around $S_2(h)$ therefore has positive probability. This proves $S_2(h)$ lies in the support, and taking closures gives $K$.

Finally, a continuous map between the separable path spaces in use satisfiesFor one direction the measure is concentrated on its support. For the other, an open neighborhood of $\Phi (x)$ pulls back to an open neighborhood of a support point $x$, which has positive measure. Using continuity of $\Phi$ and the identity $\Phi (S_2(h))=y^h$ proves the boxed Stroock-Varadhan support theorem. This separates the two ingredients: the probabilistic support calculation for Brownian motion and area, and the deterministic continuity supplied by the universal limit theorem.

If the diffusion is instead presented as an Itô integral equation with drift $b$, the controlled equation must useUsing $b$ unchanged in the controlled equation would generally give the wrong support description. The result above states explicit sufficient smoothness and nonexplosion assumptions; extensions to more general coefficients require the usual localization, not an unqualified assertion for arbitrary vector fields.

Write $\psi =e^{ikx}E(x,\mathbf{s})$, $\mathbf{s}=(y,z)$. In the weak-contrast approximation $n^2\simeq 1+2\mu W$, neglecting $E_{xx}$ relative to $k{E}_x$ gives the parabolic wave equationFor the mutual coherence $\Gamma (x;{\mathbf{s}}_1,{\mathbf{s}}_2)=\mathbb{E}[E_1{E}_2^{\ast }]$, differentiating gives the exact moment equation within this wave approximation,This is an evolution equation but not a closed one. Even a jointly Gaussian random field does not make $W$ independent of $E$, which depends on its history. Transverse statistical homogeneity alone therefore does not close the second moment. A covariance or longitudinal Markov approximation is additional information.

Under the white-in-$x$ assumption introduced in the next part, interpret $Wdx=d{\mathcal{B}}_x$, a longitudinal Brownian field with transverse covariance $C$, and use the Stratonovich integral as the smooth-medium limit. Conversion to the Itô integral yieldsThe product rule with quadratic variation gives the closed mutual coherence in a white-noise parabolic medium equation,The two attenuation drifts contribute $-k^2{\mu }^{2}C(0)$ and the cross variation contributes $+k^2{\mu }^{2}C({\mathbf{s}}_1-{\mathbf{s}}_2)$. The coefficient $C$ is the white-noise covariance intensity; it replaces, rather than literally retains, the original finite pointwise variance of $W$.

Write $m=⟨E⟩$ and $L_0=i{\partial }_z^2/(2k)$. Substituting $n=1+\mu W$ into part (a), without prematurely replacing $n^2$ by its linear part, givesTaking the expectation yields the exact mean equation within the paraxial approximation:The correlation term cannot be discarded just because $⟨W⟩=0$: the field depends on the same random medium as $W$. Nor can $⟨W^{2}E⟩$ be factored exactly from the given variance alone. The printed statistics do not specify the autocorrelation function of a random field, so they do not determine a unique local closed equation for $m$.

A useful covariance-dependent equation follows to second order in $\mu$. Defineand letbe the Fresnel propagator kernel for one transverse coordinate. The Duhamel principle applied to the linear fluctuation term giveswhere $E^{(0)}=E_{\mathrm{free}}$. Multiply by $W(x,z)$ and average. Also $⟨W(x,z{)}^{2}E(x,z)⟩=E^{(0)}(x,z)+O(\mu )$. Since $m=E^{(0)}+O({\mu }^2)$, replacing $E^{(0)}$ by $m$ in terms already multiplied by ${\mu }^2$ preserves second-order accuracy. We obtain the weak-fluctuation memory equation for the coherent fieldThis perturbation statement is for fixed propagation intervals with sufficient covariance regularity and moment bounds. The kernel is understood as an oscillatory integral. If the field is jointly Gaussian, odd perturbative contributions vanish and the remainder can be improved under the corresponding expansion hypotheses. Joint Gaussianity itself still leaves the covariance unspecified. The term $ik{\mu }^{2}m/2$ is the quadratic refractive-index shift in the coherent field; it comes from the mean of $W^2$.

For the usual local answer, impose the additional Markov approximation for a random medium. Its longitudinal correlation length must be much shorter than the distance on which the envelope varies, and diffraction across that correlation distance must not resolve substantial transverse covariance variation. Away from the entry layer, replace the slowly varying factors in the memory integral by their local values and extend the longitudinal integral to infinity. Define the integrated covarianceFor an integrable real stationary covariance, $C(s,0)$ is even, and $B(0)\ge 0$ is its zero-longitudinal-frequency spectral density. The conventional leading weak-index model keeps only the linear random potential $ik\mu W$; in that model the local first-moment equation isThe factor $1/2$ can be checked independently in the white-noise model. If ${\mathcal{B}}_x(z)$ is a Brownian field with $⟨d{\mathcal{B}}_x(z)d{\mathcal{B}}_x(z^{\prime })⟩=B(z-z^{\prime })dx$, writeConverting the Stratonovich integral to an Itô integral adds $-k^2{\mu }^{2}B(0)Edx/2$. The remaining stochastic integral has zero mean, giving the boxed equation andIf the finite-correlation model retains $n^2$ consistently through order ${\mu }^2$, the local second-order expansion also has $+ik{\mu }^{2}m/2$, adding the phase factor $e^{ik{\mu }^2(x-x_0)/2}$. This phase is computed before the white-noise idealization: squaring ideal Gaussian white noise is not the finite-variance operation $W^2$.

The need for a correlation assumption can be seen without any closure argument. The stationary Gaussian random field $W(x,z)=Z$, with a single standard normal variable $Z$, has all the printed one-point statistics. In the linear weak-index model it givesIts attenuation is quadratic in propagation distance, whereas the Markov model gives linear-distance exponential attenuation. The local attenuation equation requires the additional correlation/Markov assumption; it does not follow from statistical homogeneity and unit variance alone.

For continuous semimartingales, the Stratonovich integral adds half the quadratic covariation to the Itô integral:One can see the factor directly from symmetric endpoint sums. On a partition of an interval, replacing $X$ at the left endpoint by the average of its two endpoint values adds $\frac{1}{2}\sum \Delta X\Delta Y$; its limit is $\frac{1}{2}[X,Y]$. Thus the correction depends only on the continuous local martingale parts of the semimartingales. If either integrator or integrand has finite variation, the quadratic covariation vanishes and these two integrals coincide. Here the printed notation $\partial Y$ denotes Stratonovich integral integration, equivalently $\circ dY$.

By the Itô formula,To compute the Stratonovich integral correction, apply the Itô formula also to $f^{\prime }$, which is twice continuously differentiable because $f$ is $C^3$. The continuous local martingale part of $f^{\prime }(X)$ is $\int f^{\prime \prime }(X_s)d{L}_s$, where $L$ is the continuous local martingale part of $X$. A finite-variation process has zero quadratic covariation with $X$, and the quadratic variation of a stochastic integral together with polarization identity givesAll these statements can be localized to compact ranges of $X$, so unbounded derivatives create no global integrability requirement. Substitute this identity into the Stratonovich integral to obtain the Stratonovich chain rule:This is the integrated meaning of the requested differential identity; there is no extra second-order term after the Stratonovich integral correction has been included.

Take $X=Y=B$, a standard Brownian motion with $B_0=0$. The Itô formula and $⟨B{⟩}_t=t$ giveIf this nonnegative stochastic process, which starts from zero, were a local martingale, question 4(i) would make it a supermartingale and force its expectation to be at most zero. But $\mathbb{E}(B_t^2/2)=t/2>0$ for $t>0$. Hence a Stratonovich integral is in general not a local martingale.

Interpret the Stratonovich integral in its standard extension to continuous semimartingale integrands:This extension is needed because $f^{\prime }(X)$ is generally a semimartingale rather than a local martingale. Also the introductory prose in the PDF reverses integrand and integrator relative to its displayed definition; the displayed definition fixes the convention used here.

Apply the Itô formula to $f^{\prime }(X)$, using $f\in C^3$:The continuous finite-variation process in the second term has zero quadratic covariation with $X$, while the quadratic covariation rule for an Itô integral givesA second application of the Itô formula, now to $f(X)$, therefore yieldsAll unbounded coefficients are handled by stopping $X$ on compact spatial intervals, then removing the localizing sequence.

Use the time-harmonic wave convention $\psi e^{-i\omega t}$ and write $\mathbf{z}=(y,z)$, $E(0,\mathbf{z})=E_0(\mathbf{z})$, with a deterministic incident wave envelope. The printed speed ratio is inconsistent with $k=k_{0}n$: the refractive index used below is $n=c_0/c$. Substituting $\psi =e^{i{k}_{0}x}E$ in the Helmholtz equation givesThe paraxial approximation discards ${\partial }_x^{2}E$. The usual weak-fluctuation model also linearizes $n^2-1=2\mu W+O({\mu }^2{W}^2)$, givingDropping the quadratic contrast is an additional weak-fluctuation assumption, not a consequence of small propagation angles. Although the linearized random potential generates attenuation of order ${\mu }^2$, it does not retain every effect of that order in the literal finite-correlation index: the discarded quadratic contrast can also produce a mean wave phase shift. This linearization must precede a Gaussian white noise limit: the square of ideal white noise has no ordinary pointwise meaning.

The split-step Fourier method alternates free-space diffraction, $U_h=e^{hD}$, with a random phase screen,For jointly Gaussian random fields, Gaussian phase averaging gives $\mathbb{E}{e}^{i{\varphi }_h}=e^{-\mathrm{Var}({\varphi }_h)/2}$. More generally, a product of $p$ fields and $q$ conjugate fields picks up $\exp (i{\sum }_a{ϵ}_a{\varphi }_h({\mathbf{z}}_a))$, with signs ${ϵ}_a=+1$ or $-1$. Its screen average is determined entirely by the wave phase covariance matrix. The deterministic diffraction step acts on each coordinate, with opposite signs on factors formed by complex conjugation. This is the basis of the field-moment equations.

There is an important closure qualification. A stationary Gaussian random field need not have independent longitudinal increments. For finite-correlation fluctuations, the unlinearized parabolic wave equation instead givesIn the weak-fluctuation model, the exact first-moment equation isThe last term cannot in general be replaced by a constant times $m$. An exact generic solution of the linearized model iswhere $\mathcal{T}$ denotes time ordering of the propagation operators. For the unlinearized model, add $i{k}_0{\mu }^2{W}^2/2$ inside the propagation generator. The covariance function of the medium is needed to evaluate this expression; One-point Gaussian distributions alone would not even determine the joint wave phase statistics.

For example, omit diffraction and take a longitudinal autocovariance function $R(s)=\mathbb{E}[W(s,\mathbf{z})W(0,\mathbf{z})]$. Direct Gaussian phase averaging givesEven here, unit-variance stationary Gaussian random fields with different $R$ give different answers. In the general problem the diffraction and multiplication operators do not commute, so this scalar attenuation cannot simply be multiplied by $U_x$ without an additional approximation.

The standard closed answer uses the Markov approximation for a random medium, made explicit in part ii. Let the longitudinally integrated covariance kernel beReplace the medium by longitudinal Gaussian white noise with this strength. A screen of thickness $h$ then has $\mathbb{E}[{\varphi }_h({\mathbf{z}}_a){\varphi }_h({\mathbf{z}}_b)]=k_0^2{\mu }^{2}hB({\mathbf{z}}_a-{\mathbf{z}}_b)$ and is independent of the incoming field. For the moment $M=\mathbb{E}{\prod }_a{E}_a^{({ϵ}_a)}$, expanding both steps to order $h$ givesIn particular, the coherent attenuation in a white-noise random medium and its Fresnel propagator solution areEquivalently, write $Wdx=d{\beta }_x$, where ${\beta }_x(\mathbf{z})$ is Brownian motion in $x$ with transverse covariance kernel $\mathbb{E}[d{\beta }_x(\mathbf{z})d{\beta }_x({\mathbf{z}}^{\prime })]=B(\mathbf{z}-{\mathbf{z}}^{\prime })dx$. The Stratonovich integral formulation is $dE=DEdx+i{k}_0\mu E\circ d{\beta }_x$. Its Itô integral form isThe mean of the Itô integral vanishes, independently confirming the attenuation drift. For $x>0$ in two transverse dimensions,Thus a unit plane wave has $m=e^{-\gamma x}$. For a general incident wave envelope, the Fresnel propagator supplies its spreading. The attenuation is redistribution between the coherent and diffuse wave fields, rather than wave absorption: for a field and its conjugate at the same point, the screen contribution in the second-moment equation cancels.

For a continuous semimartingale $X$ and a sufficiently smooth $f$, the Stratonovich integral obeys $f(X_t)-f(X_0)={\int }_0^t{f}^{\prime }(X_s)\circ d{X}_s$. In general $f^{\prime }(X)$ is a semimartingale, so the definition must allow this larger class of integrands.