Next Article in Journal
Simultaneous and Multiplexed Measurement of Curvature and Strain Based on Optical Fiber Fabry-Perot Interferometric Sensors
Previous Article in Journal
Ultra-Broadband, Compact Arbitrary Ratio Power Splitters Enabled by Adiabatic Sub-Wavelength Grating
Previous Article in Special Issue
The Progress and Trend of Heterogeneous Integration Silicon/III-V Semiconductor Optical Amplifiers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Model of Dual Fabry–Perot Etalon-Based External-Cavity Tunable Laser Using Finite-Difference Traveling-Wave Method

1
State Key Laboratory of Integrated Optoelectronics, Institute of Semiconductors, Chinese Academy of Sciences, Beijing 100083, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Communication and Integrated Photonics Laboratory, Xiongan Institute of Innovation, Chinese Academy of Sciences, Xiong’an New Area 071899, China
*
Authors to whom correspondence should be addressed.
Member, IEEE.
Photonics 2023, 10(5), 579; https://doi.org/10.3390/photonics10050579
Submission received: 22 March 2023 / Revised: 22 April 2023 / Accepted: 9 May 2023 / Published: 16 May 2023
(This article belongs to the Special Issue Integrated Microwave Photonics)

Abstract

:
A physical model of an external-cavity tunable laser (ECTL) utilizing the vernier effect of a dual Fabry–Perot (FP) etalon is presented and simulated using the finite-difference traveling wave (FDTW) method. In this paper, we provide a detailed explanation of the physical principle and construction process of the model, as well as the simulation results for the laser. The model is precisely established by studying the time-dependent changes in the carrier concentration and optical field of different wavelengths inside the laser before reaching a steady state. By determining multiple parameters in the tuning region and gain region, the proposed model can calculate and predict various laser parameters, such as output power and side-mode suppression ratio (SMSR). Moreover, the FDTW method displays the change process of various parameters, such as carrier concentration and spectrum, in the convergence of various positions in the laser with femtosecond time resolution. This capability is promising for in-depth research on the inner mechanism of lasers.

1. Introduction

Semiconductor lasers have attracted significant attention since their invention due to their small size, long lifetime, narrow linewidth, high-frequency modulation, and wide tunable range. In various application scenarios, laser performance requirements can differ significantly. For instance, in high-tech fields, such as high-resolution spectroscopy and broadband communication network systems, narrow-linewidth tunable semiconductor lasers are essential. External-cavity lasers provide a mature approach to fabricate narrow-linewidth tunable lasers. Using external-cavity technology, semiconductor lasers can produce stable single-longitudinal-mode output and can also be tuned in the range of tens to hundreds of nanometers [1]. Furthermore, several other properties of semiconductor lasers have been improved, including higher output power and greater SMSR [2]. Due to these desirable features, ECTLs have found wide usage in white-light interference, wavelength-division multiplexing (WDM) systems, optical coherence tomography, coherent optical communications, gas detection, atomic physics experiments, atomic clocks, and other fields.
Currently, there are several mature schemes for ECTLs, including the thermally tunable etalon type [3], MEMS-mechanical tuning type [4], acousto-optic (AO) filter type [5], and silicon-based double microring type [6]. Although MEMS-tuned external cavities and AO filters can achieve good lasing characteristics, their output power is limited. Moreover, MEMS-based ECTLs have mechanical parts and complex configurations, while those using AO filters have a large cavity length of 50 mm. Although ECTLs based on microring resonators have been extensively adopted due to their compact size and superior performance, FP etalons remain relevant for specialized applications where their unique advantages are indispensable. FP etalons exhibit a broader free spectral range (FSR) and reduced insertion loss compared to microring resonators, rendering them apt for scenarios where these attributes are of critical importance. For instance, FP etalons facilitate a wider tuning range and increased output power, which is attributed to their diminished insertion loss. Moreover, they exhibit reduced susceptibility to fabrication imperfections and environmental influences, resulting in enhanced stability and reliability under particular circumstances. In contrast, ECTLs prepared by an etalon can change the output frequency by modulating the FSR, combined with the wavelength position of the transmission peak of the etalon, to achieve the wavelength selection function [7]. This design can achieve high SMSR, low relative intensity noise (RIN), and narrow linewidth [8,9], thereby meeting the application requirements of WDM systems. The channel accuracy of the ECTL designed by K. Sato in the C-band with a channel spacing of 50 GHz is less than 0.6 GHz [10]. In 2018, Giovanni B. de Farias proposed an ECTL that meets the C-band Integrated Tunable Laser Assembly (ITLA) specification, and 16-quadrature amplitude modulation (QAM) transmission was verified [9]. These results demonstrate the excellent properties of ECTLs, but unfortunately, they lack further theoretical analysis. Moreover, existing theories [11] are limited to steady-state solutions, making it challenging to demonstrate the process of steady-state formation and track the change process of various physical quantities of a laser before reaching a steady state. Therefore, it is difficult to analyze the reasons for various results using these theories.
J. Leuthold proposed a fitted piecewise cubic curve to describe stimulated emission gain in relation to various measurable physical parameters in the gain region of a laser [12]. This model has been extensively used in numerous publications and textbooks [13,14,15]. In this study, we integrate the aforementioned gain model with the FDTW method [16,17]. The gain region is partitioned into independent input and output systems with distinct carrier concentration and optical field distribution, while the tuning region is treated as a power phase modulation system that is time invariant. By considering the time differential perspective, we simulate and explicate the steady state of the double-etalon external-cavity tunable laser. Lastly, we present the simulations of the steady-state laser spectra obtained from a series of measurable physical parameters. Figure 1 displays the laser structure described in this article.
The structure of this paper is organized as follows: In Section 1, the transmission and reflection models for a tuning region that consists of two etalons are presented, along with a method for pre-selecting the FSR of the etalon. Section 3 describes the gain region model based on the work of J. Leuthold [12]. Section 4 introduces and discusses the simulation method used in this study. Section 5 presents the simulation results and a thorough discussion of the findings. Finally, in Section 6, the conclusions drawn from the research are presented.

2. Tuning Region

The tuning region comprises of two etalons and a reflector on the right, as depicted in Figure 1. The steady-state intensity transmission characteristics of the etalons can be explained using the current etalon filter power spectrum theory [18,19]. Since the etalons are a part of the laser resonator in our model, their impact on the optical path and phase of the cavity cannot be overlooked. The longitudinal-mode structure of the laser is determined by the longitudinal length of the resonator and the reflectivity of the reflector for the resonator formed by the etalons. Figure 2 shows the amplitude transmittance, amplitude reflectance, etalon internal refraction angle θ , etalon refractive index n , and etalon thickness d of each surface. Considering the angle θ between the etalons and the optical axis, the continuous-wave field components accumulate discretely with the inclined angle, thus suppressing parasitic oscillations within the resonant cavity while inevitably introducing tilt loss. Consequently, as the filtering components within the cavity, we do not need to focus on the reflective behavior in the gap of the etalons; instead, we should concentrate on their transmission properties. By neglecting absorption, the following relationship is satisfied i = 1 , 2 :
r i = r i r i 2 + t i t i = 1 R i = r i 2
The intensity reflectivity of the two surfaces is denoted as R i , and assuming λ is the wavelength of the input laser, the phase difference δ caused by a round trip of the beam in the etalon satisfies the following relationship:
δ = 2 n d cos θ λ 2 π
Let E 0 be the input light-field amplitude, E be the total transmitted light-field amplitude, E be the total reflected light-field amplitude, and the amplitude transmittance be t , we have the following relationship:
t = E E 0 = N = 1 + t 1 t 2 r 1 N 1 r 2 N 1 exp j N 1 δ = 1 R 1 1 R 2 1 R 1 R 2 exp j δ
The corresponding total intensity transmission coefficients are as follows:
T = t t = 1 R 1 1 R 2 1 R 1 R 2 2 + 4 R 1 R 2 sin 2 δ 2
When R 1 = R 2 = R 0 , let the fineness parameter F = 4 R 0 / 1 R 0 2 , and let the total phase difference Δ φ = t a n 1 R 0   s i n   δ / 1 R 0   c o s   δ , then
t = 1 R 0 1 R 0 exp j δ = 1 R 0 1 + R 0 2 2 cos δ exp j Δ φ
T = 1 1 + F sin 2 δ 2
The FSR and FWHM of the etalons are calculated as follows:
W F S R = c 2 n d cos θ   i n   H z 1 2 n d cos θ   i n   m 1 λ 2 2 n d cos θ   i n   m
W F W H M = 2 W F S R π F
The transmission spectrum T of the etalons can be completely determined by W F S R and W F W H M , as shown in Figure 3.
A high side-mode suppression ratio (SMSR) is essential for external-cavity tunable semiconductor lasers as it influences spectral purity, system performance, signal stability, and the minimization of crosstalk and interference. This crucial parameter directly impacts the applicability and effectiveness of such lasers across various fields, including optical communication, sensing, metrology, and spectroscopy. The SMSR of a laser is determined by the tuning region and the gain region, each based on distinct mechanisms.
For the tuning region, given that the gain spectrum is sufficiently flat, the transmission spectrum’s SMSR within the tuning region will be directly reflected in the final spectrum. The gain chip will non-linearly exacerbate this process, ultimately resulting in the laser’s SMSR being significantly larger than the SMSR of the transmission spectrum in the tuning region [20].
Equations (5)–(8), in conjunction with Figure 3, demonstrate that the filtering properties of an etalon, such as the SMSR of the transmitted spectrum, are dependent on its W F S R and W F W H M . In practical applications, it is often necessary to identify an appropriate combination of W F S R within a certain W F W H M while ensuring that the upper bound of the W F S R is not exceeded. This approach achieves the desired filtering performance when outputting specific wavelengths, such as 80 or 120 channels, within that band. The concept of finding such an FSR combination is as follows:
Expressing (4) with W F S R and W F W H M , we have
T = 1 1 + 4 W F S R 2 π 2 W F W H M 2 sin 2 ν π W F S R
which is simplified to
T Δ ν = 1 1 + 4 Δ ν 2 W F W H M 2
The SMSR of the filter is ( Δ ν 0 , Δ W F S R ) as follows:
m a x 1 1 + 4 Δ ν 2 W F W H M 1 2 × 1 1 + 4 Δ ν Δ W F S R 2 W F W H M 2 2
We traverse the possible FSR combinations and calculate the maximum value of their SMSR at the output of each band to select the most suitable FSR combination.

3. Gain Region

The gain region consists of the gain medium, which is electrically pumped to realize carrier concentration inversion, and the results in the stimulated emission optical gain for the incident light while also spontaneously emitting broad-spectrum fluorescence. The processes of spontaneous and stimulated emissions, as well as non-radiative transitions, together consume inversion carriers, and a constant current generates inversion carriers at a constant rate. The following analysis focuses on the changes in various related parameters as light passes through a thin layer in the gain region, as illustrated in Figure 4.
During the time interval d   t = d   x × n g / c , when photons travel from the left surface to the right surface of the thin layer, the reverse carrier concentration in the thin layer is pumped by the injection current in the gain region. The depletion of the inverse carrier concentration in the thin layer is caused by nonradiative recombination, spontaneous emission, Auger recombination, and stimulated emission. The changes in electron concentration and photon concentration during this time interval are given by [16]:
d N d t = η I i n j e c t q w d L A N B N 2 C N 3 + g λ , N N P λ c n g d λ
d   N P λ d   t = Γ c n g g λ , N N P λ + Γ β B λ N 2
The relationship between photon concentration and electric-field strength density is as follows:
N P λ = 2 ε 0 n e f f n g h ν E λ 2
For its rate of change in amplitude, we have the following:
d E s p o n λ = h ν 4 ε 0 n e f f n g Γ β B λ N 2 N P λ d t
d E s t i m λ = h ν 4 ε 0 n e f f n g Γ c n g g λ , N N P λ
After d x , the phase change is as follows:
d   φ = 2 π n g d   x λ
Thus, for the total light-field amplitude change, we have the following:
E o u t λ = e i φ 0 + d   φ E i n λ + d E s t i m λ + e i φ r d E s p o n λ
where φ 0 is the phase of E i n λ and φ r is a uniformly distributed random phase.
Next, we discuss the gain parameter g λ , N . We use a polynomial model for the material gain to describe g λ , N at a specific carrier concentration [12]:
g λ , N = c N λ λ z N 2 + d N λ λ z N 3 , λ < λ z N 0 , λ λ z N
where
c N = 3 g p N λ z N λ p N 2
d N = 2 g p N λ z N λ p N 3
g p N = a 0 N N 0 + a ¯ a 0 N 0   exp N N 0
λ p N = λ 0 [ b 0 N N 0 + b 1 N N 0 2 ]
λ z N = λ z 0 z 0 N N 0
Up to this point, we have determined two important parameters of the gain region: the gain coefficient g λ and the coefficient B for spontaneous radiative recombination.
To use the above model for laser simulations, several parameters of the tuning region and the gain region must be determined before the FDTW method can be applied [16,17]. See Table 1 for the details of the parameters.

4. Simulation

In this simulation, all the light-field behavior of the laser at a certain position and time is determined by the left-direction and right-direction complex-amplitude vectors at that position, which are composed of n + 1 complex numbers (where n = ( λ e n d λ s t a r t ) / Δ λ ), with each complex number representing the complex amplitude of the light field at a specific frequency.
In order to simulate the gain region, we partition it into a finite number of slices at equal intervals along the propagation direction of the light field, where the step size is denoted by Δ x , as illustrated in Figure 4. It is assumed that the injection current in the gain region remains constant until the steady state is reached. During the tiny time interval Δ t = n g Δ x / c , when the envelope of the optical field amplitude passes through the slice at a group velocity of c / n g , the carrier concentration remains constant. Any carrier concentration that is consumed or generated during this process will act on the slice after Δ t , thereby forming a gain sheet input–output system.
The gain region comprises n n = L / Δ L gain slice input and output systems connected in series. At a given time t = m Δ t ,   where   m = 0 , 1 , 2 , , the outputs of adjacent gain slice systems are connected to the inputs of the next system at the start of the next time interval Δ t . Specifically, for the leftmost gain sheet system, the left input light field is always the product of the left output light field and R l e f t , while for the rightmost gain sheet system, the input light field on the right is always the product of the output light field on the right and r t u n i n g λ . The propagation time of the light field in the tuning region is neglected, as this region is treated as passive and its physical quantities do not vary during our analysis. Moreover, the process of establishing the steady state inside the tuning region is not the focus of our study. The neglect of the propagation of the tuning region does not affect the overall steady state of the laser. The gain region, formed by the gain slices in series, is depicted in Figure 5.
At the outset, the n gain slices in the gain region exhibit no gain, and the reverse carrier concentration arises from the injection current generation in the gain region and the consumption of non-radiative recombination, spontaneous emission, and Auger recombination. This carrier concentration is determined by the steady state, which takes the following form:
η I i n j e c t q w d L A N 0 , i + B N 0 , i 2 + C N 0 , i 3 = 0
The inversion carrier concentration N 0 , i of the initial state can be determined, with the carrier concentration being the same for each sheet. The next step is to examine the input and output states of each gain slice in the initial state. As no gain exists in the initial state, the left and right outputs of each slice correspond to their spontaneous emission, with the following values:
E i n i t i a l λ = h ν w d Δ L N 0 , i 2 ε 0 n e f f n g e i φ r
In the initial state, the electric-field strength within the gain region is denoted by E i n i t i a l λ at each corresponding wavelength.
Next, we discuss the boundary conditions for iterative processes. Specifically, for the leftmost slice, the input on the left is always equal to the product of the output on the left and the reflectivity of the left mirror amplitude. This can be expressed as follows:
E 1 , t + Δ t , l , i n λ = R l e f t × E 1 , t , l , o u t λ
The initial numeral denotes the slice number from left to right. In the aforementioned equation, we are specifically referring to the leftmost slice, hence denoted by ‘1’. The second numeral corresponds to time. The third symbol signifies the position of the light field, where ‘l’ denotes the left surface of the sheet, and ‘r’ represents the right surface of the sheet. The final symbol signifies either incoming or outgoing light of the slice, where ‘in’ refers to incoming light and ‘out’ refers to outgoing light.
Regarding the rightmost slice, it can be observed that its input consistently equals the multiplication of the right output and the amplitude reflectance of the tuning region:
E 1 , t + Δ t , r , i n λ = r t λ × E 1 , t , r , o u t λ
The power of the light emitted at any time is calculated as follows:
P λ = 2 E n , t , r , o u t λ 2 ε 0 Δ t
Next, we examine the changes in the optical field and the carrier concentration for each slice per unit time Δ t during the iteration.
Within each slice, as each time interval Δ t elapses, we analyze the correlation between the output light field of all slices and the input light field in the corresponding direction using Equation (18), which can be expressed as follows:
E n , t + Δ t , l , o u t λ = f E ( E n , t , r , i n λ , N i , t )
E n , t + Δ t , r , o u t λ = f E ( E n , t , l , i n λ , N i , t )
Simultaneously, the evolution of the carrier concentration is governed by Equation (12), which can be formulated as follows:
N i , t + Δ t = f N ( N i , t , E n , t , r , i n λ , E n , t , l , i n λ )
For slices located between the leftmost and rightmost slices, during each time interval Δ t , the input light field of a given slice is determined by the output light field of its adjacent neighbor at the preceding moment. This relationship can be expressed as follows:
E n , t + Δ t , l , i n λ = E n 1 , t , r , o u t λ
E n , t + Δ t , r , i n λ = E n + 1 , t , l , o u t λ
Once the essential parameters, boundary conditions, and initial conditions have been established, the iterative process can commence. The iteration continues until the results converge or until a reasonable amount of time has passed, at which point all relevant physical parameters will have been determined.

5. Results

Compared to the simulation results obtained by other models, our proposed model provides a comprehensive visualization of the dynamic changes in various physical quantities over time, rather than just the final converged results. Additionally, our model simplifies the time-domain behavior of the two etalons in the tuning region, resulting in a faster convergence speed than the actual situation while preserving the spectral and energy changes without affecting the laser’s lasing spectrum. Using this model, we can explain and predict various physical phenomena and laws observed in the longitudinal-mode behavior of a laser, as demonstrated in Figure 6 and Figure 7, from a more microscopic perspective.
In the absence of any specific instructions regarding parameter settings in this paper, the following default values are adopted: injection current I = 0.1   A , reflectivity of the mirror on the right side of the tuning area R = 0.5 , and output wavelength λ = 1590   nm . Please refer to Table 1 for the remaining parameter values.
In Figure 6, the wavelengths studied were selected based on the maximum points of the transmission spectrum in the tuning region, thus allowing us to examine various physical quantities with minimal computational effort. As shown in Figure 6a, the side modes briefly disappear before the main mode narrows. Figure 6b indicates that the power converges more quickly than the SMSR due to the faster convergence of the average carrier concentration, while the convergence process of the carrier concentration distribution is slower, leading to the slower convergence of the secondary mode. Figure 6c shows that the selection process of the main mode is rapid, but the spectral line shaping takes longer.
The reverse carrier concentration exhibits an opposite trend to the output power. In Figure 7, we see a rapid drop in the carrier concentration, followed by a slow rise and eventual convergence to a steady state. This is due to the consumption of inversion carriers during lasing and a gain clamping effect. The gradient of the spatial distribution of carriers is relatively small compared to the change in carrier concentration over time, and carrier concentration near the output side is slightly lower due to the stronger optical field and greater carrier consumption.
Figure 8 demonstrates that under the current parameters, the laser has a flat gain spectrum from 1583 nm to 1596 nm, which is important for designing tunable lasers to achieve specific functions. The gain spectrum at different wavelengths can be calculated using our model under different parameters.
Figure 9a shows the convergence process of the output power under different injection currents, which is similar for different currents with rapid initial growth followed by slow decline and convergence. The convergence is determined by the gain spectrum at different carrier concentrations. Figure 9b indicates that the convergence results are linear for currents between 0.01 A and 0.1 A. The convergence process for different output wavelengths varies depending on its position on the gain spectrum and the injection current.
Figure 10a shows the convergence process of the carrier concentration at different injection currents, while Figure 10b indicates that the converged system gains are almost the same under different currents. In the steady state, a larger current produces more photons in the cavity, thus consuming more carriers in the gain region, which requires a larger current to supplement the consumed carriers.
Figure 11 illustrates the convergence of the laser’s SMSR over time at various currents. In comparison to Figure 3, the laser output’s SMSR is two to three orders of magnitude higher than the SMSR of the transmission spectrum. This can be attributed to the gain region discussion, where the gain chip exponentially amplifies the light at each wavelength, consequently increasing the SMSR by several orders of magnitude. With all other conditions remaining constant, raising the current will solely enhance the power of the primary mode. Therefore, a linear relationship exists between the converged SMSR and current, as evident in Figure 11b.
In summary, this paper describes a model for simulating the behavior of a tunable laser based on a distributed feedback structure. The model takes into account the interaction between the optical field and carrier concentration in the active region of the laser, as well as the effects of various physical phenomena, such as gain clamping and mode competition. The simulation results show that the model can predict the behavior of the laser with high accuracy and provide insights into the underlying physical mechanisms involved. The model can also be used to optimize the design of tunable lasers for specific applications by calculating the gain spectra under different parameters.

6. Conclusions

This paper presents a comprehensive analysis of the operation and performance of a dual FP etalon-based ECTL. The paper provides insights into the effect of the tuning region on the phase modulation of light and proposes an idea for designing etalon FSR. Additionally, the paper uses the FDTW method to discuss the changes in the differential gain index of the gain region with the carrier concentration and analyzes the influence of laser injection current on various parameters, such as power, carrier concentration distribution, and SMSR. The theoretical model presented in this paper can provide a valuable reference for the design and optimization of ECTLs.

Funding

This research was funded by the National Key Research and Development Program of China (No. 2022YFB2804501).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kobayashi, K.; Mito, I. Single frequency and tunable laser diodes. J. Light. Technol. 1988, 6, 1623–1633. [Google Scholar] [CrossRef]
  2. O’carroll, J.; Phelan, R.; Kelly, B.; Byrne, D.; Smyth, F.; Cardiff, B.; Anandarajah, P.M.; Barry, L.P. Narrow linewidth discrete mode laser diodes for advanced modulation formats. In OFC/NFOEC; Optics: Washington, DC, USA, 2012; pp. 1–3. [Google Scholar]
  3. Finot, M.; McDonald, M.; Bettman, B.; Sell, J.; Daiber, A.; Chapman, W.B.; Kozlovsky, W.J. Thermally tuned external cavity laser with micromachined silicon etalons: Design, process and reliability. In Proceedings of the 54th Electronic Components and Technology Conference (IEEE Cat. No.04CH37546), Las Vegas, NV, USA, 1–4 June 2004; Volume 1, pp. 818–823. [Google Scholar]
  4. Numata, K.; Camp, J.; Krainak, M.A.; Stolpner, L. Performance of planar-waveguide external cavity laser for precision measurements. Opt. Express 2010, 18, 22781–22788. [Google Scholar] [CrossRef] [PubMed]
  5. Takabayashi, K.; Takada, K.; Hashimoto, N.; Doi, M.; Tomabechi, S.; Nakazawa, T.; Morito, K. Widely (132 nm) wavelength tunable laser using a semiconductor optical amplifier and an acousto-optic tunable filter. Electron. Lett. 2004, 40, 1187–1188. [Google Scholar] [CrossRef]
  6. Segawa, T.; Matsuo, S.; Kakitsuka, T.; Sato, T.; Kondo, Y.; Suzuki, H. Full C-Band Tuning Operation of Semiconductor Double-Ring Resonator-Coupled Laser with Low Tuning Current. IEEE Photonics Technol. Lett. 2007, 19, 1322–1324. [Google Scholar] [CrossRef]
  7. Xue, J.; Pei, X.; Su, B.; Sun, L.; Zhao, H. Study of F-P Etalon Using in Lasers. Laser Optoelectron. Prog. 2012, 49, 104–107. [Google Scholar]
  8. De Merlier, J.; Mizutani, K.; Sudo, S.; Sato, K.; Kudo, K. Wavelength Channel Accuracy of an External Cavity Wavelength Tunable Laser with Intracavity Wavelength Reference Etalon. J. Light. Technol. 2006, 24, 3202–3209. [Google Scholar] [CrossRef]
  9. De Farias, G.B.; de Andrade, H.A.; Zanvettor, L.T.; de Carvalho, L.H.; Ubaldi, M.C.; Righetti, A.; Meli, F.; Ferrario, L.; Grasso, G. Narrow linewidth (<100 kHz) and high output power (>16 dBm) external cavity lasers for coherent optical communications. In Proceedings of the 2018 SBFoton International Optics and Photonics Conference (SBFoton IOPC), Campinas, Brazil, 8–10 October 2018; pp. 1–5. [Google Scholar]
  10. Sato, K.; Mizutani, K.; Sudo, S.; Tsuruoka, K.; Naniwae, K.; Kudo, K. Wideband External Cavity Wavelength-Tunable Laser Utilizing a Liquid-Crystal-Based Mirror and an Intracavity Etalon. J. Light. Technol. 2007, 25, 2226–2232. [Google Scholar] [CrossRef]
  11. Xu, C.; Du, X.; Jin, Y.; Chen, W.; Zhu, N. Analytical Models for External Cavity Lasers Using Tunable Etalons. IEEE Photonics J. 2022, 14, 1–8. [Google Scholar]
  12. Leuthold, J.; Mayer, M.; Eckner, J.; Guekos, G.; Melchior, H.; Zellweger, C. Material gain of bulk 1.55 μm InGaAsP/InP semiconductor optical amplifiers approximated by a polynomial model. J. Appl. Phys. 2000, 87, 618–620. [Google Scholar] [CrossRef]
  13. Schares, L.; Schubert, C.; Schmidt, C.; Weber, H.; Occhi, L.; Guekos, G. Phase dynamics of semiconductor optical amplifiers at 10–40 GHz. IEEE J. Quantum Electron. 2003, 39, 1394–1408. [Google Scholar] [CrossRef]
  14. Das, N.K.; Yamayoshi, Y.; Kawaguchi, H. Analysis of basic four-wave mixing characteristics in a semiconductor optical amplifier by the finite-difference beam propagation method. IEEE J. Quantum Electron. 2000, 36, 1184–1192. [Google Scholar]
  15. Wang, J.; Maitra, A.; Poulton, C.G.; Freude, W.; Leuthold, J. Temporal Dynamics of the Alpha Factor in Semiconductor Optical Amplifiers. J. Light. Technol. 2007, 25, 891–900. [Google Scholar] [CrossRef]
  16. Khan, F.K.; Cassidy, D.T. Analysis of broadly tunable coupled-cavity semiconductor lasers. Appl. Opt. 2011, 50, 4834–4843. [Google Scholar] [CrossRef] [PubMed]
  17. Rossetti, M.; Bardella, P.; Montrosset, I. Modeling Passive Mode-Locking in Quantum Dot Lasers: A Comparison between a Finite-Difference Traveling-Wave Model and a Delayed Differential Equation Approach. IEEE J. Quantum Electron. 2011, 47, 569–576. [Google Scholar] [CrossRef]
  18. Gunning, W. Double-cavity electrooptic Fabry-Perot tunable filter. Appl. Opt. 1982, 21, 3129–3131. [Google Scholar] [CrossRef] [PubMed]
  19. Taylor, H. Design of multireflector resonant bandpass filters for guided wave optics. J. Light. Technol. 2001, 19, 866–871. [Google Scholar] [CrossRef]
  20. Chen, Q.; Lu, Q.; Guo, W. Theory and simulation of multi-channel interference (MCI) widely tunable lasers. Opt. Express 2015, 23, 18040–18051. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Basic structure of the mentioned laser. The inner ring of the etalons is a thermistor for temperature monitoring, and the outer ring is a metal electrode for temperature adjustment.
Figure 1. Basic structure of the mentioned laser. The inner ring of the etalons is a thermistor for temperature monitoring, and the outer ring is a metal electrode for temperature adjustment.
Photonics 10 00579 g001
Figure 2. Laser propagation behavior in the etalon.
Figure 2. Laser propagation behavior in the etalon.
Photonics 10 00579 g002
Figure 3. The FSRs and FWHMs of the two etalons shown in the transmission spectrum.
Figure 3. The FSRs and FWHMs of the two etalons shown in the transmission spectrum.
Photonics 10 00579 g003
Figure 4. Analysis of gain region.
Figure 4. Analysis of gain region.
Photonics 10 00579 g004
Figure 5. The gain region changes state through Δ t .
Figure 5. The gain region changes state through Δ t .
Photonics 10 00579 g005
Figure 6. The output power of each wavelength varies with time. (a) Thermal diagram of the output power at each wavelength as a function of time. (b) Variation plot of dominant wavelength output power and SMSR with time. (c) Power spectra at five instants. The different colored points in (b) correspond to the curves of the same color in (c).
Figure 6. The output power of each wavelength varies with time. (a) Thermal diagram of the output power at each wavelength as a function of time. (b) Variation plot of dominant wavelength output power and SMSR with time. (c) Power spectra at five instants. The different colored points in (b) correspond to the curves of the same color in (c).
Photonics 10 00579 g006
Figure 7. Variation in the carrier concentration with time at various positions in the gain region. (a) Heat map of the change in carrier concentration at each location. (b) Average carrier concentration in the gain region as a function of time. (c) Carrier concentration distributions and ranges at five instants. The different colored points in (b) correspond to the curves of the same color in (c).
Figure 7. Variation in the carrier concentration with time at various positions in the gain region. (a) Heat map of the change in carrier concentration at each location. (b) Average carrier concentration in the gain region as a function of time. (c) Carrier concentration distributions and ranges at five instants. The different colored points in (b) correspond to the curves of the same color in (c).
Photonics 10 00579 g007
Figure 8. Output power at different output wavelengths. (a) The output spectra at every 0.1 nm from 1580 nm to 1602 nm, with different colors indicating the output spectra at different output wavelengths. (b) The envelope composed of the main peaks of the output spectra of different colors, showing a relatively flat output power from 1583 nm to 1596 nm.
Figure 8. Output power at different output wavelengths. (a) The output spectra at every 0.1 nm from 1580 nm to 1602 nm, with different colors indicating the output spectra at different output wavelengths. (b) The envelope composed of the main peaks of the output spectra of different colors, showing a relatively flat output power from 1583 nm to 1596 nm.
Photonics 10 00579 g008
Figure 9. Power behavior changes with current. (a) Convergence process of power under different injection currents. (b) Convergence results and fitting curves of power under different currents. The different colored points in (b) correspond to the curves of the same color in (a).
Figure 9. Power behavior changes with current. (a) Convergence process of power under different injection currents. (b) Convergence results and fitting curves of power under different currents. The different colored points in (b) correspond to the curves of the same color in (a).
Photonics 10 00579 g009
Figure 10. Carrier concentration changes with current. (a) Convergence process of carrier concentration under different injection currents. (b) Initial carrier concentration and converged carrier concentration under different currents. The different colored points in (b) correspond to the curves of the same color in (a).
Figure 10. Carrier concentration changes with current. (a) Convergence process of carrier concentration under different injection currents. (b) Initial carrier concentration and converged carrier concentration under different currents. The different colored points in (b) correspond to the curves of the same color in (a).
Photonics 10 00579 g010
Figure 11. Behavior of the SMSR at different currents. (a) Changes in SMSR over time. (b) Variation in the results with current after the SMSR converges. The different colored points in (b) correspond to the curves of the same color in (a).
Figure 11. Behavior of the SMSR at different currents. (a) Changes in SMSR over time. (b) Variation in the results with current after the SMSR converges. The different colored points in (b) correspond to the curves of the same color in (a).
Photonics 10 00579 g011
Table 1. Parameters appearing in the model.
Table 1. Parameters appearing in the model.
SymbolValueParameter
q 1.602 × 10 19   C Electron charge
c 2.998 × 10 8   m / s The speed of light in vacuum
h 6.625 × 10 34   J · s Planck’s constant
ε 0 8.854 × 10 12   F / m Vacuum dielectric constant
Δ λ 1   nm Wavelength interval of our work
w 2   μ m The width of the gain region
d 0.2   μ m The thickness of the gain region
L 1   mm The length of the gain region
Δ L 1   μ m The length of the gain slice
A 1.25 × 10 8 Unimolecular non-radiative recombination coefficient
B 0.8 × 10 16   m 3 / s Bimolecular spontaneous radiative recombination coefficient
C 0.9 × 10 42   m 6 / s Auger recombination coefficient
η 0.7 Current injection efficiency
W F S R 1 424   GHz The FSR of the first etalon at room temperature
W F S R 2 450   GHz The FSR of the second etalon at room temperature
W F W H M 1 10   GHz The FWHM of the first etalon at room temperature
W F W H M 2 10   GHz The FWHM of the second etalon at room temperature
N 0 6.5 × 10 23   m 3 Transparency carrier density
a 0 3.13 × 10 20   m 2 Differential gain factor
λ 0 1.575 × 10 6   μ m Band-edge wavelength
b 0 3.17 × 10 32   m 3 μ m Peak wavelength shift coefficient
λ z 0 1.625 × 10 6   μ m Begin-to-zero gain region
z 0 2.5 × 10 33   m 3 μ m Zero-gain wavelength shift coefficient
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Du, X.; Xu, C.; Jin, Y.; Zhai, K.; Wang, W.; Chen, W.; Zhu, N. A Model of Dual Fabry–Perot Etalon-Based External-Cavity Tunable Laser Using Finite-Difference Traveling-Wave Method. Photonics 2023, 10, 579. https://doi.org/10.3390/photonics10050579

AMA Style

Du X, Xu C, Jin Y, Zhai K, Wang W, Chen W, Zhu N. A Model of Dual Fabry–Perot Etalon-Based External-Cavity Tunable Laser Using Finite-Difference Traveling-Wave Method. Photonics. 2023; 10(5):579. https://doi.org/10.3390/photonics10050579

Chicago/Turabian Style

Du, Xinhao, Changda Xu, Ya Jin, Kunpeng Zhai, Wenting Wang, Wei Chen, and Ninghua Zhu. 2023. "A Model of Dual Fabry–Perot Etalon-Based External-Cavity Tunable Laser Using Finite-Difference Traveling-Wave Method" Photonics 10, no. 5: 579. https://doi.org/10.3390/photonics10050579

APA Style

Du, X., Xu, C., Jin, Y., Zhai, K., Wang, W., Chen, W., & Zhu, N. (2023). A Model of Dual Fabry–Perot Etalon-Based External-Cavity Tunable Laser Using Finite-Difference Traveling-Wave Method. Photonics, 10(5), 579. https://doi.org/10.3390/photonics10050579

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop