Statistical Q-Plate: Randomized OAM Conversion

• Statistical Q-plate is a space-variant birefringent device that converts spin to orbital angular momentum through controlled random phase fluctuations.
• The device leverages spatial disorder or temporal modulation to yield a probabilistic OAM spectrum, introducing mode crosstalk and spectral broadening.
• Analytical and numerical models reveal that key parameters like phase variance and correlation length crucially determine conversion fidelity and OAM bandwidth.

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" (Moriya, 4 Sep 2025), 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 (Piccirillo et al., 2010). 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 (Moriya, 4 Sep 2025). 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 (Moriya, 4 Sep 2025).

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$0-plate" (Piccirillo et al., 2010), the converted component undergoes a spin flip and an OAM shift of $q$1, 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 (Moriya, 4 Sep 2025). 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 (Piccirillo et al., 2010). 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 $q$2 acts as a unitary operator $q$3 with the transformations

$q$4

$q$5

The conversion and non-conversion probabilities are therefore

$q$6

with maximum conversion at $q$7 (Piccirillo et al., 2010).

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 $q$8, so a mode with transverse phase $q$9 is shifted to $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$0 when the polarization flips (Piccirillo et al., 2010). In the disorder-based theory, the same geometric phase mechanism is retained, but the phase $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$1 includes a random component through $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$2 (Moriya, 4 Sep 2025).

For a monochromatic, paraxial, RHCP Gaussian input beam of amplitude $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$3 and waist $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$4, passage through a thin half-wave plate $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$5 with random fast-axis $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$6 cross-polarizes the input into LHCP and adds the geometric phase $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$7. In subsequent coherence and OAM analyses, the random phase appears only through differences $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$8, because a common global random phase does not affect coherence or OAM projections (Moriya, 4 Sep 2025).

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 $l_{\mathrm{out}}=l_{\mathrm{in}}+2\sigma q,$9 and yields a binary distribution over converted and unconverted channels (Piccirillo et al., 2010). 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 $\sigma=+1$0 values (Moriya, 4 Sep 2025). These are different stochastic mechanisms.

3. Gaussian disorder model and exact ensemble-averaged coherence

The disorder model in (Moriya, 4 Sep 2025) treats $\sigma=+1$1 as a stationary, isotropic, zero-mean Gaussian random field with variance $\sigma=+1$2 and normalized spatial correlation function $\sigma=+1$3 characterized by correlation length $\sigma=+1$4. The correlation satisfies

$\sigma=+1$5

with $\sigma=+1$6 and $\sigma=+1$7. The beam-scale control parameter is the dimensionless correlation length

$\sigma=+1$8

Gaussianity permits exact averaging by the Gaussian moment theorem. Defining $\sigma=+1$9, one has $\sigma=-1$0, with $\sigma=-1$1 and $\sigma=-1$2. The resulting exact ensemble-averaged mutual coherence is

$\sigma=-1$3

This expression preserves the ideal azimuthal factor $\sigma=-1$4 and multiplies it by a disorder-dependent damping factor $\sigma=-1$5 (Moriya, 4 Sep 2025).

The physical interpretation is direct. When two points are so close that the disorder is strongly correlated, $\sigma=-1$6, the damping is weak. When they are separated beyond the correlation scale, $\sigma=-1$7 decreases and the coherence is suppressed by the phase variance $\sigma=-1$8. The spectrum is therefore controlled not only by disorder strength but by the competition between beam waist and disorder correlation length through $\sigma=-1$9.

For equal radii, the separation entering the correlation is

$(\sigma_z=+1)$0

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 (Moriya, 4 Sep 2025).

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

The ensemble-averaged OAM power in integer mode $(\sigma_z=+1)$1 is obtained by angular Fourier projection of the mutual coherence at equal radius, followed by radial integration: $(\sigma_z=+1)$2

$(\sigma_z=+1)$3

Substituting the exact coherence yields

$(\sigma_z=+1)$4

Expanding the damping factor in a uniformly convergent Taylor series,

$(\sigma_z=+1)$5

and defining the angular kernel

$(\sigma_z=+1)$6

the paper derives the exact normalized series

$(\sigma_z=+1)$7

using the normalization $(\sigma_z=+1)$8 (Moriya, 4 Sep 2025).

Two rigorous structural results accompany this spectrum. Proposition 1 establishes total power conservation and the zero-disorder limit: $(\sigma_z=+1)$9 Proposition 2 proves absolute convergence for all $\ell=+2q$0 by bounding $\ell=+2q$1 and $\ell=+2q$2, and gives a rigorous tail bound using the upper incomplete gamma function (Moriya, 4 Sep 2025).

For Gaussian spatial correlation, $\ell=+2q$3, the kernel closes analytically in terms of modified Bessel functions: $\ell=+2q$4 or, in dimensionless radius $\ell=+2q$5,

$\ell=+2q$6

The limiting regimes are sharply resolved. In the coarse-disorder limit $\ell=+2q$7, the disorder is nearly constant across the beam waist, so $\ell=+2q$8 in the high-intensity region and

$\ell=+2q$9

Coarse disorder does not broaden the ensemble-averaged OAM spectrum. In the fine-disorder limit $\alpha(r)=q\theta+\delta\alpha(r),$0, $\alpha(r)=q\theta+\delta\alpha(r),$1 for most separations and, to leading order,

$\alpha(r)=q\theta+\delta\alpha(r),$2

Power conservation then forces the remaining fraction $\alpha(r)=q\theta+\delta\alpha(r),$3 into $\alpha(r)=q\theta+\delta\alpha(r),$4 through the higher-order terms. The asymptotic scaling

$\alpha(r)=q\theta+\delta\alpha(r),$5

shows that the broadened component is localized around $\alpha(r)=q\theta+\delta\alpha(r),$6 with an effective width inversely proportional to $\alpha(r)=q\theta+\delta\alpha(r),$7 (Moriya, 4 Sep 2025).

5. Universal scaling structure and numerical corroboration

A central contribution of (Moriya, 4 Sep 2025) is a universal scaling framework that collapses the OAM spectra across disorder regimes. The “super-universal” control parameter is

$\alpha(r)=q\theta+\delta\alpha(r),$8

which maps $\alpha(r)=q\theta+\delta\alpha(r),$9 to $\delta\alpha(r)$0, with $\delta\alpha(r)$1 in the fine-disorder limit, $\delta\alpha(r)$2 in the coarse-disorder limit, and $\delta\alpha(r)$3 at $\delta\alpha(r)$4. A universal coordinate centered at the nominal mode is

$\delta\alpha(r)$5

where $\delta\alpha(r)$6, and the spectrum is written as

$\delta\alpha(r)$7

Because $\delta\alpha(r)$8, $\delta\alpha(r)$9 is a probability density sampled on the lattice $\alpha(\phi)=q\phi$0 (Moriya, 4 Sep 2025).

The limiting forms of the universal spectrum mirror the asymptotic disorder analysis. As $\alpha(\phi)=q\phi$1, $\alpha(\phi)=q\phi$2, corresponding to the ideal, unbroadened spectrum. As $\alpha(\phi)=q\phi$3, the spectrum contains a central spike of weight $\alpha(\phi)=q\phi$4 plus a broadened background of weight $\alpha(\phi)=q\phi$5 that decays rapidly for $\alpha(\phi)=q\phi$6 (Moriya, 4 Sep 2025). 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 $\alpha(\phi)=q\phi$7, the same normalization, and the same OAM projections as in the exact theory. The simulations propagate

$\alpha(\phi)=q\phi$8

and use a unitary angular FFT to ensure that the $\alpha(\phi)=q\phi$9 spectrum reduces to a single bin at $q$00 (Moriya, 4 Sep 2025).

The numerical agreement is quantitative. In the fine-disorder regime $q$01, the power-averaged nominal-mode fraction $q$02 had $q$03 slope $q$04 versus $q$05, slightly shallower than $q$06 because of incoherent background. A coherent-amplitude estimator,

$q$07

recovered $q$08, confirming the exact coherent fraction. The background width obeyed $q$09; at $q$10, the measured slope $q$11 matched the theoretical $q$12 within $q$13. The universal collapse $q$14 was validated across varying $q$15 and $q$16 with pairwise RMS distances $q$17 (Moriya, 4 Sep 2025).

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 $q$18 nm gives the Malus-like relations

$q$19

and the contrast

$q$20

enables calibration of $q$21. If $q$22 is modulated within the detector integration window $q$23, the output becomes a statistical mixture over OAM outcomes. For an input $q$24 and retardation distribution $q$25,

$q$26

Fast modulation $q$27 realizes a mixed state with tunable weights, whereas slow modulation $q$28 yields time-resolved switching between pure outcomes (Piccirillo et al., 2010).

The experimentally realized EOQPs in (Piccirillo et al., 2010) use nematic E7 liquid crystal at $q$29C between ITO-coated glass substrates, driven by $q$30 kHz AC voltage. The measured Freedericksz threshold is $q$31 V for both devices. A thin device, EOQP2, exhibited measured decay time $q$32 s, approximately $q$33 ms, allowing millisecond reprogramming of OAM probabilities. This makes the time-averaged “Statistical Q-Plate” operationally distinct from the spatially disordered plate of (Moriya, 4 Sep 2025), 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$34 case, a Gaussian input is converted to an LG mode with $q$35, and 3D PIC simulations report power conversion efficiency as high as $q$36 for $q$37 generation from $q$38 (Qu et al., 2017). 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 (Moriya, 4 Sep 2025) are device-level. The coherent fraction $q$39 in the nominal mode quantifies SAM-to-OAM fidelity under fine disorder; measuring the central spike’s coherent amplitude yields the phase variance $q$40 directly. The broadened background width scaling as $q$41 permits inference of $q$42. Together, experimental spectra can determine both $q$43 and $q$44 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 $q$45 quantifies the loss of mode purity under high-frequency internal noise (Moriya, 4 Sep 2025). 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.