Next Article in Journal
The Combined-Unified Hybrid Censored Samples from Pareto Distribution: Estimation and Properties
Next Article in Special Issue
Reduced Dimensionality Multiphysics Model for Efficient VCSEL Optimization
Previous Article in Journal
Finding Location Visiting Preference from Personal Features with Ensemble Machine Learning Techniques and Hyperparameter Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Semiconductor Laser Linewidth Theory Revisited

1
Ferdinand-Braun-Institut gGmbH, Leibniz-Institut für Höchstfrequenztechnik, Gustav-Kirchhoff-Str. 4, 12489 Berlin, Germany
2
Weierstraß-Institut für Angewandte Analysis und Stochastik, Leibniz-Institut im Forschungsverbund Berlin e.V., Mohrenstr. 39, 10117 Berlin, Germany
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(13), 6004; https://doi.org/10.3390/app11136004
Submission received: 26 April 2021 / Revised: 27 May 2021 / Accepted: 3 June 2021 / Published: 28 June 2021
(This article belongs to the Special Issue Research and Application of Semiconductor Lasers)

Abstract

:
More and more applications require semiconductor lasers distinguished not only by large modulation bandwidths or high output powers, but also by small spectral linewidths. The theoretical understanding of the root causes limiting the linewidth is therefore of great practical relevance. In this paper, we derive a general expression for the calculation of the spectral linewidth step by step in a self-contained manner. We build on the linewidth theory developed in the 1980s and 1990s but look from a modern perspective, in the sense that we choose as our starting points the time-dependent coupled-wave equations for the forward and backward propagating fields and an expansion of the fields in terms of the stationary longitudinal modes of the open cavity. As a result, we obtain rather general expressions for the longitudinal excess factor of spontaneous emission (K-factor) and the effective α -factor including the effects of nonlinear gain (gain compression) and refractive index (Kerr effect), gain dispersion, and longitudinal spatial hole burning in multi-section cavity structures. The effect of linewidth narrowing due to feedback from an external cavity often described by the so-called chirp reduction factor is also automatically included. We propose a new analytical formula for the dependence of the spontaneous emission on the carrier density avoiding the use of the population inversion factor. The presented theoretical framework is applied to a numerical study of a two-section distributed Bragg reflector laser.

1. Introduction

Many applications of semiconductor lasers utilized in miniaturized contemporary photonic integrated devices for coherent optical communication, optical atomic clocks, atom interferometry, gravitational wave detection, space-based metrology, and optical quantum sensing impose strict requirements on the coherence of the light source, which can be expressed in terms of the spectral linewidth. The theoretical understanding of the governing factors limiting the linewidth is therefore of great practical relevance. After a first prediction by Schawlow and Townes [1], the corresponding theoretical framework was developed in the 1960s [2,3,4,5,6]). Later, in the 1980s and 1990s with the rapidly advancing technology of the fabrication of semiconductor lasers and their applications the theory of the spectral linewidth was further refined. The following five milestones can be distinguished:
(i)
The first milestone is the discovery of the enhancement of the fraction of the spontaneous emission going into the lasing mode in gain-guided lasers and the derivation of a corresponding excess factor (K-factor) by Petermann [7]. Siegman recognized this effect as a general property of non-Hermitian laser cavities [8]. Later his discussion of the power-nonorthogonality of the transversal modes and its consequences was extended to the case of the power-nonorthogonality of the longitudinal modes of laser cavities [9]. In Ref. [10] it was discovered, that the longitudinal modes can become degenerate for certain parameter configurations resulting in an infinite K-factor. The occurrence of such exceptional points is not restricted to lasers but is inherent to non-Hermitian systems, see [11] for a recent review.
(ii)
The second milestone is the discovery of a linewidth enhancement in semiconductor gain materials caused by refractive index fluctuations in response to fluctuations of the carrier density. Due to gain clamping, intensity fluctuations (which have negligible direct effect on the linewidth) can cause substantial refractive index changes, which in turn lead to fluctuations of the phase. The magnitude of this amplitude-phase coupling is quantified by the linewidth enhancement factor α H or α -factor introduced by Henry [12]. Later it was found that in DFB lasers the α -factor has to be replaced by an effective factor α H , eff [13] and a general expression for α H , eff valid for distributed feedback (DFB), distributed Bragg reflector (DBR) and external cavity lasers was derived in [14].
(iii)
The third milestone is the discovery of the possibility of linewidth reduction due to optical feedback from an external cavity by Patzak et al. [15]. Later a relation to the reduction of the frequency chirp was established [16,17].
(iv)
The forth milestone is related to the deterioration of the linewidth due to charge carriers injected into a phase tuning section within the cavity. Amann and Schimpe figured out that carrier noise is the origin resulting in an additional contribution to the linewidth, if the α -factors are different in the gain and phase tuning sections [18].
(v)
The last milestone is the discovery of the enhancement of the linewidth due to fluctuations of the shape of the profile of the optical power in the cavity by Tromborg and co-workers [14], which is particularly important in the vicinity of instabilities [19,20]. The most sophisticated linewidth theory including fluctuations of the shape of the power profile was published in [21,22].
In this paper, we will derive a semi-analytical expression of the spectral linewidth from a modern perspective step by step in a self-contained manner collecting all necessary ingredients that are otherwise found only scattered in the literature. The underlying theoretical approach is the classical Langevin formalism [23] where the deterministic equations for the optical field and the carrier density are supplemented by noise sources (Langevin forces). By means of these stochastic terms, quantum field theoretical phenomena (in particular spontaneous emission), which are essential for the laser linewidth, can be adequately treated within the framework of the semi-classical theory. The correlation functions of the Langevin noise sources are determined by the fluctuation-dissipation theorem, for which we refer to Refs. [24,25,26].
While the starting point of most authors is the Helmholtz equation and an expansion of its Green’s function [14,27], it is more transparent to start from the time-dependent coupled-wave equations and to expand the forward and backward propagating fields in terms of the stationary longitudinal modes of the open cavity [10] yielding simpler expressions. The linewidth expression derived in such a manner includes automatically the findings described as milestones (i)–(iii), i.e., longitudinal K-factor, effective α -factor, and chirp reduction factor. Additionally, we propose a new analytical formula for the dependence of the spontaneous emission on the carrier density, which is independent of the commonly employed population inversion factor, i.e., the ratio between the spontaneous rate of downward band-to-band transitions and the stimulated rate of downward and upward transitions [28]. The new formula avoids problems with the singularity of the population inversion factor near the transparency density and is in good agreement with microscopic calculations, see Section 9. The theory can be easily extended to include carrier noise (iv), too, as sketched in the Outlook in Section 11.

2. Prerequisites and Basic Assumptions

We consider edge-emitting semiconductor lasers as sketched in Figure 1 consisting of an arbitrary number of sections of different functionality. The transverse cross section is uniform within each section, but may vary from section to section. The preferred propagation direction of the optical field under lasing conditions is along the cavity axis parallel to the longitudinal coordinate z. Thus laser can be well described by the traveling-wave ansatz for the main component of the electric field strength
E ( x , y , z , t ) = e ( x , y ) 2 ε 0 c n Ψ + ( z , t ) e i β 0 z + Ψ ( z , t ) e i β 0 z e i ω 0 t + c . c . ,
employing a scalar approximation and neglecting the dynamics in the transverse ( x , y ) -plane. This means that we assume index-guiding in the transverse (lateral–vertical) plane and lasing in a single transverse mode and exclude gain-guiding, as appropriate for narrow-linewidth lasers. The real-valued field distribution e ( x , y ) is obtained as a solution of a real-valued waveguide equation normalized to e 2 ( x , y ) d x d y = 1 [29]. Moreover, ε 0 is the vacuum permittivity, c the vacuum speed of light, n the reference modal index, β 0 = 2 π n / λ 0 the reference propagation constant, λ 0 the reference wavelength and ω 0 = 2 π c / λ 0 the reference frequency. The prefactor is chosen such that Ψ 2 | Ψ + | 2 + | Ψ | 2 is the optical power P with unit W. The left and right propagating fields Ψ ± vary slowly both in the longitudinal coordinate z and time t and fulfill the stochastic time-dependent coupled-wave equations [30,31]
n g ( z ) c t Ψ ( z , t ) + σ z z + i M ( z , N , Ψ 2 ) + i D ^ ( z , t ) Ψ ( z , t ) = F sp ( z , t , N )
with
Ψ ( z , t ) = Ψ + ( z , t ) Ψ ( z , t ) , σ z = 1 0 0 1 ,
and
M ( z , N , Ψ 2 ) = Δ β ( z , N , Ψ 2 ) κ + ( z ) κ ( z ) Δ β ( z , N , Ψ 2 ) ,
Δ β ( z , N , Ψ 2 ) = β ( z , N , Ψ 2 ) β 0 = 2 π λ 0 Δ n ( z ) + β N ( z , N , Ψ 2 ) ,
β N ( z , N , Ψ 2 ) = 2 π λ 0 Δ n N ( z , N , Ψ 2 ) + i 2 g ( z , N , Ψ 2 ) α ( z , N , Ψ 2 ) .
The Langevin forces
F sp ( z , t , N ) = F sp + ( z , t , N ) F sp ( z , t , N ) ,
model the spontaneous emission into the forward and backward traveling waves and are assumed to have zero mean
F sp = 0 .
Here, n g is the real-valued group index, β is the propagation factor and Δ n is the static index detuning with respect to the reference modal index. β N describes the dependence of the propagation factor on the excess carrier density N, consisting of a real part proportional to Δ n N and an imaginary part with g being the modal gain and α the optical losses. κ ± are complex coupling coefficients describing a Bragg grating (coupling forward and backward traveling waves) and D ^ ( z , t ) is the dispersion operator. The ensemble or temporal average is denoted by · . The dependence of Δ n N , g and α on the optical power Ψ 2 is due to nonlinear effects like the Kerr effect, gain compression, and two-photon absorption, which affect the stationary states only slightly but can have a significant impact on the dynamic properties.
All quantities entering (2) are obtained by weighting the original quantities entering Maxwell’s equations with the intensity distribution e 2 ( x , y ) of the transverse mode. If we choose the origin of z such that the grating given by a complex dielectric function ε ( x , y , z ) has the property ε ( x , y , z ) = ε ( x , y , z ) , then κ + = κ = κ holds. Equation (2) has to be solved subject to the usual boundary conditions
Ψ + ( 0 , t ) = r 0 Ψ ( 0 , t ) and Ψ ( L , t ) = r L e 2 i β 0 L Ψ + ( L , t ) ,
with L being the total cavity length and r 0 and r L the complex-valued reflection coefficients at the facets at z = 0 and z = L , respectively. Scattering matrices establish transition conditions on Ψ at the interfaces between different sections [22,30].
The coupled-wave Equation (2) must be supplemented with an equation for the (excess) carrier density N. Although the theoretical treatment is quite general, we will specifically treat diode lasers in this paper, where the carrier dynamics is governed actually by both electrons and holes. We assume charge neutrality and consider only the excess carriers in the active region with density N = n n 0 = p p 0 (where n and p are the electron and hole densities and n 0 , p 0 are the corresponding equilibrium densities), neglecting any transport and capture effects. Thus, the excess carrier density is assumed to obey the rate equation
N t j ( z , N ) q d + R ( z , N ) + R st ( z , N , Ψ 2 ) = F N ( z , t , N )
with the rate of spontaneous and non-radiative recombination R and the rate of stimulated recombination
R st = Re Ψ * · [ g 2 i D ^ ] Ψ d W ω 0
which can be derived by multiplying (2) from the left with [ Ψ + * , Ψ * ] , the complex conjugate of (2) with [ Ψ + , Ψ ] and adding both equations to obtain a balance equation for the electromagnetic energy in differential form (Poynting’s theorem) [30]. Here the dot product means Ψ 1 · Ψ 2 = Ψ 1 + Ψ 2 + + Ψ 1 Ψ 2 and q is the elementary charge, the reduced Planck’s constant and d and W are the thickness and width, respectively, of the active region. The current density can be related to the (sectional) applied voltage by means of the Joyce model [32,33]
j ( z , N ) = U s U F ( N ) W L s R s
with sectional applied voltage U s , resistance R s , injection current I s , length L s , and the Fermi voltage U F . The carrier-density dependent injection current density given by (7) counteracts longitudinal spatial hole burning, which substantially lowers the degradation of the side mode suppression in strongly coupled DFB lasers compared to a model assuming a constant current density [34]. The frequency modulation response, albeit not topic of the present paper, is also affected as shown in [35].
Dispersion, i.e., the frequency dependence of the dielectric function, enters (2) first via the group index n g in front of the time derivative resulting from a linearization of the real part of the dielectric function with respect to the frequency ω . Second, the dispersion of the optical gain (and associated refractive index) is described in the time domain by the operator D ^ ( z , t ) , which can be obtained, e.g., by approximating the gain spectrum by a Lorentzian (or a series of Lorentzians) and a back-transformation into the time domain. This results in auxiliary differential equations for polarization functions [36,37], see Appendix A. Due to the fact, that the spectral width of the gain in semiconductor lasers is much larger than the width of the cavity resonances, D ^ ( z , t ) can be expanded again linearly in the frequency domain. In contrast to [36], we neglect here the dependence of the dispersion operator on the carrier density N.
Besides Ψ , we consider N as the only further stochastic variable and spontaneous emission as the only source of noise, neglecting all other noise sources. The theory can be easily extended to include carrier noise as sketched in the Outlook in Section 11. We assume that κ does not depend on N, i.e., we exclude gain-coupled lasers, and neglect fluctuations of the shape of the power profile. We neglect the impact of non-lasing side modes on the linewidth of the lasing mode and employ a single-mode approximation. Thus the linewidth rebroadening caused by mode partitioning and poor side-mode suppression, cf.  [38,39], can not be not accounted for.
Regarding the properties of the noise [40], we employ the usual assumptions of ergodicity (ensemble or statistical average equals temporal average) and stationarity (i.e.,  ξ ( t ) ξ ( t ) depends only on the time difference t t ) of any stochastic variable ξ . The fluctuations are assumed to have Gaussian probability distributions and delta-correlated covariance functions in time (Markov approximation) and space, i.e., ξ ( z , t ) ζ ( z , t ) = 2 D ξ ζ ( z ) δ ( t t ) δ ( z z ) with D ξ ζ being called diffusion coefficient. In order to obtain a Lorentzian shape of the optical spectrum around the lasing frequency, some approximations have to be employed which are described in Section 4. The full width at half maximum (FWHM) of the Lorentzian is commonly called intrinsic linewidth. As we consider only spontaneous emission noise, we calculate the quantum limit of the intrinsic linewidth. If non-Markovinan noise is considered, which is beyond of the scope of the paper, the lineshape tends towards a Gaussian in the vicinity of the lasing frequency with a correspondingly larger FWHM [41,42].

3. Equation of the Field Amplitude

Throughout this work, we restrict ourselves to the analysis of the noise at continuous wave (CW) emission, i.e., at steady-state lasing. Therefore, we expand Ψ ( z , t ) in terms of the stationary eigenmodes Φ m ( z ) (longitudinal modes),
Ψ ( z , t ) = m f m ( t ) Φ m ( z ) ,
satisfying the equation
σ z z + i M z , N ¯ , Ψ 2 ¯ + i β D ( z , Ω m ) + i n g ( z ) c Ω m Φ m ( z ) = 0
with
M z , N ¯ , Ψ 2 ¯ = Δ β z , N ¯ , Ψ 2 ¯ κ ( z ) κ ( z ) Δ β z , N ¯ , Ψ 2 ¯ M ¯ ,
Δ β z , N ¯ , Ψ 2 ¯ = β z , N ¯ , Ψ 2 ¯ β 0 β ¯ β 0 = Δ β ¯ .
Here, N ¯ and Ψ 2 ¯ are the steady-state carrier density and optical power, respectively, at which the operator of the eigenvalue problem (9) is expanded ( β ¯ , M ¯ are the corresponding propagation factor and matrix M). The complex eigenvalue Ω m describes both, the frequency deviation Re ( Ω m ) = ω m ω 0 (or the wavelength deviation Δ λ m = Re ( Ω m ) d λ / d ω | λ 0 ), and the mode damping rate Im ( Ω m ) . The dispersion β D ( z , Ω m ) is given by the operator D ( z , t ) in the frequency domain. Equation (9) has to be solved subject to the boundary conditions
Φ m + ( 0 ) = r 0 Φ m ( 0 ) Φ m ( L ) = r L e 2 i β 0 L Φ m + ( L )
following from (4).
In the case of stationary eigenmodes considered here, the two-point boundary value problem (9) together with the boundary conditions (12) is in fact a quasi-linear eigenvalue problem. The point of expansion Ψ 2 ¯ , N ¯ of the nonlinear operator must be chosen self-consistently, i.e., it is required to satisfy the conditions (36) and (68) for the (mean) stationary state given below. Furthermore, in the presence of dispersion β D ( Ω m ) , the eigenvalue problem is also nonlinear in the eigenvalue. The problem is similar to the self-consistent field method in electronic structure theory [43] and can be treated in a similar way (e.g., linearization of the optical-power and the frequency dependency of the gain dispersion and fixed-point iteration, where the point of expansion is updated to the most recent stationary state in each step). Alternatively, the system can be directly propagated in the time-domain (as it has been done in Section 10), until convergence to a suitable CW state has been achieved.
Note that different modes (eigensolutions of (9)) fulfill the orthogonality relation
( Φ m , β ω Φ n ) = 0 for m n
which can be derived by multiplying (9) from the left by [ Φ n , Φ n + ] , the corresponding equation for Φ n by [ Φ m , Φ m + ] , integrating along z, using (12) and subtracting both equations. Here the short notation
( Φ , Ψ ) = 0 L ( Φ + Ψ + Φ Ψ + ) d z .
denotes an inner product, distinguished from a standard Hilbert space scalar product. At an exceptional point the expansion (8) includes also a generalized eigenmode [10]. For m = n , the factor β / ω , which can be interpreted as an inverse complex-valued group velocity, is the exact derivative
β ω = n g c + β D ω | Ω m .
In the general case of m n this holds only approximately:
β ω = n g c + β D ( Ω m ) β D ( Ω n ) Ω m Ω n n g c + β D ω | Ω m .
A system of ordinary differential equations for the amplitudes f m ( t ) can be obtained by inserting (8) into (2), using (9), multiplying from left with [ Φ n , Φ n + ] , and integrating along z. In order to exploit the orthogonality relation (13) with the approximation (16), we have again to expand
β D ( ω ˜ ) β D ( Ω m ) = β D ω | Ω m ( ω ˜ Ω m ) ,
where ω ˜ = ω ω 0 . Note that this is consistent with the slowly varying amplitude approximation. Transforming back into the time domain yields
β D ( ω ˜ ) β D ( Ω m ) f m ( ω ˜ ) e i ω ˜ t d ω ˜ = D ^ β D ( Ω m ) f m ( t )
and
β D ω | Ω m ( ω ˜ Ω m ) f m ( ω ˜ ) e i ω ˜ t d ω ˜ = β D ω | Ω m Ω m f m ( t ) + i t f m ( ω ˜ ) e i ω ˜ t d ω ˜ = β D ω | Ω m Ω m f m ( t ) + i f m t .
Therefore the amplitudes f m fulfill
( Φ m , β ω Φ m ) f m t i Ω m ( Φ m , β ω Φ m ) f m + i n ( Φ m , Δ M Φ n ) f n = ( Φ m , F sp )
with
Δ M = M M ¯ = ( β β ¯ ) 1 0 0 1 .
In what follows, we employ the single mode approximation and drop the subscript m. Then (20) becomes a single differential equation for the complex-valued amplitude f of the lasing mode
f t i Ω f + i ( Φ , Δ M Φ ) ( Φ , β ω Φ ) f = F f .
The new Langevin noise source entering (22) is given by
F f = ( Φ , F sp ) ( Φ , β ω Φ ) .
For an analysis above threshold, it is beneficial to derive equations for the modulus | f | and the phase φ of f, because the fluctuations of | f | are damped but those of φ are not. For the modulus squared | f | 2 (called intensity in the rest of the paper) it follows
| f | 2 t + 2 Im ( Ω ) | f | 2 2 Im ( Φ , Δ M Φ ) ( Φ , β ω Φ ) | f | 2 = 2 Re ( F f f * ) .
An equation for the phase can be determined from
Im f t f * Re ( Ω ) | f | 2 + Re ( Φ , Δ M Φ ) ( Φ , β ω Φ ) | f | 2 = Im ( F f f * )
and
f t = e i φ | f | t + i f φ t Im f t f * = | f | 2 φ t .
The stochastic term on the right-hand side of (24) describes intensity noise, which does not have a vanishing expectation value anymore: F f f * 0 . It is therefore convenient to rewrite the term by distilling out the expectation value and introducing two new real-valued Langevin forces F | f | 2 and F φ with zero mean as
2 F f f * = 2 F f f * + F | f | 2 + 2 i | f | 2 F φ .
It will be shown below that F f f * is real-valued, see (66). Then, the final equations for the modulus squared and phase of f read
| f | 2 t = 2 Im ( Ω ) | f | 2 + 2 Im ( Φ , Δ M Φ ) ( Φ , β ω Φ ) | f | 2 + 2 F f f * + F | f | 2
and
φ t = Re ( Ω ) Re ( Φ , Δ M Φ ) ( Φ , β ω Φ ) + F φ .
Now we consider small fluctuations
N ( z , t ) = N ( z ) + δ N ( z , t ) ,
φ ( t ) = φ ( t ) + δ φ ( t ) ,
| f ( t ) | 2 = | f | 2 + δ | f ( t ) | 2
and expand
β = β ( N , | f | 2 ) + β N | N , | f | 2 δ N + β | f | 2 | N , | f | 2 δ | f | 2
around the mean values
N ( z ) = N ¯ ( z ) , φ ( t ) = Re ( Ω ) t .
The mean intensity | f | 2 is required to satisfy
Im ( Ω ) | f | 2 = F f f * .
Due to the extreme smallness of the spontaneous emission F f f * going into the lasing mode, (35) can be replaced by
Im ( Ω ) = 0
above threshold ( | f | 2 0 ) , which will be exploited in what follows (i.e., Ω = Re ( Ω ) ). The intensity and phase fluctuations fulfill the equations
δ | f | 2 t = 2 Im ( Φ , β N δ N Φ ) + ( Φ , β | f | 2 Φ ) δ | f | 2 ( Φ , β ω Φ ) | f | 2 + F | f | 2
and
δ φ t = Re ( Φ , β N δ N Φ ) + ( Φ , β | f | 2 Φ ) δ | f | 2 ( Φ , β ω Φ ) + F φ ,
which are the basis for the derivation of the linewidth formula in the following sections.

4. Lorentzian Line Shape

The power spectral density (PSD) of the optical field at z = 0 can be calculated utilizing the Wiener–Khinchin theorem as [44]
S E ( ω ) = + E * ( x , y , 0 , t ) E ( x , y , 0 , t + τ ) e i ω τ d τ d x d y
with
E * ( x , y , 0 , t ) E ( x , y , 0 , t + τ ) = 1 | r 0 | 2 2 ε 0 c n e 2 ( x , y ) f * ( t ) f ( t + τ ) | Φ ( 0 ) | 2 e i ω 0 τ + c . c .
We do not address here technical issues regarding the non-existence of Fourier integrals of fluctuating quantities here and refer to the discussion in [44]. Hence we have to calculate the PSD
S f ( ω ) = + f * ( t ) f ( t + τ ) e i ( ω ω 0 ) τ d τ
with
f * ( t ) f ( t + τ ) = | f ( t ) | | f ( t + τ ) | e i [ φ ( t + τ ) φ ( t ) ] .
Now we employ here a first approximation: We neglect in (42) the intensity fluctuations because above threshold they are damped, in contrast to the phase fluctuations, which can been seen by comparing (37) and (38), see the discussion in [45]. Therefore, we obtain
f * ( t ) f ( t + τ ) | f | 2 e i Δ τ φ e i Ω τ
with
Δ τ φ = δ φ ( t + τ ) δ φ ( t )
and Ω = Re ( Ω ) . The random variable Δ τ φ is assumed to be Gaussian distributed [46], such that it holds
e i Δ τ φ = 1 2 π Var ( Δ τ φ ) e ( Δ τ φ ) 2 2 Var ( Δ τ φ ) e i Δ τ φ d Δ τ φ = e 1 2 Var ( Δ τ φ ) .
The variance of the phase fluctuations is
Var ( Δ τ φ ) = Δ τ φ Δ τ φ 2 = ( Δ τ φ ) 2 = δ φ ( t + τ ) 2 2 δ φ ( t ) δ φ ( t + τ ) + δ φ ( t ) 2 = 2 δ φ ( 0 ) 2 δ φ ( 0 ) δ φ ( τ ) ,
where we exploited Δ τ φ = 0 and the stationarity property. We introduce the PSD of phase fluctuations which is even in the frequency ( S δ φ ( ω ˜ ) = S δ φ ( ω ˜ ) )
S δ φ ( ω ˜ ) = + δ φ ( t ) δ φ ( t + τ ) e i ω ˜ τ d τ
and the PSD of optical frequency fluctuations δ ω ( ω ˜ ) = i ω ˜ δ φ ( ω ˜ )
S δ ω ( ω ˜ ) = + δ ω ( t ) δ ω ( t + τ ) e i ω ˜ τ d τ = ω ˜ 2 S δ φ ( ω ˜ ) .
With these spectral densities, the variance of the phase fluctuations can be written as
Var ( Δ τ φ ) = 1 π S δ φ ( ω ˜ ) 1 e i ω ˜ τ d ω ˜ = 1 π S δ φ ( ω ˜ ) 2 sin 2 ω ˜ τ 2 i sin ω ˜ τ d ω ˜ = 2 π S δ φ ( ω ˜ ) sin 2 ω ˜ τ 2 d ω ˜ = τ 2 2 π S δ ω ( ω ˜ ) sin ω ˜ τ 2 ω ˜ τ 2 2 d ω ˜ .
We carry out a second approximation by noting that the function sin 2 ω ˜ τ 2 / ( ω ˜ τ 2 ) 2 peaks strongly at ω ˜ = 0 for large | τ | , and obtain
Var ( Δ τ φ ) | τ | π S δ ω ( 0 ) sin 2 x x 2 d x = | τ | S δ ω ( 0 ) .
In this step, we have recovered a central property of Brownian motion, as we can observe that the root mean square displacement of the phase fluctuation grows as | τ | . Inserting (50) into (45) and, finally, the result into (43) yields
f * ( t ) f ( t + τ ) | f | 2 e 1 2 S δ ω ( 0 ) | τ | e i Ω τ
and
S f ( ω ) | f | 2 + e 1 2 S δ ω ( 0 ) | τ | + i ( ω 0 + Ω ω ) τ d τ = | f | 2 S δ ω ( 0 ) ω ( ω 0 + Ω ) 2 + S δ ω ( 0 ) 2 2 .
This is a Lorentzian centered at ω = ω 0 + Ω with the FWHM
Δ ω = S δ ω ( 0 )
given by the PSD of the optical frequency fluctuations taken at zero Fourier frequency.
In the following sections of this paper, we will focus on the static white noise limit, where the frequency noise PSD S δ ω ( ω ) is effectively approximated by a constant, such that the second approximation carried out on the step from (49) to (50) becomes in fact redundant.

5. Correlation Functions

The correlation function of the spontaneous emission noise F f defined in (23) is
F f * ( t ) F f ( t ) = ( Φ * , F sp * ( z , t ) ) ( Φ , F sp ( z , t ) ) ( Φ , β ω Φ ) 2 .
The numerator reads
( Φ * , F sp * ( z , t ) ) ( Φ , F sp ( z , t ) ) = 0 L 0 L Φ + * F sp * ( z , t ) + Φ * F sp + * ( z , t ) Φ + F sp ( z , t ) + Φ F sp + ( z , t ) d z d z = 0 L 0 L | Φ + | 2 F sp * ( z , t ) F sp ( z , t ) + | Φ | 2 F sp + * ( z , t ) F sp + ( z , t ) d z d z = 2 0 L Φ ( z ) 2 D sp ( z ) d z δ ( t t )
with
Φ 2 = | Φ + | 2 + | Φ | 2
taking into account the noise covariance functions [25,26,27] for index and absorption coupling
F sp ± * ( z , t ) F sp ± ( z , t ) = 2 D sp ( z ) δ ( t t ) δ ( z z ) , F sp * ( z , t ) F sp ± ( z , t ) = 0 , F sp ± ( z , t ) F sp ± ( z , t ) = 0 , F sp ( z , t ) F sp ± ( z , t ) = 0
with
2 D sp ( z , N , Ψ 2 ) = ω 0 n sp ( N ) g ( z , N , Ψ 2 ) ,
where n sp is the population inversion factor (a Bose–Einstein distribution function multiplied by 1 ). The correlation functions for gain-coupling can be found in [47].
Since Ψ has the unit W , D sp must have the unit Ws/m as it is. The diffusion coefficient D f * f defined by
F f * ( t ) F f ( t ) = 2 D f * f δ ( t t ) .
is therefore
D f * f = D f * f ( N , Ψ 2 ) = 0 L Φ ( z ) 2 D sp ( z , N , Ψ 2 ) d z ( Φ , β ω Φ ) 2 .
The correlation functions of the Langevin forces entering (37) and (38) read
F | f | 2 ( t ) F | f | 2 ( t ) = 2 D | f | 2 | f | 2 δ ( t t ) = 4 D f * f | f | 2 δ ( t t ) ,
F φ ( t ) F φ ( t ) = 2 D φ φ δ ( t t ) = D f * f | f | 2 δ ( t t ) ,
F | f | 2 ( t ) F φ ( t ) = 2 D | f | 2 φ δ ( t t ) = 0 ,
which can be derived using the transformation rules for the diffusion coefficients [23]
D | f | 2 | f | 2 = | f | 2 f | f | 2 f * + | f | 2 f * | f | 2 f D f * f , D φ φ = φ f φ f * + φ f * φ f D f * f , D | f | 2 φ = | f | 2 f φ f * + | f | 2 f * φ f D f * f ,
and the relations | f | 2 = f * f and φ = i 2 ln ( f / f * ) . Following [48], the ensemble average defined in (27) can be obtained using the formal solution of (22) in the form of
f t = f t Δ t + t Δ t t f * ( t ) t d t ,
where the time step 0 < Δ t γ R 1 is assumed to be short in comparison to the (inverse) relaxation rate γ R = Im Ω Φ , Δ M Φ / Φ , ω β Φ . Substituting (65) into F f t f * t yields
F f t f * t = F f t f * t Δ t + t Δ t t F f t f * ( t ) t d t = F f t f * t Δ t i t Δ t t Ω Φ , Δ M Φ Φ , β / ω Φ * F f t f * t d t = + t Δ t t F f t F f * t d t ,
where we have used (22) in the second step. The first term vanishes because of causality, since the fluctuations of f * t Δ t are not correlated with the future noise F f t . Consequently, the correlation function factorizes and results in zero due to the zero-mean property of the Langevin force. The second term can be dropped since it is non-zero only over a set of measure zero (i.e., at t = t ). Finally, we obtain with (59)
F f t f * t = t Δ t t F f t F f * t d t = D f * f ,
where a factor 2 is canceled by 1/2 encountered in integrating only half of the δ -function.

6. Effective Linewidth Enhancement Factor

Using the ansatz (8), we can approximate the rate of stimulated emission (6) as
R st = g + 2 Im β D ( Ω ) d W ω 0 | f | 2 Φ 2 .
The stationary mean value of the carrier density and its fluctuations satisfy
j ( N ) e d = R ( N ) + g ( N , | f | 2 Φ 2 ) + 2 Im β D ( Ω ) d W ω 0 | f | 2 Φ 2
and
δ N t + δ N τ d + R st | f | 2 δ | f | 2 = F N ,
respectively, with the inverse differential carrier lifetime (taken at N and | f | 2 )
1 τ d = 1 q d j N + R N + R st N .
The FWHM of the Lorentzian line shape is given by the spectral density of the frequency fluctuations at zero Fourier frequency, cf. (53). This is equivalent to setting the time derivatives of the fluctuations equal to zero (static limit)
δ | f | 2 t = 0
and
δ N t = 0 .
Furthermore, we omit carrier noise ( F N = 0 ), considering only spontaneous emission noise as mentioned in Section 2. With these approximations, the carrier density fluctuations are directly related to the intensity fluctuations by
δ N = τ d R st | f | 2 δ | f | 2 .
We substitute (73) into (37) and obtain
δ | f | 2 = F | f | 2 2 | f | 2 Im ( h α )
with
h α = ( Φ , τ d β N R st | f | 2 Φ ) ( Φ , β | f | 2 Φ ) ( Φ , β ω Φ ) .
Inserting (74) back into (73) yields
δ N = τ d R st | f | 2 F | f | 2 2 | f | 2 Im ( h α ) .
Finally, substituting (76) into (38) results in
δ φ t = F | f | 2 2 | f | 2 Re ( h α ) Im ( h α ) + F φ ,
which can be considered as a defining equation for the so-called effective Henry’s α -factor
α H , eff = Re ( h α ) Im ( h α ) .
If α H , eff 0 , the Langevin force resulting in fluctuations of the intensity leads also to fluctuations of the phase. By virtue of the different integrals in the numerator and denominator of (75), the carrier dependence of the propagation factor, nonlinear gain (gain compression), absorption (two-photon absorption) and index (Kerr effect) as well as gain dispersion are accounted for. The impact of longitudinal spatial hole burning and multi-section cavity structures (e.g., including Bragg gratings, passive sections or external cavities) is also correctly described. If all of these additional effects are neglected, for a Fabry–Pérot (FP) laser the α -factor is recovered
α H = Re β N Im β N = 4 π λ 0 Δ n N N g N α N .
The second equality is obtained using (3c). The minus sign in (78), resulting in a corresponding minus sign in (79), has been chosen in agreement with the usual convention ensuring positive values of α H around the gain peak [49]. The derivative α / N of the modal absorption in (79) is often neglected, cf. [49]. One should keep in mind that α H is not a constant but depends on the carrier density (and the wavelength) because modal index and gain vary differently with carrier density.

7. Spectral Linewidth

Collecting the results of the previous sections, we can now calculate the spectral linewidth. Rewriting (53) as
Δ ω = + δ φ ( t ) t δ φ ( t + τ ) t d τ
and inserting the equation of the phase fluctuations (77), we obtain
Δ ω = α H , eff 2 4 | f | 2 2 F | f | 2 ( t ) F | f | 2 ( t + τ ) d τ + F φ ( t ) F φ ( t + τ ) d τ ,
where we have taken (63) into account. Inserting (61) and (62) results in the remarkably simple expression for the spectral linewidth
Δ ν = Δ ω 2 π = D f * f 2 π | f | 2 1 + α H , eff 2 .
In the following, we will rewrite this expression in terms of the intra-cavity photon number
I ph = 1 ω 0 0 L P Re β ω d z = P 0 ω 0 0 L Φ 2 Re β ω d z ( 1 | r 0 | 2 ) | Φ ( 0 ) | 2
where
P ( z ) = | f | 2 Φ ( z ) 2
is the stationary optical power inside the cavity and
P 0 = ( 1 | r 0 | 2 ) | f | 2 | Φ ( 0 ) | 2
is the outcoupled power at z = 0 . Moreover, we define the rate of spontaneous emission into the lasing mode
R sp = 0 L Φ 2 n sp g d z 0 L Φ 2 Re β ω d z
and the longitudinal excess factor of spontaneous emission
K = 0 L Φ 2 Re β ω d z 2 ( Φ , β ω Φ ) 2 .
Now, (82) can be written as
Δ ν = K R sp 4 π I ph 1 + α H , eff 2
which is the standard form found in the literature [14,46].
The mode profile Φ ( z ) entering (83), (86), and (87) is obtained by solving (9), (36) and (68) self-consistently above threshold. Thus, Φ is a mode of the active cavity affected by spatial hole burning. The expression for the K-factor generalizes the one given in [10] for non-uniform and complex-valued group velocity ( β / ω ) 1 and is in basic correspondence with [14,50,51], but fundamentally different from the modified K-factor introduced in [52].
At an exceptional point, the mode is orthogonal to itself, i.e., ( Φ , β ω Φ ) = 0 , and the K-factor approaches infinity. However, one has to keep in mind that at an exceptional point the single-mode approximation used here fails. Instead, for the description of the dynamics in the vicinity of such a point one has use at least one eigenmode and the corresponding generalized eigenmode [10].

8. Impact of a Passive Section, Chirp Reduction Factor and the Fabry–Pérot Case

In this section, we apply the general theory developed above to special configurations that allow for further simplifications (i.e., neglect of spatial hole burning, nonlinear gain and index as well gain dispersion). Thereby, it is shown that the general framework entails some well-established formulas from the literature as special cases.

8.1. Linewidth of a Laser Consisting of a Gain Chip Subject to Feedback from an External Cavity

We consider a laser consisting of one active section for z [ 0 , l ] with length l and an arbitrary number of passive sections (including an external cavity) for z [ l , L ] with total length L l . The active section is of the Fabry–Pérot (FP) type having no Bragg grating ( κ = 0 for z [ 0 , l ] ) . The passive sections may contain Bragg gratings but no excess carriers ( N = 0 , β / N = 0 , D sp = 0 , R st = 0 for z [ l , L ] ). The reference plane is located just inside the active section so that a possible finite reflectivity at the interface between active and passive sections belongs to the passive sections. Furthermore, we neglect spatial hole burning (i.e., N / z = 0 ), nonlinear gain and index ( β / | f | 2 = 0 , R st / | f | 2 = R st / | f | 2 ) and gain dispersion ( β D = 0 , β / ω = n g / c ). See Figure 2 for an illustration of the device.
The functions Φ ± ( z ) as solution of (9) are defined on the whole cavity z [ 0 , L ] . In the active section they are given by
Φ ± ( z ) | active = Φ ± ( l ) e i ( Δ β ¯ + n g c Ω ) ( z l )
and hence
Φ + ( z ) Φ ( z ) | active = Φ + ( l ) Φ ( l )
holds. Furthermore we define the left and right reflectivities
r + = Φ ( l ) Φ + ( l ) and r = Φ + ( l ) Φ ( l ) r 0 e 2 i l Δ β ¯ + n g c Ω ,
where we exploited the fact that the Φ ± ( z ) are continuous at z = l . The derivative of r + can be expressed as (see Appendix B)
ln ( r + ) ω = 2 i l L Φ + Φ β ω d z Φ + ( l ) Φ ( l ) .
Therefore,
( Φ , β ω Φ ) = 2 0 l Φ + Φ β ω d z + 2 l L Φ + Φ β ω d z = 2 Φ + ( l ) Φ ( l ) l β ω | active + i 2 ln ( r + ) ω = 2 l Φ + ( l ) Φ ( l ) β ω | active χ
where we introduced the parameter
χ = 1 + i ln ( r + ) ω 2 l β ω = 1 + i c 2 l n g ln ( r + ) ω
as in [14]. We omit | active in what follows. Furthermore
( Φ , τ d β N R st | f | 2 Φ ) = 2 | f | 2 0 l Φ + Φ τ d β N R st d z = 2 | f | 2 Φ + ( l ) Φ ( l ) β N 0 l τ d R st d z
holds. Then from (75), (78) and (79)
α H , eff = Re β N β ω χ * Im β N β ω χ * = Re β N χ * Im β N χ * = Re α H χ + i χ Im α H χ + i χ
and
1 + α H , eff 2 = ( 1 + α H 2 ) | χ | 2 Im 2 α H χ + i χ
follow. Similarly, we obtain
D f * f = 0 l Φ 2 D sp d z 2 0 l Φ + Φ β ω d z 2 1 χ 2 .
Hence, the product of (97) and (98) is given by
D f * f 1 + α H , eff 2 = 0 l Φ 2 D sp d z 2 0 l Φ + Φ β ω d z 2 1 + α H 2 Im 2 α H χ + i χ
and the spectral linewidth can be written as
Δ ν = Δ ν FP F 2
with the chirp reduction factor
F = Im α H χ + i χ = 1 c 2 n g l Im ln ( r + ) ω + α H c 2 n g l Re ln ( r + ) ω
and
Δ ν FP = 0 l Φ 2 D sp d z 2 π | f | 2 2 0 l Φ + Φ β ω d z 2 1 + α H 2
being the linewidth of a FP laser having the cavity length l and the reflection coefficient r + ( Ω ) at the rear facet. Equation (101) agrees with [53,54] as well as [16,55] if the differently chosen harmonic time dependence ( e i ω 0 t as used here versus e i ω 0 t ) is observed. The relation to the static frequency chirp is given in Appendix C. Introducing an effective length of the passive cavity
l p = c 2 n g , p Im ln ( r + ) ω ,
the chirp reduction factor can be written as
F = 1 + n g , p l p n g l + α H c 2 n g l Re ln ( r + ) ω .
Thus the second term on the right-hand side of (101) can be interpreted as the ratio of the round-trip times in the passive and active sections. The chirp reduction factor allows an estimation of the reduction of the linewidth by a passive section or an external cavity for given group index, length and Henry’s α -factor of the active section.

8.2. Linewidth of the Fabry–Pérot Laser Cavity

We start from (88) and employ the same approximations resulting in (100) and (101). The rate of spontaneous emission into the lasing mode (86) is
R sp ( N ) = 0 l Φ ( z ) 2 n sp ( N ) g ( N ) d z 0 l Φ ( z ) 2 Re β ω d z = c n g n sp ( N ) g ( N ) .
The expression (105) should be compared with
R ˜ sp ( N ) = β sp V R rad ( N )
often used in rate equation based modeling [56], where β sp is the dimensionless spontaneous emission factor (ratio between the spontaneous emission going into the lasing mode and the spontaneous emission into all modes), V the volume of the active region and R rad = B N 2 the rate of radiative spontaneous recombination (B bimolecular recombination coefficient). Equalizing (105) and (106) at the lasing threshold yields
β sp = c n g n sp ( N th ) g ( N th ) V R rad ( N th ) ,
which is of the order 10 5 for typical edge-emitting lasers. The photon number (83) can be written as
I ph = P 0 ω 0 0 l Φ ( z ) 2 Re β ω d z ( 1 | r 0 | 2 ) | Φ ( 0 ) | 2 = n g c P 0 ω 0 ( 1 | r 0 | 2 ) | Φ ( 0 ) | 2 0 l | Φ + ( 0 ) | 2 e 2 Im ( β ¯ ) z + | Φ ( 0 ) | 2 e 2 Im ( β ¯ ) z d z = n g c P 0 ω 0 | r 0 | 2 e 2 Im ( β ¯ ) l 1 e 2 Im ( β ¯ ) l + 1 2 Im ( β ¯ ) ( 1 | r 0 | 2 ) .
The complex mode frequencies of the FP cavity (obtained from (89) applying the boundary conditions (12) at z = 0 and z = l ) are the solutions of
Ω m = c n g m π l β ¯ i 2 l ln r 0 r + ( Ω m )
with m Z . To ensure Im ( Ω ) = 0 for the lasing mode, (109) leads to the usual threshold condition
α out 1 l ln ( | r 0 r + ( Re ( Ω ) ) | ) = 2 Im ( β ¯ ) g α .
In what follows, we set r l = r + ( Ω ) (reflectivity seen by the laser at z = l ). From (110) it follows
I ph = n g c P 0 ω 0 | r 0 | 2 1 | r 0 r l | 1 | r 0 r l | + 1 2 Im ( β ¯ ) ( 1 | r 0 | 2 ) = n g c P out ω 0 α out ,
where we used the relation between outcoupled power P 0 and total output power P out
P 0 = P out 1 + 1 | r l | 2 1 | r 0 | 2 | r 0 | | r l | .
The final result is
Δ ν FP = K FP 4 π c 2 n g 2 ω 0 n sp ( α + α out ) α out P out 1 + α H 2 ,
where Petermann’s K-factor
K = 0 l Φ 2 Re β ω d z 2 2 0 l Φ + Φ β ω d z 2
of the FP cavity can be similarly derived as
K FP = ( | r 0 | + | r l | ) ( 1 | r 0 r l | ) 2 | r 0 r l | ln ( | r 0 r l | ) 2
in agreement with [9,27]. Equation (113) together with (112) was used by several authors, e.g., [12,55,57]. However, one should be aware of the approximations involved in the derivation.
The original Schawlow–Townes formula [1] is obtained by setting K FP = 1 , n sp = 1 , α = 0 , α H = 0 , and relating the outcoupling losses α out to the spectral FWHM of a cavity resonance by Δ ν cav = c / n g · α out / 2 π to
Δ ν ST = π ω 0 P out ( Δ ν cav ) 2 .
There are two reasons why this differs from the original formula by a factor of 4. First, in [1] the half widths instead of full widths are used (factor of 2). Second, in [1] the sub-threshold case is considered resulting in a further factor of 2, which can be seen as follows: Below threshold the spectral density can be directly calculated from (22) for Δ M = 0 by Fourier transformation via
f * ( ω ) f ( ω ) = 2 π S f ( ω ) δ ( ω ω )
with the result
S f ( ω ) = 2 D f * f ω Re ( Ω ) 2 + Im ( Ω ) 2 = 2 D f * f ( ω Re ( Ω ) ) 2 + D f * f | f | 2 2 ,
where we have used (35) and (66). This is a Lorentzian with the FWHM
Δ ν subthr = D f * f π | f | 2 ,
which is a factor of 2 larger than (82) (for α H , eff = 0 ). See [6] for a more thorough discussion.

9. Population Inversion Factor

The population inversion factor (sometimes misleadingly called ‘spontaneous emission factor’) introduced in (58) and used in (86) and (113) is given by [28]
n sp ( N ) = 1 1 exp ω 0 q U F ( N ) k B T ,
where U F is the Fermi voltage (spacing of quasi-Fermi potentials of holes and electrons), q the elementary charge, k B the Boltzmann constant and T the temperature. It has a singularity at the transparency density N tr , where ω 0 = U F ( N tr ) and g ( N tr ) = 0 . Therefore, the replacement
r sp ( N ) n sp ( N ) g ( N ) = n sp ( α + α out )
with constant n sp often employed is critical because for high-Q-cavities with low α and α out the spontaneous emission is underestimated. We calculated the modal gain, the spontaneous emission into the lasing mode and the inversion factor employing a numerical simulation based on a 8 × 8   k · p band structure calculation and a free carrier theory for the optical response functions with phenomenological corrections for many-body effects such as band-gap renormalization, Coulomb enhancement and transition broadening [58]. The results for a 5 nm thick InGaAs quantum well emitting around 1064 nm embedded into an AlGaAs-based waveguide structure are compared with analytical models in Figure 3. The increase of the absorption (negative gain) at small increasing carrier densities observed in the simulation is caused by band-gap narrowing, shifting the absorption edge to longer wavelengths [59].
The gain at a fixed wavelength is modeled as
g = g N tr ln max ( N , N cl ) N tr
and the modal spontaneous emission as
r sp = g N tr 2 ln 1 + N N tr 2
with differential gain g = 19 · 10 22 m 2 , transparency density N tr = 1.7· 10 24 m 3 and gain clamping density N cl = 7 · 10 23 m 3 . Equation (123) exhibits the correct asymptotic behavior r sp N 2 for N N tr and r sp = g for N . The agreement of r sp between simulation and model is remarkably good without the need of introducing any new parameters. The singularity of the inversion factor at N = N tr is clearly visible in Figure 3. If a constant value of n sp = 2 is used, the spontaneous emission is considerably under- or over-estimated (depending on the respective carrier density).

10. Numerical Results for a DBR Laser

In the simulations described below we study a simple two-section DBR laser as sketched in Figure 1 and Figure 2. It consists of a 1 mm long active gain section (without Bragg grating) and a 3 mm long passive Bragg reflector section (without active layer) with a coupling coefficient of κ = 2 cm 1 . The full set of parameters characterizing the simulated DBR laser is listed in Table 1. The usage of such DBR lasers having long reflector sections with low coupling coefficients was suggested in Ref. [16] for the first time. In fact, in Ref. [57] an intrinsic linewidth of 2 kHz at an output power of 180 mW was reported. This device had a 3 mm long gain section and a 1 mm long reflector section, where the active layer extended over the whole cavity (i.e., including the reflector section). The device studied here resembles the one presented in Ref. [60], where an intrinsic linewidth of 4 kHz at an output power of 73 mW was reported. Experimentally, the intrinsic linewidth is determined from the plateau of the frequency noise spectrum, as in the theory [61].
Within the active section, we assume the logarithmic gain model (122) supplemented with the nonlinear gain saturation factor
g ( z , N , Ψ 2 ) = g ( z ) N tr ( z ) ln max ( N , N cl ( z ) ) N tr ( z ) 1 1 + Ψ 2 / P sat ( z ) ,
and the refractive index model with square-root-like dependency on the carrier density
Δ n N = λ 0 α ˜ H g N tr 2 π max ( N , N cl , i ) N tr 1 .
Note, that at the transparency carrier density N tr , the fixed parameter α ˜ H agrees with the α -factor α H from (79), whereas Δ n N vanishes. The α -factor correspondingly increases for N > N tr . The value of α ˜ H = 1 was chosen in agreement with simulations [58] and is in basic agreement with measurements, where α -factors between 1 and 2 for a highly-strained InGaAs quantum well were obtained, depending on the detuning [62].
Other changes of the refractive index, including self-heating-induced self- and cross-heating contributions, are modeled by [30]
Δ n ( z ) = Δ n 0 ( z ) + Δ n T ( z ) , Δ n T ( z ) = ν s I , z [ 0 , l ] ν c I , z [ l , L ] ,
where I is the injection current into the active (gain) section. We employ the commonly used cubic model for the rate of non-radiative and spontaneous recombination of the carriers in the active section
R ( N ) = A N + B N 2 + C N 3 ,
is approximated by describing Shockley–Read–Hall recombination, direct band-to-band recombination and Auger recombination. Finally, the current density (7) in the gain section is approximated by
j ( z , N ) = I W l d U F d N N 1 l 0 l N d z W l R s ,
where the derivative of the Fermi voltage d U F / d N is taken at a fixed carrier density.
The numerical solution of the coupled-wave Equations (2), (4), (5) and (A1) in the time domain were performed with the software package LDSL-tool developed at the Weierstrass Institute [63].
Several characteristics of the steady states obtained in the numerical simulations of the DBR laser with an upsweep of bias current are shown in Figure 4. The almost-periodic jumps towards the steady-state defined by the adjacent longitudinal optical mode and corresponding changes of the state characteristics are induced by the self-and cross-heating modeled according to (126). While the self-heating of the active section is mainly responsible for the fast shift of the lasing wavelength to longer values and periodic jumps to the shorter-wavelength state, the cross-heating induces a reduced shift of the peak wavelength of the DBR reflectivity and, therefore, the corresponding shift of the mean lasing wavelength, see Figure 4c. The maximum emitted optical power, see Figure 4a, and the minimal (averaged) carrier density, see Figure 4e, within each period coincide well with the maximum field reflection provided by the Bragg grating shown in Figure 4g.
The overall decay of the estimated spectral linewidth, see Figure 4b, reflects its inverse proportionality to the field intensity, see (82). The decay of the linewidth within each period is consistent with the decrease of the modulus of the effective linewidth enhancement factor α H , eff , see Figure 4d, the increase of the chirp reduction factor F shown in Figure 4f, and the behavior of Petermann’s K-factor plotted in Figure 4h. The narrowest linewidth within each period is observed at the long-wavelength flank of the DBR reflectivity just before the transition to the neighboring state. The dependence of F and α H on the deviation of the lasing wavelength from the peak wavelength of the DBR reflectivity confirms earlier findings [16,64,65] that F increases and α H , eff decreases at the long-wavelength flank of the reflection spectrum.

11. Outlook

The theory presented can be easily extended to take into account carrier noise. Including the source F N 0 , we obtain instead of (73) the relation
δ N = τ d F N τ d R st | f | 2 δ | f | 2 .
Substituting this expression into (37) yields a simple modification of the Langevin force
F | f | 2 F | f | 2 + 2 Im ( Φ , β N τ d F N Φ ) ( Φ , β ω Φ ) | f | 2 ,
leading to two additional additive contributions to the linewidth due to carrier noise and the cross-correlation between carrier and intensity noise. The same procedure can be applied to the calculation of the modulation of amplitude and frequency in response to an external modulation of the current injection term in (5).
In case of gain coupling, (21) has to be replaced by
Δ M = β ¯ β κ ¯ κ κ ¯ κ β ¯ β
and (33) has to be supplemented by
κ = κ ( N ) + κ N δ N .
Note that for loss or gain coupling or for higher order gratings there are additional contributions to the imaginary part of β . Furthermore, the rate of stimulated recombination (6) has to be modified [66] and the diffusion coefficient (60) is varied because of non-vanishing F sp * ( z , t ) F sp + ( z , t ) 0 and its complex conjugate [22,47].
Fluctuations of the shape of the power profile can be accounted for by a linearization of (2) and (5) around a steady state, performing a Fourier transformation and solving for the fluctuations δ Ψ , e.g., by the Green’s function method [22] to calculate the noise spectral densities S ξ ( ω ˜ ) from
ξ * ( ω ˜ ) ξ ( ω ˜ ) = 2 π S ξ ( ω ˜ ) δ ( ω ˜ ω ˜ ) ,
where ξ is the variable of interest (e.g., phase or intensity).

12. Summary

We have derived in a self-contained manner the spectral linewidth of edge-emitting multi-section semiconductor lasers starting from the time-dependent coupled-wave equations with a Langevin noise source. For this, we have expanded the forward and backward propagating fields into the longitudinal modes of the open cavity and have obtained very general expressions for the effective linewidth enhancement factor α H and the longitudinal excess factor of spontaneous emission (K-factor), including the effects of nonlinear gain (gain compression) and index (Kerr effect), gain dispersion, longitudinal spatial hole burning for multi-section cavity structures. We have shown that the general linewidth expression contains the effect of linewidth narrowing due to an external cavity by the chirp reduction factor F. Finally, we have investigated the dependence of the population inversion factor on the carrier density and proposed a new analytical formula. Based on the derived expressions, the spectral linewidth as well as α H , K, and F have been calculated for a two-section DBR laser as functions of the injection current and the optical output power, taking into account the thermal detuning between gain and reflector sections. The mode jumps appearing with increasing current are accompanied by sudden rises of the spectral linewidth due to the rise of the modulus of α H and the drop of F.

Funding

This research was partially funded by the Federal Ministry of Education and Research (BMBF) under grant number 50RK1972 and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—The Berlin Mathematics Research Center MATH+ (EXC-2046/1, project ID: 390685689, grant AA2-13).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data can be received from the authors on request.

Acknowledgments

The authors are indebted to Hans-Jürgen Wünsche for a critical reading of the manuscript.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
CWContinuous wave
DBRDistributed Bragg reflector
DFBDistributed feedback
FPFabry–Pérot
PSDpower spectral density
FWHMFull width at half maximum
HWHMHalf width at half maximum

Appendix A. Dispersion Operator

The dispersion operator D ^ ( z , t ) in the coupled-wave equations (2) models the frequency dependence of the gain and associated refractive index. The approach followed here is based on approximating the material gain spectrum locally around the gain peak by a Lorentzian, see Figure A1, and supplementing the time-domain coupled-wave equations (2) with additional dynamical equations for auxiliary polarizations [31,36,37]
t P ( z , t ) = i ω D ( z ) γ D ( z ) P ( z , t ) + γ D ( z ) Ψ ( z , t ) ,
where ω D is the detuning between the peak gain frequency ω p and the reference frequency ω 0 , γ D is the half width at half maximum (HWHM) of the gain curve and
P ( z , t ) = P + ( z , t ) P ( z , t ) .
The corresponding dispersion operator in (2) reads
D ^ ( z , t ) Ψ ( z , t ) = i 2 g D ( z ) Ψ ( z , t ) P ( z , t ) ,
where g D describes the difference between the maximum amplitude gain g ( z , N , Ψ 2 ) , see (3c), and the gain at the reference frequency. The parameters ω D , γ D , and g D are obtained from a fit to a microscopically computed gain spectrum or experimental data.
A Fourier transformation of (A2) yields, together with the frequency domain solution of the polarization Equation (A1), an expression for the dispersion factor (complex Lorentzian)
β D z , Ω = g D ( z ) 2 Ω ω D z i Ω ω D z + γ D z .
The associated effective gain spectrum entering the rate of stimulated emission (67) reads
g eff ( z , N , Ψ 2 , Ω ) = g ( z , N , Ψ 2 ) + 2 Im β D ( z , Ω ) = g ( z , N , Ψ 2 ) g D Ω ω D z 2 Ω ω D z 2 + γ D 2 z ,
which is maximal at Ω = ω D , i.e., if the laser operates at the peak gain frequency.
Figure A1. Illustration of the Lorentzian fitted (effective) gain spectrum and labeling of the dispersion parameters ω D , γ D , and g D . The peak gain is denoted by g eff Ω = ω D = g and g eff Ω = 0 is the effective gain at the reference frequency ω 0 .
Figure A1. Illustration of the Lorentzian fitted (effective) gain spectrum and labeling of the dispersion parameters ω D , γ D , and g D . The peak gain is denoted by g eff Ω = ω D = g and g eff Ω = 0 is the effective gain at the reference frequency ω 0 .
Applsci 11 06004 g0a1
An alternative representation of (A2) is obtained with the help of the Green’s function of the polarization Equation (A1) as
D ^ z , t Ψ z , t = i 2 g D ( z ) Ψ z , t γ D ( z ) 0 t e i ω D ( z ) γ D ( z ) τ Ψ z , t τ d τ ,
where dispersion is induced via the convolution of the memory kernel with the time-delayed field amplitude (we omitted the quickly decaying initial value term here).

Appendix B. Calculation of Derivative of r+

To derive (92) we have to calculate
ω ln ( r + ) = 1 r + ω r + = 1 Φ ω Φ 1 Φ + ω Φ + ,
where we abbreviated / ω ω . The coupled-wave Equation (9) together with their derivatives with respect to frequency can be written as
z Φ + ( z ) = i β ^ ( Ω ) Φ + ( z ) i κ Φ ( z ) , z Φ ( z ) = i β ^ ( Ω ) Φ ( z ) i κ Φ + ( z ) , z ω Φ + ( z ) = i β ^ ( Ω ) ω Φ + ( z ) i κ ω Φ ( z ) i ω β ( Ω ) Φ + ( z ) , z ω Φ ( z ) = i β ^ ( Ω ) ω Φ ( z ) i κ ω Φ + ( z ) i ω β ( Ω ) Φ ( z ) ,
where β ^ ( Ω ) = Δ β ¯ ( z ) + β D ( Ω ) + n g Ω / c and ω β ( Ω ) is defined in (15). By multiplying the left- and right-hand sides of these equations by ω Φ , ω Φ + , Φ , and Φ + , respectively, we get
ω Φ z Φ + = i β ^ ( Ω ) Φ + ω Φ i κ Φ ω Φ , ω Φ + z Φ = i β ^ ( Ω ) Φ ω Φ + i κ Φ + ω Φ + , Φ z ω Φ + = i β ^ ( Ω ) Φ ω Φ + + i κ Φ ω Φ + i ω β ( Ω ) Φ Φ + , Φ + z ω Φ = i β ^ ( Ω ) Φ + ω Φ + i κ Φ + ω Φ + + i ω β ( Ω ) Φ + Φ .
Adding all these equations implies
z Φ + ω Φ Φ ω Φ + = 2 i ω β ( Ω ) Φ + Φ .
An integration over the passive sections corresponding to the coordinate interval [ l , L ] yields
Φ + ω Φ Φ ω Φ + | l L = 2 i l L ω β ( Ω ) Φ + Φ d z .
Due to the reflecting boundary condition, the expression Φ + ω Φ Φ ω Φ + vanishes at z = L and
Φ + ( l ) ω Φ ( l ) + Φ ( l ) ω Φ + ( l ) = 2 i l L ω β ( z , Ω ) Φ + ( z ) Φ ( z ) d z
follows. Dividing by Φ + ( l ) Φ ( l ) finally yields
ω log ( r + ) = 2 i l L ω β ( z , Ω ) Φ + ( z ) Φ ( z ) d z Φ + ( l ) Φ ( l ) .

Appendix C. The Relation to Static Frequency Chirp

Here we justify the naming of the factor F defined in (101). We start from the roundtrip condition
r + ( ω ˜ , g ) r ( ω ˜ , g ) = 1
following from (91) which can be transformed into equations for the modulus
Re ln ( r ( ω ˜ , g ) ) + Re ln ( r + ( ω ˜ ) ) = 0
and the phase
h ( ω ˜ , g ) Im ln ( r ( ω ˜ , g ) ) + Im ln ( r + ( ω ˜ ) ) = 2 π m
with ω ˜ = ω ω 0 and m Z . As a result of the solution of (A13) or, equivalently, (A14) and (A15), we obtain the real-valued frequency Ω = ω ˜ and the gain g of the steady lasing state.
Let the modal index in the active section vary due to a fluctuation of some parameter such as carrier density or temperature. Equation (A15) establishes a relation between the index fluctuation δ Δ n and the frequency fluctuation δ ω ,
δ h = h ω δ ω + h n δ n = 0 ,
where δ n δ Δ n and 1 / n 1 / Δ n . The ω -derivative of h is
h ω = Im ln ( r ) ω + Im ln ( r + ) ω + Im ln ( r ) g g ω .
The ω -derivative of the lasing gain g can be determined from (A14),
g ω = Re ln ( r ) ω + Re ln ( r + ) ω Re ln ( r ) g = 1 l Re ln ( r + ) ω
taking into account
ln ( r ) ω = 2 i l n g c
and
ln ( r ) g = l ( 1 + i α H ) ,
following from (91) and (11), such that
h ω = 2 l n g c + Im ln ( r + ) ω α H Re ln ( r + ) ω
is obtained. Furthermore, (91) and (11) imply
ln ( r ) n = 2 i l 2 π λ 0
and
h n = 4 π l λ 0 .
Therefore,
δ h = 2 l n g c 1 c 2 n g l Im ln ( r + ) ω + α H c 2 n g l Re ln ( r + ) ω δ ω 4 π l λ 0 δ n = 0
and
δ ω = ω 0 δ n n g F 1
is gained. For the solitary laser with r + ( ω ˜ ) taken at a fixed ω ˜ = Ω , the same analysis results in
h ω = 2 l n g c
and
δ ω s = ω 0 δ n n g .
Therefore in a laser with a passive section or an external cavity, fluctuations of the frequency δ ω due to a fluctuation of some parameter in the active section are reduced (or enhanced) by the factor
δ ω δ ω s = 1 F
compared to a the solitary laser, as stated in, e.g., [16,17]. If the active section is not a single FP cavity, the chirp reduction factor can be generalized to [67]
F = 1 + Re ( 1 i α H ) ln ( r + ) ω ln ( r ) ω .

References

  1. Schawlow, A.L.; Townes, C.H. Infrared and optical masers. Phys. Rev. 1958, 112, 1940. [Google Scholar] [CrossRef] [Green Version]
  2. Haken, H. A nonlinear theory of laser noise and coherence. I. Z. Phys. 1964, 181, 96–124. [Google Scholar] [CrossRef]
  3. Lax, M. Classical noise V: Noise in self-sustained oscillators. Phys. Rev. 1967, 160, 290. [Google Scholar] [CrossRef]
  4. Haug, H.; Haken, H. Theory of noise in semiconductor laser emission. Z. Phys. A 1967, 204, 262–275. [Google Scholar] [CrossRef]
  5. Scully, M.O.; Lamb, W.E., Jr. Quantum theory of an optical maser. I. General theory. Phys. Rev. 1967, 159, 208. [Google Scholar] [CrossRef]
  6. Risken, H. Zur Statistik des Laserlichts. Fortschr. Phys. 1968, 16, 261–308. [Google Scholar] [CrossRef]
  7. Petermann, K. Calculated spontaneous emission factor for double-heterostructure injection lasers with gain-induced waveguiding. IEEE J. Quantum Electron. 1979, 15, 566–570. [Google Scholar] [CrossRef]
  8. Siegman, A. Excess spontaneous emission in non-Hermitian optical systems. II. Laser oscillators. Phys. Rev. A 1989, 39, 1264. [Google Scholar] [CrossRef]
  9. Hamel, W.; Woerdman, J. Nonorthogonality of the longitudinal eigenmodes of a laser. Phys. Rev. A 1989, 40, 2785. [Google Scholar] [CrossRef]
  10. Wenzel, H.; Bandelow, U.; Wünsche, H.J.; Rehberg, J. Mechanisms of fast self pulsations in two-section DFB lasers. IEEE J. Quantum Electron. 1996, 32, 69–78. [Google Scholar] [CrossRef]
  11. Özdemir, Ş.K.; Rotter, S.; Nori, F.; Yang, L. Parity–time symmetry and exceptional points in photonics. Nat. Mater. 2019, 18, 783–798. [Google Scholar] [CrossRef]
  12. Henry, C.H. Theory of the linewidth of semiconductor lasers. IEEE J. Quantum Electron. 1982, 18, 259–264. [Google Scholar] [CrossRef]
  13. Amann, M.C. Linewidth enhancement in distributed-feedback semiconductor lasers. Electron. Lett. 1990, 26, 569–571. [Google Scholar] [CrossRef]
  14. Tromborg, B.; Olesen, H.; Pan, X. Theory of linewidth for multielectrode laser diodes with spatially distributed noise sources. IEEE J. Quantum Electron. 1991, 27, 178–192. [Google Scholar] [CrossRef]
  15. Patzak, E.; Sugimura, A.; Saito, S.; Mukai, T.; Olesen, H. Semiconductor laser linewidth in optical feedback configurations. Electron. Lett. 1983, 19, 1026–1027. [Google Scholar] [CrossRef]
  16. Kazarinov, R.; Henry, C.H. The relation of line narrowing and chirp reduction resulting from the coupling of a semiconductor laser to passive resonator. IEEE J. Quantum Electron. 1987, 23, 1401–1409. [Google Scholar] [CrossRef]
  17. Tromborg, B.; Olesen, H.; Pan, X.; Saito, S. Transmission line description of optical feedback and injection locking for Fabry–Perot and DFB lasers. IEEE J. Quantum Electron. 1987, 23, 1875–1889. [Google Scholar] [CrossRef]
  18. Amann, M.C.; Schimpe, R. Excess linewidth broadening in wavelength-tunable laser diodes. Electron. Lett. 1990, 26, 279–280. [Google Scholar] [CrossRef]
  19. Olesen, H.; Tromborg, B.; Lassen, H.; Pan, X. Mode instability and linewidth rebroadening in DFB lasers. Electron. Lett. 1992, 28, 444–446. [Google Scholar] [CrossRef]
  20. Schatz, R. Longitudinal spatial instability in symmetric semiconductor lasers due to spatial hole burning. IEEE J. Quantum Electron. 1992, 28, 1443–1449. [Google Scholar] [CrossRef]
  21. Olesen, H.; Tromborg, B.; Pan, X.; Lassen, H. Stability and dynamic properties of multi-electrode laser diodes using a Green’s function approach. IEEE J. Quantum Electron. 1993, 29, 2282–2301. [Google Scholar] [CrossRef]
  22. Tromborg, B.; Lassen, H.E.; Olesen, H. Traveling wave analysis of semiconductor lasers: Modulation responses, mode stability and quantum mechanical treatment of noise spectra. IEEE J. Quantum Electron. 1994, 30, 939–956. [Google Scholar] [CrossRef]
  23. Lax, M. Classical noise IV: Langevin methods. Rev. Mod. Phys 1966, 38, 541–566. [Google Scholar] [CrossRef]
  24. Lax, M. Quantum noise. IV. Quantum theory of noise sources. Phys. Rev. 1966, 145, 110. [Google Scholar] [CrossRef]
  25. Marani, R.; Lax, M. Spontaneous emission in non-Hermitian optical systems: Distributed-feedback semiconductor lasers. Phys. Rev. A 1995, 52, 2376. [Google Scholar] [CrossRef] [PubMed]
  26. Henry, C.H.; Kazarinov, R.F. Quantum noise in photonics. Rev. Mod. Phys. 1996, 68, 801. [Google Scholar] [CrossRef]
  27. Henry, C.H. Theory of spontaneous emission noise in open resonators and its application to lasers and optical amplifiers. J. Light. Technol. 1986, 4, 288–297. [Google Scholar] [CrossRef]
  28. Lasher, G.; Stern, F. Spontaneous and stimulated recombination radiation in semiconductors. Phys. Rev. 1964, 133, A553. [Google Scholar] [CrossRef]
  29. Agrawal, G.P.; Dutta, N.K. Semiconductor Lasers, 2nd ed.; Kluwer Academic Publishers: Boston, MA, USA, 1993. [Google Scholar] [CrossRef]
  30. Radziunas, M. Traveling wave modeling of nonlinear dynamics in multisection laser diodes. In Handbook of Optoelectronic Device Modeling & Simulation Volume II: Lasers, Modulators, Photodetectors, Solar Cells, and Numerical Methods; CRC Press: Boca Raton, FL, USA, 2017; pp. 153–182. [Google Scholar] [CrossRef]
  31. Sieber, J.; Bandelow, U.; Wenzel, H.; Wolfrum, M.; Wünsche, H.J. Travelling wave equations for semiconductor lasers with gain dispersion. WIAS Preprint 1998, 459, 1–13. [Google Scholar] [CrossRef]
  32. Joyce, W.B. Current-crowded carrier confinement in double-heterostructure lasers. J. Appl. Phys. 1980, 51, 2394–2401. [Google Scholar] [CrossRef]
  33. Zeghuzi, A.; Radziunas, M.; Wenzel, H.; Wünsche, H.J.; Bandelow, U.; Knigge, A. Modeling of current spreading in high-power broad-area lasers and its impact on the lateral far field divergence. Proc. SPIE 2018, 10526, 105261H. [Google Scholar] [CrossRef]
  34. Bandelow, U.; Wenzel, H.; Wünsche, H.J. Influence of inhomogeneous injection on sidemode suppression in strongly coupled DFB semiconductor lasers. Electron. Lett. 1992, 28, 1324–1326. [Google Scholar] [CrossRef]
  35. Lassen, H.E.; Wenzel, H.; Tromborg, B. Influence of series resistance on modulation responses of DFB lasers. Electron. Lett. 1993, 29, 1124–1126. [Google Scholar] [CrossRef]
  36. Ning, C.Z.; Indik, R.; Moloney, J. Effective Bloch equations for semiconductor lasers and amplifiers. IEEE J. Quantum Electron. 1997, 33, 1543–1550. [Google Scholar] [CrossRef]
  37. Bandelow, U.; Radziunas, M.; Sieber, J.; Wolfrum, M. Impact of gain dispersion on the spatio-temporal dynamics of multisection lasers. IEEE J. Quantum Electron. 2001, 37, 183–188. [Google Scholar] [CrossRef]
  38. Krüger, U.; Petermann, K. The semiconductor laser linewidth due to the presence of side modes. IEEE J. Quantum Electron. 1988, 24, 2355–2358. [Google Scholar] [CrossRef]
  39. Pan, X.; Tromborg, B.; Olesen, H. Linewidth rebroadening in DFB lasers due to weak side modes. IEEE Photonic Tech. L. 1991, 3, 112–114. [Google Scholar] [CrossRef]
  40. Gardiner, C. Stochastic Methods: A Handbook for the Natural and Social Sciences; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  41. Agrawal, G.P.; Roy, R. Effect of injection-current fluctuations on the spectral linewidth of semiconductor lasers. Phys. Rev. A 1988, 37, 2495. [Google Scholar] [CrossRef]
  42. Di Domenico, G.; Schilt, S.; Thomann, P. Simple approach to the relation between laser frequency noise and laser line shape. Appl. Opt. 2010, 49, 4801–4807. [Google Scholar] [CrossRef]
  43. Woods, N.D.; Payne, M.C.; Hasnip, P.J. Computing the self-consistent field in Kohn–Sham density functional theory. J. Phys. Condens. Matter 2019, 31, 453001. [Google Scholar] [CrossRef] [Green Version]
  44. Landau, L.D.; Lifshitz, E.M. Course of Theoretical Physics Vol. 5: Statistical Physics (Part 1); Butterworth–Heinemann: Amsterdam, The Netherlands, 1980. [Google Scholar] [CrossRef]
  45. Vahala, K.; Yariv, A. Semiclassical theory of noise in semiconductor lasers-Part II. IEEE J. Quantum Electron. 1983, 19, 1102–1109. [Google Scholar] [CrossRef]
  46. Henry, C.H. Phase noise in semiconductor lasers. J. Light. Technol. 1986, 4, 298–311. [Google Scholar] [CrossRef]
  47. Baets, R.G.; David, K.; Morthier, G. On the distinctive features of gain coupled DFB lasers and DFB lasers with second-order grating. IEEE J. Quantum Electron. 1993, 29, 1792–1798. [Google Scholar] [CrossRef]
  48. Chow, W.W.; Koch, S.W.; Sargent, M.I. Semiconductor-Laser Physics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1994. [Google Scholar] [CrossRef]
  49. Osinski, M.; Buus, J. Linewidth broadening factor in semiconductor lasers—An overview. IEEE J. Quantum Electron. 1987, 23, 9–29. [Google Scholar] [CrossRef]
  50. Champagne, Y.; McCarthy, N. Global excess spontaneous emission factor of semiconductor lasers with axially varying characteristics. IEEE J. Quantum Electron. 1992, 28, 128–135. [Google Scholar] [CrossRef]
  51. Wenzel, H.; Wünsche, H.J. An equation for the amplitudes of the modes in semiconductor lasers. IEEE J. Quantum Electron. 1994, 30, 2073–2080. [Google Scholar] [CrossRef]
  52. Pick, A.; Cerjan, A.; Liu, D.; Rodriguez, A.W.; Stone, A.D.; Chong, Y.D.; Johnson, S.G. Ab initio multimode linewidth theory for arbitrary inhomogeneous laser cavities. Phys. Rev. A 2015, 91, 063806. [Google Scholar] [CrossRef] [Green Version]
  53. Komljenovic, T.; Srinivasan, S.; Norberg, E.; Davenport, M.; Fish, G.; Bowers, J.E. Widely tunable narrow-linewidth monolithically integrated external-cavity semiconductor lasers. IEEE J. Sel. Top. Quantum Electron. 2015, 21, 214–222. [Google Scholar] [CrossRef]
  54. Tran, M.A.; Huang, D.; Bowers, J.E. Tutorial on narrow linewidth tunable semiconductor lasers using Si/III-V heterogeneous integration. APL Photonics 2019, 4, 111101. [Google Scholar] [CrossRef] [Green Version]
  55. Boller, K.J.; van Rees, A.; Fan, Y.; Mak, J.; Lammerink, R.E.; Franken, C.A.; van der Slot, P.J.; Marpaung, D.A.; Fallnich, C.; Epping, J.P.; et al. Hybrid integrated semiconductor lasers with silicon nitride feedback circuits. Photonics 2020, 7, 4. [Google Scholar] [CrossRef] [Green Version]
  56. Coldren, L.A.; Corzine, S.W.; Mashanovitch, M.L. Diode Lasers and Photonic Integrated Circuits; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef]
  57. Spießberger, S.; Schiemangk, M.; Wicht, A.; Wenzel, H.; Erbert, G.; Tränkle, G. DBR laser diodes emitting near 1064 nm with a narrow intrinsic linewidth of 2 kHz. Appl. Phys. B 2011, 104, 813. [Google Scholar] [CrossRef]
  58. Wenzel, H.; Erbert, G.; Enders, P.M. Improved theory of the refractive-index change in quantum-well lasers. IEEE J. Sel. Top. Quantum Electron. 1999, 5, 637–642. [Google Scholar] [CrossRef]
  59. Wenzel, H.; Klehr, A.; Erbert, G.; Sebastian, J.; Tränkle, G.; Pereira, M., Jr. Effect of band gap renormalization on threshold current and efficiency of a distributed Bragg reflector laser. Appl. Phys. Lett. 2000, 76, 2653–2655. [Google Scholar] [CrossRef]
  60. Brox, O.; Wenzel, H.; Fricke, J.; Della Casa, P.; Maaßdorf, A.; Matalla, M.; Wenzel, S.; Wicht, A.; Knigge, A. Novel 1064 nm DBR lasers combining active layer removal and surface gratings. Electron. Lett. 2021. [Google Scholar] [CrossRef]
  61. Schiemangk, M.; Spießberger, S.; Wicht, A.; Erbert, G.; Tränkle, G.; Peters, A. Accurate frequency noise measurement of free-running lasers. Appl. Opt. 2014, 53, 7138–7143. [Google Scholar] [CrossRef]
  62. Batrak, D.V.; Sof’ya, A.B.; Borodaenko, A.V.; Drakin, A.E.; Bogatov, A.P. Simulation of the material gain in quantum-well InGaAs layers used in 1.06-μm heterolasers. Quantum Electron. 2005, 35, 316. [Google Scholar] [CrossRef]
  63. Radziunas, M. Longitudinal Dynamics in multisection Semiconductor Lasers. Available online: https://www.wias-berlin.de/software/ldsl/ (accessed on 22 March 2021).
  64. Patzak, E.; Meissner, P.; Yevick, D. An analysis of the linewidth and spectral behavior of DBR lasers. IEEE J. Quantum Electron. 1985, 21, 1318–1325. [Google Scholar] [CrossRef]
  65. Chaciński, M.; Schatz, R. Impact of losses in the Bragg section on the dynamics of detuned loaded DBR lasers. IEEE J. Quantum Electron. 2010, 46, 1360–1367. [Google Scholar] [CrossRef]
  66. Jonsson, B.; Lowery, A.J.; Olesen, H.; Tromborg, B. Instabilities and nonlinear L-I characteristics in complex-coupled DFB lasers with antiphase gain and index gratings. IEEE J. Quantum Electron. 1996, 32, 839–850. [Google Scholar] [CrossRef]
  67. Tronciu, V.; Wenzel, H.; Wünsche, H.J. Instabilities and Bifurcations of a DFB Laser Frequency-Stabilized by a High-Finesse Resonator. IEEE J. Quantum Electron. 2017, 53, 1–9. [Google Scholar] [CrossRef]
Figure 1. Schematic view of a two-section DBR laser consisting of an electrically biased gain section and a passive Bragg reflector section, exemplifying an edge-emitting multi-section semiconductor laser.
Figure 1. Schematic view of a two-section DBR laser consisting of an electrically biased gain section and a passive Bragg reflector section, exemplifying an edge-emitting multi-section semiconductor laser.
Applsci 11 06004 g001
Figure 2. Schematic illustration of a two-section DBR laser featuring an active gain section and a passive Bragg grating as considered in Section 8.
Figure 2. Schematic illustration of a two-section DBR laser featuring an active gain section and a passive Bragg grating as considered in Section 8.
Applsci 11 06004 g002
Figure 3. (a) Modal gain (blue) and rate of spontaneous emission (green), left axis, and population inversion factor (red), right axis, versus carrier density obtained from a microscopic simulation (symbols) and the analytic models (122), (123) and n sp = r sp / g (solid lines) for a fixed wavelength λ 0 = 1064 nm . The dashed green line is the spontaneous emission computed using a constant population inversion factor n sp = 2 . See the Appendix A for the consideration of the gain at a fixed wavelength within the framework of the present model. (b) Microscopically computed gain spectra for carrier densities ranging from N = 1 · 10 18 cm 3 to 1 · 10 19 cm 3 . The position of the fixed wavelength λ 0 is shown by a blue line, the peak gain position is indicated by a red line.
Figure 3. (a) Modal gain (blue) and rate of spontaneous emission (green), left axis, and population inversion factor (red), right axis, versus carrier density obtained from a microscopic simulation (symbols) and the analytic models (122), (123) and n sp = r sp / g (solid lines) for a fixed wavelength λ 0 = 1064 nm . The dashed green line is the spontaneous emission computed using a constant population inversion factor n sp = 2 . See the Appendix A for the consideration of the gain at a fixed wavelength within the framework of the present model. (b) Microscopically computed gain spectra for carrier densities ranging from N = 1 · 10 18 cm 3 to 1 · 10 19 cm 3 . The position of the fixed wavelength λ 0 is shown by a blue line, the peak gain position is indicated by a red line.
Applsci 11 06004 g003
Figure 4. Simulated characteristics of the DBR laser as functions of the upsweeped injection current I. Dotted and dashed lines indicate mode jumps and maxima of the DBR reflectivity, respectively. (a) Output power P 0 , (b) spectral linewidth Δ ν , (c) deviation from reference wavelength λ λ 0 , (d) effective α -factor α eff , H , (e) stationary carrier density N , (f) chirp reduction factor F, (g) reflection coefficient | r + | 2 , and (h) Petermann K-factor.
Figure 4. Simulated characteristics of the DBR laser as functions of the upsweeped injection current I. Dotted and dashed lines indicate mode jumps and maxima of the DBR reflectivity, respectively. (a) Output power P 0 , (b) spectral linewidth Δ ν , (c) deviation from reference wavelength λ λ 0 , (d) effective α -factor α eff , H , (e) stationary carrier density N , (f) chirp reduction factor F, (g) reflection coefficient | r + | 2 , and (h) Petermann K-factor.
Applsci 11 06004 g004
Table 1. Parameters of the simulated DBR laser.
Table 1. Parameters of the simulated DBR laser.
ParameterSymbolUnitValueFirst Used
reference wavelength λ 0 m 1.064· 10 6
front facet reflectivity | r 0 | 2 0.3 (4)
rear facet reflectivity | r L | 2 0(4)
internal optical loss (both sections) α m−160(3c)
group index (both sections) n g 3.9(2)
built-in index detuning (both sections) Δ n 0 0(126)
Active section
lengthlm 1 · 10 3
differential gain g m2 18.62· 10 22 (122)
α -factor α ˜ H 1(79)
transparency carrier density N tr m−3 1.7· 10 24 (122)
self-heating induced index tuning ν s A−1 9.16· 10 3 (126)
gain clamping density N cl m−3 0.8· 10 24 (122)
index clamping density N cl , i m−3 0.01· 10 24 (125)
gain saturation power P sat W 5.866 (124)
dispersion peak amplitude g D m−150(A3)
dispersion peak frequency detuning ω D rad s−10(A3)
dispersion HWHM γ D rad s−1 83.25· 10 12 (A3)
thickness of active regiondm 4 · 10 6 (5)
width of active regionWm 10 · 10 9 (6)
series resistance R s Ω 1(7)
Fermi voltage derivative d U F d N Vm3 3 · 10 26 (128)
defect recombination coefficientAs−1 4 · 10 8 (127)
bimolecular recombination coefficientBm3 s−1 1 · 10 16 (127)
Auger recombination coefficientCm6 s−1 4 · 10 42 (127)
injection currentIA [ 0 , 300 ] · 10 3 (7)
Passive section
length L l m 3 · 10 3
coupling coefficient κ m−1200(10)
cross-heating induced index tuning ν c A−1 1.83· 10 3 (126)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wenzel, H.; Kantner, M.; Radziunas, M.; Bandelow, U. Semiconductor Laser Linewidth Theory Revisited. Appl. Sci. 2021, 11, 6004. https://doi.org/10.3390/app11136004

AMA Style

Wenzel H, Kantner M, Radziunas M, Bandelow U. Semiconductor Laser Linewidth Theory Revisited. Applied Sciences. 2021; 11(13):6004. https://doi.org/10.3390/app11136004

Chicago/Turabian Style

Wenzel, Hans, Markus Kantner, Mindaugas Radziunas, and Uwe Bandelow. 2021. "Semiconductor Laser Linewidth Theory Revisited" Applied Sciences 11, no. 13: 6004. https://doi.org/10.3390/app11136004

APA Style

Wenzel, H., Kantner, M., Radziunas, M., & Bandelow, U. (2021). Semiconductor Laser Linewidth Theory Revisited. Applied Sciences, 11(13), 6004. https://doi.org/10.3390/app11136004

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