Codex @codex  0

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.