1. Introduction
Atmospheric turbulence is typically classified by its intensity into strong turbulence, weak turbulence, and moderate turbulence. Strong turbulence refers to turbulence with high intensity accompanied by significant fluctuations. In regions of strong turbulence, a beam experiences significant phase and logarithmic amplitude fluctuations induced by atmospheric disturbances. The turbulent atmosphere fragments the beam’s wavefront into numerous small speckles, creating a discontinuous phase surface marked by numerous phase singularities. This phenomenon leads to pronounced fluctuations in light intensity, causing flickering. During such conditions, the wavefront sensor struggles to accurately measure the complete distribution of the wavefront phase, resulting in a notable reduction in the efficacy of conjugate phase correction [
1,
2]. In contrast, wavefront sensor-less adaptive optics do not require a wavefront sensor to measure the distortion information of the signal light. Instead, it utilizes the control signal as an optimization parameter. Furthermore, parameters, such as imaging clarity, received optical energy, and other system performance indicators, can directly serve as objective functions for optimization algorithms. Leveraging these algorithms can yield nearly optimal correction outcomes. Hence, this technology is suitable for rectifying wavefront aberrations in free-space laser communication systems operating amidst strong turbulence conditions [
3,
4]. In a free-space optical communication system, when the selected signal wavelength has poor penetration or severe attenuation in the atmospheric channel, the system employs beacon light with a wavelength different from signal light, in conjunction with the Acquisition, Pointing, and Tracking (APT) system, to achieve capture, aiming, and tracking, thereby establishing the communication link [
5,
6]. However, because of the different wavelengths of signal light and beacon light, atmospheric dispersion and beam diffraction effects may lead to differences in wavefront phase aberrations after atmospheric turbulence transmission [
7,
8]. When the correction system adjusts the signal light by measuring the beam quality of the beacon light emitted from the same optical path as the signal light, such as the Strehl Ratio (SR), an optimization parameter for wavefront correction algorithms, a divergence emerges between the beam SR value captured by the CCD camera, which measures the beacon light, and the actual SR value of the signal light. This discrepancy ultimately leads to a certain degree of residual aberration within the transmission even after the correction system has been applied.
In recent years, numerous studies have delved into the disparities in wavefront aberrations between signaling light and beacon light within dual-wavelength free-space optical communication systems. In 2007, Li discovered that the wavelength of a beam, alongside atmospheric turbulence and receiving aperture characteristics, influences the performance of adaptive optics systems [
9]. Subsequently, in 2008, Nicholas identified variations in wavefront aberrations among different wavelengths transmitted through atmospheric turbulence, underscoring their significant impact on the correction of higher-order aberrations by adaptive optics systems [
10]. In 2015, Gorelaya conducted an indoor space laser communication experimental link featuring the co-optical transmission of signal and beacon lights. Using a Shack–Hartmann wavefront sensor at the receiving end, Gorelaya measured the wavefront phase aberrations of optical signals at wavelengths of 530 nm and 1060 nm [
11]. Additionally, in 2020, Xu scrutinized the polarization chromatic aberration between 1300 nm and 589 nm beams arising from the surface coating of the optical original in the Cassegrain optical antenna. The findings revealed that this type of polarization chromatic aberration encompasses out-of-focus aberration and a small amount of spherical aberration, with the out-of-focus aberration potentially reaching up to 1.14 rad [
12]. Further contributing to the understanding of wavefront variances, Ke derived a formula for the overall undulation variance in the wavefront for beams of different wavelengths, highlighting differences in wavefront variance under identical transmission conditions [
13]. Finally, in 2023, Wu discussed the wavefront distortion differences among beams induced by wavelength disparities in weak turbulence regions [
14].
In general, wavefront aberration disparities in beams across weak turbulence, moderate turbulence, and strong turbulence regions are well-documented, albeit with less emphasis on the intricacies of the strong turbulence domain. Moreover, the challenge of effectively addressing aberration differences among beams of varying wavelengths in the strong turbulence region remains largely unresolved. In this article, we meticulously analyze the variances in wavefront phase aberrations within the received cross-section subsequent to the transmission of signal and beacon lights of disparate wavelengths through the strong turbulence region. Additionally, we propose a corrective coefficient method employing wavefront-free correction. This method relies on the relationships of Zernike polynomial coefficients corresponding to beams of different wavelengths. Ultimately, the efficacy of our proposed method is validated through simulation and experimentation.
2. Differential Analysis of Beam Wavefront Aberrations at Different Wavelengths in Regions of Strong Turbulence
Signal light with wavelength and beacon light with wavelength are Gaussian beams, which, after long-distance transmission through the atmospheric channel, arrive at the surface of the receiving antenna. This process can be approximated as plane wave incidence. In the region of strong turbulence, the phase of the beam is torn into a discontinuous scattering distribution across the receiving cross-section. According to the Markov approximation, the following holds [
15].
where
denotes the refractive index of the air at a point on the cross-section of the transmission path,
and
represent any position vector on the transmission path cross-section, respectively, and
denotes the correlation function between different points. Furthermore,
denotes the proportionality coefficient, where
and
represent the position vector on the receiving cross section, and its value in the strong turbulence region is tentatively set to 1 to guarantee uniform homogeneity.
is the shock response function that describes the response function to the effect of the turbulent phase, where z represents the transmission distance along the propagation path. The wavefront phase functions of the signal beam and the beacon beam at the receiving end can be approximated as response functions with phases modulated by turbulence. These can be expressed as follows [
15]:
where
and
denoted the wavefront phase functions for wavelengths
and
, respectively.
L denotes the length of the transmission link, and
denotes the position vector in the receiver cross-section. Furthermore,
denotes the shock response function;
and
denote the wave numbers corresponding to wavelengths
and
, respectively. According to Equations (2) and (3), it is evident that the phase function response to atmospheric turbulence at a certain point on the receiving cross-section is related to the shock response function
and the wave number
. Therefore, the wavefront phase aberration of beams with different wavelengths after turbulence will differ. This dispersion non-equal halo error between the phases of beams of different wavelengths is represented by
and can be expressed as:
Based on Equation (4), it can be observed that
is affected by the wavelength difference, transmission distance, distance between the two points in the received cross-section, and atmospheric turbulence. In regions of strong turbulence, the wavefront phase changes rapidly, which in turn leads to rapid and random changes in the aberration difference between the two wavefronts. These changes do not satisfy the conditions of a strictly smooth stochastic process [
16]. The average value can only be maintained constant over a fairly short period. To solve this problem, instead of studying the difference value
directly, we examine the rate of change in the difference value, denoted as
. We consider the relative rate of change to satisfy the conditions of a smooth stochastic process, and the values that change slowly can be described by the structure function
of the phase difference as follows:
where
denotes the position vector between two points on the receiving cross-section and
represents ensemble average. From Equation (5), it can be observed that the total amount of the phase structure function is composed of a stabilizing quantity and small portion that exhibits slight variations. The spatial distribution of these variations in the receiving cross-section can be expressed as follows:
where
denotes the mutual coherence function. The phase response function is defined as the response to a unit impulse that is uniformly distributed across the receiving surface. At this point, the impact of atmospheric turbulence on the optical phase front can be described by the distribution of refractive index variations, which establishes a relationship between the structure function of phase difference and structure function for the random distribution of refractive index. This relationship is denoted as follows:
where
denotes the phase structure function of atmospheric turbulence. Based on the relationship between the turbulent power spectrum and phase structure function, the following relationship for
can be obtained [
17] as follows:
where
represents the number of interrelationships between two points on the spot,
denotes the number of spatial waves, and
denotes the spatial power spectrum of the turbulence. The modified von-Kármán power spectrum is used here, and it can be expressed as:
In the above equation, the relationship between the parameter
and outer scale
of turbulence can be expressed as
, where
denotes a proportionality parameter. Typically, its value is
in the region of strong turbulence. Furthermore,
denotes the atmospheric refractive index structure constant, which characterizes the strength of the turbulence and is a constant in horizontal transport. When
, turbulence intensity is weak;
indicates strong turbulence intensity; and
corresponds to moderate turbulence intensity [
17].
denotes a parameter that is related to
and can be expressed as
, where
denotes the internal scale of turbulence. The expansion
in Equation (8) can be expressed using the Maclaurin series as [
17]:
By substituting Equations (9) and (10) into Equation (8) and integrating, Equation (8) becomes [
17]:
By substituting Equation (11) into Equation (7) and transforming, Equation (7) can be expressed as:
where function
is a type of confluent hypergeometric function. It can serve as an approximation for large parameters when the distance
between two points on the receiver’s cross-section satisfies the condition
. This approximation for large parameters is provided in [
17]:
where
denotes the gamma function.
When distance
between two points on the receiving cross-section satisfies the condition of
. The following small parameter approximation can be used [
17]
When the above-mentioned equation is combined with Equation (12), it can be observed that the beam wavefront aberration difference is not only related to the beam wavelength but also closely related to the inner and outer turbulence scales.
The wavefront sensor fails to capture the continuous phase distribution within the strong turbulence region, rendering the phase relationship unusable for correcting wavefront aberrations. Additionally, conventional geometric optics analysis proves insufficient for analyzing wavefront phase across different wavelengths. Consequently, we opted for the statistical optics method, also termed as the phase structure function method. In statistical optics, the coefficient average, also referred to as the ensemble average, can be considered equal to the time average of a single stochastic process under the preconditions of a generalized stochastic smooth process. Leveraging the definitions of the coefficient average and root mean square (RMS), we derive the RMS of the error between beams of varying wavelengths by squaring the mean square value of Equation (12), which can be expressed as:
It is assumed that the values of the outer scale of turbulence
correspond to 20 m, 30 m, 50 m, 60 m, or 80 m, the turbulence inner scale parameter
is 0.1 m, and the atmospheric refractive index structural parameter value is
. When signal light with a wavelength of 1550 nm and beacon light with a wavelength of 632.8 nm are transmitted in the atmospheric channel in a co-optical path with
L = 10 km, the variation curves of the difference RMS in the wavefront phase aberration between different points on the received cross-section are shown in
Figure 1.
Figure 1 illustrates the beam wavefront aberration difference at the receiving end cross-section with the augmentation of the outer scale of turbulence when the distance between the two points is held constant. Notably, beyond an outer scale of 40 m, the rate of change in wavefront aberration diminishes, suggesting a gradual reduction in the impact of the turbulence outer scale on beam modulation. Furthermore, when the outer scale of turbulence is fixed, an increase in the distance between the two points at the receiver does lead to a higher wavefront aberration at the receiver’s cross-section; however, this overall variation magnitude remains relatively small. This observation underscores that, while different distances between two points influence the wavefront aberration at the receiver’s cross-section to some extent, the effect is not substantial.
It is assumed that the turbulence inner scale
values range from 0.01 m to 0.15 m and from 0.2 m to 0.4 m, with a turbulence outer scale parameter of
= 40 m and an atmospheric refractive index structural constant value of
. When signal light with a wavelength of 1550 nm and beacon light with a wavelength of 632.8 nm are transmitted through the atmospheric channel on a co-optical path with
L = 10 km, the influence of the turbulence inner scale on the RMS difference in wavefront phase aberration across various points on the received cross-section is shown in
Figure 2. Upon comparing
Figure 2a,b, it becomes apparent that when the separation between the two points is smaller than the inner scale of turbulence, there is a discernible effect on the wavefront phase aberration difference. As the distance between the points increases and approaches the inner scale of turbulence, this effect gradually diminishes, and the curves begin to overlap significantly. Subsequently, when the distance between the points surpasses the inner scale of turbulence, the wavefront phase aberration differences stabilize and the curves’ slopes exhibit no significant changes. This outcome may be attributed to beam diffraction within small-scale turbulence. However, as the distance between the points extends beyond the inner scale of turbulence, the light undergoes refraction once more, leading to a diminished influence of small-scale turbulence on the beam.
It is assumed that the turbulence inner scale
is 0.10 m, turbulence outer scale
is 50 m, and atmospheric refractive index structure constant corresponds to
. Signal light with a wavelength of 1550 nm and beacon lights with wavelengths of 632.8 nm, 530 nm, 850 nm, and 950 nm traverse the atmospheric channel along a shared optical path spanning 10 km. Variations in wavefront phase aberration across different points on the received cross-section are influenced by these wavelength discrepancies, as depicted in
Figure 3. In
Figure 3a, within the receiving aperture cross-section, the distance between two points exhibits minimal impact on the wavefront aberration difference between beams of distinct wavelengths, consistently maintaining a stable level. Notably, when the beacon light wavelength is 1280 nm, the disparity between the signal light and beacon light wavefront aberration remains below 1, rendering it practically negligible.
Figure 3b highlights that as the wavelength disparity between the signal light and the beacon light diminishes, the wavefront aberration difference progressively decreases. When the wavelengths of the signal light and the beacon light align, the wavefront aberration difference reaches 0. Thus, only when the signal light and beacon light wavelengths closely match, the accuracy of adaptive optics for wavefront correction can be maximized. Consequently, in practical applications, the wavelength difference between beacon light and signal light should be considered a critical factor when selecting the beacon light wavelength.