Sub-megahertz linewidth 780 . 24 nm distributed feedback laser for 87 Rb applications

A distributed feedback GaAs-based semiconductor laser with a laterally coupled grating is demonstrated at a wavelength of 780.24 nm with up to 60 mW power. A mode expander and aluminum-free active layers have been used to reduce the linewidth to 612 kHz while maintaining high output power. The laser demonstrates over 40 dB side-mode suppression ratio with >0.3 nm of tuning suitable for atom cooling experiments with the D2 87Rb atomic transition. This laser has substantial potential to be integrated into miniaturized cold atom systems.

Emerging technologies based on cold atoms such as atomic clocks [1,2] rely on accurate measurements through locking a laser to a stable atomic transition. These high-accuracy devices, however, require stringent laser properties in terms of linewidth, stability, and power output; therefore, the laser systems are typically bulky, expensive, and power-hungry [1,2]. Semiconductor lasers offer good chip-scale integrability [3][4][5], low production cost, and reduced power consumption [6], which are all desirable properties for the realization of compact clock systems [7]. Of particular relevance is the GaAs/AlGaAs material platform because it can cover a wide wavelength emission range from 700 nm to 900 nm, which includes important cold atom transition lines such as 87 Rb and 133 Cs at 780.24 nm and 852 nm, respectively [8].
One of the critical issues for the development of high-power lasers on a GaAs/AlGaAs material platform is the aluminumoxide formation on the air-exposed surfaces. Not only does such oxide formation decrease the power density threshold for catastrophic optical damage (COD) and mirror damage (COMD) [9,10], but it also degrades the laser output in terms of carrier lifetime, propagation losses, and self-heating [11][12][13]. Moreover, the aluminum fraction in the epilayer material poses several issues for epitaxially regrown structures where the oxygen contamination of the aluminum-containing layers can substantially decrease the fabrication yield. One of the possible solutions to mitigate the aluminum oxidation is the use of an aluminum-free active area [14]. From an optical gain point of view, the majority of the structures reported in the literature are designed with a large overlap between the optical mode and the gain region to minimize the laser threshold. Such design strategy is, however, not ideal for high-power and narrow-linewidth lasers as structures with low power density and low internal losses are desirable to both extend the power range with singlemode operation and reduce the linewidth in accordance to the Schawlow-Townes formula [15].
In this work, a distributed feedback (DFB) laser design is proposed to reduce the linewidth at the emission wavelength of 780.24 nm. An aluminum-free active area was designed with InGaAsP quaternary layers to reduce the risk of COD and COMD. In addition, a mode expander layer was introduced in the n-doped AlGaAs cladding to stretch the optical mode vertically, thereby decreasing the average power density and the internal propagation losses. Finally, a long cavity of 4 mm was employed to further decrease the average power density and reduce the laser linewidth. The fabricated DFB lasers exhibited stable emission wavelength of 780.24 nm with powers up to 60 mW and linewidths down to 612 kHz, a value suitable for the 87 Rb atom cooling at the D2 transition (<1 MHz [16]).
The first aspect we considered in the laser optimization was the epilayer design of the active region. In order to limit the interaction of the optical mode with the aluminum-containing layers, InGaAsP quaternary compounds were chosen for the quantum wells (QWs) and barriers, and these layers were simulated using the Schrödinger-Poisson current solver Nextnano++ [17]. Two compressively strained 4-nm-thick In 0.413 Ga 0.587 As 0.44 P 0.56 QWs inserted between 10-nm-thick In 0.27 Ga 0.73 As 0.44 P 0.56 barriers provide laser emission around 780.24 nm. The active area is enclosed between two graded separate confinement heterostructure (GRINSCH) layers of AlGaAs with aluminum fraction gradually increasing from 0.4 to 0.55. Finally, the GRINSCH layers are sandwiched between two doped Al 0.55 Ga 0.45 As cladding layers. The use of an aluminum-free active area reduces the optical confinement factor of the optical mode within aluminum-containing layers by 10% in comparison to standard GaAs/AlGaAs designs [18] and gives more tolerance in defining the DFB gratings as the active area can be etched without introducing additional scattering losses or decreasing the power density threshold for COD and COMD [9,10]. The Schawlow-Townes formula is an effective tool to analyze the impact of the epilayer design on the laser linewidth [15]. For semiconductor DFB lasers, the formula can be expressed as Equation 1 indicates that lasers with low ν require low propagation losses α i , equivalent mirror losses for DFB lasers α m , mode confinement in the active QW region QW , and Henry factor α H . Also, linewidth narrowing can be obtained by increasing the output power P out of the laser. It is worth noting that α i and QW are predominantly related to the design of the gain material, while α m depends on the grating design; hence, these contributions to the linewidth can be independently optimized. The linewidth enhancement factor α H has a substantial impact on nonlinear effects such as spatial hole burning (SHB), and is critical for narrow-linewidth semiconductor lasers as it increases for high values of P out [19], which explains the linewidth broadening experienced by high-power lasers. The other terms (i.e., hν photon energy and r sp /R sp spontaneous emission fraction in the lasing mode) effecting the linewidth in Eq. (1) can be assumed as fixed terms in this treatment and will not be considered in the following.
As reported in previous work for 1550 nm lasers [20], the introduction of a mode expander layer in the n-doped cladding decreases several parameters that play a role in defining the laser linewidth in Eq. (1): (i) the mode confinement QW ; (ii) the propagation losses α i ; and (iii) the average power density in the transversal direction P tran . The mode expander was implemented by introducing a layer in the lower cladding with a higher refractive index (i.e., Al 0.35 Ga 0.65 As) in comparison with the surrounding layer (i.e., Al 0.55 Ga 0.45 As), so as to pull the optical mode further into the n-doped cladding and expand the mode size. The magnitude of this effect depends on the mode expander thickness and on its distance from the active area; hence, the laser parameters QW and α i can be engineered as a function of the mode expander design. The vertical mode profiles with and without the mode expander, the band structure, and the refractive index profile around the active area are reported in Fig. 1. The distance of the mode expander from the active area is a critical parameter because the lower energy bandgap of the mode expander layer can create an asymmetrical confinement for electrons and holes in the active area. From simulation of the energy levels and carrier confinements in the active area by Nextnano++ software, a distance between the active and the mode expander larger than 100 nm was found to be sufficient to ensure a symmetrical band structure in the active area. On the other hand, a larger distance would reduce the effect of the mode expander on the active area, which would demand a wider mode expander to achieve an equivalent mode attraction. The lower energy bandgap of the mode expander, however, acts as a potential barrier for electrons and holes; hence, it is preferable to limit the mode expander thickness to minimize the value of an additional series resistance through the epilayer. Based on these considerations, a mode expander distance of 100 nm was chosen for the epilayer design.
The fraction of the modal power overlapping with the active region QW and the propagation losses α i , given by the sum of the propagation losses in the active region α QW = QW k p , n-doped α n = n k n , and p-doped cladding α p = p k p [20], are simulated through a finite difference eigenmode (FDE) 1Dsimulation tool (Lumerical Mode solution [21]). The factors k QW,n, p represent the normalized absorption coefficients of the QWs and claddings for doping levels n= 10 18 cm −3 and p = 6 × 10 17 cm −3 , respectively, and are obtained considering the free carrier absorption and the intervalence band absorption as the main loss mechanisms [22][23][24]. Both QW and α i are reported in Fig. 2 as a function of the mode expander thickness.
From Eq. (1), a lower value of QW provides a narrower laser linewidth; however, the reduction of QW also effects other laser parameters such as the threshold current and the output power P out [25]. Therefore, a trade-off value of 1% was selected as a design parameter as this QW value is sufficient to achieve a threshold current density as low as 1 kA/cm −2 [26] while it narrows the linewidth by a factor 2.4 in comparison to an epilayer design without the mode expander layer. For a mode expander thickness of 250 nm, which corresponds to a QW value of 1%, Fig. 2 also indicates a decrease of the propagation losses from α i from 5.4 cm −1 to values below 4.1 cm −1 . This is a consequence of the larger mobility of electrons [23] compared to holes [24] in AlGaAs compounds, which translates into lower absorption in n-doped materials compared to p-doped materials. The mode expander design also has a good tolerance in terms of thickness and material composition. An inaccuracy of 10 nm in layer thickness or 1% in aluminum fraction, which is significant compared to the tolerances for QW growths, translates in a variation of α i and QW lower than 0.1 cm −1 and 0.1%, respectively. The decrease of α i entails a further reduction of the linewidth according to Eq. (1). Moreover, the mode expander decreases the average power density in the transversal direction P tran , which is linked to the power output by the equation where L is the cavity length and A is the modal cross-sectional area (i.e., defined as the area contained between the points where the field intensity decreases by a factor 1/e over the maximum). The area A depends on the transversal geometry of the waveguide, but the vertical and horizontal contributions can be considered separately as the former is mainly effected by the material epilayer and the latter is shaped by the ridge waveguide geometry. The vertical mode size d , which represents the projection of A in the vertical line including the maximum of the power intensity, is increased from 0.4 µm to 0.48 µm by the mode expander layer. As a consequence, the mode size expansion reduces the average transversal power density P tran by about 20% for the same output power.
Regarding the Bragg grating design, the two key parameters to be optimized are the coupling factor κ and the cavity length L. The product κ L is the total grating coupling and determines the modal gain threshold α m in Eq. (1). It should be noted that although the grating coupling depends only on the product κ L, gratings with the same values of κ L but with different relative values of κ and L have different impacts on the laser properties. In fact, for the same product κ L, a longer cavity translates into a lower power density and, hence, a reduction of the nonlinear effects, which could broaden the linewidth. A longer cavity, however, has higher propagation losses per round trip α i L, which increases the lasing threshold and decreases P out ; hence, its value has to be carefully optimized by considering all the laser design parameters simultaneously.
The coupling factor κ was simulated through a FDE 2Dsimulator with a longitudinal approximation and confirmed through a 3D-simulator EigenMode Expansion (EME) by Lumerical MODE software for a shallow-etched sidewall Bragg grating, with a width of 1.2 µm and 2.2 µm for the two sections of the ridge waveguide, as described in a previous work [27]. The considered etching depth for the sidewall grating is 1.6 µm, which corresponds to stopping the etch 150 nm above the active area. Although the coupling coefficient κ is reduced by the effect of the mode expander layer, which pulls the optical mode out of the grating etched in the p-doped cladding, a lower value of κ is functional for achieving narrow linewidth, provided that the total grating coupling κ L is above 1 [28]. In order to have the Bragg condition at the target wavelength of 780.24 nm, the grating period at the first order is ∼117 nm , which gives very low tolerances to fabrication inaccuracies. Hence, a third-order grating geometry with a 50% duty cycle and a 351 nm period was selected to relax the fabrication tolerances, providing a coupling factor of approximately 4 cm −1 [27], which is consistent with the coupling coefficient measured from the stop band in the DFB spectrum. The cavity length was chosen to be 4 mm, which gives a κ L product of 1.6, a value optimal to ensure stable DFB lasing operation under single-mode operation [28]. Although a longer cavity increases the threshold current and the total propagation losses, such an increase is less prominent in lasers with a small number of QWs and low mode confinement QW , as for the epilayer design in this work [25]. Longer cavities are also expected to improve the thermal behavior of the laser as the dissipated heat can be spread over a longer length.
The epitaxial layers were grown by a commercial supplier on 650-µm-thick GaAs substrates. The sidewall DFB grating was patterned on a 500-nm-thick hydrogen silsesquioxane resist by electron beam lithography using a Vistec VB6 tool. The pattern was transferred to the sample by reactive ion etching using an Oxford Instrument ICP180 tool with a Cl 2 /BCl 3 /Ar/N 2 chemistry to etch through the upper cladding with a high aspect ratio, as shown in Fig. 3. The sample was then covered with a SiO 2 passivation layer, and a narrow window on the top of the waveguide was etched to allow for current injection. The top (TiPtAu) and bottom (AuGeNi) contacts were evaporated after the sample was thinned down to a thickness of approximately 250 µm. Finally, the device was cleaved before antireflection (AR) and high-reflection (HR) coatings were deposited on the cavity facets to suppress any parasitic Fabry-Perot effects.
The fabricated devices were fully characterized for lightvoltage-current (LIV), optical spectrum, and linewidth characteristics under continuum-wave operation and at a fixed temperature of 20 • C. The LIV characteristics shown in Fig. 4 indicate a threshold current of 140 mA with a voltage bias lower than 2 V up to an injection current of 250 mA, which demonstrates a good electrical conductivity of the epilayer. The DFB laser has kink-free emission between the threshold and an injection current of 240 mA, albeit the power output linearly increases until the injection current of 275 mA exceeds 60 mW. The laser spectrum was measured by an optical spectrum analyzer with 0.01 nm resolution for increasing values of the injected current in order to confirm single-mode operation and emission wavelength tuning (see Fig. 5). The DFB laser clearly exhibits single-mode behavior with no evidence of mode hopping and a side-mode suppression ratio (SMSR) exceeding 40 dB within the 180-250 mA current range and with a peak wavelength tunability of 0.3 nm. Finally, the laser linewidth was characterized through a heterodyne detection technique [29]. The laser signal was beaten with a commercial 50 kHz Ti:sapphire laser source at 780 nm wavelength and coupled to a fast photodetector to generate a radio frequency (RF) signal. The generated RF beat note was then resolved by a RF spectrum analyzer in order to measure the full width at half-maximum (FWHM) of the convolved signals. A Voigt line shape, the convolution of a Gaussian and a Lorentzian line shape, was employed to fit the RF beat note. The natural Lorentzian-shaped laser linewidth was extrapolated from the fitted Voigt linewidth. The measured RF signal and its Voigt-shaped fit are shown in Fig. 6 with a measurement sweep time of 25 µs. The RF fitting indicates that the laser has a linewidth of 612 kHz, a value that is suitable for 87 Rb atomic cooling applications [16].
In conclusion, an optimized, aluminum-free active region DFB laser with a mode expander and a long-cavity laterally coupled Bragg grating has been demonstrated at 780.24 nm wavelength with a 60 mW emission. The measured SMSR of 40 dB and linewidth of 612 kHz match the requirements for 87 Rb D2 transitions for atom cooling applications.
Disclosures. The authors declare no conflicts of interest.