Spatio-temporal theory of lasing action in optically-pumped rotationally excited molecular gases

doi:10.1364/OE.19.007513

Spatio-temporal theory of lasing action in optically-pumped rotationally excited molecular gases

Song-Liang Chua, Christine A. Caccamise, Dane J. Phillips, John D. Joannopoulos, Marin Soljačić, Henry O. Everitt, and Jorge Bravo-Abad »View Author Affiliations

Optics Express, Vol. 19, Issue 8, pp. 7513-7529 (2011)
http://dx.doi.org/10.1364/OE.19.007513

View Full Text Article

Acrobat PDF (5770 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Theoretical framework
3. Results and discussion
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

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

Abstract

We investigate laser emission from optically-pumped rotationally excited molecular gases confined in a metallic cavity. To this end, we have developed a theoretical framework able to accurately describe, both in the spatial and temporal domains, the molecular collisional and diffusion processes characterizing the operation of this class of lasers. The effect on the main lasing features of the spatial variation of the electric field intensity and the ohmic losses associated to each cavity mode are also included in our analysis. Our simulations show that, for the exemplary case of methyl fluoride gas confined in a cylindrical copper cavity, the region of maximum population inversion is located near the cavity walls. Based on this fact, our calculations show that the lowest lasing threshold intensity corresponds to the cavity mode that, while maximizing the spatial overlap between the corresponding population inversion and electric-field intensity distributions, simultaneously minimizes the absorption losses occurring at the cavity walls. The dependence of the lasing threshold intensity on both the gas pressure and the cavity radius is also analyzed and compared with experiment. We find that as the cavity size is varied, the interplay between the overall gain of the system and the corresponding ohmic losses allows for the existence of an optimal cavity radius which minimizes the intensity threshold for a large range of gas pressures. The theoretical analysis presented in this work expands the current understanding of lasing action in optically-pumped far-infrared lasers and, thus, could contribute to the development of a new class of compact far-infrared and terahertz sources able to operate efficiently at room temperature.

1. Introduction

It has been forty years since the pioneering work of Chang and Bridges [1

T. Y. Chang and T. J. Bridges, “Laser actions at 452, 496, and 541 μm in optically pumped CH₃F,” Opt. Commun. 1, 423–426 (1970). [CrossRef]

] demonstrated generation of submillimeter radiation via stimulated emission from optically-pumped methyl fluoride (CH₃F) gas. Since then, a rich spectrum of laser lines through the 50 μm – 2 mm wave-length region have been obtained from a variety of different molecular gases [2

T. K. Plant, L. A. Newman, E. J. Danielewicz, T. A. DeTemple, and P. D. Coleman, “High power optically pumped far infrared lasers,” IEEE Trans. Microwave Theory Tech. MT22, 988–990 (1980).

–10

R. L. Crownover, H. O. Everitt, D. D. Skatrud, and F. C. DeLucia, “Frequency stability and reproductibility of optically pumped far-infrared lasers,” Appl. Phys. Lett. 57, 2882–2884 (1990). [CrossRef]

]. These developments have made this class of lasers—usually referred to as optically-pumped far infrared (OPFIR) lasers—a central resource for many scientific and technological applications. The physical mechanism underlying the operation of conventional OPFIR lasers has also been extensively studied in past decades [2

T. K. Plant, L. A. Newman, E. J. Danielewicz, T. A. DeTemple, and P. D. Coleman, “High power optically pumped far infrared lasers,” IEEE Trans. Microwave Theory Tech. MT22, 988–990 (1980).

–15

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.

]. In essence, the gain in these systems derives from an optically-driven population inversion between two rotational states in an excited vibrational level of a gas molecule; the laser emission frequency is determined by the frequency of the transition between these two rotational states. The feedback mechanism enabling laser action is often provided by designing the cavity, in which the molecular gas is confined, to act as a one-dimensional Fabry-Perot resonator for the emitted laser frequency.

However, although research on conventional OPFIR lasers is generally considered as a mature field, some central challenges remain for the ultimate application of this technology, such as low efficiency, poor tunability, or limited capability of integration into a compact platform. More than two decades ago it was shown [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

] that some of the aforementioned drawbacks can be overcome by considering the active role that high-energy vibrational states play in the laser operation (previous approaches assumed that the populations of those vibrational states were negligible). These findings demonstrated the feasibility of developing OPFIR lasers able to operate at large pressures (larger than the pressure cut-off determined by the so-called vibrational bottleneck [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).

]) and featuring volumes up to four orders of magnitude smaller than those characterizing typical OPFIR lasers. The advent of quantum cascade lasers to replace CO₂ pump lasers may combine to produce truly compact submillimeter-wave laser sources in the near future.

Surprisingly, despite their great promise, up to date there has been limited interest in further exploring the potential of these discoveries [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

] for the development of a novel class of compact far-infrared and terahertz (THz) sources able to operate efficiently at room temperature. In order to do that, it is critical to develop a complete theoretical understanding of the physical mechanisms underlying lasing action in the considered systems. The purpose of this paper is to present such a theoretical analysis. Specifically, we introduce a theoretical framework which generalizes previous works in this area to allow for a realistic description of the spatio-temporal dynamics characterizing the molecular collisional and diffusion processes occurring in the considered OPFIR lasers. Importantly, the effect on the main laser features of both the electric-field (E-field) intensity distribution inside the metallic cell and the ohmic losses associated with the penetration of the E-field inside the metallic walls are also included in our theoretical approach. As we show in this work, both effects are particularly important for a realistic description of the operation of compact OPFIR lasers based on a metallic cavity which waveguides the emitted laser frequency.

Using this theoretical framework, we perform a detailed numerical analysis of laser emission from optically pumped methyl fluoride (CH₃F) gas confined in a cylindrical copper cavity, and find that the model compares well against experimental results. Our calculations show that, for this exemplary system, the region of maximum population inversion concentrates near the cavity walls. The outcome of these calculations also implies that, as a consequence of this population inversion distribution, the lowest threshold intensity corresponds to the cavity mode whose E-field intensity profile presents an optimal spatial overlap with the population inversion distribution, while at the same time, the ohmic losses associated to that cavity mode are minimized. In the case of the considered cylindrical cavity both features are satisfied by the TE₀₁ mode. We also analyze how the threshold intensity depends on the gas pressure and the cavity radius. The results from these simulations are explained in terms of the corresponding variation of the overall gain of the system and the ohmic losses associated with the lowest threshold mode. Remarkably, we predict that as the cavity size is varied, the interplay between both factors yields an optimal value of the cavity radius for which the threshold intensity is minimized for a large range of gas pressures.

The paper is organized as follows. Section 2 discusses the theoretical model used throughout this paper. Section 3 presents a detailed analysis of lasing action in the considered molecular gas lasers. Finally, in Section 4 we provide a set of conclusions of this work.

2. Theoretical framework

We start by introducing the theoretical model used to describe the molecular dynamics of the considered optically-pumped polyatomic molecular gases. Fig. 1 illustrates the energy level diagram and the different energy transfer mechanisms that form the basic structure of our model. In the following, we describe each of the processes sketched in Fig. 1 and show how to quantify them using experimentally measured parameters. Before starting this analysis, we point out that the theoretical approach presented here is an extension of previous studies on methyl fluoride CH₃F [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

, 7

R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).

] in order to account for the spatial variations and diffusion of the collisional processes. We also point out that although for demonstration the model is applied to the case of methyl-halide gases, our theoretical analysis can be easily generalized to most OPFIR systems, both pulsed and CW, and operating over a wide range of experimental conditions utilizing an appropriate gas that satisfied the desired pumping and lasing transition frequencies [15

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.

Fig. 1 Schematic diagram of the general model used to describe the dynamics of a OPFIR molecular gas. Details on each process labeled are described in the text.

The experimentally-motivated thermal pool concept has been previously introduced in OPFIR laser models [3

W. H. Matteson and F. C. De Lucia, “Millimeter wave spectroscopic studies of collision-induced energy transfer processes in the ¹³CH₃F laser,” IEEE J. Quantum Electron. 19, 1284–1293 (1983). [CrossRef]

, 7

R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).

, 8

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]

, 15

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.

]. In the case of molecules, a thermal pool is observed to be a collection of rotational states that are in thermal equilibrium, and thus, are populated according to a Boltzmann distribution at the operating temperature. Such equilibrium may be attained, for instance, due to the fast relaxation rates among the levels. However, the total population of these pools varies through the flow of energy governed by the collisional and diffusive processes shown in Fig. 1. This fortuitous consolidation of numerous rotational levels into a small number of pools makes numerical simulation tractable, and more importantly, makes rotational energy transfer studies possible in the first place. In the particular case of methyl halides, there exist two thermal pools for each vibrational state, corresponding to the two symmetry types (A – and E – type), but without loss of generality they will be combined here to describe CW OPFIR laser operation at low pressure [8

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]

]. In addition, all other higher vibrational states that are not connected by the pump are also modeled as pools. This inclusion is important at high pressure and pump intensity since it provides an additional channel through which molecules from the excited vibrational state may be transferred out, and hence, reduces the absorption on the corresponding FIR lasing transition (represented as the pink arrow in the excited state of Fig. 1).

At the relatively low pressure regime considered in this paper, molecules are returned to the ground thermal pool from the pools of the excited and higher excited states, shown as green arrows in Fig. 1, primarily by colliding with the walls of the cavity in which the molecular gas is confined. Unlike molecule-molecule de-excitation processes, these involve diffusion through the gas towards the wall where their rates of occurrences decrease with pressure, and can be defined as k _wall = vλ _MFP/R ² [16

L. E. Reichl, A Modern Course in Statistical Physics (John Wiley & Sons Inc., 1998).

]. λ _MFP is the pressure-dependent mean free path and v is the average relative velocity between molecules. At very low pressure (or equivalently, within λ _MFP from the cavity walls), the transport mechanism changes from diffusion to ballistic and the rate becomes pressure-independent so that k _direct = v _abs/R, with v _abs defined to be the average absolute velocity of a molecule. It is clear from these expressions that the wall rates are very different across the cavity, and we explicitly account for this spatial variation in our numerical simulations. On the other hand, the vibrational relaxation between the excited and higher excited pools can be dominated by molecule-molecule collisions, in which case, the transition rates between them are k _exc = N _tot vσ _exc. Here, N _tot is the total density of molecules and σ _exc is the collision cross section measured from experiments.

Non-equilibrium among rotational and vibrational states in the model is solely introduced by the pump (i.e. the blue arrow in Fig. 1), and the resulting nonthermal distribution of population in some of the rotational levels can no longer be considered part of the thermal pool but must be modeled individually. These nonthermal states, illustrated in Fig. 1 as level J ₁, J ₂. . . etc. in each vibrational state, only exist in states that are of the same symmetry type and vibrational state as those connected by the pump. The total population of any state J_i (i.e. N _tot,J_i ) is simply the sum of its nonthermal and thermal parts, N _{tot,J
_i} = N_{J_i} + f_{J_i} N _pool, where f_{J_i} is the fraction of the total pool population N _pool in state J_i , and N_{J_i} is the amount the state population differs from rotational thermal equilibrium because of the pump. Note that N_{J_i} can be positive (excited vibrational state) or negative (ground vibrational state).

In general, the pump will create significant non-equilibrium not only in the states that it directly connects, but also in those adjacent to it. To account for this, our numerical model allows one to include as many rotational states as are necessary to model such a pump-induced non-equilibrium accurately. Generally, the number of J states to include depends on the ratio of the dipole-dipole rate to the K-swap rate (see end of paragraph for the description of these two processes) — the slower the K-swap rate, the more J levels one needs to consider. As already pointed out, the total population of any nonthermal state in this model is the combined population from both the nonthermal and thermal contributions, but since the thermal pools are in rotational equilibrium, the rotational relaxations, shown as black arrows in Fig. 1, are only modeled for the nonthermal population. In methyl fluoride, the rotational nonthermal processes are primarily dipole-dipole collisions whose rates are computed from the dipole-dipole collision cross section σ _DD measured from experiments [9

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in CH₃F: The ΔJ = n, ΔK = 0 processes,” J. Chem. Phys. 92, 6480–6491 (1990). [CrossRef]

]. Note that their allowed transitions are restricted to ΔJ = ±1 at the same K = K_i as the pumping transition. Lastly, the K-swap process [8

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]

] (i.e. the red arrows in Fig. 1) allows population to be transferred out of nonthermal states into the rest of the rotational levels within the same thermal pool, and its rate constant can again be measured experimentally. The K-swap collisional processes (which encompass the ΔK=3n process and the V-swap process [8

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]

, 15

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.

] when the A – and E – symmetry thermal pools are combined as they are here) convert the nonthermal molecules to thermal molecules but not vice versa, and thus, allowing the nonthermal states to return to zero when the pump is turned off.

All of the aforementioned energy transfer mechanisms except K-swap include reverse processes which account for the appropriate degeneracy of the upper and lower levels so that microscopic reversibility is satisfied. Moreover, spatial diffusion of the molecules within the cavity, though a relatively slow process, is also included in order to account for the molecular gas dynamics within the cavity. The pump transition rate, assuming a cylindrical waveguide cavity, is determined from R _pump = αP _pump/πR ² N _tot f _low hν _IR, where P _pump is the pump power, f _low = f _{J
₁} is the fraction of molecules in the lower ro-vibrational state J ₁ (as depicted in Fig. 1), ν _IR is the pump infrared frequency, and α, which takes into account both Doppler broadening and mismatch between pump and absorption lines, measures the gas’ ability to absorb the infrared radiation. Hence, pump transition rate may be increased by reducing R, increasing P _pump, or decreasing ν _IR. The last of these, however, is constrained by the need to match the frequency of the laser line with that of a ro-vibrational transition in the molecules, and in general, cannot be varied at will. Using such an implementation, our model accommodates both CW and pulsed pumping with arbitrary spatial profiles.

Lastly, we describe the various processes in the model via rate equations, within the rate-equation approximation [17

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

]. The most general rate description of a rotational level would include rotational relaxation terms due to dipole-dipole collisions among the non-thermal J levels, terms representing the pumping mechanism, K-swapping terms that allow flow of molecules out of the non-thermal levels (with equilibrating effects on non-thermal levels), and the stimulated emission terms at the FIR wavelength. Since we are interested in the steady state distribution of molecules across the cavity, spatial diffusion terms are also added to the rate equations describing each levels. Stimulated and spontaneous transitions at the lasing wavelength can, however, be neglected in the rate equations when operated near the threshold if the spontaneous emission lifetime, τ _spont, is relatively longer than all other transition lifetimes in the considered system. In methyl fluoride, for instance, τ _spont is on the order of 10s or more. In contrast, the next longest lifetime in the model occurs for the vibrational relaxation terms and yet, has lifetime on the order of 1 × 10⁻³ s or less. The relative magnitude of these radiative terms involving τ _spont in the rate equations are further diminished by the fact that only a fraction of these transitions contributes to the desired laser mode, especially when the cavity is large. With that, the rate equation describing the molecular density of J ₂ level in the excited vibrational state (i.e.

N_{J_{2}}^{exc}

) of Fig. 1 is given by

\begin{array}{l} \frac{\partial N_{J_{2}}^{exc} (r, t)}{\partial t} & = & - R_{pump} (r, t) Δ N_{pump} (r, t) - N_{J_{2}}^{exc} (r, t) \sum_{i = 1}^{3} \frac{1}{τ_{J_{2}, i}} + \sum_{j = 1}^{2} \frac{N_{j}^{exc} (r, t)}{τ_{j, J_{2}}} \\ + D \nabla^{2} N_{J_{2}}^{exc} (r, t) \end{array}

(1)

where

Δ N_{pump} (r, t) = N_{tot, J_{2}}^{exc} (r, t) - g_{J_{2}}^{exc} N_{tot, J_{1}}^{grd} (r, t) / g_{J_{1}}^{grd}

with

N_{tot, J_{1}}^{grd}

being the total molecular density of the level J ₁ in the ground vibrational state after accounting for the nonthermal density and the contribution from the corresponding fraction of the thermal pool, and

g_{J}^{grd} = 2 J + 1

is the degeneracy of the ground state J^th level. Similar definitions hold for

N_{tot, J_{2}}^{exc}

and

g_{J}^{exc}

. τ_{J
₂, i} two adjacent rotational levels and also the thermal for i = 1,2,3 is the decay lifetime out of

N_{J_{2}}^{exc}

to the two adjacent rotational levels and also the thermal pool. τ_{j,J
₂} for j = 1,2 represents the decay into

N_{J_{2}}^{exc}

from its two adjacent levels

N_{J_{1}}^{exc}

and

N_{J_{3}}^{exc}

. Finally, D is the diffusion constant given by the product of the molecule velocity and mean free path. In the same way, the thermal pool equation has terms from the K-swap and spatial diffusion processes, as well as vibrational relaxations due to molecule-molecule collisions and diffusion to wall. The equation describing the excited vibrational pool in Fig. 1 is then given by

\begin{array}{l} \frac{\partial N_{pool}^{exc} (r, t)}{\partial t} & = & - N_{pool}^{exc} (r, t) \sum_{k = 1}^{2} \frac{1}{τ_{exc, k} (r)} + \sum_{n = 1}^{2} \frac{N_{pool}^{n} (r, t)}{τ_{n, exc} (r)} + \sum_{m = 1}^{m_{0}} \frac{N_{J_{m}}^{exc} (r, t)}{τ_{J_{m}}^{swap}} \\ + D \nabla^{2} N_{pool}^{exc} (r, t) \end{array}

(2)

The first two terms relate to the vibrational relaxations among the thermal pools, where τ _exc, _k (r) for k = 1,2 is the decay lifetime from the excited thermal pool to the ground and higher excited pools, while τ_n _,exc(r) is the reverse. As noted earlier, the values of these diffusive lifetimes depend on their distances from the cell wall.

τ_{J_{m}}^{swap}

for m = 1,2,...,m ₀ in the third term is the lifetime of the transition from each of the m ₀ rotational levels to the excited pool due to the K-swap molecule-molecule collision process. As in Eq. (1), the last term represents spatial diffusion. Similar equations to Eq. (1) and Eq. (2) are used to represent every nonthermal level or thermal pool that contributes to the lasing action, which in turn depends on their degree of non-equilibrium caused by the pump (see Appendix). This way, the complete set of rate equations describing our laser system obeys the diffusion equation ∂N _tot(r,t)/∂t = D∇² N _tot(r,t), and so, conserves the total molecular density.

In our simulations, we consider only the immediate nonthermal levels adjacent to the pumped levels and initialize the density of molecules in the thermal pools according to a Boltzmann-distributed partition based on their respective vibrational energies. The non-thermal J_i levels are set equal to zero before any pumping occurs. The set of equations is then evolved in time until steady state is reached, allowing the temporal development of the lasing action to be tracked.

Laser action in this class of OPFIR molecular gas lasers is obtained as follows: Pumping is achieved with a line-tunable laser (e.g. CO₂ laser) which excites molecules from a specific rotational level in the ground vibrational state into a specific rotational level in an excited vibrational state. Because the excited vibrational level is comparatively empty, these photo-excited molecules create a population inversion between the pumped rotational state and the one immediately below it. This leads directly to stimulated emission between rotational states in the excited vibrational state. These photo-excited non-equilibrium molecules are subsequently rotationally and vibrationally relaxed to the ground state through collisions with other molecules and with the chamber walls, respectively.

The population inversion responsible for lasing exists as long as the pump-induced difference in nonthermal population between the upper and lower state N_{J_i} – N_{J_i} _–1 is greater than the difference in the corresponding pool-contributed thermal populations (f_{J_i} – f_{J_i} _–1)N _pool for those states. As the pressure grows, the pool molecules have an increasingly difficult time reaching the walls and de-exciting; thus, the thermal contribution from the pool increasingly quenches the non-thermal population inversion. Pool quenching creates the vibrational bottleneck that was once believed to limit OPFIR laser operation to low pressures regardless of pump intensity. However, it is now understood that as the excited vibrational state thermal pool fills, near resonant collisions among its constituents can excite molecules to even higher-lying vibrational levels, thereby providing an alternate relaxation pathway [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

]. Therefore, the pressure at which vibrational bottleneck occurs for a given laser geometry depends on the pump intensity. By concentrating on low pressure operation in this manuscript, these additional pathways may be ignored to simplify the analysis without affecting the generality of the results.

When molecules are photo-excited from the ground vibrational state to an excited state, two types of population inversion can be induced, and lasing arises from each. The first is the primary lasing transition already described between the pumped rotational state in the excited vibrational level and the rotational state below it. The second is the refilling transition between the IR pump depleted rotational state in the ground vibrational level and the rotational state above it. When the pump is eventually turned off, the K-swap process allows equilibrium conditions in the molecules to be restored. In the following illustration of this theoretical approach, without loss of generality, we present results only for the primary lasing transition.

3. Results and discussion

Figure 2(a) illustrates the setup of the OPFIR laser system considered in this section. A cylindrical waveguide resonator of length L _cell and radius R is filled with a suitable gas that lases at the desired THz frequency, while pump power at a much higher IR frequency enters the system from the front window. Depending on the reflectivities of the front and back windows at THz, lasing output power can then escape from the cavity via both channels. In order to compare our numerical calculations with past experimental results, the cavity modeled in the following simulation has length L _cell = 12 cm and front window reflectivity Γ_front = 0.96. The back window is opaque (i.e. Γ_back = 1). Furthermore, R is varied across a range of commercially available values from the smallest possible before modal cutoff at 0.08 cm to 1 cm. The cavity is then filled with an isotopic isomer of methyl fluoride gas, ¹³CH₃F, where a total of six rotational levels, J_ν _₀ = 3,4,5 in the ground vibrational state and J_ν _₃ = 4,5,6 in the excited state, are included (see detailed rate equations in the Appendix). Here, the K-swap processes occur at roughly 1–2 times the rates of the rotational transitions, and through this process, the lower and higher adjacent J levels can be modeled as part of their respective pool. For low pressure operation, we include the doubly degenerate ν ₆ state as our only higher excited vibrational state. As we shall see next, these relatively small number of rotational-vibrational levels are already sufficient to attain reasonable match between numerical predictions and experimental data. The system is uniformly pumped at room temperature (T = 300 K) with CW CO₂ laser at 31 THz in which molecules are excited between J_ν _₀ = 4 and J_ν _₃ = 5 to produce primary lasing transition at 0.245 THz in J_ν _₃ = 4 – 5 of the ν ₃ vibrational state, with τ _spont roughly chosen to be 15s [19

I. Shamah and G. Flynn, “Vibrational relaxation induced population inversions in laser pumped polyatomic molecules,” Chem. Phys. 55, 103–115 (1981). [CrossRef]

]. It is also noted that the occupancy of ν ₃ is comparatively low (f = e ^{−E_ν
/kT} ≈ 1/150), making the condition favorable for the creation of population inversion between J_ν _₃ = 5 and the level immediately below. Lastly, K_i in Fig. 1 equals 3 for 13CH₃F and we have approximated the K-swap process as the faster ΔK=3n process where the V-swap rate is slow and primarily equilibrates A and E states in the pool. Other molecules with a slower or absent ΔK=3n process will need to approximate the K-swap process by the V-swap process. The rate constants and other salient molecular parameters used are provided in the Appendix.

Fig. 2 (a) The OPFIR laser system considered in our numerical modeling: A cylindrical waveguide resonator of length L _cell and radius R is filled with a suitable gas that lases at the desired THz frequency, while pump power at a much higher IR frequency enters the system from the front window. In our example, ¹³CH₃F is pumped with CO₂ laser at 31 THz to produce lasing at 0.245 THz. (b) Frequency dependence of the ohmic losses for the five lowest order modes (i.e. TE₁₁, TM₀₁, TE₂₁, TE₀₁, TM₁₁) and three higher order ones with low losses (i.e. TE₀₂, TE₁₂, TE₂₂). Plot is for a copper cavity with R = 0.26 cm (and assuming L _cell ≫ R) so that several modes (not all shown in plot) have cut-offs below the THz lasing frequency. Despite this, only the lowest loss mode TE₀₁ exist in the cavity near the lasing threshold. Vertical green line shows the band of frequency at which our system operates. (c) Intensity profile of the three lowest loss modes supported by the R = 0.26 cm waveguide cavity.

Figure 2(b) shows a plot of the ohmic losses in a R = 0.26 cm copper waveguide cavity as a function of frequency. This loss prediction is made assuming L _cell ≫ R so that only the circular plane of the cavity matters and that the copper conductivity is high enough that the wall currents flow uniformly within a skin depth of the surface [18

R. Bansal (ed.), Handbook of Engineering Electromagnetics (Marcel Dekker, Inc., 2004). [CrossRef]

]. At the operating frequency, shown as a green vertical band in figure, the cavity supports multiple modes but only the five lowest order modes and three higher order ones that have low ohmic losses are included in the plot. The cut-off frequency of each of these modes can also be deduced from the figure. The lowest loss mode, TE₀₁, will be the first to lase upon threshold and remains the only mode present for pump power near the threshold, despite the cavity’s ability to support multiple modes at THz. Alternatively, one can ensure single-mode operation even for large pump power by selecting the smallest possible cavity size with radius ∼ λ _THz/2, where λ _THz is the wavelength of the THz lasing output. For our purpose, it suffices to simply operate near the threshold for single-mode operation since this also facilitates independent design optimizations of our system based on molecular gas physics and photonics considerations. Figure 2(c) depicts the intensity profiles of the three lowest loss modes of the cavity. It may be noted that the magnitude of ohmic loss scales roughly with the fraction of modal intensity residing near the cell wall, which in turn relates to the amount of penetration, to within a skin depth, into the wall. At low enough pump rate, only losses of TE₀₁ mode (both ohmic losses and leakage through the end mirrors) can be compensated by the gain produced in ¹³CH₃F while the other modes remain suppressed. Thus, we shall consider TE₀₁ to be the lowest-threshold lasing mode.

Figure 3 plots typical spatio-temporal results of the numerical simulation. In particular, the evolution of the laser transition ΔN(λ _THz) and the ν ₃ thermal pool to steady state are presented and compared at 250 and 350 mTorr, assuming P _pump to be 10 watts in CW operation. A few initial observations are in order: (i) the radial spatial profile of ΔN(r) is the inverse of that for ν ₃ pool. This is the pool quenching that leads to the vibrational bottleneck mentioned earlier. Note the rapid removal of molecules from ν ₃ near the cell wall is critical to maintaining the inversion. (ii) The cell center takes a long time to attain steady state given that the wall rate is slow far away from the wall (i.e. molecules have to travel a long distance to reach the cell boundary). (iii) Similarly, a longer time is required to reach steady state at a higher pressure due to the decreased k _wall near the cell center. In Fig. 3, steady state is attained after 100 μs at 250 mTorr while it takes 300 μs to do so at 350 mTorr. (iv) The magnitude of the variation across the radial direction of the cavity is greater at higher pressure, as will be explained shortly. While we mostly deal with CW systems in this paper, the ability to track the full temporal development of the laser action allows pulsed systems to be studied as well.

Fig. 3 Time dependent properties of laser system with CW pumping at P _pump = 10 W. (a) Time evolution of the inversion in ν ₃ state at 250 mTorr across the radial axis of the cavity until steady state behavior is observed. (b) Time evolution of the ν ₃ thermal pool population at 250 mTorr across the radial axis of the cavity until steady state behavior is observed. (c) Same as in (a) except operated at 350 mTorr (using the same magnitude range for the color bar). (d) Same as in (b) except operated at 350 mTorr (using the same magnitude range for the color bar).

We next verify our numerical model against experimental results for a R = 0.26 cm copper cavity ¹³CH₃F laser. The experimental configuration was identical to that described in [6

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]

], with a 1 mm diameter output coupler hole in the front reflector through which the IR pump and THz laser power were transmitted [Fig. 2(a)]. The slope efficiency was measured as the ratio of the THz laser power emitted to the corresponding IR pump intensity, and the threshold pump intensity was the extrapolation of this slope to zero THz power. Numerical predictions and experimental measurements of both the pump power and slope efficiency at threshold are compared at two pressures: 100 mTorr and 300 mTorr. Here, we adopt a semi-analytic approach for these comparisons, where knowledge of the losses in the OPFIR system and population inversions calculated with our model are used simultaneously. In the following, we first discuss the match for the threshold power intensity before addressing that of the slope efficiency.

The blue lines in Figs. 4(a) and 4(b) depict the unsaturated effective gain values of the system, as predicted by our model, for a range of pump intensity at 100 and 300 mTorr respectively. The gain for the J = 4 – 5 lasing transition in ν ₃ state of 13CH₃F may be written as [20

A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, 2007).

, 21

S. L. 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–1562 (2011). [CrossRef] [PubMed]

]

〈 γ (λ_{THz}) 〉 = 〈 Δ N 〉 \frac{λ_{THz}^{2}}{4 π^{2} n^{2} τ_{spont} Δ ν (λ_{THz})}

(3)

where 〈ΔN〉 is the pressure-dependent steady-state effective inversion defined as

〈 Δ N (t_{\infty}) 〉 = \frac{\int_{V_{ACT}} d r | E_{0} (r) |^{2} Δ N (r, t_{\infty})}{\int_{V_{ACT}} d r | E_{0} (r) |^{2}} .

(4)

E ₀(r) is the normalized mode profile of the passive cavity (∫ d r ε ₀ n ²(r) |E ₀(r)|² = 1) and V _ACT being the volume of active region of the considered structure. In our case, V _ACT is the entire cavity (note that this implies that the gas’s behavior is less relevant in regions where the modal intensity is small but plays an important role in regions of high field intensity). In the considered case,

Δ N (r) = N_{tot, J = 5}^{ν_{3}} (r) - g_{J = 4}^{ν_{3}} N_{tot, J = 4}^{ν_{3}} (r) / g_{J = 5}^{ν_{3}}

with

N_{tot, J}^{ν_{3}}

being the total molecular density of the upper (J = 5) and lower (J = 4) lasing transition level in ν ₃ state after accounting for the fractional contribution from the excited thermal pool, and

g_{J}^{ν_{3}} = 2 J + 1

is the degeneracy of the excited state’s J^th level. E ₀(r) is the normalized field of the TE₀₁ mode and n ≈ 1 for the gas. The width of the transition line, Δν, is broadened due to a combination of Doppler effects and molecular collisions and has a strong dependence on both pressure and transition frequency. The former is associated with the Gaussian lineshape while the latter, which leads to pressure broadening, is Lorentzian. Hence,

Δ ν = {(Δ ν_{D}^{2} + Δ ν_{P}^{2})}^{1 / 2}

where the two terms are the Doppler and pressure broadening widths respectively, and the overall lineshape function used is the one that corresponds to the larger broadening parameter of the two. It is noted that for the same half width, pressure-broadened lines are stronger far away from the line-center than the Doppler-broadened ones. The gain defined in Eq. (3) accounts for the number of modes only in the frequency interval ν _THz ± Δν/2 by assuming a large optical resonator whose dimension is several wavelengths so that its exact shape becomes less significant. We find that this applies in our case since the cylindrical copper waveguide is almost 100λ _THz long and a few λ _THzs in diameter. Moreover, the inversion density, and hence gain, is not saturated by stimulated emission for the entire range of input intensity plotted. As we are only interested in the threshold value in this part of the work, such approach remains valid.

Fig. 4 (a) Method of threshold predictions via gain-loss balancing used in plotting (c) for every pressure point. Here, we operate at 100 mTorr with R = 0.26 cm. The blue line is the unsaturated effective gain value for a range of pump intensity predicted by the aforementioned model while the red line is the magnitude of the loss coefficient due to absorption by the metallic cavity and energy leakage from the front window. Intersection between the two lines corresponds to the lasing threshold. (b) Same as in (a) except that the pressure is at 300 mTorr. (c) Left axis depicts the relationship between the threshold intensity and pressure for a range of R values from 0.08 cm to 1 cm as predicted from numerical model. The corresponding linewidth, Δν, of the lasing transition (over which it may be tuned) is also illustrated in the top axis. A general trend exists such that the lasing threshold increases with pressure. The right axis shows the mean free path, λ _MFP, within the gaseous system which at a fixed temperature, is inversely proportional to pressure. Inset is the ohmic loss [see Fig. 2(b)] of the TE₀₁ mode at 0.245 THz, with the radii of interest marked as square markers.

The threshold for lasing may be approximated (neglecting spontaneous emission) to occur when the gain balances the total losses of the OPFIR system. Losses in our setup are assumed to be due to either ohmic losses from the E-field penetration into the copper wall or to energy leakage from the front mirror. For L _cell ≫ R, we can treat the system as one-dimensional along the cylinder axis for loss analysis and quantify its coefficient as 2α _ohmic – ln(Γ_frontΓ_back)/L _cell, where α _ohmic is the amplitude loss coefficient shown in Fig. 2(b). The magnitudes of these loss coefficients are indicated as horizontal red lines in Figs. 4(a) and 4(b). As such, the input threshold intensity predicted from our model is determined at the intersection between the red and blue lines, which indicates the point at which losses are compensated by the gain at λ _THz so that any further increase in pump intensity saturates the gain and is thereby channeled into coherent output power that escapes from the front window. From Figs. 4(a) and 4(b), the numerical thresholds attained are 0.78 and 2.42 W/cm² for the laser operated at 100 mTorr and 300 mTorr, respectively.

To make realistic comparisons with experiment, the decreasing absorption of the IR pump beam as it propagates along the cavity has to be considered. The pressure-dependent IR absorption coefficient, α′, for ¹³CH₃F is 71 m⁻¹Torr⁻¹ [14

H. O. Everitt, “Collisional Energy Transfer in Methyl Halides,” PhD Thesis (Department of Physics, Duke University, 1990).

] so that the ability of the pump to excite the gas decreases with increasing pressure. The effect is to reduce the amount of IR power effectively pumping the gas and can be captured by defining an incident pump power

P_{pump}^{inc} = P_{pump} / κ

, where κ is a pressure-dependent factor that accounts for the net fraction of incident IR power available for excitation in the L _cell long cavity. This factor is simply given by

(1 / L_{cell}) \times \int_{0}^{L_{cell}} e^{- α^{'} x} d x

. Hence, κ = 0.67 and 0.36 for our 100 mTorr and 300 mTorr setups resulting in actual pump thresholds (

P_{pump}^{inc} / π R^{2}

) of 1.2 and 6.7 W/cm², respectively. These results compare well with experimentally measured thresholds of 1.1±0.4 and 5.5±1.5 W/cm². Possible sources of discrepancy in the comparison include instabilities that make measuring the threshold difficult, deviation from spatially uniform pumping in the experiment, the one-dimensional treatment in balancing gain-loss, and the possibility that the emission was from the refilling transition in the ν ₀ state. Even so, the current treatment provides adequately accurate predictions of the lasing thresholds. More importantly, this comparison acts to verify our claim that the aforementioned model of the OPFIR is physically robust and rigorous, yet numerically tractable.

Next, we present a similar comparison for the slope efficiency at threshold using a semi-analytic approach. The set of rate equations in Eq. (1), with the inclusion of the induced lasing transition rate, are solved analytically at steady-state to obtain the following approximate relationship between output power P _out and input power P _in above threshold:

P_{out} = f_{out} P_{sat} (\frac{P_{in}}{P_{in}^{th}} - 1)

(5)

where f _out accounts for the fraction of emitted power that escapes via the front window, and

P_{in}^{th}

is the IR power at threshold. P _sat is the saturation power at THz in which the gain in the system is half its unsaturated value. In this particular system, P _sat = h̄ω _THz ∫ _V _{_ACT} d rΔN _th(r)/τ _eff where ΔN _th is the inversion at threshold and τ _eff is an effective lifetime that accounts for the various rates between the rotational levels in the lasing ν ₃ state, and is found to be

τ_{eff} = \frac{1 / τ_{J_{1}}^{swap} - 1 / τ^{'}}{1 / τ_{J_{1}}^{swap} τ_{J_{2}, J_{1}} - 1 / τ^{'} τ_{J_{1}}}

(6)

with

\begin{matrix} \frac{1}{τ^{'}} = \frac{τ_{J_{3}} - τ_{J_{3}, J_{2}} - τ_{J_{2}, J_{3}} τ_{J_{3}, J_{2}} / τ_{J_{2}}^{swap}}{τ_{J_{2}, J_{3}} τ_{J_{3}, J_{2}}}, \\ \frac{1}{τ_{J_{1}}} = \frac{1}{τ_{J_{1}, J_{2}}} + \frac{1}{τ_{J_{1}}^{swap}} and \frac{1}{τ_{J_{3}}} = \frac{1}{τ_{J_{3}, J_{2}}} + \frac{1}{τ_{J_{3}}^{swap}} . \end{matrix}

The lifetimes in Eq. (6) describe the transition rates in the excited ν ₃ state where lasing occurs and are defined in accordance to their definitions in Fig. 1 and Eq. (1). With that, the slope dP _out/dP _in is

f_{out} P_{sat} / P_{in}^{th}

. Using numerical results of ΔN _th and

P_{in}^{th}

found from the full model and accounting for the pressure-dependent IR absorption discussed in the previous paragraph, we obtain a threshold slope efficiency of 24.5×10⁻⁶ and 39.2×10⁻⁶ at 100 mTorr and 300 mTorr, respectively. Note that the best slope efficiency for the considered ¹³CH₃F laser is 8 × 10⁻³, which is the Manley-Rowe limit when the quantum efficiency is 1. These calculated values compare reasonably well with experimental slope values of (44 ± 22) × 10⁻⁶ and (48 ± 24) × 10⁻⁶, for which the greatest source of uncertainty is the calibration of the the Golay cell used to detect the THz radiation. Given this, these results further validate our model.

We extend the threshold determination method of Figs. 4(a) and 4(b) to study how the threshold intensity behaves as both cell radius and pressure are varied in Fig. 4(c). A general trend exists such that the lasing threshold increases with pressure for every radius plotted. This can primarily be attributed to the fact that the linewidth, Δν (λ _THz), broadens with increasing pressure so that 〈γ(λ _THz)〉 reduces [see Eq. (3)]. Note that the variation of the linewidth with pressure is captured in the top axis of the figure. Since the total loss of the system is pressure-independent, such gain reduction increases the threshold linearly with pressure for each value of R. Second order effects that further raise the threshold at high pressure originate from the vibrational bottleneck: A higher pressure impedes the molecular diffusion to the cell wall and the subsequent decrease in wall collision rates reduces the flow of molecules out of the the excited thermal pool, leading to enhanced absorption and hence, a weaker inversion 〈ΔN〉 and gain. Equivalently, the same effects on 〈ΔN〉 can be understood by recognizing that λ _MFP progressively makes up a smaller region near the cell wall as the pressure increases such that a greater proportion of the molecules are in the non-ballistic regime, which also leads to slower k _wall. Red dashed line plotted against the right axis depicts the pressure dependence of λ _MFP. Over the range of low pressures considered, the vibrational bottleneck is only observed for large radius cells (i.e. R ≥ 0.26 cm in the current setup) while for smaller cavities, there is no vibrational bottleneck and the threshold dependence on pressure remains approximately linear.

From Fig. 4(c), it can also be seen that the threshold intensity increases with radius for R values larger than 0.26 cm but decreases with radius for R values smaller than 0.26 cm. To understand why R = 0.26 cm is nearly optimal, the reader is referred to the inset that illustrates the reduction in ohmic losses of the TE₀₁ mode at 0.245 THz when R grows from the waveguide cut-off value, with the radii of interest examined in the main figure marked as square markers. For small radius, the ohmic loss is high and hence, possesses great influence on the direction of threshold shift. In this regime, when the radius increases, the corresponding fall in ohmic loss has the tendency to lower the threshold. On the other hand, at larger R values where the ohmic losses are low, the vibrational bottleneck increasingly raises or prevents threshold with increasing R and pressure. However, at pressures less than 100 mTorr in large cells (i.e. R ≥ 0.26 cm), the pressure-dependent threshold lines in Fig. 4(c) converge to the linear behavior observed for cells with R ≤ 0.26 cm. As we shall see next, this follows directly from the fact that λ _MFP grows large at low pressures, so for a constant intensity and pump rate R _pump, the inversion 〈ΔN〉 and gain exhibit little radial dependence. As a result, variation of R has a much smaller effect on the wall collision rate, and correspondingly on 〈γ〉, at lower pressures. So, for cell radii where the ohmic losses are also low, the threshold remains almost unchanged as R is varied.

To understand the spatial effects on laser performance more clearly, the radial dependence of the inversion and gain in a R = 0.5 cm cavity pumped by CW CO₂ laser at P _pump = 100 W are studied in Fig. 5. From the intensity-weighted averaging scheme used in Eq. (4), we have already noted that it is most desirable to match the peak of the cavity’s eigenmode intensity to regions where ΔN(r) or γ(r) is maximum. Such consideration plays an important role in designing the cavity itself. Figure 5(a) illustrates the distribution of ΔN(r) within the cavity for values of pressure ranging from 50 mTorr to 450 mTorr. In general, strong inversions are favored near the cell wall where the wall collision rates are the highest, and so, allow rapid depopulation of the excited and higher excited thermal pools to reduce pool quenching. This directly implies an advantage in exciting a lasing mode whose intensity peaks near the cell wall. However, such approach will also have to be weighted against the higher ohmic losses that arise from the field’s enhanced interaction with the cell, in a similar trade-off effects to what was already discussed in Fig. 4(c). The details and treatment of the cavity design will be left for future work. We further note that the magnitude of radial variation in the cell increases with pressure. For instance, at 50 mTorr, a good proportion of the cavity is ballistic [see λ _MFP variation in Fig. 4(c)] so that ΔN is approximately uniform, i.e. ΔN(r _cen) = 0.96ΔN(r _wall) where r _cen and r _wall indicate the position at the cell center and wall, respectively. On the other hand, at 450 mTorr, most of the molecules are in the non-ballistic regime where k _wall ∝1/r ² as opposed to k _direct ∝1/r. This leads to ΔN(r _cen) to be only 14% of ΔN(r _wall). Indeed, it was found that the maximum inversion for the set of parameters used in this example increases with pressure because of the corresponding increase in N _tot. Figure 5(b) depicts a similar plot for the gain. The same characteristics are observed among the two except that γ(r) is also influenced by the laser transition linewidth which broadens as pressure increases. This decrease of gain with pressure is again consistent with the discussion presented for Fig. 4(c), explaining why most OPFIR lasers use large diameter cavities operating at low pressure.

Fig. 5 (a) Radial spatial variation of the excited inversion in ν ₃ state of ¹³CH₃F for a range of operating pressure from 50 to 450 mTorr. 0 cm marks the center of the cylindrical cavity while 0.5 cm marks the copper wall. All four panels in this figure assumed CW pumping at P _pump = 100 W. Near the cell wall, the diffusion is ballistic allowing significant reduction of molecules from the ν ₃ pool, and hence, results in a higher population inversion. (b) Same as in (a) except that the transition width at 0.245 THz is also factored in to study the gain. (c) Gain dependence on pressure for a fixed set of pump parameters. From a molecular gas physics standpoint, results clearly indicate that small R cavities are favored in terms of gain magnitude and pressure cut-off. (d) Data extracted from (c) where the left axis and right axis respectively depict the optimum operating pressure and the corresponding gain, as a function of cell radius R. Again, results here favored small sized cavity for high pressure operation.

However, this design misses a critical point. Since R _pump increases when the same pump power is more tightly confined, the gain tends to grow as cavity radii shrink until it eventually saturates (because N _tot is a conserved quantity). Figures 5(c) and 5(d) track the optimum radius size for lasing based on the numerical model developed. An optimum operating pressure exists for each radius where the gain is maximized. Beyond that, vibrational bottleneck begins to play a dominant role, explaining the fall of the gain as pressure rises. At high enough pressure, the lasing action will be quenched. For example, in Fig. 5(c) the cut-off pressure occurs at 380 mTorr for R = 1 cm. Because of the role of higher lying vibrational levels, the cutoff pressure is pump intensity dependent such that the greater the pumping intensity, the higher the cut-off pressure. Thus, the optimum pressure and gain results summarized in Fig. 5(d) based on the molecular gas physics greatly favor small-sized cavities, especially if the effects of waveguide ohmic losses can be overcome.

4. Conclusions

We have presented a comprehensive theoretical analysis of lasing action in optically-pumped rotationally excited molecular gases confined in a metallic cavity. In order to do that, we have developed a theoretical framework that provides a realistic spatio-temporal description of the molecular de-excitation and diffusion processes occurring inside the cavity. Our approach also accounts for the influence on the main laser features of the electric-field intensity distribution and the ohmic losses associated to the different cavity modes. To illustrate this theoretical framework, we have analyzed in detail the case of lasing action from optically-pumped methyl fluoride gas confined in a cylindrical copper cavity. Our simulations have shown that for this system, the region of maximum population inversion is located near the cavity walls. The outcome of our calculations has also shown that the lowest lasing threshold mode corresponds to the cavity mode whose electric-field intensity profile presents an optimal spatial overlap with the corresponding population inversion distribution inside the cavity, and at the same time, minimizes the ohmic losses at the lasing frequency. Furthermore, we have predicted the existence of an optimal radius, determined by the interplay between the overall gain in the system and the corresponding ohmic losses, which minimizes the lasing threshold intensity for a large interval of gas pressures. We believe the results reported in this paper could contribute to the development of a novel class of room-temperature compact sources of far-infrared and THz radiation.

Appendices

Appendix: Rate equations and rate constants of ¹³CH₃F gas lasers

In this appendix, we first present the set of rate equations used to model the ¹³CH₃F gas lasers. Figure 6 and Eq. (7) to Eq. (15) refer to a specific example of a methyl fluoride gas discussed in Sec. 3, adopted from the general model presented in Fig. 1, Eq. (1) and Eq. (2) of Sec. 2. We have labeled the rotational levels from N ₁ to N ₆, and the thermal pools as N_A , N_B and N_C . Hence, τ ₁₂, τ ₂₃, τ ₄₅, τ ₅₆, and their reverse processes represent rotational transitions due to dipole-dipole collisions, τ ₁ _A , τ ₄ _B , and similar processes represent K-swap transition lifetimes from the rotational levels (i.e. N ₁ and N ₄) into their respective thermal pools (i.e. N_A and N_B ), and τ_AB , τ_AC , and their reverse processes are vibrational transitions which arise due to collisions with the cell wall. Readers should note that all notations in the main body of the paper are kept consistent with those used in Fig. 1, Eq. (1) and Eq. (2).

Fig. 6 ¹³CH₃F OPFIR laser model used in this paper: results of diagnostic studies and theoretical considerations restrict the degrees of freedom so that a numerically tractable, yet physically accurate model is attained. Details on each process labeled are described in the text.

\frac{\partial N_{1} (r, t)}{\partial t} = - N_{1} (r, t) (\frac{1}{τ_{12}} + \frac{1}{τ_{1 A}}) + \frac{N_{2} (r, t)}{τ_{21}} + D \nabla^{2} N_{1} (r, t)

(7)

\begin{array}{l} \frac{\partial N_{2} (r, t)}{\partial t} & = & - N_{2} (r, t) (\frac{1}{τ_{21}} + \frac{1}{τ_{23}} + \frac{1}{τ_{2 A}}) + \frac{N_{1} (r, t)}{τ_{12}} + \frac{N_{3} (r, t)}{τ_{32}} \\ + R_{pump} (r, t) Δ N_{pump} (r, t) + D \nabla^{2} N_{2} (r, t) \end{array}

(8)

\frac{\partial N_{3} (r, t)}{\partial t} = - N_{3} (r, t) (\frac{1}{τ_{32}} + \frac{1}{τ_{3 A}}) + \frac{N_{2} (r, t)}{τ_{23}} + D \nabla^{2} N_{3} (r, t)

(9)

\frac{\partial N_{4} (r, t)}{\partial t} = - N_{4} (r, t) (\frac{1}{τ_{45}} + \frac{1}{τ_{4 B}}) + \frac{N_{5} (r, t)}{τ_{54}} + D \nabla^{2} N_{4} (r, t)

(10)

\begin{array}{l} \frac{\partial N_{5} (r, t)}{\partial t} & = & - N_{5} (r, t) (\frac{1}{τ_{54}} + \frac{1}{τ_{56}} + \frac{1}{τ_{5 B}}) + \frac{N_{4} (r, t)}{τ_{45}} + \frac{N_{6} (r, t)}{τ_{65}} \\ - R_{pump} (r, t) Δ N_{pump} (r, t) + D \nabla^{2} N_{5} (r, t) \end{array}

(11)

\frac{\partial N_{6} (r, t)}{\partial t} = - N_{6} (r, t) (\frac{1}{τ_{65}} + \frac{1}{τ_{6 B}}) + \frac{N_{5} (r, t)}{τ_{56}} + D \nabla^{2} N_{6} (r, t)

(12)

\begin{array}{l} \frac{\partial N_{A} (r, t)}{\partial t} & = & - N_{A} (r, t) (\frac{1}{τ_{A B} (r)} + \frac{1}{τ_{A C} (r)}) + \frac{N_{B} (r, t)}{τ_{B A} (r)} + \frac{N_{C} (r, t)}{τ_{C A} (r)} \\ + \frac{N_{1} (r, t)}{τ_{1 A}} + \frac{N_{2} (r, t)}{τ_{2 A}} + \frac{N_{3} (r, t)}{τ_{3 A}} + D \nabla^{2} N_{A} (r, t) \end{array}

(13)

\begin{array}{l} \frac{\partial N_{B} (r, t)}{\partial t} & = & - N_{B} (r, t) (\frac{1}{τ_{B A} (r)} + \frac{1}{τ_{B C}}) + \frac{N_{A} (r, t)}{τ_{A B} (r)} + \frac{N_{C} (r, t)}{τ_{C B}} \\ + \frac{N_{4} (r, t)}{τ_{4 B}} + \frac{N_{5} (r, t)}{τ_{5 B}} + \frac{N_{6} (r, t)}{τ_{6 B}} + D \nabla^{2} N_{B} (r, t) \end{array}

(14)

\begin{array}{l} \frac{\partial N_{C} (r, t)}{\partial t} & = & - N_{C} (r, t) (\frac{1}{τ_{C A} (r)} + \frac{1}{τ_{C B}}) + \frac{N_{A} (r, t)}{τ_{A C} (r)} + \frac{N_{B} (r, t)}{τ_{B C}} \\ + D \nabla^{2} N_{C} (r, t) \end{array}

(15)

where

Δ N_{pump} (r, t) = N_{5}^{tot} (r, t) - g_{5} N_{2}^{tot} (r, t) / g_{2}

with

N_{2}^{tot}

(r,t) being the total molecular density of J = 4 nonthermal level in ν ₀ state after accounting for the fractional contribution from thermal pool A, and g ₂ = 9 is its degeneracy. Similar definitions hold for

N_{5}^{tot}

and g ₅ = 11.

Table 1 provides the rate constants and relevant molecular parameters used in our calculations [8

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]

, 14

H. O. Everitt, “Collisional Energy Transfer in Methyl Halides,” PhD Thesis (Department of Physics, Duke University, 1990).

, 15

H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.

Table 1 Translational Cross Sections (in Å ²) of ¹³CH₃F

Mass (AMU)	35
Gas kinetic collision cross section	44
ν ₃ → ν ₆ cross section at 300 K, σ _exc	2.65
Dipole-dipole cross section at 300 K, σ _DD	320
K-swap cross section at 300 K, σ _K-swap	137

Acknowledgments

We acknowledge helpful discussions with Dr. Ofer Shapira and Prof. Steven Johnson. This work was supported in part by the MRSEC Program of the National Science Foundation under award number DMR-0819762, the Army Research Office through the Institute for Soldier Nanotechnologies under Contract No. W911NF-07-D0004, and the Army’s competitive in-house innovative laboratory research program. J.B.-A. acknowledges financial support by Spanish MICINN through Ramón y Cajal program, grant no. RyC-2009-05489. Song-Liang Chua acknowledges the financial support from the DSO National Laboratories, Singapore.

References and links

1.	T. Y. Chang and T. J. Bridges, “Laser actions at 452, 496, and 541 μm in optically pumped CH₃F,” Opt. Commun. 1, 423–426 (1970). [CrossRef]
2.	T. K. Plant, L. A. Newman, E. J. Danielewicz, T. A. DeTemple, and P. D. Coleman, “High power optically pumped far infrared lasers,” IEEE Trans. Microwave Theory Tech. MT22, 988–990 (1980).
3.	W. H. Matteson and F. C. De Lucia, “Millimeter wave spectroscopic studies of collision-induced energy transfer processes in the ¹³CH₃F laser,” IEEE J. Quantum Electron. 19, 1284–1293 (1983). [CrossRef]
4.	M. S. Tobin, “A review of optically pumped NMMW lasers,” Proc. IEEE 73, 61–85 (1985). [CrossRef]
5.	P. K. Cheo (ed.), Handbook of Molecular Lasers (Marcel Dekker, Inc., 1987), pp. 495–569.
6.	H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]
7.	R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).
8.	H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of ¹³CH₃F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]
9.	H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in CH₃F: The ΔJ = n, ΔK = 0 processes,” J. Chem. Phys. 92, 6480–6491 (1990). [CrossRef]
10.	R. L. Crownover, H. O. Everitt, D. D. Skatrud, and F. C. DeLucia, “Frequency stability and reproductibility of optically pumped far-infrared lasers,” Appl. Phys. Lett. 57, 2882–2884 (1990). [CrossRef]
11.	D. Dangoisse, P. Glorieux, and J. Wascat, “Diffusion and vibrational bottleneck in optically pumped submillimeter laser,” Int. J. Infrared Milimeter Waves 2, 215–229 (1981). [CrossRef]
12.	J. O. Henningsen and H. G. Jensen, “The optically pumped far-infrared laser: Rate equations and diagnostic experiments,” IEEE J. Quantum Electron. QE-11, 248–252 (1975). [CrossRef]
13.	R. J. Temkins and D. R. Cohn, “Rate equations for an optically pumped, far-infrared laser,” Opt. Commun. 16, 213–217 (1976). [CrossRef]
14.	H. O. Everitt, “Collisional Energy Transfer in Methyl Halides,” PhD Thesis (Department of Physics, Duke University, 1990).
15.	H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.
16.	L. E. Reichl, A Modern Course in Statistical Physics (John Wiley & Sons Inc., 1998).
17.	A. E. Siegman, Lasers (Univ. Science Books, 1986).
18.	R. Bansal (ed.), Handbook of Engineering Electromagnetics (Marcel Dekker, Inc., 2004). [CrossRef]
19.	I. Shamah and G. Flynn, “Vibrational relaxation induced population inversions in laser pumped polyatomic molecules,” Chem. Phys. 55, 103–115 (1981). [CrossRef]
20.	A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, 2007).
21.	S. L. 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–1562 (2011). [CrossRef] [PubMed]

OCIS Codes
(140.3070) Lasers and laser optics : Infrared and far-infrared lasers
(140.3460) Lasers and laser optics : Lasers
(140.4130) Lasers and laser optics : Molecular gas lasers

ToC Category:
Lasers and Laser Optics

History
Original Manuscript: February 10, 2011
Manuscript Accepted: March 21, 2011
Published: April 4, 2011

Citation
Song-Liang Chua, Christine A. Caccamise, Dane J. Phillips, John D. Joannopoulos, Marin Soljačić, Henry O. Everitt, and Jorge Bravo-Abad, "Spatio-temporal theory of lasing action in optically-pumped rotationally excited molecular gases," Opt. Express 19, 7513-7529 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-8-7513

Sort: Author | Year | Journal | Reset

References

T. Y. Chang and T. J. Bridges, “Laser actions at 452, 496, and 541 μm in optically pumped CH3F,” Opt. Commun. 1, 423–426 (1970). [CrossRef]
T. K. Plant, L. A. Newman, E. J. Danielewicz, T. A. DeTemple, and P. D. Coleman, “High power optically pumped far infrared lasers,” IEEE Trans. Microwave Theory Tech. MT22, 988–990 (1980).
W. H. Matteson and F. C. De Lucia, “Millimeter wave spectroscopic studies of collision-induced energy transfer processes in the 13CH3F laser,” IEEE J. Quantum Electron. 19, 1284–1293 (1983). [CrossRef]
M. S. Tobin, “A review of optically pumped NMMW lasers,” Proc. IEEE 73, 61–85 (1985). [CrossRef]
P. K. Cheo (ed.), Handbook of Molecular Lasers (Marcel Dekker, Inc., 1987), pp. 495–569.
H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]
R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).
H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of 13CH3F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]
H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in CH3F: The ΔJ = n, ΔK = 0 processes,” J. Chem. Phys. 92, 6480–6491 (1990). [CrossRef]
R. L. Crownover, H. O. Everitt, D. D. Skatrud, and F. C. DeLucia, “Frequency stability and reproductibility of optically pumped far-infrared lasers,” Appl. Phys. Lett. 57, 2882–2884 (1990). [CrossRef]
D. Dangoisse, P. Glorieux, and J. Wascat, “Diffusion and vibrational bottleneck in optically pumped submillimeter laser,” Int. J. Infrared Milimeter Waves 2, 215–229 (1981). [CrossRef]
J. O. Henningsen and H. G. Jensen, “The optically pumped far-infrared laser: Rate equations and diagnostic experiments,” IEEE J. Quantum Electron. QE-11, 248–252 (1975). [CrossRef]
R. J. Temkins and D. R. Cohn, “Rate equations for an optically pumped, far-infrared laser,” Opt. Commun. 16, 213–217 (1976). [CrossRef]
H. O. Everitt, “Collisional Energy Transfer in Methyl Halides,” PhD Thesis (Department of Physics, Duke University, 1990).
H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.
L. E. Reichl, A Modern Course in Statistical Physics (John Wiley & Sons Inc., 1998).
A. E. Siegman, Lasers (Univ. Science Books, 1986).
R. Bansal (ed.), Handbook of Engineering Electromagnetics (Marcel Dekker, Inc., 2004). [CrossRef]
I. Shamah and G. Flynn, “Vibrational relaxation induced population inversions in laser pumped polyatomic molecules,” Chem. Phys. 55, 103–115 (1981). [CrossRef]
A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, 2007).
S. L. 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–1562 (2011). [CrossRef] [PubMed]

Bravo-Abad, J.
Bridges, T. J.
Chang, T. Y.
Chong, Y. D.
Chua, S. L.
Cohn, D. R.
Coleman, P. D.
Crownover, R. L.
Dangoisse, D.
Danielewicz, E. J.
De Lucia, F. C.
DeLucia, F. C.
DeTemple, T. A.
Everitt, H. O.
Flynn, G.
Glorieux, P.
Henningsen, J. O.
Jensen, H. G.
Matteson, W. H.
McCormick, R. I.
Newman, L. A.
Plant, T. K.
Reichl, L. E.
Shamah, I.
Siegman, A. E.
Skatrud, D. D.
Soljac?ic, M.
Stone, A. D.
Temkins, R. J.
Tobin, M. S.
Wascat, J.
Yariv, A.
Yeh, P.

Appl. Phys. Lett.

H. O. Everitt, D. D. Skatrud, and F. C. De Lucia, “Dynamics and tunability of a small optically pumped CW far-infrared laser,” Appl. Phys. Lett. 49, 995–997 (1986). [CrossRef]
R. L. Crownover, H. O. Everitt, D. D. Skatrud, and F. C. DeLucia, “Frequency stability and reproductibility of optically pumped far-infrared lasers,” Appl. Phys. Lett. 57, 2882–2884 (1990). [CrossRef]

Chem. Phys.

I. Shamah and G. Flynn, “Vibrational relaxation induced population inversions in laser pumped polyatomic molecules,” Chem. Phys. 55, 103–115 (1981). [CrossRef]

IEEE J. Quantum Electron.

J. O. Henningsen and H. G. Jensen, “The optically pumped far-infrared laser: Rate equations and diagnostic experiments,” IEEE J. Quantum Electron. QE-11, 248–252 (1975). [CrossRef]
R. I. McCormick, H. O. Everitt, F. C. De Lucia, and D. D. Skatrud, “Collisional energy transfer in optically pumped far-infrared lasers,” IEEE J. Quantum Electron. QE-23, 2069–2077 (1989).
W. H. Matteson and F. C. De Lucia, “Millimeter wave spectroscopic studies of collision-induced energy transfer processes in the 13CH3F laser,” IEEE J. Quantum Electron. 19, 1284–1293 (1983). [CrossRef]

IEEE Trans. Microwave Theory Tech.

T. K. Plant, L. A. Newman, E. J. Danielewicz, T. A. DeTemple, and P. D. Coleman, “High power optically pumped far infrared lasers,” IEEE Trans. Microwave Theory Tech. MT22, 988–990 (1980).

Int. J. Infrared Milimeter Waves

D. Dangoisse, P. Glorieux, and J. Wascat, “Diffusion and vibrational bottleneck in optically pumped submillimeter laser,” Int. J. Infrared Milimeter Waves 2, 215–229 (1981). [CrossRef]

J. Chem. Phys.

H. O. Everitt and F. C. De Lucia, “A time-resolved study of rotational energy transfer into A and E symmetry species of 13CH3F,” J. Chem. Phys. 90, 3520–3527 (1989). [CrossRef]
H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in CH3F: The ΔJ = n, ΔK = 0 processes,” J. Chem. Phys. 92, 6480–6491 (1990). [CrossRef]

Opt. Commun.

T. Y. Chang and T. J. Bridges, “Laser actions at 452, 496, and 541 μm in optically pumped CH3F,” Opt. Commun. 1, 423–426 (1970). [CrossRef]
R. J. Temkins and D. R. Cohn, “Rate equations for an optically pumped, far-infrared laser,” Opt. Commun. 16, 213–217 (1976). [CrossRef]

Opt. Express

S. L. 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–1562 (2011). [CrossRef] [PubMed]

Proc. IEEE

M. S. Tobin, “A review of optically pumped NMMW lasers,” Proc. IEEE 73, 61–85 (1985). [CrossRef]

Other

P. K. Cheo (ed.), Handbook of Molecular Lasers (Marcel Dekker, Inc., 1987), pp. 495–569.
H. O. Everitt, “Collisional Energy Transfer in Methyl Halides,” PhD Thesis (Department of Physics, Duke University, 1990).
H. O. Everitt and F. C. De Lucia, “Rotational energy transfer in small polyatomic molecules,” in Advances in Atomic and Molecular Physics (Academic Press, 1995), Vol. 35, pp. 331–400.
L. E. Reichl, A Modern Course in Statistical Physics (John Wiley & Sons Inc., 1998).
A. E. Siegman, Lasers (Univ. Science Books, 1986).
R. Bansal (ed.), Handbook of Engineering Electromagnetics (Marcel Dekker, Inc., 2004). [CrossRef]
A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, 2007).

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