# Vector Rayleigh Diffraction of High-Power Laser Diode Beam in Optical Communication

Authors: Qiang Xu; Renxian Li; Yuanyuan Zhang; Yiping Han; Zhensen Wu
Journal: International Journal of Optics
 (2020)

---

## Abstract

Laser diodes (LDs) are widely used in optical wireless communication (OWC) and optical networks, and proper theoretical models are needed to precisely describe the complicated beam field of LDs. A novel mathematical model is proposed to describe the vectorial field of nonparaxial LD beams. Laser beam propagation is studied using the vector Rayleigh diffraction integrals, and the stationary phase method is used to find the asymptotic expansion of diffraction integral. The far-field distribution of the LD beam in the plane parallel and perpendicular to the junction is considered in detail, and the computed intensity distributions of the theory are compared with the corresponding measurements. This model is precise for single transverse model beam of LDs and can be applied to describe the LD beams in OWC and optical networks.


---

## Body

## 1. Introduction

Considering that laser diodes (LDs) are the efficient light source and easy to integrate, LD-enabled optical wireless communication (OWC) is an emerging technology for realizing high-confidentiality and high-speed point-to-point (PtP), vehicle-to-vehicle, and white-lighting data access links in free-space communication [1–6], indoor communication [7, 8], underwater communication [9–11], and optical networks [12–16].However, the output beam quality of LDs is relatively poor, such as astigmatism, high beam asymmetry, and large beam divergence [17–19], in many applications, and proper theoretical models are needed to precisely describe the optical field distribution of LDs.The problem of laser propagation is mainly dealt through paraxial approximation. However, the output facet of LDs is extremely small, and their beams are divergent and asymmetrical. The rigid optical field distributions cannot be calculated from the paraxial approximation, and the longitudinal component in beam propagation direction should be considered. Thus, the vector theory for nonparaxial beams should be used to precisely describe beam fields of LDs. Several models, such as exponential Gaussian function [20–22], Hermite–Gaussian model [23], nonparaxial diffraction of vectorial Gaussian wave [24, 25], plane waves with a small aperture [26], propagation of LD beams in the optical system [27–29], and polarization of LD beams [30], are used to describe the beam field of LDs. However, no theoretical model is used for all cases because of the complicated beam field of LDs. Thus, a novel model should be developed to precisely describe the output field of LDs, which is the aim of this paper.


## 2. Vectorial Electric Field of LD Beam

Considering that transverse electric modes are usually excited in LD,Ex,y is identified with the component Eyx,y of the electric field vector, and a source beam is linearly polarized at the plane z = 0:(1)Exx0,y0=0,(2)Eyx0,y0=E0exp−px03/2−qy02,where p and q are related to the waveguide structure of LDs, 1/p=1.22λn/dx×10−3 and 1/q=1.22λn/dy×10−3, in which λn is the beam wavelength in the active layer of LDs, dx is the waveguide width in the x direction, and dy is the waveguide width in the y direction, and E0 is a constant.Beam propagation is governed by the vector Rayleigh diffraction integrals that provide the field expression in the entire half-spacez>0. When the boundary condition at the plane z = 0 is given, the field takes the following form [24, 26]:(3)Exr=−12π∫Exr0∂Gr,r0∂zds,Eyr=−12π∫Eyr0∂Gr,r0∂zds,Ezr=12π∫Exr0∂Gr,r0∂x+Eyr0∂Gr,r0∂yds,where r0=x0i+y0j (r0 is the vector in beam output plane), r=xi+yj+zκ (r is the beam propagation vector), i,j,κ are the unit vectors in the x-, y-, and z-directions, respectively, and(4)Gr,r0=expikr−r0r−r0,in which k is the wavenumber related to wavelength λ by k=2π/λ. Substituting equation (4) into equation (3) yields [24, 26](5)Exx,y,z=−12π∫−∞+∞∫−∞+∞Exx0,y0zikr−r0−1r−r03expikr−r0dx0dy0,Eyx,y,z=−12π∫−∞+∞∫−∞+∞Eyx0,y0zikr−r0−1r−r03expikr−r0dx0dy0,Ezx,y,z=12π∫−∞+∞∫−∞+∞Exx0,y0x−x0ikr−r0−1r−r03expikr−r0+Eyx0,y0y−y0ikr−r0−1r−r03expikr−r0dx0dy0.We expandexpikr−r0 into a series, keeping the first-, second-, and third-order series expansions [31]:(6)expikr−r0≈expikr−ikxx0+yy0r+iky2+z2x02+x2+z2y02−2xyx0y02r3,where r=x2+y2+z2, replace expikr−r0 in equation (5) by equation (6), and replace r−r0 in equation (5) by r:(7)Exx,y,z=0,Eyx,y,z=−12πzikr−1r3expikr∫−∞+∞∫−∞+∞Eyx0,y0×exp−ikxx0+yy0r+iky2+z2x02+x2+z2y02−2xyx0y02r3dx0dy0,Ezx,y,z=12πikr−1r3expikr∫−∞+∞∫−∞+∞Exx0,y0x−x0+Eyx0,y0y−y0×exp−ikxx0+yy0r+iky2+z2x02+x2+z2y02−2xyx0y02r3dx0dy0.For largek (104mm−1), exp−ikxx0+yy0/r+iky2+x2x02+x2+z2y02−2xyx0y0/2r3 rapidly oscillates, and such rapid oscillations over the range of integration indicate that the integrand averages to approximately zero, except near the stationary phase. Thus, the stationary phase method is used to find the asymptotic expansion of the diffraction integral.The corresponding diffraction integral is approximated by [32](8)Ux,y,z=∬Dfx0,y0expikgx0,y0dx0dy0≈2πσkHfxs,ysexpikgxs,ys.fx0,y0=Eyx0,y0 for Eyx,y,z, and fx0,y0=Exx0,y0x−x0+Eyx0,y0y−y0 for Ezx,y,z.where(9)H=∂2g∂x02∂2g∂y02−∂2g∂x0∂y02,(10)σ=1,ifH<0,i,ifH>0,∂2g∂x02xs,ys>0,−i,ifH>0,∂2g∂x02xs,ys<0,where xs and ys are the stationary phase points, and we have(11)gx0,y0=−xx0+yy0r+y2+z2x02+x2+z2y02−2xyx0y02r3.Letting(12)∂g∂x0=0,∂g∂y0=0,we find the stationary phase points(13)xs=r2z2x,ys=r2z2y,(14)gxs,ys=−r2z2x2+y2,and ∂2g/∂x02=y2+z2/r3, ∂2g/∂y02=x2+z2/r3, and ∂2g/∂x0∂y0=xy/r3.Thus,(15)H=∂2g∂x02∂2g∂y02−∂2g∂x0∂y02=z2r4,and H>0,∂2g/∂x02xs,ys>0,(16)σ=i.Substituting equations (13)–(16) into equation (7) yields(17)Exx,y,z=0,Eyx,y,z=−iikr−1krE0⋅exp−pr2z2x3/2−qr2z2y2×expik−r2z2x2+y2z2expikr,Ezx,y,z=iikr−1kr⋅yzx2+y2z2E0⋅exp−pr2z2x3/2−qr2z2y2×expik−r2z2x2+y2expikr.Equation (17) represents the expression of vector theory for nonparaxial LD beam.The intensity profiles can be given by(18)Ixx,y,z=0,Iyx,y,z=1+k2r2k2r2E02⋅exp−2pr2z2x3/2−2qr2z2y2,Izx,y,z=1+k2r2k2r2y2z2x2+y2z22E02⋅exp−2pr2z2x3/2−2qr2z2y2,and the total intensity can be expressed as(19)Ix,y,z=Ixx,y,z+Iyx,y,z+Izx,y,z=1+k2r2k2r21+y2z2x2+y2z22E02exp−2pr2z2x3/2−2qr2z2y2.The intensity of LD beams can be investigated in two vertical planes. In the plane perpendicular to the junction (i.e.,y = 0), as shown in Figure 1, the substitution of y = 0 into equation (19) yields(20)Ix,0,z=1+k2x2+z2k2x2+z2E02exp−2px2+z2z2x3/2.Figure 1
Facet of a LD chip and the related coordinate system.In the plane parallel to the junction (i.e.,x = 0), as shown in Figure 1, the total intensity can be expressed as follows:(21)I0,y,z=1+k2y2+z2k2y2+z21+y6z6E02exp−2qy2+z2z2y2.


## 3. Experimental Procedure

The experiments were performed to examine the theoretical results using three high-power LDs (USHIO HL63391DC, TOSHIBA TOLD9441MC, and USHIO HL63290HD). The parameters are shown in Table1.Table 1
Parameter of LDs.
LDsWave modelOperating current (mA)Wavelength (nm)ParameterpParameterqUSHIO HL63391DCSingle transverse model2256391.55 × 1032.20 × 104TOSHIBA TOLD9441MCSingle transverse model506503353.05 × 104USHIO HL63290HDMultitransverse model2.46382561.55 × 104As shown in Figure2, the intensity profiles of laser beam were measured through a pinhole scan (radius is 100 μm) and a photodiode (LSGSPD-UL0.25, 0.25 mm visible light PIN photodiode, wavelength 500–880 nm, and 0.25 mm active diameter) behind the hole. The photodiode moved along the straight lines parallel to the output facet of the LDs’ chip in the x–z and y–z planes, where z = 50 mm. The uncertainty of measurements is less than 1%.Figure 2
Experimental setup.Figure3 shows the measurements and theoretical beam profiles of HL63391DC, and the intensity curve of the theory agrees with the experimental data in most portions. Figure 4 shows the light intensity profiles of TOLD9441MC. The calculated profiles agree well with the experimental data in most portions, except for the discrepancies in the low-intensity value regions in the x–z plane. The theoretical curve agrees well with the experimental data in the y–z plane. Figure 5 shows the light intensity profiles of HL63290HD, and the discrepancies of theory and measurement of this LD are greater than those of HL63391DC and TOLD9441MC because only single transverse mode exists in HL63391DC and TOLD9441MC, whereas multitransverse modes exist in HL63290HD. Thus, the output light field of the latter is more complicated, the shape intensity of two lobes appears in the y–z plane, and the theoretical curve does not fit the measurement.Figure 3
Measured and theoretical beam intensity profiles of HL63391DC. (a)x–z plane. (b) y–z plane.
(a)(b)Figure 4
Measured and theoretical beam intensity profiles of TOLD9441MC. (a)x–z plane. (b) y–z plane.
(a)(b)Figure 5
Measured and theoretical beam intensity profiles of HL63290HD. (a)x–z plane. (b) y–z plane.
(a)(b)Compared with the previous models of LD beam, including Hermite–Gaussian model [23], Gaussian model [25, 27, 28], elliptical Gaussian model [24, 29], and negative exponential Gaussian model [22, 30], the novel output model Eyx0,y0=E0exp−px03/2−qy02 in this article is more precise for single transverse model beam. For the calculation of the vector Rayleigh diffraction integrals, we expand expikr−r0 into a series by keeping the first, second, and third series expansions. The calculations make the diffraction integral of beam distribution with large divergence more reliable compared with the first two expansions in the article [24, 26].


## 4. Conclusion

A novel theoretical model for the nonparaxial vectorial field of high-power LDs was proposed, and the beam parameters were related to the structure of LDs’ waveguide. High-order approximations of the diffraction integral were calculated on the basis of the vector Rayleigh diffraction integrals, the fields parallel and perpendicular to beam propagation direction were considered, and the beam intensities of three high-power LDs beam were measured. The mathematical model provided a good fit to the experimental data of single transverse model beam of LDs. This mathematical model can be used to describe the beam propagation and shape of LDs in OWC and optical networks.



---

*Source: 7059638-2020-07-18.xml*