Next Article in Journal
Selenium Content of Goose Breast Meat Depending on the Type of Heat Processing
Previous Article in Journal
A Deep Learning Approach to Predict Supply Chain Delivery Delay Risk Based on Macroeconomic Indicators: A Case Study in the Automotive Sector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Difference in Wavefront Distortion on Beams Caused by Wavelength Differences in the Strong Turbulence Region

1
Faculty of Automation and Information Engineering, Xi’an University of Technology, Xi’an 710048, China
2
Shaanxi Civil-Military Integration Key Laboratory of Intelligence Collaborative Networks, Xi’an 710126, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(11), 4692; https://doi.org/10.3390/app14114692
Submission received: 1 April 2024 / Revised: 12 May 2024 / Accepted: 27 May 2024 / Published: 29 May 2024
(This article belongs to the Section Optics and Lasers)

Abstract

:
In free-space optical communication, the transmission of signal light and beacon light of differing wavelengths through the same atmospheric channel encounters variations in how the atmospheric refractive index absorbs and scatters light. This leads to distinct degrees of wavefront aberrations between the signal and beacon lights. In this study, we employed statistical optics to derive wavefront phase structure functions for both signal and beacon lights under conditions of strong turbulence. We explored how wavefront distortion varies among beams of different wavelengths after propagation through such turbulent conditions. Our findings revealed that as the turbulence outer scale escalates, the difference in wavefront distortion between signal and beacon lights stabilizes after an initial increase, assuming constant wavelengths. Furthermore, we observed significant changes in the relative wavefront aberrations when the inner scale of turbulence surpasses the separation between two points on the receiving apertures. As the disparity in wavelength decreases, so does the difference in wavefront aberrations. Finally, we propose a method for correcting wavefront aberrations based on coefficients of Zernike polynomials corresponding to beams with different wavelengths. This approach is validated through simulation and experimentation, demonstrating an 11% enhancement in the signal-to-optical Strehl ratio and a 0.072 increase in spot energy after the addition of correction coefficients compared with before their inclusion. These results solidify the efficacy of our method in improving adaptive optics correction accuracy.

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].
n ( R 1 ) n ( R 2 ) = B n ( R 1 R 2 ) = δ ( z ) A ( r 1 r 2 )
where n denotes the refractive index of the air at a point on the cross-section of the transmission path, R 1 and R 2 represent any position vector on the transmission path cross-section, respectively, and B n denotes the correlation function between different points. Furthermore, A r 1 r 2 denotes the proportionality coefficient, where r 1 and r 2 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. δ z 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]:
ϕ λ 1 = k 1 0 L δ n ( r , z ) d z
ϕ λ 2 = k 2 0 L δ n ( r , z ) d z
where ϕ λ 1 and ϕ λ 2 denoted the wavefront phase functions for wavelengths λ 1 and λ 2 , respectively. L denotes the length of the transmission link, and r denotes the position vector in the receiver cross-section. Furthermore, δ n ( r , z ) denotes the shock response function; k 1 and k 2 denote the wave numbers corresponding to wavelengths λ 1 and λ 2 , 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 δ n ( r , z ) and the wave number k = 2 π / λ . 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:
Δ ϕ r = ϕ λ 1 ϕ λ 2 = k 1 0 L δ n ( r , z ) d z k 2 0 L δ n ( r , z ) d z
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 Δ ϕ r directly, we examine the rate of change in the difference value, denoted as Δ ϕ r Δ ϕ r , Δ t . 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 D Δ ϕ ρ of the phase difference as follows:
D Δ ϕ ρ = Δ ϕ ρ + r Δ ϕ r 2
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:
D Δ ϕ ρ = Δ ϕ ρ + r Δ ϕ r 2 = 2 B n ρ B n r
where B n 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:
D Δ ϕ r = k 1 0 L δ n r 1 , z d z k 2 0 L δ n r 2 , z d z 2 = k 1 k 2 2 0 L δ n r 1 , z d z 0 L δ n r 2 , z d z 2 = k 1 k 2 2 D n r
where D n r 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 D n r can be obtained [17] as follows:
D n ρ = 2 B n r + ρ B n r = 8 π 0 κ 2 Φ κ 1 sin k ρ k ρ d κ
where B n r 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:
Φ κ = 0.033 C n 2 exp κ 2 / κ m 2 κ 2 + κ 0 2 11 / 6 , 0 κ <
In the above equation, the relationship between the parameter κ 0 2 and outer scale L 0 of turbulence can be expressed as κ 0 = c 0 / L 0 , where c 0 denotes a proportionality parameter. Typically, its value is 8 π in the region of strong turbulence. Furthermore, C n 2 denotes the atmospheric refractive index structure constant, which characterizes the strength of the turbulence and is a constant in horizontal transport. When C n 2 < 10 17   m 2 / 3 , turbulence intensity is weak; C n 2 > 10 13   m 2 / 3 indicates strong turbulence intensity; and 10 17   m 2 / 3 < C n 2 < 10 13   m 2 / 3 corresponds to moderate turbulence intensity [17]. κ m denotes a parameter that is related to l 0 and can be expressed as κ m = 5.92 / l 0 , where l 0 denotes the internal scale of turbulence. The expansion 1 sin κ ρ κ ρ in Equation (8) can be expressed using the Maclaurin series as [17]:
1 sin κ ρ κ ρ = n = 1 1 n 1 2 n + 1 ! κ ρ 2 n
By substituting Equations (9) and (10) into Equation (8) and integrating, Equation (8) becomes [17]:
D n ρ = 1.685 C n 2 κ m 2 / 3 F 1 1 1 3 ; 3 2 ; κ m 2 ρ 2 4 1 + 1.050 C n 2 κ 0 2 / 3 1 F 1 0 ; 2 3 ; κ 0 2 ρ 2 4 , κ 0 κ m
By substituting Equation (11) into Equation (7) and transforming, Equation (7) can be expressed as:
D Δ ϕ r = 1.685 C n 2 k 1 k 2 2 κ m 2 / 3 F 1 1 1 3 ; 3 2 ; κ m 2 ρ 2 4 1 + 1.050 C n 2 k 1 k 2 2 κ 0 2 / 3 1 F 1 1 1 ; 2 3 ; κ 0 2 ρ 2 4 , κ 0 κ m
where function F 1 1 a ; c ; x 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 l 0 ρ L 0 . This approximation for large parameters is provided in [17]:
F 1 1 a ; c ; x Γ c Γ c a x a
where Γ denotes the gamma function.
When distance ρ between two points on the receiving cross-section satisfies the condition of ρ l 0 . The following small parameter approximation can be used [17]
F 1 1 a ; c ; x 1 a c x
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:
RMS = D Δ ϕ r
It is assumed that the values of the outer scale of turbulence L 0 correspond to 20 m, 30 m, 50 m, 60 m, or 80 m, the turbulence inner scale parameter l 0 is 0.1 m, and the atmospheric refractive index structural parameter value is C n 2 = 1 × 10 12   m 2 / 3 . 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 l 0 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 L 0 = 40 m and an atmospheric refractive index structural constant value of C n 2 = 1 × 10 12   m 2 / 3 . 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 l 0 is 0.10 m, turbulence outer scale L 0 is 50 m, and atmospheric refractive index structure constant corresponds to C n 2 = 1 × 10 12   m 2 / 3 . 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.

3. Signal Optical Wavefront Aberration Correction for Dual-Wavelength Free-Space Optical Communication Systems

3.1. Theory

Zernike polynomials are a set of mutually orthogonal polynomials defined within a unit circle. Their coefficients, which are linear combinations of basis functions, can accurately describe different types of wavefront aberrations. Wavefront aberration, ϕ ρ , θ , can be expressed as a linear combination of Zernike polynomials [17]:
ϕ ρ = i = 1 N a i z i ρ , θ
where N denotes the total aberration order, a i denotes the coefficient of the Zernike polynomial of order i, ρ , θ denotes the polar coordinate, and Z i ρ , θ denotes the expression of the Zernike polynomial of order i in polar coordinates as follows:
Z i = n + 1 R n m ρ 2 cos m θ     m 0   and   m   i s   e v e n n + 1 R n m ρ 2 sin m θ     m 0   and   m   i s   o d d n + 1 R n m ρ                         m = 0
where
R n m ρ = s = 0 n m 2 1 s n s ! s ! n + m / 2 s ! n m / 2 s ! ρ n 2 s
where m and n denote the angular and radial numbers in polar coordinates of the Zernike polynomials, respectively. The essence of employing Zernike polynomials to model wavefront aberrations lies in determining both the order of the polynomials and the coefficients associated with each order. Essentially, higher-order Zernike polynomials encompass a broader spectrum of aberrations, thereby yielding greater fitting accuracy. Following the fitting of Zernike polynomials to the wavefront aberration phase for signal and beacon lights of varying wavelengths, discrepancies emerge in the coefficient matrices of the Zernike polynomials corresponding to the two beams within the same Zernike mode order.
The coefficients of the Zernike polynomials, which represent wavefront aberrations for signal light at wavelength λ 1 and for beacon light at wavelength λ 2 , can be expressed as follows [17]:
a 1 λ 1 , a 2 λ 1 , a 3 λ 1 , , a j λ 1 a 1 λ 2 , a 2 λ 2 , a 3 λ 2 , , a j λ 2
The numerical relationship between the coefficients of Zernike polynomials of the same order can be established through the wavelength function, which can be expressed as:
a 1 λ 1 = f 1 a 1 λ 2 , λ 1 a 2 λ 1 = f 2 a 2 λ 2 , λ 1     a j λ 1 = f j a j λ 2 , λ 1
Based on the actual transmission process of laser communication, we established a mode analysis model to determine the functional relationship among polynomial coefficients. We selected an infinitesimal region, denoted as Δ δ , on the wavefront phase plane and directed the signal light and beacon light into a common optical path. Assuming perfect alignment, the corresponding region Δ δ experiences the same transmission path through the atmospheric channel and the same variations in refractive index. Based on “Taylor’s” frozen-flow hypothesis [18], we selected a particular instant as the time for analysis. At this moment, the combined fluctuations in the atmospheric refractive index along the path of the beam can be approximated by a lens with an inhomogeneous refractive index distribution. This approximation can be analyzed using the formula for refractive index dispersion, which is commonly used for engineering materials. The relationship between the coefficients of the Zernike polynomials and wavelength can be expressed as follows [19]:
a i = A i + B i λ + C i λ 3.5
where a i denotes the coefficient of the Zernike polynomial of order i, and A i , B i , and C i denote the coefficients of the Conrady–Zernike polynomial. When the wavelengths differ, the relationship between their corresponding coefficients can be expressed as follows:
a i λ 1 = A i + B i λ 1 + C i λ 1 3.5 / A i + B i λ 2 + C i λ 2 3.5 a i λ 2
The proportionality correction factor between the ith order Zernike polynomial coefficients corresponding to the signal light and beacon light is as follows:
Z i = ( A i + B i λ 1 + C i λ 1 3.5 ) / ( A i + B i λ 2 + C i λ 2 3.5 )
Considering the 15th-order Zernike polynomial as an example, the Zernike polynomial coefficient vector corresponding to beacon light can be expressed as:
A 15 λ 2 = a 1 , a 2 , a 3 , , a 15 Τ
After multiplying the coefficient vector of the beacon light with the corresponding correction vector, we obtain:
A 15 λ 2 = z 1 , 0 , 0 , , 0 0 , z 2 , 0 , , 0 0 , 0 , 0 , , z 15 a 1 , a 2 , a 3 , , a 15 Τ = z 1 a 1 , z 2 a 2 , z 3 a 3 , , z 15 a 15 Τ
Furthermore, the mutual covariance matrix of the modified polynomial coefficients can be expressed as [20]:
C = E A 15 λ 2 A 15 Τ λ 2 = E z 1 a 1 , z 1 a 1 E z 1 a 1 , z 2 a 2 E z 1 a 1 , z N a N E z 2 a 2 , z 1 a 1 E z 2 a 2 , z 2 a 2 E z 2 a 2 , z N a N E z N a N , z 1 a 1 E z N a N , z 2 a 2 E z N a N , z N a N
The covariance matrix C is known as a Hermitian matrix [20]. Hence, there exists a unitary matrix X such that it is a diagonal array. We perform singular value decomposition of the covariance matrix C as follows:
C = X S X Τ
It is assumed that the coefficient vector of another linear combination B of wavefront Zernike polynomials of the same dimension as A is B = b 1 , b 2 , , b N Τ . Then, the Zernike polynomial coefficients of the wavefront of the corresponding signal optical distortion are:
A = X B

3.2. Simulation Verification

The wavefront correction system uses a charge-coupled device (CCD) camera to measure the spot energy, Strehl ratio, and other beam parameters of the laser beam, which are used as optimization indicators for the control program and converted into control signals for the deformable mirror [20,21]. The working principle of the system is shown in Figure 4, which mainly includes the following three parts: wavefront controller, wavefront corrector (deformable mirror, DM), and imaging detector (CCD). The parallel light generated by the light source forms distorted beams carrying aberrations after undergoing atmospheric turbulence transmission. This distorted beam is projected onto the deformable mirror, which corrects the distortion for the first time and reflects the remaining distorted beam onto the CCD. Next, the wavefront controller drives the intelligent algorithm to recalculate the control instructions of the deformable mirror based on the system performance indicators collected by the CCD. The deformable mirror corrects the wavefront distortion of the beacon light by adjusting its surface morphology. This process is cyclic, and after each round of correction, the deformable mirror further adjusts its shape based on the newly obtained feedback information of the light spot, gradually approaching the ideal correction state and achieving multiple closed-loop corrections of the distorted beam.
The control system adopts the stochastic parallel gradient descent (SPGD) algorithm. The algorithm first generates the basic voltage command for the deformable mirror and superimposes a random perturbation voltage signal that conforms to the Bernoulli distribution. On this basic voltage, forward and reverse voltage disturbances are applied separately, and then the change in the system objective function under each disturbance voltage is calculated. Finally, based on the calculation principles shown in Figure 5, the voltage signal that should be applied during the iteration cycle is determined. This process continues to loop until the system objective function reaches the preset optimization conditions.
In conditions characterized by an atmospheric coherence length r 0 = 0.01   m , a disturbance amplitude δ = 0.006 , a high gain coefficient of γ = 1.6 / ( 1 + J ) (where J denotes the value of the system’s objective function), and 100 iterations, numerical simulations were conducted to investigate the performance of dual-wavelength wavefront-free sensing adaptive optics technology in correcting far-field spot SR under atmospheric turbulence. The simulations employed the multi-phase screen method to model laser beam propagation through a turbulent atmosphere, with phase screens generated using Zernike polynomials [22]. The specific parameters for the simulation included a signal wavelength of 1550 nm, a beacon wavelength of 632.8 nm, a transmission distance of 10 km, a phase screen interval of Δ z = 1500 m, and a phase screen size of D = 0.4 m. The Zernike coefficient matrices for the beacon light, both pre- and post-correction, were integrated into the adaptive optics control system to rectify the wavefront aberration of the signal light. The resulting changes in SR for the signal light are depicted in Figure 6. Analysis of the figure reveals that over 100 iterations of correction, the SR value of the signal beam within the correction system lacking additional correction coefficients exhibits slow convergence, reaching a limit of 0.89, accompanied by notable jitter during the correction process. Conversely, the system augmented with correction coefficients realizes a higher SR value limit of 0.982, marking an 11% enhancement over the initial convergence target of 0.89. Furthermore, it attains convergence more rapidly, stabilizing at 70 iterations with significantly reduced jitter.

3.3. Experimental Research

The wavefront-free adaptive optics correction system is constructed indoors, with the experimental design principle illustrated in Figure 7. Both the signal light and the beacon light, emanating from a shared optical path, traverse through polarizer P1 before entering a 4F system comprising lenses L1 and L2 for beam expansion and collimation. Subsequently, the collimated beam proceeds to a spatial light modulator after being split by beam splitter BS1. Loaded with a phase screen, the spatial light modulator modulates the collimated light, and the modulated beam is reflected by mirror PM. Following reflection, the beam is again divided by BS2 and directed to the deformable mirror (DM). Here, the DM undertakes the initial correction of the distorted light, which is then converged through lens L3 and directed to the CCD. Furthermore, the CCD captures the intensity of the distorted beam and transmits this information to the computer. Utilizing optimization algorithms, the computer generates control signals for the deformable mirror, thereby executing multiple closed-loop corrections of the wavefront aberrations. The detailed equipment parameters employed in the experiment are outlined in Table 1.
In the experiments, we selected wavelengths of 1550 nm as the signal light and 632.8 nm as the beacon light, respectively. The turbulent phase screen based on the Zernike polynomial method was loaded into the liquid crystal spatial light modulator (LC-LSM) to cause wavefront distortion of the beam [23,24,25]. The SPGD algorithm based on the variable gain coefficient was used to correct the wavefront distortion. We used the average grey scale value of 10 × 10 pixel area centered on the spot’s center of mass, as captured by the CCD camera, to evaluate the calibration effect in terms of spot brightness. The size of the phase screen was D = 0.4 m. For the SPGD algorithm, we selected a perturbation amplitude value of δ = 0.006 and a gain coefficient value of γ = 1.6 / ( 1 + J ) . In order to simulate strong turbulence conditions, we set the atmospheric coherence length to r 0 = 0.01 m [26,27]. We conducted closed-loop correction experiments, and the results are shown in Figure 8a,b. By comparing the two figures, we can observe that the system without wavefront correction has a significant correction effect on the wavefront aberration in regions of strong turbulence. The SR value of the spot prior to correction stands at 0.112. Post-correction, it notably increases to 0.74, signifying a substantial enhancement in spot energy convergence. Upon employing the SPGD algorithm to generate the adaptive mirror voltage control signal, correction coefficients are incorporated to address the aberrations of the signal light. As depicted in Figure 8c, the results illustrate heightened spot energy compared with Figure 8b, where the dark spot is rectified. Notably, the beam’s SR value escalates to 0.812, underscoring the effectiveness of this approach in further refining the system’s correction accuracy without wavefront correction.

4. Conclusions

The discrepancy in wavefront distortion aberrations within an adaptive optical system represents a systematic error stemming from the disparity in wavefront lengths between the system’s signal and beacon lights. By adjusting system parameters and adopting advanced correction algorithms, the impact of turbulence on beam transmission can be effectively reduced. In our study, theoretical and experimental results show that when the turbulence external scale is greater than 40 m, the speed of wavefront aberration change decreases. Wavelength matching optimization can reduce phase difference, that is, the closer the wavelengths of the beacon light and signal light are, the closer the wavefront aberration is to 0. In numerical simulation, the introduction of correction coefficients in the wavefront correction system increased the SR value from the basic 0.89 to 0.982, by about 11%, and accelerated the convergence speed. In the experiment, the SR value of the spot before and after the implementation of wavefront correction technology jumped from 0.112 to 0.74. After further adding the correction factor, the SR value reached 0.812, an increase of about 10% compared with the initial correction, clearly demonstrating the additional enhancement effect of secondary correction on system accuracy. The further improvement in spot energy and the SR value verifies that this method can effectively improve correction accuracy.

Author Contributions

Conceptualization, M.H. and X.K.; methodology, M.H. and X.K.; software, M.H.; validation, M.H., X.K. and J.L.; formal analysis, M.H.; investigation, M.H.; resources, J.L.; data curation, M.H., X.K. and J.L.; writing—original draft preparation, M.H.; writing—review and editing, M.H. and X.K.; supervision, X.K.; project administration, M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Key Industrial Innovation Chain Project of Shaanxi Province [grant number 2017ZDCXLGY-06-01] and the General Project of the National Natural Science Foundation of China [grant number 61377080].

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the limitations of our patent work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mackey, R.; Dainty, C. Wavefront sensing and adaptive optics in strong turbulence. In Proceedings of the Opto-Ireland 2005: Photonic Engineering, Dublin, Ireland, 4–6 April 2005; Volume 5827, pp. 23–29. [Google Scholar]
  2. Wu, J.L.; Ke, X.Z. Development of adaptive optical correction and polarization control modules for 10-km free-space coherent optical communications. J. Mod. Opt. 2019, 67, 1–7. [Google Scholar]
  3. Yang, H.Z.; Cai, D.M.; Chen, B. Analysis of adaptive optics techniques without a wave-front sensor and its application in atmospheric laser communications. Chin. J. Lasers 2008, 35, 680–684. [Google Scholar] [CrossRef]
  4. Dong, B.; Yu, J. Hybrid approach used for extended image-based wavefront sensor-less adaptive optics. Chin. Opt. Lett. 2015, 13, 21–25. [Google Scholar] [CrossRef]
  5. Jiang, H.L.; Jiang, L.; Song, Y.S.; Meng, L.X.; Fu, Q.; Hu, Y.; Zhang, L.Z.; Yu, X.N. Research of Optical and APT Technology in One-Point to Multi-Point Simultaneous Space Laser Communication System. Chin. J. Lasers 2015, 42, 150–158. [Google Scholar]
  6. Ma, Y.F. Research on APT Spot Detection Algorithm in Space Optical Communication. Digit. Technol. Appl. 2017, 11, 115–116. [Google Scholar]
  7. Li, G.; He, Y.L.; Zeng, B. Spot positioning algorithm for laser radar based on triangulation. Chin. J. Sens. Actuators 2016, 29, 1692–1697. [Google Scholar]
  8. Zhang, J.Z.; Zhang, F.Z.; Li, Y.K.; Wu, Y. Theoretic analysis of expanded beacon’s anisoplanatic error. High Power Laser Part. Beams 2014, 26, 89–95. [Google Scholar]
  9. Li, F.; Jian, R.; Zhong, X.; Ding, X. Study on correcting angle-of-arrival fluctuations of space optical communication with AO. High Power Laser Part. Beams 2007, 6795, 141–146. [Google Scholar]
  10. Nicholas, D.; Alexander, G.; Christopher, D. Chromatic effects of the atmosphere on astronomical adaptive optics. Appl. Opt. 2008, 47, 1072–1081. [Google Scholar]
  11. Gorelaya, A.V.; Shubenkova, E.V.; Dmitriev, D.I.; Dmitrieva, A.D.; Kudryashov, A.V.; Lovchiy, I.L.; Shalymov, E.V.; Sheldakova, Y.V.; Tsvetkov, A.D.; Venediktov, D.V.; et al. Investigation of dual-wavelength laser beam propagation along the in-door atmospheric path. In Proceedings of the Conference on Optics in Atmospheric Propagation and Adaptive Systems, Toulouse, France, 22 September 2015. [Google Scholar]
  12. Xu, B.W.; Hu, S.J.; Li, J.; Fan, J.Q.; Tan, Y. Impact of Polarization Chromatic Aberration on Adaptive Optical Systems. Laser Optoelectron. Prog. 2020, 57, 44–51. [Google Scholar]
  13. Ke, X.Z.; Chen, X.Z. Correcting wavefront distortion of dual-wavelength beams due to atmospheric turbulence with a correction coefficient. Opt. Photonics J. 2020, 10, 64–77. [Google Scholar] [CrossRef]
  14. Wu, J.L.; Ke, X.Z.; Kang, W.L. Study on the difference of wavefront distortion on beams caused by wavelength-differences in weak turbulence region. Photonics 2023, 10, 725–739. [Google Scholar] [CrossRef]
  15. Norris, J.R. Markov Chains; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  16. Aleksandrovich, K.L.; Rytov, S.M. Elements of Random Process Theory; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  17. Andrews, L.C.; Phillips, R. Laser Beam Propagation through Random Media; SPIE Press: Bellingham, WA, USA, 2005. [Google Scholar]
  18. Zhang, A.C. Modern Meteorological Observation; Peking University Press: Beijing, China, 2000. [Google Scholar]
  19. Zhang, Q.; Han, S.; Tang, S.H. Study on functional relationship between transmitted wavefront Zernike coefficients and wavelengths. Acta Opt. Sin. 2018, 38, 0212002. [Google Scholar] [CrossRef]
  20. Weyrauch, T.; Vorontsov, M.A. Free-space laser communications with adaptive optics: Atmospheric compensation experiments. J. Opt. Fiber Commun. Rep. 2004, 1, 355–379. [Google Scholar] [CrossRef]
  21. Song, H.; Fraanje, R.; Schitter, G.; Kroese, H.; Vdovin, G.; Verhaegen, M. Model-based aberration correction in a closed-loop wavefront-sensor-less adaptive optics system. Opt. Express 2010, 18, 24070–24084. [Google Scholar] [CrossRef] [PubMed]
  22. Roddier, N.A. Atmospheric wavefront simulation using Zernike polynomials. Opt. Eng. 1990, 29, 1174–1180. [Google Scholar] [CrossRef]
  23. Phillips, J.D.; Goda, M.E.; Schmidt, J. Atmospheric turbulence simulation using liquid crystal spatial light modulators. In Advanced Wavefront Control: Methods, Devices, and Applications III; SPIE: Bellingham, WA, USA, 2005; pp. 57–67. [Google Scholar]
  24. Shan, X.; Liu, M.; Zhang, N.; Ai, Y. Laboratory simulation of non-Kolmogorov turbulence based on a liquid crystal spatial light modulator. Opt. Eng. 2017, 56, 026111. [Google Scholar] [CrossRef]
  25. Liu, Y.J.; Hu, L.F.; Cao, Z.L.; Li, D.Y.; Mu, Q.; Lu, X.H.; Xuan, L. Liquid crystal atmosphere turbulence simulator. Acta Photonica Sin. 2006, 35, 1960–1963. [Google Scholar]
  26. Ai, Y.; Duan, M.Y.; Xu, J.; Shan, X.; Chen, J.; Xiong, Z.; Jiang, R. LC-SLM laser atmospheric transmission turbulence simulation and communication test. Infrared Laser Eng. 2015, 44, 3103–3109. [Google Scholar]
  27. Liu, C.; Chen, S.; Li, X.; Xian, H. Performance evaluation of adaptive optics for atmospheric coherent laser communications. Opt. Express 2014, 22, 15554–15563. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The influence of the distance between two points on the wavefront phase aberration difference for beams with different wavelengths.
Figure 1. The influence of the distance between two points on the wavefront phase aberration difference for beams with different wavelengths.
Applsci 14 04692 g001
Figure 2. Influence of the turbulence inner scale on the difference in wavefront aberrations of beams with different wavelengths: (a) l0 = 0.01 m ~ l0 = 0.15 m and (b) l0 = 0.20 m ~ l0 = 0.40 m.
Figure 2. Influence of the turbulence inner scale on the difference in wavefront aberrations of beams with different wavelengths: (a) l0 = 0.01 m ~ l0 = 0.15 m and (b) l0 = 0.20 m ~ l0 = 0.40 m.
Applsci 14 04692 g002
Figure 3. The difference in beam wavefront phase distortion between different points on the receiving cross-section. (a) Different wavelengths. (b) Difference between different wavelengths.
Figure 3. The difference in beam wavefront phase distortion between different points on the receiving cross-section. (a) Different wavelengths. (b) Difference between different wavelengths.
Applsci 14 04692 g003
Figure 4. The working principle of the wavefront correction system.
Figure 4. The working principle of the wavefront correction system.
Applsci 14 04692 g004
Figure 5. Schematic diagram of the correction of wavefront aberration differences without wavefront correction.
Figure 5. Schematic diagram of the correction of wavefront aberration differences without wavefront correction.
Applsci 14 04692 g005
Figure 6. Signal light SR value variation curve with the number of iterations. (a) Without the wavefront correction coefficient. (b) With the wavefront correction coefficient.
Figure 6. Signal light SR value variation curve with the number of iterations. (a) Without the wavefront correction coefficient. (b) With the wavefront correction coefficient.
Applsci 14 04692 g006
Figure 7. Experimental schematic diagram of wavefront sensor-less adaptive optics.
Figure 7. Experimental schematic diagram of wavefront sensor-less adaptive optics.
Applsci 14 04692 g007
Figure 8. The change in the light pattern in the picture correction experiment. (a) Before correction; (b) correction without the correction factor, and (c) after adding the correction factor correction.
Figure 8. The change in the light pattern in the picture correction experiment. (a) Before correction; (b) correction without the correction factor, and (c) after adding the correction factor correction.
Applsci 14 04692 g008
Table 1. Parameters of experimental equipment.
Table 1. Parameters of experimental equipment.
Equipment NameTypeEquipment Parameter
Liquid crystal spatial light modulatorRL-SLM-R2Target size: 0.78″
Pixel size: 12.3 mm
Operating wavelength: 400~1620 nm
Phase modulation capability: 0~2 π @532 nm
Deformable mirrorALPAO-DMNumber of actuators: 69
Pupil diameter: 10.5 mm
Oblique wavefront modulation: 60 mm
Bandwidth: >750 Hz
Setting time: 800 μs
Analog response time: <10 μs
CCDMW-GX-650Pixel Size: 5.5 mm × 5.5 mm
Frame rate: 8 fps
Exposure time: 300–1,000,000 ms
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Han, M.; Ke, X.; Liang, J. Study on the Difference in Wavefront Distortion on Beams Caused by Wavelength Differences in the Strong Turbulence Region. Appl. Sci. 2024, 14, 4692. https://doi.org/10.3390/app14114692

AMA Style

Han M, Ke X, Liang J. Study on the Difference in Wavefront Distortion on Beams Caused by Wavelength Differences in the Strong Turbulence Region. Applied Sciences. 2024; 14(11):4692. https://doi.org/10.3390/app14114692

Chicago/Turabian Style

Han, Meimiao, Xizheng Ke, and Jingyuan Liang. 2024. "Study on the Difference in Wavefront Distortion on Beams Caused by Wavelength Differences in the Strong Turbulence Region" Applied Sciences 14, no. 11: 4692. https://doi.org/10.3390/app14114692

APA Style

Han, M., Ke, X., & Liang, J. (2024). Study on the Difference in Wavefront Distortion on Beams Caused by Wavelength Differences in the Strong Turbulence Region. Applied Sciences, 14(11), 4692. https://doi.org/10.3390/app14114692

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop