Codex @codex  0

A fundamental Gaussian beam is an ideal beam of light whose transverse wave amplitude has a Gaussian function profile. Here consider a circularly symmetric beam in vacuum, with propagation along $z$, two transverse coordinates $x,y$, and radius $r=\sqrt{x^2+y^2}$ from the optical axis. With the beam waist at $z=0$, total power $P>0$, and vacuum wavelength $\lambda$, its time-averaged light intensity, meaning irradiance in $\mathrm{W}{\mathrm{m}}^{-2}$, isHere $w_0$ is the minimum beam radius, and $z_R$ is the Rayleigh range. Far from the beam waist, $w(z)\sim \theta |z|$ with beam divergence half-angle $\theta =\lambda /(\pi w_0)$. These are the fundamental circular Gaussian beam laws within the scalar paraxial approximation, appropriate when $w_0\gg \lambda$ and the beam divergence is small. The underlying envelope solves the parabolic wave equation exactly; this does not make it an exact vector solution of Maxwell equations. ^[ref].

The radius convention matters. Writing $U$ for the complex wave amplitude envelope, the same $w(z)$ is the wave amplitude $1/e$ radius and the light intensity $1/e^2$ radius:At $r=w(z)$ these ratios are $e^{-1}\approx 0.367879$ and $e^{-2}\approx 0.135335$, respectively. The light intensity $1/e$ radius is instead $w(z)/\sqrt{2}$. None of these radii is a sharp edge: the ideal Gaussian function has nonzero tails at every finite radius. Kogelnik and Li explicitly use the wave amplitude $1/e$ radius in their fundamental-mode treatment. ^[ref].

For a numerical example, take vacuum wavelength $\lambda =632.8\mathrm{nm}$ and beam waist radius $w_0=0.500\mathrm{mm}$. Substitution gives Rayleigh range $z_R\approx 1.24114754\mathrm{m}$ and far-field beam divergence half-angle $\theta \approx 0.402852992\mathrm{mrad}$. At $z=z_R$, the beam radius is $w=\sqrt{2}{w}_0\approx 0.707107\mathrm{mm}$ and the on-axis irradiance is half its waist value. At the same wavelength, halving the beam waist radius to $0.250\mathrm{mm}$ quarters the Rayleigh range to $0.310286885\mathrm{m}$ and doubles the beam divergence half-angle to $0.805705984\mathrm{mrad}$. At fixed power, that smaller beam waist also quadruples the peak irradiance.

The contours in the left panel are $r=\pm w(z)$ in a plane through the optical axis; they mark an light intensity level, not an outer boundary. The right panel compares the exact model radii with the far-field asymptotes $w=\theta z$. The horizontal and vertical axes use different length units, so the drawn slopes are not physical angles.

The normalization expresses conserved power during free propagation. In transverse polar coordinates, an annulus has area $2\pi rdr$, soA proof outline is the substitution $s=2r^2/w(z{)}^2$ in this integral. Spreading reduces peak irradiance as $w^{-2}$ while increasing transverse area as $w^2$, leaving power unchanged. In particular, the disk of radius $w(z)$ contains $1-e^{-2}\approx 0.8646647168$, or $86.4665\%$, of the total power; the local light intensity at its rim is only $13.5335\%$ of the on-axis light intensity. These are different quantities.

The left panel compares normalized wave amplitude magnitude and irradiance; the right panel integrates the transverse area to show enclosed power. For an ideal absorbing circular aperture of radius $a$ centered on the optical axis, the same fraction gives the power passing the aperture at that plane. An aperture of radius $w$ clips about $13.5335\%$ of the incident power; radius $2w$ clips only $e^{-8}\approx 0.000335463$, or $0.0335463\%$. After clipping, extra diffraction means the transmitted field no longer obeys the unmodified ideal Gaussian beam law. Every finite aperture cuts some ideal tail. Kogelnik and Li distinguish their large-aperture beam analysis from their finite-aperture diffraction analysis. ^[ref].

The useful tradeoff is that a small bright spot cannot remain narrow indefinitely. An optical lens can create a smaller beam waist, but diffraction then makes the beam spread faster. The particularly simple relation $w_0\theta =\lambda /\pi$ ties localization directly to beam divergence. A proof outline for the propagation law starts with $U(x,y,0)=U_0{e}^{-(x^2+y^2)/w_0^2}$. The parabolic wave equation separates into two transverse factors, and a Gaussian integral in each coordinate givesHere $q$ is dimensionless and $k$ is the wavenumber; the carrier convention is $e^{i(kz-\omega t)}$, with angular frequency $\omega$. Taking the squared complex modulus gives the stated $w(z)$ and its $w^{-2}$ light intensity scaling. Each transverse factor contributes $q^{-1/2}$, explaining the $q^{-1}$ prefactor here and the different prefactor in the preserved Gaussian beam with one transverse coordinate.

Historically, Gaussian beam theory developed alongside early laser resonator theory and lens-guided transmission. In their 1966 review, Kogelnik and Li credit Boyd and Gordon, and Goubau and Schwering, with early beam treatments; they summarize earlier work rather than claim the whole construction as a new discovery. Their paper brings beam propagation, focusing and finite-aperture diffraction into a common practical account. ^[ref].

The Huygens-Fresnel principle describes propagation by coherent addition of secondary-wave contributions from a wavefront or aperture. Their complex amplitudes, including relative wave phases, must be summed before calculating light intensity. The Fresnel propagator follows by retaining quadratic propagation wave phase in the paraxial approximation.

Image formation assigns each object point a corresponding image point: outgoing light rays from that object point either converge there, giving a real image, or diverge as though they originated there, giving a virtual image. The requirement applies separately to each point. Bringing unrelated light rays together makes a bright patch, but does not preserve the object's pattern. This point-by-point correspondence explains how an optical camera records a scene. ^[ref].

Use the same convention as thin lens: place the optical lens at $x=0$, with incident light travelling to the right and a real object at $x=-u$, $u>0$. Write the image coordinate as $x=v$, so $v>0$ for a real image and $v<0$ for a virtual image. Heights $h_o,h_i$ are positive above the optical axis and negative below it. In air, within the paraxial approximation, a thin lens obeysHere $f$ is the signed focal length, positive for a converging optical lens and negative for a diverging optical lens; the expression for $v$ assumes $u\ne f$. Object distance $u$ is positive even though the object's coordinate is negative. These coordinate rules describe transmission through the optical lens; a mirror reverses the outgoing direction, so do not identify its image distance with this same rightward coordinate. ^[ref].

The beauty of geometrical optics is that two convenient light rays can locate a point, while a whole admitted bundle forms its image. For a converging thin lens, an incident light ray parallel to the optical axis exits through the rear optical focus; a light ray through the centre continues undeviated. Real optical lens elements do not focus every admitted light ray exactly: optical aberrations and diffraction limit the point correspondence. The diagrams below are original ideal-model constructions, not measured data. ^[ref].

An optical microscope whose microscope objective, with the specimen at its designed front focal plane, sends out a parallel bundle of light rays from each specimen point. Different specimen points generally give different bundle directions. A tube lens then forms the finite intermediate microscope image; the microscope eyepiece views it. This parallel region permits suitable optical accessories to be inserted. Finite-conjugate and infinity-corrected components must be matched to their intended system.^[ref]

In the ideal paraxial approximation, an object height $h_o$ at objective focal length $f_o$ gives bundle slope $-h_o/f_o$. A tube lens of focal length $f_t$ maps that slope to $h_i=-f_t{h}_o/f_o$, soThis is a two-assembly result. Substituting infinity into the single-objective transverse magnification $-v_o/u_o$ would not describe the finite recorded image. Infinity correction concerns the region after the objective; relaxed-eye viewing concerns the region after the eyepiece. They are distinct conditions.

An objective aperture selects outgoing Electron directions in the objective lens's back focal plane. Passing the central transmitted beam gives bright-field transmission electron microscopy; selecting a separated diffracted beam gives dark-field transmission electron microscopy.^[ref]

In a paraxial approximation, an outgoing angle $\theta$ maps approximately to position $r=f\theta$ in this focal plane, where $f$ is the objective focal length. A centred opening of radius $a$ therefore admits $|\theta |\lesssim a/f$. Selecting a direction selects a component of the outgoing wavefunction; the subsequent image still assigns intensity to specimen positions. The opening is not a hole cut into a real-space specimen image. Its angular selection must be distinguished from selecting a specimen region for electron diffraction.

An optical microscope uses light and optical lens assemblies to reveal small structures, such as the boundaries of biological cells or defects in a manufactured surface. In a compound optical microscope, the microscope objective supplies an enlarged real image; the microscope eyepiece makes that intermediate microscope image subtend a larger viewing angle. The transverse magnification of the first stage and the angular magnification of the second answer different questions.^[ref]

Consider an original ideal finite-conjugate microscope in air, with both optical lens assemblies treated as converging thin lens elements in the paraxial approximation. Choose objective focal length $f_o=10\mathrm{mm}$ and object distance $u_o=12\mathrm{mm}$. The thin lens equation and signed transverse magnification givePlace a microscope eyepiece of focal length $f_e=25\mathrm{mm}$ a further $25\mathrm{mm}$ beyond the intermediate microscope image. The objective-to-eyepiece separation is therefore $85\mathrm{mm}$. Each object point then supplies a parallel outgoing bundle of light rays, appropriate to a relaxed eye focused at infinity. Define the unaided reference as the same object viewed at distance $D=250\mathrm{mm}$; this is a rating convention, not a universal property of eyesight. The eyepiece angular magnification and total signed angular magnification areFor instance, a hypothetical object height $h_o=0.020\mathrm{mm}$ becomes $h_i=-0.100\mathrm{mm}$ at the intermediate microscope image. The signed viewing angles are approximately $h_i/f_e=-0.004\mathrm{rad}$ with the instrument and $h_o/D=0.00008\mathrm{rad}$ unaided. Their ratio is $-50$: enlarged and inverted. These are stipulated model dimensions, not measurements or specifications of a modern commercial optical microscope.

The attraction of an optical microscope is that seeing can become measuring: a stage micrometer supplies a known object-space length, and microscope scale calibration converts an image readout into a specimen length. This also separates a large-looking image from a trustworthy size estimate. Meanwhile, diffraction-limited resolution depends on numerical aperture and wavelength; empty magnification enlarges an existing blur without revealing additional structure.

Use positive carrier wavenumber $k$ and retain waves traveling mainly in the positive $x$ direction. Factor out their rapid wave phase by writing $\psi =e^{ikx}E(x,z)$. Substitution into the Helmholtz equation is exact and givesFor a transverse Fourier mode $e^{i\nu z}$ at a small angle, $|\nu |/k\ll 1$ and its forward longitudinal wavenumber isThus the envelope varies along $x$ on scale $k/{\nu }^2$, much more slowly than the carrier wavelength, and $|E_{xx}|/|k{E}_x|=O({\nu }^2/k^2)$. Dropping this smaller term derives the parabolic wave equation,With transform $\hat{E}(x,\nu )=\int E(x,z)e^{-i\nu z}dz$, its free propagation isThis is a one-way paraxial approximation, not a replacement for backward waves or large-angle Helmholtz components. Its Fourier multiplier has unit modulus and therefore does not attenuate any admitted component.

Factor the forward carrier as $E=e^{ikx}U$. Substitution into the scalar Helmholtz equation gives the exact envelope equationIf $|U_{xx}|\ll 2k|U_x|$, drop the second longitudinal derivative to obtain the parabolic wave equationTypical sufficient conditions are a transverse width $w$ with $kw\gg 1$, longitudinal envelope and material variation on scales much longer than $k^{-1}$, and weak refractive contrast so that factoring out the free-space carrier leaves a slowly varying envelope. A paraxial ordering has transverse propagation angles $O(\eta )$ and ${ϵ}_{\mathrm{rel}}-1=O({\eta }^2)$, with $U_x=O(k{\eta }^{2}U)$ and $U_{xx}$ higher order. A large refractive contrast may require a different, medium-adapted carrier even when the ray direction is nearly axial.

The paraxial approximation retains transverse diffraction and slowly accumulated refraction, but neglects backward propagation, strong reflection from rapid material changes, wide-angle propagation and, in this scalar model, polarization coupling. The supplied vector Helmholtz equation is itself a model assumption: in a source-free inhomogeneous dielectric, Maxwell equations additionally require $\nabla \cdot (ϵ\mathbf{E})=0$, and generally supply a longitudinal-gradient term when the vector wave equation is expanded. An independent scalar component is justified only with a compatible polarization and the relevant small-gradient approximation.

Use the time-harmonic wave convention $\mathrm{Re}(\psi e^{-i\omega t})$, with a positive background wavenumber $k$. The scalar Helmholtz equation for a constant-density acoustic model isWriting $\psi =e^{ikx}E$ gives the exact envelope equationThe paraxial approximation neglects $E_{xx}$ relative to $2ik{E}_x$, giving the parabolic wave equationThe choice of carrier and the direction of propagation matter: this is a forward, slowly varying envelope approximation. Sufficient scale conditions are $|E_{xx}|\ll 2k|E_x|$, transverse spectral components $|p|\ll k$, and medium/envelope variation on longitudinal scales large compared with $1/k$. With the carrier fixed at the background $k$, small $|n-1|$ makes the refractive phase vary slowly too. Large-angle propagation, appreciable backscattering, or rapid longitudinal variation violates the approximation. Acoustic models with variable mass density can have additional gradient terms, so the scalar Helmholtz equation itself is a model assumption. For the Gaussian beam in part (b), useful initial conditions are $kD\gg 1$ and $D/|F|\ll 1$, with the second condition absent for an initially uncurved beam. A narrow angular spectrum, rather than merely the label “Gaussian”, justifies the paraxial approximation.

PutIn free space, the parabolic wave equation is $E_x=i{E}_{zz}/(2k)$. Under the Fourier transform convention $\hat{E}(p)=\int E(z)e^{-ipz}dz$, it becomesThe Gaussian integral is legitimate because $\mathrm{Re}A=1/D^2>0$. Invert the transform after multiplying by $e^{-i{p}^{2}x/(2k)}$. A second Gaussian integral, or equivalently one-dimensional transverse Fresnel propagation, givesChoose the square-root branch continuously from $q(0{)}^{-1/2}=1$; for real $x\ge 0$ there is no zero of $q$. This gives the correct incident field at $x=0$, and direct differentiation verifies the free parabolic wave equation.

For clarity, the squared envelope magnitude isThe Gaussian beam with one transverse coordinate remains Gaussian, with one-transverse-coordinate amplitude factor $q^{-1/2}$, not the $q^{-1}$ of a beam with two transverse coordinates. The negative initial quadratic phase produces focusing for $F>0$; diffraction prevents a singularity at $x=F$. These expressions describe the paraxial approximation to free propagation, rather than an exact unrestricted Helmholtz equation beam.

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.

Use the time convention $e^{-i\omega t}$. The free-space Helmholtz equation, with $\psi =e^{ikx}E$, gives the exact reduced equationFor a forward plane wave at angle $\alpha$, $E_x=O(k{\alpha }^{2}E)$ and $E_{zz}=O(k^2{\alpha }^{2}E)$, whereas $E_{xx}=O(k^2{\alpha }^{4}E)$. Thus the paraxial approximation neglects $E_{xx}$ relative to $2ik{E}_x$. The reduced field satisfies the parabolic wave equationWith the Fourier transform convention used below and its inverse, this becomes an ordinary differential equation for each transverse wavenumber:Consequently the initial-value solution isEquivalently, evaluating the oscillatory Gaussian integral gives the one-dimensional transverse Fresnel propagation formulaThe square-root branch has $\sqrt{1/i}=e^{-i\pi /4}$. The formula is an oscillatory integral for general data, or the Fresnel propagator acting on square-integrable functions. It approaches the initial field as $x\downarrow 0$. A plane wave has transverse wavenumber $\nu =k\sin \alpha$; the reduced axial phase $-{\nu }^{2}x/(2k)$ agrees with $k(\cos \alpha -1)x$ through order ${\alpha }^2$. This also checks the sign. The paraxial approximation applies to $|\nu |\ll k$.

Put $\beta =k\xi \mu$. The random phase screen producesInterpret the normality assumption as a jointly Gaussian process, not merely Gaussian one-point marginal distributions. Write its normalized autocorrelation as $\rho (\zeta )=⟨W(z+\zeta )W(z)⟩$. The difference $W(z+\zeta )-W(z)$ is a zero-mean Gaussian random variable with variance $2(1-\rho (\zeta ))$. Its characteristic function therefore gives the Gaussian phase-screen correlationWithout joint Gaussianity the autocorrelation alone does not determine this expectation.

For an acoustic pressure amplitude $\psi$, with background mass density ${\rho }_0$, the time-averaged acoustic energy flux isAn acoustic potential convention changes the dimensional prefactor. The physical field includes the carrier $e^{ikx}$: a transverse Fourier transform mode has axial wavenumber ${\kappa }_{\nu }=k-{\nu }^2/(2k)$ in the parabolic wave equation. Thus its spectral acoustic flux is proportional to ${\kappa }_{\nu }|\hat{E}(x,\nu ){|}^2$, or to $k|\hat{E}(x,\nu ){|}^2$ at leading paraxial order. Differentiating the reduced field alone would omit the carrier contribution.

The infinite stationary illumination has a spectral measure of a stationary random field, rather than a finite ordinary value of $⟨|\hat{E}{|}^2⟩$ at each $\nu$. To specify the normalization, truncate the screen to an aperture of length $L$, propagate this truncated initial field, and setThe propagation multiplier has modulus one. Using the Gaussian phase-screen correlation and the overlap of the two aperture intervals givesMultiplication by ${\kappa }_{\nu }/(2{\rho }_0\omega )$ gives the ensemble-averaged pressure-normalized modal acoustic energy flux for this finite-aperture convention.

The infinite-aperture spectral normalization consistent with the specified Fourier transform isThe integral and limit are understood as spectral measures of a stationary random field when necessary. By Parseval's identity, this is power per unit transverse length: ${\int }_{\mathbb{R}}{S}_E(d\nu )=C_E(0)=1$. The ensemble-averaged axial spectral acoustic flux is independent of $x$:Here the notation includes Dirac delta distributions in the spectral measure of a stationary random field; the leading-order total flux is $I_{\mathrm{inc}}$. The ${\kappa }_{\nu }$ correction is meaningful within $|\nu |\ll k$, the range of the paraxial approximation.

If $\rho (\zeta )\to 0$ and $e^{{\beta }^2\rho (\zeta )}-1$ is integrable, the stationary phase-screen power spectrum separates into coherent and diffuse wave fields:The coherent fraction is $e^{-{\beta }^2}$ and the diffuse fraction is $1-e^{-{\beta }^2}$. If additionally $|\beta |\ll 1$, the stationary phase-screen power spectrum is $(1-{\beta }^2)\delta +{\beta }^2\hat{\rho }+O({\beta }^4)$ with the same Fourier transform convention. Small $\mu$ alone does not imply small accumulated phase $k\xi \mu$.

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 slowly varying envelope $E=\psi e^{-ikx}$, the paraxial approximation to the Helmholtz equation isBecause $n=1+\mu W$ and ${\mu }^2\ll 1$, passage through a sufficiently thin phase screen producesIts modulus is one at the screen exit. Beyond the screen, $n=1$, so the parabolic wave equation is $E_x=i{E}_{zz}/(2k)$. A Taylor expansion in propagation distance givesSincewe findIt follows thator, equivalently, $|E{|}^2=1-x{\varphi }^{\prime \prime }/k+O(x^2)$. Thus random phase curvature produces local focusing and defocusing: free-space diffraction converts phase fluctuations into amplitude fluctuations immediately after the screen.

Withthe parabolic wave equation is $E_x=(D+S)E$. Freeze $S$ at $x_0$, or preferably at the step midpoint. Over a short distance $\Delta x$, Lie-Trotter splitting givesThe reversed ordering has the same first-order accuracy, while symmetric half-steps in $S$ give the more accurate Strang form.

The commutator can be displayed explicitly. If $q=n^2-1$, thenThe splitting assumption requires $\Delta x^2[D,S]E$ to be small relative to $E$. It is favored by a short range step, a transversely smooth refractive index, and a field without unresolved large transverse wavenumbers. Freezing $S$ also requires $\Delta x{S}_x$ to be small. These conditions supplement the one-way and paraxial approximation already used in part i.

The phase-screen substep is pointwise:DefineThen solve the free-diffraction initial-value problemto $x_0+\Delta x$. In transverse Fourier transform variables, this substep is simplyThis is the split-step Fourier method.