Short-exposure turbulence effects on the modulation transfer function of annular aperture optical systems

Abstract. The results of recent analytical work performed to characterize the short-exposure blurring effects of path integrated optical turbulence for annular aperture systems is described using the modulation transfer function. Results presented here extend the circular case to the annular aperture problem involving a circularly symmetric central obscuration associated with telescope optics secondary mirrors. To enable this extension, adjustments to procedures that compute the tilt-variance and phase-tilt correlation effects are described, and the results are illustrated. A four-dimensionally parameterized function of turbulence strength, diffraction influence, angular frequency, and fractional central obscuration is given.


Introduction
The standard theory for modeling atmospheric turbulence impacts in systems performance models [1][2][3] has been Fried's 4 short-exposure atmospheric modulation transfer function (MTF). Reanalysis of Fried's approach has shown that when diffraction influences on the phase structure function are included, 5 the tilt-phase (TP) correlation term is numerically evaluated, 6 and path-dependent turbulence strength variations are considered. 7 However, like the original Fried work, only circular aperture geometries have been considered. This omits a significant class of telescope optics involving a central obscuration produced by reflector optics' circularly symmetric secondary mirror. The current paper summarizes efforts to close that gap.
The closest parallel works in the area of optical turbulence effects on annular aperture systems appear to be limited. 8,9 However, these primarily considered the modeling of phase perturbation terms and not the full short-exposure MTF.
The extension of the theory to incorporate this additional effect involves the introduction of a free parameter, here designated as C, the ratio of the inner obscuration diameter, D 1 , to the outer system aperture diameter, D 2 , such that 0 ≤ C < 1. The limit C ¼ 0 corresponds to the circular solution.
The revised theory is addressed in general in Sec. 2. Thereafter, individual components of the calculation are considered separately, including the revised system MTF, the tilt-variance calculation, and the TP correlation calculation technique.

Short-Exposure Turbulence Modulation
Transfer Function The short-exposure turbulence MTF (SEMTF), based on Fried's approach, evaluates the impact of statistically averaged phase and amplitude fluctuations on image quality produced by an optical system. The statistics of the turbulence field is characterized by the wave and phase structure functions, denoted by DðrÞ and D ϕ ðrÞ, respectively. To generate averaged tilt-removed effects, the angle-of-arrival linear phase is mathematically removed. The linear tilt effects are deemed to shift the position of the centroid of the arriving pattern of photons but to not affect image quality. Only the remaining higher order phase and amplitude perturbations produce blur, resulting in MTF degradation.
The engine for the calculations performed involves the Rytov approximation's (RA) single-scattering model. This method was compared to the reportedly more robust method of Charnotskii 10 in Ref. 6. That paper also discussed limits of the RA based on Dashen's criteria. 11 However, here, for completeness, this discussion will be continued in Sec. 6.

Dimensionless Variables Used
The general optical problem encountered is illustrated in Fig. 1. An extended source object plane to be imaged is positioned at z ¼ −L. The signal emitted is an incoherent light source consisting of uncorrelated photons, some of which pass through the receiver aperture entrance pupil at z ¼ 0. The aperture consists of an annular ring whose maximum diameter is D 2 and whose inner obscuration diameter is D 1 . The light studied is monochromatic with wavelength λ and wavenumber k ¼ 2π∕λ, and it passes through an atmosphere characterized by turbulence strength parameter C 2 n (m −2∕3 ) that in this paper is independent of path position.
From these variables, we can define a set of dimensionless variables that parameterize the optical problem. First, C ¼ D 1 ∕D 2 characterizes the aperture geometry. Second, Q ¼ S ¼ D 2 ∕ρ o will also be used. This model ignores innerand outer-scale of turbulence effects, whose impacts are minimal for incoherent imaging problems due to the mitigating factors of aperture averaging on the inner scale and aperture size (D 2 < L o ) on the outer scale. Impacts of pathvarying turbulence were considered in Ref. 7.
The final governing variable is the dimensionless angular frequency ω, expressed in radians of frequency per radian of angle: Given a system aperture of outer diameter D 2 and wavelength λ, the maximum angular frequency is denoted Ω o ¼ D 2 ∕λ. This is the upper limit of variable Ω ¼ Rf developed by Goodman 12 to aggregate the effects of the system focal length R and spatial frequency f in the system's image plane. Here, Goodman's Ω∕Ω 0 is further reduced to the dimensionless frequency ω ¼ Ω∕Ω o , with limits 0 ≤ ω ≤ 1.

Wave and Phase Structure Functions
Let us next consider the form of the integrated wave and phase structure functions, given in terms of the above dimensionless parameters that describe the aggregate atmospheric turbulence properties. These functions were expressed in integral form for spherical wave propagation by Lawrence and Strohbehn 13 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 6 3 ; 3 3 5 DðvÞ ¼ D l ðvÞ þ D ϕ ðvÞ for the wave structure function, and E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 6 3 ; 2 6 3 ; (2) for the phase structure function. Here, l is a log-amplitude perturbation and ϕ is a phase perturbation at the system entrance pupil. J 0 ðxÞ is the Bessel function of the first kind, order zero, and c ¼ ðL þ zÞ∕L is a dimensionless path position, equal to zero at the object plane (z ¼ −L) and unity at the system aperture (z ¼ 0). The structure functions represent expectations of squared differences, e.g., D ϕ ðvÞ ¼ h½ϕðxÞ − ϕðx þ vÞ 2 i. The factor Φ n ðκÞ is the refractive index power spectrum. In this analysis, a Kolmogorov spectrum is considered such that Φ n ðκÞ ¼ βC 2 n ðzÞκ −11∕3 , where β ≈ 0.033 and κ (m −1 ) is the spatial frequency of atmospheric turbulence fluctuations.
Use of the Kolmogorov spectrum implies that the turbulence structure is homogeneous and isotropic. Note that the Bessel function terms involving J 0 produce a low frequency cutoff that further restricts outer scale influences.
The v variable in Eqs. (1) and (2) denotes a distance of spatial separation in a transverse plane between two observation points. In its truest sense, this variable must be treated as a vector quantity, but due to the choice of the isotropic/ homogeneous Kolmogorov spectrum it loses its directional dependence. Hereafter, v will be a transverse distance measured in the system entrance pupil. In this case, it is also common to normalize v as a dimensionless variable (e.g., u ¼ v∕D 2 ) such that the above integrals can be reduced to the forms: 5 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 3 2 6 ; 5 9 8 and E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 3 2 6 ; 5 3 7 Here, αðxÞ is a transitional function that varies between 1/2 at x ¼ 10 −8 , ∼3∕4 at x ¼ 10 −1 , and ∼1 beyond about x ¼ 10 (e.g., Fig. 1 of Ref. 5). In the numerical processing performed, this function is tabulated and interpolated as necessary. Also, JðxÞ encapsulates the behavior of the wave structure function, while YðxÞ separates the u dependence of the phase structure function from the S dependence. This separation will be critical in Sec. 4's TP analysis.
The final form on the right side of Eq. (4) emphasizes that the turbulence-related S factor's influence has been decoupled from the u dependence in YðQuÞ. This will allow turbulence-related effects to be factored out of the integral in Eq. (50) of Sec. 4.

Short-Exposure Turbulence Modulation Transfer Function Decomposition
Given these preliminaries, the SEMTF can now be evaluated. In its compact form, the integral expression defining this function is given as 4,5 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 ; 3 2 6 ; 2 7 8 is the area of the system entrance pupil, WðvÞ are window functions over this pupil, considered as binary functions (0 or 1), lðvÞ is a propagated amplitude perturbation, and ψðvÞ is a tilt-corrected phase perturbation. The wavefront tilt over the aperture is characterized by a, such that the tilt-corrected phase function in the system aperture can be expressed as ψðvÞ ¼ ϕðvÞ − a · v. Fried 4 showed in detail that l and ψ are independent random variables. Angle brackets are used to denote the expectation operator used to create the statistics for the net mean SEMTF.
The factor A imposes a normalization on the SEMTF integral: At zero frequency, ω ¼ 0, we find that This result also necessarily assumes hexpð2lÞi ¼ 1 [see, e.g., Equation (20) of Ref. 5. The principal difference between the short-and long-exposure imaging problems is the presence of ψ that replaces ϕ from the long-exposure case. The key element of the analysis of the SEMTF then reduces to the analysis of the expectation functions in Eq. (5). First, note that due to turbulence homogeneity, the relative starting point of these functions is indeterminate. Hence E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 7 ; 6 3 ; 6 1 9 hexpflðvÞ þ lðv − D 2 ωÞgi ≡ hexpflð0Þ þ lð−D 2 ωÞgi: Also, because the turbulence is isotropic, the resulting functions will only depend on magnitude Dω. We further define variables χðuÞ ¼ lð0Þ þ lðD 2 uÞ and ζðuÞ ¼ ψð0Þ − ψðD 2 uÞ. Prior analysis 4 showed that Eq.
Variables χ and ζ are both Gaussian random variates that are independent of one another. Variable ζ is also a zero-mean quantity.
The annular aperture window function WðvÞ can be constructed from a pair of functions W o ðv∕DÞ ¼ Cylðv∕DÞ as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 9 ; 6 3 ; 4 2 2 Wðv∕D 2 ; where CylðrÞ ¼ 1, if r < 1∕2 and zero when r ≥ 1∕2. The open portion of the aperture integrates to A ¼ Seemingly, the remaining step is to numerically evaluate the integrals with D l ðvÞ and D ψ ðvÞ. However, upon decomposing D ψ ðD 2 ωÞ, we find three terms E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 0 ; 6 3 ; 3 1 4 To assess these components, note that the first term on the right is just the phase structure function, D ϕ ðD 2 ωÞ. Combining this term with the amplitude structure function of Eq. (8) yields a wave structure function. This combined term represents the long-exposure atmospheric turbulence MTF kernel. The remaining terms correspond to the shortexposure corrections. The first of these extra terms (center term on the right) is due to the tilt-variance. However, this term is positive in Eq. (10), which turns negative in Eq. (8), resulting in an overall degradation influence. Therefore, it is the final term on the right in Eq. (10), designated here as the TP correlation, that supplies the counteracting short-exposure correction to long-exposure turbulent blur. The remaining task involves evaluating these remaining two terms, as well as assessing the net system MTF, defined as M 0 ðωÞ. To facilitate these computations, let us introduce a centered geometry E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 1 ; 3 2 6 ; 7 3 0 where TV is the tilt-variance term and TP is the tilt-phase correlation effect.
In the limit of negligible turbulence, the exponential term in Eq. (11) approaches unity, leaving the system MTF, M 0 ðω; CÞ, involving the integral of two shifted Eq. (9) window functions. The shifted overlapping windows are illustrated in Fig. 2. The frequency ω becomes the shift magnitude in the normalized aperture picture (D 2 becomes 1/2).
Function M 0 appears in integral form as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 2 ; 3 2 6 ; 5 7 1 To evaluate this integral, one must simply evaluate the areas of the two outer overlapping circles (the sum of areas S 1 through S 10 ) plus the overlapping inner circles (areas S 3 and S 4 ), then subtract off twice the overlap between an inner and an outer circle (e.g., areas S 1 through S 4 ). To evaluate these areas, consider a circle of radius R centered on the origin and integrated from Let us write this result as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 3 ; 3 2 6 ; 4 2 3Ũ where −R ≤ p ≤ þR. Using this result, for the overlapping outer circles, the net area is Fig. 2 Intersection of annular apertures displaced center-to-center from one another by amount ω of unit outer diameter and inner diameter C ¼ 9∕14. Subregions are labeled as S 1 through S 10 . Cross marks center of coordinate system (origin of v vector space).
Downloaded From: https://www.spiedigitallibrary.org/journals/Optical-Engineering on 17 Apr 2022 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 4 ; 6 3 ; 7 5 2 S 1 þ : : : þ S 10 ¼ 2ð1∕2Þ 2 U½ðω∕2Þ∕ð1∕2Þ ¼ UðωÞ∕2: The corresponding overlapping inner circle area is E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 5 ; 6 3 ; 7 0 9 S 3 þ S 4 ¼ 2ðC∕2Þ 2 U½ðω∕2Þ∕ðC∕2Þ ¼ C 2 Uðω∕CÞ∕2; (15) (assuming an overlap exists). The overlap between inner and outer circles is also evaluated using the U function but requires that arguments l 1 and l 2 be determined. From Fig. 2, if rightward-directed arrows are considered positive valued, leftward arrows as negative valued, and ω as positive, then E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 6 ; 6 3 ; 5 9 2 l 1 þ l 2 ¼ ω: To evaluate l 1 and l 2 , we also note that right triangles A-B-D and A-C-D share the common side of length h. Side B-D is C∕2 long and side C-D is 1/2 long, leading to the following equations: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 7 ; 6 3 ; 5 1 7 l 2 1 þ h 2 ¼ ðC∕2Þ 2 and l 2 2 þ h 2 ¼ ð1∕2Þ 2 : Subtracting these equations and using symbol e ¼ ; t e m p : i n t r a l i n k -; e 0 1 8 ; 6 3 ; 4 6 3 l 2 2 − l 2 1 ¼ e: Equation (16) can be used to eliminate either l 1 or l 2 from Eq. (18), resulting in E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 9 ; 6 3 ; 4 0 9 These results are directly comparable to the approach of Figs. 9-4 of Ref. 14. These expressions are just the arguments of the related U functions needed to evaluate these areas. Doubling the size of a single inner-outer overlapped area to account for S 1 -S 4 plus S 3 -S 6 , produces the result E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 0 ; 6 3 ; 3 0 4 Cylinder functions were added to handle the windowing of the function components. For shifts ω∕2 < 1 − C, the doubled inner circles lie completely inside the outer circles, whose net contributions are 2πC 2 ∕4. Combining normalized area terms and dividing byÂ yields the complete system MTF function, written as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 1 ; 6 3 ; 1 6 4 M 0 ðω;CÞ¼G 1 ðω;CÞ Various versions of this function for different choices of parameter C are plotted in Fig. 3. While the basic function is normalized such that M 0 ð0; CÞ ¼ 1, the curves plotted in Fig. 3 have been multiplied by (1 − C 2 ) to preserve a sense of the relative entrance pupil in terms of its light-gathering capacity.
The next two sections cover the calculation of the tiltvariance (TV) term in Sec. 3 and the TP correlation term in Sec. 4.

Annular Aperture Tilt-Variance Calculations
In this section, we consider the TV term in Eq. (11). In general, the TV is found to be aperture geometry dependent, but because the turbulence is assumed isotropic and the apertures are circularly symmetric, the tilt statistics is also isotropic (ha · ai ¼ ha 2 x i þ ha 2 Components a x and a y are also decorrelated (ha x a y i ¼ 0). Therefore, we can write E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 2 ; 3 2 6 ; 3 2 8 It remains to evaluate ha · ai. To do so, accounting for annular aperture influences, we return to an analysis of Fried's 15 phase perturbation model. The instantaneous phase perturbation function ϕðvÞ can be decomposed into an infinite series of weighted orthonormal basis functions E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 3 ; 3 2 6 ; 2 0 8 This expression may also be dissected in the form E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 4 ; 3 2 6 ; 1 4 7 where the lower order terms are resolved and the higher order remainder, ℋ N , is unresolved. Fried ordered the F μ basis functions such that the first few would involve the piston (index 0), tip (1), and tilt (2) terms, such that E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 5 ; 6 3 ; 7 1 9 ϕðvÞ ¼ Φ 3 ðvÞ þ ℋ 3 ðvÞ; where Φ 3 ðvÞ ¼ p 0 F 0 ðvÞ þ p x F x ðvÞ þ p y F y ðvÞ: Also, Φ 0 ¼ 0 (no expansion terms), such that ℋ 0 ðvÞ ¼ ϕðvÞ. In general, the total unresolved instantaneous phase variance averaged over the system aperture can be evaluated via the integral E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 6 ; 6 3 ; 6 2 5 In incoherent imaging, once a photon appears on the system image plane, knowledge of the mean phase is lost. Thus, the lowest order measureable phase error is Δ 2 1 . Short-exposure imaging further reduces the phase error by eliminating tip and tilt effects since a centroid shift will not reduce pointsource image quality. Hence, we focus on Δ 2 3 . The difference, y Þ∕A, is proportional to the TV. Evaluating this difference involves calculating p x and p y coefficients based on integrations using the F j basis functions. Fried required that these follow the normalization rule E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 7 ; 6 3 ; 4 3 7 The piston, tip, and tilt terms are given by E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 8 ; 6 3 ; 3 8 2 F 0 ðvÞ ¼ E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 9 ; 6 3 ; 3 3 7 where W T ¼ 8∕½πD 4 2 ð1 − C 4 Þ 1∕2 has dimensions of inverse area. These factors yield orthonormal annular functions (see, e.g., Ref. 9) as opposed to Fried's circular aperture factors.
The expansion coefficients p μ are then evaluated using E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 0 ; 6 3 ; 2 6 6 p μ ¼ Z dxWðjxjÞϕðxÞF μ ðxÞ: Dimensionally speaking, the p μ 's are the lengths, the F μ 's are the inverse lengths, and their products are dimensionless. Fried's 4 a · v product was also dimensionless but used lengths for v and inverse lengths for a, although their analogs F μ and p μ , respectively, were of opposing dimensional senses. The tilt vector components a μ are therefore related to quantities p μ as a μ v μ ¼ p μ F μ ¼ ðp μ W T Þv m u, such that E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 1 ; 6 3 ; 1 4 6 ha · ai ¼ hW 2 To evaluate Δ 1 1 − Δ 2 3 , let us follow Fried's development. We first write this difference as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 2 ; 3 2 6 ; 7 5 2 We then evaluate the two forms from Eq. (32) using the p μ forms of Eqs. (30) in (26) E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 3 ; 3 2 6 ; 7 0 7 Now, let us introduce Fried's 15 "trick" involving the rewriting of Δ 2 0 in a form resembling the bracketed final section of Eq. (33). First, E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 4 ; 3 2 6 ; 5 3 0 The trick involves recognizing that the bracketed region only produces a result for its F 0 ðxÞ term in Eq. (34). All other components integrate to zero due to orthogonality. For the zeroth component, the integral evaluates to A 1∕2 , which cancels with F 0 ¼ A −1∕2 .
Symmetrizing Eq. (34) into the form, E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 5 ; 3 2 6 ; 3 6 7 Equations (33) Optical Engineering 015105-5 January 2017 • Vol. 56 (1) Tofsted: Short-exposure turbulence effects on the modulation transfer function. . . using E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 9 ; 6 3 ; 7 4 1 Terms related to s can now be collected and integrated separately. Due to circular symmetry of the entrance pupil, the resulting functions depend only on magnitude r ¼ jrj. (40) where the C parameter denotes the additional influence of the central obscuration and the F μ functions involving dimensional arguments v (m) have been converted to functions G μ featuring dimensionless spatial distances u ¼ r∕D 2 . At this point, given available space, just the results of further manipulations will be shown. Calculation details were discussed in Ref. 16.
For the circular aperture case, the tilt-variance can be written as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 4 1 ; 6 3 ; 4 7 9 ha · ai ¼ 2. Σ P ðx; A; B 1 ; B 2 ; C 1 ; C 2 ; DÞ By requiring B 1 C 1 ¼ B 2 C 2 , both the value of the function and its slope will be contiguous at x ¼ D. This function approximates the curve plotted in Fig. 2 of Ref. 5. The current form has been tailored for the range −2 <q < þ2. The GðQÞ function is parameterized such that it models the main diffraction effects on the tilt-variance while the relative strength of the effect has been aggregated in the leading constant.
The addition of the central obscuration modulates Eq. (41) to the following form: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 4 6 ; 3 2 6 ; 6 1 5 The third function,GðQ; CÞ, has been similarly formulated as GðQÞ so that it approaches a limiting value of unity for large Q but also reflects the influences of C. It can be expressed as a modulated version of the original GðQÞ function  16], which is consistent with the concept that tilt effects are dominated by low wavenumber turbulence, with wavelengths wider than the system aperture.

Annular Aperture Tilt-Phase Correlation
Calculations The remaining impact to be quantified from Eq. (11) is the influence of the TP correlation TP. But, in doing so, we recognize, based on Eqs. (24) and (46), that the TV term is independent of the integration variable v in Eq. (11), as is the wave structure function [first term in the exponential in Eq. (11)]. Therefore, these terms can be factored out of the integral. Conversely, the TP term is dependent on v, and we must write E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 4 9 ; 3 2 6 ; 2 2 7 where the aperture effect has been normalized such that the entrance pupil area can be written usingÃ ¼ πð1 − C 2 Þ∕4.
Integration variable v has also been replaced by u ¼ v∕D 2 .
The full transition to a dimensionless form for TP requires several additional steps. Introducing a ¼ W T p and the Eq. (30) definition of the components into Eq. (10), then using Eq. (10) in the remnants of D ψ in Eq. (8) (gives us a 1/2 factor), and applying isotropic and homogeneous conditions for the structure functions Introducing dimensionless u ¼ v∕D 2 and rewriting the two phase structure function elements based on Eq. (4)'s decomposition into J and Y functions E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 5 1 ; 6 3 ; 5 8 4 The resulting integral is now only a function of u, Q, and C.
We also use the normalizedW window function. The results of this integral can be tabulated and interpolated separately for a more general recombining with multiple values of the leading turbulence strength factor. The results are also independent of the vector nature of ω due to rotational aperture symmetry. Thus, the tabulated results need only be fivedimensional, facilitated by assuming ω ¼ ωx is oriented along the x-axis.
To evaluate the integral in Eq. (51), let us divide the window functionW into an inner and an outer portion using the form E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 5 2 ; 6 3 ; 3 9 1 F X ðQ; C; u; ωÞ ¼ such that E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 5 3 ; 6 3 ; 3 0 1 TP ¼ 32JðS∕QÞ πð1 − C 4 Þ ½F X ðQ; 1; u; ωÞ − F X ðQ; C; u; ωÞ: Finally, to evaluate M TP , we use the following N-factor adjustment E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 5 4 ; 6 3 ; 2 3 9 M TP ðωÞ ¼ expðNÞ where variable N ensures that the integral can be evaluated within the numerical precision of the machine in the presence of large S (strong turbulence). Use of F X means only two sets of tabulated results need to simultaneously reside in the computer memory to evaluate results for any given C value. (57) where the TV term in Eq. (11) has been expanded using Eq. (22) and ha · ai from Eq. (44). The remaining M TP effect must be evaluated numerically using Eq. (54). However, rather than leave the results only in this tabulated format, an integrated short-exposure correction term has been developed as a perturbation of the original solution suggested by Fried in treating the circular case where X ¼ 2.1S replaces S.
To tabulate these effects, function Vðω; X; Q; CÞ was evaluated at 103 ω values, 41 logarithmically increasing X values ranging from 0.1 to 10 3 , 9 logarithmically increasing Q values ranging from 1/16 to 16, and 12 C values ranging from 0.00 to 0.95, for a total of 456,084 calculations.
To parameterize the C variability, a parameterized C dependence is constructed for each parameter E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 6 0 ; 3 2 6 ; 1 8 8 UðCÞ The original circular aperture solution was handled as an 80-parameter model. Using the above expression, each of the 80 parameters would need to be characterized by six separate parameters, resulting in a 480-parameter model.
The forms taken by the V at different C values are plotted in a series of graphs below. Figure 4 illustrates samples from the original circular aperture data set. The two groups of lines denote Q ¼ 1∕16 and Q ¼ 16 behaviors at [essentially] the extents of low Q (high diffraction effects) and high Q (farfield) imaging conditions. (Telescope type systems typically operate in the range of Q ¼ 2: : : 4.) One finds that the behavior of V does not change significantly above Q ¼ 2.
The significant feature of the V curves at high Q values is that they exceed unity over a range of frequencies ω but fall below unity as ω → 1. By falling below 1.0 at high frequency, the resulting forms obey diffraction limits, as expected. Conversely, when V > 1 at intermediate frequencies, one sees a midrange frequency correction of long-exposure turbulence blurring. This recovery is limited to high Q system configurations.
Extending the results illustrated in Fig. 4 to noncircular apertures, Figs. 5 and 6 illustrate behaviors under conditions of moderate (C ¼ 0.50) and large (C ¼ 0.95) obstructions. As shown, the main impact of increasing central obscuration is to affect midrange spatial frequencies. V values at these frequencies become increasingly attenuated as C increases. Nonetheless, interestingly, the floor values of the functions (at ω ¼ 0 and ω ¼ 1 frequencies) remain unaffected.
To more directly compare these results, Fig. 7 shows V values from a moderately strong turbulence scenario (X ¼ 10) for a moderately large aperture (Q ¼ 4) over a range of C parameter values. Here, C varies from 0.00 to 0.95, but most of the dynamic effects (degradation of the MTF in the high C range) are not significant before C ¼ 0.30 and are only severe beyond C ¼ 0.50. Thus, the majority of systems operating at C ≤ 0.50 experience turbulence effects nearly identical to that of circular apertures, differing primarily in terms of the form of the system MTF.
To further examine the variability of the V function, let us next note from a comparison of Figs. 5 and 6 how similar the evolution of a given family of V curves is for a given Q and C combination as X (S) and ω are varied. The primary elements that appear to change dramatically are the height of the baseline value of V along with the height of the peak above the baseline. The former is termed V 0 ðQ; CÞ, and the latter is designated as V M ðQ; CÞ. Figures 8 and 9 illustrate these functional behaviors.
As illustrated in Fig. 8, there are minimal differences in V 0 as C varies. The primary driving factor appears to be due to diffraction, as parameterized by its Q variations.
The V M function appears to exhibit both variations in C and Q. As opposed to the V 0 function, the influence of the central obscuration can significantly modulate this height. Note, particularly, that increasing C has a degrading effect on the midrange frequency enhancement noted for the circular aperture case. Nonetheless, we find that for relatively    large apertures (Q > 2), there is in fact no value of C that will not produce the needed enhancement to result in V > 1 at some range of frequencies.
On the other hand, if Q is small, e.g., Q < 1∕16, there is no C value that can produce the frequency enhancement region. Also, we see again that for C < 0.3 and Q > 2 one finds a plateau region where the behavior of the annular system approximates that of the circular aperture system.

Applicability
Dashen's 11 analysis of the applicability of the RA is now reviewed. Dashen proposed two parameters, Φ and Θ, to describe the influences of diffraction and turbulence on the propagation environment. According to Dashen, if Θ ≫ Φ the RA should apply, where Φ represents the rms phase fluctuation using first-order geometric optics to measure turbulence effects and Θ represents the diffraction effects due to a characteristic scale size of the turbulence itself. These two parameters were studied in Ref. 6, resulting in the relation Using typical values, this condition is easily achieved (e.g., L o ≥ 0.5 m, C 2 n < 10 −12 m −2∕3 , L ≤ 10 km). Potvin 17 has further shown recently that RA is consistent with beam steering based on Fermat rays (the Eikonal) affected by large-scale turbulent perturbations.
Moreover, from this author's perspective, the RA is under-rated, often being conflated with the Born approximation. The alternative, the mutual coherence function (MCF) approach, measures the degree to which an advancing wavefront loses coherence. Incoherent light involves the propagation of individual decorrelated photons. It is difficult to imagine how any given photon could become decorrelated from itself. Hence, MCF does not appear to be examining the same problem, and it may be that the RA is closer to the correct approach for handling incoherent radiation.

Conclusions
The complete model of the annular aperture SEMTF involves 480 parameter values consistent with the 80-parameter circular aperture model expanded to the annular aperture through use of the polynomial construction of Eq. (61), which multiplies the total number of parameters by a factor of six. The full model is available upon request and approval through appropriate channels.
In light of the degrading behavior of the SEMTF as C increases, one wonders whether additional terms that would significantly improve the performance could be relatively easily introduced. Hufnagel 18 addressed higher order phase corrections, but only for circular apertures. Perhaps the full annular aperture problem could supply added correction capabilities at high turbulence levels that are not achievable for circular apertures.
The development and analysis process described have been designed to provide system modelers and developers a tool for simulating a broad spectrum of system conditions and behaviors. In general, the behavior of annular aperture systems do not deviate significantly from circular aperture systems until C > 0.4. Study results may also be of use to adaptive optics designers by providing a baseline of performance when no correction is applied.  Plots of difference between the maximum V ðω; X ; Q; CÞ values (captured from the X ¼ 1000 results) above the minimum V values (V 0 ) for various Q and C combinations.