---
title: Circular-Beam Approximation
url: https://www.emergentmind.com/topics/circular-beam-approximation
type: topic
---

# Circular-Beam Approximation

Searching arXiv for recent and relevant papers on circular-beam approximation and related usages.
Calling arXiv search tool with targeted queries.
In recent arXiv literature, the expression *circular-beam approximation* designates more than one construction. One use models atmospheric free-space channels by assuming that the instantaneous beam profile at the receiver remains a circular Gaussian with random beam-spot size and wandering but no ellipticity; another uses a circular-beam approximation in eigen-emittance tailoring to obtain a simple relation between transverse eigen-emittance variation and the beam phase integral; and a separate paraxial-optics literature studies Circular-Beams as normalized monochromatic solutions carrying Orbital Angular Momentum [2507.12947] [2104.06818] [1501.07062]. A plausible implication is that the term should be interpreted from context rather than treated as a single standardized formalism.

## 1. Scope and terminological usage

The principal uses of the term in the cited literature can be organized as follows.

| Context | Object modeled | Defining statement |
|---|---|---|
| Paraxial optics | Circular-Beams (CiBs) | “a very general solution of the paraxial wave equation carrying Orbital Angular Momentum” |
| Atmospheric quantum channels | PDT of transmittance | “a simplified alternative” to the elliptic-beam approximation |
| Eigen-emittance tailoring | Coupled transverse beam dynamics | a relation between eigen-emittance difference and beam phase integral |

For CiBs, the object is not an approximation in the narrow sense but a beam family with explicit normalization, Laguerre–Gaussian decomposition, far-field divergence, and experimentally relevant special cases [1501.07062]. For turbulent channels, the approximation is explicitly a model reduction: the elliptic-beam description tracks two spot sizes plus orientation, whereas the circular-beam model tracks only a scalar beam-spot variable \(S\), reducing integrals over beam-shape parameters from \(4\)D down to \(1\)D [2507.12947]. For eigen-emittance tailoring, the approximation applies under experimentally relevant conditions and yields the working formula \(\Delta\varepsilon(s)=|\mathcal W(s)|\), with \(\mathcal W\) the beam phase integral or vorticity [2104.06818].

## 2. Circular-Beams as paraxial OAM solutions

Circular-Beams are defined in cylindrical coordinates \((r,\phi,z)\), with time dependence \(e^{-i\omega t}\) understood. Two complex beam parameters are introduced,
\[
q_0=i z_0-d_0,\qquad q_1=i z_1-d_1,
\]
with \(z_0=\pi W_0^2/\lambda\), and then
\[
q(z)=z+q_0,\qquad \widetilde q(z)=z+q_1,
\]
\[
\xi=\frac{q_1-q_0}{q_0^*-q_1},\qquad \frac{1}{\chi^2(z)}=\frac{i k}{2}\Bigl[\frac{1}{q(z)}-\frac{1}{\widetilde q(z)}\Bigr],
\]
where \(k=2\pi/\lambda\). The normalized, monochromatic Circular–Beam carrying orbital index \(m\in\mathbb Z\) and radial parameter \(p\in\mathbb C\) is [1501.07062]
\[
\mathrm{CiB}_{p,m}^{(q_0,q_1)}(r,\phi,z)
=
N_{p,m}
\Bigl[(1+\xi)\frac{\widetilde q(z)}{q(z)}\Bigr]^{p/2}
\Bigl[\frac{\sqrt2\,r}{W_0}\Bigr]^{|m|}
\,{}_1F_1\!\Bigl(-\frac p2,|m|+1;\frac{r^2}{\chi^2(z)}\Bigr)
e^{i m\phi}
e^{-\frac{i k r^2}{2q(z)}-i(p/2)\arg(1+\xi)+i(|m|+1)\zeta(z)},
\]
with \(\zeta(z)=\arg(q(z))\).

Normalization is imposed by
\[
1=\int_0^{2\pi}d\phi\int_0^\infty |\mathrm{CiB}_{p,m}|^2\, r\,dr,
\]
and the closed form is
\[
|N_{p,m}|^{-2}
=
{}_2F_1\!\Bigl(-\frac p2,-\frac{p^*}{2};1+|m|;|\xi|^2\Bigr).
\]
Convergence, and hence square-integrability, holds under the conditions listed in eq. (8) of the paper, in particular whenever \(|\xi|<1\) or \(p=2\ell\in 2\mathbb N\) if \(|\xi|>1\), etc. [1501.07062].

Because \(\{\mathrm{LG}_{n,m}(r,\phi,z)\}_{n=0}^\infty\) is an orthonormal basis for the paraxial wave equation, the CiB admits the expansion
\[
\mathrm{CiB}_{p,m}^{(q_0,q_1)}(r,\phi,z)
=
\sum_{n=0}^{\infty} C_{n,m}(z)\,\mathrm{LG}_{n,m}(r,\phi,z-d_0),
\]
with
\[
C_{n,m}(z)
=
N_{p,m}^{-1}\,\xi^n\,
\frac{\Gamma(n-\frac p2)}{\Gamma(-\frac p2)}
\sqrt{\frac{|m|!}{n!(|m|+n)!}}.
\]
The coefficients are independent of \(z\) except for the common Gouy-phase in the \(\mathrm{LG}\) modes. This establishes CiBs as a parametrized OAM-carrying family with an explicit modal representation in the standard Laguerre–Gaussian basis.

## 3. Asymptotics, divergence, and experimental generation of CiBs

Near the axis, \(r\ll W(z)\), one has
\[
{}_1F_1\!\Bigl(-\frac p2,|m|+1;\frac{r^2}{\chi^2}\Bigr)\approx 1,
\]
so that \(\psi\sim r^{|m|}e^{i m\phi}\) and
\[
I(r,z)\approx |N_{p,m}|^2\Bigl|\,(1+\xi)\frac{\widetilde q}{q}\Bigr|^p
\Bigl[\frac{\sqrt2\,r}{W_0}\Bigr]^{2|m|}.
\]
Hence \(I\propto r^{2|m|}\) as \(r\to 0\) [1501.07062].

In the far field, \(z\gg z_0,z_1\), the asymptotic forms \(q(z)\approx z\), \(\widetilde q(z)\approx z\), \(\chi^2(z)\approx -2iz/k\), and \(W(z)\approx (W_0/z_0)\,z\) imply that the CiB acquires an effective Gaussian envelope in angle \(\theta=r/z\). Equivalently, from the LG expansion,
\[
I(r,z)\simeq \frac{2}{\pi W^2(z)}\exp\!\Bigl[-\frac{2r^2}{W^2(z)}\Bigr]\,R\!\Bigl(\frac{r}{W(z)}\Bigr),
\]
where \(R(u)\) is a slowly varying polynomial in \(u\) determined by the coefficients \(C_{n,m}\). The overall divergence half-angle (rms) is
\[
\theta_{\rm rms}
=
\theta_0\sqrt{\,1+|m|+\Phi_{p,m}^{(\xi)}-\Re\{\xi(p-\Phi_{p,m}^{(\xi)})\}\,},
\qquad
\theta_0=\frac{W_0}{\sqrt2\,z_0},
\]
with
\[
\Phi_{p,m}^{(\xi)}
=
\frac{|p\xi|^2}{2+2|m|}
\frac{{}_2F_1(1-\frac p2,1-\frac{p^*}{2};2+|m|;|\xi|^2)}
{{}_2F_1(-\frac p2,-\frac{p^*}{2};1+|m|;|\xi|^2)}.
\]
For the special hypergeometric case \(|\xi|=1\), a simpler closed form is available [1501.07062].

The experimentally relevant special cases are explicit. For phase-only OAM modes, \(p=-|m|\) and \(|\xi|=1\), one merely needs a Gaussian TEM\(_{00}\) input passed through a phase-only element such as a q-plate or spiral phase plate \(\exp(i m\phi)\); no amplitude mask is needed. For higher-order radial shapes, \(p>-|m|\) and \(|\xi|=1\), an amplitude mask with transmittance \(\propto r^{p+|m|}\) plus the same \(\exp(i m\phi)\) yields the full CiB of given \((p,m)\). Since a CiB is specified by \((q_0,q_1,p,m)\), standard Gaussian-beam ABCD laws apply separately to \(q_0\) and \(q_1\), and \(\xi\) is invariant. One may thus design lenses or telescopes to adjust divergence or focus [1501.07062].

## 4. Circular-beam approximation for turbulent atmospheric quantum channels

In a free-space atmospheric quantum channel, the input–output relation for the Glauber–Sudarshan \(P\)-function is
\[
P_{\rm out}(\alpha)=\int_0^1 d\eta\;\mathcal P(\eta)\;\frac{1}{\eta}\;
P_{\rm in}\!\Bigl(\frac{\alpha}{\sqrt\eta}\Bigr),
\]
where \(\eta\in[0,1]\) is the random channel transmittance and \(\mathcal P(\eta)\) is the PDT. All statistics of the fluctuating channel enter via \(\mathcal P(\eta)\); once \(\mathcal P(\eta)\) is known, one can compute output photon-number distributions, quadrature variances, and nonclassicality witnesses such as Mandel \(Q\) via averages over \(\eta\) [2507.12947].

The circular-beam approximation assumes that the instantaneous beam profile at the receiver remains a circular Gaussian with random beam-spot size and wandering but no ellipticity. Three fluctuating random variables characterize the beam: the centroid position \(\mathbf r_0=(x_0,y_0)\) with \(\mathrm{var}(x_0)=\sigma_{\rm bw}^2\); the instantaneous squared beam-spot radius
\[
S=4\int_{\mathbb R^2} d^2r\;(x-x_0)^2\,I(r;L),
\]
modeled by a log-normal distribution; and the transmittance
\[
\eta=\int_{\mathcal A} d^2r\;I(r;L),
\]
where \(\mathcal A\) is the circular receiver aperture of radius \(a\). The key simplification versus the elliptic-beam model is that one need only track scalar \(S\), instead of two principal axes and orientation, reducing integrals over beam-shape parameters from \(4\)D down to \(1\)D [2507.12947].

For fixed \(S\), the conditional PDT is obtained from the beam-wandering model. For a circular Gaussian of spot size \(\sqrt{S/2}\) centered at random \(\mathbf r_0\sim\mathcal N(0,\sigma_{\rm bw}^2)\),
\[
\mathcal P(\eta\mid S)
=
\frac{R^2(S)}{\sigma_{\rm bw}^2\,\lambda(S)}
\biggl(\ln\frac{\eta_0(S)}{\eta}\biggr)^{\frac{2}{\lambda(S)}-1}
\exp\!\Bigl[-\frac{R^2(S)}{2\sigma_{\rm bw}^2}
\bigl(\ln\tfrac{\eta_0(S)}{\eta}\bigr)^{2/\lambda(S)}\Bigr]
\]
for \(0\le \eta\le \eta_0(S)\), with
\[
\eta_0(S)=1-\exp(-a^2/S),
\]
and explicit functions \(\lambda(S)\) and \(R(S)\) given in terms of modified Bessel functions \(I_0\) and \(I_1\). Marginalizing over the log-normal \(P(S)\) yields the full PDT through a one-dimensional integration over \(S\) [2507.12947].

## 5. Moment matching, validity range, and quantum observables

The first two moments of transmittance are
\[
\langle \eta\rangle
=
\int_{\mathcal A} d^2r\;\langle I(r;L)\rangle
=
\int_0^1 d\eta\;\eta\,\mathcal P(\eta),
\]
\[
\langle \eta^2\rangle
=
\int_{\mathcal A} d^2r_1 d^2r_2\;\langle I(r_1;L)I(r_2;L)\rangle
=
\int_0^1 d\eta\;\eta^2\,\mathcal P(\eta).
\]
Two procedures are given to determine the log-normal parameters \((\mu,\sigma)\) for \(P(S)\). One matches beam-spot moments \(\langle S\rangle\) and \(\langle S^2\rangle\), then uses
\[
\mu=\ln\!\frac{\langle S\rangle^2}{\sqrt{\langle S^2\rangle}},\qquad
\sigma^2=\ln\!\frac{\langle S^2\rangle}{\langle S\rangle^2}.
\]
The other matches transmittance moments by solving numerically for \((\mu,\sigma)\) from the system
\[
\langle\eta\rangle=\int_0^\infty dS\,P(S\mid \mu,\sigma)\,\langle\eta\rangle_S,\qquad
\langle\eta^2\rangle=\int_0^\infty dS\,P(S\mid \mu,\sigma)\,\langle\eta^2\rangle_S,
\]
where explicit closed forms for \(\langle\eta\rangle_S\) and \(\langle\eta^2\rangle_S\) exist in terms of Marcum–\(Q\) functions [2507.12947].

The comparison with the elliptic-beam approximation is specific. The elliptic-beam model matches experiments best near \(a\approx 1.25\,W_{\rm LT}\) for moderate turbulence, but circular-beam with transmittance-moment matching extends good agreement to a wider range of \(a/W_{\rm LT}\). By matching \(\langle\eta\rangle\) and \(\langle\eta^2\rangle\) directly, the circular-beam PDT reproduces the true first two moments for all \(a\), even outside the narrow elliptic-beam regime. The limitation is equally explicit: the method requires accurate inputs for \(\langle\eta\rangle\) and \(\langle\eta^2\rangle\), and errors in analytical moment formulas propagate into PDT biases, especially for bright inputs [2507.12947].

The model is used for several quantum observables. For sub-Poissonian light,
\[
Q_{\rm out}
=
\frac{\langle\eta^2\rangle}{\langle\eta\rangle}\,Q_{\rm in}
+
\frac{\langle\Delta\eta^2\rangle}{\langle\eta\rangle}\,\langle n\rangle_{\rm in}.
\]
For click-detector nonclassicality,
\[
Q_N
=
N\,\frac{\langle\Delta n^2\rangle}{\langle n\rangle\,(N-\langle n\rangle)}-1.
\]
For quadrature squeezing,
\[
\langle:\Delta x^2:\rangle_{\rm out}
=
\langle\eta\rangle\,\langle:\Delta x^2:\rangle_{\rm in}
+
\langle\Delta T^2\rangle\,\langle x\rangle_{\rm in}^2,
\qquad T=\sqrt\eta.
\]
Using the circular-beam PDT with moment matching, \(Q_{\rm out}\) is found in excellent agreement with phase-screen simulation; the model reproduces simulated \(Q_7\) within a few percent; and accurate modeling of \(\langle\Delta T^2\rangle\) allows prediction of loss-induced anti-squeezing and verification of operating regimes for squeezed-state QKD over free-space links. The summary given in the source attributes to the method a simple \(1\)D-integral PDT model, extended range of validity across beam-aperture ratios and moderate turbulence, sufficient accuracy for predicting nonclassicality transfer in quantum communication, and orders-of-magnitude reduction in computational effort versus the full elliptic-beam model [2507.12947].

## 6. Circular-beam approximation in eigen-emittance tailoring

In beam dynamics, the circular-beam approximation concerns transverse eigen-emittances rather than optical field envelopes. Let \(f_b(x,x',y,y')\) be the normalized \(4\)D beam distribution, with second moments \(\langle uv\rangle\), and define the rms area
\[
A=\sqrt{\langle x^2\rangle\langle y^2\rangle-\langle xy\rangle^2}.
\]
The beam phase integral, also called the beam vorticity, is
\[
\mathcal W
=
\int_A [\nabla\times \bar{\mathbf r}'(x,y)]\cdot d\mathbf A
=
\oint_{\mathcal C}\bar{\mathbf r}'\cdot d\mathbf C,
\]
with centroid-slope field
\[
\bar x'(x,y)=\frac{\langle x'x\rangle}{\langle x^2\rangle}x+\frac{\langle x'y\rangle}{\langle y^2\rangle}y,\qquad
\bar y'(x,y)=\frac{\langle y'x\rangle}{\langle x^2\rangle}x+\frac{\langle y'y\rangle}{\langle y^2\rangle}y.
\]
Its algebraic form is
\[
A\mathcal W
=
\langle y^2\rangle\langle x y'\rangle
-
\langle x^2\rangle\langle y x'\rangle
+
\langle xy\rangle(\langle xx'\rangle-\langle yy'\rangle).
\]
Under experimentally relevant conditions, the paper introduces an approximation providing a very simple and powerful relation between transverse eigen-emittance variation and the beam phase integral [2104.06818].

The derivation starts from the extended Busch theorem. For short coupling elements, such as solenoid fringes, one may assume: the beam rms area is effectively constant, \(\delta A=0\); behind the element all inter-plane correlations vanish, so \(\mathcal W\to 0\); and the change of \(\mathcal W\) is dominated by the change of field, \(\delta\mathcal W=-A\,\delta(B_s/(B\rho))\). Under these approximations,
\[
\frac{A B_s}{(B\rho)}+\mathcal W=\text{constant},
\qquad
(\Delta\varepsilon)^2-\mathcal W^2=\text{constant}.
\]
If the beam enters the coupling region uncorrelated, \((\Delta\varepsilon_i=0,\mathcal W_i=0)\), then
\[
\Delta\varepsilon(s)=|\mathcal W(s)|
\]
becomes the working formula. Within these conditions—short coupling sections, fully decoupled output, no substantial evolution of \(\mathcal W\) in subsequent drifts, and constant beam rigidity—one may use \(\Delta\varepsilon=\mathcal W\) without recourse to the full \(4\times 4\) moment-matrix diagonalization [2104.06818].

The interpretation is that any beam can be viewed as the superposition of two circular modes rotating with opposite angular velocities. Following Brinkmann et al. (1999), the inverse round-to-flat projector maps the beam into two rigidly rotating \(2\)D objects at frequencies \(\pm\omega\), and the two \(4\)D eigen-emittances coincide with the areas of these two oppositely rotating circular modes. The paper states that the difference of eigen-emittances is given by the beam phase integral or vorticity rather than by angular momentum. For perfectly round beams, \(\mathcal W\) reduces to twice the rms angular momentum and one recovers the Kim formula \(\Delta\epsilon=2\mathcal L\); the result thus extends the angular-momentum picture to arbitrary non-cylindrical beams by replacing angular momentum with the more general phase integral [2104.06818].

The practical examples given are a flat-beam photoinjector, where \(\Delta\epsilon_f=\mathcal W_f=(A_0B_0)/(B\rho)\), and the EMTEX emittance-transfer experiment, where \(\Delta\epsilon_{\rm final}=\sqrt{\Delta\epsilon_{\rm initial}^2+\mathcal W_{\rm after}^2}\), in excellent agreement with measured \(\sim 2.2\) mm mrad. Simple model objects include rigid rotation at constant \(\omega\), giving \(\Delta\epsilon=2|\omega|A\), and pure shear \(x'=a\,y\), giving \(\Delta\epsilon=|a|A\) [2104.06818].

## 7. Related circular-symmetry approximations and common distinctions

A separate paraxial approximation concerns beams passing apertures with sharp boundaries rather than the two circular-beam approximations just described. For a uniform circular aperture of radius \(R\), with initial field \(U(r,0)=A\) for \(r\le R\) and zero outside, the sharp-boundary method gives
\[
U(r,z)\approx
A\,e^{ikz}\,\frac12\,\operatorname{erfc}\!\Biggl(\frac{r-R}{\sqrt{2\,i\,z/k}}\Biggr).
\]
No other fitting parameters appear; the edge-thickness parameter is taken \(\to 0\) for a truly sharp boundary [1605.04540].

Its regime of validity is also explicit: paraxial propagation, \(z\gg \lambda\), with the physical transition from inside to outside occurring over a transverse scale \(\Delta r\ll R\). The approximation is most accurate in the region \(r\approx R\pm O(\sqrt{z/k})\) near the edge. It captures the onset of Fresnel diffraction as edge-blurring, with the initially sharp jump at \(r=R\) smoothed over a radial width \(O(\sqrt{z/k})\), but it does not generate the full oscillatory ring structure of the Airy pattern [1605.04540].

This distinction matters conceptually. In the atmospheric-channel model, “circular” refers to the assumed instantaneous beam shape at the receiver; in eigen-emittance tailoring, it refers to the decomposition into two circular modes rotating with opposite angular velocities; and in the sharp-boundary diffraction problem, it refers to the geometry of the aperture. This suggests that the phrase *circular-beam approximation* is best read as a family resemblance across circular symmetry reductions, rather than as a unique formalism.

Source: https://www.emergentmind.com/topics/circular-beam-approximation