Von Kármán Phase Screens

• Von Kármán phase screens are numerical models representing optical phase perturbations derived from a bounded turbulence spectrum, crucial for simulating atmospheric effects.
• They utilize methods like FFT-based simulation, subharmonic augmentation, and randomized spectral sampling to accurately capture both low- and high-frequency turbulence characteristics.
• Advanced techniques, including autoregressive Fourier updates and 3D Karhunen–Loève constructions, extend static phase screens to time-evolving and multi-wavelength applications.

Von Kármán phase screens are numerical representations of optical phase perturbations whose second-order statistics are prescribed by the von Kármán turbulence spectrum. In optical-turbulence theory, the phase power spectral density is commonly written as

$$W_\varphi(f)=0.023\,r_0^{-5/3}\,\bigl(f^2+f_0^2\bigr)^{-11/6}, \qquad f_0=\frac{1}{L_0},$$

or, equivalently in spatial-frequency notation,

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$

with $r_0$ the Fried coherence length and $L_0$ the outer scale. In the Kolmogorov limit $L_0\to\infty$, the spectrum reduces to the familiar $-11/3$ inertial-range law. These screens are used to model atmospheric turbulence in ground-based astronomy, adaptive optics, laser propagation, and atmospheric wave-optics simulations, both as static realizations and as time-dependent fields (Mathar, 14 Oct 2025, Chhabra et al., 2021, Srinath et al., 2015).

1. Statistical model and structure function

The defining feature of a von Kármán phase screen is that its covariance and structure function are derived from a bounded turbulence spectrum with finite outer scale. For a phase screen $\varphi(\mathbf r)$, the two-point structure function is

$$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$

and in the Kolmogorov limit it follows the $5/3$-law

$D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$

More generally, the Wiener–Khinchin relation gives

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$0

with covariance

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$1

For finite $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$2, a closed-form expression used in FFT-based analysis is

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$3

which tends to $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$4 as $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$5 (Mathar, 14 Oct 2025, Chhabra et al., 2021).

The same model can be written at the refractive-index level. In three dimensions, the von Kármán refractive-index spectrum is given as

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$6

while alternative formulations use

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$7

and, when an inner scale is included,

$$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$8

These equivalent parameterizations underpin different simulation frameworks, including single-layer phase screens, volumetric models, and multi-wavelength cross-spectral synthesis (Mathar, 2011, IV et al., 1 May 2025).

2. Classical frozen-screen construction

Time-dependent phase screens in ground-based astronomy are typically simulated in the frozen-screen approximation by establishing a static phase screen on a large pupil and dragging an aperture equivalent to the size of the actual input pupil across this oversized phase screen. The speed of this motion sweeping through the large phase screen is equivalent to a wind speed that changes the phase screen as a function of time (Mathar, 14 Oct 2025).

The standard FFT-based generator discretizes a square grid of physical size $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$9 on $r_0$0 pixels, with sampling interval $r_0$1, Fourier spacing $r_0$2, and Nyquist frequency $r_0$3. On the discrete frequency grid,

$r_0$4

the DC term $r_0$5 is set to zero to remove piston, complex Gaussian coefficients

$r_0$6

are generated with $r_0$7, and an inverse FFT yields a real-space screen $r_0$8 (Chhabra et al., 2021).

This construction is computationally convenient, but the literature identifies specific failure modes. FFT-based phase-screen simulations give accurate results only when the screen size $r_0$9 is much larger than the outer scale parameter $L_0$0; otherwise, they fall short in correctly predicting both the low and high frequency behaviours of turbulence induced phase distortions. Traditional Fourier methods also suffer from screen periodicity and low spatial frequency power deficiency. In comparison studies, pure DFT misses low $L_0$1 entirely, and subharmonic-augmented DFT still leaves residual bias at low frequencies (Chhabra et al., 2021, Srinath et al., 2015, Charnotskii, 2019).

3. Low-frequency compensation and alternative sampling strategies

A large part of the modern literature on von Kármán phase screens is concerned with recovering low-frequency content without excessive computational cost. The classic subharmonic method adds extra sinusoids at fractional frequencies $L_0$2, $L_0$3, but reported limitations include unequal sampling of low-$L_0$4 annuli and the need for empirical patch-normalization factors that depend on $L_0$5. A generalized approach replaces fixed low-frequency weights with a Gaussian phase autocorrelation matrix: after FFT and subharmonics, one computes a residual

$L_0$6

fits that residual with a Gaussian-kernel model, and adds the corresponding autocorrelation correction to restore missing variance and correlations. The same framework also performs high-frequency compensation by removing piston/tip/tilt, zero-clamping negative power, and fitting a second Gaussian smoothing operator in the high-$L_0$7 band. For $L_0$8 as small as $L_0$9, grid sizes up to $L_0\to\infty$0, and $L_0\to\infty$1 up to $L_0\to\infty$2, the maximum $L_0\to\infty$3 is reported as $L_0\to\infty$4 (Chhabra et al., 2021).

Randomized Spectral Sampling (RSS) replaces uniform sampling of the PSD on the discrete Fourier lattice by a randomized shift of the entire sampling grid. Each realization draws offsets $L_0\to\infty$5 uniformly in the fundamental Fourier cell, samples the spectrum at $L_0\to\infty$6, performs a single inverse FFT, and multiplies by a compensating phase ramp. Because the sampling grid is shifted away from the origin, low-frequency content such as tip, tilt, and focus appears without additional subharmonic subgrids. For bounded von Kármán models with finite outer scale, RSS alone yields RMS structure-function error of $L_0\to\infty$7 for $L_0\to\infty$8 on $L_0\to\infty$9, $-11/3$0, and $-11/3$1 grids, whereas the traditional uniform grid has errors up to $-11/3$2. The same work reports that RSS-generated screens match Cheon’s closed-form angle-of-arrival theory substantially better than traditional FFT screens (Paulson et al., 2019).

A separate line of work uses sparse trigonometric expansions rather than FFT grids. Charnotskii compares subharmonic-augmented DFT, randomized DFT, Sparse Spectrum (SS), and Sparse Spectrum with Uniform segments (SU), all targeted at the same continuous von Kármán PSD (Charnotskii, 2019).

Method	Mechanism	Reported behavior
DFT+SH	FFT grid plus low-frequency subharmonics	residual bias $-11/3$3 remains even with $-11/3$4
PWD / randomized DFT	jittered Cartesian grid, optionally with SH	unbiased in expectation; $-11/3$5 slower than pure FFT
SS	direct sparse cosine or complex-exponential sum	$-11/3$6 structure-function error once $-11/3$7
SU	sparse spectrum with uniform sector sampling	$-11/3$8 error for $-11/3$9 after sufficient averaging

The comparison shows that SS and SU generate unbiased samples for screens with bounded phase variance and show superior computational effectiveness for one megapixel and larger screens. A plausible implication is that the choice of generator is governed less by the turbulence model itself than by how faithfully a numerical scheme tiles the low-$\varphi(\mathbf r)$0 region and how efficiently it handles large domains (Charnotskii, 2019).

4. Time evolution beyond static translation

The most direct time-dependent generalization of a von Kármán phase screen is the autoregressive Fourier-domain update

$\varphi(\mathbf r)$1

where $\varphi(\mathbf r)$2 is the target von Kármán PSD, $\varphi(\mathbf r)$3 is complex white noise, and

$\varphi(\mathbf r)$4

recovers pure frozen flow. Allowing $\varphi(\mathbf r)$5 injects stochastic refresh, or “boiling,” while preserving the desired stationary spectrum $\varphi(\mathbf r)$6. This method was introduced to address limitations of Fourier-based methods such as screen periodicity and low spatial frequency power content, and AO simulator comparisons revealed significantly elevated residual closed-loop temporal power for small increases in added stochastic content at each time step (Srinath et al., 2015).

An alternative is the ergodic three-dimensional Karhunen–Loève construction. Instead of a single large frozen screen, one constructs a realization of $\varphi(\mathbf r)$7 inside a 3D sphere of radius $\varphi(\mathbf r)$8, defines a homogeneous and isotropic covariance

$\varphi(\mathbf r)$9

and solves the Mercer eigenproblem

$$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$0

The resulting random field is expanded as

$$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$1

with iid standard Gaussian $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$2. Two-dimensional phase screens are then obtained as planar cuts through the 3D volume, displaced by $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$3. Because the field is statistically stationary and isotropic in 3D, any planar slice—of any orientation or offset—automatically has the correct 2D von Kármán statistics, and moving the plane at constant speed $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$4 makes successive frames obey Taylor’s frozen-flow hypothesis in the mean. For implementation, Mathar proposes expansion in 3D Zernike functions to exploit spherical symmetry and reduce the eigenproblem to smaller dense blocks. The real-time slicing cost is stated as $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$5 for $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$6 frames, $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$7 pixels, and $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$8 modes, with $$D_\varphi(\Delta\mathbf r) = \bigl\langle\bigl|\varphi(\mathbf r+\Delta\mathbf r)-\varphi(\mathbf r)\bigr|^2\bigr\rangle,$$9 and $5/3$0 described as manageable on modern workstations (Mathar, 14 Oct 2025).

5. Finite-thickness and multi-wavelength generalizations

The standard phase-screen formalism usually assumes that the phase structure function is proportional to the path length through a turbulent layer. For a von Kármán refractive-index spectrum, Mathar reworked the integral exactly and found that the linear dependence on layer thickness $5/3$1 is only approximate. In the canonical form,

$5/3$2

the multiplicative factor $5/3$3 tends to $5/3$4 for $5/3$5, recovering the standard $5/3$6 law, but for $5/3$7 the asymptotics soften to

$5/3$8

The correction is more important for large than for small outer scales, and the numerical examples in the paper show reductions from about $5/3$9 to $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$0 relative to the linear-$D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$1 model, depending on $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$2, $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$3, and $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$4. The stated practical implication is that in tomographic AO or long-baseline interferometry, where one models turbulence in narrow altitude slices and separations up to tens of metres, omitting the correction leads to an overestimate of wavefront variance and an underestimate of the coherence radius $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$5 recovered from phase-screen statistics (Mathar, 2011).

Multi-wavelength phase screens extend the same logic to a set of coupled optical wavenumbers $D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$6. Hyde IV et al. use the modified von Kármán refractive-index spectrum

$D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$7

with single-wavelength phase spectrum

$D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$8

and two-wavelength cross-spectrum

$D_\varphi(\Delta r)=2.91\,(\Delta r/r_0)^{5/3}.$9

At each spectral point, a $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$00 correlation matrix is formed from normalized cross-spectra, Cholesky-factorized, and used to generate correlated Gaussian Fourier coefficients before inverse FFT. The method is validated by comparing the theoretical two-wavelength optical-path-length structure function to simulated results, which are reported to be in excellent agreement (IV et al., 1 May 2025).

6. Uses, caveats, and recurrent misconceptions

Von Kármán phase screens are applied across several simulation regimes: ground-based astronomy under frozen flow, long-range laser propagation, adaptive-optics telemetry fitting, beam-wander modeling over multiple screens, tomographic AO, long-baseline interferometry, and multi-wavelength wave-optics propagation (Mathar, 14 Oct 2025, Paulson et al., 2019, Srinath et al., 2015, IV et al., 1 May 2025). Their common role is to provide statistically controlled phase perturbations while accommodating finite outer scale and, in some formulations, inner scale, layer thickness, and wavelength coupling.

Several recurring misconceptions are corrected in the cited literature. One is that frozen translation of a static screen is the only principled way to generate time dependence; the AR method and the ergodic 3D-KL construction show that time evolution can instead be generated by controlled stochastic refresh or by slicing a stationary 3D random field (Srinath et al., 2015, Mathar, 14 Oct 2025). A second is that low-frequency correction alone is sufficient for FFT-based generators; the generalized autocorrelation-matrix method explicitly treats both low- and high-frequency errors, while RSS identifies small-scale accuracy issues for non-Kolmogorov power laws and proposes a white-noise remedy for spectral content outside the sampled band (Chhabra et al., 2021, Paulson et al., 2019). A third is that all “von Kármán screens” are statistically equivalent once they share the same PSD; comparison studies show that discretization strategy materially affects bias, computational cost, periodicity, and the reproduction of derived observables such as the structure function or angle-of-arrival variance (Charnotskii, 2019, Paulson et al., 2019).

The literature also distinguishes bounded and unbounded spectra. For a bounded von Kármán model with finite $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$01, core RSS is often sufficient; for spectra diverging at $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$02, such as pure Kolmogorov $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$03 with no inner or outer bounds, pure RSS may underrepresent the structure function at large scales and requires subharmonics. RSS further notes that absolute phase distribution can become log-normal for $$W_\phi(k)=0.023\,r_0^{-5/3}\,\bigl(k^2+L_0^{-2}\bigr)^{-11/6},$$04 domain, whereas relative-phase statistics remain Gaussian (Paulson et al., 2019). This suggests that, even within a single nominal turbulence law, the screen generator, the finite computational domain, and the observable of interest together determine whether a simulation is merely convenient or statistically faithful.