Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Basic phase-locking, noise, and modulation properties of optically mutual-injected terahertz quantum cascade lasers

Open Access Open Access

Abstract

The phase-locking, noise, and modulation properties of two face-to-face optically mutual-injected terahertz (THz) quantum cascade lasers (QCLs) are analyzed theoretically. In the phase-locking range, the two THz QCLs are in stationary states working at the same frequency. Outside the phase-locking range, the amplitude and the instantaneous frequency of the optical field oscillate with time, and the power spectrum shows a series of discrete peaks. For strong mutual injection, the optical field of the THz QCL array also exhibits oscillatory behavior. Coherent collapse or chaotic behavior is not observed within the range of the parameters used in this simulation. The spontaneous emission noise of phase-locked THz QCLs is higher than that of THz QCLs at free-running operation, and mutual injection may introduce additional modulation peaks in the noise spectrum. The modulation response of the mutual-injected THz QCLs to an individual modulation is investigated. It is found that the modulation bandwidth and the phase difference are significantly dependent on the modulation parameters. These results are helpful for further understanding the nonlinear dynamic behaviors of THz QCLs under optical injection and provide theoretical support for the development of THz QCL phase-locked arrays.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

In the last two decades, quantum cascade lasers (QCLs) have become important sources in the mid-infrared and terahertz spectral range [1, 2]. QCL is a new type of unipolar laser based on intersubband transition and resonant tunneling in a repeated stack of semiconductor multiple quantum well heterostructures. Compact size, widely tunable frequency range, and high coherence make THz QCLs useful in imaging [3, 4], material analysis [4, 5], gas monitoring [6], and THz-communication [7].

In applications, high-performance THz QCLs are highly desired. High power THz QCLs have been achieved through design optimization and fabrication improvements [8–10]. The output power of pulsed devices can reach Watt level [8], and the continuous wave operation can reach hundreds of milliwatts [9]. However, these devices usually work with multi-modes, which may limit their application in some areas where spectral purity is highly required. The power of single-mode THz QCLs is still low, and the beam quality of facet-emitting THz QCLs is bad. In order to improve the single-mode output power and beam quality of THz QCLs, many efforts have been made, such as using graded photonic heterostructures [10], master-oscillator power-amplifier [11], tapered waveguide [12], and symmetric wafer bonded active regions [13].

It is an important method by using optically mutual-injected phase-locked laser array to improve the output power and beam quality of single-mode devices. The dynamical behavior of phase-locked lasers has long been concerned due to the requirement of phase locking between independent laser sources, and the diode laser arrays have been extensively studied [14, 15]. The symmetry-breaking [16, 17], stability condition [18, 19], and chaos synchronization [20, 21] in laser arrays have been investigated in detail. It has been found that the phase-locking properties of mutual-injected lasers are strongly dependent on the delay time, the coupling strength and the detuning frequency [22–24].

For THz-QCLs, great efforts have been made in developing Y-coupled devices [25], global antenna mutual coupling [26], phase-locked surface-emitting arrays [27], external cavity coherent quantum cascade laser array [28], and ring resonator arrays [29], et al. However, compared with the developed mutual-injected diode laser arrays, the studies on mutual-injected THz QCLs are still limited. The optical mutual injection could be regarded as the basic mechanism for face-to-face injected, lateral coupled or even some other kinds of phase-locked QCL arrays. Ultra-short carrier lifetime and small linewidth enhancement factor may make THz QCL arrays different from traditional semiconductor lasers. The lack of understanding for the basic dynamics of optically mutual-injected QCLs is an obstacle for developing corresponding devices. In this work, we will investigate the dynamics of two face-to-face mutual-injected THz QCLs to reveal the basic phase-locking property of THz QCL arrays.

Noise is a key factor affecting the performance of THz QCLs in potential applications [30–32]. Noises in THz QCLs come from spontaneous emission, non-radiative losses of carriers, interactions between photons and carriers, fluctuations of carriers, and so on. The spontaneous emission is an intrinsic noise source for lasers. The spontaneous emission noise of free-running and master-slave injection locking mid-infrared QCL has been investigated [33–35]. For mutual-injected THz QCLs, an important problem is how mutual injection affects the noise level. So we will characterize the spontaneous emission noise of mutual-injected THz QCLs.

The frequency response and phase characteristics of THz QCLs under modulation are closely related to applications in high-speed data transmission. Although mid-infrared wireless communication has been implemented by employing the mid-infrared QCLs as transmitter [36], it is still a long-standing challenge to realize terahertz wireless transmission. At present, the research on modulation properties of THz QCLs has attracted much attention [37–39]. For the mutual-injected THz QCL array, an interesting problem is whether the whole arrays can respond to an individual modulation. We will discuss this problem and fundamental modulation properties for THz QCL arrays.

The paper is organized as follows: a theoretical model of optically mutual-injected THz QCLs based on the rate equations is described in Section 2. In Section 3, we firstly study the phase locking boundary and basic working properties of optically mutual-injected THz QCLs within and out of the phase locking range. Then we investigate the spontaneous emission noise of the mutual-injected THz QCLs under different injection current and coupling strength. Finally, we perform a small signal analysis of the rate equations to discuss the frequency response and phase difference of mutual-injected THz QCLs, followed with a conclusion in Section 4.

2. Theoretical model of optically mutual-injected terahertz quantum cascade lasers

The basic scheme of two face-to-face optically mutual-injected THz QCLs system is shown in the inset of Fig. 1. We assume that the two QCLs are single-mode ones and have identical parameters but differ in free-running frequencies. The rate equations describing the dynamics of the system are written as:

dεA(t)dt=iΩA0εA(t)+12(1+iα){Zg[NA(t)NT]1τp}εA(t)+κcτleiCpεB(tτext),
dεB(t)dt=iΩB0εB(t)+12(1+iα){Zg[NB(t)NT]1τp}εB(t)+κcτleiCpεA(tτext),
dNA(t)dt=ηIin,AqγeNA(t)g[NA(t)NT]|εA(t)|2,
dNB(t)dt=ηIin,BqγeNB(t)g[NB(t)NT]|εB(t)|2,
where εA(B) is the complex electric field of THz QCL A (B). NA(B) is the carrier number of one single period of active region in THz QCL A (B). ΩA0(B0) is the angular frequency of the free-running THz QCL A(B). ΔΩ = ΩB0 − ΩA0 is the detuning frequency. Z is the number of stages of the THz QCLs. g is the differential gain coefficient of the active region. NT is the carrier number of one single active region at transparency and τp is the lifetime of photon in the cavity. κc and Cp represent the coupling coefficient and coupling phase of the mutual injection, respectively. κc=ζ(1R)/R where R is the reflectivity of the emitting facet of the laser, ζ is the fraction of the injected light coupled into the lasing mode. In general, Cp does not qualitatively influence the results [40]. So in the latter calculation Cp is taken as 0.1π for simplification. It is worthwhile to note that we ignore the multiple longitudinal-mode effects caused by the cavity length L in our simulations. Then the cavity length only affects the round-trip time τl = 2nef f L/c, where nef f is the effective refractive index of the active region, and c is the speed of light in vacuum. The carrier decay rate is expressed as γe = 1/τe where τe is the lifetime of carrier. α is the linewidth enhancement factor. α is taken as 0.1 in simulations, which is a typical value for QCLs[41, 42].

 figure: Fig. 1

Fig. 1 The phase-locking range of face-to-face optically mutual-injected THz QCLs. The shadow region is the phase-locked range of the mutual-injected THz QCLs. Inset: the schematic of two THz QCLs with mutual injection. The red dashed line indicates the limitation of κ due to stability requirements.

Download Full Size | PDF

The delay time is expressed as τext = d/c, where d is the distance between the light-emitting facets of laser A and B. The value of d significantly affects the dynamics of laser array. It has been shown that for a very long d, coupled lasers exhibit periodic modulation [40] and low-frequency fluctuation [43]. In our simulation, d is fixed at 0.1m, which is small enough to avoid the complicated dynamical behaviors. This value is typical for an experimental setup of mutual-injected lasers on separated cryostats. For monolithic QCL arrays, the value of d may be further smaller. But even a smaller d will not change the dynamics of mutual-injected QCLs qualitatively.

Iin, A and Iin,B are the corresponding injection currents. η denotes the current injection efficiency into the lasing levels. q is the electronic charge. In our model, the processes of carrier’s injection are not discussed in detail. The phenomenon that the injection efficiency η decreases with the increase of injection current in realistic QCL is not considered here. so η is kept asa constant in the simulation. In general, the injection current can be expressed as (1 + a)Ith , where a is the pump parameter [33]. So in the later simulations, Ith is used as the unit of injection current directly.

For the convenience of numerical simulation, the use of dimensionless parameters is necessary. Denoting the instantaneous frequencies of the two THz QCLs by ΩA and ΩB, the slowly-varying complex electric fields can be expressed as EA(t) = |εA| exp[−iA − ΩA0)t] and EB(t) = |εB | exp[−iB − ΩB0)t]. We introduce the dimensionless time s and normalize the time to the photon lifetime by s = t/τp. Then the variables and parameters in Eqs. (1)(4) are nondimensionalized in accordance with the following rules:

EA(s)=εAexpiΩA0tEc,nA(s)=NANc0Nc,
EB(s)=εBexpiΩB0tEc,nB(s)=NBNc0Nc,
where Ec=γc/g is the normalization factor of optical field. Nc = γ/Zg is the normalization factor of carrier number where γ = 1/τp is the decay rate of phonon. Nc0 = NT + γ/Zg is a constant shift to ensure n = 0 at the laser threshold, i.e., nA and nB are the excess carrier number of THz QCLs A and B, respectively. The dimensionless rate equations for EA, EB and the excess carrier number nA, nB can be written as:
dEA(s)ds=12(1+iα)nA(s)EA(s)+κeiCpeiωB0θEB(sθ)eiΔωs,
dEB(s)ds=12(1+iα)nB(s)EB(s)+κeiCpeiωA0θEA(sθ)eiΔωs,
TdnA(s)ds=pAnA(s)[1+nA(s)]|EA(s)|2,
TdnB(s)ds=pBnB(s)[1+nB(s)]|EB(s)|2.

Where Tτe/τp = γ/γe is defined as the ratio of carrier lifetime and photon lifetime. κ = κc(τp/τl) is defined as the coupling strength. It should be noted that κ is controllable because ζ and R can be changed by light path and experimental design. κ is quite similar to the feedback parameter C in Lang-Kobayashi equation [44] of self-mixing [45]. The coupling strength κ may be also calibrated in the similar manner of the calibration of feedback parameter C in self-mixing experiments [46]. So in this paper, κ is directly used as a simulation parameter for simplicity. θτext /τp is the dimensionless external cavity one-trip time. pA(B) = (Zg/γ)(ηIin, A(B)/eNT)− 1 is the pump rate of THz QCL A (B). ωA0B0, and Δω are the dimensionless parameters of ΩA0, ΩB0 and ΔΩ respectively. The meaning of all parameters and typical values for the simulation of THz QCLs are listed in Table 1.

Tables Icon

Table 1. Parameters values used in the simulations.

3. Results and discussion

3.1. Stability of optically mutual-injected terahertz quantum cascade lasers

In order to study the stability properties of the system, we first determine the phase-locking range by analyzing the stationary solutions of Eqs. (7)(10). When the QCLs are phase-locked, their fields and excess carrier numbers are assumed to have the form as follows:

EAL=FALei(ωLωA0)s+iϕL,nAL=pFAL2FAL2+1,
EBL=FBLei(ωLωB0)s,nBL=pFBL2FBL2+1,
where ωL is the locked frequency of the coupled laser system, and ϕL is the locked phase of the field EAL (s) relative to EBL (s). FAL(BL) and nAL(BL) are the steady-state solutions of the corresponding amplitude of electric field and excess carrier number of THz QCL A(B). Iin, A = Iin,B = Iin and pA = pB = p are assumed in order to simplify the analysis. By inserting Eqs. (11)(12) into Eqs. (7)(10), we have:
FAL=pnAL1+nAL,nAL=2κFBLFALcos(ωLθ+Cp+ϕL),
FBL=pnBL1+nBL,nBL=2κFALFBLcos(ωLθ+CpϕL),
ωLΩA0=κFBLFAL1+α2sin(ωLθ+Cp+arctanα+ϕL),
ωLΩB0=κFALFBL1+α2sin(ωLθ+Cp+arctanαϕL).

Although those coupled equations cannot be solved analytically, they provide some useful information about the condition of the phase-locking.

We further assume that the amplitude of the electric field satisfies such a relationship that FAL /FBL ≈ 1. This approximation is reasonable when the two phase-locked THz QCLs are identical except their free-running frequencies. Based on the Eqs. (15) and (16), the detuning frequency can be expressed as

Δω=2κ1+α2cos(ωLθ+Cp+arctanα)sin(ϕL).

According to the boundedness of the trigonometric function, the boundary of the phase-locked detuning frequency can be obtained by:

|Δω|max=2κ1+α2.

The phase-locking range is shown as the shaded area in Fig. 1. With the increase of the coupling strength κ, the model allows larger frequency detuning for phase locking.

The working properties of mutual-injected THz QCLs depend on the coupling strength κ and detuning frequency ΔΩ/2π. In order to reveal the dynamic properties of optically mutual-injected THz QCLs, we calculate the time evolution of |EA|, |EB |, instantaneous frequency, and power spectral density from Eqs. (7)(10) for different detuning frequency ΔΩ/2π and coupling strength κ. As shown in Figs. 2(a)2(c), the equations are solved numerically with ΔΩ/2π = 0.5GHz and κ = 9.87 × 10−3. The coupling strength is moderate and the frequency detuning is within the locking range. In this case for κ = 9.87 × 10−3, the maximum allowable locking range is 0.547 GHz based on Eq. (18). So the two THz QCLs with ΔΩ/2π = 0.5GHz are phase-locked. From the time evolution of optical fields in Fig. 2(a), it can be confirmed that the lasers are rapidly stabilized after a short period of interaction. The sharp peak at ΔΩ/2π = 0 in the power spectral density curve in Fig. 2(b) indicates that the DC component of the output power is dominant. The instantaneous frequencies of the two THz QCLs are equal and stable.

 figure: Fig. 2

Fig. 2 Time evolution of |EA|, |EB | (first column), the corresponding power spectral density (second column), and instantaneous frequency (third column) with different coupling strength κ and detuning frequency ΔΩ/2π. The coupling strength in the first three rows are set as κ = 9.87 × 10−3, which is the case of moderate coupling. (a)–(c) within the phase-locking regime, ΔΩ/2π = 0.5GHz, (d)–(f) out of the pahse-locking regime, ΔΩ/2π = 5GHz, (g)–(i) out of the pahse-locking regime, ΔΩ/2π = 0.55GHz. The fourth row is the case of strong coupling with κ = 0.247. (j)–(l) out of the pahse-locking regime, ΔΩ/2π = 0.2GHz. The effective injection current is 1.5Ith in all simulations.

Download Full Size | PDF

In Figs. 2(d)2(f), we keep the coupling strength κ unchanged and increase the detuning frequency ΔΩ/2π to 5.0GHz which is far greater than the locking boundary. It is found that the electric fields oscillate with time. Since |EA|≈|EB |, the power spectral density of the two THz QCLs is almost equal as shown in Fig. 2(d). In addition, the power spectral density has sidebands at n · ΔΩ/2π(n = 1, 2, 3…). It indicates that the coupling of THz QCLs is equivalent to mutual frequency modulation of frequency ΔΩ/2π if the detuning frequency is large enough. In this case, the instantaneous frequencies of the two lasers exhibit periodic oscillations near their free-running frequencies.

Then, as shown in Figs. 2(g)2(i), we still keep κ unchanged and set the detuning frequency ΔΩ/2π to 0.55GHz which is a little greater than the locking boundary determined by Eq. (18). The optical field cannot reach a stable state, but remains a bounded oscillation even after a long time evolution. In the frequency domain, the power spectral density of both THz QCLs is split into many separated lines. These peaks do not correspond to detuning frequencies n · ΔΩ/2π(n = 1, 2, 3…). The frequencies of the two lasers are still not locked. The instantaneous frequencies of the two THz QCLs show periodic oscillations.

Figures 2(j)2(l) show the strong coupling case with κ = 0.247 and ΔΩ/2π = 0.2GHz. The lasers are in the phase-locking range according to Eq. (18), but unlike the case in Figs. 2(a)2(c), the two lasers cannot achieve their steady states with such a large coupling strength. In time domain, the amplitude of optical field still has a very small oscillation with time, as shown in Fig. 2(j). Moreover, although the instantaneous frequencies of the two lasers are very close, they also behave as small oscillations and cannot be locked. Compared with Fig. 2(h), the power spectral density of the non-DC component is much lower, and the sidebands do not correspond to n · ΔΩ/2π(n = 1, 2, 3…). Conventional semiconductor lasers exhibit coherent collapse or chaotic behavior under strong light injection. Similar behavior is not found in the simulation of mutual-injected THz QCL array. By comparison, the stability of the THz QCL array under strong injection is better than that of diode laser arrays.

3.2. Spontaneous emission noise of optically mutual-injected terahertz quantum cascade lasers

After discussing the phase-locking property of mutual-injected THz QCLs, we continue to analyze the noise property of the model in the phase-locking range. An intrinsic source of noise in THz QCLs is the fluctuations in the complex field EA, EB caused by the spontaneous emission of photons. This kind of noise is very important for dynamics and can not be ignored in general. The equations of electric fields with noise are expressed as:

dEAds=12(1+iα)nAEA+κeiCpeiωB0θEB(sθ)eiΔωs+XA(s),
dEBds=12(1+iα)nBEB+κeiCpeiωA0θEA(sθ)eiΔωs+XB(s).

Where XA(B) denotes the spontaneous emission noise of THz QCL A which follows the Langevin force formulation. Time evolutions of the real and imaginary parts of XA(B) are independent random processes. XA and XB have zero mean < XA >=< XB >= 0 and the following correlations:

<XA(s)XA(s)¯>=β(nA+nA0)δ(ss),
<XB(s)XB(s)¯>=β(nB+nB0)δ(ss),
where β = 10−6 is spontaneous emission factor representing the fraction of spontaneous emission processes contribute to the lasing mode. nA0 = nB0 = gNT /γ is the dimensionless carrier number at threshold. Considering the symmetry of the detuning frequency in the model, we only calculate the relative intensity noise (RIN) of THz QCL A. PA is the normalized photon number so that PA = |EA|2. The RIN spectrum can be calculated by:
RINA(ω)=limT1T|0T[PA(s)<PA(s)>]eiωsds|2<PA(s)>2.

The RINs are obtained by averaging over 20 independent calculations in order to avoid fluctuation caused by randomness and limited time evolution.

In order to reveal the effect of mutual injection on the spontaneous emission noise of devices, we compare the RIN between the mutual-injected array and the free-running laser. The RIN of a free-running QCL can be solved analytically by a semiclassical noise model [33, 34]. The RIN of the single free-running THz QCL A is calculated and shown in Fig. 3. The RIN spectrum is flat at low frequencies and exhibits white noise, which is consistent with the results in [33]. For κ = 9.87 × 10−3 and ΔΩ/2π = 0.2GHz, the RIN of the optically mutual-injected THz QCLs is also given in Fig. 3. The RIN spectrum of phase-locked mutual-injected THz QCLs is still frequency-independent in low-frequency region. There is no characteristic peak corresponding to relaxation oscillations on the RIN spectrum. The overdamping of relaxation oscillations is attributed to the ultra-short electron lifetime in THz QCLs. Compared with the free-running case, the mutual injection scheme makes the noise level of the laser increase by 57dB/Hz.

 figure: Fig. 3

Fig. 3 The comparison of RINs for mutual-injected THz QCLs and the free-running one. Black line: moderate coupling within the pahse-locking regime with κ = 9.87 × 10−3 and ΔΩ/2π = 0.2GHz; Red line: moderate coupling out of the phase-locking regime with κ = 9.87 × 10−3 and ΔΩ/2π = 20GHz; Blue line: the RIN of a free-running THz QCL. And in all cases, Iin = 1.5Ith.

Download Full Size | PDF

According to the numerical calculation in Section 3.1, the two THz QCLs cannot be locked steadily when the detuning frequency is too large, such as κ = 9.87 × 10−3 and ΔΩ/2π = 20GHz. The feedback term κeiCpeiωB0θEB(sθ)eiΔΩs in Eq. (19) can be regarded as a modulation signal of 20GHz mixed in the noise XA(s). Then the coupled lasers are approximately free-running with the noise containing modulated signals. We calculate the RIN of the THz QCL A in this case, and compare it to the RIN with κ = 9.87 × 10−3 and ΔΩ/2π = 0.2GHz where the laser can be locked steadily. As shown in Fig. 3, in these two cases, the RINs are approximately equal below several GHz and show the characteristics of white noise. There are peaks at 20, 40, 60GHz in the case κ = 9.87 × 10−3 , ΔΩ/2π = 20GHz resulted from the frequency modulation.

The RINs of THz QCL A with coupling coefficients κ = 9.87 × 10−4 and κ = 9.87 × 10−2 are shown in Fig. 4(a). The RINs increase with coupling coefficient. Although increasing the coupling strength κ makes it easier for the array to achieve phase locking, the spontaneous emission noise increases.

 figure: Fig. 4

Fig. 4 RINs of two mutual-injected THz QCLs with different coupling strengths and injection currents. (a) κ = 9.87 × 10−4 , 9.87 × 10−2, with the effective injection current and detuning frequency set as Iin = 1.5Ith and ΔΩ/2π = 0.01GHz, respectively. (b) Iin = 1.5Ith, 2Ith, with the coupling strength and detuning frequency set as κ = 9.87 × 10−3, ΔΩ/2π = 0.2GHz, respectively.

Download Full Size | PDF

It is known that the spontaneous emission noises of THz QCLs are sensitive to current changes. Fig. 4(b) shows the RIN of THz QCL A with injection currents 1.5Ith and 2Ith with κ = 9.87 × 10−3 , ΔΩ/2π = 0.2GHz. It is found that an appropriate increase of current can reduce the spontaneous emission noise of phase-locked THz QCLs.

3.3. Modulation response characteristics of optically mutual-injected terahertz quantum cascade lasers

The modulation response of THz QCLs is an important aspect of laser performance. Under the condition of high-frequency modulation, ensuring high output signal amplitude is essential to high-speed data transmission. Theoretical and experimental researches on the modulation responses of free-running and injection-locked mid-infrared QCLs have been discussed in [47, 50]. The optically mutual-injected THz QCLs model is one of the most typical delay-coupled systems, and its frequency response characteristics are of great concern. In order to investigate the modulation response of mutual-injected THz QCLs, a small-signal analysis is carried out to Eqs. (7)(10). ΔFA, FB, ΔϕA, ΔϕB, ΔnA and ΔnB are small variances around the steady state values, i.e.,

{FA=ΔFA+FAL,FB=ΔFB+FBLϕA=ΔϕA+ϕAL,ϕB=ΔϕB+ϕBLnA=ΔnA+nAL,nB=ΔnB+nBL,
the rate equations for small signal modulation can be derived by substituting Eq. (24) into Eqs. (7)(10).
dΔFA(s)ds=12nALΔFA(s)+12ΔnA(s)FALκFBL[ΔϕA(s)ΔϕB(sθ)]sin(ωLθ+Cp+ϕL)+κΔFB(sθ)cos(ωLθ+Cp+ϕL),
dΔFB(s)ds=12nBLΔFB(s)+12ΔnB(s)FBLκFAL[ΔϕB(s)ΔϕA(sθ)]sin(ωLθ+Cp+ϕL)+κΔFA(sθ)cos(ωLθ+Cp+ϕL),
dΔϕA(s)dsFAL=(ωLωA0)ΔFA(s)+12αnALΔFA(s)κFBL[ΔϕA(s)ΔϕB(sθ)]cos(ωLθ+Cp+ϕL)κΔFB(sθ)sin(ωLθ+Cp+ϕL),
dΔϕB(s)dsFBL=(ωLωB0)ΔFB(s)+12αnBLΔFB(s)κFAL[ΔϕB(s)ΔϕA(sθ)]cos(ωLθ+CpϕL)κΔFA(sθ)sin(ωLθ+CpϕL),
dΔnA(s)dsTA=pA2FALΔFA(s)(nAL+1)ΔnA(s)(FAL2+1),
dΔnB(s)dsTB=pB2FBLΔFB(s)(nBL+1)ΔnB(s)(FBL2+1).

After taking the Laplace transform and rewriting them into a matrix form we can get:

(Q11Q12Q13Q14Q15Q16Q21Q22Q23Q24Q25Q26Q31Q32Q33Q34Q35Q36Q41Q42Q43Q44Q45Q46Q51Q52Q53Q54Q55Q56Q61Q62Q63Q64Q65Q66)(ΔFAΔFBΔϕAΔϕBΔnAΔnB)=Zgqγγe(0000ΔJAΔJB).
where the matrix elements are expressed as
{Q11=12nALλQ12=κeλθcosψALQ13=κFBLsin(ψAL)Q14=eλθκFBLsin(ψAL)Q15=12FALQ16=0Q21=κeλθcosψBLQ22=12nBLλQ23=eλθκFALsin(ψBL)Q24=κFALsin(ψBL)Q25=0Q26=12FBLQ31=12αnAL(ωLωA0)Q32=κeλθsin(ψAL)Q33=FBLλFALQ34=eλθκFBLQ35=12αFALQ36=0Q41=κeλθsin(ψBL)Q42=12αnBL(ωLωB0)Q43=eλθκFALQ44=κFALλFBLQ45=0Q46=12αFBLQ51=2FAL(nAL+1)Q52=0Q53=0Q54=0Q55=(FAL2+1)TλQ56=0Q61=0Q62=2FBL(nBL+1)Q63=0Q64=0Q65=0Q66=(FBL2+1)Tλ
with ψAL = ωL θ + Cp + ϕL and ψBL = ωL θ + CpϕL . Considering the symmetry of our model, we define the output signal as the electric field amplitude of THz QCL A. When only adding the modulation signal ΔJA to the laser A, the transfer function is written as
HAA(λ)=qγγeZgΔFAΔJA=Q1(000010).

Then the frequency response is expressed as |HAA()| and the phase difference between the output signal and the modulation signal is ∠HAA(). Due to the injection coupling of the two THz QCLs, there is frequency response in THz QCL A even when the modulation signal is only added to THz QCL B, the corresponding transfer function is

HAB(λ)=qγγeZgΔFAΔJB=Q1(000001).

And the corresponding modulation response and phase difference are |HAB()| and ∠HAB(), respectively.

The transfer function method is widely used to describe the frequency response of linear time-invariant systems. In our case, because the matrix Q contains delay term eλθ , the transfer function contains eλθ to e−3λθ terms which make the system complex. Fortunately, we do not particularly need the analytic zeros and poles of the transfer function. We are more concerned about the common nature of this kind of coupled system under modulation. So the numerical method is applied to calculate the modulation response of mutual-injected THz QCLs. The frequency responses of the lasers under different coupling strength κ and injection current Iin are calculated. The obtained Bode diagrams are presented in Fig. 5.

 figure: Fig. 5

Fig. 5 Within the phase-locking regime, the normalized frequency response |HAA| and |HAB | with different coupling strengths and injection currents. The frequency detunning ΔΩ/2π is set as 0.01GHz. (a) Iin = 1.5Ith, different coupling strength with κ = 9.87×10−4, 9.87 × 10−2. (b) κ = 9.87 × 10−3, different injection current Iin = 1.5Ith,2.0Ith . The −3 dB modulation renponse is indicated by the blue dashed line.

Download Full Size | PDF

The mutual-modulation frequency response |HAB()| is always less than the self-modulation frequency response |HAA()|. When the modulation is added to THz QCL A, we notice that ΔJA(s) linearly affects ΔnA(s) in Eq. (29) through pA = (Zg/γ)(ηIin, A/eNT)− 1. So ΔFA(s) is linearly dependent on ΔJA(s) through the term 1/2ΔnA(s)FAL in Eq. (25). According to the symmetry of the equations, the relationship between ΔFB(s) and ΔJB(s) is similar. However, when the modulation signal is added to THz QCL B, ΔJB(s) changes ΔFA(s) through the term κΔFB(sθ)cos(ωL θ + Cp + ϕL) in Eq. (25). Because κ is very small, it causes 1/2ΔnA(s)FALΔFB(sθ)cos(ωL θ + Cp + ϕL) in general. This is the main reason for |HAB()| < |HAA()|.

The 3dB bandwidth f3dB is defined as the frequency at which the modulation response |H(2πif3dB)| = 1/2|H(0)|. It is found that increasing the coupling strength and injection current can increase the modulation bandwidth of both modulation method. The modulation bandwidth of conventional semiconductor laser is restricted by the relaxation oscillation resulting from the interaction of the carriers and the photons. The normalized frequency response of the mutual-injected THz QCL array does not contain a characteristic peak corresponding to relaxation oscillation due to the ultra-short carrier lifetime of THz QCLs. Then for phase-locking arrays of THz QCLs, the absence of relaxation oscillations makes them very promising in high-frequency applications.

Since the model is a delayed coupling system, it is necessary to study the phase difference between input and output signals. ∠HAA() and ∠HAB() are shown in Figs. 6(a) and 6(b) respectively. As shown in Fig. 6, in both cases the low-frequency limit of phase difference is π, which comes from the negative coefficient −qγγe/Zg in Eqs. (33) and (34). It means that the phase of ΔFA is completely synchronized with ΔJA and ΔJB by low-frequency modulation. However, with the increase of modulation frequency, phase asynchrony is generated between response and modulation, and the phase difference decreases with the increase of the modulation frequency. It is found that ∠HAB() drops faster under the same delay time parameter. The results show that when d = 0.05m and the modulation frequency is 20 GHz, ∠HAA decreases by less than 30 degrees, but ∠HAB has fallen by more than 360 degrees.

 figure: Fig. 6

Fig. 6 The phase difference between input and output signals of the mutual-injected THz-QCLs under frequency modulation with d = 0.05m, 0.1m, 0.15m. (a) self-modulation ∠HAA(), (b) mutual-modulation ∠HAB(). The coupling strength, detuning frequency, and effective injection current are set as κ = 9.87 × 10−3, ΔΩ/2π = 0.01GHz, and 1.5Ith , respectively.

Download Full Size | PDF

Figure 6 also shows that ∠HAA() is independent of delay time, since 1/2ΔnA(s)FAL in Eq. (25) does not contain the delay time θ. On the contrary, ∠HAB() shows the dependence of delay time, and drops faster as the delay time increases. The reason behind this is that the instantaneous value of κΔFB(sθ)cos(ωL θ + Cp + ϕL) in Eq. (25) is related to the history of ΔFB, a phase difference is inevitable. It is found that ∠HAA() gradually tends to 180 degrees as the modulation frequency increases, which is similar to the properties of free-running and master-slave injection arrays. In contrast, ∠HAB() still depends on the modulation frequency.

4. Conclusion

In conclusion, the theoretical investigation of two face-to-face optically mutual-injected THz QCLs is carried out. It is shown that coupling coefficient κ and linewidth enhancement factor α determine the maximum allowable frequency detuning of phase-locking operation. According to the working parameters, the mutual-injected array can operate in the mode of stable phase locking or periodic oscillation. Particularly when the detuning frequency is far away from the locking range boundary, the effect of mutual injection on each laser is equivalent to a modulation of the detuning frequency. By studying the spontaneous emission noise in the model, it is found that the mutual-injected THz QCLs show a higher white noise level than free-running operation within several GHz frequency range. And a series of peaks corresponding to n · ΔΩ/2π(n = 1, 2, 3…) appear on the noise spectrum when the array is working far away from the locking region. By the method of small signal analysis, we get the modulation response property of the phase-locked THz QCLs. Increasing the injection current and coupling strength can increase the modulation bandwidth of the whole array. The modulation and response of the two THz QCLs at low-frequency limit are completely synchronous. With the increase of modulation frequency, phase difference is generated within the array. These results may help to further understand the nonlinear dynamic behavior and noise characteristics of THz QCLs under light injection, and provide theoretical support for the development of THz QCL phase-locked array.

Funding

National Natural Science Foundation of China (NSAF) Joint Fund (U1730246); National Natural Science Foundation of China (NSFC-RS) (F040302); Foundation of Director of IAPCM (ZYSZ1518-16).

References

1. J. Faist, F. Capasso, D. L. Sivco, C. Sirtori, A. L. Hutchinson, and A. Y. Cho, “Quantum cascade laser,” Science 264(5158), 553–556 (1994). [CrossRef]   [PubMed]  

2. M. Razeghi, Q. Y. Lu, N. Bandyopadhyay, W. Zhou, D. Heydari, Y. Bai, and S. Slivken, “Quantum cascade lasers: from tool to product,” Opt. Express 23(7), 8462–8475 (2015). [CrossRef]   [PubMed]  

3. A. W. M. Lee, Q. Qin, S. Kumar, B. S. Williams, Q. Hu, and J. L. Reno, “Real-time terahertz imaging over a standoff distance (> 25 meters),” Appl. Phys. Lett. 89(14), 141125 (2006). [CrossRef]  

4. A. D. Rakić, T. Taimre, K. Bertling, Y. L. Lim, P. Dean, D. Indjin, Z. Ikonić, P. Harrison, A. Valavanis, S. P. Khanna, M. Lachab, S. J. Wilson, E. H. Linfield, and A. G. Davies, “Swept-frequency feedback interferometry using terahertz frequency QCLs: a method for imaging and materials analysis,” Opt. Express 21(19), 22194–22205 (2013). [CrossRef]  

5. A. G. Davies, A. D. Burnett, W. Fan, E. H. Linfield, and J. E. Cunningham, “Terahertz spectroscopy of explosives and drugs,” Mater. Today 11(3), 18–26 (2008). [CrossRef]  

6. G. N. Rao and A. Karpf, “External cavity tunable quantum cascade lasers and their applications to trace gas monitoring,” Appl. Opt. 50(4), A100–A115 (2011). [CrossRef]   [PubMed]  

7. C. Wang, J. Cao, L. Gu, Q. Wu, and Z. Tan, “20 Mbps wireless communication demonstration using terahertz quantum devices,” Chin. Opt. Lett. 13(8), 81402 (2015). [CrossRef]  

8. L. Li, L. Chen, J. Zhu, J. Freeman, P. Dean, A. Valavanis, A. G. Davies, and E. H. Linfield, “Terahertz quantum cascade lasers with >1 W output powers,” Electron. Lett. 50(4), 309–311 (2014). [CrossRef]  

9. X. Wang, C. Shen, T. Jiang, Z. Zhan, Q. Deng, W. Li, W. Wu, N. Yang, W. Chu, and S. Duan, “High-power terahertz quantum cascade lasers with ~ 0.23 W in continuous wave mode,” AIP Adv. 6(7), 075210 (2016). [CrossRef]  

10. G. Xu, R. Colombelli, S. P. Khanna, A. Belarouci, X. Letartre, L. Li, E. H. Linfield, A. G. Davies, H. E. Beere, and D. A. Ritchie, “Efficient power extraction in surface-emitting semiconductor lasers using graded photonic heterostructures,” Nat. Commun. 3, 952 (2012). [CrossRef]   [PubMed]  

11. H. Zhu, F. Wang, Q. Yan, C. Yu, J. Chen, G. Xu, L. He, L. Li, L. Chen, A. G. Davies, E. H. Linfield, J. Hao, P. B. Vigneron, and R. Colombelli, “Terahertz master-oscillator power-amplifier quantum cascade lasers,” Appl. Phys. Lett. 109(23), 231105 (2016). [CrossRef]  

12. Y. F. Li, J. Wang, N. Yang, J. Liu, T. Wang, F. Liu, Z. Wang, W. Chu, and S. Duan, “The output power and beam divergence behaviors of tapered terahertz quantum cascade lasers,” Opt. Express 21(13), 15998–16006 (2013). [CrossRef]   [PubMed]  

13. M. Brandstetter, C. Deutsch, M. Krall, H. Detz, D. C. MacFarland, T. Zederbauer, A. M. Andrews, W. Schrenk, G. Strasser, and K. Unterrainer, “High power terahertz quantum cascade lasers with symmetric wafer bonded active regions,” Appl. Phys. Lett. 103(17), 171113 (2013). [CrossRef]  

14. E. Kapon, J. Katz, and A. Yariv, “Supermode analysis of phase-locked arrays of semiconductor lasers,” Opt. Lett. 9(4), 125–127 (1984). [CrossRef]   [PubMed]  

15. M. V. Romalis, “Narrowing of high power diode laser arrays using reflection feedback from an etalon,” Appl. Phys. Lett. 77(8), 1080–1081 (2000). [CrossRef]  

16. F. Rogister and J. García-Ojalvo, “Symmetry breaking and high-frequency periodic oscillations in mutually coupled laser diodes,” Opt. Lett. 28(14), 1176–1178 (2003). [CrossRef]   [PubMed]  

17. X. Qi, I. L. Garanovich, A. A. Sukhorukov, W. Krolikowski, A. Mitchell, G. Zhang, D. N. Neshev, and Y. S. Kivshar, “Polychromatic solitons and symmetry breaking in curved waveguide arrays,” Opt. Lett. 35(9), 1371–1373 (2010). [CrossRef]   [PubMed]  

18. P. K. Jakobsen, R. A. Indik, J. V. Moloney, A. C. Newell, H. G. Winful, and L. Raman, “Diode-laser array modes: discrete and continuous models and their stability,” J. Opt. Soc. Am. B 8(8), 1674–1680 (1991). [CrossRef]  

19. S. Zhao, Z. Wang, and Z. Wang, “Scalability and stability of a coherent fiber laser array using fiber couplers,” J. Opt. Soc. Am. B 27(5), 1131–1138 (2010). [CrossRef]  

20. L. Junges, A. Gavrielides, and J. A. C. Gallas, “Synchronization properties of two mutually delay-coupled semiconductor lasers,” J. Opt. Soc. Am. B 33(7), C65–C71 (2016). [CrossRef]  

21. N. Jiang, W. Pan, L. Yan, B. Luo, W. Zhang, S. Xiang, L. Yang, and D. Zheng, “Chaos synchronization and communication in mutually coupled semiconductor lasers driven by a third laser,” J. Lightwave Technol. 28(13), 1978–1986 (2010). [CrossRef]  

22. Y. Xie, Y. Li, J. Wang, N. Yang, W. Chu, and S. Duan, Chapter 5 in Quantum Cascade Lasers (IntechOpen, 2017).

23. D. Liu, C. Sun, B. Xiong, and Y. Luo, “Nonlinear dynamics in integrated coupled DFB lasers with ultra-short delay,” Opt. Express 22(5), 5614–5622 (2014). [CrossRef]   [PubMed]  

24. E. Wille, M. Peil, I. Fischer, and W. Elsäßer, “Dynamical scenarios of mutually delay-coupled semiconductor lasers in the short coupling regime,” Proc. SPIE 5452, 41–50 (2004). [CrossRef]  

25. O. P. Marshall, S. Chakraborty, M. Khairuzzaman, H. E. Beere, and D. A. Ritchie, “Reversible mode switching in Y-coupled terahertz lasers,” Appl. Phys. Lett. 96(102), 111105 (2013). [CrossRef]  

26. T.-Y. Kao, J. L. Reno, and Q. Hu, “Phase-locked laser arrays through global antenna mutual coupling,” Nat. Photonics 10(8), 541–546 (2016). [CrossRef]  

27. T. Y. Kao, Q. Hu, and J. L. Reno, “Phase-locked arrays of surface-emitting terahertz quantum-cascade lasers,” Appl. Phys. Lett. 96(10), 101106 (2010). [CrossRef]  

28. R. Vallon, B. Parvitte, L. Bizet, G. M. D. Naurois, B. Simozrag, G. Maisons, M. Carras, and V. Zeninari, “External cavity coherent quantum cascade laser array,” Infrared Phys. Technol. 76, 415–420(2016). [CrossRef]  

29. Y. Halioua, G. Xu, S. Moumdji, L. Li, J. Zhu, E. H. Linfield, A. G. Davies, H. E. Beere, D. A. Ritchie, and R. Colombelli, “Phase-locked arrays of surface-emitting graded-photonic-heterostructure terahertz semiconductor lasers,” Opt. Express 23(5), 6915–6923 (2015). [CrossRef]   [PubMed]  

30. M. Ravaro, V. Jagtap, C. Manquest, P. Gellie, G. Santarelli, C. Sirtori, S. P. Khanna, E. H. Linfield, and S. Barbieri, “Spectral properties of THz quantum-cascade lasers: frequency noise, phase-locking and absolute frequency measurement,” J. Infrared Milli. Terahz. Waves 34(5–6), 342–356 (2013). [CrossRef]  

31. M. Yamanishi, “Theory of intrinsic linewidth based on fluctuation-dissipation balance for thermal photons in THz quantum-cascade lasers,” Opt. Express 20(27), 28465–28478 (2012). [CrossRef]   [PubMed]  

32. T. Liu and Q. J. Wang, “Fundamental frequency noise and linewidth broadening caused by intrinsic temperature fluctuations in quantum cascade lasers,” Phys. Rev. B 84(12), 125322 (2011). [CrossRef]  

33. T. Gensty and W. Elsäßer, “Semiclassical model for the relative intensity noise of intersubband quantum cascade lasers,” Opt. Commun. 256(1), 171–183 (2005). [CrossRef]  

34. H. Simos, A. Bogris, D. Syvridis, and W. Elsäßer, “Intensity noise properties of mid-infrared injection locked quantum cascade lasers: I. modeling,” IEEE J. Quantum Electron. 50(2), 98–105 (2014). [CrossRef]  

35. M. M. Sheikhey, M. Goudarzi, R. Yadipour, and H. Baghban, “Analytical investigation of relative intensity noise properties of injection-locked mid-IR quantum cascade lasers,” J. Opt. Soc. Am. B 33(11), D57–D64 (2016). [CrossRef]  

36. F. Capasso, R. Paiella, R. Martini, R. Colombelli, C. Gmachl, T. L. Myers, M. S. Taubman, R. M. Williams, C. G. Bethea, K. Unterrainer, H. Y. Hwang, D. L. Sivco, A. Y. Cho, A. M. Sergent, H. C. Liu, and E. A. Whittaker, “Quantum cascade lasers: ultrahigh-speed operation, optical wireless communication, narrow linewidth, and far-infrared emission,” IEEE J. Quantum Electron. 38(6), 511–532 (2002). [CrossRef]  

37. F. Wang, X. G. Guo, C. Wang, and J. C. Cao, “Ultrafast population dynamics in electrically modulated terahertz quantum cascade lasers,” New J. Phys. 15(15), 075009 (2013). [CrossRef]  

38. W. Maineult, L. Ding, P. Gellie, P. Filloux, C. Sirtori, S. Barbieri, T. Akalin, J. F. Lampin, I. Sagnes, H. E. Beere, and D. A. Ritchie, “Microwave modulation of THz quantum cascade lasers: a transmission-line approach,” Appl. Phys. Lett. 96(2), 021108 (2010). [CrossRef]  

39. J. F. Webb, K. S. C. Yong, and M. K. Haldar, “Harmonic, intermodulation and cross-modulation distortion in directly modulated quantum cascade lasers,” IOP Conf. Ser.: Mater. Sci. Eng. 131(1), 012025 (2016). [CrossRef]  

40. E. A. Viktorov, A. M. Yacomotti, and P. Mandel, “Semiconductor lasers coupled face-to-face,” J. Opt. B: Quantum Semiclassical Opt. 6(2), L9–L12 (2004). [CrossRef]  

41. R. P. Green, J. Xu, L. Mahler, A. Tredicucci, F. Beltran, G. Giuliani, H. E. Beere, and D. A. Ritchie, “Linewidth enhancement factor of terahertz quantum cascade lasers,” Appl. Phys. Lett. 92(7), 071106 (2008). [CrossRef]  

42. C. Wang, F. Grillot, V. Kovanis, and J. Even, “Rate equation analysis of injection-locked quantum cascade lasers,” J. Appl. Phys. 113(6), 063104 (2013). [CrossRef]  

43. C. Q. Quiroz, J. T. Alsina, J. Roma, M. C. Torrent, and C. Masoller, “Quantitative identification of dynamical transitions in a semiconductor laser with optical feedback,” Sci Rep 6, 37510 (2016). [CrossRef]  

44. R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” IEEE J. Quantum Electron. 16(3), 347–355 (1980). [CrossRef]  

45. X. Qi, G. Agnew, I. Kundu, T. Taimre, Y. Lim, K. Bertling, P. Dean, A. Grier, A. Valavanis, E. Linfield, A. Giles Davies, D. Indjin, and A. Rakić, “Multi-spectral terahertz sensing: proposal for a coupled-cavity quantum cascade laser based optical feedback interferometer,” Opt. Express 20(9), 10153–10165 (2017). [CrossRef]  

46. J. Xi, Y. Yu, J. F. Chicharo, and T. Bosch, “Estimating the parameters of semiconductor lasers based on weak optical feedback self-mixing interferometry,” IEEE J. Quantum Electron. 41(8), 1058–1064 (2005). [CrossRef]  

47. B. Meng and Q. J. Wang, “Theoretical investigation of injection-locked high modulation bandwidth quantum cascade lasers,” Opt. Express 20(2), 1450–1464 (2012). [CrossRef]   [PubMed]  

48. S. Kumar, B. S. Williams, S. Kohen, Q. Hu, and J. L. Reno, “Continuous-wave operation of terahertz quantumcascade lasers above liquid-nitrogen temperature,” Appl. Phys. Lett. 84(14), 2494 (2004). [CrossRef]  

49. B. Williams, S. Kumar, Q. Hu, and J. Reno, “Operation of terahertz quantum-cascade lasers at 164 K in pulsed mode and at 117 K in continuous-wave mode,” Opt. Express 13(9), 3331–3339 (2005). [CrossRef]   [PubMed]  

50. R. Paiella, R. Martini, F. Capasso, C. Gmachl, H. Y. Hwang, D. L. Sivco, J. N. Baillargeon, A. Y. Cho, E. A. Whittaker, and H. C. Liu, “High-frequency modulation without the relaxation oscillation resonance in quantum cascade lasers,” Appl. Phys. Lett. 79(16), 2526–2528 (2001). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 The phase-locking range of face-to-face optically mutual-injected THz QCLs. The shadow region is the phase-locked range of the mutual-injected THz QCLs. Inset: the schematic of two THz QCLs with mutual injection. The red dashed line indicates the limitation of κ due to stability requirements.
Fig. 2
Fig. 2 Time evolution of |EA|, |EB | (first column), the corresponding power spectral density (second column), and instantaneous frequency (third column) with different coupling strength κ and detuning frequency ΔΩ/2π. The coupling strength in the first three rows are set as κ = 9.87 × 10−3, which is the case of moderate coupling. (a)–(c) within the phase-locking regime, ΔΩ/2π = 0.5GHz, (d)–(f) out of the pahse-locking regime, ΔΩ/2π = 5GHz, (g)–(i) out of the pahse-locking regime, ΔΩ/2π = 0.55GHz. The fourth row is the case of strong coupling with κ = 0.247. (j)–(l) out of the pahse-locking regime, ΔΩ/2π = 0.2GHz. The effective injection current is 1.5Ith in all simulations.
Fig. 3
Fig. 3 The comparison of RINs for mutual-injected THz QCLs and the free-running one. Black line: moderate coupling within the pahse-locking regime with κ = 9.87 × 10−3 and ΔΩ/2π = 0.2GHz; Red line: moderate coupling out of the phase-locking regime with κ = 9.87 × 10−3 and ΔΩ/2π = 20GHz; Blue line: the RIN of a free-running THz QCL. And in all cases, Iin = 1.5Ith.
Fig. 4
Fig. 4 RINs of two mutual-injected THz QCLs with different coupling strengths and injection currents. (a) κ = 9.87 × 10−4 , 9.87 × 10−2, with the effective injection current and detuning frequency set as Iin = 1.5Ith and ΔΩ/2π = 0.01GHz, respectively. (b) Iin = 1.5Ith, 2Ith, with the coupling strength and detuning frequency set as κ = 9.87 × 10−3, ΔΩ/2π = 0.2GHz, respectively.
Fig. 5
Fig. 5 Within the phase-locking regime, the normalized frequency response |HAA| and |HAB | with different coupling strengths and injection currents. The frequency detunning ΔΩ/2π is set as 0.01GHz. (a) Iin = 1.5Ith, different coupling strength with κ = 9.87×10−4, 9.87 × 10−2. (b) κ = 9.87 × 10−3, different injection current Iin = 1.5Ith,2.0Ith . The −3 dB modulation renponse is indicated by the blue dashed line.
Fig. 6
Fig. 6 The phase difference between input and output signals of the mutual-injected THz-QCLs under frequency modulation with d = 0.05m, 0.1m, 0.15m. (a) self-modulation ∠HAA(), (b) mutual-modulation ∠HAB(). The coupling strength, detuning frequency, and effective injection current are set as κ = 9.87 × 10−3, ΔΩ/2π = 0.01GHz, and 1.5Ith , respectively.

Tables (1)

Tables Icon

Table 1 Parameters values used in the simulations.

Equations (34)

Equations on this page are rendered with MathJax. Learn more.

d ε A ( t ) d t = i Ω A 0 ε A ( t ) + 1 2 ( 1 + i α ) { Z g [ N A ( t ) N T ] 1 τ p } ε A ( t ) + κ c τ l e i C p ε B ( t τ e x t ) ,
d ε B ( t ) d t = i Ω B 0 ε B ( t ) + 1 2 ( 1 + i α ) { Z g [ N B ( t ) N T ] 1 τ p } ε B ( t ) + κ c τ l e i C p ε A ( t τ e x t ) ,
d N A ( t ) d t = η I i n , A q γ e N A ( t ) g [ N A ( t ) N T ] | ε A ( t ) | 2 ,
d N B ( t ) d t = η I i n , B q γ e N B ( t ) g [ N B ( t ) N T ] | ε B ( t ) | 2 ,
E A ( s ) = ε A exp i Ω A 0 t E c , n A ( s ) = N A N c 0 N c ,
E B ( s ) = ε B exp i Ω B 0 t E c , n B ( s ) = N B N c 0 N c ,
d E A ( s ) d s = 1 2 ( 1 + i α ) n A ( s ) E A ( s ) + κ e i C p e i ω B 0 θ E B ( s θ ) e i Δ ω s ,
d E B ( s ) d s = 1 2 ( 1 + i α ) n B ( s ) E B ( s ) + κ e i C p e i ω A 0 θ E A ( s θ ) e i Δ ω s ,
T d n A ( s ) d s = p A n A ( s ) [ 1 + n A ( s ) ] | E A ( s ) | 2 ,
T d n B ( s ) d s = p B n B ( s ) [ 1 + n B ( s ) ] | E B ( s ) | 2 .
E A L = F A L e i ( ω L ω A 0 ) s + i ϕ L , n A L = p F A L 2 F A L 2 + 1 ,
E B L = F B L e i ( ω L ω B 0 ) s , n B L = p F B L 2 F B L 2 + 1 ,
F A L = p n A L 1 + n A L , n A L = 2 κ F B L F A L cos ( ω L θ + C p + ϕ L ) ,
F B L = p n B L 1 + n B L , n B L = 2 κ F A L F B L cos ( ω L θ + C p ϕ L ) ,
ω L Ω A 0 = κ F B L F A L 1 + α 2 sin ( ω L θ + C p + arctan α + ϕ L ) ,
ω L Ω B 0 = κ F A L F B L 1 + α 2 sin ( ω L θ + C p + arctan α ϕ L ) .
Δ ω = 2 κ 1 + α 2 cos ( ω L θ + C p + arctan α ) sin ( ϕ L ) .
| Δ ω | max = 2 κ 1 + α 2 .
d E A d s = 1 2 ( 1 + i α ) n A E A + κ e i C p e i ω B 0 θ E B ( s θ ) e i Δ ω s + X A ( s ) ,
d E B d s = 1 2 ( 1 + i α ) n B E B + κ e i C p e i ω A 0 θ E A ( s θ ) e i Δ ω s + X B ( s ) .
< X A ( s ) X A ( s ) ¯ > = β ( n A + n A 0 ) δ ( s s ) ,
< X B ( s ) X B ( s ) ¯ > = β ( n B + n B 0 ) δ ( s s ) ,
R I N A ( ω ) = lim T 1 T | 0 T [ P A ( s ) < P A ( s ) > ] e i ω s d s | 2 < P A ( s ) > 2 .
{ F A = Δ F A + F A L , F B = Δ F B + F B L ϕ A = Δ ϕ A + ϕ A L , ϕ B = Δ ϕ B + ϕ B L n A = Δ n A + n A L , n B = Δ n B + n B L ,
d Δ F A ( s ) d s = 1 2 n A L Δ F A ( s ) + 1 2 Δ n A ( s ) F A L κ F B L [ Δ ϕ A ( s ) Δ ϕ B ( s θ ) ] sin ( ω L θ + C p + ϕ L ) + κ Δ F B ( s θ ) cos ( ω L θ + C p + ϕ L ) ,
d Δ F B ( s ) d s = 1 2 n B L Δ F B ( s ) + 1 2 Δ n B ( s ) F B L κ F A L [ Δ ϕ B ( s ) Δ ϕ A ( s θ ) ] sin ( ω L θ + C p + ϕ L ) + κ Δ F A ( s θ ) cos ( ω L θ + C p + ϕ L ) ,
d Δ ϕ A ( s ) d s F A L = ( ω L ω A 0 ) Δ F A ( s ) + 1 2 α n A L Δ F A ( s ) κ F B L [ Δ ϕ A ( s ) Δ ϕ B ( s θ ) ] cos ( ω L θ + C p + ϕ L ) κ Δ F B ( s θ ) sin ( ω L θ + C p + ϕ L ) ,
d Δ ϕ B ( s ) d s F B L = ( ω L ω B 0 ) Δ F B ( s ) + 1 2 α n B L Δ F B ( s ) κ F A L [ Δ ϕ B ( s ) Δ ϕ A ( s θ ) ] cos ( ω L θ + C p ϕ L ) κ Δ F A ( s θ ) sin ( ω L θ + C p ϕ L ) ,
d Δ n A ( s ) d s T A = p A 2 F A L Δ F A ( s ) ( n A L + 1 ) Δ n A ( s ) ( F A L 2 + 1 ) ,
d Δ n B ( s ) d s T B = p B 2 F B L Δ F B ( s ) ( n B L + 1 ) Δ n B ( s ) ( F B L 2 + 1 ) .
( Q 11 Q 12 Q 13 Q 14 Q 15 Q 16 Q 21 Q 22 Q 23 Q 24 Q 25 Q 26 Q 31 Q 32 Q 33 Q 34 Q 35 Q 36 Q 41 Q 42 Q 43 Q 44 Q 45 Q 46 Q 51 Q 52 Q 53 Q 54 Q 55 Q 56 Q 61 Q 62 Q 63 Q 64 Q 65 Q 66 ) ( Δ F A Δ F B Δ ϕ A Δ ϕ B Δ n A Δ n B ) = Z g q γ γ e ( 0 0 0 0 Δ J A Δ J B ) .
{ Q 11 = 1 2 n A L λ Q 12 = κ e λ θ cos ψ A L Q 13 = κ F B L sin ( ψ A L ) Q 14 = e λ θ κ F B L sin ( ψ A L ) Q 15 = 1 2 F A L Q 16 = 0 Q 21 = κ e λ θ cos ψ B L Q 22 = 1 2 n B L λ Q 23 = e λ θ κ F A L sin ( ψ B L ) Q 24 = κ F A L sin ( ψ B L ) Q 25 = 0 Q 26 = 1 2 F B L Q 31 = 1 2 α n A L ( ω L ω A 0 ) Q 32 = κ e λ θ sin ( ψ A L ) Q 33 = F B L λ F A L Q 34 = e λ θ κ F B L Q 35 = 1 2 α F A L Q 36 = 0 Q 41 = κ e λ θ sin ( ψ B L ) Q 42 = 1 2 α n B L ( ω L ω B 0 ) Q 43 = e λ θ κ F A L Q 44 = κ F A L λ F B L Q 45 = 0 Q 46 = 1 2 α F B L Q 51 = 2 F A L ( n A L + 1 ) Q 52 = 0 Q 53 = 0 Q 54 = 0 Q 55 = ( F A L 2 + 1 ) T λ Q 56 = 0 Q 61 = 0 Q 62 = 2 F B L ( n B L + 1 ) Q 63 = 0 Q 64 = 0 Q 65 = 0 Q 66 = ( F B L 2 + 1 ) T λ
H A A ( λ ) = q γ γ e Z g Δ F A Δ J A = Q 1 ( 000010 ) .
H A B ( λ ) = q γ γ e Z g Δ F A Δ J B = Q 1 ( 000001 ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.