Steady-state ab initio laser theory for N-level lasers

doi:10.1364/OE.20.000474

« Show journal navigation

Steady-state ab initio laser theory for N-level lasers

Alexander Cerjan, Yidong Chong, Li Ge, and A. Douglas Stone »View Author Affiliations

Optics Express, Vol. 20, Issue 1, pp. 474-488 (2012)
http://dx.doi.org/10.1364/OE.20.000474

View Full Text Article

Acrobat PDF (1049 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (24)
Cited By
Figures (6)
Metrics

Abstract

We show that Steady-state Ab initio Laser Theory (SALT) can be applied to find the stationary multimode lasing properties of an N-level laser. This is achieved by mapping the N-level rate equations to an effective two-level model of the type solved by the SALT algorithm. This mapping yields excellent agreement with more computationally demanding N-level time domain solutions for the steady state.

1. Introduction

Semiclassical laser theory, which neglects the quantum fluctuations of the electromagnetic field, is widely used to describe and simulate lasers [1

H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.

W. E. Lamb, “Theory of an optical maser,” Phys. Rev. 134, A1429 (1964). [CrossRef]

]. In principle, it correctly describes the laser thresholds and frequencies, the spatial pattern of the lasing modes, and the laser output power, including all classical non-linear effects, such as spatial hole-burning, gain saturation, and mode and phase locking. Essentially the theory describes Maxwell’s equations in an open cavity, coupled to the non-linear polarization of the gain medium. The gain polarization can be described using either a classical non-linear oscillator model [3

A. E. Siegman, Lasers (University Science Books, 1986).

], or a quantum-mechanical model of N atomic levels in which the polarization and level populations obey the equations of motion of the quantum density matrix. The simplest version of the theory, used widely in textbooks, is the two-level Maxwell-Bloch (MB) model [1

H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.

]; however, most design and characterization simulations of lasers use models with N = 3 or more levels. In addition, most theoretical solutions for the semiclassical laser equations employ a large number of simplifying assumptions in order to make them analytically tractable, most notably neglecting the openness of the cavity and/or treating only simple one-dimensional (1D) or ring cavities, as well as approximating the nonlinear interactions to cubic order. The results are typically not useful for quantitative modeling. Until recently, the only useful way to obtain quantitative results for non-trivial laser structures was to integrate the semiclassical laser equations in space and time. For novel and interesting modern laser structures with non-trivial 2D and 3D cavity geometries, such simulations are at the limits of computational feasibility, making it difficult to study a large parameter space or ensemble of designs.

In the past five years a new approach to finding the stationary solutions of the semiclassical laser equations has been developed, known as Steady state Ab initio Laser Theory (SALT) [6

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]

–9

L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]

]. SALT treats the openness of the cavity exactly, and the multimode non-linear interactions to infinite order (with two approximations, to be discussed below). It is applicable to cavities of arbitrary complexity in 2D and 3D, although we will discuss here only the scalar wave equation coupled to a gain medium. Importantly, this approach eliminates the need to perform a time integration to steady-state, dramatically reducing the computational effort and allowing one to study cavities with high spatial complexity, such as 2D random lasers [8

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

, 10

H. Cao, “Review on the latest developments in random lasers with coherent feedback,” J. Phys. A 38, 10497–10535 (2005). [CrossRef]

], photonic crystal lasers [11

O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science 284, 1819–1821 (1999). [CrossRef] [PubMed]

, 12

S. Chua, Y. D. Chong, A. D. Stone, M Soljac̆ić, and J. Bravo-Abad, “Low-threshold lasing action in photonic crystal slabs enabled by Fano resonances,” Opt. Express 19, 1539 (2011). [CrossRef] [PubMed]

] and chaotic disk microlasers [13

C. Gmachl, F. Capasso, E. E. Narimanov, J. U. Nöckel, A. D. Stone, J. Faist, D. L. Sivco, and A. Y. Cho, “High-power directional emission from microlasers with chaotic resonators,” Science 280, 1556–1564 (1998). [CrossRef] [PubMed]

, 14

Ge Li, Yale PhD thesis, 2010.

SALT was originally formulated to solve the steady-state two-level Maxwell-Bloch equations [6

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]

, 7

H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A 76, 013813 (2007). [CrossRef]

] in the standard slowly-varying envelope approximation (SVEA), and an iterative algorithm was developed [7

H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A 76, 013813 (2007). [CrossRef]

] to solve the resulting SALT equations. Subsequently it was realized that the SVEA afforded no advantage in numerical solutions of the SALT equations, so this approximation was dropped, leading to slightly different SALT equations [14

Ge Li, Yale PhD thesis, 2010.

, 15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

], which were used in all subsequent work [8

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

, 9

L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]

, 12

]. Recently an important generalization of the SALT solution algorithm was developed, improving its performance for laser cavities with complex spatial index variation and inhomogeneous pumping [9

L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]

]. In its present form, SALT only employs two approximations, the stationary inversion approximation (SIA) and the rotating wave approximation (RWA), which are both well satisfied for most lasers of interest. In Ref. [15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

], the results of SALT were compared to the full time-dependent solution of the MB equations for a simple 1D cavity in the multimode regime, well above threshold. Excellent agreement was found in the parameter regime for which the SIA holds. This was, to our knowledge, the first demonstration of a frequency-domain method which agrees with exact time-domain methods for above-threshold multimode lasing.

Previous applications of SALT have focused on two-level gain media. In order for SALT to be a useful modeling tool, it is necessary to demonstrate that it can be applied to N-level lasers. In the current work, we show analytically that the steady-state equations for an N-level laser can be reduced to those for an effective two-level system, and hence solved using the efficient SALT algorithm with essentially the same degree of computational effort. We also explore how this effective two-level system differs from the ordinary two-level laser. Next, we present a numerical comparison between the results of SALT calculations and exact N-level finite-difference time-domain (FDTD) calculations, for the same simple 1D laser studied in [15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

], as well as for a 1D random laser. We note that a similar comparison between SALT and FDTD has been performed for a four-level high-Q single-mode photonic crystal laser in Ref. [12

], with good agreement found. Here, we test SALT’s accuracy in treating the more challenging case of multimode lasing in low-Q and random cavities.

2. Effective two-level systems

2.1. Four level system analysis

We illustrate the approach using the semi-classical laser equations [1

H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.

] for the four-level atomic gain medium shown in Fig. 1:

4 π {\ddot{P}}^{+} = c^{2} \nabla^{2} E^{+} - ε_{c} (r) {\ddot{E}}^{+}

(1)

{\dot{P}}^{+} = - (i ω_{a} + γ_{⊥}) P^{+} + \frac{g^{2}}{i \bar{h}} E^{+} (ρ_{22} - ρ_{11})

(2)

{\dot{ρ}}_{33} = 𝒫 (ρ_{00} - ρ_{33}) - γ_{23} ρ_{33}

(3)

{\dot{ρ}}_{22} = γ_{23} ρ_{33} - γ_{12} ρ_{22} - \frac{1}{i \bar{h}} E^{+} \cdot ({(P^{+})}^{*} - P^{+})

(4)

{\dot{ρ}}_{11} = γ_{12} ρ_{22} - γ_{01} ρ_{11} + \frac{1}{i \bar{h}} E^{+} \cdot ({(P^{+})}^{*} - P^{+})

(5)

{\dot{ρ}}_{00} = γ_{01} ρ_{11} - 𝒫 (ρ_{00} - ρ_{33}) .

(6)

The RWA has already been made, and we have assumed a lasing structure with one or two directions of translational symmetry, so that TM and TE polarizations are conserved and Maxwell’s equations reduce to a scalar Helmholtz equation [16

The equations are written for the TM case, the modifications for TE are straightforward.

]. E⁺ and P⁺ are the positive frequency components of the scalar electric and polarization fields respectively. ρ_ii is the population density of level |i〉, ω_a is the frequency of the gain center, γ_⊥ is the gain width (polarization dephasing rate), g is the dipole matrix element, 𝒫 is the pump rate, and γ_ij is the decay rate from level |j〉 to level |i〉. The four levels are labelled from 0 – 3 in order of increasing energy (Fig. 1). The polarization equation, Eq. (2), is obtained from the four-level density matrix equation of motion, assuming that only the level 2 → 1 transition will be inverted and lase. Often in FDTD calculations a real classical oscillating dipole equation is used to describe the polarization [12

]; in Appendix A, we show that this yields essentially the same results, with an appropriate identification of parameters. In the rate equations, Eqs. (3)–(6), the pump is coherently acting between levels |0〉 and |3〉.

Fig. 1 Schematic of a four-level gain medium.

The polarization equation, Eq. (2), incorporates the inversion D(x,t) ≡ ρ₂₂(x,t) – ρ₁₁(x,t), similar to the polarization equation for a two-level laser. By assuming that the non-lasing populations are stationary, ρ̇₀₀ = ρ̇₃₃ = 0, we can show that D(x,t) obeys

\dot{D} = - γ_{| |}^{'} (D - D_{0}^{'}) - \frac{2}{i \bar{h}} E^{+} \cdot ({(P^{+})}^{*} - P^{+}),

(7)

which is precisely the form of the inversion equation for the two-level medium [1

H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.

, 18

H. Fu and H. Haken, “Multifrequency operations in a short-cavity standing-wave laser,” Phys. Rev. A 43, 2446–2454 (1991). [CrossRef] [PubMed]

]. The parameters D′₀ and γ′_|| serve as an effective equilibrium inversion and inversion relaxation rate respectively, and are given by [14

Ge Li, Yale PhD thesis, 2010.

γ_{| |}^{'} = 2 γ_{12} (1 + \frac{S}{2 + \frac{γ_{01}}{𝒫} + 2 \frac{γ_{01}}{γ_{23}}})

(8)

D_{0}^{'} = \frac{S 𝒫 n}{γ_{01} + (S + 2 + 2 \frac{γ_{01}}{γ_{23}}) 𝒫},

(9)

where S = (γ₀₁ – γ₁₂)/γ₁₂ and n = ∑_iρ_ii is the total density of gain atoms. Equation (9) for the inversion in the absence of laser emission (which here acts as the effective pump parameter) has been discussed by Siegman [3

A. E. Siegman, Lasers (University Science Books, 1986).

], while Eq. (8) for the effective relaxation rate has been derived for a special case by Khanin [19

Y. I. Khanin, Principles of Laser Dynamics (Elsevier, 1995).

]. These expressions have not been used previously to solve the four-level lasing equations in terms of the two-level solutions, as we do here. If we use an incoherent pump, the 𝒫(ρ₀₀ – ρ₃₃) terms in Eqs. (3) and (6) would be replaced by 𝒫ρ₀₀; the effect of this would be the removal of the factor of 2 preceding the ratio γ₀₁/γ₂₃ in the denominators of Eqs. (8) and (9). For typical laser systems, each denominator is dominated by the γ₀₁/𝒫 term, so the difference between coherently and incoherently pumping the system is negligible [17

The observation that coherent and incoherent pumping are nearly equivalent for most systems, is invalid when the coherent pumping is supplied at a similar frequency to the atomic lasing transition and thus interactions between the lasing field and pumping field must be taken into account.

Given the parameters {γ₀₁,γ₁₂,γ₂₃, 𝒫} describing the four-level medium of Fig. 2, we can calculate the effective pump and relaxation rates, and feed those (together with the cavity dielectric function, lasing transition frequency ω_a, and gain linewidth γ_⊥) into the SALT algorithm. That will yield the steady-state lasing properties of this four-level laser.

Fig. 2 Modal intensities as functions of the normalized equilibrium inversion D′₀/D_c (effective pump) in a 1D microcavity edge emitting laser (schematic inset). The cavity is bounded on one side by a perfect mirror and on the other side by air, and has uniform refractive index n = 1.5. Solid lines are results obtained by the time-independent SALT method; open circles are results of FDTD simulations with a coherently pumped four-level medium (Fig. 1); solid triangles are results of FDTD simulations with a coherently pumped six-level medium with a lasing transition between |3〉 and |1〉. Simulation parameters are given in Appendix D. Both the four-level and six-level media are chosen to satisfy SIA. The dephas-ing rate is γ_⊥ = 4.0. The four-level system is in the linear regime described by Eqs. (16) and (17). The six-level system is calculated using the formula in the appendix B, but is in the non-linear regime described by Eqs. (18) and (19). The spectra at D₀/D_c = 0.488, and the gain curve, are shown in the upper left inset.

2.2. Arbitrary number of levels

A general N-level laser can be treated via an analogous procedure. Suppose we have a gain medium with an arbitrary number of levels, N. We assume that there is only a single lasing transition, between two levels which we denote by |u〉 and |l〉. The population of each of the N – 2 non-lasing levels obeys

{\dot{ρ}}_{i i} = \sum_{j} γ_{i j} ρ_{j j} - \sum_{j} γ_{j i} ρ_{i i} + γ_{i u} ρ_{u u} + γ_{i l} ρ_{l l} - (γ_{u i} + γ_{l i}) ρ_{i i},

(10)

where the sums are taken over all of the non-lasing states and γ_ij is the rate at which atoms transition from state |j〉 to |i〉 which is either a decay or pump rate depending on the relative energies of the states. In this way, we can incorporate decay and pump processes between any levels. We again assume that ρ̇_ii = 0 for all non-lasing levels. Then

\sum_{j} [(s_{i} + γ_{u i} + γ_{l i}) δ_{i j} - γ_{i j}] ρ_{i j} = γ_{i u} ρ_{u u} + γ_{i l} ρ_{l l},

(11)

where s_i = ∑_jγ_ji and δ_ij is the Kronecker delta. The term in brackets on the left hand side corresponds to an (N – 2)×(N – 2) matrix, which we denote as R. Upon inverting Eq. (11) and substituting it into the equations of motion for the lasing levels, we obtain

γ_{| |}^{'} = \frac{B_{l} T_{u} - B_{u} T_{l}}{T_{u} + T_{l}}, D_{0}^{'} = \frac{B_{u} + B_{l}}{T_{u} + T_{l}} \frac{n}{γ_{| |}^{'}},

(12)

where

T_{u / l} = 1 + \sum_{i j} {[R^{- 1}]}_{i j} γ_{j, u / l}

(13)

B_{u} = - s_{u} + \sum_{i j} (γ_{u i} - γ_{l i}) {[R^{- 1}]}_{i j} γ_{j u},

(14)

B_{l} = s_{l} + \sum_{i j} (γ_{u i} - γ_{l i}) {[R^{- 1}]}_{i j} γ_{j l} .

(15)

The details of this calculation are given in Appendix B.

3. Physical limits of interest

Returning to the typical four-level case, we take note of two important physical regimes. The first, the linear regime, is γ₂₃ ∼ γ₀₁ ≫ γ₁₂ ≫ 𝒫, for which one recovers the expected behavior that the equilibrium inversion increases linearly with the pump and that γ′_|| is a constant:

γ_{| |}^{'} \approx 2 γ_{12},

(16)

D_{0}^{'} \approx \frac{𝒫}{γ_{12}} n .

(17)

In this case, varying the equilibrium inversion and the pump strength are essentially equivalent.

The second regime of interest, the non-linear regime, is γ₂₃ ∼ γ₀₁ ≫ γ₁₂ ∼ 𝒫, i.e. when the slow decay rate between the lasing levels is on the same order as the pump rate. In this regime, γ′_|| increases with increasing pump and D′₀ saturates with increasing pump:

γ_{| |}^{'} \approx 2 (γ_{12} + 𝒫),

(18)

D_{0}^{'} \approx \frac{1}{1 + \frac{𝒫}{γ_{12}}} (\frac{𝒫}{γ_{12}}) n .

(19)

This regime is also interesting from the viewpoint of SALT. As γ′_|| increases with 𝒫, a laser could satisfy the inequality γ′_|| ≪ γ_⊥ near threshold, leading to stationary inversion and an accurate solution via SALT, but fail to satisfy the inequality as the pump becomes stronger, leading to a decrease in the accuracy of SALT.

For a system with an arbitrary number of levels, the first regime always occurs at sufficiently small pump values; the second regime is obtainable if electrons in the upper lasing level are relatively long-lived compared to electrons in other levels.

4. Brief summary of SALT

For completeness we briefly outline SALT. The E⁺ and P⁺ fields are assumed to obey a multimode ansatz

\begin{array}{l} E^{+} (\vec{r}, t) = \sum_{μ = 1}^{M} Ψ_{μ} (\vec{r}) e^{- i ω_{μ} t}, \\ P^{+} (\vec{r}, t) = \sum_{μ = 1}^{M} p_{μ} (\vec{r}) e^{- i ω_{μ} t}, \end{array}

(20)

where the indices μ = 1,2,···, M label the different lasing modes, and the field and polarization are now explicitly scalar quantities. The total number of modes, M, is not given, but increases in unit steps from zero as we increase the pump strength D₀. The values of D₀ at which each step occurs are the (interacting) modal thresholds, to be determined self-consistently from the theory. The real numbers ω_μ are the lasing frequencies of the modes (henceforth c = 1), which will also be determined self-consistently.

We insert the ansatz Eq. (20) into the two-level laser equations, and employ the stationary inversion approximation (SIA) Ḋ = 0. The result is a set of coupled nonlinear differential equations, which are the fundamental equations of SALT [9

L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]

[\nabla^{2} + (ε_{c} (\vec{r}) + \frac{γ_{⊥} D (\vec{r})}{k_{μ} - k_{a} + i γ_{⊥}}) k_{μ}^{2}] Ψ_{μ} (\vec{r}) = 0,

(21)

D (\vec{r}) = D_{0} (\vec{r}) {[1 + \sum_{ν = 1}^{M} Γ_{ν} {| Ψ_{ν} (\vec{r}) |}^{2}]}^{- 1} .

(22)

Ψ and D are now dimensionless, measured in their natural units

E_{c} = \bar{h} \sqrt{γ_{| |} γ_{⊥}} / (2 g)

and D_c = h̄γ_⊥/(4πg²), and

Γ_{ν} \equiv γ_{⊥}^{2} / (γ_{⊥}^{2} + {(ω_{ν} - ω_{a})}^{2})

is the Lorentzian gain curve evaluated at frequency ω_ν. Note that these equations are time-independent; Eq. (21) is a stationary wave equation for the electric field mode Ψ_μ, with an effective dielectric function consisting of both the “passive” contribution ε_c(r⃗) and an “active” contribution from the gain medium. The latter is frequency-dependent, and has both a real part and a negative (amplifying) imaginary part. It also includes infinite-order nonlinear “hole-burning” modal interactions, seen in the |Ψ_ν|² dependence of Eq. (22). In addition, we make the key requirement that Ψ_μ must be purely outgoing outside the cavity; it is this condition that makes the problem non-Hermitian. It is worth noting that the SIA is not needed until at least two modes are above threshold, so Eq. (22) is exact for single-mode lasing up to and including the second threshold (aside from the well-obeyed RWA).

These equations are solved efficiently by projecting them onto a complete biorthogonal set of purely outgoing states with external wavevectors, k_μ, equal to the lasing frequencies. We refer to these states as the threshold constant flux (TCF) states, because one member of the basis set is always equal to the (non-interacting) threshold lasing mode, leading to very rapid convergence of the basis expansion above threshold [9

L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]

]. The major computational effort in solving the SALT equations by this approach is in calculating the (linear) TCF states. The SALT solutions are obtained for successive values of the pump increasing from the first threshold; at each step, the coefficients of the lasing modes in the TCF basis are obtained using a standard non-linear solver, with the coefficients from the previous step as an initial guess (which is never far from the correct solution). Unlike FDTD, in the current version of SALT one cannot simply directly solve at a fixed pump value, well above threshold. However, even with this limitation, SALT is much more efficient, and provides substantial physical insight, as we will discuss below.

5. Numerical comparison

To perform a well-controlled comparison between SALT and the four-level laser equations, Eqs. (1)–(6), as well as N-level generalizations, we studied 1D microcavity lasers for which the FDTD calculations are tractable and fast enough to generate extensive steady-state data. We first consider the same simple edge emitting uniform-index laser treated in Refs. [6

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]

, 7

H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A 76, 013813 (2007). [CrossRef]

, 15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

], with a perfect mirror at the origin, active region of length L terminating abruptly in air (see schematic, Fig. 2). The simulations were carried out using standard FDTD for the electromagnetic field, and Crank-Nicholson discretization for the polarization and rate equations based on the method of Bidégaray [20

B. Bidégaray, “Time discretizations for Maxwell-Bloch equations,” Numer. Meth. Partial Differential Equations 19, 284–300 (2003). [CrossRef]

] (in which the polarization and inversion are spatially aligned with the electric field but updated at the same time steps as the magnetic field). The reported modal intensities are calculated by Fourier transforming the electric field at the cavity boundary after the simulation has reached steady state (see §6 for a discussion of the steady-state criterion). The lasing transition frequency ω_a is chosen so that n₀k_aL = 60, corresponding to roughly ten wavelengths of radiation within the cavity. Physical quantities are reported in terms of their natural scales, D_c = h̄γ_⊥/(4πg²) and

E_{c} = (\bar{h} / 2 g) \sqrt{γ_{⊥} γ_{| |}^{'}}

. In addition, we take c = h̄ = 1 and measure rates in dimensionless units, i.e. γ_meas = γ_realL/c. We note that the parameters chosen accurately reflect those of real microcavities at optical frequencies [12

]; the complete set of simulation parameters is given in Appendix D.

As shown in Fig. 2, we find close agreement between SALT calculations and FDTD simulations. At a representative pump strength D′₀ = 0.488D_c, the mode intensities produced by SALT differ from those of the four- and six-level FDTD simulations by ∼ 1%, while the frequencies differ by < 0.1%. The difference in mode frequencies between SALT and FDTD also exists at the first lasing threshold, for which an analytical value can be calculated. There, we find that the FDTD simulation has a 0.2% error in the first mode frequency, while SALT has a 0.08% error; this error arises from the spatial discretization of the cavity employed in both approaches [21

For any 1D cavity which is uniformly pumped the TCF states for solving SALT can also be found using a transfer matrix method which does not require discretizing space. We use a more general TCF solver in the calculations presented here which does discretize space.

]. It is worth emphasizing that SALT treats the non-linearity to infinite order; in the earlier work on the Maxwell-Bloch model [15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

] it was shown that for this same cavity the common cubic approximation for the non-linearity fails both quantitatively and qualitatively.

These results demonstrate that so long as the system satisfies SIA, the mapping between systems with an arbitrary number of levels to an effective two-level system is nearly exact, and SALT is able to very accurately determine the steady state properties of the cavity. If two cavities, each with an arbitrary number of levels, have the same effective parameters D′₀ and γ′_||, and otherwise have the same polarization relaxation rate and atomic transition frequency, the cavities are equivalent from the electromagnetic point of view, and will have identical lasing properties.

The six-level simulations shown in Fig. 2 occupy the non-linear parameter regime of Eqs. (18) and (19), i.e. γ′_|| is a linear function and D′₀ a non-linear function of 𝒫. However, the unscaled modal intensity leaving the cavity is still, to leading order, linear in 𝒫. This can be seen by rearranging Eq. (22), inserting the expressions for γ′_|| and D′₀, and noting that at the end of the cavity the inversion is roughly independent of the pump strength. This result is discussed further in Appendix C.

The mapping between the N-level laser and two-level SALT breaks down at large pump strengths, when the condition γ′_|| ≪ γ_⊥ is violated due to the increase of γ′_|| with 𝒫. In Ref. [15

L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]

], following an argument by Haken [18

H. Fu and H. Haken, “Multifrequency operations in a short-cavity standing-wave laser,” Phys. Rev. A 43, 2446–2454 (1991). [CrossRef] [PubMed]

], it was demonstrated that violating this condition causes the SIA for the two-level model to break down. This effect can be seen in the four-level laser data in Fig. 3, where γ′_|| = 0.1, γ_⊥ = 4.0 and accuracy is already lost for the third lasing mode. For the six-level data of Fig. 3, which is in the non-linear parameter regime, SIA is satisfied and SALT agrees with the FDTD simulations for small values of the normalized equilibrium inversion; for larger values of D′₀, the SALT and FDTD results begin to diverge.

Fig. 3 Breakdown of the equivalence between SALT and FDTD when SIA is not valid is shown here in two different ways. Here, modal intensities as a function of the normalized equilibrium inversion D′₀/D_c (effective pump) are shown for a 1D microcavity edge emitting laser with γ_⊥ = 4.0 and n = 1.5. Solid lines again represent results obtained from SALT, while open circles represent FDTD simulations of a simple four-level system with γ′_|| = 0.1. Triangles represent FDTD simulations of a six-level system in the non-linear parameter regime in which γ′_|| ∼ 0.001 for D′₀ ≤ 0.1, and thus satisfying SIA, but γ′_|| ∼ 0.01 for D′₀ ≥ 0.45, and consequently no longer satisfying SIA.

Finally, to demonstrate that the mapping to an effective two-level model works equally well for a complex laser cavity, Fig. 4 shows a comparison between SALT and FDTD simulations for a four-level gain medium in a 1D random dielectric structure. A number of studies have been published on random lasers using such simulations [22

X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70 (2000). [CrossRef] [PubMed]

,23

X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B 69, 104202 (2004). [CrossRef]

]; SALT provides a much more efficient method for such studies, which often require generating a statistical ensemble of lasers. Here, the passive cavity dielectric function contains ∼ 31 layers, alternating randomly between regions with refractive indices n₁ = 1.25 and n₂ = 1. Each random layer was generated according to the formula d_1,2 = 〈d_1,2〉(1+ηζ) where 〈d₁〉 = (1/3)(L/30) and 〈d₂〉 = (2/3)(L/30) are the average thicknesses of the layers, η = 0.9 represents the degree of randomness of the cavity, and ζ ∈ [−1,1] is a randomly generated number. The gain medium was added uniformly to the entire cavity, and the coherent pump was likewise uniform. The transition frequency was chosen such that n₁k_aL = 120, corresponding to roughly 20 wavelengths inside of the cavity. We find only small discrepancies between the SALT and FDTD results, with ∼ 1.1% difference in the modal intensities. These differences did not vary significantly between different realizations of the random laser.

Fig. 4 SALT and FDTD results for a 1D random laser. Modal intensities are plotted against the normalized equilibrium inversion D′₀/D_c (effective pump). Solid lines represent SALT results, and circles represent FDTD simulations for a four-level system with γ′_|| = 0.001. The refractive index distribution of the edge emitting random laser is described in the text. The gain medium has γ_⊥ = 4.0 and is in the regime described by Eqs. (16) and (17). Left inset: log-log plot of the indicated region where three modes turn on in close proximity. Right inset: schematic of the cavity structure.

6. Computational efficiency of SALT for N-level systems

In this section we present a set of benchmarks comparing the computational efficiency of SALT to FDTD. SALT calculations enjoy three main advantages over FDTD simulations of the semi-classical laser equations. First and foremost, SALT directly finds the steady-state solutions, so no time integration is involved, which substantially decreases computational effort. Second, SALT unambiguously determines how many modes are lasing at a given pump, whereas it can be difficult to determine, especially for multimode lasing, when an FDTD simulation has reached the steady-state with all modes that will lase “on”. Third, within SALT, with minimal additional computational effort, it is possible to monitor modes which are below threshold via a modified threshold matrix [8

H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]

], and hence to ascertain if more modes are likely to turn on in some interval of pump.

Structurally SALT has one disadvantage with respect to FDTD. The convergence time in FDTD is determined by the longest time scale in the problem, which is the greater of beating period between consecutive modes and the relaxation oscillation time. These time scales are relatively independent of the number of modes lasing in the cavity, so the efficiency is largely independent of the pump in the multimode regime. This is not the case for SALT, as the computational time increases as N² where N is the number of lasing modes. As SALT was implemented in this and previous works, it automatically calculates the entire lasing fields and spectrum up to a specific pump level, using the results from the lower pump values iteratively to expedite convergence. Thus it is not optimized to produce numerical data at a single given pump level, well above threshold, as one can do easily with FDTD. However studies of the convergence of SALT with an initial guess far away from the final solution have shown that SALT is rather robust and flows to the correct stable solution [14

Ge Li, Yale PhD thesis, 2010.

], and codes optimized for this different type of calculation should be possible As we see in Fig. 5, even a non-optimal implementation of SALT is substantially more efficient than FDTD even when calculating the steady state of a single pump value.

Fig. 5 Comparison of SALT and FDTD run-times. Modal intensities are shown as a function of the run-time for SALT (squares) and four-level FDTD simulations (circles), using the parameters of Fig. 2. FDTD simulations that have not begun to lase are marked as crosses. Plot (a) shows data for D₀/D_c = 0.071, just above the first lasing threshold. SALT determined the steady-state single modal intensity in under three minutes, while the FDTD required ∼ 5000 minutes to reach steady state. Plot (b) shows data for D₀/D_c = 0.486, well above the third lasing threshold. SALT calculated all data up to and including this pump value in under 90 minutes, whereas FDTD required > 500 minutes for the first two modes to reach steady-state, with the third mode intensity (green circles) still fluctuating after 5000 minutes (not shown).

Calculating full modal intensity/frequency curves as a function of the pump strength, such as in Fig. 2, is generally much more efficient using SALT. For example, in order to generate the curves seen in Fig. 2, SALT ran for a little under 2 processor hours. To generate all of the FDTD data for the four level simulations took 267 processor days. If one is attempting to explore a large parameter space of designs or system parameters, SALT may make studies feasible which are simply impractical using FDTD, particularly in more realistic 2D and 3D structures.

As mentioned before, the bulk of the computational effort required for the SALT algorithm, especially in higher dimensions, is in solving for the TCF states. The TCF generalized eigenvalue problem is essentially identical in computational complexity to solving for a portion of the linear resonance spectrum of a cavity, something which is challenging in higher dimensional complex structures but feasible, particularly if one can use the scalar wave equation as in 2D. In practice the number of TCF states needed is often ∼ 10 – 20 and should not exceed ∼ 100 in cases of interest. Advanced finite element codes can be adapted to solve the TCF problem and promise to increase its efficiency in the future. It should be noted that the TCF basis library needs to be generated at a few dozen k values, which does increase the full computation compared to a single resonance calculation. Once one has a TCF basis library, using the SALT algorithm to iterate above threshold does not directly scale with the dimensionality of the system.

Finally, the current implementations of SALT assumes the electric field or magnetic field is perpendicular to the direction of wave propagation, leading to a scalar wave equation; the vectorial generalization of SALT is under investigation and will complicate the TCF calculation to the same extent that they do for more traditional linear Maxwell solvers.

7. Conclusion

We have found that using stationarity conditions on the non-lasing level populations, the rate equations for an N-level laser can be mapped to an effective two-level model, for which the steady-state multimode lasing properties are efficiently solvable using Steady-state Ab initio Laser Theory (SALT). Using this mapping, we found excellent agreement between SALT and FDTD simulations for the modal frequencies, thresholds and above-threshold intensities, in the expected domain of validity. SALT is typically several orders of magnitude more efficient computationally than time domain solution of the laser-rate equations, assuming only steady-state properties are needed.

Appendices

Appendix A: Classical polarization in the laser-rate equations

Throughout this paper, we have used the density matrix equations of motion for the polarization field. However, much of the literature uses the classical oscillating dipole equation, in which the gain atoms are assumed to be dipoles undergoing harmonic oscillations and the quantum density matrix is neglected. This appendix briefly demonstrates that the two models are equivalent, and derives the relevant parameter redefinitions. A more thorough discussion has been given by Boyd [24

R. W. Boyd, Nonlinear Optics (Academic Press, 2008).

The polarization equation used here,

{\dot{P}}^{+} = - (i ω_{a} + γ_{⊥}) P^{+} + \frac{g^{2}}{i \bar{h}} E^{+} (ρ_{u u} - ρ_{l l})

(A1)

is derived from the density matrix equation

{\dot{ρ}}_{u l} = (i ω_{a} + γ_{⊥}) ρ_{u l} - \frac{i}{\bar{h}} g E (ρ_{u u} - ρ_{l l})

(A2)

where |u〉 and |l〉 are the upper and lower lasing states, ρ_ij is the density matrix element ij and g_ul = g_lu = g is the coupling constant of the lasing states to the electric field, which allows for the definition P⁺ = gρ_ul. Alternatively, following the notation of Boyd [24

R. W. Boyd, Nonlinear Optics (Academic Press, 2008).

], we could consider the equation of motion for M = gρ_ul + c.c. = P⁺ + P⁻, which is the expectation value of the dipole moment induced by the applied field, i.e. the classical oscillating dipole field. Thus

\dot{M} = g {\dot{ρ}}_{u l} + c.c. = (i ω_{a} + γ_{⊥}) ρ_{u l} - \frac{i}{\bar{h}} g E D + c.c.,

(A3)

using Eq. (A2), and

\ddot{M} = (- ω_{a}^{2} + 2 i ω_{a} γ_{⊥} + γ_{⊥}^{2}) g ρ_{u l} - \frac{ω_{a}}{\bar{h}} g^{2} E D + c.c. .

(A4)

This can be rewritten as

\ddot{M} + 2 γ_{⊥} \dot{M} + ω_{a}^{2} M = - γ_{⊥}^{2} M - \frac{2 ω_{a} g^{2}}{\bar{h}} E D,

(A5)

where D = ρ_uu – ρ_ll is the inversion density of the lasing states. For

ω_{a}^{2} ≫ γ_{⊥}^{2}

, we can discard the term

γ_{⊥}^{2} M

, resulting in the traditional form of the classical oscillating dipole polarization field equation. This gives the classical coupling constant σ = 2ω_ag²/h̄ [3

A. E. Siegman, Lasers (University Science Books, 1986).

A similar analysis can used to show that the inversion equation can be rewritten as

\dot{D} = - γ_{| |} (D - D_{0}) + \frac{2 E}{\bar{h} ω_{a}} \dot{M},

(A6)

which demonstrates how the change in the inversion is dependent upon the classical polarization field. Therefore, so long as

ω_{a}^{2} ≫ γ_{⊥}^{2}

, a condition that is usually satisfied, these two formulations of the polarization dynamics are equivalent.

Appendix B: Effective two-level parameters from N-level rate equations

This appendix derives the effective equilibrium inversion and relaxation rate for a gain medium with an arbitrary number of levels. We allow decays between any two levels, even if they are not adjacent. We assume there is a single lasing transition, between levels |u〉 and |l〉, which need not be adjacent. The rate equation for an arbitrary non-lasing level in the system is

{\dot{ρ}}_{i i} = \sum_{j} γ_{i j} ρ_{j j} - \sum_{j} γ_{j i} ρ_{i i} + γ_{i u} ρ_{u u} + γ_{i l} ρ_{l l} - (γ_{u i} + γ_{l i}) ρ_{i i},

(B1)

where the summations are taken over all non-lasing levels. Here we do not distinguish between decay rates and pumping rates; γ_ij is simply interpreted as the rate at which level |j〉 transitions into level |i〉, regardless of the energies of those states. If the populations of all the non-lasing transitions are stationary, i.e. ρ̇_ii = 0, then we can rewrite Eq. (B1) as

\sum_{j} R_{i j} ρ_{j j} = γ_{i u} ρ_{u u} + γ_{i l} ρ_{l l},

(B2)

R_{i j} \equiv (s_{i} + γ_{u i} + γ_{l i}) δ_{i j} - γ_{i j} .

(B3)

Here, s_i ≡ ∑_jγ_ji and δ_ij is the Kronecker delta. Inverting Eq. (B2) gives

ρ_{i i} = \sum_{j} {[R^{- 1}]}_{i j} (γ_{j u} ρ_{u u} + γ_{j l} ρ_{l l}) .

(B4)

Hence, we can express the total number density of gain atoms as

n = \sum_{i} ρ_{i i} + ρ_{u u} + ρ_{l l}

(B5)

= T_{u} ρ_{u u} + T_{l} ρ_{l l}

(B6)

where

T_{u} = 1 + \sum_{i j} {[R^{- 1}]}_{i j} γ_{j u},

(B7)

T_{l} = 1 + \sum_{i j} {[R^{- 1}]}_{i j} γ_{j l} .

(B8)

Noting that D = ρ_uu – ρ_ll, we can write the populations of the lasing states as

ρ_{u u} = \frac{n + T_{l} D}{T_{l} + T_{u}},

(B9)

ρ_{l l} = \frac{n - T_{u} D}{T_{l} + T_{u}} .

(B10)

From the equations of motion for the lasing levels, we have the inversion equation

\dot{D} = {\dot{ρ}}_{u u} - {\dot{ρ}}_{l l} = \sum_{i} (γ_{u i} - γ_{l i}) ρ_{i i} - s_{u} ρ_{u u} + s_{l} ρ_{l l} - \frac{2}{i \bar{h}} E^{+} \cdot ({(P^{+})}^{*} - P^{+}),

(B11)

where s_u = ∑_jγ_ju + γ_lu and s_l is defined similarly. Inserting Eq. (B4) into this equation gives

\dot{D} = B_{u} ρ_{u, u} + B_{l} ρ_{l, l} - \frac{2}{i \bar{h}} E^{+} \cdot ({(P^{+})}^{*} - P^{+}),

(B12)

where

B_{u} \equiv - s_{u} + \sum_{i j} (γ_{u i} - γ_{l i}) {[R^{- 1}]}_{i j} γ_{j u},

(B13)

B_{l} \equiv s_{l} + \sum_{i j} (γ_{u i} - γ_{l i}) {[R^{- 1}]}_{i j} γ_{j l} .

(B14)

Plugging in Eqs. (B9) and (B10) now yields the inversion equation in the desired form Eq. (7), with

γ_{| |}^{'} = \frac{B_{l} T_{u} - B_{u} T_{l}}{T_{u} + T_{l}},

(B15)

D_{0}^{'} = \frac{B_{u} + B_{l}}{B_{l} T_{u} - B_{u} T_{l}} n .

(B16)

Appendix C: Inversion as a function of the pump

In this appendix we discuss the surprising result that the unscaled modal intensities of the six-level simulations discussed in Fig. 2 as a function of the pump are approximately linear, as shown in Fig. 6, even though these six-level simulations are in the non-linear parameter regime. In simple treatments of lasers [3

A. E. Siegman, Lasers (University Science Books, 1986).

] the inversion is assumed to be clamped after the first lasing threshold. However, a more detailed analysis lead to the inclusion of spatial hole-burning effects which forces the inversion in the cavity to change beyond the first lasing threshold in a non-linear manner. This result then, that the unscaled modal intensity is linear in the pump rate, can be understood from the observation that at the end of the cavity, r⃗ = L, all of the lasing modes have a maximum in their fields, and thus the inversion is effectively clamped at this point beyond the first lasing threshold, which can be seen in Fig. 6(b).

Fig. 6 (a) Unscaled modal intensity of the six-level simulations from Fig. 2 as a function of the pump. A cross sectional area of 1m² is assumed to calculate the power. (b) Inversion as a function of the pump at the cavity boundary. Dashed lines in plots a and b correspond to the pump values shown in plot c. (c) Inversion as a function of position within the cavity for three different pump values, cyan corresponds with 𝒫 = 3.75 × 10⁸s⁻¹, magenta with 𝒫 = 1.65 × 10⁹s⁻¹, and orange with 𝒫 = 4.85 × 10⁹s⁻¹ to show the evolution of the inversion within the cavity as a function of the pump strength.

To understand this formally, we begin by rewriting Eq. (22) and removing the scaling factors E_c and D_c gives

\frac{2 g^{2}}{{\bar{h}}^{2}} \sum_{ν = 1}^{N} Γ_{ν} {| Ψ_{ν} (\vec{r}) |}^{2} = γ_{| |}^{'} (\frac{D_{0}^{'}}{D (\vec{r})} - 1) .

(C1)

Substituting in Eqs. (18) and (19), which are valid for this simulation, gives

\frac{g^{2}}{{\bar{h}}^{2}} \sum_{ν = 1}^{N} Γ_{ν} {| Ψ_{ν} (\vec{r}) |}^{2} = 𝒫 (\frac{N}{D (\vec{r})} - 1) - γ_{12} .

(C2)

The inversion D is a function of both position and the pump. However, for r⃗ = L corresponding to the cavity edge, D should be mostly independent of the pump, as at this location every mode is at its maximum intensity and the effect of spatial hole-burning is most pronounced. The FDTD simulation results, shown in Fig. 6, demonstrate that D indeed varies very weakly with 𝒫 at the cavity edge. Since the left hand side of Eq. (C2) evaluated at r = L is proportional to the total output power, and this must hold for any number, N, of lasing modes, pump-independence of D(r) implies a linear behavior of the unscaled output power of each mode with pump, as observed in Fig. 6(a).

Appendix D: Simulation constants

In this appendix we list the parameters used in each of the FDTD simulations described above. These constants are matrices, with γ_ij denoting the decay rate from |j〉 to |i〉. These values are given in their dimensionless form, i.e. γ_meas = γ_realL/c. Unlisted entries are zero. We also note that throughout this paper |0〉 denotes the ground state, so these matrices are 0 indexed.

For the four-level simulations in Fig. 2,

γ_{4 lv, fig 2} = (\begin{matrix} \cdot & 0.8 & \cdot & \cdot \\ \cdot & \cdot & 5 \times 10^{- 4} & \cdot \\ \cdot & \cdot & \cdot & 0.8 \\ \cdot & \cdot & \cdot & \cdot \end{matrix}) .

(D1)

The dipole matrix element is g = 2.3 · 10⁻¹²m^3/2, and the number of gain atoms is n = 5 · 10²³m⁻³. The pump 𝒫 was varied between 3 · 10⁻⁶ and 3 · 10⁻⁵.

Thus, for an optical wavelength of λ = 628nm, the requirement in Fig. 2 that n₀kL = 60 means that L = 4μm. Using this length, the decay rates can be converted to their unit-full values as γ_⊥ = 3 · 10¹⁴s⁻¹, γ₂₃ = γ₀₁ = 6 · 10¹³s⁻¹, γ₁₂ = 3.75 · 10¹⁰s⁻¹, and the pump at threshold is 𝒫 = 3 · 10⁸s⁻¹. Similarly, the dipole matrix element also acquires units of inverse time, and can be expressed as g²/h̄ = 3.98 · 10⁻⁹m³/s, which corresponds to a coupling constant in the classical oscillating dipole picture of σ = 10⁻⁴C²/kg. These constants can be seen to be similar to those used in other studies of optical microcavities [4

A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]

, 12

For the six-level simulations in Fig. 2,

γ_{6 lv, fig 2} = (\begin{matrix} \cdot & 0.8 & 10^{- 5} & 10^{- 5} & 10^{- 5} & 10^{- 5} \\ \cdot & \cdot & 0.8 & 10^{- 5} & 10^{- 5} & 10^{- 4} \\ \cdot & \cdot & \cdot & 5 \times 10^{- 5} & 10^{- 5} & 10^{- 5} \\ \cdot & \cdot & \cdot & \cdot & 0.8 & 10^{- 5} \\ \cdot & \cdot & \cdot & \cdot & \cdot & 0.8 \\ \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \end{matrix}) .

(D2)

Furthermore, γ₁₅ = 10⁻⁴, and the lasing transition is between levels |3〉 and |1〉 (where the ground state is again |0〉 and the states are numbered in order of increasing energy).

For the four-level simulations in Fig. 3,

γ_{4 lv, fig 3} = (\begin{matrix} \cdot & 0.8 & \cdot & \cdot \\ \cdot & \cdot & 5 \times 10^{- 2} & \cdot \\ \cdot & \cdot & \cdot & 0.8 \\ \cdot & \cdot & \cdot & \cdot \end{matrix}) .

(D3)

For the six-level simulations in Fig. 3,

γ_{6 lv, fig 3} = (\begin{matrix} \cdot & 0.8 & 10^{- 5} & 10^{- 5} & 10^{- 5} & 10^{- 5} \\ \cdot & \cdot & 0.8 & 10^{- 5} & 10^{- 5} & 10^{- 4} \\ \cdot & \cdot & \cdot & 5 \times 10^{- 5} & 10^{- 5} & 10^{- 5} \\ \cdot & \cdot & \cdot & \cdot & 0.8 & 10^{- 5} \\ \cdot & \cdot & \cdot & \cdot & \cdot & 0.8 \\ \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \end{matrix}) .

(D4)

The four-level simulations of the random cavity in Fig. 4 used the same parameters as the four-level simulations in Fig. 2.

Acknowledgments

We thank Robert J. Tandy, Hui Cao, and Peter Bermel for helpful discussions. This work was partially supported by NSF grant No. DMR-0908437.

References and links

1.	H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.
2.	W. E. Lamb, “Theory of an optical maser,” Phys. Rev. 134, A1429 (1964). [CrossRef]
3.	A. E. Siegman, Lasers (University Science Books, 1986).
4.	A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]
5.	K. S. Yee, “Numerical solution of the initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans Antennas Propag . 14, 302–307 (1966). [CrossRef]
6.	H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]
7.	H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A 76, 013813 (2007). [CrossRef]
8.	H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]
9.	L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A 82, 063824 (2010). [CrossRef]
10.	H. Cao, “Review on the latest developments in random lasers with coherent feedback,” J. Phys. A 38, 10497–10535 (2005). [CrossRef]
11.	O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science 284, 1819–1821 (1999). [CrossRef] [PubMed]
12.	S. Chua, Y. D. Chong, A. D. Stone, M Soljac̆ić, and J. Bravo-Abad, “Low-threshold lasing action in photonic crystal slabs enabled by Fano resonances,” Opt. Express 19, 1539 (2011). [CrossRef] [PubMed]
13.	C. Gmachl, F. Capasso, E. E. Narimanov, J. U. Nöckel, A. D. Stone, J. Faist, D. L. Sivco, and A. Y. Cho, “High-power directional emission from microlasers with chaotic resonators,” Science 280, 1556–1564 (1998). [CrossRef] [PubMed]
14.	Ge Li, Yale PhD thesis, 2010.
15.	L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]
16.	The equations are written for the TM case, the modifications for TE are straightforward.
17.	The observation that coherent and incoherent pumping are nearly equivalent for most systems, is invalid when the coherent pumping is supplied at a similar frequency to the atomic lasing transition and thus interactions between the lasing field and pumping field must be taken into account.
18.	H. Fu and H. Haken, “Multifrequency operations in a short-cavity standing-wave laser,” Phys. Rev. A 43, 2446–2454 (1991). [CrossRef] [PubMed]
19.	Y. I. Khanin, Principles of Laser Dynamics (Elsevier, 1995).
20.	B. Bidégaray, “Time discretizations for Maxwell-Bloch equations,” Numer. Meth. Partial Differential Equations 19, 284–300 (2003). [CrossRef]
21.	For any 1D cavity which is uniformly pumped the TCF states for solving SALT can also be found using a transfer matrix method which does not require discretizing space. We use a more general TCF solver in the calculations presented here which does discretize space.
22.	X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70 (2000). [CrossRef] [PubMed]
23.	X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B 69, 104202 (2004). [CrossRef]
24.	R. W. Boyd, Nonlinear Optics (Academic Press, 2008).

OCIS Codes
(140.3430) Lasers and laser optics : Laser theory
(140.3945) Lasers and laser optics : Microcavities

ToC Category:
Lasers and Laser Optics

History
Original Manuscript: September 27, 2011
Revised Manuscript: November 21, 2011
Manuscript Accepted: November 22, 2011
Published: December 21, 2011

Citation
Alexander Cerjan, Yidong Chong, Li Ge, and A. Douglas Stone, "Steady-state ab initio laser theory for N-level lasers," Opt. Express 20, 474-488 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-1-474

Sort: Author | Year | Journal | Reset

References

H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.
W. E. Lamb, “Theory of an optical maser,” Phys. Rev.134, A1429 (1964). [CrossRef]
A. E. Siegman, Lasers (University Science Books, 1986).
A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag.46, 334–340 (1998). [CrossRef]
K. S. Yee, “Numerical solution of the initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans Antennas Propag. 14, 302–307 (1966). [CrossRef]
H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A74, 043822 (2006). [CrossRef]
H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A76, 013813 (2007). [CrossRef]
H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science320, 643–646 (2008). [CrossRef] [PubMed]
L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A82, 063824 (2010). [CrossRef]
H. Cao, “Review on the latest developments in random lasers with coherent feedback,” J. Phys. A38, 10497–10535 (2005). [CrossRef]
O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science284, 1819–1821 (1999). [CrossRef] [PubMed]
S. Chua, Y. D. Chong, A. D. Stone, M Soljac̆ić, and J. Bravo-Abad, “Low-threshold lasing action in photonic crystal slabs enabled by Fano resonances,” Opt. Express19, 1539 (2011). [CrossRef] [PubMed]
C. Gmachl, F. Capasso, E. E. Narimanov, J. U. Nöckel, A. D. Stone, J. Faist, D. L. Sivco, and A. Y. Cho, “High-power directional emission from microlasers with chaotic resonators,” Science280, 1556–1564 (1998). [CrossRef] [PubMed]
Ge Li, Yale PhD thesis, 2010.
L. Ge, R. J. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative verification of ab initio self-consistent laser theory,” Opt. Express16, 16895 (2008). [CrossRef] [PubMed]
The equations are written for the TM case, the modifications for TE are straightforward.
The observation that coherent and incoherent pumping are nearly equivalent for most systems, is invalid when the coherent pumping is supplied at a similar frequency to the atomic lasing transition and thus interactions between the lasing field and pumping field must be taken into account.
H. Fu and H. Haken, “Multifrequency operations in a short-cavity standing-wave laser,” Phys. Rev. A43, 2446–2454 (1991). [CrossRef] [PubMed]
Y. I. Khanin, Principles of Laser Dynamics (Elsevier, 1995).
B. Bidégaray, “Time discretizations for Maxwell-Bloch equations,” Numer. Meth. Partial Differential Equations19, 284–300 (2003). [CrossRef]
For any 1D cavity which is uniformly pumped the TCF states for solving SALT can also be found using a transfer matrix method which does not require discretizing space. We use a more general TCF solver in the calculations presented here which does discretize space.
X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett.85, 70 (2000). [CrossRef] [PubMed]
X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B69, 104202 (2004). [CrossRef]
R. W. Boyd, Nonlinear Optics (Academic Press, 2008).

IEEE Trans Antennas Propag

K. S. Yee, “Numerical solution of the initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans Antennas Propag. 14, 302–307 (1966). [CrossRef]

IEEE Trans. Antennas Propag.

A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag.46, 334–340 (1998). [CrossRef]

J. Phys. A

H. Cao, “Review on the latest developments in random lasers with coherent feedback,” J. Phys. A38, 10497–10535 (2005). [CrossRef]

Numer. Meth. Partial Differential Equations

B. Bidégaray, “Time discretizations for Maxwell-Bloch equations,” Numer. Meth. Partial Differential Equations19, 284–300 (2003). [CrossRef]

Opt. Express

Phys. Rev.

W. E. Lamb, “Theory of an optical maser,” Phys. Rev.134, A1429 (1964). [CrossRef]

Phys. Rev. A

H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A74, 043822 (2006). [CrossRef]
H. E. Türeci, A. D. Stone, and L. Ge, “Theory of the spatial structure of nonlinear lasing modes,” Phys. Rev. A76, 013813 (2007). [CrossRef]
H. Fu and H. Haken, “Multifrequency operations in a short-cavity standing-wave laser,” Phys. Rev. A43, 2446–2454 (1991). [CrossRef] [PubMed]
L. Ge, Y. D. Chong, and A. D. Stone, “Steady-state ab initio laser theory: generalizations and analytic results,” Phys. Rev. A82, 063824 (2010). [CrossRef]

Phys. Rev. B

X. Jiang, S. Feng, C. M. Soukoulis, J. Zi, J. D. Joannopoulos, and H. Cao, “Coupling, competition, and stability of modes in random lasers,” Phys. Rev. B69, 104202 (2004). [CrossRef]

Phys. Rev. Lett.

X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett.85, 70 (2000). [CrossRef] [PubMed]

Science

C. Gmachl, F. Capasso, E. E. Narimanov, J. U. Nöckel, A. D. Stone, J. Faist, D. L. Sivco, and A. Y. Cho, “High-power directional emission from microlasers with chaotic resonators,” Science280, 1556–1564 (1998). [CrossRef] [PubMed]
O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science284, 1819–1821 (1999). [CrossRef] [PubMed]
H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science320, 643–646 (2008). [CrossRef] [PubMed]

Other

A. E. Siegman, Lasers (University Science Books, 1986).
H. Haken, Light: Laser Dynamics (North-Holland Phys. Publishing, 1985), Vol. 2.
Ge Li, Yale PhD thesis, 2010.
Y. I. Khanin, Principles of Laser Dynamics (Elsevier, 1995).
The equations are written for the TM case, the modifications for TE are straightforward.
The observation that coherent and incoherent pumping are nearly equivalent for most systems, is invalid when the coherent pumping is supplied at a similar frequency to the atomic lasing transition and thus interactions between the lasing field and pumping field must be taken into account.
For any 1D cavity which is uniformly pumped the TCF states for solving SALT can also be found using a transfer matrix method which does not require discretizing space. We use a more general TCF solver in the calculations presented here which does discretize space.
R. W. Boyd, Nonlinear Optics (Academic Press, 2008).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper