BER Analysis of Coherent Free Space Optical Communication Systems with Holographic Modal Wavefront Sensor
 Author: Liu Wei, Yao Kainan, Huang Danian, Cao Jingtai, Wang Liang, Gu Haijun
 Publish: Current Optics and Photonics Volume 1, Issue1, p1~6, 25 Feb 2017

ABSTRACT
Degradation of biterrorrate (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 [14]. 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 ShackHartmann 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 biterrorrate (BER) performance of the FSO links in weak nonKolmogorov 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 lengthr_{0} (D /r_{0} ) was large enough, the Zernike mode was accurate [13, 14]. Li evaluated the performance of the coherent FSO system employing quadrature array phaseshift 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 127subaperture SHWFS, under differentD /r_{0} [1619]. 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 computergenerated holographic elements to wavefront sensing and corrected the response curve [2125]. 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 [2628].In this paper, basedon 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
A _{min}Z_{i} and a reference beam focused at pointA 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 aberrationA _{max}Z_{i} and a reference beam focused at a different pointB , 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,
A_{i} , which is between the maximum and minimum values, The beams will be focused at pointA and pointB , 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
U _{S,0} , the phase function isϕ (r_{s} ) , then the complex amplitude of the signal beam can be expressed as:And the complex amplitude of the reference beam can be expressed as:
Where,
U _{R,0} is the amplitude of the reference beam, andϕ (r_{R} ) 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]:
Where
c.c is a constant, is the conjugate ofU_{S} (r ) ,U_{P} (r_{P} ,r ) is the replay beam, and theG (r_{d} ,r ) is the Green’s Function. In free space,G (r_{d} ,r ) can be approximately expressed as:Where,
k_{d} is the wave vector of laser beam (U_{d} )r_{d} and (x_{d} ,y_{d} ,z_{d} ), (x ,y ,z ) are space coordinates.Assume the Fresnel approximation conditions are satisfied, the thickness of the hologram is ignored, and
λ =λ_{S} =λ_{P} =λ_{d} , (U_{d} )r_{d} can be expressed as [29]:Where, Π is the pupil function, and the range of integration is limited within the pupil.
a andb 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]:
where,
U_{S} (r ,φ ) andU_{LO} (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:where
H_{T} is the magnetic field distribution after frequency mixing,E_{T} ×H_{T} 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]:
where
B is the noise bandwidth of the detector, andI_{LO} 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]:where
P_{S} is the laser signal power. Specially, whenω_{IF} = 0, which is considered as homodyne detection, we can obtain mixing efficiency as [16]: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]:
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:where
N_{p} 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]: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 offaxis 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 peaktovalley (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
N_{p} = 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 rootmeansquare (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.

[FIG. 1.] The theory of HMWFS.

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[FIG. 2.] Computergenerated holograph.

[FIG. 3.] The spots on holographic image plane.

[FIG. 4.] The intensity distribution on focal plane before correction.

[FIG. 5.] The intensity distribution on focal plane before correction.

[TABLE 1.] The BER before and after aberrations correction

[FIG. 6.] The results of the six groups simulation.

[TABLE 2.] The BER of the six groups simulation results