Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Stretchable Strain Sensor with Small but Sufficient Adhesion to Skin
Next Article in Special Issue
Revealing Local Temporal Profile of Laser Pulses of Intensity above 1014 W/cm2
Previous Article in Journal
Hardware-in-the-Loop Test Bench for Simulation of Catenary–Pantograph Interaction (CPI) with Linear Camera Measurement
Previous Article in Special Issue
Sensing Trace-Level Metal Elements in Water Using Chirped Femtosecond Laser Pulses in the Filamentation Regime
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Air Turbulence on Vortex Beams in Nonlinear Propagation

1
School of Electronics and Information Engineering, Tiangong University, Tianjin 300387, China
2
Tianjin Key Laboratory of Optoelectronic Detection Technology and Systems, Tianjin 300387, China
3
School of Physical Science and Technology, Tiangong University, Tianjin 300387, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(4), 1772; https://doi.org/10.3390/s23041772
Submission received: 14 January 2023 / Revised: 1 February 2023 / Accepted: 3 February 2023 / Published: 4 February 2023
(This article belongs to the Special Issue Sensing with Femtosecond Laser Filamentation)

Abstract

:
Vortex beams with orthogonality can be widely used in atmospheric applications. We numerically analyzed the statistical regularities of vortex beams propagating through a lens or an axicon with different series of turbulent air phase screens. The simulative results revealed that the distortion of the transverse intensity was sensitive to the location and the structure constant of the turbulence screen. In addition, the axicon can be regarded as a very useful optical device, since it can not only suppress the turbulence but also maintain a stable beam pattern. We further confirmed that a vortex beam with a large topological charge can suppress the influence of air turbulence. Our outcomes are valuable for many applications in the atmospheric air, especially for optical communication and remote sensing.

1. Introduction

Ultrafast laser filamentation has attracted more and more interest due to its potential applications, such as THz generation [1,2], pulse compression [3,4], air lasing [5,6] and remote sensing [7,8]. Ultrafast laser filamentation is a unique nonlinear optical phenomenon that occurs in transparent optical media when the laser reaches the critical power (Pcr ≈ 5 GW in air) [9,10,11]. It is mainly caused by the dynamic balance between the optical Kerr-effect-induced self-focusing and defocusing effects of the plasma generated by a high-intensity laser or higher-order Kerr effect (HOKE) [12,13]. Fruitful nonlinear optical phenomena can be observed during the ultrafast laser filamentation evolution, including self-focusing, group velocity dispersion (GVD), self-steepening, self-phase modulation (SPM), multi-photon and tunneling ionization (MPI/TI) and stimulated amplification [9], etc. Recently, several studies demonstrated that ultrafast laser filamentation was able to achieve free space optical (FSO) communication through clouds and fog in the atmosphere [14,15]. These reports open a new and exciting prospect for the atmospheric application of optical communication.
As is known, a vortex beam is considered to have many promising applications since its special helical wavefront can generate the orbital angular momentum (OAM). A vortex beam is characterized by the phase factor exp( i l θ ), where θ is the azimuthal coordinate and l, termed the topological charger, refers to an integer counting the number of intertwined helices, which represents the number of 2π phase shifts along the circle [16]. The mode of OAM state is equal to this number l. A vector normal to a helical phase front follows a spiral trajectory around the propagation axis. The different OAM modes of the vortex beam are orthogonal with each other [17]. More importantly, it is found that a vortex beam with orthogonality can be widely used in FSO communication and remote sensing systems [18,19]. In the past decade, several methods have been proposed by which vortex beams carrying OAM can intensively increase the information rate and capacity of FSO communication systems [20,21,22,23,24]. Although a vortex beam usually has a ring structure in transverse modulation, it still undergoes the process of self-focusing and defocusing effects when the peak intensity of the laser power is high enough to exceed the critical power. Therefore, a vortex beam carrying OAM during nonlinear propagation is also considered to be laser filamentation.
These types of atmospheric applications, such as optical communication [25,26,27,28], remote sensing and LIDAR technology (light detection and ranging) [29], etc., require the long-distance propagation of the laser beam through the atmosphere. Turbulent air is the main obstacle during the propagation of ultrafast filamentation through the atmosphere. It generates optical perturbation in the air density, resulting in fluctuations in the refractive index in the atmosphere and instabilities in the optical transverse modulation [25,26]. Some studies were devoted to the influence of turbulent air on the filament distance [30,31], optical pulse broadening [32] and transverse filament wandering [33,34], etc. Due to the existence of turbulence, the vortex beam will suffer from serious beam expansion, drift and distortion. Consequently, the code information that is carried by the vortex beam will be easily lost, and the blurring and scintillation that are generated by phase-front distortions will severely degrade the system’s stability [35,36]. Therefore, it is very significant to deeply understand the influence of turbulent air on the process of filamentation.
In this paper, we present the effect of air turbulence on the transverse beam wandering when vortex beams propagate through a lens or an axicon during nonlinear propagation. In order to obtain the statistical regularities of vortex beam propagation in turbulent air, we used 20 independent simulations with different series of stochastic turbulence phase screens for each condition. The statistical analyses found that beam wandering was affected by the location and the structure constant of the turbulence screens. Meanwhile, it was found that a nonlinear effect is better in suppressing the air turbulence during propagation. In addition, in our simulations, the axicon could be regarded as a very useful optical device in not only suppressing the turbulence, but also maintaining a stable beam pattern. We further confirmed that a large topological charge can reduce the influence of air turbulence during propagation. Our results are meaningful to remote sensing and optical communication in the atmosphere.

2. Numerical Simulation Model

The numerical simulation model of vortex beam propagation through air turbulence consists of two parts, including the filamentation model and the turbulent phase screen model. Our theoretical study mainly relies on the numerical simulations of a 2D + 1 [A(x,y,z)] nonlinear wave equation [33,37]:
2 i k 0 A z + ( 2 x 2 + 2 y 2 ) A + 2 k 0 2 Δ n A + 2 k 0 2 n ˜ ( x , y , z ) A = 0 ,
where A denotes the amplitude of a continuous wave beam, and k0 represents the wave number of the beam, whose central wavelength is 800 nm. In Equation (1), Δn corresponds to the intensity-dependent refractive index, which is determined by the optical Kerr effect ( Δ n ker = n 2 I , n2 = 2 × 10−19 cm2/W) and the plasma defocusing effect ( Δ n p l a s m a = σ I m )-induced nonlinear refractive index. Here, σ is an empirical parameter determined from our previous work, which gives rise to a clamped intensity of 5 × 1013 W/cm2 [38] in air. Term m is chosen to be equal to 8, which is approximately the reported effective nonlinearity order in air by a near-infrared femtosecond laser [39]. Since we focus mainly on the spatial distribution of the vortex beam, the temporal aspects of the nonlinear propagation are not considered in Equation (1). The validity of this type of simplified simulation model has been clarified by previous studies on multi-filamentation [40,41].
Term n ˜ ( x , y , z ) of Equation (1) corresponds to the fluctuations in the air refractive index caused by the air turbulence. Corresponding to the reality of optical perturbation in the atmosphere, we use the classical model of the modified Karman spectrum describing the stationary, isotropic and homogeneous turbulent air [34]:
ϕ n ( κ x , κ y , κ z ) = 0.033 C n 2 ( κ 2 + κ 0 2 ) 11 / 6 exp [ ( κ 2 / κ m 2 ) ] , 0 κ < ,
where Cn2 denotes the refractive index structure constant, whose values vary from 10−15 cm−2/3 in relative weak air turbulence to approximately 10−11 cm−2/3 in relatively strong turbulence. Term κ refers to the spatial wave number, which is given by κ = κ x 2 + κ y 2 + κ z 2 , while κ m = 5.92 / l 0 and κ 0 = 2 π / L 0 . Terms l 0 and L 0 are the inner and outer scales of turbulent air, respectively. The scale of the air turbulence air refractive index fluctuations varies from the inner scale l 0 , which is approximately 0.1~1 cm, to the outer scale L 0 , which can be several meters [42].
In practice, the typical air turbulence phase image consists of a number of superimposed turbulence spectral densities distributed over a distance. During the propagation, we divide the cumulative modified Karman spectrum into several parts, and each part is regarded as a complete “sub-phase screen”, which has a unified amplitude. The spectral density of phase fluctuations on each part Δ z has the form of
F ϕ ( κ x , κ y ) = 2 π k 0 2 Δ z ϕ n ( κ x , κ y , 0 ) .
According to the method described in previous studies [42,43], the spatial phase fluctuations are reconstructed by the fast Fourier transform (FFT) of the spatial spectrum density indicated by Equation (3). The random complex 2D field of air turbulence phase screen ϕ ( x , y ) is given by
ϕ ( x , y ) = 1 N M p = N 2 N 2 1 q = M 2 M 2 1 α p q ( ε p q + i η p q ) exp [ i 2 π ( p x N + q m M ) ]
where N and M are the numbers of the numerical simulation grids. Terms ε p q and η p q are statistically independent random numbers, and α p q is determined by the spectral density of phase fluctuations as follows:
α p q 2 = F ϕ ( p Δ k x , q Δ k y ) Δ k x Δ k y
where Δ k x and Δ k y are connected with the transverse dimensions of the phase screen Lx and Ly, respectively, as
Δ k x = 2 π L x , Δ k y = 2 π L y
Finally, the air refractive index fluctuations n ˜ ( x , y , z ) of Equation (1) can be expressed as follows:
n ˜ ( x , y , z ) = ϕ ( x , y ) k 0 Δ z .
We chose the numerical simulation grids (1200 × 1200) with a unit grid size of 5 μm for both the vortex beam propagation model and the turbulent phase screen model. Because of the random and stochastic characteristic, as an illustration, Figure 1 shows the representative simulated turbulent phase screen with the structure constant Cn2 = 5.6 × 10−13 cm−2/3 and the distance Δz = 1 m. The color scale reflects the values of phase fluctuations in radians. It is worth noting that there are several models used to generate turbulent phase screens. Currently, the most common methods for turbulent phase screen generation use an FFT of random spectral density with statistics matching the phase turbulence to a certain extent. The FFT phase is complemented by subharmonic terms that compensate for the lack of low-frequency FFT components. However, it is only capable of generating phase screens that are restricted to a rectangular spatial grid [44]. The validity of the turbulent phase screen model generated by the FFT method has been proven in our previous studies [42].

3. Results and Discussion

In our numerical simulation, we set topological charge l as 10, since it has been reported that a large topological charge can reduce the influence of air turbulence during propagation [45]. Recently, the axicon has become a promising optical device in elongating the length of a single filament [41,46,47]. Therefore, a vortex beam with a topological charge of 10, which propagates through a lens or an axicon, has been considered, respectively, during our work.
Figure 2a,b demonstrate the simulated vortex beam’s intensity distribution along the axis of x = 0 during nonlinear propagation when the laser powers are set as 20 times the critical power for self-focusing by using a lens and an axicon without air turbulence, respectively. Figure 2a shows that the vortex beam is focused by the lens of f = 50 cm, which corresponds to one half of the whole propagation distance, 100 cm. Curved filamentation is formed by the focusing lens, which can be explained by the strong energy confinement along a bent trajectory caused by vortex phase modulation introduced into the initial beam [48]. Figure 2c,d indicate the corresponding cross-section patterns of laser intensity at the propagation distance z = 60 cm. There are some ring structures shown in Figure 2c,d. Figure 2b,d shows that after being focused by an axicon, a vortex beam can be converted into a high-order Bessel beam. Moreover, the order of the Bessel beam is equal to the topological charge number of the vortex beam. The bottom angle of the axicon was set as 0.28° in the simulation, so that the non-diffraction zone of the axicon coincided with the propagation distance of 100 cm.
In order to perform the statistical regularities of vortex beam propagation in turbulent air, we used 20 independent simulations with different series of stochastic turbulence phase screens for each condition. The turbulent phase screens were introduced at the positions of z = 0 cm, 25 cm, 50 cm and 75 cm, corresponding to the beginning position of the initial vortex beam, the onset, the intermediate stage and the end of the filamentation process, respectively. The receiving screen was placed at the location of 100 cm to record the transverse intensity patterns during the propagation process of the laser beams.
Figure 3 demonstrates the simulated results of the statistical characteristics of vortex beam wandering focused by a lens when the air turbulence with the structure constant Cn2 = 5.6 × 10−13 cm−2/3 operates in the nonlinear propagation for every 25 cm of distance, respectively. Because of the ring structure of the transverse intensity distribution, as opposed to a single spot, the central coordinates (x0, y0) of the vortex beam were calculated with the following formula:
x 0 = x I ( x , y ) d x d y I ( x , y ) d x d y , y 0 = y I ( x , y ) d x d y I ( x , y ) d x d y
where I (x, y) refers to the ring structure intensity generated by the vortex beam. Figure 3 implies that the resulting vortex beam wandering was random due to the stochastic characteristic of turbulence. The red dots correspond to the central coordinates (x0, y0) of vortex beams. The simulations have been repeated 20 times for each position where the turbulent phase screens were introduced; thus, there are 20 red dots in each of Figure 3a–d. The average values of the transverse intensity central displacement <R0>, which is defined as R 0 = x 0 2 + y 0 2 , are approximately 128 μm, 101 μm, 71 μm and 40 μm, corresponding to the locations of 0 cm, 25 cm, 50 cm and 75 cm, where the turbulent phase screens were placed. The turbulence phase screens will lead to phase perturbation and provoke the instability of the wavefront continuously; thus, the influence of turbulent air on beam wandering is stronger when the turbulence occurs prior to the onset of propagation as compared to when it takes place at the middle or the end of propagation. Therefore, we can conclude that the air turbulence occurring earlier can result in a greater influence on the beam wandering. In the meantime, we numerically simulated a vortex beam propagating through a lens in turbulent air using 20 different series of stochastic turbulence phase screens during linear propagation. The average values <R0> calculated were 141 μm, 119 μm, 91 μm and 48 μm, when the phase screens were placed at 0 cm, 25 cm, 50 cm and 75 cm, respectively. The average values of spatial displacement become 110%, 118%, 128% and 120% of that during nonlinear propagation at the same location. It is well established that the turbulence effect’s contribution to the refractive index gradient is negligible compared with the nonlinear effect, which was induced by the Kerr effect within the filamentation. Previous work reported that filamentation was very resistant to turbulence once it was formed [34]. Thus, we included the consideration that the nonlinear evolution of the vortex beam can suppress the air turbulence better than linear evolution.
As a comparison, Figure 4 illustrates the simulated results of vortex beam wandering focused by an axicon when the air turbulence with the structure constant Cn2 = 5.6 × 10−13 cm−2/3 operates at different positions of 0 cm, 25 cm, 50 cm and 75 cm during nonlinear propagation. The numerically transverse intensity patterns were recorded at 100 cm. The average values of the transverse intensity central displacement <R0> are approximately 101 μm, 74 μm, 55μm and 33 μm, which become 79%, 73%, 77% and 78% of that focused by a lens at the same location. Using the same methods mentioned above, vortex beam wandering focused by an axicon was simulated in linear propagation. The average values <R0> calculated are 109 μm, 89 μm, 65 μm and 36 μm, which are all approximately 75% more or less than that focused by a lens during linear propagation. Because of the non-diffraction propagation, it is indicated that the axicon, as the focusing optical device, can suppress the optical perturbation induced by air turbulence.
The angle-of-arrival (AOA) fluctuation plays a significant role in FSO communication and atmospheric remote sensing [49]. The AOA fluctuation of an optical wavefront of the receiver aperture is associated with distortion in the focal plane of an imaging system. Under the geometrical optics method, the definition of AOA has the following form [50]:
β = arctan Δ S f ( r a d i a n ) ,
where ΔS represents the phase shift, and f’ is the focal length of the receiving antenna; we set f’ = 1 m for convenience in the simulations. The AOA fluctuation is defined as the standard deviation of β.
Figure 5 shows the simulated average values of AOA of vortex beams focused by a lens and an axicon when air turbulence was introduced at the various positions (z = 0, 25, 50 and 75 cm), respectively, during linear and nonlinear propagation. The error bar of Figure 5 denotes the AOA fluctuation. We can see that the average values of AOA decreased gradually when the air turbulence screen was placed from the beginning to the end of the propagation. Meanwhile, it can be observed that there are no obvious differences among the various positions where turbulent screens were installed. Therefore, the AOA fluctuation is not significantly affected by the different stages of propagation compared with the results of beam wandering as shown in Figure 3 and Figure 4. Meanwhile, the average values of AOA in nonlinear propagation are slightly smaller than that in the linear condition. As a comparison, we discovered, surprisingly, that the average values of AOA of the vortex beam focused by an axicon are nearly an order of magnitude lower than that focused by a lens. The results further confirm that a laser beam propagating through an axicon with a relatively large bottom angle (generally ≥ 0.2°) not only can suppress the turbulence, but can also maintain a more stable beam pattern.
In addition, we discussed the degree of transverse intensity patterns’ distortion. Figure 2b,d illustrate that the ring structure created by the vortex beam is usually of a circular pattern. However, due to the optical perturbation in the atmosphere, the circular-shaped ring structure is easily transformed into an elliptic or other irregular oval shape. Here, distortion factor δ is expressed by the following formula:
δ = Δ a Δ b ,
where Δ a is the maximum distance between the edges of the ring structure along the x-axis and Δ b is along the y-axis.
Figure 6 depicts the distortion factor δ of the ring structure generated by a lens and an axicon during nonlinear propagation when 20 different turbulence screens were located at z = 0 cm, 25 cm, 50 cm, 75 cm, respectively. Figure 6 shows that the variation in distortion factors focused by an axicon are more random and irregular than that focused by a lens during nonlinear propagation. It can be interpreted that the ring transverse intensity of vortex beams focused by a lens is usually larger than that focused by an axicon due to the divergence; therefore, small changes in Δa and Δb can lead to great fluctuations in distortion factor δ according to Equation (10).
In order to simulate the real effect of air turbulence in the atmosphere, the first turbulent phase screen was placed at the beginning of the propagation of the initial vortex beam and another three phase screens of the same series were introduced along the propagation path at every 25 cm simultaneously. The receiving screen was placed at the location of 100 cm to record the transverse intensity patterns. In other words, in each numerical simulation, a total of four phase screens were set up along the entire propagation path to fit the actual atmospheric conditions. The simulations have been repeated 20 times using a different series of air turbulence screens with the same structure constant of Cn2 = 5.6 × 10−13 cm−2/3. Figure 7 demonstrates the representative simulated vortex beam intensity distribution focused by a lens and an axicon with four turbulent screens simultaneously located at four different positions during the whole nonlinear propagation process when the laser powers are set as 20 times the critical power for self-focusing. Figure 7a,b depict the vortex beam intensity distribution along the axis of x = 0, while Figure 7c,d indicate the corresponding cross-section patterns of laser intensity at the propagation distance z = 70 cm. Compared with Figure 2, the vortex beam intensity distribution does not show cylindrical symmetry due to the effect of air turbulence. There are several “hot” spots along the ring structure shown in Figure 7c, which may lead to several single filaments. Meanwhile, it can be observed in Figure 7b that the high-order Bessel beam generated by the vortex beam propagating across the axicon does not transmit in a straight manner along the central axis. As shown in Figure 7d, the transverse intensity distribution of the Bessel beam is converted into an inhomogeneous elliptical pattern.
Although there are some other parameters in atmospheric air, such as temperature, relative humidity and wind speed, these effects eventually lead to the variation in turbulent structure constant Cn2 [51]. The optical effects of atmospheric turbulence depend primarily on the structure parameter Cn2. Consequently, we performed a statistical analysis of the transverse spatial displacements as a function of turbulent structure constant Cn2. As shown in Figure 8, by varying structure constant Cn2 from 10−13 to 10−11 cm−2/3, we simulated the nonlinear propagation process of a vortex beam focused by a lens and an axicon with four turbulent screens simultaneously. In the meantime, we set the topological charge as 4 in addition to 10 to investigate the properties of another mode of OAM state carried by vortex beam propagation in turbulent air. Each process of simulation was repeated 20 times with a different series of air turbulence screens. Figure 8 describes the average spatial displacements <R0> as a function of turbulent structure constant Cn2 when vortex beams carrying 4 and 10 modes of OAM state propagate through a lens and an axicon during nonlinear propagation (P = 20Pcr). The average values of spatial displacements <R0> become gradually larger with the increasing structure constant Cn2. According to Equations (1)–(7), it can be easily understood that the turbulent structure constant contributes to the refractive index gradient more efficiently than any other parameter. Thus, the larger the structure constant is, the more intensively the air turbulence affects the beam wandering. Comparing the topological charge of 10 with that of 4 in Figure 8, the relatively larger topological charge behaves better with respect to suppressing the effect of air turbulence, which is consistent with the previously reported conclusion [45].

4. Conclusions

In summary, we have investigated the statistical regularities of a vortex beam propagating through a lens or an axicon with different series of turbulent air phase screens. We mainly simulated the evolution process of beam wandering owing to turbulent phase screens with various parameters. The results revealed that the distortion of the spatial displacement was sensitive to the location and the structure constant Cn2 of the turbulent screens. The effect of the optical perturbation owing to air turbulence is more intensive when the phase screen is located at the early stages of the propagation. Turbulent structure constant Cn2, which can affect the refractive index gradient efficiently, is a significant parameter in the atmospheric air. Compared with linear propagation, we find that the nonlinear process of vortex beams can effectively suppress the effect of air turbulence. In addition, the axicon can be regarded as a very useful optical device, since it can not only suppress the turbulence but also maintain a stable beam pattern. Furthermore, we simulated vortex beams carrying different OAM modes propagating in turbulent air. It can be concluded that a large topological charge can reduce the influence of air turbulence during propagation. Our outcomes are valuable for many applications in atmospheric air, especially for optical communication and remote sensing.

Author Contributions

Conceptualization, methodology, writing—review and editing, D.Z.; data curation, writing—original draft preparation, C.L.; supervision, project administration, X.S.; investigation, Y.L.; visualization, Y.Z.; resources, software, H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 61605144).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.L.; Chen, Y.P.; Cui, S.; Feng, H.; Chen, M.; Zhang, Z.; Yu, J.; Chen, L.M.; Sheng, Z.M.; Zhang, J. Manipulation of polarizations for broadband terahertz waves emitted from laser plasma filaments. Nat. Photonics 2018, 12, 554–559. [Google Scholar] [CrossRef]
  2. Dey, I.; Jana, K.; Fedorov, V.Y.; Koulouklidis, A.D.; Mondal, A.; Shaikh, M.; Sarkar, D.; Lad, A.D.; Tzortzakis, S.; Couairon, A.; et al. Highly efficient broadband terahertz generation from ultrashort laser filamentation in liquids. Nat. Commun. 2017, 8, 1–7. [Google Scholar] [CrossRef]
  3. Varela, O.; Alonso, B.; Sola, I.J.; Román, S.; Zaïr, A.; Méndez, C.; Roso, L. Self-compression controlled by the chirp of the input pulse. Opt. Lett. 2010, 35, 3649–3651. [Google Scholar] [CrossRef]
  4. Adachi, S.; Suzuki, T. Self-compression of femtosecond deep-ultraviolet pulses by filamentation in krypton. Opt. Lett. 2017, 42, 1883–1886. [Google Scholar] [CrossRef]
  5. Xu, H.L.; Lötstedt, E.; Lwasaki, A.; Yamanouchi, K. Sub-10-fs population inversion in N2+ in air lasing through multiple state coupling. Nat. Commun. 2015, 6, 8347. [Google Scholar] [CrossRef]
  6. Yao, J.P.; Jiang, S.C.; Chu, W.; Zeng, B.; Wu, C.Y.; Lu, R.F.; Li, Z.T.; Xie, H.Q.; Li, G.H.; Wang, Z.S.; et al. Population Redistribution Among Multiple Electronic States of Molecular Nitrogen Ions in Strong Laser Fields. Phys. Rev. Lett. 2016, 116, 143007. [Google Scholar] [CrossRef]
  7. Xu, H.L.; Cheng, Y.; Chin, S.-L.; Sun, H.-B. Femtosecond laser ionization and fragmentation of molecules for environmental sensing. Laser Photonics Rev. 2015, 9, 275–293. [Google Scholar] [CrossRef]
  8. Odhner, J.; Levis, R. Optical Spectroscopy Using Gas-Phase Femtosecond Laser Filamentation. Rev. Phys. Chem. 2014, 65, 605–628. [Google Scholar] [CrossRef]
  9. Couairon, A.; Mysyrowicz, A. Femtosecond filamentation in transparent media. Phys. Rep. 2007, 441, 47–189. [Google Scholar] [CrossRef]
  10. Bergé, L.; Skupin, S.; Kasparian, J.; Wolf, J.-P. Ultrashort filaments of light in weakly ionized, optically transparent media. Rep. Prog. Phys. 2007, 70, 1633–1713. [Google Scholar] [CrossRef] [Green Version]
  11. Kasparian, J.; Wolf, J.-P. Physics and applications of atmospheric nonlinear optics and filamentation. Opt. Express 2008, 16, 466–493. [Google Scholar] [CrossRef] [PubMed]
  12. Béjot, P.; Kasparian, J.; Henin, S.; Loriot, V.; Vieillard, T.; Hertz, E.; Faucher, O.; Lavorel, B.; Wolf, J.-P. Higher-Order Kerr Terms Allow Ionization-Free Filamentation in Gases. Phys. Rev. Lett. 2010, 104, 103903. [Google Scholar] [CrossRef] [PubMed]
  13. Béjot, P.; Kasparian, J.; Lavorel, B.; Wolf, J.-P.; Faucher, O. Transition from Plasma-Driven to Kerr-Driven Laser Filamentation. Phys. Rev. Lett. 2011, 106, 243902. [Google Scholar] [CrossRef] [PubMed]
  14. Schimmel, G.; Produit, T.; Mongin, D.; Kasparian, J.; Wolf, J.-P. Free space laser telecommunication through fog. Optica 2018, 5, 1338–1341. [Google Scholar] [CrossRef]
  15. Schroeder, M.C.; Andral, U.; Wolf, J.-P. Opto-mechanical expulsion of individual micro-particles by laser-induced shockwave in air. AIP Adv. 2022, 12, 095119. [Google Scholar] [CrossRef]
  16. Willner, A.E.; Huang, H.; Yan, Y.; Ren, Y.; Ahmed, N.; Xie, G.; Bao, C.; Li, L.; Cao, Y.; Zhao, Z.; et al. Optical communications using orbital angular momentum beams. Adv. Opt. Photonics 2015, 7, 66–106. [Google Scholar] [CrossRef]
  17. Kim, H.-C.; Lee, Y.H. Hermite–Gaussian and Laguerre–Gaussian beams beyond the paraxial approximation. Opt. Commun. 1999, 169, 9–16. [Google Scholar] [CrossRef]
  18. Wang, A.D.; Zhu, L.; Deng, M.L.; Lu, B.; Guo, X.J. Experimental demonstration of OAM-basedtransmitter mode diversity data transmission under atmosphere turbulence. Opt. Express 2021, 29, 13171–13182. [Google Scholar] [CrossRef]
  19. Zhu, L.; Wang, A.D.; Deng, M.L.; Lu, B.; Guo, X.J. Free-space optical communication with quasi-ring Airy vortex beam under limited-size receiving aperture and atmospheric turbulence. Opt. Express 2021, 29, 32580–32590. [Google Scholar] [CrossRef]
  20. Roux, F.S.; Wellens, T.; Shatokhin, V.N. Entanglement evolution of twisted photons in strong atmospheric turbulence. Phys. Rev. A 2015, 92, 012326. [Google Scholar] [CrossRef]
  21. Zhang, Y.W.; Agnew, M.; Roger, T.; Roux, F.S.; Konrad, T.; Faccio, D.; Leach, J.; Forbes, A. Simultaneous entanglement swapping of multiple orbital angular momentum states of light. Nat. Commun. 2017, 8, 1–7. [Google Scholar] [CrossRef]
  22. Leonhard, N.D.; Shatokhin, V.N.; Buchleitner, A. Universal entanglement decay of photonic-orbital-angular-momentum qubit states in atmospheric turbulence. Phys. Rev. A 2015, 91, 012345. [Google Scholar] [CrossRef]
  23. Wang, J.; Yang, J.Y.; Fazal, I.M.; Ahmed, N.; Yan, Y.; Huang, H.; Ren, Y.X.; Yue, Y.; Dolinar, S.; Tur, M.; et al. Terabit free-space data transmission employing orbital angular momentum multiplexing. Nat. Photonics 2012, 6, 488–496. [Google Scholar] [CrossRef]
  24. Li, J.; Chen, X.; McDuffie, S.; Najjar, M.A.M.; Rafsanjani, S.M.H.; Korotkova, O. Mitigation of atmospheric turbulence with random light carrying OAM. Opt. Commun. 2019, 446, 178–185. [Google Scholar] [CrossRef]
  25. Ayeni, O.O.; Ojo, J.S.; Owolawi, P.A.; Ajewole, M.O. Performance comparison of diversity techniques with Addictive White Gaussian Channel (AWGC) in free space optical communication (FSO) under atmospheric turbulence scenario. Int. J. Phys. Sci. 2023, 18, 12–18. [Google Scholar] [CrossRef]
  26. Walsh, S.M.; Karpathakis, S.F.E.; McCann, A.S.; Dix-matthews, B.P.; Frost, A.M.; Gozzard, D.R.; Gravestock, C.T.; Schediwy, S.W. Demonstration of 100 Gbps coherent free-space optical communications at LEO tracking rates. Sci. Rep. 2022, 12, 18345. [Google Scholar] [CrossRef]
  27. Ebrahimi, F.; Ghassemlooy, Z.; Olyaee, S. Investigation of a hybrid OFDM-PWM/PPM visible light communications system. Opt. Commun. 2018, 429, 65–71. [Google Scholar] [CrossRef]
  28. Eso, E.; Ghassemlooy, Z.; Zvanovec, S.; Sathian, J.; Abadi, M.M.; Younus, O.I. Performance of Vehicular Visible Light Communications underthe Effects of Atmospheric Turbulence with Aperture Averaging. Sensors 2021, 21, 2751. [Google Scholar] [CrossRef]
  29. Kasparian, J.; Rodriguez, M.; Méjean, G.; Yu, J.; Salmon, E.; Wille, H.; Bourayou, R.; Frey, S.; André, Y.-B.; Mysyrowicz, A.; et al. White-Light Filaments for Atmospheric Analysis. Science 2003, 301, 61–64. [Google Scholar] [CrossRef]
  30. Kandidov, V.P.; Kosareva, O.G.; Tamarov, M.P.; Brodeur, A.; Chin, S.L. Nucleation and random movement of filaments in the propagation of high-power laser radiation in a turbulent atmosphere. Quantum Electron. 1999, 29, 911–915. [Google Scholar] [CrossRef]
  31. Peñano, J.R.; Sprangle, P.; Hafizi, B.; Ting, A.; Gordon, D.F.; Kapetanakos, C.A. Propagation of ultra-short, intense laser pulses in air. Phys. Plasmas 2004, 11, 2865–2874. [Google Scholar] [CrossRef]
  32. Sunilkumar, K.; Anand, N.; Moorthy, K.K.; Ilavazhagan, G. Enhanced optical pulse broadening in free-space optical links due to the radiative effects of atmospheric aerosols. Opt. Express 2021, 29, 865–876. [Google Scholar] [CrossRef]
  33. Chin, S.L.; Talebpour, A.; Yang, J.; Petit, S.; Kandidov, V.P.; Kosareva, O.G.; Tamarov, M.P. Filamentation of femtosecond laser pulses in turbulent air. Appl. Phys. B-Lasers Opt. 2002, 74, 67–76. [Google Scholar] [CrossRef]
  34. Houard, A.; Franco, M.; Prade, B.; Durécu, A.; Lombard, L.; Bourdon, P.; Vasseur, O.; Fleury, B.; Robert, C.; Michau, V.; et al. Femtosecond filamentation in turbulent air. Phys. Rev. A 2008, 78, 033804. [Google Scholar] [CrossRef]
  35. Gill, H.S.; Singh, M.L.; Singh, M.; Priyanka; Kaur, S.; Kaur, H. Analysis of full reference quality metrics for image transmission over a MIMO OWC channel under varying turbulent conditions. Int. J. Commun. Syst. 2023, 36, e5426. [Google Scholar]
  36. Naghshvarianjahromi, M.; Kumar, S.; Deen, M.J. Free Space Ground to Satellite Optical Communications Using Kramers-Kronig Transceiver in the Presence of Atmospheric Turbulence. Sensors 2022, 22, 3435. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, W.W.; Chin, S.L. Abnormal wavelength dependence of the self-cleaning phenomenon during femtosecond-laser-pulse filamentation. Phys. Rev. A 2007, 76, 013826. [Google Scholar] [CrossRef]
  38. Sun, X.D.; Gao, H.; Zhang, S.W.; Liu, W.W. Numerical simulation of the generation of multiple laser filaments by an axicon array. J. Mod. Opt. 2013, 60, 1637–1643. [Google Scholar] [CrossRef]
  39. Talebpour, A.; Yang, J.; Chin, S.L. Semi-empirical model for the rate of tunnel ionization of N2 and O2 molecule in an intense Ti:sapphire laser pulse. Opt. Commun. 1999, 163, 29–32. [Google Scholar] [CrossRef]
  40. Dubietis, A.; Tamošauskas, G.; Fibich, G.; Ilan, B. Multiple filamentation induced by input-beam ellipticity. Opt. Lett. 2004, 29, 1126–1128. [Google Scholar] [CrossRef]
  41. Sun, X.D.; Gao, H.; Zeng, B.; Xu, S.Q.; Liu, W.W.; Cheng, Y.; Xu, Z.Z.; Mu, G.G. Multiple filamentation generated by focusing femtosecond laser with axicon. Opt. Lett. 2012, 37, 857–859. [Google Scholar] [CrossRef] [PubMed]
  42. Zeng, T.; Gao, H.; Sun, X.D.; Liu, W.W. Turbulence-induced beam wandering during femtosecond laser filamentation. Chin. Opt. lett. 2015, 13, 070008. [Google Scholar] [CrossRef]
  43. Johansson, E.M.; Gavel, D.T. Simulation of stellar speckle imaging. Proc. SPIE 1994, 2200, 372–383. [Google Scholar]
  44. Charnotskii, M. Comparison of four techniques for turbulent phase screens simulation. J. Opt. Soc. Am. A 2020, 37, 738–747. [Google Scholar] [CrossRef]
  45. Yuan, Y.S.; Lei, T.; Li, Z.H.; Li, Y.J.; Gao, S.C.; Xie, Z.W.; Yuan, X.C. Beam wander relieved orbital angular momentum communication in turbulent atmosphere using Bessel beams. Sci. Rep. 2017, 7, 2045–2322. [Google Scholar] [CrossRef]
  46. Polynkin, P.; Kolesik, M.; Roberts, A.; Faccio, D.; Trapani, P.D.; Moloney, J. Generation of extended plasma channels in air using femtosecond Bessel beams. Opt. Express 2008, 16, 15733–15740. [Google Scholar] [CrossRef]
  47. Akturk, S.; Zhou, B.; Franco, M.; Couairon, A.; Mysyrowicz, A. Generation of long plasma channels in air by focusing ultrashort laser pulses with an axicon. Opt. Commun. 2009, 282, 129–134. [Google Scholar] [CrossRef]
  48. Papazoglou, D.G.; Suntsov, S.; Abdollahpour, D.; Tzortzakis, S. Tunable intense Airy beams and tailored femtosecond laser filaments. Phys. Rev. A 2010, 81, 061807. [Google Scholar] [CrossRef]
  49. Liao, T.H.; Gao, Q. Image recovery from double amplitudes in fractional Fourier domain. Chin. Phys. 2006, 15, 347–352. [Google Scholar]
  50. Liu, Y.; Song, X.; Ni, X.; Qi, J.; Liu, Z. Measurement of angle-of-arrival fluctuations over a real atmospheric turbulent path. In Proceedings of the 2015 International Conference on Optoelectronics and Microelectronics, Changchun China, 16–18 July 2015. [Google Scholar]
  51. Frederickson, P.; Davidson, K.; Zeisse, C.; Bendall, I. A Comparison of Near-surface Bulk and Scintillation Cn2 Measure-ments during EOPACE. Proc. SPIE 1998, 3433, 77–88. [Google Scholar]
Figure 1. Representative turbulent phase screen with the structure constant Cn2 = 5.6 × 10−13 cm−2/3 and the distance Δz = 1 m.
Figure 1. Representative turbulent phase screen with the structure constant Cn2 = 5.6 × 10−13 cm−2/3 and the distance Δz = 1 m.
Sensors 23 01772 g001
Figure 2. Simulated intensity distribution without air turbulence by using (a) a lens and (b) an axicon longitudinal distribution during nonlinear propagation when the ultrafast laser powers are set to 20 times the critical power, respectively. (c,d) transverse distribution at z = 60 cm during nonlinear propagation corresponding to (a,b).
Figure 2. Simulated intensity distribution without air turbulence by using (a) a lens and (b) an axicon longitudinal distribution during nonlinear propagation when the ultrafast laser powers are set to 20 times the critical power, respectively. (c,d) transverse distribution at z = 60 cm during nonlinear propagation corresponding to (a,b).
Sensors 23 01772 g002
Figure 3. Spatial displacements focused by a lens (f = 50 cm) when the air turbulence (Cn2 = 5.6 × 10−13 cm−2/3) operates at z = (a) 0 cm, (b) 25 cm, (c) 50 cm and (d) 75 cm during nonlinear propagation.
Figure 3. Spatial displacements focused by a lens (f = 50 cm) when the air turbulence (Cn2 = 5.6 × 10−13 cm−2/3) operates at z = (a) 0 cm, (b) 25 cm, (c) 50 cm and (d) 75 cm during nonlinear propagation.
Sensors 23 01772 g003
Figure 4. Spatial displacements focused by an axicon (the bottom angle is 0.28°) when the air turbulence (Cn2 = 5.6 × 10−13 cm−2/3) operates at z = (a) 0 cm, (b) 25 cm, (c) 50 cm and (d) 75 cm during nonlinear propagation.
Figure 4. Spatial displacements focused by an axicon (the bottom angle is 0.28°) when the air turbulence (Cn2 = 5.6 × 10−13 cm−2/3) operates at z = (a) 0 cm, (b) 25 cm, (c) 50 cm and (d) 75 cm during nonlinear propagation.
Sensors 23 01772 g004
Figure 5. The simulated average values of AOA of vortex beams focused by (a) a lens and (b) an axicon when turbulence was introduced at various displacements, respectively, during linear (blue) and nonlinear (red) propagation.
Figure 5. The simulated average values of AOA of vortex beams focused by (a) a lens and (b) an axicon when turbulence was introduced at various displacements, respectively, during linear (blue) and nonlinear (red) propagation.
Sensors 23 01772 g005
Figure 6. The distortion factor δ of the ring structure generated by (a) a lens and (b) an axicon during nonlinear propagation with 20 different turbulence screens located at various positions, respectively.
Figure 6. The distortion factor δ of the ring structure generated by (a) a lens and (b) an axicon during nonlinear propagation with 20 different turbulence screens located at various positions, respectively.
Sensors 23 01772 g006
Figure 7. Representative simulated intensity distribution with 4 turbulent phase screens in the whole propagation process by using (a) a lens and (b) an axicon, showing longitudinal distribution during nonlinear propagation, respectively; (c,d) transverse distribution at z = 70 cm during nonlinear propagation corresponding to (a,b).
Figure 7. Representative simulated intensity distribution with 4 turbulent phase screens in the whole propagation process by using (a) a lens and (b) an axicon, showing longitudinal distribution during nonlinear propagation, respectively; (c,d) transverse distribution at z = 70 cm during nonlinear propagation corresponding to (a,b).
Sensors 23 01772 g007
Figure 8. The average spatial displacements <R0> as a function of turbulent structure constant Cn2 when vortex beams carrying 4 (blue) and 10 (red) modes of OAM state propagate through (a) a lens and (b) an axicon during nonlinear propagation.
Figure 8. The average spatial displacements <R0> as a function of turbulent structure constant Cn2 when vortex beams carrying 4 (blue) and 10 (red) modes of OAM state propagate through (a) a lens and (b) an axicon during nonlinear propagation.
Sensors 23 01772 g008
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

Zhu, D.; Li, C.; Sun, X.; Liu, Y.; Zhang, Y.; Gao, H. The Effect of Air Turbulence on Vortex Beams in Nonlinear Propagation. Sensors 2023, 23, 1772. https://doi.org/10.3390/s23041772

AMA Style

Zhu D, Li C, Sun X, Liu Y, Zhang Y, Gao H. The Effect of Air Turbulence on Vortex Beams in Nonlinear Propagation. Sensors. 2023; 23(4):1772. https://doi.org/10.3390/s23041772

Chicago/Turabian Style

Zhu, Di, Chunhua Li, Xiaodong Sun, Yali Liu, Yuqi Zhang, and Hui Gao. 2023. "The Effect of Air Turbulence on Vortex Beams in Nonlinear Propagation" Sensors 23, no. 4: 1772. https://doi.org/10.3390/s23041772

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