Next Article in Journal
Dual-Branch Feature Fusion Network for Salient Object Detection
Next Article in Special Issue
Nonlinear Dynamics of a Single-Mode Semiconductor Laser with Long Delayed Optical Feedback: A Modern Experimental Characterization Approach
Previous Article in Journal
High-Flexibility Control of Structured Light with Combined Adaptive Optical Systems
Previous Article in Special Issue
Effects of Asymmetric Coupling Strength on Nonlinear Dynamics of Two Mutually Long-Delay-Coupled Semiconductor Lasers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Semiconductor Lasers under External Optical Feedback from Both Sides of the Laser Cavity

by
Mónica Far Brusatori
* and
Nicolas Volet
Department of Electrical and Computer Engineering, Aarhus University, 8200 Aarhus, Denmark
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(1), 43; https://doi.org/10.3390/photonics9010043
Submission received: 30 November 2021 / Revised: 7 January 2022 / Accepted: 12 January 2022 / Published: 14 January 2022
(This article belongs to the Special Issue Nonlinear Dynamics of Semiconductor Lasers and Their Applications)

Abstract

:
To increase the spectral efficiency of coherent communication systems, lasers with ever-narrower linewidths are required as they enable higher-order modulation formats with lower bit-error rates. In particular, semiconductor lasers are a key component due to their compactness, low power consumption, and potential for mass production. In field-testing scenarios their output is coupled to a fiber, making them susceptible to external optical feedback (EOF). This has a detrimental effect on its stability, thus it is traditionally countered by employing, for example, optical isolators and angled output waveguides. In this work, EOF is explored in a novel way with the aim to reduce and stabilize the laser linewidth. EOF has been traditionally studied in the case where it is applied to only one side of the laser cavity. In contrast, this work gives a generalization to the case of feedback on both sides. It is implemented using photonic components available via generic foundry platforms, thus creating a path towards devices with high technology-readiness level. Numerical results shows an improvement in performance of the double-feedback case with respect to the single-feedback case. In particularly, by appropriately selecting the phase of the feedback from both sides, a broad stability regime is discovered. This work paves the way towards low-cost, integrated and stable narrow-linewidth integrated lasers.

1. Introduction

The effect of external optical feedback (EOF) on diode laser dynamics has been extensively studied for the past half century [1,2]. EOF has been proven to affect laser performance, showing regimes that can aid in linewidth reduction [3,4,5,6], as well as others responsible for highly unstable behavior, from mode hopping to the case of coherence collapse [7,8,9,10,11,12]. Methods to improve laser stability thus need to take EOF into account, as even weak feedback can be detrimental. A traditional approach to mitigate its effects is to include an off-chip isolator at the laser output. Yet, this component negatively impacts the dimensions of packaged devices as well as fabrication times and costs. As such, research is ongoing to develop an integrated solution that can minimize the negative effects of EOF. Efforts include adjusting the feedback phase to tune into line-narrowing regimes [13,14], using unidirectional phase modulators [15,16], reducing the linewidth enhancement factor, e.g., using quantum dots [17,18,19], employing electromagnetic effects [20,21], harnessing the mode propagation properties of ring lasers [22,23], or the extended cavity approach [24,25,26,27].
Double external feedback on the same side of the laser cavity has been previously explored for chaos stabilization in Reference [28]. Its results numerically show that, by introducing an additional feedback term, a chaotic regime can be stabilized in terms of output power, where a robust stable region for a wide parameter range can be seen. Furthermore, linewidth is numerically shown to be reduced with respect to the single feedback case. However, it lacks an analytical expression for the linewidth, and it does not explore the effect of the added term on other optical feedback regimes, which limits the scope of the analysis [7]. Dual external cavities have also been explored to maximize the sensitivity of self-mixing interferometers for sensing applications [29,30,31]. These studies focus on improving the interferometric signal in terms of output power characteristics, yet they do not study the effect EOF has on linewidth, which is crucial for other applications such as coherent communications [32] or frequency metrology [33]. Furthermore, in all works exploring dual external cavities, both external reflections return to the same side of the laser cavity and thus have the same propagation direction. Moreover, phase shifts in the external phase cavity are not accounted for, thus ignoring, e.g., possible phase shifts at the mirrors.
The established line of thought relies on the key assumption that feedback is introduced from only one side of the laser cavity. Current integration technologies make this assumption obsolete, as they allow for arbitrarily complex design geometries with a variety of functionalities, such as tunability and modulation, while maintaining narrow-linewidth performance [34,35,36,37,38]. Consequently, this work aims to extend the theoretical foundations of EOF to the case of feedback coupling into the laser cavity from both sides. This system is studied to obtain and analyze its dynamic rate equations. Furthermore, the frequency noise power spectral density and subsequently the intrinsic linewidth’s dependence on feedback is computed. The Lang–Kobayashi approach is used [39], where an additional term to account for the extra feedback term is included. In a similar fashion, to obtain analytical solutions both small-signal and weak feedback conditions are assumed. The obtained equations are then numerically solved. Results show that feedback-insensitivity can be achieved by tuning the feedback parameters. In contrast with previous works, this can be accomplished in a monolithic platform by including a single active component in a straight-forward geometric configuration. As such, design complexity is reduced which enables devices with a compact form factor. Furthermore, existing methods for laser fabrication are suitable for realizing the proposed device. In particular, the proposed laser system can be made with mature photonic integrated components available in generic foundry platforms [40,41], thus creating a path towards devices with high technology-readiness level while maintaining low cost and size.

2. Rate Equations Model

This section includes a derivation of the dynamic equations of a laser cavity with EOF coupled into the laser cavity from both sides. Starting from the Lang–Kobayashi model [39], the lasing frequency and threshold gain shifts due to feedback are obtained, as well as an analytical expression for the frequency noise power spectral density from which the change of intrinsic linewidth can be computed.
This work proposes a revised laser system, as shown in Figure 1. The laser cavity of length L is delimited by two mirrors with complex reflection coefficients ρ 1 and ρ 2 , respectively. Assuming two interfaces at each side of the main cavity, two additional back-reflections ( ρ 1 , ext and ρ 2 , ext ) are considered and accounted for by computing effective reflection coefficients. The case of weak feedback is discussed here, for which:
| ρ j , ext | 1 .
The following parameters are introduced:
κ j 1 ρ j 2 ρ j ρ j , ext t cav
ϕ j ω FB t j + ϕ m j
δ ω ω FB ω ref
t cav = 2 L / v g ,
with j = 1 , 2 , where κ j is the coupling coefficient; t cav is the cavity round-trip time, with the group velocity v g ; ϕ j is the phase delay due to the external cavities determined by the external round-trip time t j , the lasing frequency in the presence of feedback ω FB , and a phase shift at the external mirrors ϕ m j ; and ω ref is the free running laser frequency.
The first step is extracting the lasing conditions of the proposed laser system.

2.1. Lasing Conditions

By analyzing the slowly-varying amplitude of the complex electric field, effective reflection coefficients ρ j eff can be obtained:
ρ j eff ρ j 1 + κ j t cav e ± i ϕ j ,
where the plus and minus signs corresponds to j = 1 (left mirror) and j = 2 (right mirror), respectively. The additional reflection influences the lasing condition, as shown in Appendix A, and thus the threshold gain and lasing frequency. To calculate them, the following definitions are convenient:
G FB Γ g FB v g
G ref Γ g ref v g
Δ ϕ m ϕ m 1 ϕ m 2 ,
γ H 1 + α H 2
θ H arctan α H
where Γ is the confinement factor; g FB and g ref are the threshold gain coefficient with and without feedback, respectively; and α H is the linewidth enhancement factor [42]. From Appendix B, it is possible to obtain:
G FB G ref δ G ( A 16 ) 2 κ 2 cos ( ϕ 2 ) 2 κ 1 cos ( ϕ 1 )
δ ω ( A 23 ) γ H κ 2 sin ( ϕ 2 + θ H ) + κ 1 sin ( θ H ϕ 1 ) .
The relation between the right-hand terms determines the shift in threshold gain and lasing frequency with feedback. When comparing with previous literature regarding dual external cavity lasers, a sign difference is observed with respect to the lasing frequency shift equations found in References [28,29], arising from the different propagation direction of the additional feedback term. A clearer contrast is present with the lasing frequency shift shown in ref. [30], where a single term accounts for both external reflections. Furthermore, a phase shift at the external mirrors ( ϕ m j ) is not considered in the mentioned sources, thus the feedback terms included in both the excess gain and lasing frequency variations are distinct with respect to those in Equation (5).
Given the transcendental form of Equation (5b), a numerical analysis under different feedback conditions is studied in Section 3. Nevertheless, an analytical solution for lasing frequency stability can be found for the condition:
δ ω = 0 .
Under this condition, Equation (5b) becomes:
κ 2 sin ( ϕ 2 + θ H ) = ( 6 ) κ 1 sin ( θ H ϕ 1 ) ,
which is determined by the feedback parameters κ j and ϕ j , the latter being dependent on the time delay t j as well as the lasing frequency. Finding a stable solution that does not depend on the lasing frequency is of particular interest, as it can be advantageous for tunable lasers and their numerous applications. With the following assumption:
κ 2 = κ 1 κ ,
Equation (7) can be rewritten as:
sin ( ϕ 2 + θ H ) = sin ( θ H ϕ 1 ) = sin ( ϕ 1 θ H ) ϕ 2 + θ H = ϕ 1 θ H + 2 m π .
Without loss of generality, the parameter m is set to m = 0 , thus:
2 θ H = ϕ 1 ϕ 2 = ( 2 b ) = ( 4 c ) ω FB ( t 1 t 2 ) + Δ ϕ m 2 θ H Δ ϕ m = ω FB ( t 1 t 2 ) .
Choosing an equal time delay in both external cavities:
t 2 = t 1 t ext ,
sets the right hand term of Equation (10) to zero, so that:
2 θ H = Δ ϕ m .
This result suggests that by tuning the phase in the external cavities, so that condition (12) is met, it is possible to obtain a feedback-insensitive lasing frequency. An active method to tune the phase is however required as α H is dependent on laser parameters, such as carrier density and wavelength [43]. This can be managed by, for example, phase shifters, which are mature and widely used components that can be included on-chip in a laser. The shown stable solution thus requires meeting the conditions (8), (11) and (12), which constrain the feedback parameters of one side of the cavity with respect to those of the other side, but do not restrict their absolute value. Nevertheless, if conditions (8) and (11) are not met, solutions for stable performance become frequency dependent. Such a case would thus only be satisfied for certain lasing frequency values for a given set of feedback parameters, which can potentially yield unstable solutions for other frequencies.
Regarding the variations in threshold gain, under conditions (8) and (11) Equation (5a) becomes:
δ G = 2 κ cos ( ω ref t ext + ϕ m 1 ) + cos ( ω ref t ext + ϕ m 2 ) ) ,
which vanishes if:
cos ( ω ref t ext + ϕ m 1 ) = cos ( ω ref t ext + ϕ m 2 ) cos ( ω ref t ext + ϕ m 1 ) = cos ( ω ref t ext + ϕ m 2 π ) ω ref t ext + ϕ m 1 = ω ref t ext + ϕ m 2 π + 2 m π m Z m = 0 Δ ϕ m = π .
This condition, while different than condition (12), also yields stability regardless of lasing frequency in this case for the threshold gain. Both cases are studied numerically in Section 3.

2.2. Rate Equations

In order to obtain the frequency noise (FN) power spectral density (PSD), and from it the laser intrinsic linewidth, the laser rate equations for the intensity and phase, as well as one for the carrier number, need to be studied. The former two can be extracted from the dynamic equations for the field inside the laser cavity, following the Lang-Kobayashi [39] approach. Its full derivation is shown in Appendix C. Furthermore, Langevin noise terms are included to account for shot noise fluctuations. The following definitions are useful to simplify notation:
A ( t ) = S ( t ) e i ϕ ( t )
S j ± κ j S ( t ± t j )
ϕ t j ± ϕ ( t ± t j )
Δ Φ 1 + = ( 15 c ) ϕ ( t ) ϕ t 1 + ϕ 1
Δ Φ 2 = ( 15 c ) ϕ ( t ) ϕ t 2 + ϕ 2
Δ G G FB τ ph 1 ,
with j = 1 , 2 relating to the EOF components from the left and right, respectively, and τ ph is the photon decay time, which accounts for cavity and mirror losses. The field amplitude A is assumed to be slowly varying, where S is the photon number inside the laser cavity and ϕ is the phase of the field. With these definitions, the rate equations of the system can be written as:
S ˙ = ( A 27 ) S Δ G + 2 S 2 S cos ( Δ Φ 2 ) + 2 S 1 + S cos ( Δ Φ 1 + ) + R sp + F S
ϕ ˙ = ( A 28 ) α H Δ G 2 δ ω S 2 S sin Δ Φ 2 S 1 + S sin Δ Φ 1 + + F ϕ
N ˙ = [ 3 ] I G S ( t ) N τ sp 1 + F N ,
where I is the effective rate of injected current (in electrons), τ sp is the carrier lifetime, and R sp is the spontaneous recombination rate. This system is compatible with those presented in refs. [28,31], by accounting for the sign change in the delayed feedback term corresponding to the left external mirror, due to it having the opposite propagation direction. The steady state solution of Equation (16) can be seen in Equation (A33). The Langevin noise sources F S ( t ) , F ϕ ( t ) and F N ( t ) satisfy [44]:
F i ( t ) = 0
F i ( t 1 ) F j ( t 2 ) = 2 D i j δ ( t 1 t 2 ) with i , j = S , ϕ or N ,
where:
D S S = R sp S ; D ϕ ϕ = R sp 4 S ; D N N = R sp S + N τ sp 1 ; D S N = R sp S ,
are standard diffusion coefficients. A usual approach for solving the system from Equation (16) involves a small-signal analysis. Small deviations from a steady-state value are assumed:
S S 0 + S Δ = S 0 + e i Ω t S 0 p ( Ω ) d Ω with S 0 S Δ
ϕ ϕ Δ = e i Ω t ϕ 0 p ( Ω ) d Ω
N N 0 + N Δ = N 0 + e i Ω t N 0 p ( Ω ) d Ω with N 0 N Δ ,
where the steady state value of the phase is assumed to be zero. The full linearization of the rate equations is shown in Appendix D, which uses the following definitions:
κ j c κ j t j cos ϕ j
κ j s κ j t j sin ϕ j
K s κ 2 s + κ 1 s
K c 1 + κ 2 c κ 1 c
ζ S R sp / S 0
a g = Γ v g a
G i a g ( N i N tr )
τ e 1 a g S 0 + τ sp 1 ,
where a linear approximation for the gain has been introduced, with a the differential gain coefficient and N tr the number of electrons at transparency. Applying the Fourier transform to Equation (A36a–c), the following system of equations is obtained in the frequency domain:
i Ω K c S 0 p = ( A 36 a ) a g S 0 N 0 p ζ S S 0 p 2 i Ω S 0 K s ϕ 0 p + F ^ S
2 i Ω K c ϕ 0 p = ( A 36 b ) α H a g N 0 p + i Ω K s S 0 S 0 p + 2 F ^ ϕ
i Ω N 0 p = ( A 36 c ) τ e 1 N 0 p G 0 S 0 p + F ^ N ,
where the unknowns S 0 p , ϕ 0 p , N 0 p , and F ^ S , F ^ ϕ and F ^ N depend on the Fourier angular frequency Ω . These equations are the first step to obtain the FN PSD.

2.3. Power Spectral Density and Laser Intrinsic Linewidth

The next step is to calculate the FN PSD, and from it the laser intrinsic linewidth. By solving the system in Equation (21) it is possible to find an expression for the PSD [45]:
S f ( 1 ) ( Ω ) = Ω 2 2 π 2 | ϕ 0 p ( Ω ) | 2 ,
which, using the following definitions:
F 0 ζ S 2 + K c 2 τ e 2 2 K c a g G 0 S 0
F 1 K s 2 + K c 2
Λ 4 4 F 1 D ϕ ϕ
F 2 K c α H + K s
F 3 K c α H K s
Δ 4 F 1 2 ,
and:
Λ 2 a g 2 F 2 2 D N N + 4 D ϕ ϕ K s 2 τ e 2 + F 0 2 K s S 0 D S S a g τ e 1 F 2 ζ S α H α H G 0
Λ 0 α H 2 a g 2 D S S ζ S 2 + G 0 2 + 2 ζ S G 0 + ζ S N τ sp 1 + τ e 1 ζ S + a g G 0 S 0 2 4 D ϕ ϕ
Δ 2 K c 2 ζ S 2 + τ e 2 F 1 2 2 F 1 a g G 0 S 0 F 3
Δ 0 a g G 0 S 0 F 3 + K c ζ S τ e 1 2 ,
can be written as:
4 π 2 S f ( 1 ) = ( 23 ) = ( A 46 ) Λ 4 Ω 4 + Λ 2 Ω 2 + Λ 0 Δ 4 Ω 4 + Δ 2 Ω 2 + Δ 0 .
This is an analytical solution for the FN PSD of the isolated laser system. From the following expression [46]:
S f ( 1 ) ( Ω 0 ) = 2 π Δ f ,
which is valid for a Lorentzian lineshape, the intrinsic linewidth can be obtained. As shown in Appendix F:
F Δ f Δ f 0 1 + α H 2 = ( A 49 ) 1 + γ H κ 2 t 2 cos ϕ 2 + θ H γ H κ 1 t 1 cos ϕ 1 θ H 2 ,
where Δ f 0 is the Schawlow–Townes linewidth [47]. The expression found for the intrinsic linewidth has two feedback terms, one contribution from each side, with a sign that depends on ϕ 1 and ϕ 2 . Recalling from Equation (2b) that these quantities are a function of t j and ϕ m j , a proper design of the laser can yield linewidth stability or a reduction of the intrinsic linewidth with respect to the case of one-sided feedback. This is further explored using a numerical analysis in Section 3. Additionally, it is possible to find an analytical expression for Equation (26) under the conditions for frequency stability, namely conditions (8), (11) and (12). Using this assumptions in Equation (26):
F = 1 + γ H κ t ext cos ϕ 2 + θ H cos ϕ 1 θ H 2 .
Taking a closer look at the feedback terms yields:
cos ϕ 2 + θ H = ( 12 ) ω ref t ext + ϕ m 1 2 θ H + θ H = cos ϕ 1 θ H .
Using Equation (28) in Equation (27) yields a value of F = 1 which indicates that, under the assumed conditions, the intrinsic linewidth is insensitive to feedback. This result is significant as under the same condition the frequency is also feedback-insensitive, as shown in Section 2.1, regardless of lasing frequency. It is worth noting however that weak feedback is assumed in this analysis, with which the upper bound for feedback strength under which these equation are valid is not established. Nevertheless, achieving stability in the full frequency domain even under this condition is an improvement with respect to the single feedback case.

3. Numerical Study

Laser stability is studied by numerically evaluating the equations for the shift in lasing frequency, threshold gain and intrinsic linewidth under the revised EOF conditions, namely Equations (5) and (26), under different feedback parameters. Particular attention is given to the previously analyzed case under conditions (8), (11) and (12), which shows feedback-insensitive solutions for the lasing frequency and intrinsic linewidth. System tolerances to these conditions are explored by varying each while keeping the other two fixed. The simulated equations are plotted as a function of the unperturbed laser frequency multiplied by t 2 , which represents the phase delay in the right external cavity for the free running laser frequency. It is kept between 0 and 1 (i.e., ω ref t 2 ( 0 , 2 π ] ) given the periodicity of the numerically solved functions. Furthermore, without loss of generality, ϕ m 2 is kept fixed at zero so that the value of Δ ϕ m is selected by varying ϕ m 1 . Additionally, simulations assume α H = 3 . This value is compatible with measurements for semiconductor lasers [48], and thus meeting condition (12) requires that Δ ϕ m = 2 θ H 2.5 . A summary of the values used in the numerical solutions is given in Table 1.
Results are compared with the case with feedback from a single side, where:
κ 1 = 0 .
In this case, as shown in [7], as feedback strength increases solutions for the lasing frequency become multi-valued. This gives rise to instabilities in the system such as mode hopping or coherence collapse regimes. The separation between single-valued and multi-valued solutions is related to the coefficient:
C = γ H κ 2 t 2 ,
where C = 1 is the critical value that separates both behaviours.
Numerical solutions of the proposed system under conditions (8) and (11) are thus compared to the single feedback case for three cases:
Case 1: C = 0.5 . This represents the single-feedback case with a single solution, and results for various values of Δ ϕ m are shown in Figure 2. Column A, B and C show the numerical solutions for the variations in lasing frequency, threshold gain and intrinsic linewidth, given by Equations (5) and (26), respectively. By plotting the logarithm of the latter, negative values correspond to linewidth narrowing. The blue and orange plots represent the double-feedback and the single-feedback case, respectively, and these labels are maintained throughout the document.
The upper row shows the case where Δ ϕ m = 0 , i.e., additional phase shifts at, for example, the external mirrors ( ϕ m j ) are kept at zero. Linewidth narrowing for a wide range of frequencies can be observed, evidenced, as mentioned, by l n ( F ) < 1 . As an example, when f ref t 2 = 0.75 , a 74% and a 42% reduction in linewidth is seen with respect to the free running laser and the single feedback case, respectively. Additionally, the signal is singled valued in the full domain, yet it is close to the critical point where multi-valued solutions arise. As Δ ϕ m increases, the amplitude of the lasing frequency shift is reduced until becoming zero for all values when condition (12) is met, as expected from previous analysis. Under this condition the intrinsic linewidth does not experience fluctuations either, and the amplitude of the threshold gain fluctuations is lower than in the single feedback case, indicating better stability across the three analyzed parameters with respect to the single feedback case. For the case of Δ ϕ m = π , the threshold gain shows no fluctuations as predicted by Equation (14), and while the lasing frequency and intrinsic linewidth fluctuations are no longer zero, they are less pronounced than in the single feedback case. Further increases in Δ ϕ m show an increase in the fluctuations across all functions, and for ϕ m > 6 multi-valued solutions arise.
Case 2: C = 1 . This represents the limiting case between single and multivalued solutions in the single-feedback case. Results for various values of Δ ϕ m are shown in Figure 3. As expected, the threshold gain is stable for Δ ϕ m = π , and meeting condition (12) results in a stable lasing frequency and intrinsic linewidth. As Δ ϕ m deviates from these optimal points in either direction, the amplitude of fluctuations increase until reaching multi-valued solutions for Δ ϕ m < 1.5 and Δ ϕ m > 3.5 . Comparing these results with the previous case shows that as feedback increases, the single valued solutions become more sensitive to the value of Δ ϕ m .
Case 3: C = 1.3 . This represents the single-feedback case with multi-valued solutions. Results for various values of Δ ϕ m are shown in Figure 4. In the single feedback case, multi-valued solutions are present in a given frequency range, and this span increases with increasing feedback strength. The multi-valued characteristics are evidenced experimentally with unstable regimes characterized by mode hopping and eventually coherence collapse for sufficiently high feedback. In contrast, the system proposed in this work shows that by tuning the value of Δ ϕ m to meet condition (12), even with increasing feedback it is possible to achieve stable performance regardless of frequency. In the case shown in Figure 4 for C = 1.3 , single valued solutions can be found for Δ ϕ m ( 1.5 , 3.3 ) which is equivalent to a phase variation of more than 90 . Still, comparing with previous cases it is possible to see that as feedback increases, the single valued solutions tolerance with respect to Δ ϕ m decreases. Nevertheless, it is an improvement with respect to the single feedback case which shows no single value solutions across all frequencies for C > 1 .
Taking all cases into account, it can be seen that linewidth narrowing regions are present for all cases of Δ ϕ m analyzed. The only exception are the stable cases when meeting the three conditions (8), (11) and (12), yet the feedback-insensitivity provided by this case is also beneficial. Selecting an appropriate Δ ϕ m can thus be used to harness linewidth narrowing properties at a desired frequency.
Furthermore, system tolerances to conditions (8) and (11) are studied, while maintaining condition (12). Results for up to a 20% deviation are shown in Figure A2,Figure A3 and Figure A4 for cases 1, 2 and 3, respectively, where single valued solutions are observed in all cases in the full frequency span. Results show a high tolerance with respect to feedback strength. Looking at case 3, single valued solutions are obtained for κ 1 / κ 2 ( 0.2 , 1.7 ) . While the system is no longer feedback-insensitive, results evidence single-valued solutions that are robust with respect to condition (8). The system is more sensitive with respect to (11), with single valued solutions achieved within t 1 / t 2 ( 0.47 , 1.2 ) . Nevertheless, tolerances become once again stricter as feedback increases for both parameters, thus laser design is of paramount importance and should focus on meeting the discussed stability conditions. In particular, choosing equal lengths for both external cavities should suffice to meet condition (11). Fabrication tolerances in foundry processes are the main limiting factor for time delay accuracy. To meet condition (8), a possible approach is to merge the output of the two external cavities into a single one using a coupler, which can be included in the laser on-chip. Finally, as mentioned before, condition (12) can be met by using a phase shifter, which is a mature component in active platforms.
All in all, results demonstrate that, with proper design of the laser cavity, conditions (8), (11) and (12) can be met, with which it is possible to obtain lasing frequency and intrinsic linewidth insensitivity to feedback.

4. Discussion

The existing mature photonic integration fabrication processes are very flexible with respect to device geometry. They are however limited by a lack of commercially available on-chip isolators, and thus new approaches are required to minimize the effect feedback has on laser stability. The current work proposes a theoretical extension of laser dynamics under EOF by considering two external reflections, one from each side of the cavity. The proposed analysis yields new laser dynamic equations. These are numerically solved, which show the existence of a stable regime with high feedback tolerances. Similar results have been seen in previous double external cavity schemes [28], with EOF on the same side of the laser cavity, which showed robust stabilization of the chaotic regime by tuning the external feedback parameters. Yet other feedback regimes are not explored, and neither analytical stability conditions nor dependence with free laser running frequency are reported. The present analysis shows that feedback-insensitivity is achievable for a wide range of feedback strengths, under conditions that can be met with current laser fabrication processes and components, such as phase shifters for meeting condition (12) and couplers for condition (8).
It has previously been shown, for the single feedback case, that tuning the feedback phase can result in linewidth narrowing [13], however this still requires low feedback levels and a precise phase shift which can suffer variations due to external parameters, such as temperature or driving currents. In the approach proposed in this work, the stability conditions do not require specific values, instead relating the feedback parameters from one side to those from the other side. For example, condition (11) only implies equal round-trip times at the external cavities, regardless of value. This allows for additional flexibility in the feedback parameters and gives versatility to the device. Furthermore, this method allows for feedback-insensitivity across the full spectra, which is not seen in the single feedback case. This is of particular importance for tunable lasers, as all lasing frequencies are thus equally affected. Additionally, it relaxes the need for an isolator, reducing the cost and size of packaging processes. Another significant improvement of the proposed method with respect to the single feedback case is the increase in feedback tolerance: by choosing feedback parameters close to the stability conditions of Equations (11) and (12), higher levels of feedback strength are allowed without seeing multi-valued solutions, which results in experimentally seen mode-hopping. As a weak feedback approximation is used, the upper bound for feedback tolerance cannot be extracted from this analysis. Despite this, even under this approximation, tolerances are higher than that of the single feedback case. Furthermore, this system has a high tolerance to deviations from the optimal stability conditions as analyzed in the previous section. Linewidth narrowing can be achieved in these cases for certain frequency values, which can be tuned by selecting appropriate feedback parameters, as was the case for single feedback conditions, while maintaining stable solutions.
Finally, while the dynamics under consideration are complex, the laser system itself involves a straight-forward configuration using widely used on-chip components, which are available in generic foundry platforms. Previously studied methods to reduce feedback sensitivity include resourceful yet intricate designs. The proposed system is, in contrast, potentially easier to design, fabricate and characterize. An experimental study of this laser system is essential to validate the obtained results, and more importantly to explore the limitations of the model, and is the next step for a more comprehensive understanding of the proposed system. It is worth noting that, for the single feedback case, the Lang-Kobayashi approach yields results compatible with experimental data [49,50,51], which suggests that future experimental realization of the method here presented will be compatible with its theoretical predictions.

5. Conclusions

This work explores an extension of the theoretical background of EOF. By assuming that feedback couples into the laser cavity from both sides, new dynamic equations are found for the lasing frequency, the threshold gain and the intrinsic linewidth. These are numerically evaluated to analyze laser stability. Solutions with linewidth reduction are observed, where tuning of the feedback parameters yielded a 74% and a 42% reduction with respect of the free running laser and the single feedback case, respectively. Results also show the existence of a stable solution, with feedback-insensitive lasing frequency and intrinsic linewidth, regardless of the lasing frequency. This case is obtained by tuning the phase of the feedback field, for external cavities with equal lengths and coupling factors. Furthermore, the feedback-insensitive case exists regardless of the feedback strength, within a weak feedback approximation, which is an major improvement with respect to the single feedback case. Additionally, the stability conditions show good tolerances with respect to all feedback parameters, albeit they become stricter as the feedback strenght increases. Choosing feedback parameters close to the feedback-insensitive conditions ensures stable solutions that are feedback tolerant. Finally, the proposed system relies on few components in straight-forward configurations, and the stable conditions can be met with mature components available in generic foundry platforms. This enables close to market, low cost, feedback-tolerant semiconductor lasers which have direct applications in multiple fields that rely on stable laser sources, such as coherent communications and spectroscopy.

Author Contributions

Formal analysis, M.F.B.; Funding acquisition, N.V.; Project administration, N.V.; Resources, N.V.; Supervision, N.V.; Writing—original draft, M.F.B.; Writing—review & editing, M.F.B. and N.V. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge support from Independent Research Fund Denmark.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Amplitude and Phase Conditions

To obtain the revised lasing conditions resulting from the additional feedback component, extracting the amplitude and phase of the effective reflection coefficients is needed. In polar notation:
ρ j eff = ρ j eff e i φ j ,
where the magnitude is computed as:
ρ j eff 2 / ρ j 2 = ( 3 ) 1 + κ j t cav cos ( ϕ j ) 2 + κ j t cav sin ( ϕ j ) 2 = 1 + 2 κ j t cav cos ( ϕ j ) + κ j 2 t cav 2 .
Assuming condition (1) of weak external feedback, the last term on the right-hand side of Equation (A2) is neglected, resulting in:
ρ j eff ( ω FB ) ρ i = 1 + 2 κ j t cav cos ( ϕ j ) ( 1 ) 1 + κ j t cav cos ( ϕ j ) .
The phase of the effective reflection coefficient is extracted from:
φ j = arctan ρ j eff / ρ j eff = ( 3 ) arctan ± κ j t cav sin ( ϕ j ) 1 + κ j t cav cos ( ϕ j ) ( 1 ) arctan [ ± κ j t cav sin ( ϕ j ) ] ( 1 ) ± κ j t cav sin ( ϕ j ) .
The fields traveling forward and backward in the laser cavity, E f and E b shown in Figure A1, can now be related by the effective reflection coefficients:
E f ( z = 0 ) = ρ 1 eff E b ( z = 0 )
E b ( z = L ) = ρ 2 eff E f ( z = L ) .
Using the propagation constant:
β n ω / c ,
with n the effective refractive index of the lasing mode and c the speed of light in vacuum, the fields can be written as:
E f = A f e i β z + 1 2 ( Γ g α ) z
E b = A b e i β ( L z ) + 1 2 ( Γ g α ) ( L z ) ,
where g is the gain coefficient and α is the attenuation coefficient. Replacing Equation (A5) into Equation (A7):
E f 0 = ( A 1 ) ρ 2 eff e i φ 2 A b e i β L + ( Γ g α ) L / 2
E b 0 = ( A 1 ) ρ 1 eff e i φ 1 A f e i β L + ( Γ g α ) L / 2 ,
and inserting Equation (A8a) into Equation (A8b) results in:
1 = ρ 1 eff e i φ 1 ρ 2 eff e i φ 2 e 2 i β L + ( Γ g α ) L = ρ 1 eff ρ 2 eff e i ( 2 β L φ 1 φ 2 ) e ( Γ g α ) L .
Once lasing has been established, the gain assumes its threshold value:
g = g FB ,
where g FB is the threshold gain with feedback. Thus, Equation (A9) yields a lasing condition for the amplitude:
1 = ρ 1 eff ρ 2 eff e Γ g FB α L = ( A 10 ) ( A 3 ) ρ 2 ρ 1 1 + κ 2 t cav cos ( ϕ 2 ) 1 + κ 1 t cav cos ( ϕ 1 ) e Γ g FB α L ,
and the phase:
2 π m = 2 β L φ 2 φ 1 = ( A 4 ) 2 β L + κ 2 t cav sin ( ϕ 2 ) κ 1 t cav sin ( ϕ 1 ) , m Z .
where the influence of feedback gives rise to two terms in Equations (A11) and (A12), one from each side. The new lasing conditions result in a variation of the lasing frequency and threshold gain of the system, and thus have to be studied to determine the laser dynamics.
Figure A1. Schematic of the effective cavity of the laser, resulting from calculating effective reflection coeffitients.
Figure A1. Schematic of the effective cavity of the laser, resulting from calculating effective reflection coeffitients.
Photonics 09 00043 g0a1

Appendix B. Threshold Gain Reduction and Lasing Frequency Shift

Under feedback from both sides of the laser cavity, new lasing conditions are found which subsequently result in a shift of the laser threshold gain and lasing frequency with respect to the case without feedback, in which:
κ j = 0 .
In this case, the amplitude condition from Equation (A11) becomes:
1 = ( A 13 ) ρ 2 ρ 1 e Γ g th α L .
Using the expansion:
ln 1 + x x ,
the threshold gain reduction due to feedback can be found by computing the ratio between Equations (A11) and (A14):
1 = ρ 2 ρ 1 1 + κ 2 t cav cos ( ϕ 2 ) 1 + κ 1 t cav cos ( ϕ 1 ) e Γ g FB α L ρ 2 ρ 1 e Γ g ref α L = 1 + κ 2 t cav cos ( ϕ 2 ) 1 + κ 1 t cav cos ( ϕ 1 ) e Γ g FB g ref L ( 2 ) G FB G ref 2 κ 2 cos ( ϕ 2 ) 2 κ 1 cos ( ϕ 1 ) .
The relation between the right hand terms determines the threshold gain reduction, as discussed in Section 2.1.
The phase lasing condition from Equation (A12) yields the lasing frequency shift equation. Consider the following definitions related to the effective refractive index [52]:
n n + i n
n g n + ω n ω
n c G 2 ω v g ,
α H Δ n / Δ n n N = n n n N
= ( A 17 a ) = ( A 17 c ) = ( A 17 d ) G N α H c 2 ω v g ,
where n g is the group refractive index. To find the lasing frequency shift equation, calculating the change in β is first needed:
c δ β = ( A 6 ) δ ( n ω ) = ω δ n + n δ ω = ( 2 c ) ω n N N N th + n ω δ ω + n δ ω = ( A 17 b ) = ( A 17 e ) G FB G th 2 v g α H c + n g δ ω = ( 5 a ) κ 2 cos ( ϕ 2 ) + κ 1 cos ( ϕ 1 ) α H c v g + n g δ ω ,
with which:
2 L δ β = ( 5 a ) κ 2 cos ( ϕ 2 ) + κ 1 cos ( ϕ 1 ) α H t cav + t cav δ ω .
Furthermore, using:
sin arctan ( x ) = x 1 + x 2
cos arctan ( x ) = 1 1 + x 2
sin ( x ± y ) = sin ( x ) cos ( y ) ± cos ( x ) sin ( y )
cos ( θ H ) = ( 4 e ) = ( A 20 b ) = ( 4 d ) 1 / γ H ,
the following can be computed:
α H cos ( ϕ j ) = ( A 20 a ) = ( 4 d ) = ( 4 e ) γ H sin θ H cos ( ϕ j ) = ( A 20 c ) γ H sin ( θ H ± ϕ j ) cos ( θ H ) sin ( ϕ j ) = ( A 20 d ) γ H sin θ H ± ϕ j sin ( ϕ j ) / γ H α H cos ( ϕ j ) ± sin ( ϕ j ) = γ H sin θ H ± ϕ j .
Finally, using the phase condition in Equation (A12), and assuming without generality loss that m = 0 , it is possible to compute:
2 π m = ( A 19 ) t cav δ ω + κ 2 t cav α H cos ( ϕ 2 ) + sin ( ϕ 2 ) + κ 1 t cav α H cos ( ϕ 1 ) sin ( ϕ 1 ) m = 0 δ ω = κ 2 α H cos ( ϕ 2 ) + sin ( ϕ 2 ) + κ 1 α H cos ( ϕ 1 ) sin ( ϕ 1 ) δ ω = ( A 21 ) γ H κ 2 sin ( ϕ 2 + θ H ) + κ 1 sin ( ϕ 1 θ H ) ,
which describes the lasing frequency shift as a function of the feedback parameters κ 1 , κ 2 , ϕ 1 and ϕ 2 , Similarly to the threshold gain reduction, the interaction between the two right hand terms determines the lasing frequency stability which is discussed in Section 2.1.

Appendix C. Deriving the Rate Equation for the Intensity and Phase

To further inspect the laser dynamics, the rate equations for the intensity and phase must be studied. Considering Equation (15), and a slowly varying electric field given by:
E ( t ) = A ( t ) e i ω FB t ,
and following the approach from Lang and Kobayashi [39], the laser field equation which considers EOF from both sides of the laser cavity can be written as:
E ˙ = i ω ref + Δ G 1 i α H 2 E ( t ) + κ 2 E ( t t 2 ) + κ 1 E ( t + t 1 ) ( A 23 ) d A ( t ) e i ω FB t d t = ( 15 ) e i ω FB t i ω ref + Δ G 1 i α H 2 A ( t ) + κ 2 A ( t t 2 ) e i ϕ 2 + κ 1 A ( t + t 2 ) e i ϕ 1 A ˙ ( t ) = ( 2 c ) = ( 2 b ) i δ ω + Δ G 1 i α H 2 A ( t ) + κ 2 A ( t t 2 ) e i ϕ 2 + κ 1 A ( t + t 2 ) e i ϕ 1 .
The last two right-hand terms appear as a result of the imposed feedback conditions, each term to account for feedback on each side of the cavity. In the case without feedback the lasing frequency become ω FB = ω ref and κ j = 0 , thus recovering the no-feedback field Equation [42]. The slowly varying field amplitude A can be modeled as Equation (15a), and thus the rate equations for the photon number S and phase ϕ can be found using:
S ˙ = d [ A A * ] d t = A A * ˙ + A * A ˙
ϕ ˙ = 1 S A * A ˙ .
Replacing Equations (15a) and (A24) into Equation (A25a) the photon rate equation reads:
S ˙ = ( 15 d ) = ( 15 e ) S Δ G + 2 S 2 S cos ( Δ Φ 2 ) + 2 S 1 + S cos ( Δ Φ 1 + ) .
In the case of the phase, its rate equation comes from replacing Equations (15a) and (A24) into Equation (A25b):
ϕ ˙ = ( 15 ) Δ G 2 α H δ ω S 2 S sin Δ Φ 2 S 1 + S sin Δ Φ 1 + .
Equations (A26) and (A27) are the amplitude and phase rate equations for the laser system proposed in this work. These are the starting point to compute the frequency noise PSD, and extract the intrinsic linewidth.

Appendix D. Small-Signal Analysis

To find the FN PSD, the system shown in Equation (16) is to be solved. This is done using the small-signal analysis proposed in Equation (19). Assuming a narrow-linewidth laser, i.e., a long coherence time with respect to the external cavity lenghts:
t ext < t coh ,
the following approximation is valid:
Ω t ext 1 .
By linearizing the following expressions:
S ( t ± t j ) S ( t ) ( A 28 ) = 1 ± S ˙ S t j ( 19 a ) 1 ± i Ω S Δ S 0 t j ( A 29 ) 1 ± i Ω S Δ 2 S 0 t j
ϕ ϕ ( t ± t j ) ( A 28 ) ϕ ϕ t j ϕ ˙ = ( 19 b ) i t j Ω ϕ Δ ,
it is possible to rewrite Equation (16) as:
i Ω S Δ = ( A 30 ) S Δ G + R sp + 2 κ 2 t cav S 1 i Ω S Δ 2 S 0 t 2 cos i Ω ϕ Δ t 2 + ϕ 2 + 2 κ 1 S 1 + i Ω S Δ 2 S 0 t 1 cos i Ω ϕ Δ t 1 + ϕ 1 + F S
i Ω ϕ Δ = ( A 30 ) α H Δ G 2 δ ω κ 2 t cav 1 i Ω S Δ 2 S 0 t 2 sin i Ω ϕ Δ t 2 + ϕ 2
+ κ 1 1 + i Ω S Δ 2 S 0 t 1 sin i Ω ϕ Δ t 1 + ϕ 1 + F ϕ
i Ω N Δ = ( 19 ) = ( A 30 ) I G FB S N τ sp 1 + F N .
Solving Equation (A31) requires the steady-state solution of Equation (16). Under stationary conditions:
S ˙ = 0 S ( t ) = S ( t ± t j ) = S 0
ϕ ˙ = 0 ϕ ( t ) = ϕ ( t ± t j )
N ˙ = 0 N ( t ) = N ( t ± t j ) = N 0 ,
the steady-state equations are:
τ ph 1 = G 0 + 2 κ 2 cos ϕ 2 + 2 κ 1 cos ϕ 1 + R sp S 0
δ ω = α H G 0 τ ph 1 2 κ 2 sin ϕ 2 + κ 1 sin ϕ 1
I = G 0 S 0 + N 0 τ sp 1 ,
where the Langevin noise terms are not included as their mean value is zero. Next, using Equation (A33) and the following expansions:
sin ( x + Δ ) sin ( x ) + Δ cos x
cos ( x + Δ ) cos ( x ) Δ sin x ,
Equation (A31a) can be rewritten as:
i Ω S Δ = ( 19 ) = ( A 33 ) ( 20 ) G 0 + a g N Δ G 0 + 2 κ 2 cos ϕ 2 + κ 1 cos ϕ 1 + R sp S 0 S 0 + S Δ + 2 { κ 2 S 0 + S Δ cos ( ϕ 2 ) i Ω ϕ Δ t 2 sin ϕ 2 κ 2 S 0 + S Δ i Ω S Δ 2 S 0 t 2 cos ( ϕ 2 ) i Ω ϕ Δ t 2 sin ϕ 2 + κ 1 S 0 + S Δ cos ( ϕ 1 ) i Ω ϕ Δ t 1 sin ϕ 1 + κ 1 S 0 + S Δ i Ω S Δ 2 S 0 t 1 cos ( ϕ 1 ) i Ω ϕ Δ t 1 sin ϕ 1 } + R sp + F S .
Simplifying this equation, and neglecting the quadratic terms yields:
i Ω S Δ ( 20 ) a g N Δ S 0 ζ S S Δ i Ω ϕ Δ 2 K s S 0 + i Ω S Δ ( κ 1 c κ 2 c ) + F S .
In a similar way, Equation (A31b) can be rewritten using Equations (A33) and (20):
i Ω ϕ Δ α H 2 G 0 + a g N Δ τ ph 1 α H G 0 τ ph 1 2 κ 2 sin ϕ 2 + κ 1 sin ϕ 1 κ 2 1 i Ω S Δ 2 S 0 t 2 sin ϕ 2 + i Ω ϕ Δ t 2 cos ϕ 2 + κ 1 1 + i Ω S Δ 2 S 0 t 1 sin ϕ 1 + i Ω ϕ Δ t 1 cos ϕ 1 + F ϕ 2 i Ω ϕ Δ = ( 20 a ) = ( 20 b ) α H a g N Δ + 2 ( κ 1 c κ 2 c ) i Ω ϕ Δ + i Ω S Δ S 0 K s + 2 F ϕ .
Finally, Equation (A31c) can be rewritten as:
i Ω N Δ = ( 19 ) = ( A 33 ) ( 20 ) G 0 S 0 G 0 + a g N Δ S 0 G 0 + a g N Δ S Δ τ sp 1 N Δ + F N i Ω N Δ = ( 20 h ) τ e 1 N Δ G 0 S Δ + F N .
The linearized rate equations of the laser under study are thus Equations (A35a)–(A35c), from which the power spectral density, and subsequently the intrinsic linewidth, can be computed.

Appendix E. Power Spectral Density

The next step is to find an expression for ϕ 0 p from which the FN PSD, and thus the laser intrinsic linewidth, can be computed. The following definitions are convenient:
A ϕ i Ω + τ e 1 i Ω K c + a g G 0 i Ω + τ e 1 S 0 + ζ S
2 A S i Ω K s i Ω + τ e 1 S 0 α H a g G 0
2 A N α H A ϕ + 2 S 0 A S i Ω + τ e 1 a g
B ϕ K c A ϕ α H a g G 0 S 0 K s + i Ω K s 2 i Ω + τ e 1 .
First, N 0 p is extracted from Equation (21c):
N 0 p = F ^ N G 0 S 0 p i Ω + τ e 1 .
Next, replacing Equation (A37) into Equation (21a) yields the expression for S 0 p :
i Ω K c + ζ S S 0 p = a g F ^ N G 0 S 0 p i Ω + τ e 1 S 0 2 i Ω S 0 K s ϕ 0 p + F ^ S A ϕ i Ω + τ e 1 S 0 p = a g S 0 F ^ N i Ω + τ e 1 2 i Ω S 0 K s ϕ 0 p + F ^ S S 0 p = ( A 36 a ) a g S 0 F ^ N 2 i S 0 Ω K s i Ω + τ e 1 ϕ 0 p + i Ω + τ e 1 F ^ S A ϕ .
Finally, inserting Equation (A37) and (A38) into Equation (21b) and grouping the terms with ϕ 0 p , F ^ S , F ^ ϕ and F ^ N yields:
2 i Ω B ϕ ϕ 0 p = ( A 36 ) α H a g G 0 S 0 + α H A ϕ + i Ω K s i Ω + τ e 1 a g i Ω + τ e 1 F ^ N + 2 A S F ^ S + 2 A ϕ F ^ ϕ = ( A 36 b ) α H A ϕ + 2 S 0 A S a g i Ω + τ e 1 F ^ N + A S F ^ S + 2 A ϕ F ^ ϕ ϕ 0 p = ( A 36 c ) A N F ^ N + A S F ^ S + A ϕ F ^ ϕ i Ω B ϕ .
With Equation (A39) it is possible to calculate an expression for Equation (22):
2 π 2 | B ϕ | 2 S f ( 1 ) = ( 22 ) Ω 2 | B ϕ | 2 ϕ 0 p ( Ω ) ϕ 0 p * ( Ω ) = ( A 39 ) Ω 2 | B ϕ | 2 A N F ^ N + A S F ^ S + A ϕ F ^ ϕ i Ω B ϕ A N F ^ N + A S F ^ S + A ϕ F ^ ϕ i Ω B ϕ * = A N F ^ N + A S F ^ S + A ϕ F ^ ϕ A N F ^ N + A S F ^ S + A ϕ F ^ ϕ * .
It can be seen from Equation (A36) that the coefficients A i are independent of time, and assuming an ergodic process they can be taken out of the average in Equation (A40), obtaining:
π 2 | B ϕ | 2 S f ( 1 ) = ( 17 b ) | A N | 2 D N N + | A S | 2 D S S + | A ϕ | 2 D ϕ ϕ + A S A N * + A N A S * D S N
4 π 2 | B ϕ | 2 S f ( 1 ) = ( A 36 ) = ( 23 ) ζ S α H a g 2 + Ω 2 a g 2 F 2 2 D N N + Ω 4 K c 2 + Ω 2 F 0 + τ e 1 ζ S + a g G 0 S 0 2 4 D ϕ ϕ + Ω 2 τ e 1 K s S 0 2 + 2 α H a g G 0 K s S 0 + Ω 4 K s S 0 2 + α H a g G 0 2 D S S 2 α H 2 a g 2 G 0 ζ S + Ω 2 K s S 0 a g τ e 1 F 2 ζ S α H D S S 2 π 2 S f ( 1 ) = ( 23 ) Λ 4 Ω 4 + Λ 2 Ω 2 + Λ 0 2 | B ϕ | 2 .
Furthermore, using the following definitions:
Δ 4 F 1 2
Δ 2 K c 2 ζ S 2 + τ e 2 F 1 2 2 F 1 a g G 0 S 0 F 3
Δ 0 a g G 0 S 0 F 3 + K c ζ S τ e 1 2 ,
the expression for B ϕ Equation (A36d) can be rewritten as:
B ϕ = a g G 0 S 0 F 3 Ω 2 F 1 + i Ω K c ζ S + τ e 1 F 1 + K c ζ S τ e 1 | B ϕ | 2 = Ω 2 K c ζ S + τ e 1 F 1 2 + a g G 0 S 0 F 3 Ω 2 F 1 + K c ζ S τ e 1 2 = ( A 43 c ) Ω 4 Δ 4 + Ω 2 K c ζ S + τ e 1 F 1 2 2 F 1 a g G 0 S 0 F 3 2 F 1 K c ζ S τ e 1 + Δ 0
= ( A 43 c ) Ω 4 Δ 4 + Ω 2 Δ 2 + Δ 0 ,
with which the expression of the FN PSD becomes:
4 π 2 S f ( 1 ) = ( 23 ) Λ 4 Ω 4 + Λ 2 Ω 2 + Λ 0 Δ 4 Ω 4 + Δ 2 Ω 2 + Δ 0 .

Appendix F. Expression for the Intrinsic Linewidth

The laser intrinsic linewidth can be found from Equation (25). Defining:
β A g τ e 1 ζ S a g G 0 S 0 ; δ A g ζ S G 0 ; Δ f 0 = R sp 4 π S 0 ,
where, following [3]:
δ A g 0 β A g ,
as δ A g < 10 2 , which accounts for shot noise in the generation and recombination of minority carriers, and β A g is inversely proportional to the laser power which above threshold becomes negligible. Starting from Equation (A42) and setting Ω = 0 as required by Equation (25):
4 π Δ f = Λ 0 / Δ 0 = ( A 43 c ) = ( 23 ) ζ S α H a g 2 R sp S 0 1 + G 0 2 ζ S 2 + 2 G 0 ζ S + N τ sp 1 + τ e 1 ζ S + a g G 0 S 0 2 R sp S 0 a g G 0 S 0 F 3 + K c ζ S τ e 1 2 Δ f Δ f 0 = α H G 0 2 ζ S 2 + G 0 2 + 2 G 0 ζ S + ζ S 2 N τ sp 1 R sp S 0 + β A g + 1 2 K c β A g + F 3 2 = β A g + 1 2 + α H 2 1 + 2 δ A g + δ A g 2 1 + N τ sp 1 R sp S 0 K c β A g + F 3 2
( A 47 ) Δ f Δ f 0 1 + α H 2 F 3 2 = ( 23 ) = ( 20 ) 1 + κ 2 t 2 cos ϕ 2 α H sin ϕ 2 κ 1 t 1 cos ϕ 1 + α H sin ϕ 1 2 = ( A 20 a ) = ( A 20 b ) 1 + γ H κ 2 t 2 cos ϕ 2 + θ H γ H κ 1 t 1 cos ϕ 1 θ H 2 .
The presence of EOF from both sides of the laser cavity results in two terms in the linewidth expression, one for each side, as was seen in the threshold gain reduction and lasing frequency shift due to feedback. This result is discussed in Section 2.3.

Appendix G. Suplementary Images: Tolerances

Figure A2. Numerical solutions for C = 0.5 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure A2. Numerical solutions for C = 0.5 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g0a2
Figure A3. Numerical solutions for C = 1 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure A3. Numerical solutions for C = 1 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g0a3
Figure A4. Numerical solutions for C = 1.3 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure A4. Numerical solutions for C = 1.3 under condition (12) for a ± 20 % variation of conditions (8) and (11). Double feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g0a4

References

  1. Sciamanna, M.; Shore, K.A. Physics and applications of laser diode chaos. Nat. Photonics 2015, 9, 151–162. [Google Scholar] [CrossRef] [Green Version]
  2. Soriano, M.C.; García-Ojalvo, J.; Mirasso, C.R.; Fischer, I. Complex photonics: Dynamics and applications of delay-coupled semiconductors lasers. Rev. Modern Phys. 2013, 85, 421–470. [Google Scholar] [CrossRef] [Green Version]
  3. Agrawal, G.P. Line Narrowing in a Single-Mode Injection Laser Due to External Optical Feedback. IEEE J. Quantum Electron. 1984, 20, 468–471. [Google Scholar] [CrossRef]
  4. Olesen, H.; Osmundsen, J.H.; Tromborg, B. Nonlinear dynamics and spectral behaviour for an external cavity laser. IEEE J. Quantum Electron. 1986, 22, 762–773. [Google Scholar] [CrossRef]
  5. Petermann, K. External optical feedback phenomena in semiconductor lasers. IEEE J. Sel. Top. Quantum Electron. 1995, 1, 480–489. [Google Scholar] [CrossRef]
  6. Patzak, E.; Olesen, H.; Sugimura, A.; Saito, S.; Mukai, T. Spectral linewidth reduction in semiconductor lasers by an external cavity with weak optical feedback. Electron. Lett. 1983, 19, 938–940. [Google Scholar] [CrossRef]
  7. Tkach, R.W.; Chraplyvy, A.R. Regimes of Feedback Effects in 1.5-µm Distributed Feedback Lasers. J. Lightw. Technol. 1986, 4, 1655–1661. [Google Scholar] [CrossRef]
  8. Lenstra, D.; Verbeek, B.; Boef, A.D. Coherence collapse in single-mode semiconductor lasers due to optical feedback. IEEE J. Quantum Electron. 1985, 21, 674–679. [Google Scholar] [CrossRef]
  9. Henry, C.; Kazarinov, R. Instability of semiconductor lasers due to optical feedback from distant reflectors. IEEE J. Quantum Electron. 1986, 22, 294–301. [Google Scholar] [CrossRef]
  10. Ebisawa, S.; Komatsu, S. Orbital Instability of Chaotic Laser Diode with Optical Injection and Electronically Applied Chaotic Signal. Photonics 2020, 7, 25. [Google Scholar] [CrossRef] [Green Version]
  11. Locquet, A. Routes to Chaos of a Semiconductor Laser Subjected to External Optical Feedback: A Review. Photonics 2020, 7, 22. [Google Scholar] [CrossRef] [Green Version]
  12. Gomez, S.; Huang, H.; Duan, J.; Combrié, S.; Shen, A.; Baili, G.; de Rossi, A.; Grillot, F. High coherence collapse of a hybrid III–V/Si semiconductor laser with a large quality factor. J. Phys. Photonics 2020, 2, 025005. [Google Scholar] [CrossRef]
  13. Zhao, D.; Andreou, S.; Yao, W.; Lenstra, D.; Williams, K.; Leijtens, X. Monolithically Integrated Multiwavelength Laser With Optical Feedback: Damped Relaxation Oscillation Dynamics and Narrowed Linewidth. IEEE Photonics J. 2018, 10, 6602108. [Google Scholar] [CrossRef]
  14. Aoyama, K.; Yokota, N.; Yasaka, H. Strategy of optical negative feedback for narrow linewidth semiconductor lasers. Opt. Express 2018, 26, 21159–21169. [Google Scholar] [CrossRef]
  15. Doerr, C.R.; Dupuis, N.; Zhang, L. Optical isolator using two tandem phase modulators. Opt. Lett. 2011, 36, 4293–4295. [Google Scholar] [CrossRef]
  16. Van Schaijk, T.T.M.; Lenstra, D.; Williams, K.A.; Bente, E.A.J.M. Model and experimental validation of a unidirectional phase modulator. Opt. Express 2018, 26, 32388–32403. [Google Scholar] [CrossRef] [PubMed]
  17. Septon, T.; Becker, A.; Gosh, S.; Shtendel, G.; Sichkovskyi, V.; Schnabel, F.; Sengül, A.; Bjelica, M.; Witzigmann, B.; Reithmaier, J.P.; et al. Large linewidth reduction in semiconductor lasers based on atom-like gain material. Optica 2019, 6, 1071–1077. [Google Scholar] [CrossRef]
  18. Duan, J.; Huang, H.; Dong, B.; Jung, D.; Norman, J.C.; Bowers, J.E.; Grillot, F. 1.3-µ m Reflection Insensitive InAs/GaAs Quantum Dot Lasers Directly Grown on Silicon. IEEE Photonics Technol. Lett. 2019, 31, 345–348. [Google Scholar] [CrossRef]
  19. Yamasaki, K.; Kanno, K.; Matsumoto, A.; Akahane, K.; Yamamoto, N.; Naruse, M.; Uchida, A. Fast dynamics of low-frequency fluctuations in a quantum-dot laser with optical feedback. Opt. Express 2021, 29, 17962–17975. [Google Scholar] [CrossRef] [PubMed]
  20. Huang, D.; Pintus, P.; Bowers, J.E. Towards heterogeneous integration of optical isolators and circulators with lasers on silicon. Opt. Mat. Express 2018, 8, 2471–2483. [Google Scholar] [CrossRef] [Green Version]
  21. Yan, W.; Yang, Y.; Liu, S.; Zhang, Y.; Xia, S.; Kang, T.; Yang, W.; Qin, J.; Deng, L.; Bi, L. Waveguide-integrated high-performance magneto-optical isolators and circulators on silicon nitride platforms. Optica 2020, 7, 1555–1562. [Google Scholar] [CrossRef]
  22. Lenstra, D.; van Schaijk, T.T.M.; Williams, K.A. Toward a Feedback-Insensitive Semiconductor Laser. IEEE J. Sel. Top. Quantum Electron. 2019, 25, 1–13. [Google Scholar] [CrossRef]
  23. Khoder, M.; der Sande, G.V.; Danckaert, J.; Verschaffelt, G. Effect of External Optical Feedback on Tunable Micro-Ring Lasers Using On-Chip Filtered Feedback. IEEE Photonics Technol. Lett. 2016, 28, 959–962. [Google Scholar] [CrossRef]
  24. Komljenovic, T.; Liang, L.; Chao, R.L.; Hulme, J.; Srinivasan, S.; Davenport, M.; Bowers, J.E. Widely-Tunable Ring-Resonator Semiconductor Lasers. Appl. Sci. 2017, 7, 732. [Google Scholar] [CrossRef]
  25. Kasai, K.; Nakazawa, M.; Ishikawa, M.; Ishii, H. 8 kHz linewidth, 50 mW output, full C-band wavelength tunable DFB LD array with self-optical feedback. Opt. Express 2018, 26, 5675–5685. [Google Scholar] [CrossRef]
  26. Morton, P.A.; Morton, M.J. High-Power, Ultra-Low Noise Hybrid Lasers for Microwave Photonics and Optical Sensing. J. Lightw. Technol. 2018, 36, 5048–5057. [Google Scholar] [CrossRef]
  27. Xiang, C.; Morton, P.A.; Bowers, J.E. Ultra-narrow linewidth laser based on a semiconductor gain chip and extended Si3N4 Bragg grating. Opt. Lett. 2019, 44, 3825–3828. [Google Scholar] [CrossRef]
  28. Liu, Y.; Ohtsubo, J. Dynamics and chaos stabilization of semiconductor lasers with optical feedback from an interferometer. IEEE J. Quantum. Elect. 1997, 33, 1163–1169. [Google Scholar] [CrossRef]
  29. Mezzapesa, F.P.; Columbo, L.L.; Dabbicco, M.; Brambilla, M.; Scamarcio, G. QCL-based nonlinear sensing of independent targets dynamics. Opt. Express 2014, 22, 5867. [Google Scholar] [CrossRef] [PubMed]
  30. Zhu, W.; Chen, Q.; Wang, Y.; Luo, H.; Wu, H.; Ma, B. Improvement on vibration measurement performance of laser self-mixing interference by using a pre-feedback mirror. Opt. Laser Eng. 2018, 105, 150–158. [Google Scholar] [CrossRef]
  31. Ruan, Y.; Liu, B.; Yu, Y.; Xi, J.; Guo, Q.; Tong, J. High sensitive sensing by a laser diode with dual optical feedback operating at period-one oscillation. Appl. Phys. Lett. 2019, 115, 011102. [Google Scholar] [CrossRef]
  32. Seimetz, M. Laser Linewidth Limitations for Optical Systems with High-Order Modulation Employing Feed Forward Digital Carrier Phase Estimation; Paper OTuM2; OFC: Auckland, New Zealand, 2008. [Google Scholar]
  33. Liang, W.; Ilchenko, V.S.; Eliyahu, D.; Dale, E.; Savchenkov, A.A.; Seidel, D.; Matsko, A.B.; Maleki, L. Compact stabilized semiconductor laser for frequency metrology. Appl. Opt. 2015, 54, 3353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Barton, J.S.; Skogen, E.J.; Mašanović, M.L.; Denbaars, S.P.; Coldren, L.A. A widely tunable high-speed transmitter using an integrated SGDBR laser-semiconductor optical amplifier and Mach–Zehnder modulator. IEEE J. Sel. Top. Quantum Electron. 2003, 9, 1113–1117. [Google Scholar] [CrossRef]
  35. Huang, D.; Tran, M.A.; Guo, J.; Peters, J.; Komljenovic, T.; Malik, A.; Morton, P.A.; Bowers, J.E. High-power sub-kHz linewidth lasers fully integrated on silicon. Optica 2019, 5, 745–752. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, H.; Kim, D.; Harfouche, M.; Santis, C.T.; Satyan, N.; Rakuljic, G.; Yariv, A. Narrow-Linewidth Oxide-Confined Heterogeneously Integrated Si III-V Semiconductor Lasers. IEEE Photonics Technol. Lett. 2017, 29, 2199–2202. [Google Scholar] [CrossRef]
  37. Lin, Y.; Browning, C.; Timens, R.B.; Geuzebroek, D.H.; Roeloffzen, C.G.H.; Hoekman, M.; Geskus, D.; Oldenbeuving, R.M.; Heideman, R.G.; Fan, Y.; et al. Characterization of Hybrid InP-TriPleX Photonic Integrated Tunable Lasers Based on Silicon Nitride (Si3N4/SiO2) Microring Resonators for Optical Coherent System. IEEE Photonics J. 2018, 10, 1400108. [Google Scholar] [CrossRef]
  38. Li, B.; Jin, W.; Wu, L.; Chang, L.; Wang, H.; Shen, B.; Yuan, Z.; Feshali, A.; Paniccia, M.; Vahala, K.J.; et al. Reaching fiber-laser coherence in integrated photonics. Opt. Lett. 2021, 46, 5201–5204. [Google Scholar] [CrossRef]
  39. Lang, R.; Kobayashi, K. External Optical Feedback Effects on Semiconductor Injection Laser Properties. IEEE J. Quantum Electron. 1980, 16, 347–355. [Google Scholar] [CrossRef]
  40. Smit, M.; Leijtens, X.; Ambrosius, H.; Bente, E.; van der Tol, J.; Smalbrugge, B.; de Vries, T.; Geluk, E.J.; Bolk, J.; van Veldhoven, R.; et al. An introduction to InP-based generic integration technology. Semicond. Sci. Technol. 2014, 29, 083001. [Google Scholar] [CrossRef]
  41. Augustin, L.M.; Santos, R.; den Haan, E.; Kleijn, S.; Thijs, P.J.A.; Latkowski, S.; Zhao, D.; Yao, W.; Bolk, J.; Ambrosius, H.; et al. InP-Based Generic Foundry Platform for Photonic Integrated Circuits. IEEE J. Sel. Top. Quantum Electron. 2018, 24, 6100210. [Google Scholar] [CrossRef]
  42. Henry, C.H. Theory of the Linewidth of Semiconductor Lasers. IEEE J. Quantum Electron. 1982, 18, 259–264. [Google Scholar] [CrossRef]
  43. Osinski, M.; Buus, J. Linewidth broadening factor in semiconductor lasers–An overview. IEEE J. Quantum Electron. 1987, 23, 9–29. [Google Scholar] [CrossRef]
  44. Lax, M.; Noise, V.C. Noise in Self-Sustained Oscillators. Phys. Rev. 1967, 160, 290–307. [Google Scholar] [CrossRef]
  45. Rowe, H.E. Signals and Noise in Communication Systems, 1st ed.; D. Van Nostrand Company. Inc.: Princeton, NJ, USA, 1965. [Google Scholar]
  46. Spano, P.; Piazzolla, S.; Tamburrini, M. Theory of noise in semiconductor lasers in the presence of optical feedback. IEEE J. Quantum Electron. 1984, 20, 350–357. [Google Scholar] [CrossRef]
  47. Schawlow, A.L.; Townes, C.H. Infrared and Optical Masers. Phys. Rev. 1958, 112, 1940–1949. [Google Scholar] [CrossRef] [Green Version]
  48. Yu, Y.; Giuliani, G.; Donati, S. Measurement of the Linewidth Enhancement Factor of Semiconductor Lasers Based on the Optical Feedback Self-Mixing Effect. IEEE Photonics Technol. Lett. 2004, 16, 990–992. [Google Scholar] [CrossRef]
  49. Schunk, N.; Petermann, K. Numerical analysis of the feedback regimes for a single-mode semiconductor laser with external feedback. IEEE J. Quantum Electron. 1988, 24, 1242–1247. [Google Scholar] [CrossRef]
  50. Panozzo, M.; Quintero-Quiroz, C.; Tiana-Alsina, J.; Torrent, M.C.; Masoller, C. Experimental characterization of the transition to coherence collapse in a semiconductor laser with optical feedback. Chaos 2017, 27, 114315. [Google Scholar] [CrossRef] [Green Version]
  51. Happach, M.; Felipe, D.D.; Friedhoff, V.N.; Kresse, M.; Irmscher, G.; Kleinert, M.; Zawadzki, C.; Rehbein, W.; Brinker, W.; Mohrle, M.; et al. Influence of Integrated Optical Feedback on Tunable Lasers. IEEE J. Quantum Electron. 2020, 56. [Google Scholar] [CrossRef]
  52. Petermann, K. Laser Diode Modulation and Noise (Advances in Optoelectronics 3); Kluwer Academic Publishers: Dordrecht, The Netherlands, 1988; Chapter 2. [Google Scholar]
Figure 1. Schematic of a laser cavity affected by external optical feedback from both sides.
Figure 1. Schematic of a laser cavity affected by external optical feedback from both sides.
Photonics 09 00043 g001
Figure 2. Numerical solutions for C = 0.5 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure 2. Numerical solutions for C = 0.5 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g002
Figure 3. Numerical solutions for C = 1 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure 3. Numerical solutions for C = 1 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g003
Figure 4. Numerical solutions for C = 1.3 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Figure 4. Numerical solutions for C = 1.3 under conditions (8) and (11) for different Δ ϕ m . Double-feedback case shown in blue, single feedback case shown in orange. Column (A) shows the lasing frequency shift results. Column (B) shows the threshold gain shift. Column (C) shows the intrinsic linewidth variations.
Photonics 09 00043 g004
Table 1. Summary of values used in numerical solutions of Equations (5) and (26).
Table 1. Summary of values used in numerical solutions of Equations (5) and (26).
ParameterValue
ω ref t ext (0,2 π )
α H 3
0.1581 (case 1)
κ 2 t 2 0.3162 (case 2)
0.4111 (case 3)
ϕ m 2 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Far Brusatori, M.; Volet, N. Dynamics of Semiconductor Lasers under External Optical Feedback from Both Sides of the Laser Cavity. Photonics 2022, 9, 43. https://doi.org/10.3390/photonics9010043

AMA Style

Far Brusatori M, Volet N. Dynamics of Semiconductor Lasers under External Optical Feedback from Both Sides of the Laser Cavity. Photonics. 2022; 9(1):43. https://doi.org/10.3390/photonics9010043

Chicago/Turabian Style

Far Brusatori, Mónica, and Nicolas Volet. 2022. "Dynamics of Semiconductor Lasers under External Optical Feedback from Both Sides of the Laser Cavity" Photonics 9, no. 1: 43. https://doi.org/10.3390/photonics9010043

APA Style

Far Brusatori, M., & Volet, N. (2022). Dynamics of Semiconductor Lasers under External Optical Feedback from Both Sides of the Laser Cavity. Photonics, 9(1), 43. https://doi.org/10.3390/photonics9010043

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