Turbulence-Resistant All Optical Relaying Based on Few-Mode EDFA in Free-Space Optical Systems
Shanyong Cai Zhiguo Zhang Xue Chen


























Journal of  Class Files, Vol. xx, No. xx, xxx xxxx
Shell et al.: Bare Demo of IEEEtran.cls for IEEE Journals


























Turbulence-Resistant All Optical Relaying Based on Few-Mode EDFA in Free-Space Optical Systems
    Shanyong Cai, Member, OSA,
        Zhiguo Zhang,
        and Xue Chen
Shanyong Cai, Zhiguo Zhang and Xue Chen are with State Key Laboratory of Information Photonics and Optical Communications, Beijing University of Posts and Telecommunications, Beijing, 100876, China. e-mail: (caishanyong@bupt.edu.cn, zhangzhiguo@bupt.edu.cn, xuechen@bupt.edu.cn).

Manuscript received xxx xx, xxxx; revised xxx xx, xxxx. This work is supported by National Natural Science Foundation of China (No. 61671076 and No. 61801043) and “the Fundamental Research Funds for the Central Universities”.
    December 30, 2023
====================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================






We investigate the communication performance of a few-mode EDFA based all-optical relaying system for atmospheric channels in this paper. A dual-hop free space optical communication model based on the relay with two-mode EDFA is derived. The BER performance is numerically calculated. Compared with all-optical relaying system with single-mode EDFA, the power budget is increased by 4 dB, 7.5 dB and 11.5 dB at BER = 1×10^-4 under the refractive index structure constant C_n^2 = 2×10^-14, 5×10^-14 and 1×10^-13 respectively when a few mode fiber supporting 4 modes is utilized as the receiving fiber at the destination. The optimal relay location  is slightly backward from the middle of the link. The BER performance is the best when mode-dependent gain of FM-EDFA is zero.





free space optical communication, relay, few mode fiber, EDFA.







§ INTRODUCTION


Recently, free-space optical (FSO) communication systems have attracted a lot of attention for its high bandwidth, unlicensed spectrum and high security <cit.>. However, in FSO system, atmospheric turbulence results in reduced transmission distance and degraded bit error rate (BER). A variety of techniques have been proposed to mitigate atmospheric turbulence fading including spatial diversity <cit.>, wavefront correction <cit.> and relay-assisted optical communication system <cit.>. Relay-assisted transmission system relays the data signal from the source to destination through intermediate nodes (relays). It exploits the fact that the turbulence induced fading variance is distance-dependent and yields improved transmission distance and BER by taking advantage of the shorter hop distances <cit.>. Moreover, relaying systems support optical connection between source and destination which do not have a line of sight.

The current literature on relaying system including serial and parallel electrical relaying and all-optical relaying <cit.>. Compared with electrical relaying where the relays employ optical-to-electrical (OE) and electrical-to-optical (EO) convertors, all-optical relaying avoids the requirement for high-speed electronics and electro-optics. The state-of-art all-optical relaying system is based on single-mode (SM) erbium-doped fiber amplification (EDFA) which only receives and amplifies the power of spatial fundamental mode (Gaussian) <cit.>. Thus, the turbulence induced fading variance at the relay is high which limits the hop distance. In fact, the distorted wavefront after atmospheric channel is a superposition of spatial fundamental mode and spatial high-order modes (e.g. Laguerre-Gaussian). If multiple spatial modes are received and amplified, the turbulence induced fading at the relay will be improved. Recently, spatial-division multiplexing (SDM) optical fiber communication system has attracted a lot of attention <cit.> for its ability to increase system capacity. Few mode fiber (FMF) components like FMFs <cit.> and few-mode erbium-doped fiber amplifications (FM-EDFAs) <cit.> have been fabricated and implemented in SDM optical fiber communication system. FMF components have also been paid attention to in FSO systems <cit.>. For example, experimental results show that the relative standard deviation (RSD) of turbulence fading can be significantly improved by using a FMF receiver compared with single mode fiber (SMF) receiver <cit.>. Moreover, FM-EDFAs which could receive and amplify several spatial modes simultaneously have been used in preamplified receiver in FSO system at ECOC'17 to resist turbulence <cit.>.

In this paper, a scheme of FM-EDFA based all-optical relaying system to resist turbulence for atmospheric channels is proposed which is the extension of our earlier work on CLEO'18 <cit.>. The relay is composed of a two-mode EDFA with two mode fiber (TMF) pigtails. At the relay, the TMF pigtail is utilized as a receiver leading to the improvement of the turbulence induced fading. After FM-EDFA, the mixed spatial modes field is equally amplified and forward to the destination. Compared with previous SM-EDFA based relaying system, the simulation resuls show that system BER is optimized for both fixed and variable gain amplification under atmospheric turbulence.



§ DUAL-HOP FSO MODEL BASED ON FM-EDFA

In this section, we develop model for dual-hop FSO system based on FM-EDFA. Firstly, the scheme is given. Secondly, the signals and noises at the destination terminal for dual-hop FSO communication systems are obtained according to the general approach introduced in [10].





 §.§ Scheme 

The scheme of FM-EDFA based dual-hop transmission system is shown in Fig. 1(a). The relay is composed of receiving lens, FM-EDFA and transmitting lens (Fig. 1(a,c)). At the input end of relay, the signal and background light is coupled into FMF pigtail of FM-EDFA. In the remainder of this paper, assuming FM-EDFA only supports two modes (LP_01 mode and degenerate LP_11 modes). Assuming the mode dependent gain of FM-EDFA is zero, the noisy spatial mode fields are equally amplified by FM-EDFA and forwarded to the destination. A FMF which supports N modes (N=1, 2 or 4) is placed on the destination to collect light field. The optical links between source and relay and between relay and destination are atmospheric turbulence channels. The light field is distorted by atmospheric turbulence in the links.

The fiber specifications and link parameters are shown in Fig. 1(d-e). The source-destination link length is 5.0 km. The relay is placed in the middle of the link (d_1 = d_2). SMF at Source/Relay or FMF at Relay/Destination is placed on the focal point of transmiting/receiving lens. SMF is commercially available SMF-28 fiber. The mode field diameter of SMF-28 is 10.4 μ m  (Fig. 1(d)). FMF is taken to be graded-index fiber with profile parameter ρ=2, thus the spatial distribution of core mode in the FMF is given by

    ψ_m,n( r,θ) = C*L_n - 1^m( r^2/ω ^2)*( r/ω)^m*exp( - r^2/2ω ^2)
    {[ cos(mθ ); sin(mθ ) ].,

where C is the amplitude constant, L_n^m are the associated Laguerre polynomials, and ω is a constant relating modal diameter. Here, ω = 3.89 μ m is set, so the mode field diameter of LP_01 mode is 11 μ m which is equal to that of OFS' FMFs  (Fig. 1(d)). For relaying system with FM-EDFA, the transmitting lens focal lengths of the first and second hop are f_t^1 = 20 cm and f_t^2 = 21.15 cm respectively (Fig. 1(a,e)) to ensure that the diameter of the output beams are the same. The receiving lens focal lengths of the first and second hop are equal to 40 cm (f_r^1 = f_r^2 =40 cm) (Fig. 1(e)). In addition, in relaying system with SM-EDFA, the receiving lens focal lengths of the first hop is set to 38 cm for fair comparison with relaying system with FM-EDFA (Fig. 1(e)). At this time, the coupling efficiencies from free space to LP_01 mode in SMF/FMF are the same and all equal to 44% without atmospheric turbulence.

Compared with SM-EDFA based all-optical relay (Fig. 1(b)), FM-EDFA based relay characterizes lower turbulence fading RSD. In analogy with aperture averaging, the reason can be attributed to mode averaging where multiple spatial modes contained in the turbulence induced distorted light field are received simultaneously. In addition, the average receiving power at the relay with FM-EDFA is higher than that with SM-EDFA. In moderate to strong turbulence, the average receiving power can be increased by 4 dB for TMF compared with SMF <cit.>. After SM-EDFA or FM-EDFA, the fixed amplification gain of EDFA keeps the average relay output power constant at P_r.



 §.§ Model




  §.§.§ Channel Model


The considered system operates at wavelength λ = 1550 nm with on-off keying (OOK) and direct detection (DD). The temporal random atmospheric channel impairments are denoted by h_i(t) for the ith hop. The channel impairments h_i(t) usually can be formulated as h_i(t)≈√(L_i) h_p^iη_i. L_i = α_attnd_i denotes the atmospheric attenuation caused by absorption and scattering effects where α_attn and d_i represent the attenuation coefficient and hop distance, respectively. It can be considered constant since the weather changes slowly over a long time. h_p^i denotes impairment caused by the spatial spread of the optical beam and also can be considered constant for a given transmitting beam diameter and link distance. η_i denotes turbulence-induced fading factor and is a random coefficients. In this paper, to obtain h_i when FMF is utilized as receiving fiber in the two hops, atmospheric turbulence is simulated by phase plates with appropriate randomness.

Turbulence phase plates are generated via Fourier Transform method <cit.>. The basic idea is using atmospheric phase spectrum according to atmospheric turbulence theory to filter Gaussian random noises. The phase screen in x-space is calculated by

    φ(x,y)=IFFT(C·σ(k_x,k_y))


    σ^2(k_x,k_y)=0.023r_0^-5/3[k_x^2+k_y^2]^-11/6,

where σ^2(k_x,k_y) is the variance of phase spectrum, and C is a matrix of complex Gaussian random numbers. r_0 is the atmospheric coherence length, which denotes the turbulence strength and is a function of refractive index structure constant, link distance and optical beam. For Gaussian beam, r_0 can be approximated by the simple expression <cit.>

    r_0 = [8/3(a+0.62Λ^11/6)]^3/5(0.423C_n^2k^2d)^-3/5,    l_0 ≪ r_0 ≪ L_0,

where C_n^2 is the refractive index structure constant, k is the wave number, d is link distance, a and Λ are the parameters related to the beam waist, divergence angle of transmitting Gaussian beam and transmission distance (d). The detailed derivation of a and Λ can be referred to reference <cit.>. In this paper, D/r_0 is used to characterize the effect of turbulence on FSO beams where D=38 mm is the transmitting beam diameter in our scheme. Source-destination length d_sd is used to calculate D/r_0.





  §.§.§ Dual-Hop Transmission Model Based on FM-EDFA


A diagram of dual-hop model with FM-EDFA is shown in Fig. 2. For relay, at the input port of FM-EDFA, the total received electric field is composed of three spatial modes (including two degenerate modes) and is the sum of the signal electric field and the background radiation electric field

    E_r(r,θ,t)=√(2xP_t)∑_m=1^3{h_1,1me^jϕ_1,1me^jω_0tψ_m(r,θ)}
     +√(2N_bδν)∑_m=1^3∑_l=-M^M{e^jϕ_l,me^jω_ltψ_m(r,θ)},

where P_t and x∈{0,2} denote the transmit power of the source and the transmitted OOK symbol. h_1,1me^jϕ_1,1m represent fading coefficients from transmitting LP_01 mode to the mth fiber mode of FM-EDFA in the first hop (m=1 corresponds to LP_01, m=2 corresponds to LP_11a and m=3 corresponds to LP_11b, respectively.). N_b = P_b/B_o denotes the background radiation power spectral density, where P_b and B_o are the background radiation power of single  mode and the optical bandwidth of the system, respectively. The background radiation noise is modeled as a white Gaussian noise in space, thus, values of N_b corresponding to three spatial modes are the same. Following the approach in <cit.>, the electric field corresponding to the background radiation in (5) is written as a sum of 2M+1 (M = B_o/(2δν)) cosine terms at frequencies ω_l = ω_0 +2π lδν with center frequency ω_0 =2πν_0 (ν_0: optical center frequency), where δν denotes the spacing of the considered frequencies. ϕ_l,m denotes a random phase corresponding to the mth mode. ψ_m(r,θ) is the spatial complex amplitude of the mth fiber mode.

FM-EDFA is followed to amplify the received electric field. ASE noise is introduced after FM-EDFA. For amplification, a fixed gain or a variable gain can be employed. The fixed amplification gain keeps the average output power of the relay limited at P_r and is expressed as:

    G_fixed=P_r+3n_sphν_0B_o/P_tE[α_1]+3P_b+3n_sphν_0B_o,

where E[α_1] is the average channel impairment of first hop with α_1 = ∑_m=1^3h_1,1m^2. 3P_b represents the total background radiation power of three modes. n_sp is the amplifier spontaneous emission factor. h denotes Planck`s constant. The variable amplification gain keeps the output power of relay stabilized at constant value P_r at all times and is given by

    G_var=P_r+3n_sphν_0B_o/P_tα_1+3P_b+3n_sphν_0B_o.

Channel impairment α_1, required in the above equation can be estimated at the relay based on pilot symbols emitted by the source. For the sake of brevity, in the remainder of this paper, G_EDFA represents both G_fixed and G_var of the relay  and supposing P_r = P_t.

The total electric field at the destination is given by

    E_d(r,θ,t) = ∑_nE_d,n(t)ψ_n(r,θ),

which is the superposition of multiple spatial mode fields (ψ_n(r,θ)) supported by the receiving FMF in the destination. The time dependent component E_d,n(t) is given by

    E_d,n(t) = E_s(t) + E_b,r(t) + E_ASE(t) + E_b,d(t)
     = √(2xP_t)|H_1HH_2^T|cos(ω_0t + ϕ_s) 
     + √(2N_bδν G_EDFA)∑_m=1^3∑_l=-M^Mh_2,mncos(ω_lt + ϕ_r,l,mn)
     + √(2N_0δν)∑_m=1^3∑_l=-M^Mh_2,mncos(ω_lt + ϕ_ASE,l,mn)
    +  √(2N_bδν)∑_l=-M^Mcos(ω_lt + ϕ_d,l,n),

which is the add of signal light electric field E_s, background radiation electric field at the relay E_b,r, ASE noise after FM-EDFA E_ASE and background radiation electric field at the destination E_b,d. H_1 = [h_1,11e^jϕ_1,11, h_1,12e^jϕ_1,12, h_1,13e^jϕ_1,13], H_2=[h_2,1ne^jϕ_2,1n, h_2,2ne^jϕ_2,2n, h_2,3ne^jϕ_2,3n] and H_1HH_2^T=∑_m=1^3{h_1,1me^jϕ_1,1m√(G_EDFA)h_2,mne^jϕ_2,mn} where
h_2,mne^jϕ_2,mn denote the fading coefficients from the mth spatial mode at relay to the nth spatial mode at the destination for the second hop (Fig. 2). N_0=n_sp(G_EDFA-1)hν_0 is the spectral density of the ASE noise. ϕ_r,l,mn, ϕ_ASE,l,mn and ϕ_d,l,n in (9) are mutual independent random phase. At the destination, after photodetector (PD), the optical signal is converted into electrical signal. The expression of the photocurrent is as follows:
 
    i_photo(t)=R∑_nE_d,n^2(t),

where R=ρ q/hν_0 denotes the responsivity of the PD, ρ and q=1.6× 10^-19 denote the efficiency of PD and the charge of an electron, respectively. The useful signal terms and unwanted noise terms can be obtained by inserting (9) into (10). Following the approach in reference <cit.>, the expressions of these terms are as follows.

Signal Direct Current term :

    I_s = ∑_nRxP_t|H_1HH_2^T|^2.


Background-Background Direct Current terms :

    I_b× b,r = ∑_nRN_bB_oG_EDFA∑_m=1^3h_2,mn^2,



    I_b× b,d = ∑_nRN_bB_o,

where I_b× b,r and I_b× b,d denote the background-background direct current at relay and destination, respectively.

Background-Background Beat Noise:

    σ_b× b,r^2 = ∑_nR^2N_b^2G_EDFA^2(2B_eB_o-B_e^2)
    (∑_m=1^3h_2,mn^4 + 4∑_i=1^3∑_j=1,j≠ i^3h_2,in^2h_2,jn^2),



    σ_b× b,d^2 = ∑_nR^2N_b^2(2B_eB_o-B_e^2),



    σ_b,r× d^2 = ∑_n4R^2N_b^2G_EDFA(∑_m=1^3h_2,mn^2)(2B_eB_o-B_e^2),

where σ_b× b, r^2, σ_b× b, d^2 and σ_b,r× d^2 denote the variance of the relay background-background radiation beat noise, the destination background-background radiation beat noise and the relay-destination background beat noise, respectively. B_e is the electronic bandwidth of the system. The total background-background radiation beat noise variance is given by

    σ_b× b^2 = σ_b× b, r^2+σ_b× b, d^2+σ_b,r× d^2.


ASE-ASE Direct Current term :

    I_ASE× ASE = ∑_nRN_0B_o(∑_m=1^3h_2,mn^2).


ASE-ASE Beat Noise:

    σ_ASE× ASE^2 = ∑_nR^2N_0^2(2B_eB_o-B_e^2)
    (∑_m=1^3h_2,mn^4 + 4∑_i=1^3∑_j=1,j≠ i^3h_2,in^2h_2,jn^2).


Signal-Background Beat Noise:

    σ_s× b,r^2 = ∑_n 4R^2xP_t|H_1HH_2^T|^2G_EDFA(∑_m=1^3h_2,mn^2)N_bB_e,



    σ_s× b,d^2 = ∑_n 4R^2xP_t|H_1HH_2^T|^2N_bB_e,

where σ_s× b,r^2 and σ_s× b,d^2 denote the variance of beat noise between signal and relay background radiation and that of beat noise between signal and destination background radiation, respectively. 

Signal-ASE Beat Noise:

    σ_s× ASE^2 = ∑_n 4R^2xP_t|H_1HH_2^T|^2(∑_m=1^3h_2,mn^2)N_0B_e.


Background-ASE Beat Noise:

    σ_b× ASE^2 = ∑_n 4R^2N_bN_0(2B_oB_e-B_e^2)(∑_m=1^3h_2,mn^2)
    (1+G_EDFA∑_m=1^3h_2,mn^2),

which is composed of relay background-ASE beat noise and destination background-ASE beat noise.

Shot Noise:

    σ_shot, off^2 = 2q(I_b× b,r+I_b× b,d+I_ASE× ASE)B_e,


    σ_shot, on^2 = σ_shot, off^2 + 2qI_sB_e,

where σ_shot, off^2 and σ_shot, on^2 denote the shot noise variances for transmit symbol x=0 and x=2, respectively.

In addition to the noise terms discussed above, the impact of thermal noise of PD at the destination is also important. Finally, the total noise variance of relaying system is

    σ_off^2 = σ_th^2 + σ_shot, off^2 + σ_b× b^2 +σ_ASE× ASE^2 +σ_b× ASE^2,

for transmit symbol x=0, and

    σ_on^2 = σ_off^2 + 2qI_sB_e + σ_s× b,r^2 + σ_s× b,d^2 + σ_s× ASE^2,

for transmit symbol x=2 with thermal noise variance σ_th^2 = 4KTB_e/R_L. K is the Boltzmann constant, T is the temperature in Kelvin, and R_L is the photodetector load resistance. In general, the model for dual-hop transmission system based on FM-EDFA is composed of signal current I_s, in (11) and the noise variances in (26) and (27). In the remainder of this paper, the performance comparison between dual-hop transmission system with FM-EDFA and dual-hop transmission system with SM-EDFA will be given where the model for dual-hop transmission system based on SM-EDFA can be obtained when only LP_01 mode is considered at the relay.



§ NUMERICAL RESULTS

To evaluate the system performance, we present simulation results for BER of dual-hop all-optical relaying system with FM-EDFA. At first, the system parameters of the relaying system is given in Table 1. The system operates at bit-rates (BR) 2 Gpbs. The BER is averaged over 1000 different fading states with 2 × 10^6 bits transmitted per fading state.



In the presence of atmospheric turbulence, atmospheric fading coefficients h_i are time-varying. With the aid of the 1000 turbulence screens generated by Fourier Transform method, the 1000 values of h_i are calculated over N=1000 fading states. Considering the refractive index structure constant C_n^2=1× 10^-13 and no atmospheric attenuation, the statistical distribution of 1000 values of α_1 (α_1 = ∑_m=1^3h_1,1m^2 denotes the atmospheric fading at the relay) is displayed in Fig. 3. Obviously, the average value of fading at the relay with FM-EDFA (E(α_1, FM)=9.84 dB) is smaller than that of the relay with SM-EDFA (E(α_1, SM)=12.5 dB). The atmospheric fading at the relay with FM-EDFA shows a narrower distribution (RSD is restrained by 69%) which is expected to bring improved BER for all optical relaying system.














A FMF which supports N modes (N=1, 2 or 4) is placed on the destination. For N=1, the BER for dual-hop relaying system with SM-EDFA or with FM-EDFA is given in Fig. 4-6. Fig. 4 shows the BER vs transmitting power (P_t) under fixed and variable gain with C_n^2=2× 10^-14 (D/r_0 = 1.093) and α_attn=0.43dB/km. With the increase of transmitting power, the BER of fixed gain relaying system with FM-EDFA decreases faster than that of relaying system with SM-EDFA. The power budget is increased by 0.6 dB at BER = 1× 10^-4 compared with fixed gain relaying system with SM-EDFA. The reason is mainly due to the increase of signal-dependent noise. When signal-dependent noise is the major noise, the SNR of two relaying system is optimized. At this time, for the turbulence-resistant effect of FM-EDFA, BER performance of fixed gain relaying system with FM-EDFA is better.

With the increase of turbulence level, the growth of power budget for fixed gain relaying system with FM-EDFA is more obvious compared to the fixed gain relaying system with SM-EDFA. Fig. 5 shows the BER performance considering the refractive index structure constant C_n^2=5× 10^-14 (D/r_0 = 1.895) and attenuation coefficient α_attn=4.2 dB/km. Obviously, the BER of relaying system with FM-EDFA is always better. The power budget is increased by 1.8 dB at BER = 1× 10^-4. Furthermore, the power budget is increased by 2 dB at BER = 1× 10^-4 under C_n^2=1× 10^-13 (D/r_0 = 2.872) and α_attn=0.43 dB/km (Fig. 6). From Fig. 4 to Fig. 6, with the increasing intensity of turbulence, the turbulence compensation effect of FM-EDFA is becoming more and more obviously. In general, the stronger the turbulence is, the better the performance of FM-EDFA based relaying system is than the relaying system with SM-EDFA.

Fig. 4-6 also depict the BER of variable gain relaying system with SM-EDFA or with FM-EDFA as a function of transmitting power. The BER performance of the variable amplification gain is always outperform that of the fixed amplification gain. Different from the fixed gain system, the BER of variable gain relaying system with SM-EDFA is superior to that of variable gain relaying system with FM-EDFA under weak to moderate turbulence (D/r_0 = 1.093, D/r_0 = 1.895 and D/r_0 = 2.872)(Fig. 4-6).

To comment on whether the mode-dependent gain (MDG) of FM-EDFA have an impact on the system performance, MDG of FM-EDFA is taken into account for fixed gain amplification relaying system in Fig. 7. The performance is the best when MDG = 0 and a little bit worse when MDG = 3 dB (the gain coefficient of LP_01 mode is twice that of LP_11 mode).





To analyze the optimal relay placement of dual-hop relaying system with FM-EDFA. Fig. 8 shows the BER performance corresponding to different relay locations d_1 = 1 km, 2.5 km, 3 km and 4 km. d_1 is the link lengths of the first hop.  The total source-destination length is 5 km (d_1 + d_2 = 5 km). The best BER performance occurs at d_1 = 3 km under atmospheric turbulence and haze (C_n^2=5× 10^-14, α_attn=4.2dB/km) when N=1. For fair comparison, the focal lengths of receiving lenses of relaying system with different relay locations are optimized to achieve the maximum coupling efficiency.


Using FMF supporting multiple spatial modes at the destination can further increase the system performance. Fig. 9 and Fig. 10 give the BER for dual-hop relaying system for N=2 and N=4, respectively. Similar to the relaying system with N=1, the BER performance of the fixed gain relaying system with FM-EDFA is superior to that of the fixed gain relaying system with SM-EDFA. For N=2, the power budget is increased by 2.2 dB, 6.3 dB and 12 dB at BER = 1× 10^-4 under the refractive index structure constant C_n^2 = 2× 10^-14, 5× 10^-14 and 1× 10^-13, respectively. For N=4, the power budget is increased by 4 dB, 7.5 dB and 11.5 dB at BER = 1× 10^-4 under the refractive index structure constant C_n^2 = 2× 10^-14, 5× 10^-14 and 1× 10^-13, respectively. Under the same C_n^2, at BER = 1× 10^-4, the increase of power budget for FM-EDFA based relaying systems with N=2 or N=4 are larger than that with N=1. Moreover, when a four mode fiber is placed on the destination, the BER of the variable gain system with FM-EDFA is also always better than that of the variable gain system with SM-EDFA (Fig. 10) which is different from the relaying system with N=1 (Fig. 4-6). In addition, under the same turbulence level, the bigger the N-value, the smaller the transmitting power required for the same BER.



§ CONCLUSION

A FM-EDFA based all-optical relaying system is proposed in this paper to resist atmospheric turbulence. Atmospheric turbulence is modeled with turbulence phase screen. At the relay, using FM-EDFA as the receiver leads to the improvement of the turbulence induced fading which is expected to bring improved BER for all optical relaying system. A model for dual-hop FSO system with FM-EDFA is derived and the BER performance of FM-EDFA based relaying system is given. Compared with relaying system based on SM-EDFA, the BER of FM-EDFA based relaying system is always superior to that of SM-EDFA based relaying system for the fixed gain amplification. When a FMF supporting 2 or 4 spatial modes is placed at the destination, the BER performance of relaying system will be further improved. Under moderate turbulence D/r_0 = 2.872 (C_n^2 = 1× 10^-13), at BER = 1× 10^-4, the power budget of the fixed gain relaying system with FM-EDFA is increased more than 10 dB than that of the fixed gain relaying system with SM-EDFA. In general, the proposed FM-EDFA based all-optical relaying systems have better communication capacity for its turbulence suppression ability and have potential to be used in future last-mile FSO systems.













































1

khalighi-icst-2014 M. A. Khalighi and M. Uysal, “Survey on free space optical communication: A communication theory perspective,” IEEE Commun. Surv. Tut., vol. 16, pp. 2231-2258, 2014.

zhu-IEEtc-2002 Z. Xiaoming and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun., vol. 50, pp. 1293-1300, 2002.

navidpour-IEEETWC-2007 S. M. Navidpour, M. Uysal, and M. Kavehrad, “BER performance of free-space optical transmission with spatial diversity,” IEEE Trans. Wireless Commun., vol. 6, pp. 2813-2819, 2007.

love-ao-1997 G. D. Love, “Wave-front correction and production of Zernike modes with a liquid-crystal spatial light modulator,” App. Opt., vol. 36, pp. 1517-1524, 1997.

safari-IEEETWC-2008 M. Safari and M. Uysal, “Relay-assisted free-space optical communication,” IEEE Trans. Wireless Commun., vol. 7, pp. 5441-5449, 2008.

karimi-jlt-2009 M. Karimi and M. Nasiri-Kenari, “BER analysis of cooperative systems in free-space optical networks,” J. Lightwave Technol., vol. 27, pp. 5639-5647, 2009.

rjeily-IEEETC-2011 C. Abou-Rjeily and A. Slim, “Cooperative Diversity for Free-Space Optical Communications: Transceiver Design and Performance Analysis,” IEEE Trans. Commun., vol. 59, pp. 658-663, 2011.

karimi-jlt-2011 M. Karimi and M. Nasiri-Kenari, “Free space optical communications via optical amplify-and-forward relaying,”  J. Lightwave Technol., vol. 29, pp. 242-248, 2011.

kazemlou-jlt-2011 S. Kazemlou, S. Hranilovic, and S. Kumar, "All-optical multihop free-space optical communication systems," J. Lightwave Technol., vol. 29, pp. 2663-2669, 2011.

bayaki-IEEETC-2012 E. Bayaki, D. S. Michalopoulos, and R. Schober, "EDFA-based all-optical relaying in free-space optical systems," IEEE Trans. Commun., vol. 60, pp. 3797-3807, 2012.

kashani-ICL-2012 M. A. Kashani, M. M. Rad, M. Safari, and M. Uysal, “All-optical amplify-and-forward relaying system for atmospheric channels,” IEEE Commun. Lett., vol. 16, pp. 1684-1687, 2012.

nor-jlt-2017 N. A. M. Nor, Z. Ghassemlooy, J. Bohata, P. Saxena, M. Komanec, S. Zvanovec, et al., “Experimental investigation of all-optical relay-assisted 10 Gb/s FSO link over the atmospheric turbulence channel,” J. Lightwave Technol., vol. 35, pp. 45-53, 2017.

kashani-jocn-2013 M. A. Kashani, M. Safari, and M. Uysal, “Optimal relay placement and diversity analysis of relay-assisted free-space optical communication systems,” J. Opt. Commun. Netw., vol. 5, pp. 37-47, 2013.

ruiz-jocn-2018 R. Boluda-Ruiz, A. Garca-Zambrana, C. Castillo-Vzquez, B. Castillo-Vzquez, and S. Hranilovic, “Amplify-and-forward strategy using MRC reception over FSO channels with pointing errors,” J. Opt. Commun. Netw., vol. 10, pp. 545-552, 2018.

richardson-nature-2013 D. Richardson, J. Fini, and L. Nelson, “Space-division multiplexing in optical fibres,” Nat. Photonics, vol. 7, pp. 354-362, 2013.

li-aop-2014 G. Li, N. Bai, N. Zhao, and C. Xia, “Space-division multiplexing: the next frontier in optical communication,” Adv. Opt. Photonics, vol. 6, pp. 413-487, 2014.

sillard-ofc-2016 P. Sillard, “Few-mode fibers for space division multiplexing,” in Optical Fiber Communication Conference, Anaheim, California, 2016, p. Th1J.1.

bai-oe-2012 N. Bai, E. Ip, Y.-K. Huang, E. Mateo, F. Yaman, M.-J. Li, et al., “Mode-division multiplexed transmission with inline few-mode fiber amplifier,” Opt. Express, vol. 20, pp. 2668-2680, 2012.

ozdur-ol-2013 I. Ozdur, P. Toliver, A. Agarwal, and T. K. Woodward, “Free-space to single-mode collection efficiency enhancement using photonic lanterns,” Opt. Lett., vol. 38, pp. 3554-3557, 2013.

arikawa-ecoc-2016 M. Arikawa, T. Ishikawa, K. Hosokawa, S. Takahashi, Y. Ono, and T. Ito, “Mitigation of fading caused by atmospheric turbulence with FMF coupling and maximum ratio combining used in 320-m free-space optical transmission of 10 Gb/s BPSK,” in ECOC, 2016, pp. 1205-1207.

zheng-oe-2016 D. Zheng, Y. Li, E. Chen, B. Li, D. Kong, W. Li, et al., “Free-space to few-mode-fiber coupling under atmospheric turbulence,” Opt. Express, vol. 24, pp. 18739-18744, 2016.

huang-ecoc-2017 Bin Huang et al., “Turbulence-resistant free-space optical communication using few-mode preamplified receivers,” in ECOC, 2017, paper Tu.2.E.4.

geisler-spie-2018 D. J. Geisler, T. M. Yarnall, G. Lund, C. M. Schieler, M. L. Stevens, N. K. Fontaine, et al., “Experimental comparison of 3-mode and single-mode coupling over a 1.6-km free-space link,” in Proc. SPIE, Free-Space Laser Communication and Atmospheric Propagation, Vol. 10524, 105240H, 2018.

zheng-oe-2018 D. Zheng, Y. Li, H. Zhou, Y. Bian, C. Yang, W. Li, et al., “Performance enhancement of free-space optical communications under atmospheric turbulence using modes diversity coherent receipt,” Opt. Express, vol. 26, pp. 28879-28890, 2018.

cai-cleo-2018 Shanyong Cai, Zhiguo Zhang, Xue Chen, “Few-mode EDFA based all-optical relaying system for atmospheric channels,” in CLEO, 2018, San Jose, California, p. JTu2A.37.

xu-cleo-2014 X. Tong et al., “Simulating atmospheric turbulence using a spatial light modulator based on fourier transform,” in CLEO, 2014, San Jose, California, p. SM4J.3.

andrews-2005 L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media, 2nd ed.1em plus 0.5em minus 0.4emWA Bellingham:SPIE 2005.

olsson-jlt-1989 N. A. Olsson, “Lightwave systems with optical amplifiers,” J. Lightwave Technol., vol. 7, pp. 1071-1082, 1989.



















Shanyong Cai
received his Ph.D. degree in information and communication engineering from Beijing University of Posts and Telecommunications (BUPT), Beijing, China, in 2017 and the M.S. degree in optics from South China Normal University, Guangzhou, China, in 2013. He currently holds a postdoctoral position with the Institute of Information Photonics and Optical Communications, BUPT. His research interests include space-division multiplexing (SDM), free space optical communication (FSO) and orthogonal frequency-division multiple access (OFDMA) PON.


Zhiguo Zhang
received the B.S. degree from Shandong University, Shandong, China, in 2002 and the Ph.D. degree from Beijing University of Posts and Telecommunications (BUPT), Beijing, China, in 2007. He is now an associate professor of Institute of Information Photonics and Optical Communications, BUPT. His main research interests include free space optical communication, optical fiber sensor and optical access networks.


Xue Chen
received the B.S. degree from Dalian University of Technology, Dalian, China, in 1982 and the M.S. degree from Beijing University of Posts and Telecommunications (BUPT), Beijing, China, in 1985. She is now a professor of Institute of Information Photonics and Optical Communications, BUPT. Her main research interests focus on backbone optical transmission and optical access networks.

















