"Hidden phase" in two-wavelength adaptive optics.

Two-wavelength adaptive optics (AO), where sensing and correcting (from a beacon) are performed at one wavelength λ B and compensation and observation (after transmission through the atmosphere) are performed at another λ T , has historically been analyzed and practiced assuming negligible irradiance fluctuations (i.e., weak scintillation). Under these conditions, the phase corrections measured at λ B are robust over a relatively large range of wavelengths, resulting in a negligible decrease in AO performance. In weak-to-moderate scintillation conditions, which result from distributed-volume atmospheric aberrations, the pupil-phase function becomes discontinuous, producing what Fried called the "hidden phase" because it is not sensed by traditional least-squares phase reconstructors or unwrappers. Neglecting the hidden phase has a significant negative impact on AO performance even with perfect least-squares phase compensation. To the authors' knowledge, the hidden phase has not been studied in the context of two-wavelength AO. In particular, how does the hidden phase sensed at λ B relate to the compensation (or observation) wavelength λ T ? If the hidden phase is highly correlated across λ B and λ T , like the least-squares phase, it is worth sensing and correcting; otherwise, it is not. Through a series of wave optics simulations, we find an approximate expression for the hidden-phase correlation coefficient as a function of λ B , λ T , and the scintillation strength. In contrast to the least-squares phase, we determine that the hidden phase (when present) is correlated over a small band of wavelengths centered on λ T . Over the range λ B ,λ T ∈[1,3]µm and in weak-to-moderate scintillation conditions (spherical-wave log-amplitude variance σ χ2∈[0.1,0.5]), we find the average hidden-phase correlation linewidth to be approximately 0.35 µm. Consequently, for |λ B -λ T | greater than this linewidth, including the hidden phase does not significantly improve AO performance over least-squares phase compensation.

Active-illumination extension to the Priest and Meier pBRDF.

This paper develops a 3D vector solution for the scattering of partially coherent laser-beam illumination from statistically rough surfaces. Such a solution enables a rigorous comparison to the well-known Priest and Meier polarimetric bidirectional reflectance distribution function (pBRDF) [Opt. Eng.41(5), 988 (2002)10.1117/1.1467360]. Overall, the comparison shows excellent agreement for the normalized spectral density and the degree of polarization. Based on this agreement, the 3D vector solution also enables an extension to the Priest and Meier pBRDF that accounts for the effects of active illumination. In particular, the 3D vector solution enables the development of a closed-form expression for the spectral degree of coherence. This expression provides a gauge for the average speckle size based on the spatial-coherence properties of the laser source. Such an extension is of broad interest to long-range applications that deal with speckle phenomena.

Pulse-quality metric for nonstationary partially coherent fields.

This paper generalizes a pulse-quality metric referred to as P2, i.e., the time analogue of Siegman's beam quality factor M2, to include pulsed (nonstationary) random fields of any state of coherence. The analysis begins with the derivation of a general P2 relation, which we then specialize to the important cases of coherent and Schell-model pulsed beams. As examples, we derive the P2 for two stochastic sources: (1) a cosine Gaussian-correlated Schell-model pulsed beam and (2) a nonuniformly correlated pulsed beam. For both of these sources, we generate (in simulation) random instances of each and compare the simulated (Monte Carlo) P2, i.e., computed directly from its definition, to the theoretical quantity. The agreement is excellent, thereby validating our P2 analysis.

Multi-Gaussian random variables for modeling optical phenomena.

A generalization of the classic Gaussian random variable to the family of multi-Gaussian (MG) random variables characterized by shape parameter M > 0, in addition to the mean and the standard deviation, is introduced. The probability density function (PDF) of the MG family members is an alternating series of Gaussian functions with suitably chosen heights and widths. In particular, for integer values of M, the series has a finite number of terms and leads to flattened profiles, while reducing to the classic Gaussian PDF for M = 1. For non-integer, positive values of M, a convergent infinite series of Gaussian functions is obtained that can be truncated in practical problems. For all M > 1, the MG PDF has flattened profiles, while for 0 < M < 1, the MG PDF has cusped profiles. Moreover, the multivariate extension of the MG random variable is obtained and the log-multi-Gaussian random variable is introduced. In order to illustrate the usefulness of these new random variables for optics, the application of MG random variables to the characterization of novel speckle fields is discussed, both theoretically and via numerical simulations.

Twisted space-frequency and space-time partially coherent beams.

We present partially coherent sources that are statistically twisted in the space-frequency and space-time domains. Beginning with the superposition rule for genuine partially coherent sources, we derive source plane expressions for the cross-spectral density (CSD) and mutual coherence functions (MCFs) for twisted space-frequency and space-time Gaussian Schell-model (GSM) beams. Using the Fresnel approximation to the free-space Green's function, we then paraxially propagate the CSD and MCF to any plane [Formula: see text]. We discuss the beams' behavior as they propagate, with particular emphasis on how the beam shape rotates or tumbles versus z. To validate our analysis, we simulate the generation and subsequent propagation of twisted space-frequency and space-time GSM beams. We compare the simulated moments to the corresponding theoretical predictions and find them to be in excellent agreement. Lastly, we describe how to physically synthesize twisted space-frequency and space-time partially coherent sources.

Stochastic complex transmittance screens for synthesizing general partially coherent sources.

We develop a method to synthesize any partially coherent source (PCS) with a genuine cross-spectral density (CSD) function using complex transmittance screens. Prior work concerning PCS synthesis with complex transmittance screens has focused on generating Schell-model (uniformly correlated) sources. Here, using the necessary and sufficient condition for a genuine CSD function, we derive an expression, in the form of a superposition integral, that produces stochastic complex screen realizations. The sample autocorrelation of the screens is equal to the complex correlation function of the desired PCS. We validate our work by generating, in simulation, three PCSs from the literature-none has ever been synthesized using stochastic screens before. Examining planar slices through the four-dimensional CSD functions, we find the simulated results to be in excellent agreement with theory, implying successful realization of all three PCSs. The technique presented herein adds to the existing literature concerning the generation of PCSs and can be physically implemented using a simple optical setup consisting of a laser, spatial light modulator, and spatial filter.

Generating electromagnetic nonuniformly correlated beams.

We develop a method to generate electromagnetic nonuniformly correlated (ENUC) sources from vector Gaussian Schell-model (GSM) beams. Having spatially varying correlation properties, ENUC sources are more difficult to synthesize than their Schell-model counterparts (which can be generated by filtering circular complex Gaussian random numbers) and, in past work, have only been realized using Cholesky decomposition-a computationally intensive procedure. Here we transform electromagnetic GSM field instances directly into ENUC instances, thereby avoiding computing Cholesky factors resulting in significant savings in time and computing resources. We validate our method by generating (via simulation) an ENUC beam with desired parameters. We find the simulated results to be in excellent agreement with the theoretical predictions. This new method for generating ENUC sources can be directly implemented on existing spatial-light-modulator-based vector beam generators and will be useful in applications where nonuniformly correlated beams have shown promise, e.g., free-space/underwater optical communications.

Generating dark and antidark beams using the genuine cross-spectral density function criterion.

In this work, we demonstrate how to generate dark and antidark beams-diffraction-free partially coherent sources-using the genuine cross-spectral density function criterion. These beams have been realized in prior work using the source's coherent-mode representation and by transforming a J0-Bessel correlated partially coherent source using a wavefront-folding interferometer. We generalize the traditional dark and antidark beams to produce higher-order sources, which have not been realized. We simulate the generation of these beams and compare the results to the corresponding theoretical predictions. The simulated results are found to be in excellent agreement with theory, thus validating our analysis. We discuss the pros and cons of our synthesis approach vis-à-vis the prior coherent modes work. Lastly, we conclude this paper with a brief summary, and a discussion of how to physically realize these beams and potential applications.

Partially coherent sources generated from the incoherent sum of fields containing random-width Bessel functions.

Using the criterion for a genuine cross-spectral density function, we demonstrate the realization of an Im-Bessel correlated source, which has only recently been achieved using the source's coherent-mode representation. In addition, with just a simple change, we create a whole new class of partially coherent sources that have not been realized. We simulate the generation of these sources and compare the results to theoretical predictions to validate our analysis. The partially coherent sources described herein can easily be synthesized using spatial light modulators, and the approach presented in this Letter can be used to design sources for optical trapping, optical tweezers, and other related applications.

Array tilt in the atmosphere and its effect on optical phased array performance.

This paper investigates atmospheric array tilt and its effect on target-in-the-loop optical phased array (OPA) performance. Assuming a direct-solve, piston-only-phase-compensation OPA, two expressions for the atmospheric array tilt variance are derived using Mellin transform techniques. The first-the "full" array tilt variance-is germane when the OPA is sensitive to atmospheric tilt and is shown to significantly impact OPA target-plane intensity. The second-the Zernike-tilt-removed array tilt variance-is relevant when a separate system compensates for atmospheric tilt (the more likely scenario) and is shown to negligibly affect OPA performance. To show how atmospheric array tilt errors affect target-plane intensity, moments of the far-zone (or focused) array intensity, as functions of the array tilt variance, are derived and discussed. Lastly, Monte Carlo simulation results are presented to validate the theoretical array tilt variance expressions.

Behavior of tiled-aperture arrays fed by vector partially coherent sources.

This paper explores the far-zone behavior of partially coherent arrays. We derive an expression for the far-zone spectral density valid for any array composed of circular elements and fed by fields with Schell-model cross-spectral density functions. This expression is written as the sum of convolution integrals, making it easy to physically interpret. We discuss this expression at length and present examples. Lastly, we validate our analysis by comparing Monte Carlo averages from wave-optics simulations with theory. We conclude this paper with a brief summary of the results and potential uses of our work.

Polychromatic wave-optics models for image-plane speckle. 1. Well-resolved objects.

Polychromatic laser light can reduce speckle noise in wavefront-sensing and imaging applications that use direct-detection schemes. To help quantify the achievable reduction in speckle, this paper investigates the accuracy and numerical efficiency of three separate wave-optics methods. Each method simulates the active illumination of extended objects with polychromatic laser light. In turn, this paper uses the Monte Carlo method, the depth-slicing method, and the spectral-slicing method, respectively, to simulate the laser-object interaction. The limitations and sampling requirements of all three methods are discussed. Further, the numerical efficiencies of the methods are compared over a range of conditions. The Monte Carlo method is found to be the most efficient, while spectral slicing is more efficient than depth slicing for well-resolved objects. Finally, Hu's theory is used to quantify method accuracy when possible (i.e., for well-resolved objects). In general, the theory compares favorably to the simulation methods.

Polychromatic wave-optics models for image-plane speckle. 2. Unresolved objects.

Polychromatic laser light can reduce speckle noise in many wavefront-sensing and imaging applications. To help quantify the achievable reduction in speckle noise, this study investigates the accuracy of three polychromatic wave-optics models under the specific conditions of an unresolved object. Because existing theory assumes a well-resolved object, laboratory experiments are used to evaluate model accuracy. The three models use Monte-Carlo averaging, depth slicing, and spectral slicing, respectively, to simulate the laser-object interaction. The experiments involve spoiling the temporal coherence of laser light via a fiber-based, electro-optic modulator. After the light scatters off of the rough object, speckle statistics are measured. The Monte-Carlo method is found to be highly inaccurate, while depth-slicing error peaks at 7.8% but is generally much lower in comparison. The spectral-slicing method is the most accurate, always producing results within the error bounds of the experiment.

Monte Carlo simulations of three-dimensional electromagnetic Gaussian Schell-model sources.

This article presents a method to simulate a three-dimensional (3D) electromagnetic Gaussian-Schell model (EGSM) source with desired characteristics. Using the complex screen method, originally developed for the synthesis of two-dimensional stochastic electromagnetic fields, a set of equations is derived which relate the desired 3D source characteristics to those of the statistics of the random complex screen. From these equations and the 3D EGSM source realizability conditions, a single criterion is derived, which when satisfied guarantees both the realizability and simulatability of the desired 3D EGSM source. Lastly, a 3D EGSM source, with specified properties, is simulated; the Monte Carlo simulation results are compared to the theoretical expressions to validate the method.

Partially coherent sources with circular coherence: comment.

In [Opt. Lett.42, 1512 (2017)OPLEDP0146-959210.1364/OL.42.001512], the authors present a new class of non-uniformly correlated sources with circular coherence. They also describe a basic experimental setup for synthesizing this class of sources, which uses the Van Cittert-Zernike theorem. Here, we present an alternative way to analyze these sources and a different way to generate them.

Modeling random screens for predefined electromagnetic Gaussian-Schell model sources.

In a previous paper [Opt. Express22, 31691 (2014)] two different wave optics methodologies (phase screen and complex screen) were introduced to generate electromagnetic Gaussian Schell-model sources. A numerical optimization approach based on theoretical realizability conditions was used to determine the screen parameters. In this work we describe a practical modeling approach for the two methodologies that employs a common numerical recipe for generating correlated Gaussian random sequences and establish exact relationships between the screen simulation parameters and the source parameters. Both methodologies are demonstrated in a wave-optics simulation framework for an example source. The two methodologies are found to have some differing features, for example, the phase screen method is more flexible than the complex screen in terms of the range of combinations of beam parameter values that can be modeled. This work supports numerical wave optics simulations or laboratory experiments involving electromagnetic Gaussian Schell-model sources.

Temporal coherence effects on target-based phasing of laser arrays.

This paper studies how temporal coherence (in particular, linewidth broadening introduced to suppress stimulated Brillouin scattering) affects target-based phasing of fiber laser arrays. A radio-frequency modulated array whose elements are fed by a broadband laser source phasing on a remote step target is theoretically analyzed. An expression for the detector plane irradiance, ultimately used to phase the array on the target, is derived and discussed in detail. Simulation results of a seven-element hexagonal array phasing on a distant step target with scattering surfaces separated by many coherence lengths are presented to validate the theoretical findings.

Two spatial light modulator system for laboratory simulation of random beam propagation in random media: comment.

In a recent paper [Appl. Opt.55, 1112 (2016)APOPAI0003-693510.1364/AO.55.001112], the authors present a procedure for synthesizing Gaussian phase screens to generate Schell-model sources. Their synthesis method is analyzed, revealing interesting phenomena not fully described by the authors.

Two-dimensional simulation of optical wave propagation through atmospheric turbulence.

A methodology for the two-dimensional simulation of optical wave propagation through atmospheric turbulence is presented. The derivations of common statistical field moments in two dimensions, required for performing and validating simulations, are presented and compared with their traditional three-dimensional counterparts. Wave optics simulations are performed to validate the two-dimensional moments and to demonstrate the utility of performing two-dimensional wave optics simulations so that the results may be scaled to those of computationally prohibitive 3D scenarios. Discussions of the benefits and limitations of two-dimensional atmospheric turbulence simulations are provided throughout.

LSPV+7, a branch-point-tolerant reconstructor for strong turbulence adaptive optics.

Optical wave propagation through long paths of extended turbulence presents unique challenges to adaptive optics (AO) systems. As scintillation and branch points develop in the beacon phase, challenges arise in accurately unwrapping the received wavefront and optimizing the reconstructed phase with respect to branch cut placement on a continuous facesheet deformable mirror. Several applications are currently restricted by these capability limits: laser communication, laser weapons, remote sensing, and ground-based astronomy. This paper presents a set of temporally evolving AO simulations comparing traditional least-squares reconstruction techniques to a complex-exponential reconstructor and several other reconstructors derived from the postprocessing congruence operation. The reconstructors' behavior in closed-loop operation is compared and discussed, providing several insights into the fundamental strengths and limitations of each reconstructor type. This research utilizes a self-referencing interferometer (SRI) as the high-order wavefront sensor, driving a traditional linear control law in conjunction with a cooperative point source beacon. The SRI model includes practical optical considerations and frame-by-frame fiber coupling effects to allow for realistic noise modeling. The "LSPV+7" reconstructor is shown to offer the best performance in terms of Strehl ratio and correction stability-outperforming the traditional least-squares reconstructed system by an average of 120% in the studied scenarios. Utilizing a continuous facesheet deformable mirror, these reconstructors offer significant AO performance improvements in strong turbulence applications without the need for segmented deformable mirrors.