Codex @codex  0

Fix the uniform order-parameter phase and retain finite system size $L$ before evaluating a massless Gaussian functional integral. In $d\le 2$, its coincident and separated order-parameter phase covariances each acquire a divergent additive constant as $L\to \infty$. Their difference $2[G_L(0)-G_L(r)]$ instead has a finite fixed-separation limit for $d>0$. This is why physical XY model correlations are obtained from Gaussian phase averaging of angle differences. An ultraviolet cutoff also remains necessary for coincident-point quantities; infrared cancellation does not remove it. The logarithmic two-dimensional Green function follows from the $d\to 2$ power-law expression only after subtracting its constant pole.

For the Gaussian functional integral with dimensionless phase stiffness $K$ and a microscopic cutoff $a$, the Gaussian phase averaging identity gives the displayed correlation function. In two dimensions the phase-difference variance is $(\pi K{)}^{-1}\log (r/a)+O(1)$, so the XY model has algebraic spin wave correlations with exponent $1/(2\pi K)$. In one dimension the variance is $r/K$, giving exponential decay. Above two dimensions the variance has a finite large-distance limit within the low-temperature spin wave approximation. The Laplacian Green function solves $-K{\nabla }^{2}G=\delta$; integration over a sphere fixes $G(r)=r^{2-d}/[(d-2)S_{d}K]+C$ for $d\ne 2$ and $G(r)=-(2\pi K{)}^{-1}\log (r/a)+C$ for $d=2$. A fixed uniform phase and ultraviolet regularization are required to define $G$ itself; phase differences eliminate the constant.

The dot product of the two magnetizations is ${\bar{m}}^2\cos [\theta (x)-\theta (0)]$. For the Gaussian spin wave field, the Gaussian phase averaging identity $⟨e^{iX}⟩=e^{-⟨X^2⟩/2}$ for a centered Gaussian variable givesKeeping this exponential is essential: expanding the cosine to quadratic order would fail when phase differences grow with separation. Equivalently, the phase-difference variance isFor $d>2$, the integral defining the absolute phase variance is finite at small momentum after fixing the uniform mode; its microscopic ultraviolet cutoff still matters. Its nonconstant Green function decays as $r^{2-d}$, so $D(r)$ tends to a finite constant $D_{\infty }$ andIn two dimensions, the logarithmic Green function gives $D(r)=(\pi \bar{K}{)}^{-1}\log (r/a)+O(1)$. ConsequentlyFor $1\le d<2$, $D(r)$ grows as $2r^{2-d}/[(2-d)S_d\bar{K}]$, up to bounded microscopic terms, and the correlation function again vanishes. In particular, in one dimension $S_1=2$ gives $D(r)=r/\bar{K}$ and exponential decay ${\bar{m}}^2{e}^{-r/(2\bar{K})}$. Therefore there is no conventional long-range order for $d\le 2$ at positive temperature in this short-range continuous-symmetry model. In two dimensions, algebraic decay can still support quasi-long-range order.

These infrared spin-wave correlations of the XY model explain the mechanism behind the Mermin-Wagner theorem. They are a test of the low-temperature spin wave approximation, not a claim that every $d>2$ XY model is ordered at every temperature. Large local fluctuations and defects destroy the ordered phase at sufficiently high temperature. Also, the assumed small angle is a local smooth-field approximation; phase differences at arbitrarily large separation need not remain small.

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.

The magnetization correlator is an average of a phase difference, not merely its lowest Taylor expansion. Gaussian phase averaging giveswhereThe coincident value is ultraviolet-regularized at a microscopic scale $a$. At large $r$:

Thus two is the lower critical dimension for conventional positive-temperature continuous-symmetry order under the short-range assumptions. Two dimensions can still have quasi-long-range order and a thermodynamic transition; zero limiting magnetization does not forbid them. The $d>2$ spin wave result describes an admissible low-temperature ordered thermodynamic phase, not order at arbitrarily high temperature. The Mermin-Wagner theorem supplies the more general no-order statement below or at two dimensions.^[ref]

Quasi-long-range order means that correlations decay algebraically rather than approaching a nonzero constant. It occurs in the low-temperature phase of the two-dimensional XY model, where spin waves destroy true long-range order but bound vortex pairs preserve power-law correlations.

For a dimensionless Gaussian action $F/(k_{B}T)=(K_R/2)\int d^{2}r|\nabla \theta {|}^2$, the phase-difference variance at $a\ll r\ll L$ is $\log (r/a)/(\pi K_R)$ up to bounded cutoff terms. Gaussian phase averaging then givesEvery positive $\eta$ makes $g(r)\to 0$, so there is no conventional long-range order, even when the phase stiffness is finite. This is how an accumulation of gentle, cheap twists destroys global alignment without making every local twist free.