Abstract
Applying the angular spectrum theory, we derive the expression of a new Hermite-Gaussian (HG) vortex beam. Based on the new Hermite-Gaussian (HG) vortex beam, we establish the model of the received probability density of orbital angular momentum (OAM) modes of this beam propagating through a turbulent ocean of anisotropy. By numerical simulation, we investigate the influence of oceanic turbulence and beam parameters on the received probability density of signal OAM modes and crosstalk OAM modes of the HG vortex beam. The results show that the influence of oceanic turbulence of anisotropy on the received probability of signal OAM modes is smaller than isotropic oceanic turbulence under the same condition, and the effect of salinity fluctuation on the received probability of the signal OAM modes is larger than the effect of temperature fluctuation. In the strong dissipation of kinetic energy per unit mass of fluid and the weak dissipation rate of temperature variance, we can decrease the effects of turbulence on the received probability of signal OAM modes by selecting a long wavelength and a larger transverse size of the HG vortex beam in the source’s plane. In long distance propagation, the HG vortex beam is superior to the Laguerre-Gaussian beam for resisting the destruction of oceanic turbulence.
© 2017 Optical Society of America
1. Introduction
In recent years, optical beams with orbital angular momentum (OAM) propagating in random media have attracted much attention due to the applications in optical wireless communication. The main motivation is that beams carrying OAM can realize arbitrary base-N quantum digits in principle [1,2], which offers the additional degrees of freedom for information coding and increases information capacity of the communication link. Due to the presence of various constituents in the ocean and the fluctuations of the refractive index of water, the propagation of optical signal through an oceanic medium will suffer attenuation and wavefront distortion [3,4]. It is well known that oceanic turbulence will distort OAM modes of optical beams, which introduces intermodal crosstalk between channels. After Nikishov et al. [5] presented the power spectrum of oceanic turbulence, there are great attentions paid to understand the laws of OAM modes propagation in the oceanic turbulence [6–10]. Huang et al. [9] analyzed the influence of beam parameters, the rate of dissipation of turbulent and temperature salinity fluctuations on the average intensity of Gaussian Schell-model vortex beam propagating through oceanic turbulence, and showed that the partially coherent vortex beams have robust turbulent resistance. Liu et al. [10] studied the evolution properties of partially coherent flat-topped vortex hollow beams propagating through oceanic turbulence, and revealed that this beam evolves into the Gaussian-like beam for the long distance propagation and this beam with higher order M will lose its initial dark hollow center slower. The effects of oceanic turbulence and light source parameters on the received probability of partially coherent Laguerre–Gaussian (LG) vortex beam was investigated [6], and showed that the partially coherent vortex beam provides worse performance than fully coherent vortex beam. The channel capacity of Bessel–Gauss Beams propagating through turbulent ocean was investigated [7], which revealed that BG beams are better than LG beams to mitigate the effects of oceanic turbulence because of the non-diffraction and self-healing characteristics of BG beams.
The experimental and theoretical results demonstrated that turbulence can be anisotropic [11–17]. Due to the earth's rotation, the turbulence of ocean will also be anisotropic [18,19]. The average polarizability of Gaussian Schell-model quantized beams and the spreading and wandering of Gaussian–Schell model laser beams in anisotropic oceanic turbulence were studied respectively [18,19]. However, the energy distribution of OAM modes of LG beams decreases with the increasing quantum number of OAM modes, and the crosstalk between OAM modes increases with increasing the quantum number of OAM modes [2,7]. We cannot make use of the high order OAM modes of LG beams to improve the capacity of the free-space optical links. To improve capacity of the free-space optical links with the high order OAM modes, a more suitable vortex beam should be found. To the best of our knowledge, the transmission properties of Hermite-Gaussian (HG) vortex beams carrying OAM in anisotropic turbulent ocean have not been reported yet.
In this paper, we first develop a new angular spectrum of the vortex beam and derive a new HG vortex beam based on the angular spectrum theory. Then, we study the propagation of the new HG vortex beam in anisotropic turbulent ocean, and establish the model of the received probability distribution of OAM modes. In the end, we discuss the influences of the HG vortex beam parameters and the anisotropic turbulent ocean on probability density of signal OAM modes and crosstalk OAM modes in weak horizontal oceanic turbulent channels.
2. Anisotropic oceanic turbulence
In this paper, we study the HG vortex beam propagation in the turbulence of the pure seawater [20]. Due to the earth’s rotation, the turbulence of the ocean will be anisotropic [11]. For the sake of simplicity, here we employ the method that the anisotropic power spectrum has been proposed in atmospheric turbulence [15–17] and consider the turbulent anisotropy only existing along the propagation direction of the beam as the discussion of Ref [18]. We can present the anisotropic power spectrum of oceanic turbulence by an effective anisotropic factor . The anisotropic spectrum of oceanic turbulence can be expressed as
where is the spatial frequency of turbulent fluctuations, . is the temperature structure constant for given [20] that characters the temperature fluctuation strength of oceanic turbulence, is the rate of dissipation of kinetic energy per unit mass of fluid ranging from to , is the dissipation rate of temperature variance and has the range from to , is the inner scale of oceanic turbulence, , , , , is the ratio of temperature and salinity contributions to the refractive index spectrum, which in the seawaters can vary in the interval, with and the corresponding to dominating temperature-induced and salinity-induced optical turbulence, respectively; .3. Propagation of HG vortex beams in the absence of turbulence
The propagation of a monochromatic scalar light field in free space is represented by the angular spectrum for the half-space [21,22]:
where is the angular frequency, is the wavenumber in free space, is the wavelength, and represent the components of the wave vector along the axis and axis respectively, is transmission distance, represents the angular spectrum.For the lowest-order Gaussian beam field, the product of the divergence and transverse spread reaches a minimum and can be represented by [22]:
The expressions of the divergence and transverse spread for the symmetry can be defined as
where represents the partial derivative.The product of the divergence and transverse spread in or direction meets the minimum respectively. In direction, the product can be expressed by [22]:
Applying the Cauchy-Schwartz inequality
where and stand for the arbitrary functions, the superscript represents the complex conjugate.Equation (7) can reach its minimum in the condition that and are proportional to each other. The angular spectrum for the Gaussian beams can be represented by
It is known that the angular spectrum of the HG beam is given by [21,22]
Inserting the vortex factor into the angular spectrum of the HG beam, we can obtain a new angular spectrum. We call this new angular spectrum as the angular spectrum of the HG vortex beam and express it as
where and are integers.Substituting Eq. (10) into Eq. (2), the propagation of the HG vortex beam in free space of absence turbulence can be investigated by modulating the angular spectrum
whereand , is the initial radius of beams in the source plane.According to these relations
Equation (11) can be rewritten byBy setting the new variables and , the partial derivative becomes . Therefore, Eq. (16) can be rewritten by
With the help of the integral expression [21–23]
the HG vortex beam is given byNeglecting the constant factor, we can obtain the HG vortex beams by using the following mathematical formulae and as
where ,, and . The Fig. 1 presents the intensity distribution for HG vortex beams in the planes .
Fig. 1 Intensity profiles of HG vortex beams in the absence of turbulence (a) , (b) ,(c) . When in the plane, the intensity of the HG vortex beam decreases with deviating from the optical axis center. In addition, the intensity forms annular shapes along the optical axis center with increasing quantum number .
In particular, when , the new HG vortex beam can degrade into HG beam as
where represents Hermite polynomials.Moreover, let and , the new HG vortex beam can be reduced as beam [24] and can be given by
where represents Laguerre polynomials.4. Probability density of OAM modes in oceanic turbulence
The HG vortex beam is a superposition of eigenstates and the conditional probability of obtaining a measurement of the OAM of a photon is obtained by summing the probabilities associated with that eigenvalue [2],
where is the quantum number of the OAM modes of beams in turbulence, and is the superposition coefficients.The aberrations are random, so by taking the ensemble average, the probability that a measurement of the OAM will yield a value is found as follows,
where denotes average over the ensemble of oceanic turbulence.By the superposition theory of plane waves with phase , the can be written as
where represents the field amplitude of the HG vortex beam propagation in oceanic turbulence.The complex amplitude of the HG vortex beam is given by [25]
where is the phase distortion caused by oceanic turbulence.The ensemble average of the superposition coefficients is represented as
where denotes average over the ensemble of the source and turbulence, denotes average over the ensemble of the light source and denotes average over the ensemble of oceanic turbulence. represents the coherence function of the phase, can be expressed aswhere , ,, and , and is given bySubstituting Eq. (28) into Eq. (24), Eq. (24) can be rewritten as
Since the aberrations introduced by oceanic turbulence are random variables, the coherence function of the phase can be expressed as
where is the phase structure function of turbulent aberrations. For Kolmogorov turbulence phase aberrations, the phase structure function of turbulent aberrations is given by [25,26]where represents the spatial coherence radius of a spherical wave propagating in the oceanic turbulence.Make use of the relations and , the coherence function of the phase is given by
In Markov approximation, i.e. , the spatial coherence radius can be written as
After a long but straightforward calculation [4], we can obtain the lateral coherence length of the spherical wave in anisotropic oceanic turbulence. The Eq. (34) is reduced to
Based on the Eq. (30) and Eq. (33), we can obtain the received probability [16]
whereFor , is the received probability of transmitting OAM modes . In addition, when , is the received probability of transmitting OAM modes with index followed by detection of OAM mode l.
5. Numerical analysis
We now investigate the statistical properties of OAM modes of the HG vortex beam propagating through oceanic turbulence by numerical simulation. In the following analysis, the calculation parameters of HG vortex beams are , , , , , , , , , , , unless other variable parameters are specified in calculation.
Figure 2 displays the received probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence for several values of parameters and . It can be seen that the received probability of signal OAM modes decreases monotonously with increasing propagation distance. Figure 2(a) indicates that the received probability of signal OAM modes of the beam is smaller than the HG vortex beam, when the transmission distance is greater. It means that the HG vortex beam has a stronger resistance to the oceanic turbulence than the beam, and the resistance increases with the increase of transmission distance. In addition, from Fig. 2(a-b), the received probability of signal OAM modes decreases with increasing OAM mode number for a given propagation distance. Therefore, the HG vortex beam with smaller is less affected by oceanic turbulence than the beam from the aspect of the signal OAM probability in long distance transmission.

Fig. 2 Probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence for several values of parameter (a) , (b) respectively.
The discussion of Fig. 2 is related to the received probability of signal OAM modes. The intermodel crosstalk between channels can be caused by optical OAM beams propagation through oceanic turbulence, which represents the OAM energy level difference . Figure 3 displays the difference between the HG vortex beams and the conventional LG beam propagation through oceanic turbulence causing the received probability of crosstalk OAM modes. We find that the received probability of crosstalk OAM modes of the conventional LG beam is obviously larger than that of the HG vortex beam. That is to say, the effects of oceanic turbulence on the HG vortex beam compared to the conventional LG beam are significantly less. The results will be very useful for underwater optical communications.

Fig. 3 Probability of crosstalk OAM modes of the HG vortex beam propagation through oceanic turbulence versus propagation distance for different modes .
For the HG vortex beam propagation in oceanic turbulence, the received probability of signal OAM modes will be improved in anisotropic turbulence in Fig. 4. Figure 4 also shows that the received probability of signal OAM modes of the HG vortex beam as functions of propagation distance for different anisotropic factor . For a given value of propagation distance , when the anisotropic factor increases, the signal probability of HG vortex beams increases.

Fig. 4 Probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence versus propagation distance for different anisotropic factor .
To comprehend the regularity of the received probability of signal OAM modes of the HG vortex beam, we calculate the influence of wavelength for different inner scale in Fig. 5. We deduce from Fig. 5 that, for any given inner scale , there is an increase of the received probability of signal OAM modes with increasing . The result shows that the HG vortex beam with longer wavelength is less affected by oceanic turbulence. In addition, Fig. 5 also shows that the received probability of signal OAM modes increases with the increase of for any given . The phenomenon can be explained that the larger inner scale of turbulent eddy results in small light scattering, and then the small light scattering produces a small energy loss. Thus, the received probability of signal OAM modes of the HG vortex beam can be improved under the greater inner scale of oceanic turbulence.

Fig. 5 Probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence versus wavelength for different inner scale .
We also analyze the influence of the dissipation rate of temperature variance and the source’s transverse size on the received probability of signal OAM modes for the HG vortex beam. From Fig. 6, we can see that, as the source’s transverse size increases, the received probability of signal OAM modes of the HG vortex beam increases. This is caused by the larger source’s transverse size to make a smaller beam spread and beam wander which are caused by the diffraction and refringence of the vacuum and oceanic turbulence, respectively. The smaller beam spread and beam wander lead to the smaller transformations of the space distribution of beams propagation in the link. In addition, as would be expected, the received probability of signal OAM modes of the HG vortex beam decreases with increasing dissipation rate of temperature variance from Fig. 6, especially in the larger source’s transverse size .

Fig. 6 Probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence versus the dissipation rate of temperature variance for different the source’s transverse size .
To further comprehend the influence of oceanic turbulence parameters, in particular the rate of dissipation of kinetic energy per unit mass of fluid and the ratio of temperature and salinity contributions to the refractive index spectrum on the received probability of the signal OAM mode of the HG vortex beam. Figure 7 shows that the received probability of the signal OAM mode is a function of for different . It is evident from Fig. 7 that the larger value of will obtain the larger received probability of signal OAM modes. This implies that the salinity fluctuations dominated turbulence affects the HG vortex beam propagation more than the temperature fluctuations dominated turbulence. Figure 7 also shows that the larger rate of dissipation of kinetic energy per unit mass of fluid enhances the detection probability of the signal OAM mode of the HG vortex beam. This conclusion comes from the larger rate of dissipation of kinetic energy per unit mass of fluid reducing oceanic turbulence.

Fig. 7 Probability of signal OAM modes of the HG vortex beam propagation through oceanic turbulence versus the rate of dissipation of kinetic energy per unit mass of fluid for different the ratio of temperature and salinity contributions to the refractive index spectrum .
6. Conclusion
We develop a new angular spectrum of the vortex beam based on the angular spectrum of the HG beam and the vortex factor , which evolves a new HG vortex beam. Then we establish the model of the received probabilities of OAM modes of the HG vortex beam propagation through anisotropic oceanic turbulence. The influence of the HG vortex beam on the received probabilities of signal and crosstalk OAM modes shows that the HG vortex beam is superior to the LG beam propagation in anisotropic oceanic turbulence for long-distance transmission. We also find that the HG vortex beam has robust turbulent resistance under the condition of anisotropic oceanic turbulence, longer , , , , , or smaller . The HG vortex beam can realize the goal for enhancing the link capacity of underwater optical communication.
Funding
Fundamental Research Funds for the Central Universities (Grant No.JUSRP51716A).
References and links
1. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef] [PubMed]
2. C. Paterson, “Atmospheric turbulence and orbital angular momentum of single photons for optical communication,” Phys. Rev. Lett. 94(15), 153901 (2005). [CrossRef] [PubMed]
3. Y. Baykal, “Intensity fluctuations of multimode laser beams in underwater medium,” J. Opt. Soc. Am. A 32(4), 593–598 (2015). [CrossRef] [PubMed]
4. Y. Wu, Y. Zhang, Y. Li, and Z. Hu, “Beam wander of Gaussian-Schell model beams propagating through oceanic turbulence,” Opt. Commun. 371, 59–66 (2016). [CrossRef]
5. V. V. Nikishov and V. I. Nikishov, “Spectrum of turbulent fluctuations of the sea-water refraction index,” Int. J. Fluid Mech. Res. 27(1), 82–98 (2000). [CrossRef]
6. M. Cheng, L. Guo, J. Li, Q. Huang, Q. Cheng, and D. Zhang, “Propagation of an optical vortex carried by a partially coherent Laguerre-Gaussian beam in turbulent ocean,” Appl. Opt. 55(17), 4642–4648 (2016). [CrossRef] [PubMed]
7. M. Cheng, L. Guo, J. Li, and Y. Zhang, “Channel capacity of the OAM based free-space optical communication links with Bessel–Gauss beams in turbulent ocean,” IEEE Photonics J. 8, 7901411 (2016).
8. J. Xu and D. Zhao, “Propagation of a stochastic electromagnetic vortex beam in the oceanic turbulence,” Opt. Laser Technol. 57, 189–193 (2014). [CrossRef]
9. Y. Huang, B. Zhang, Z. Gao, G. Zhao, and Z. Duan, “Evolution behavior of Gaussian Schell-model vortex beams propagating through oceanic turbulence,” Opt. Express 22(15), 17723–17734 (2014). [CrossRef] [PubMed]
10. D. Liu, Y. Wang, and H. Yin, “Evolution properties of partially coherent flat-topped vortex hollow beam in oceanic turbulence,” Appl. Opt. 54(35), 10510–10516 (2015). [CrossRef] [PubMed]
11. B. Galperin, S. Sukoriansky, N. Dikovskaya, P. L. Read, Y. H. Yamazaki, and R. Wordsworth, “Anisotropic turbulence and zonal jets in rotating flows with a β-effect,” Nonlinear Process. Geophys. 13(1), 83–98 (2006). [CrossRef]
12. M. S. Belen’kii, J. D. Barchers, S. J. Karis, C. L. Osmon, J. M. Brown II, and R. Q. Fugate, “Preliminary experimental evidence of anisotropy of turbulence and the effect of non-Kolmogorov turbulence on wavefront tilt statistics,” Proc. SPIE 3762, 396–406 (1999). [CrossRef]
13. M. S. Belen’kii, E. Cuellar, K. A. Hughes, and V. A. Rye, “Experimental study of spatial structure of turbulence at Maui Space Surveillance Site (MSSS),” Proc. SPIE 6304, 63040U (2006). [CrossRef]
14. V. P. Lukin, “Investigation of some peculiarities in the structure of large scale atmospheric turbulence,” Proc. SPIE 2200, 384–395 (1994). [CrossRef]
15. V. P. Lukin, “Investigation of the anisotropy of the atmospheric turbulence spectrum in the low-frequency range,” Proc. SPIE 2471, 347–355 (1995). [CrossRef]
16. Y. Zhu, M. Chen, Y. Zhang, and Y. Li, “Propagation of the OAM mode carried by partially coherent modified Bessel-Gaussian beams in an anisotropic non-Kolmogorov marine atmosphere,” J. Opt. Soc. Am. A 33(12), 2277–2283 (2016). [CrossRef] [PubMed]
17. I. Toselli, “Introducing the concept of anisotropy at different scales for modeling optical turbulence,” J. Opt. Soc. Am. A 31(8), 1868–1875 (2014). [CrossRef] [PubMed]
18. Y. Li, Y. Zhang, Y. Zhu, and M. Chen, “Effects of anisotropic turbulence on average polarizability of Gaussian Schell-model quantized beams through ocean link,” Appl. Opt. 55(19), 5234–5239 (2016). [CrossRef] [PubMed]
19. Y. Wu, Y. Zhang, Y. Zhu, and Z. Hu, “Spreading and wandering of Gaussian-Schell model laser beams in an anisotropic turbulent ocean,” Laser Phys. 26(9), 095001 (2016). [CrossRef]
20. W. Lu, L. Liu, and J. F. Sun, “Influence of temperature and salinity fluctuations on propagation behaviour of partially coherent beams in oceanic turbulence,” J. Opt. A, Pure Appl. Opt. 8(12), 1052–1058 (2006). [CrossRef]
21. F. Pampaloni and J. Enderlein, “Gaussian, Hermite-Gaussian, and Laguerre-Gaussian beams: A primer,” arXiv:physics/0410021 (2004)
22. O. Korotkova, Random Light Beams Theory and Applications (CRC, 2014).
23. I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, Series and Products, 6th Ed. (Academic, 2000).
24. V. A. Sennikov, P. A. Konyaev, and V. P. Lukin, “Computer simulation of scalar vortex and annular beams LG0L in time-varying random inhomogeneous media,” Proc. SPIE 10035, 10035–10237 (2016).
25. L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media (Washington: SPIE, 2005).
26. X. Sheng, Y. Zhang, F. Zhao, L. Zhang, and Y. Zhu, “Effects of low-order atmosphere-turbulence aberrations on the entangled orbital angular momentum states,” Opt. Lett. 37(13), 2607–2609 (2012). [CrossRef] [PubMed]