BER Analysis of Coherent Free Space Optical Communication Systems with Holographic Modal Wavefront Sensor

  • cc icon
  • ABSTRACT

    Degradation of bit-error-rate (BER), caused by atmospheric turbulence, seriously hinders the performance of coherent Free Space Optical (FSO) communication systems. An adaptive optics system proves to be effective in suppressing the atmospheric turbulence. The holographic modal wavefront sensor (HMWFS) proposed in our previous work, noted for its fast detecting rates and insensitivity to beam scintillation, is applied to the coherent FSO communication systems. In this paper, based on our previous work, we first introduce the principle of the HMWFS in brief and give the BER of the coherent FSO with homodyne detection in theory, and then analyze the improvement of BER for a coherent FSO system based on our previous simulation works. The results show that the wavefront sensor we propose is better for weak atmospheric turbulence. The most obvious advantages of HMWFS are fast detecting rates and insensitivity to beam scintillation.


  • KEYWORD

    BER , Coherent FSO communication , Holographic modal wavefront sensor , Aberrations correction

  • I. INTRODUCTION

    For the free space optical communication (FSO) system, the performance of the coherent detection scheme is better than that of the intensity modulation direct detection, with higher spectral efficiencies and data rates, a greater ability to decrease both background and thermal noise, and more sensitive receivers [1-4]. However, the atmospheric turbulence greatly degrades the performance of the coherent FSO systems [5]. The wavefront distortion, scintillation, beam wandering and spreading caused by a turbulent atmosphere will not only degrade the entrance efficiency of the receiving antenna but also cause mismatch of the field of the signal beam and the local oscillator (LO) [6].

    An adaptive optics (AO) system is identified as an effective way to compensate atmospheric turbulence in a coherent FSO system [7, 8]. And as the main component of an AO system, the wavefront sensor is attracting extensive attention [9]. The Shack-Hartmann wavefront sensor (SHWFS) is widely used for turbulence compensation for a coherent FSO system [10]. Belmonte put emphasis on elucidating how the addition of AO to the transmitter or receiver can reduce the effects of atmospheric propagation and on quantifying the improvement on the performance of FSO systems regarding coherent detection [11, 12]. Zuo investigated the bit-error-rate (BER) performance of the FSO links in weak non-Kolmogorov turbulence and showed that BER decreased sharply with AO corrections. Considering the influence of both the amplitude fluctuation and spatial phase aberrations, when the ratio of receiving aperture diameter D to the coherent length r0 (D/r0) was large enough, the Zernike mode was accurate [13, 14]. Li evaluated the performance of the coherent FSO system employing quadrature array phase-shift keying modulation over the maritime atmosphere with AO system based on SHWFS for atmospheric turbulence compensation [15]. C. Liu and J. Huang analyzed the improvement of BER for the coherent FSO system with AO system based on SHWFS through numerical simulation and experimental data of the 1.8 m telescope with a 127-subaperture SHWFS, under different D/r0 [16-19]. However, the SHWFS was so sensitive to laser scintillation that the performance was limited seriously for a coherent FSO system. The focal plane wavefront sensor solved this problem. It had higher power usage ratio, but it was much more costly in time [20]. In 2000, Ghebremichael et al. applied computer-generated holographic elements to wavefront sensing and corrected the response curve [21-25]. Changhai et al. used the plane wave in combination with the Fourier lens for the replacement of spherical waves to improve the efficiency of energy utilization [26-28].

    In this paper, based-on our previous work, HMWFS is applied to the coherent FSO systems to evaluate the BER performance improvement. First, the principle of HMWFS and the BER of the coherent FSO communication are given in brief in Section 2. Then, the numerical simulation is used to verify the feasibility of the HMWFS for the coherent FSO system, and the BER improvement for the coherent FSO system with HMWFS based on our previous works are presented in Section 3. Finally, we summarize and present our conclusions in Section 4.

    II. THEORETICAL ANALYSIS

       2.1 The holographic modal wavefront sensing

    The theoretical layout of HMWFS is shown in Fig. 1 [29].

    An object beam with the minimum expected amplitude of the aberration AminZi and a reference beam focused at point A on a distant detector record one hologram, as shown in Fig. 1(a). Another hologram is recorded by an object beam with the maximum amplitude of the aberration AmaxZi and a reference beam focused at a different point B , as shown in Fig. 1(b). Then a multiplexed hologram is generated with the pair of holograms [29].

    The hologram will be reconstructed with an input beam with some amplitude, Ai , which is between the maximum and minimum values, The beams will be focused at point A and point B , as shown in Fig. 1(c). The relative intensity of at these points is related to the coefficient of Zernike item. Thus, after a calibration, the amplitude can be obtained simply by sensing the spot intensities.

    In addition, the multiplexed hologram with more pairs of holograms is needed to increase the number of aberration modes detected. Considering the pairs of spots are spatially separate, the amplitudes of these spots can be read out in parallel by a sensor array. Thus, the Zernike modal amplitude can be obtained through comparing the relative intensities of pairs of spots [29].

    According to Fourier optics, suppose the amplitude of the signal beam is US,0 , the phase function is ϕ(rs ) , then the complex amplitude of the signal beam can be expressed as:

    image

    And the complex amplitude of the reference beam can be expressed as:

    image

    Where, UR,0 is the amplitude of the reference beam, and ϕ(rR ) is the phase function of the reference beam.

    The two waves interfere on the recording medium, and generate a holographic pattern. When a replay beam with random aberrations illuminates the holographic pattern, the wavefront of the signal beam can be reconstructed on the focal plane. The diffraction field distribution on the focal plane shall be expressed as [29]:

    image

    Where c.c is a constant, is the conjugate of US( r) , UP(rP, r) is the replay beam, and the G(rd, r) is the Green’s Function. In free space, G(rd, r) can be approximately expressed as:

    image

    Where, kd is the wave vector of laser beam (Ud)rd and ( xd , yd , zd ), ( x, y, z) are space coordinates.

    Assume the Fresnel approximation conditions are satisfied, the thickness of the hologram is ignored, and λ =λS =λP =λd , (Ud)rd can be expressed as [29]:

    image

    Where, Π is the pupil function, and the range of integration is limited within the pupil. a and b are respectively the length and width of the hologram. F{.} is the Fourier transform. Thus, the Fresnel integral is simplified into a Fourier transform with thin hologram conditions. Thus, the distribution of complex amplitude on the holographic plane can be solved by a numerical method.

       2.2 The BER of the coherent FSO system

    As mentioned in our previous work, for a coherent FSO system, the received beam is mixed with an LO signal. The received beam and the LO signal are given by [16]:

    image
    image

    where, US(r, φ) and ULO(r, φ) are the amplitudes of received beam and LO signal. ψS(r, φ) and ψLO(r, φ) are the variable phases related to spatial position of the laser carrier and LO signal, respectively. ωS and ωLO are their angular frequencies, respectively. The current after frequency mixing is given by:

    image

    where HT is the magnetic field distribution after frequency mixing, ET × HT is the Poynting vector, η is the quantum efficiency of the detector, c is constant, and σ is the area of the detector.

    Since the intensity of the LO signal is much larger than the intensity of the received beam, the detector noise is mainly shot noise. Thus, the mean square value of the noise current is given by [16]:

    image

    where B is the noise bandwidth of the detector, and ILO is the direct current coursed by the LO signal. The signal to noise ratio (SNR) is defined by the ratio of the mean square value of the intermediate frequency signal current and the mean square value of the noise current [16]:

    image

    where PS is the laser signal power. Specially, when ωIF = 0, which is considered as homodyne detection, we can obtain mixing efficiency as [16]:

    image

    According to the Eq. (12), the mixing efficiency is related to both the amplitude of received beam and LO signal, and spatial phase distribution difference. Assuming the received beam and the LO signal are plane waves with uniform intensity, the mixing efficiency is influenced directly by the spatial phase distribution difference. The AO technology is generally regarded as an effective way to decrease the spatial phase distribution difference between the carrier signal and the LO signal.

    In the coherent detection system, the BER can be given by [16]:

    image

    where erfc is the complementary error function, . For synchronous binary phase shift keying (BPSK) modulation, the optical power of the received beam is expressed as:

    image

    where Np is the number of photons received within a single bit. For the BPSK modulation, the BER of the homodyne detection can be given as [16]:

    image

    III. NUMERICAL SIMULATION

    Based on our previous work, in terms of the computergenerated holograph, the center wavelength of the laser is λ = 632.8 nm, the number of sampling points of the holographic grating is 600×600 pixels, the aperture of holographs is 60 mm, the distance between image plane and holographic plane is 1500 mm, the off-axis distance of the spot on focal plane is larger than 1 cm. The low eight order Zernike aberrations (Z3, Z4, …, Z10) are coded in this simulation. Where the peak-to-valley (PV) value of Z3, Z4, Z5 is λ / 2 , and the PV value of Z6, Z7, Z8, Z9, Z10 is λ / 4 . The holograph we eventually obtained is shown in Fig. 2 [29].

    According to the theory of the HMWFS, the relative intensity distribution on the holographic image plane is shown in Fig. 3 [29]. By comparing the relative intensity on the holographic image plane, the Zernike aberrations are obtained quantitatively.

    The mixing efficiency is approximate to the Strehl ratio of the far field of the received beam for homodyne detection, and it is defined by the ratio of the far field encircled energy and the diffraction limited encircled energy [18]. The intensity distribution on the focal plane before correction (detail view) is shown in Fig. 4.

    Thus, it can be seen that the holographic modal wavefront sensor has smaller residual in low eight orders Zernike aberrations detection. The intensity distribution on focal plane after aberrations correction is shown in Fig. 5.

    Then, we analyze the BER before and after aberrations correction based on Eq. (14) when Np = 10 and η = 1. The results are shown in Table 1.

    We make another five groups numerical simulation by changing the introduced aberrations. In our simulation, we introduce the Zernike mode aberrations into the laser beam before correction. The root-mean-square (RMS) values of the aberrations before correction are between 0.15λ and 0.25λ . And RMS values of the residual aberrations after correction are nearly below 0.1λ [29]. Thus, we get the five images before and after correction as shown in Fig. 6.

    Where, the Fig. 6(a) shows the images before correction, and the Fig. 6(b) shows the images after correction. The BER before and after correction are shown in Table 2.

    As shown in Table 2, according to homodyne detection coherent FSO system, we can see that before correction, the BER is above 10−4, and after correction, the BER is below 10−5. According to the BER before correction, we can estimate the atmospheric turbulence intensity. Under some weak atmospheric turbulence conditions, the BER is about 10-4 before correction. And after AO correction based on HMWFS, the BER can reach below 10−7. In the meanwhile, under some strong atmospheric turbulence conditions, the BER is about 10−2 before correction. And after correction, the BER can achieve below 10−5. Above all, the AO system based on HMWFS is an effective way to improve the BER performance of a coherent FSO system with homodyne detection conditions. According to Table 2, we can also verify that the wavefront sensing we proposed is better for weak atmospheric turbulence with fast detecting rates and insensitivity to beam scintillation.

    IV. CONCLUSION

    In this paper, we analyzed the principle of HMWFS in theory and gave the BER expression of the coherent FSO system, and then first applied HMWFS to a coherent FSO system, and then we analyze the BER performance improvement through numerical simulation based on HMWFS. Finally, we discussed the impact of HMWFS on BER of the coherent FSO system before and after the aberrations correction based on HMWFS. According to the results of numerical simulation, the wavefront sensor we proposed is better for weak atmospheric turbulence. The most obvious advantage of HMWFS is fast detecting rates and insensitivity to beam scintillation.

    The disturbance strength is only described by the value of the introduced aberrations in our work. The relationship between the introduced aberrations and the structure constant of refractive index will be carried out in the near future.

  • 1. Wei L., Kainan Y., Danian H., Xudong L., Liang W., Yaowen L. (2016) “Performance evaluation of coherent free space optical communications with a double-stage fast-steering-mirror adaptive optics system depending on the Greenwood frequency,” [Opt. Express] Vol.24 P.13288-13302 google doi
  • 2. Ma J., Li K., Tan L., Yu S., Cao Y. (2015) “Performance analysis of satellite-to-ground downlink coherent optical communications with spatial diversity over Gamma-Gamma atmospheric turbulence,” [Appl. Opt.] Vol.54 P.7575-7585 google doi
  • 3. Geng D., Du P., Wang W., Gao G., Wang T., Gong M. (2014) “Single laser free-space duplex communication system with adaptive threshold technique and BER analysis in weak turbulent atmosphere,” [Opt. Letters] Vol.39 P.3950-3953 google doi
  • 4. Rennie R. M., Goorskey D., Whiteley M. R., Jumper E. J. (2012) “Wavefront measurements of a laser-induced break-down spark in still air,” [Appl. Opt.] Vol.51 P.2306-2314 google doi
  • 5. Liu W., Shi W., Cao J., Lv Y., Yao K., Wang S., Wang J., Chi X. (2014) “Bit error rate analysis with real-time pointing errors correction in free space optical communication system,” [Optik] Vol.12 P.324-328 google
  • 6. Dongyiel S., Such H. Y., Woo C. J. (2000) “4×10 Gb/s terrestrial optical free space transmission over 1.2 km using an EDFA preamplifier with 100 GHz channel spacing,” [Opt Express] Vol.7 P.280-284 google doi
  • 7. Yang Q., Zhao J., Wang M., Jia J. (2015) “Wavefront sensor-less adaptive opitcs based on the trust region method,” [Opt. Letters] Vol.40 P.1235-1237 google doi
  • 8. Cao J., Zhao X., Li Z., Liu W., Song Y. (2014) “Stochastic parallel gradient descent laser beam control algorithm for atmospheric compensation in free space optical communication,” [Optik] Vol.125 P.6142-6147 google doi
  • 9. Linhai H., Rao C. (2011) “Wavefront sensor-less adaptive optics: a general model-based approach,” [Opt. Express] Vol.19 P.371-379 google doi
  • 10. Ma J., Zhao F., Tan L., Yu S., Yang Y. (2010) “Degradation of single-mode fiber coupling efficiency due to localized wavefront abrerrations in free-space laser communications,” [Opt. Eng.] Vol.49 P.1-6 google
  • 11. Belmonte A., Rodriguez A., Dios F., Comeon A. (2007) “Phase compensation considerations on coherent, free-space laser communications system,” [Proc. SPIE] Vol.6736 P.A-11 google
  • 12. Belmonte A. (2008) “Influence of atmospheric phase compensation on optical heterodyne power measurements,” [Opt. Express] Vol.16 P.6756-6767 google doi
  • 13. Zuo L., Ren Y., Dang A., Hong G. (2010) “Performance of coherent BPSK systems using phase compensation and diversity techniques,” [Proceedings of IEEE Conference on Global Telecommunications] P.1-5 google
  • 14. Zuo L., Dang A., Ren Y., Guo H. (2011) “Performance of phase compensated coherent free space optical communications through non-Kolmogorov turbulence,” [Opt. Commun.] Vol.28 P.1491-1495 google
  • 15. Li M., Cvijetic M. (2015) “Coherent free space optics communications over the maritime atmosphere with use of adaptive optics for beam wavefront correction,” [Appl. Opt.] Vol.54 P.1453-1462 google doi
  • 16. Liu C., Chen S., Li X., Xian H. (2014) “Performance evaluation of adaptive optics for atmospheric coherent laser communications,” [Opt. Express] Vol.22 P.15554-15563 google doi
  • 17. Liu C., Chen M., Chen S., Xian H. (2016) “Adaptive optics for the free-space coherent optical communications,” [Opt. Commun.] Vol.361 P.21-24 google doi
  • 18. Jian H., Ke D., Chao L., Peng Z., Dagang J., Zhoushi Y. (2014) “Effectiveness of adaptive optics system in satellite-to-ground coherent optical communication,” [Opt. Express] Vol.22 P.16000-16007 google doi
  • 19. Huang J., Mei H., Deng K., Kang L., Zhu W., Yao Z. (2015) “Signal to noise ratio of free space homodyne coherent optical communication after adaptive optics compensation,” [Opt. Commun.] Vol.356 P.574-577 google doi
  • 20. Liu W., Shi W., Wang B., Yao K., Lv Y., Wang J. (2013) “Free space optical communication performance analysis with focal plane based wavefront measurement,” [Opt. Commun.] Vol.309 P.212-220 google doi
  • 21. Ghebremichael F., Andersen G. P., Gurley Kenneth S. (2008) “Holography-based wave-front sensing,” [Appl. Opt.] Vol.47 P.A62-A69 google
  • 22. Andersen G. P., Dussan L., Ghebremichael F., Chen K. (2009) “Holographic wave-front sensor,” [Opt. Eng.] Vol.48 P.085801 google
  • 23. Andersen G. P., Ghebremichael F., Gurley K. S. (2007) “Fast Computing-Free Wave-front Sensing,” [Adaptive Optics: Methods, Analysis and Applications] Vol.18 P.AWC4 google
  • 24. Mishra S. K., Bhatt R., Mohan D., Gupta A. K., Sharma A. (2009) “Differential modal Zernike wavefront sensor employing a computer-generated hologram: a proposal,” [Appl. Opt.] Vol.48 P.6458-6465 google doi
  • 25. Bhatt R., Mishra S. K., Mohan D., Gupta A. K. (2008) “Direct amplitude detection of Zernike modes by computer-generated holographic wavefront sensor: Modeling and simulation” [Optics and Lasers in Engineering] Vol.46 P.428-439 google doi
  • 26. Changhai L., Zongfu J. (2009) “Holographic Modal Wave-front Sensor: Theoretical Analysis and Simulation,” [Chinese Journal of Lasers] Vol.36 P.147-152 google
  • 27. Changhai L., Fengjie X., Haotong M., Shengyang H., Zongfu J. (2010) “Modal wave-front sensor based on binary phase-only multiplexed computer-generated hologram,” [Appl. Opt.] Vol.49 P.5117-5124 google doi
  • 28. Changhai L., Fengjie X., Shengyang H., Zongfu J. (2011) “Performance analysis of multiplexed phase computer-generated hologram for modal wave-front sensing,” [Appl. Opt.] Vol.50 P.1631-1639 google
  • 29. Liu W., Shi W., Yao K., Cao J., Wu P., Chi X. (2014) “Fiber Coupling efficiency analysis of free space optical communication systems with holographic modal wavefront sensor,” [Opt. Laser Technol.] Vol.60 P.116-123 google doi
  • [FIG. 1.] The theory of HMWFS.
    The theory of HMWFS.
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [] 
  • [FIG. 2.] Computer-generated holograph.
    Computer-generated holograph.
  • [FIG. 3.] The spots on holographic image plane.
    The spots on holographic image plane.
  • [FIG. 4.] The intensity distribution on focal plane before correction.
    The intensity distribution on focal plane before correction.
  • [FIG. 5.] The intensity distribution on focal plane before correction.
    The intensity distribution on focal plane before correction.
  • [TABLE 1.] The BER before and after aberrations correction
    The BER before and after aberrations correction
  • [FIG. 6.] The results of the six groups simulation.
    The results of the six groups simulation.
  • [TABLE 2.] The BER of the six groups simulation results
    The BER of the six groups simulation results