---
title: 'Statistical Q-Plate: Randomized OAM Conversion'
url: https://www.emergentmind.com/topics/statistical-q-plate
type: topic
---

# Statistical Q-Plate: Randomized OAM Conversion

A statistical Q-plate is a q-plate whose spin-to-orbital angular momentum conversion is governed by randomness rather than by a perfectly deterministic geometric-phase pattern. In the exact disorder-theoretic formulation of "Normalized Ensemble-Averaged OAM Spectrum in Disordered Statistical Q-Plates" [2509.03920], it is a thin, birefringent geometric-phase element whose fast-axis orientation contains random spatial fluctuations around the ideal azimuthal law. In an earlier control-oriented usage, an electrically tunable liquid-crystal q-plate becomes a "Statistical Q-Plate" when its retardation is modulated within the detector integration window, so that the measured output is a statistical mixture over OAM outcomes [1010.4473]. The common core is that a q-plate continues to mediate SAM-to-OAM conversion, but the output OAM content must be described probabilistically rather than as a single deterministic mode.

## 1. Definition and conceptual scope

In the disorder-based formulation, the ideal optic-axis distribution is  
\[
\alpha(r,\theta)=q\theta,
\]
where \(q\) is the topological charge. A circularly polarized input then undergoes SAM-to-OAM conversion according to the nominal rule  
\[
l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,
\]
under the thin-element and paraxial approximations, with \(\sigma=+1\) for RHCP and \(\sigma=-1\) for LHCP in the notation of [2509.03920]. For an RHCP beam \((\sigma_z=+1)\), the ideal device outputs LHCP with OAM \(\ell=+2q\). A statistical Q-plate replaces the ideal axis by  
\[
\alpha(r)=q\theta+\delta\alpha(r),
\]
where \(\delta\alpha(r)\) is a random, zero-mean fluctuation field. These microscopic birefringent fluctuations generate macroscopic OAM crosstalk and spectral broadening [2509.03920].

A q-plate more generally is a space-variant birefringent plate whose optic axis winds azimuthally as \(\alpha(\phi)=q\phi\), thereby coupling circular polarization to OAM. In the liquid-crystal formulation of "Photon spin-to-orbital angular momentum conversion via an electrically tunable \(q\)-plate" [1010.4473], the converted component undergoes a spin flip and an OAM shift of \(\pm 2q\), while the unconverted component remains in its original OAM state. This broader q-plate framework supplies the deterministic baseline from which statistical variants are defined.

The term therefore spans at least two distinct constructions. One is **spatially disordered**, in which randomness is internal to the device microstructure and is modeled as a random field over the transverse plane [2509.03920]. The other is **temporally randomized**, in which the device is itself deterministic at each instant but is driven so that measurement averages over different retardations and hence over different conversion probabilities [1010.4473]. This suggests that “statistical Q-plate” is best understood as a family of stochastic q-plate models rather than a single implementation class.

## 2. Ideal q-plate action and the origin of statistical behavior

In the circular polarization and OAM basis, an ideal q-plate with retardation \(\delta\) acts as a unitary operator \(U_q(\delta)\) with the transformations  
\[
|L, l\rangle \xrightarrow{U_q(\delta)} \cos\!\left(\frac{\delta}{2}\right)|L, l\rangle
+ i\,\sin\!\left(\frac{\delta}{2}\right)|R, l+2q\rangle,
\]
\[
|R, l\rangle \xrightarrow{U_q(\delta)} \cos\!\left(\frac{\delta}{2}\right)|R, l\rangle
+ i\,\sin\!\left(\frac{\delta}{2}\right)|L, l-2q\rangle.
\]
The conversion and non-conversion probabilities are therefore  
\[
P_{\text{conv}}(\delta)=\sin^2\!\left(\frac{\delta}{2}\right),\qquad
P_{\text{no conv}}(\delta)=\cos^2\!\left(\frac{\delta}{2}\right),
\]
with maximum conversion at \(\delta=\pi\) [1010.4473].

The local Jones description makes the geometric origin of OAM transfer explicit. In the circular basis, the off-diagonal terms in the space-variant Jones matrix carry helical phase factors \(e^{\pm i2q\phi}\), so a mode with transverse phase \(e^{il\phi}\) is shifted to \(l\pm 2q\) when the polarization flips [1010.4473]. In the disorder-based theory, the same geometric phase mechanism is retained, but the phase \(2\alpha(r)\) includes a random component through \(\delta\alpha(r)\) [2509.03920].

For a monochromatic, paraxial, RHCP Gaussian input beam of amplitude \(E_0\) and waist \(w_0\), passage through a thin half-wave plate \((\Gamma=\pi)\) with random fast-axis \(\alpha(r)=q\theta+\delta\alpha(r)\) cross-polarizes the input into LHCP and adds the geometric phase \(\exp(i2\alpha)\). In subsequent coherence and OAM analyses, the random phase appears only through differences \(\delta\alpha(r_1)-\delta\alpha(r_2)\), because a common global random phase does not affect coherence or OAM projections [2509.03920].

A frequent source of confusion is the relation between incomplete conversion and disorder-induced broadening. In the electrically tunable setting, incomplete conversion is governed by the deterministic retardation \(\delta\) and yields a binary distribution over converted and unconverted channels [1010.4473]. In the disordered setting, the device is locally half-wave, but random axis fluctuations broaden the converted OAM spectrum itself and redistribute power among nearby \(\ell\) values [2509.03920]. These are different stochastic mechanisms.

## 3. Gaussian disorder model and exact ensemble-averaged coherence

The disorder model in [2509.03920] treats \(\delta\alpha(r)\) as a stationary, isotropic, zero-mean Gaussian random field with variance \(V\) and normalized spatial correlation function \(g(s)\) characterized by correlation length \(L_c\). The correlation satisfies  
\[
\langle \delta\alpha(r)\,\delta\alpha(r')\rangle = V\,g(|r-r'|/L_c),
\]
with \(g(0)=1\) and \(|g(s)|\le 1\). The beam-scale control parameter is the dimensionless correlation length  
\[
\xi=\frac{L_c}{w_0}.
\]

Gaussianity permits exact averaging by the Gaussian moment theorem. Defining \(Y=2(\delta\alpha(r_2)-\delta\alpha(r_1))\), one has \(\langle e^{iY}\rangle=e^{-\langle Y^2\rangle/2}\), with \(\mathrm{Var}(Y)=8V[1-g(s)]\) and \(s=|r_2-r_1|/L_c\). The resulting exact ensemble-averaged mutual coherence is  
\[
= E_0^2 \exp\left(-\frac{r_1^2+r_2^2}{w_0^2}\right)
\exp[i2q(\theta_2-\theta_1)]
\exp[-4V(1-g(s))].
\tag{1}
\]
This expression preserves the ideal azimuthal factor \(\exp[i2q(\theta_2-\theta_1)]\) and multiplies it by a disorder-dependent damping factor \(\exp[-4V(1-g)]\) [2509.03920].

The physical interpretation is direct. When two points are so close that the disorder is strongly correlated, \(g(s)\approx 1\), the damping is weak. When they are separated beyond the correlation scale, \(g(s)\) decreases and the coherence is suppressed by the phase variance \(V\). The spectrum is therefore controlled not only by disorder strength but by the competition between beam waist and disorder correlation length through \(\xi=L_c/w_0\).

For equal radii, the separation entering the correlation is  
\[
s_{r,\phi}=r\sqrt{2-2\cos\phi}/L_c,
\]
so the entire OAM problem reduces to angular Fourier analysis of the damped azimuthal coherence kernel. This exact reduction is central because it connects microscopic birefringent randomness to a measurable macroscopic OAM spectrum without invoking perturbative truncations [2509.03920].

## 4. Exact OAM spectrum, normalization, and limiting disorder regimes

The ensemble-averaged OAM power in integer mode \(\ell\) is obtained by angular Fourier projection of the mutual coherence at equal radius, followed by radial integration:
\[
\Gamma_\ell(r)\equiv \frac{1}{2\pi}\int_{0}^{2\pi} (r,\theta;\,r,\theta+\phi)\,e^{-i\ell\phi}\,d\phi,
\tag{2}
\]
\[
P(\ell;\xi)\equiv 2\pi\int_{0}^{\infty}\Gamma_\ell(r)\,r\,dr.
\tag{3}
\]
Substituting the exact coherence yields
\[
\Gamma_\ell(r)=E_0^2\exp\left(-\frac{2r^2}{w_0^2}\right)\frac{1}{2\pi}
\int_{0}^{2\pi}\exp[i(2q-\ell)\phi]\exp[-4V(1-g(s_{r,\phi}))]\,d\phi.
\tag{4}
\]
Expanding the damping factor in a uniformly convergent Taylor series,
\[
\exp[-4V(1-g(s))]=e^{-4V}\sum_{n=0}^{\infty}\frac{(4V)^n}{n!}[g(s)]^n,
\tag{5}
\]
and defining the angular kernel
\[
C_{n,m}(r;\xi)=\frac{1}{2\pi}\int_{0}^{2\pi}\cos(m\phi)
\left[g\left(\frac{r\sqrt{2-2\cos\phi}}{L_c}\right)\right]^n\,d\phi,
\tag{6}
\]
the paper derives the exact normalized series
\[
P(\ell;\xi)=4e^{-4V}\sum_{n=0}^{\infty}\frac{(4V)^n}{n!}K_{n,|\ell-2q|}(\xi),
\tag{8}
\]
using the normalization \(E_0^2w_0^2=2/\pi\) [2509.03920].

Two rigorous structural results accompany this spectrum. Proposition 1 establishes total power conservation and the zero-disorder limit:
\[
\sum_\ell P(\ell;\xi)=1,\qquad P(\ell;0)=\delta_{\ell,2q}.
\]
Proposition 2 proves absolute convergence for all \(V\ge 0\) by bounding \(|C_{n,m}|\le 1\) and \(|K_{n,m}|\le 1/4\), and gives a rigorous tail bound using the upper incomplete gamma function [2509.03920].

For Gaussian spatial correlation, \(g(s)=e^{-s^2}\), the kernel closes analytically in terms of modified Bessel functions:
\[
C_{n,m}(r;\xi)=\exp\left(-n\frac{2r^2}{L_c^2}\right)
I_m\left(n\frac{2r^2}{L_c^2}\right),
\tag{9}
\]
or, in dimensionless radius \(p=r/w_0\),
\[
C_{n,m}(w_0p;\xi)=\exp\left(-n\frac{2p^2}{\xi^2}\right)
I_m\left(n\frac{2p^2}{\xi^2}\right).
\tag{10}
\]

The limiting regimes are sharply resolved. In the **coarse-disorder** limit \(\xi\gg 1\), the disorder is nearly constant across the beam waist, so \(\delta\alpha(r_2)-\delta\alpha(r_1)\approx 0\) in the high-intensity region and
\[
P(\ell;\xi\gg 1)=\delta_{\ell,2q}.
\tag{11}
\]
Coarse disorder does not broaden the ensemble-averaged OAM spectrum. In the **fine-disorder** limit \(\xi\ll 1\), \(g(s)\approx 0\) for most separations and, to leading order,
\[
P(\ell;\xi\ll 1,\text{ leading order})=e^{-4V}\delta_{\ell,2q}.
\tag{12}
\]
Power conservation then forces the remaining fraction \(1-e^{-4V}\) into \(\ell\neq 2q\) through the higher-order terms. The asymptotic scaling
\[
K_{n,m}(\xi)\sim \xi f_n(m\xi)
\tag{13}
\]
shows that the broadened component is localized around \(\ell=2q\) with an effective width inversely proportional to \(\xi\) [2509.03920].

## 5. Universal scaling structure and numerical corroboration

A central contribution of [2509.03920] is a universal scaling framework that collapses the OAM spectra across disorder regimes. The “super-universal” control parameter is
\[
U=\frac{\xi^2}{1+\xi^2},
\tag{14}
\]
which maps \(\xi\in[0,\infty)\) to \(U\in[0,1]\), with \(U\to 0\) in the fine-disorder limit, \(U\to 1\) in the coarse-disorder limit, and \(U=1/2\) at \(\xi=1\). A universal coordinate centered at the nominal mode is
\[
y=\frac{\ell-2q}{\sigma_{\text{total}}},
\tag{15}
\]
where \(\sigma_{\text{total}}=2\mu_1\), and the spectrum is written as
\[
P(\ell;\xi)=G(y,U).
\tag{17}
\]
Because \(\sum_\ell P(\ell;\xi)=1\), \(G(y,U)\) is a probability density sampled on the lattice \(y_\ell=(\ell-2q)/\sigma_{\text{total}}\) [2509.03920].

The limiting forms of the universal spectrum mirror the asymptotic disorder analysis. As \(U\to 1\), \(G(y,U)\to \delta(y)\), corresponding to the ideal, unbroadened spectrum. As \(U\to 0\), the spectrum contains a central spike of weight \(e^{-4V}\) plus a broadened background of weight \(1-e^{-4V}\) that decays rapidly for \(|y|\gtrsim 1\) [2509.03920]. The conceptual significance is that distinct microscopic disorder strengths and correlation lengths can map to the same universal spectral morphology after rescaling.

Monte Carlo simulations were performed with the same Gaussian correlation model \(g(s)=\exp[-(s/L_c)^2]\), the same normalization, and the same OAM projections as in the exact theory. The simulations propagate
\[
\psi(r,\theta)=E_0e^{-r^2/w_0^2}e^{i(2q\theta+2\delta\alpha(r))}
\]
and use a unitary angular FFT to ensure that the \(V=0\) spectrum reduces to a single bin at \(\ell=2q\) [2509.03920].

The numerical agreement is quantitative. In the fine-disorder regime \((\xi=0.03)\), the power-averaged nominal-mode fraction \(f_0=P(\ell=2q)\) had \(\ln(f_0)\) slope \(-3.288\) versus \(V\), slightly shallower than \(-4\) because of incoherent background. A coherent-amplitude estimator,
\[
f^{\text{coh}}_0=\left|\left\langle\int a_{2q}(r)\,r\,dr\right\rangle\right|^2,
\]
recovered \(f^{\text{coh}}_0=e^{-4V}\), confirming the exact coherent fraction. The background width obeyed \(\sigma_\ell\propto 1/\xi\); at \(V=0.30\), the measured slope \(a_{\text{fit}}=1.30706\) matched the theoretical \(a_{\text{bg}}=1.31043\) within \(0.25\%\). The universal collapse \(P(\ell;\xi)=G(y,U)\) was validated across varying \(V\) and \(\xi\) with pairwise RMS distances \(\lesssim 0.08\) [2509.03920].

## 6. Alternative realizations, related platforms, and application domains

In the electrically tunable liquid-crystal implementation, the statistical character does not arise from spatial disorder but from time-dependent retardation. Measuring converted and unconverted powers at wavelength \(532\) nm gives the Malus-like relations
\[
P_u=P_0\cos^2(\delta/2),\qquad P_c=P_0\sin^2(\delta/2),
\]
and the contrast
\[
\rho(v)=\frac{P_u-P_c}{P_u+P_c}=\cos\delta(v)
\]
enables calibration of \(\delta(V)\). If \(\delta\) is modulated within the detector integration window \(T\), the output becomes a statistical mixture over OAM outcomes. For an input \(|L,l=0\rangle\) and retardation distribution \(f(\delta)\),
\[
P(l=0)=\int f(\delta)\cos^2\!\left(\frac{\delta}{2}\right)d\delta,\qquad
P(l=+2q)=\int f(\delta)\sin^2\!\left(\frac{\delta}{2}\right)d\delta.
\]
Fast modulation \((\tau_{\text{mod}}\ll T)\) realizes a mixed state with tunable weights, whereas slow modulation \((\tau_{\text{mod}}\gg T)\) yields time-resolved switching between pure outcomes [1010.4473].

The experimentally realized EOQPs in [1010.4473] use nematic E7 liquid crystal at \(30^\circ\)C between ITO-coated glass substrates, driven by \(1\) kHz AC voltage. The measured Freedericksz threshold is \(1.85\pm 0.01\) V for both devices. A thin device, EOQP2, exhibited measured decay time \((0.80\pm 0.01)\times 10^{-2}\) s, approximately \(8\) ms, allowing millisecond reprogramming of OAM probabilities. This makes the time-averaged “Statistical Q-Plate” operationally distinct from the spatially disordered plate of [2509.03920], even though both are probabilistic OAM converters.

A related but non-liquid-crystal platform is the plasma q-plate, where a magnetized plasma implements the rotating optic axis and converts a circularly polarized Gaussian beam to a twisted beam. In the demonstrated \(q=1\) case, a Gaussian input is converted to an LG mode with \(\ell=\pm 2\), and 3D PIC simulations report power conversion efficiency as high as \(83\%\) for \(\ell=\pm 2\) generation from \(\ell_{\text{in}}=0\) [1709.07780]. That work also proposes a statistical robustness description in which local retardation fluctuations are modeled as a random field, identifying transverse inhomogeneity as the principal source of mode mixing and crosstalk. This suggests that statistical-q-plate analysis may extend beyond liquid-crystal disorder to plasma and other anisotropic media.

The practical implications stated in [2509.03920] are device-level. The coherent fraction \(e^{-4V}\) in the nominal mode quantifies SAM-to-OAM fidelity under fine disorder; measuring the central spike’s coherent amplitude yields the phase variance \(V\) directly. The broadened background width scaling as \(1/\xi\) permits inference of \(L_c/w_0\). Together, experimental spectra can determine both \(V\) and \(L_c\) by fitting the normalized exact series and by applying the universal rescaling. In OAM communications and quantum state preparation, coarse-disorder protection implies tolerance to slowly varying device imperfections, while the fine-disorder law \(f_0=e^{-4V}\) quantifies the loss of mode purity under high-frequency internal noise [2509.03920]. A plausible implication is that statistical q-plates form a common language for treating stochastic OAM generation across fabrication disorder, active temporal modulation, and complex anisotropic media.

Source: https://www.emergentmind.com/topics/statistical-q-plate