# Patent

## Title
Passively mode locked quantum cascade lasers

## Abstract

This invention relates to a self-induced transparency mode-locked quantum cascade laser having an active section comprising a plurality of quantum well layers deposited in alternating layers on a plurality of quantum barrier layers and form a sequence of alternating gain and absorbing periods, said alternating gain and absorbing periods interleaved along the growth axis of the active section.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
The application claims priority benefit under 35 USC 119(e) to U.S. provisional patent application Ser. No. 61/144,178 filed 13 Jan. 2009, the contents of which are incorporated by reference herein in their entirety.


STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT
Federal government funds were used in researching or developing this invention under National Science Foundation Grant No. EEC-0540832 under auspices of the MIRTHE Engineering Research Center.


NAMES OF PARTIES TO A JOINT RESEARCH AGREEMENT
None.
BACKGROUND
1. Field of the Invention
This invention relates to lasers that emit sub-picosecond pulses, and particularly to semiconductor lasers called quantum cascade lasers that emit in the mid- to far-infrared portion of the electromagnetic spectrum. These pulses have particular importance to applications where materials must be processed while minimizing heat transfer and damage to the materials, for example tissue removal in medicine and dentistry, and micromachining at the pico and femto scales. These pulses also have importance to applications involving spectroscopy, such as detection of trace compounds, clock synchronization in computer networks and telecommunication networks, and in light detection and ranging (lidar).
2. Background of the Invention
Although interband lasers can be high-power lasers, it is difficult to obtain any power at all in interband lasers in the mid-IR. Quantum cascade lasers are semiconductor lasers that emit in the mid- to far-infrared portion of the electromagnetic spectrum. Compared to other semiconductor diode lasers that emit in the mid-far IR, QC lasers have higher output power since laser emission is achieved from the transition of an electron through periodic thin layers of material forming a superlattice that introduces an electric potential over the device. Unlike other lasers, the lattice allows the electron to emit multiple photons as it traverses, or cascades, from one period of the lattice to the next. Semiconductor lattices in QC lasers can be made from layers of crystalline aluminum indium arsenide alternating with indium gallium arsenide, which create structures called quantum wells. One of the interesting features of QC lasers is that they operate at room temperature, without the need for cooling that might be found in non-QC lasers.
A technique for creating extremely short pulses of laser light is called mode-locking. By pumping a laser into a laser cavity consisting of two or more mirrors, the normal random oscillations of the light waves can be made synchronous and thus increase (constructive interference), called mode-locking, or may be made to interfere and damp (destructive interference). By using an electronic device that modulates the light intensity, an electrical signal can be used to establish a short pulse within the laser cavity in a technique called active mode-locking. However, in active mode-locking, pulses that are less than about one picosecond have not been obtained. An actively modelocked QCL that produces 3 ps pulses at 6.2 μm was reported by Wang et al.
In contrast, passive mode-locking does not require an external signal to maintain a continuous stream of pulses. Some examples of passive mode-locking use an intracavity absorber to absorb low-intensity light and preferentially amplify high-intensity spikes. Certain dyes in solution can act as saturable absorbers, but graphite-lattice structures, and the use of an aperture have also been used to focus and amplify the high-intensity light and attenuate the low-intensity light. Aperture transmission of high-intensity light and attenuation of low-intensity is called Kerr lens mode-locking, sometimes called “self mode-locking” However, no semiconductor laser has been passively mode-locked in the mid-IR range.
Laser light can have an interesting relationship to the medium or material that it impacts. McCall and Hahn observed a nonlinear optical propagation effect in resonant medium. According to their observation, a short pulse of coherent light above a critical input energy, for a given pulse width τ<T2, can pass through a saturable resonant medium as though the medium were transparent. However, below the critical energy this self-induced transparency (SIT) phenomenon cannot happen, rather, the pulse energy is absorbed.
Because QCLs typically have relatively short values of T1 (1-10 ps) and relatively long values of T2 ( ̃100 fs) compared to other semiconductor lasers, standard passive modelocking cannot be achieved. However, the conditions are ideal for any approach, denoted SIT modelocking. After the original observation of passive modelocking, it was speculated that passive modelocking was SIT modelocking, but this speculation was later proved to be incorrect. Still later, Kozlov studied the possibility of obtaining SIT modelocking by using gas lasers with two separate tubes in which the gases in the two different tubes have the same resonance lines and the gas in one tube has twice the dipole moment of the other. However, this design cannot be implemented in practice. Applications for mode-locked QC lasers include three-dimensional micro-, or nano-fabrication or machining, investigative and diagnostic spectroscopy, optical communication and electronic devices such as limiters or storage, and medical and dental probes and surgical devices.
BRIEF SUMMARY OF THE INVENTION
To address limitations in the prior art, there is provided in one preferred embodiment a self-induced transparency mode-locked quantum cascade laser, comprising:
(i) an active section comprising a plurality of quantum well layers deposited in alternating layers on a plurality of quantum barrier layers and form a sequence of alternating gain and absorbing periods, said alternating gain and absorbing periods interleaved along the growth axis of the active section, wherein the absorbing periods have a dipole moment of about twice that of the gain periods;

(ii) an optical cavity that houses the active section and permits amplified light to escape;

(iii) an externally supplied seed pulse; and

(iv) current injectors structured and arranged to apply an electric control field to the active section.

In another preferred embodiment, there is also provided wherein the quantum barrier layers are made using Indium-Aluminum-Arsenide, and the quantum well layers are made using Indium-Gallium-Arsenide.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser having a wavelength of between about 3 micrometers and about 12 micrometers, and in some embodiments between about 8 and about 12 micrometers.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide coherence times  gain recovery times  round trip times.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide: coherence times of at least 50 to about 200 femtoseconds, and in one embodiment about 100 femtoseconds; gain recovery times of at least 1 to about 10 picoseconds, and in one embodiment about 1 picosecond; and round trip times of at least 40-60 picoseconds, in one embodiment 50 picoseconds.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide a heterostructure having a ratio of 4 gain periods:1 absorbing period. In other embodiments, the heterostructure may include a ratio of 3 gain periods to 1 absorbing period, or 5 gain periods to 1 absorbing period.
In another preferred embodiment, there is also provided wherein the input pulse is a it pulse in the gain medium.
In another preferred embodiment, there is also provided wherein each gain period and each absorbing period have over 16 quantum layers.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to suppress continuous waves and eliminate spatial hole burning.
In another preferred embodiment, there is also provided wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser having pulse length less than 100 femtoseconds long.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a Schematic of a QCL structure with gain and absorbing periods. On the left, we show a cutaway view of the QCL structure. The active region is shown as a filled-in rectangle. We are looking in the direction along which light would propagate. Electrodes would be affixed to the top and bottom so that electrons flow vertically. The heterostructure would also be stacked vertically as shown on the right. We show one absorbing period for every four gain periods, corresponding schematically to the case in which the electron density in the gain medium Ng≅4× the electron density in the absorbing medium (Na), and we show absorbing periods that are twice as large as gain periods to indicate schematically that the dipole moment in the absorbing medium (μa)≅2× the dipole moment in the gain medium (μg).
FIG. 2. Schematic of the (a) gain and (b) absorbing media. Black straight-line arrows indicate the direction of electron flow. Red wavy arrows indicate radiative transitions.
FIG. 3. Conduction-band diagram and moduli-squared wave functions for one gain and one absorbing period of the (a) 12 μm, (b) 8 μm mode-locking QCL structures. The sequence of layer dimensions is (in Å, starting from left): (a) 37, 36, 10, 35, 10, 34, 11, 34, 12, 35, 39, 37, 12, 62, 14, 58, 28, 42, 12, 40, 13, 37, 15, 34, 15, 34, 34, 45, 11, 65, 6, 69; (b) 42, 34, 9, 33, 12, 30, 13, 28, 16, 28, 41, 27, 18, 62, 14, 58, 28, 42, 12, 40, 13, 37, 13, 34, 16, 34, 34, 9, 31, 50, 5, 84. The numbers in bold type indicate In0.52Al0.48As barrier layers and those in roman type indicate In0.53Ga0.47As well layers. Red wavy arrows indicate radiative transitions.
FIG. 4. Stability limits of the normalized absorption (ā) vs the normalized gain ( g) coefficients with different Δg0 and Δa0. The ratio T1/T2 is infinity in all cases. For a given τ and ā, the required g increases as Δg0 and |Δa0| decrease.
FIG. 5. Stability limits of the normalized absorption (ā) vs the normalized gain ( g) coefficients with different T2a/T2g. We set T1g=T1a=∞ in all cases. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0, and the absorbing medium is completely uninverted, i.e., Δa0=−1.0.
FIG. 6
Stability limits of the normalized absorption vs the normalized gain coefficients with different T1/T2. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0 and the absorbing medium is completely uninverted, i.e., Δa0=−1.0. In each bundle of dashed lines, corresponding to a fixed value of τ, T1/T2 decreases from left to right.
FIG. 7
Pulse evolution in the system. (a) Gain and absorption coefficients are in the stable regime, g=4.0, ā=3.5. (b) Gain and absorption coefficients are in the regime in which continuous waves are unstable, g=4.0, ā=1.0. (c) Gain and absorption coefficients are in the regime where any initial pulse attenuates, g=4.0, ā=7.8. The ratio T1/T2 equals 10 in all cases. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0, and the absorbing medium is completely uninverted, i.e., Δa0=−1.0.
FIG. 8
Stability limits of the normalized absorption (ā) vs the normalized gain ( g) coefficients with different values of T1a/T1g. We set T2g=T2a and T1g/T2g=10 in all cases. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0, and the absorbing medium is completely uninverted, i.e., Δa0=−1.0
FIG. 9
Stability limits of the ratio of the dipole moments in the absorbing and gain media ( μ) vs the normalized gain coefficient ( g) for three cases of normalized absorption (ā). The ratio T1/T2 is 10 in all cases. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0, and the absorbing medium is completely uninverted, i.e., Δa0=−1.0
FIG. 10
Input pulse energy limits vs normalized input pulse duration (τi/T2) for two different cases of g and ā. In both the cases, we set T1/T2=10. In equilibrium, the gain medium is completely inverted, i.e., Δg0=1.0, and the absorbing medium is completely uninverted, i.e., Δa0=−1.0.



DETAILED DESCRIPTION OF THE INVENTION
Quantum cascade lasers (QCLs) are important light sources in the midinfrared range, about 3 μm to about 12 μm. The light is generated by a transition between two subbands in the conduction band, in contrast to interband semiconductor lasers. As a consequence, the subbands have narrow linewidths and long coherence times T2 compared to interband semiconductor lasers. Values of T2 on the order of 100 fs are achievable. Another important feature of the QCLs is their rapid gain recovery times T1 compared to interband semiconductor lasers due to fast LO-phonon relaxation. Typical values of T1 are in the range 1-10 ps, which is short compared to Trt, the round-trip time in the cavity. Typical values of Trt are on the order of 50 ps. The narrow linewidths and fast recovery times of QCLs make it difficult to achieve conventional passive mode locking. Gain bandwidths that are significantly larger than the pulse bandwidths are required, and that is hard to obtain when the linewidths are narrow, as in QCLs. A saturable gain, with a recovery time that is long compared to Trt, is also required for conventional passive mode locking in order to suppress continuous waves, and the typical gain recovery times in QCLs are shorter than the round-trip times. Thus, conventional passive mode locking cannot work in QCLs that operate in a standard parameter regime. However, QCLs are an ideal tool for creating the long-predicted, never-observed SIT mode locking. The relationship T2T1Trt, which is typical for QCLs, is precisely what is needed for SIT mode locking, and the ease of band-gap engineering in QCLs makes it possible to interleave gain and absorbing periods as needed. Conversely, SIT mode locking of QCLs makes it possible to obtain mode-locked pulses from a mid-infrared laser that are less than 100 fs in duration. Using SIT mode-locking in a QC laser also addresses the stability problems relating to the growth of continuous waves and spatial hole burning. In SIT, continuous waves are suppressed and mode-locked pulses propagate in one direction or the other at a time so that spatial hole burning is not an issue.
Referring now to the FIGURES, specific non-limiting examples and preferred embodiments of SIT mode-locking structures for QC lasers are described.
In order to obtain SIT mode locking, it is necessary to have two highly coherent resonant media with nearly equal resonant frequencies. In one medium, denoted the gain medium, electrons should be injected into the upper lasing state so that the resonant states are nearly inverted. In the other medium, denoted the absorbing medium, electrons should be injected into the lower state so that the resonant states are not inverted. Also, the dipole strength in the absorbing medium should be nearly equal to twice the dipole strength in the gain medium. At the same time, the ratio of the gain per unit length to the absorption per unit length should be small enough so that the growth of continuous waves is suppressed, but large enough so that a mode-locked pulse can stably exist. It is possible to simultaneously satisfy all these conditions by interleaving gain and absorbing periods that have the required dipole strengths as shown schematically in FIG. 1.
FIG. 1 shows a schematic of a QCL structure with gain and absorbing periods. On the left, we show a cutaway view of the QCL structure. The active region is shown as a filled-in rectangle. We are looking in the direction along which light would propagate. Electrodes would be affixed to the top and bottom so that electrons flow vertically. The heterostructure would also be stacked vertically as shown on the right. We show one absorbing period for every four gain periods, corresponding schematically to the case in which the electron density in the gain medium (Ng)≅4× the electron density in the absorbing medium (Na), and we show absorbing periods that are twice as large as gain periods to indicate schematically that the dipole moment in the absorbing medium (μa)≅2× the dipole moment in the gain medium (μg). By appropriately choosing the number of gain periods and the number of absorbing periods, one can obtain any desired ratio for the gain and absorption per unit length. As long as there are many periods within the transverse wavelength of the lasing mode, the gain and absorbing periods will experience the same light intensity.
In FIG. 2, we show simplified two-level resonant structures for the gain and absorbing media. FIG. 2 shows a schematic of the (a) gain and (b) absorbing media. Black straight-line arrows indicate the direction of electron flow. Red wavy arrows indicate radiative transitions. In the gain medium, electrons are injected into level 2g and are extracted from level 1g. The carrier lifetime in 2g should be longer than the mode-locked pulse duration and the equilibrium population inversion should be nearly complete. When an optical pulse with a photon energy equal to the resonant energy impinges on the gain medium with its polarization oriented in the direction perpendicular to the layers, electrons scatter to level 1g and photons are emitted. Then, the electrons are nonradiatively extracted from level 1g. In the absorbing medium, electrons are injected into the lower level 1a. The lifetime of state 1a should again be longer than the pulse width. When a light pulse of the appropriate wavelength and polarization impinges on the absorbing medium, photons are absorbed and electrons jump to level 2a. If a light pulse has enough intensity, then photons are re-emitted with no overall loss in one Rabi oscillation time. In order for these processes in the gain medium and absorbing medium to occur simultaneously, the energy spacing between the resonant levels should be nearly the same in both media.
In the theory of resonant two-level media, both π pulses and 2π pulses play an important role. A π pulse is a pulse with sufficient energy to exactly invert the lower state population of a two-level medium if the medium is initially uninverted, or, conversely, to uninvert the upper state if the medium is initially inverted. In the former case, the pulse experiences loss and rapidly attenuates, but, in the latter case, the pulse experiences gain. The pulse duration is approximately half a Rabi oscillation period. If a pulse lasts a longer time than required to drive the population from the upper level to the lower, then the medium will amplify the first part of the pulse and attenuate the latter part, in a way that shortens the pulse. Conversely, if a pulse is initially too short, it is lengthened. Because a π pulse experiences gain, it is natural that shortly after the initial observations of passive mode locking in lasers, it was proposed that the pulses in these lasers are actually SIT-induced π pulses. However, these pulses are not suitable for use on their own as passively mode-locked laser pulses. Where one π pulse can exist, there is nothing to prevent continuous waves from generating multiple pulses leading to chaos rather than a single stable pulse oscillating in a laser cavity.
One can circumvent this difficulty by combining a gain medium in which the optical pulse is a π pulse with an absorbing medium in which the optical pulse is a 2π pulse. A 2π pulse is a pulse with sufficient energy so that in an uninverted medium the lower state population is first inverted and then returned to the lower state in approximately one Rabi oscillation. A 2π pulse like a π pulse is stable. If its initial duration is too long, the duration decreases, and, if its initial duration is too short, the duration increases. The 2π pulse propagates through the medium without loss, in contrast to continuous waves at the resonant optical frequency, which experience loss. This remarkable property is what led to the name “self-induced transparency.” Because of this property, the absorbing medium acts like the saturable loss in a conventional passively modelocked system, suppressing the growth of continuous waves, while allowing a short pulse to propagate.
It is evidently important that both the gain medium and the absorbing medium act on the optical pulse simultaneously. We may achieve this simultaneous interaction by designing a QCL structure that has the gain and absorbing periods interleaved along the growth axis of the structure, as shown in FIG. 1. By making the dipole moment in the absorbing periods twice that of the gain periods, a π pulse in the gain periods is a 2π pulse in the absorbing periods. Therefore, an injected π pulse completely depletes the gain medium as it makes its way through the QCL, whereas, the absorbing medium becomes transparent. We will show that by controlling the amount of gain and absorption per unit length in the gain and absorbing media, pulse durations can be controlled.
In order to suppress spatial hole burning, the RNGH instability, or the growth of multiple pulses, we do not want continuous waves to grow in an SIT mode-locked laser. The absorption parameter should be large enough to keep the laser operating below the threshold for the growth of continuous waves. Therefore, the laser cannot self-start and it is necessary to use some means to start the mode locking. Essentially, we need to create an optical pulse that has sufficient energy and a duration on the order of T1. This pulse can be created by using an electrical shock, electrical pulse, a mechanical vibration or an external optical seed. We suggest two optical approaches. First, we can seed the pulse from an external source by injection locking, or, second, we can use active mode locking to generate an initial pulse that will have a suitable energy and initial duration for SIT mode locking. It may also be possible to use a mechanical or an electrical impulse to start the mode locking.
Quantum Cascade Laser Structures
We have designed QCL gain and absorbing periods that fulfill the SIT mode locking requirements at two different wavelengths, 8 and 12 μm. Similar structures can be designed over a broad range of mid-IR wavelengths, from about 3 μm to about 12 μm. We use the In0.52Al0.48As/In0.53Ga0.47As material system for the active region in our design since most of the demonstrated QCLs to date have been based on this material system. However, mode-locking structures operating at less than 8 μm will be difficult to design using this material system. Since electrons are injected into the lower state in the absorbing medium, the upper state approaches close enough to the conduction-band edge at wavelengths below 8 μm to lead to a large increase in the carrier loss due to scattering to the continuum-like states. To design a mode-locking structure at shorter wavelengths, it should be possible to use a strain-balanced InAlAs/InGaAs material system.
The gain periods in our design are typical QCL periods. We design a three-quantum-well active region for the gain periods that has a diagonal transition, which lowers the dipole moment relative to designs that have a vertical transition. This choice simplifies the design of the absorbing periods. Population inversion is achieved by confining the resonant states in separate quantum wells and depopulating the lower state through phonon relaxation to another state sitting below the lower resonant state. While the dipole moment in the gain periods is not high, the upper state lifetime is larger than in the case of vertical transitions so that the gain remains high. The design of absorbing periods is different from the design of gain periods. The combined requirements of carrier injection into and extraction from the lower resonant level and a dipole moment twice that of the gain periods make it difficult to design the absorbing periods. To achieve a large dipole moment, the transition should be between two excited states. The carrier lifetime is made high by reducing scattering through phonon relaxation and reducing the carrier tunneling from the lower resonant state into the succeeding injector states. The injector regions have different designs when the electrons are tunneling into a gain or absorbing active region due to the different quantum electronic structures of gain and absorbing active regions.
The structure in FIG. 3(a) emits light at 12 μm. FIG. 3. Conduction-band diagram and moduli-squared wave functions for one gain and one absorbing period of the (a) 12 μm, (b) 8 μm mode-locking QCL structures. The sequence of layer dimensions is (in Å, starting from left): (a) 37, 36, 10, 35, 10, 34, 11, 34, 12, 35, 39, 37, 12, 62, 14, 58, 28, 42, 12, 40, 13, 37, 15, 34, 15, 34, 34, 45, 11, 65, 6, 69; (b) 42, 34, 9, 33, 12, 30, 13, 28, 16, 28, 41, 27, 18, 62, 14, 58, 28, 42, 12, 40, 13, 37, 13, 34, 16, 34, 34, 9, 31, 50, 5, 84. The numbers in bold type indicate In0.52Al0.48As barrier layers and those in roman type indicate In0.53Ga0.47As well layers. Red wavy arrows indicate radiative transitions.
Electrons are injected into level 3g, and the gain transition is between levels 3g and 2g. The dipole moment between the resonant levels is given by μg/e=1.81 nm. Level 3g has a lifetime of  ̃2 ps. Level 1g is positioned approximately at phonon resonance with 2g. Level 2g has a lifetime of  ̃0.5 ps, so that the population inversion is high. The transition energy is 101 meV. In the absorbing periods, electrons are injected into level 4a and they jump to level 5a by absorbing photons. The lifetime of level 5a is  ̃0.83 ps. The absorbing levels are separated by 101 meV and have a dipole moment μa/e=3.65 nm. After electrons are injected into level 3g, they have a long lifetime of 2 ps due to the low scattering rate. When the optical pulse arrives, the population in level 3g decreases to level 2g. Level 2g is depopulated quickly to level 1g through phonon interactions and the electrons then transit to the following injector stage. Population inversion is restored before the optical pulse makes a round trip in the laser cavity. Electrons sit in level 4a after being injected by the preceding injector stage. When an optical pulse reaches the absorbing periods, electrons from level 4a move to level 5a by absorbing photons. Since the pulse is a 2π pulse for the absorbing medium, level 5a is depleted by making a Rabi oscillation, during which photons are emitted.
In the structure given in FIG. 3(b), the gain transition is between levels 3g and 2g. The dipole moment between the gain levels is 1.55 nm and the lifetime of level 3g is −3 ps. The resonant transition energy is 150 meV in both the gain and absorbing periods. In this structure, electrons are injected into level 3a in the absorbing periods and the absorbing transition is between levels 3a and 4a. Level 4a has a lifetime of  ̃0.8 ps. The absorbing levels have a dipole moment μa/e=2.95 nm.
In a QCL, the time T1 is determined mainly by the LO-phonon relaxation rate. The LO-phonon relaxation rate depends mainly on the energy spacing between the levels and the overlap of the corresponding wave functions, so that T1 depends on the details of the band structure and can vary greatly from design to design. Indeed, strictly speaking there is not a single T1 in either the gain or the absorbing medium since there are more than two levels involved in the dynamics in both media. In the gain stage of the QCL structure given in FIG. 3, the upper state is confined in a separate quantum well from the other states in the active region, so that the phonon relaxation rate is smaller than when all the active states are in the same well and the lifetime is higher. However, the absorbing stage is designed such that the upper and lower state wave functions have a large overlap, which makes the dipole moment higher, but decreases the lifetime. Therefore, in practical designs, we find T1g>T1a.
Generally, if the gain and absorbing media are grown from the same material system, it is reasonable to assume T2g≈T2a. In the QCL structure that we propose, an In0.52Al0.48As/In0.53Ga0.47As material system is used for both the gain and absorbing periods, so that T2a should not vary much from T2g. However, there has yet to be a detailed theoretical calculation of these coherence times. In a QCL, the value of T2 depends mainly on electron-electron scattering, electron-LO phonon scattering, and interface scattering. Therefore, the values of T2g and T2a may differ somewhat, depending on the details of the design.
Maxwell-Bloch Equations
Wang et al. and Gordon et al. have observed the RNGH instability in QCLs with only gain periods. They showed evidence for Rabi oscillations and demonstrated that the two-level Maxwell-Bloch equations apply to QCLs in some parameter regimes, although they also showed that saturable absorption affects the behavior quantitatively, significantly reducing the RNGH threshold. Gordon et al. attributed the saturable absorption to Kerr lensing that increases the mode overlap with the active region and reduces the wall losses. These effects depend sensitively on the details of the QCL geometry. They also observed that spatial hole burning due to the interaction of forward- and backward-propagating waves has an important effect on the pulse spectrum. They did not find it necessary to include chromatic dispersion or other nonlinearities. Motivated by these results, we use a simple two-level model based on the standard one-dimensional Maxwell-Bloch equations. The Maxwell-Bloch equations that describe the light propagation and light-matter interaction in QCL having interleaved gain and absorbing periods may be written as




  
    
      
        
          
            
              n
              c
            
            ⁢
            
              
                ∂
                E
              
              
                ∂
                t
              
            
          
          =
          
            
              -
              
                
                  ∂
                  E
                
                
                  ∂
                  z
                
              
            
            -
            
              i
              ⁢
              
                
                  
                    kN
                    g
                  
                  ⁢
                  
                    Γ
                    g
                  
                  ⁢
                  
                    μ
                    g
                  
                
                
                  2
                  ⁢
                  
                    ɛ
                    0
                  
                  ⁢
                  
                    n
                    2
                  
                
              
              ⁢
              
                η
                g
              
            
            -
            
              i
              ⁢
              
                
                  
                    kN
                    a
                  
                  ⁢
                  
                    Γ
                    a
                  
                  ⁢
                  
                    μ
                    a
                  
                
                
                  2
                  ⁢
                  
                    ɛ
                    0
                  
                  ⁢
                  
                    n
                    2
                  
                
              
              ⁢
              
                η
                a
              
            
            -
            
              
                1
                2
              
              ⁢
              lE
            
          
        
        ,
      
    
    
      
        (
        
          1
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            ∂
            
              η
              g
            
          
          
            ∂
            t
          
        
        =
        
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  g
                
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            
              Δ
              g
            
            ⁢
            E
          
          -
          
            
              η
              g
            
            
              T
              
                2
                ⁢
                g
              
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            ∂
            
              Δ
              g
            
          
          
            ∂
            t
          
        
        =
        
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  g
                
              
              ħ
            
            ⁢
            
              η
              g
            
            ⁢
            
              E
              *
            
          
          -
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  g
                
              
              ħ
            
            ⁢
            
              η
              g
              *
            
            ⁢
            E
          
          +
          
            
              
                Δ
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              -
              
                Δ
                g
              
            
            
              T
              
                1
                ⁢
                g
              
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            ∂
            
              η
              a
            
          
          
            ∂
            t
          
        
        =
        
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  a
                
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            
              Δ
              a
            
            ⁢
            E
          
          -
          
            
              η
              a
            
            
              T
              
                2
                ⁢
                a
              
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            ∂
            
              Δ
              a
            
          
          
            ∂
            t
          
        
        =
        
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  a
                
              
              ħ
            
            ⁢
            
              η
              a
            
            ⁢
            
              E
              *
            
          
          -
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  a
                
              
              ħ
            
            ⁢
            
              η
              a
              *
            
            ⁢
            E
          
          +
          
            
              
                Δ
                
                  a
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              -
              
                Δ
                a
              
            
            
              T
              
                1
                ⁢
                a
              
            
          
        
      
    
    
      
        (
        
          1
          ⁢
          e
        
        )
      
    
  




where the subscripts g and a in Eq. (1) represent gain and absorption, respectively. The independent variables z and t are length along the light-propagation axis of the QCL and time. The dependent variables E,ηg,a, and Δg,a refer to the envelopes of the electric field, gain polarization, and gain inversion. The parameters Δg0 and Δa0 refer to the equilibrium inversion away from the mode-locked pulse. The parameters n and c denote the index of refraction and the speed of light. The parameters Ng,aΓg,a, denote the effective electron density multiplied by the overlap factor. The parameters k, l, ∈0, and h denote the wave number in the active region, the linear loss including the mirror loss, the vacuum dielectric permittivity, and Planck's constant. The notation closely follows that of Wang et al., with the differences that we have an absorbing as well as a gain medium, and we are considering unidirectional propagation, as is appropriate for a mode-locked pulse.
In order to achieve SIT mode locking, the growth of continuous waves must be suppressed. At this point, two observations are made. First, because continuous waves are suppressed, forward- and backward-propagating waves cannot interact, and spatial hole burning will not occur. Second, we did not include saturable absorption in Eq. (1) because we are not certain how to do so. This contribution was added phenomenologically to the Maxwell-Bloch equations by Wang et al. and Gordon et al., based on experimental observations in particular QCLs and was attributed to effects that depend sensitively on the geometry of those QCLs.
In order to suppress continuous waves, the gain must be below threshold. To derive this condition, we set Δg=Δg0 and Δa=Δa0 in steady state, where there is no evolution in z. We then find from Eqs. (1b) and (1d),




  
    
      
        
          
            η
            g
          
          =
          
            i
            ⁢
            
              
                μ
                g
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            
              T
              
                2
                ⁢
                g
              
            
            ⁢
            
              Δ
              
                g
                ⁢
                
                  
                
                ⁢
                0
              
            
            ⁢
            E
          
        
        ,
        
          

        
        ⁢
        
          
            η
            a
          
          =
          
            i
            ⁢
            
              
                μ
                a
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            
              T
              
                2
                ⁢
                a
              
            
            ⁢
            
              Δ
              
                a
                ⁢
                
                  
                
                ⁢
                0
              
            
            ⁢
            E
          
        
        ,
      
    
    
      
        (
        2
        )
      
    
  




where we are considering continuous waves, so that there is no dependence on t and the t-derivatives vanish. After substitution in Eq. 1a, we obtain in steady state




  
    
      
        
          
            
              
                
                  kN
                  g
                
                ⁢
                
                  Γ
                  g
                
                ⁢
                
                  μ
                  g
                  2
                
                ⁢
                
                  T
                  
                    2
                    ⁢
                    g
                  
                
                ⁢
                
                  Δ
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
              
              
                2
                ⁢
                
                  ɛ
                  0
                
                ⁢
                
                  n
                  2
                
                ⁢
                ħ
              
            
            +
            
              
                
                  kN
                  a
                
                ⁢
                
                  Γ
                  a
                
                ⁢
                
                  μ
                  a
                  2
                
                ⁢
                
                  T
                  
                    2
                    ⁢
                    a
                  
                
                ⁢
                
                  Δ
                  
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
              
              
                2
                ⁢
                
                  ɛ
                  0
                
                ⁢
                
                  n
                  2
                
                ⁢
                ħ
              
            
            -
            l
          
          =
          0
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




which may also be written as gΔg0+aΔa0−l=0, where





  
    
      
        
          g
          =
          
            
              
                kN
                g
              
              ⁢
              
                Γ
                g
              
              ⁢
              
                μ
                g
                2
              
              ⁢
              
                T
                
                  2
                  ⁢
                  g
                
              
            
            
              2
              ⁢
              
                ɛ
                0
              
              ⁢
              
                n
                2
              
              ⁢
              ħ
            
          
        
        ,
        
          
        
        ⁢
        
          a
          =
          
            
              
                
                  kN
                  a
                
                ⁢
                
                  Γ
                  a
                
                ⁢
                
                  μ
                  a
                  2
                
                ⁢
                
                  T
                  
                    2
                    ⁢
                    a
                  
                
              
              
                2
                ⁢
                
                  ɛ
                  0
                
                ⁢
                
                  n
                  2
                
                ⁢
                ħ
              
            
            .
          
        
      
    
    
      
        (
        4
        )
      
    
  




Physically, the parameter g corresponds to the gain per unit length from the gain periods of the QCL and the parameter a corresponds to the absorption per unit length from the absorbing periods. The condition for the linear gain to remain below threshold is gΔg0+aΔa0−l<0. In the case of a fully inverted gain medium, so that Δg=Δg0=1 and a fully uninverted absorbing medium so that Δa=Δa0=−1, the condition to suppress continuous waves becomes g−a−1<0.
Assuming that T1g and T1a are large enough so that they may be set equal to co in Eq. 1, and focusing on the special case in which μa=2μg and the pulse is a it pulse in the gain medium, Eq. 1 has an exact analytical solution that we may write as




  
    
      
        
          E
          =
          
            
              ħ
              
                
                  μ
                  g
                
                ⁢
                τ
              
            
            ⁢
            sech
            ⁢
            
              
            
            ⁢
            x
          
        
        ,
      
    
    
      
        (
        
          5
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            η
            g
          
          =
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  B
                  g
                
              
              2
            
            ⁢
            
              Δ
              
                g
                ⁢
                
                  
                
                ⁢
                0
              
            
            ⁢
            sech
            ⁢
            
              
            
            ⁢
            x
          
        
        ,
      
    
    
      
        (
        
          5
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            Δ
            g
          
          =
          
            
              B
              g
            
            ⁢
            
              
                Δ
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              ⁡
              
                (
                
                  
                    τ
                    
                      T
                      
                        2
                        ⁢
                        
                          
                        
                        ⁢
                        g
                      
                    
                  
                  -
                  
                    tanh
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          5
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            η
            a
          
          =
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  B
                  a
                
              
              2
            
            ⁢
            
              
                Δ
                
                  a
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              ⁡
              
                (
                
                  
                    
                      -
                      sech
                    
                    ⁢
                    
                      
                    
                    ⁢
                    x
                    ⁢
                    
                      
                    
                    ⁢
                    tanh
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                  +
                  
                    
                      τ
                      
                        3
                        ⁢
                        
                          T
                          
                            2
                            ⁢
                            a
                          
                        
                      
                    
                    ⁢
                    sech
                    ⁢
                    
                      
                    
                    ⁢
                    x
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          5
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            Δ
            a
          
          =
          
            
              
                
                  B
                  a
                
                2
              
              ⁢
              
                
                  Δ
                  
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
                ⁡
                
                  (
                  
                    1
                    +
                    
                      
                        τ
                        2
                      
                      
                        3
                        ⁢
                        
                          T
                          
                            2
                            ⁢
                            a
                          
                          2
                        
                      
                    
                  
                  )
                
              
            
            -
            
              
                B
                a
              
              ⁢
              
                
                  Δ
                  
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
                ⁡
                
                  (
                  
                    
                      
                        sech
                        2
                      
                      ⁢
                      x
                    
                    +
                    
                      
                        
                          2
                          ⁢
                          τ
                        
                        
                          3
                          ⁢
                          
                            T
                            
                              2
                              ⁢
                              a
                            
                          
                        
                      
                      ⁢
                      tanh
                      ⁢
                      
                        
                      
                      ⁢
                      x
                    
                  
                  )
                
              
            
          
        
        ,
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          5
          ⁢
          
            
          
          ⁢
          e
        
        )
      
    
  
  
    
      
        
          x
          =
          
            
              t
              τ
            
            -
            
              z
              
                v
                ⁢
                
                  
                
                ⁢
                τ
              
            
          
        
        ,
        
          

        
        ⁢
        while
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            B
            g
          
          =
          
            1
            
              1
              +
              
                τ
                /
                
                  T
                  
                    2
                    ⁢
                    g
                  
                
              
            
          
        
        ,
        
          
        
        ⁢
        
          
            and
            ⁢
            
              
            
            ⁢
            
              B
              a
            
          
          =
          
            2
            
              
                (
                
                  1
                  +
                  
                    
                      τ
                      /
                      3
                    
                    ⁢
                    
                      T
                      
                        2
                        ⁢
                        a
                      
                    
                  
                
                )
              
              ⁢
              
                (
                
                  1
                  +
                  
                    τ
                    /
                    
                      T
                      
                        2
                        ⁢
                        a
                      
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




are chosen such that Δg→Δg0 and Δa→Δa0 as t→−∞, so that the equilibrium population completely recovers on every pass of the pulse through the laser. Hence, Eq. 5 shows that stable mode-locked operation can be achieved in the proposed structure. The parameters τ and ν that correspond to the pulse duration and the pulse velocity are determined by the equations




  
    
      
        
          
            
              
                
                  τ
                  /
                  
                    T
                    
                      2
                      ⁢
                      g
                    
                  
                
                
                  1
                  +
                  
                    τ
                    /
                    
                      T
                      
                        2
                        ⁢
                        g
                      
                    
                  
                
              
              ⁢
              g
              ⁢
              
                
              
              ⁢
              
                Δ
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
            
            +
            
              
                
                  
                    
                      τ
                      2
                    
                    /
                    3
                  
                  ⁢
                  
                    T
                    
                      2
                      ⁢
                      a
                    
                    2
                  
                
                
                  
                    (
                    
                      1
                      +
                      
                        τ
                        /
                        
                          T
                          
                            2
                            ⁢
                            a
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    (
                    
                      1
                      +
                      
                        
                          τ
                          /
                          3
                        
                        ⁢
                        
                          T
                          
                            2
                            ⁢
                            a
                          
                        
                      
                    
                    )
                  
                
              
              ⁢
              a
              ⁢
              
                
              
              ⁢
              
                Δ
                
                  a
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
            
            -
            l
          
          =
          0
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          1
          v
        
        =
        
          
            n
            c
          
          -
          
            a
            ⁢
            
              
            
            ⁢
            τ
            ⁢
            
              
                
                  τ
                  /
                  2
                
                ⁢
                
                  T
                  
                    2
                    ⁢
                    g
                  
                
              
              
                
                  (
                  
                    1
                    +
                    
                      τ
                      /
                      
                        T
                        
                          2
                          ⁢
                          a
                        
                      
                    
                  
                  )
                
                ⁢
                
                  (
                  
                    1
                    +
                    
                      
                        τ
                        /
                        3
                      
                      ⁢
                      
                        T
                        
                          2
                          ⁢
                          a
                        
                      
                    
                  
                  )
                
              
            
            ⁢
            
              
                Δ
                
                  a
                  ⁢
                  
                    
                  
                  ⁢
                  0
                
              
              .
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




The full-width half-maximum pulse duration τFWHM equals 1.763τ. This solution was previously presented in the special case Δg0=1.0 and Δa0=−1.0. We now consider in more detail the special case T2g=T2a≡T2. Writing g=g/l, ā=a/l, and τ=τ/T2, we find that the condition to suppress the growth of continuous waves becomes gΔg0+āΔa0−1<0, and the equation for the pulse duration becomes 10, and 11




  
    
      
        
          3
          
            τ
            _
          
        
        =
        
          
            
              
                3
                ⁢
                
                  g
                  _
                
                ⁢
                
                  
                
                ⁢
                
                  Δ
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
              
              -
              4
            
            2
          
          +
          
            
              
                [
                
                  
                    
                      (
                      
                        
                          
                            3
                            ⁢
                            
                              g
                              _
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              Δ
                              
                                g
                                ⁢
                                
                                  
                                
                                ⁢
                                0
                              
                            
                          
                          -
                          4
                        
                        2
                      
                      )
                    
                    2
                  
                  +
                  
                    3
                    ⁢
                    
                      (
                      
                        
                          
                            g
                            _
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            Δ
                            
                              g
                              ⁢
                              
                                
                              
                              ⁢
                              0
                            
                          
                        
                        +
                        
                          
                            a
                            _
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            Δ
                            
                              a
                              ⁢
                              
                                
                              
                              ⁢
                              0
                            
                          
                        
                        -
                        1
                      
                      )
                    
                  
                
                ]
              
              
                1
                /
                2
              
            
            .
          
        
      
    
    
      
        (
        10
        )
      
    
  




Equation (10) only has a solution when ā<(3 gΔg0−2)2/12|Δa0|, where we note that Δa0<0. The conditions for stability may be summarized as





  
    
      
        
          
            (
            
              
                
                  g
                  _
                
                ⁢
                
                  
                
                ⁢
                
                  Δ
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
              
              -
              1
            
            )
          
          
            
            
              Δ
              
                a
                ⁢
                
                  
                
                ⁢
                0
              
            
            
          
        
        <
        
          a
          _
        
        <
        
          
            
              
                (
                
                  
                    3
                    ⁢
                    
                      g
                      _
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      Δ
                      
                        g
                        ⁢
                        
                          
                        
                        ⁢
                        0
                      
                    
                  
                  -
                  2
                
                )
              
              2
            
            
              12
              ⁢
              
                
                
                  Δ
                  
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
                
              
            
          
          .
        
      
    
    
      
        (
        11
        )
      
    
  




When ā is above the upper limit in Eq. (11), we have found by solving Eq. (1) computationally that an initial pulse damps away. When ā is below the lower limit, continuous waves grow. We have computationally found that multiple pulses are generated in this case.
Equation (11) defines a parameter regime in which stable mode-locked operation is possible. In FIG. 4, we present the stability limits when the population inversion in the gain and absorbing periods vary. In all cases, the lower lines indicate the limiting values for ā, below which continuous waves grow, and the upper lines indicate the limiting values for ā, above which initial pulses damp. FIG. 4 shows that the minimum value of g that is required for stable operation increases when Δg0 decreases and Δa0 increases by the same amount. There is also a slight decrease in the lower limit for ā and a larger decrease in the upper limit. Since the upper limit drops more than the lower limit, the stable parameter region becomes smaller. We also show contours of the pulse duration, normalized by the coherence time T2g, denoted τ, as given by Eq. (10). Pulse durations are approximately on the order of T2g when g≈2.5 and ā≈2.0. We also note that pulses become shorter as g and ā increase. However, both g and ā are directly proportional to the current; so, to increase the gain and absorption in a fixed ratio, one must increase the current. At the same time, we note that g and ā are directly proportional to T2. Hence, it is possible to reduce the required current by increasing T2.
We have studied what happens to the stability limits if T2a/T2g vary and we show the results in FIG. 5. In FIG. 5, we have varied T2a keeping T2g constant. In a QCL, typical values of T2g and T2a are on the order of 100 fs. A change in T2a affects the stability limits more than does a change in T2g as is evident from Eq. 8. When T2a/T2g increases, the upper stability limits increase. When g=4.0, we find that the upper limit for ā varies from 3.75 to 8.3 to 24 as T2a/T2g varies from 0.5 to 1.0 to 2.0. The lower limit for ā remains unchanged.
We now derive an energy-balance equation that describes the energy input limits for stable operation when τT2. We define Θ(z,t)=∫−∞tE(z,t′)dt′. Then, Eqs. (1b) and (1c) can be written as




  
    
      
        
          
            
              ∂
              
                
                  η
                  g
                
                ⁡
                
                  (
                  
                    z
                    ,
                    t
                  
                  )
                
              
            
            
              ∂
              t
            
          
          =
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  g
                
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            
              
                Δ
                g
              
              ⁡
              
                (
                
                  z
                  ,
                  t
                
                )
              
            
            ⁢
            
              
                ∂
                
                  Θ
                  ⁡
                  
                    (
                    
                      z
                      ,
                      t
                    
                    )
                  
                
              
              
                ∂
                t
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        12
        )
      
    
  
  
    
      
        
          
            ∂
            
              
                Δ
                g
              
              ⁡
              
                (
                
                  z
                  ,
                  t
                
                )
              
            
          
          
            ∂
            t
          
        
        =
        
          2
          ⁢
          
            
              i
              ⁢
              
                
              
              ⁢
              
                μ
                g
              
            
            ħ
          
          ⁢
          
            
              η
              g
            
            ⁡
            
              (
              
                z
                ,
                t
              
              )
            
          
          ⁢
          
            
              
                ∂
                
                  Θ
                  ⁡
                  
                    (
                    
                      z
                      ,
                      t
                    
                    )
                  
                
              
              
                ∂
                t
              
            
            .
          
        
      
    
    
      
        (
        13
        )
      
    
  




In the gain medium, the polarization and population inversion can be written in terms of a single angle α as 2iηg=cos α and Δg=sin α. We integrate both sides of Eq. (12) or Eq. (13), after substituting these expressions for ηg and Δg and assuming that Δg(z,t→−∞)=1. We then obtain α(z,t)=π/2+(μg/ h)Θ(z,t). We may similarly write 2iηa=cos β and Δa=sin β in the absorbing medium, and we then find β(z,t)=−π/2+(μa/ h)Θ(z,t), where we have set Δa(z,t→−∞)=−1. We now consider Eq. (1a) in steady state, where there is no evolution in z, and in the limit t→∞, where there is no evolution in t. We also define a normalized pulse energy

 θ(z)=(μg/ h)Θ(z,t→∞).  (14)

Equation (1a) now becomes





  
    
      
        
          
            g
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                [
                
                  
                    θ
                    _
                  
                  ⁡
                  
                    (
                    z
                    )
                  
                
                ]
              
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              
                μ
                g
              
              
                μ
                a
              
            
            ⁢
            
              
                T
                
                  2
                  ⁢
                  g
                
              
              
                T
                
                  2
                  ⁢
                  a
                
              
            
            ⁢
            
              sin
              ⁡
              
                [
                
                  
                    
                      μ
                      a
                    
                    
                      μ
                      g
                    
                  
                  ⁢
                  
                    
                      θ
                      _
                    
                    ⁡
                    
                      (
                      z
                      )
                    
                  
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        15
        )
      
    
  




where we note that the linear loss may be neglected in the limit τ<<T2. In the special case μa=2μg, and T2g=T2a≡T2, we find a(μg/μa)(T2g/T2a)sin [(μa/μg) θ]=a cos θ sin θ, so that Eq. (15) becomes cos θ=g/a, which defines the limits of the input energy that is required to generate a single pulse,

cos−1(g/a)< θ<2π−cos−1(g/a).  (16)

When the initial value of θ is within these limits, a single pulse with a final value of θ=π is generated. When the initial value of θ is below this value, the lower limit in Eq. (16), the initial pulse damps. When the initial value of θ is above the upper limit, the initial pulse splits into multiple pulses.
In the analysis up to now it has been assumed that the central carrier frequency of the light pulse and the transition frequency in both the gain and absorbing media are the same. Since the frequency of the light is largely determined by the gain medium, it is reasonable to assume that there is no detuning between the light and the gain medium. Even if the mode locking is seeded by injection locking, analogous to the experiment of Choi et al. the injection-locking laser can be tuned to the gain resonance. There may be a small detuning between the gain and absorbing media due to design or growth issues; however, it is possible to design the gain and absorbing media so that detuning is nearly absent. QCLs are currently being grown with high accuracy and experimentally observed wavelengths agree closely with the designed values. If there is a detuning of Δω between the gain and the absorbing periods and the light pulses are tuned to the gain periods, Eq. 1d becomes




  
    
      
        
          
            ∂
            
              η
              a
            
          
          
            ∂
            t
          
        
        =
        
          
            
              
                i
                ⁢
                
                  
                
                ⁢
                
                  μ
                  a
                
              
              
                2
                ⁢
                ħ
              
            
            ⁢
            E
            ⁢
            
              
            
            ⁢
            
              Δ
              a
            
          
          -
          
            
              (
              
                
                  1
                  
                    T
                    2
                  
                
                -
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  
                    Δ
                    ω
                  
                
              
              )
            
            ⁢
            
              
                η
                a
              
              .
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




Then the solutions for ηa and Δa change. Analytical solutions for ηa and Δa may be found in the presence of detuning Δω when τT2, so that the term proportional to 1/T2 may be neglected in the polarization equation. Then, the solutions for ηa and Δa become




  
    
      
        
          
            η
            a
          
          =
          
            
              
                
                  
                    Δ
                    ω
                  
                  ⁢
                  τ
                
                
                  1
                  +
                  
                    
                      (
                      
                        
                          Δ
                          ω
                        
                        ⁢
                        τ
                      
                      )
                    
                    2
                  
                
              
              ⁢
              sech
              ⁢
              
                
              
              ⁢
              x
            
            +
            
              i
              ⁢
              
                1
                
                  1
                  +
                  
                    
                      (
                      
                        
                          Δ
                          ω
                        
                        ⁢
                        τ
                      
                      )
                    
                    2
                  
                
              
              ⁢
              sech
              ⁢
              
                
              
              ⁢
              x
              ⁢
              
                
              
              ⁢
              tanh
              ⁢
              
                
              
              ⁢
              x
            
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            Δ
            a
          
          =
          
            
              -
              1
            
            +
            
              
                2
                
                  1
                  +
                  
                    
                      (
                      
                        
                          Δ
                          ω
                        
                        ⁢
                        τ
                      
                      )
                    
                    2
                  
                
              
              ⁢
              sech
              ⁢
              
                
              
              ⁢
              x
            
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          b
        
        )
      
    
  




where x=t/τ−z/ντ and Δa0=−1 at t→∞.

On physical grounds, it is apparent that the criterion for acceptable detuning is that Δω≦1/T2, since τ≦T2 and the bandwidth of the pulse in angular frequency is  ̃τ−1. If T2 is 100 fs, and we demand conservatively that Δω<1/T2, then Δω<1012 s−1, corresponding to an allowed frequency detuning of 1.6×1011 Hz, which is 2% of the carrier frequency of 8 μm light and is not overly demanding.
In order for the solution reported in Eq. (5) to be of any practical interest, it must be robust when μa differs from 2μg, when T1g and T1a are on the order of a picosecond or less, when an initial pulse that is long compared to its final, stable duration is injected into the medium, and when the initial pulse area differs from the ideal value of π in the gain medium and 2π in the absorbing medium. The Maxwell-Bloch equations must be solved computationally to determine what happens under these conditions. For computational analysis, we normalize Eq. (1). We define Ē=(μgT2g/ h)E and we introduce the retarded time t′=t−(n/c)z, the normalized time t=t/T2g, and the normalized distance z=lz, so that Eq. (1) becomes




  
    
      
        
          
            
              ∂
              
                E
                _
              
            
            
              ∂
              
                z
                _
              
            
          
          =
          
            
              
                -
                i
              
              ⁢
              
                g
                _
              
              ⁢
              
                
              
              ⁢
              
                η
                g
              
            
            -
            
              i
              ⁢
              
                
                  a
                  _
                
                
                  
                    (
                    
                      
                        T
                        
                          2
                          ⁢
                          a
                        
                      
                      /
                      
                        T
                        
                          2
                          ⁢
                          g
                        
                      
                    
                    )
                  
                  ⁢
                  
                    μ
                    _
                  
                
              
              ⁢
              
                η
                a
              
            
            -
            
              
                1
                2
              
              ⁢
              
                E
                _
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              ∂
              
                η
                g
              
            
            
              ∂
              
                t
                _
              
            
          
          =
          
            
              
                i
                2
              
              ⁢
              
                Δ
                g
              
              ⁢
              
                E
                _
              
            
            -
            
              η
              g
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              ∂
              
                Δ
                g
              
            
            
              ∂
              
                t
                _
              
            
          
          =
          
            
              i
              ⁡
              
                (
                
                  
                    
                      η
                      g
                    
                    ⁢
                    
                      
                        E
                        _
                      
                      *
                    
                  
                  -
                  
                    
                      η
                      g
                      *
                    
                    ⁢
                    
                      E
                      _
                    
                  
                
                )
              
            
            +
            
              
                
                  Δ
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
                -
                
                  Δ
                  g
                
              
              
                
                  T
                  
                    1
                    ⁢
                    
                      
                    
                    ⁢
                    g
                  
                
                /
                
                  T
                  
                    2
                    ⁢
                    g
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            
              ∂
              
                η
                a
              
            
            
              ∂
              
                t
                _
              
            
          
          =
          
            
              
                i
                2
              
              ⁢
              
                μ
                _
              
              ⁢
              
                
              
              ⁢
              
                Δ
                a
              
              ⁢
              
                E
                _
              
            
            -
            
              
                η
                a
              
              
                
                  T
                  
                    2
                    ⁢
                    a
                  
                
                /
                
                  T
                  
                    2
                    ⁢
                    g
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            
              ∂
              
                Δ
                a
              
            
            
              ∂
              
                t
                _
              
            
          
          =
          
            
              i
              ⁢
              
                
              
              ⁢
              
                
                  μ
                  _
                
                ⁡
                
                  (
                  
                    
                      
                        η
                        a
                      
                      ⁢
                      
                        
                          E
                          _
                        
                        *
                      
                    
                    -
                    
                      
                        η
                        a
                        *
                      
                      ⁢
                      
                        E
                        _
                      
                    
                  
                  )
                
              
            
            +
            
              
                
                  Δ
                  
                    a
                    ⁢
                    
                      
                    
                    ⁢
                    0
                  
                
                -
                
                  Δ
                  a
                
              
              
                
                  T
                  
                    1
                    ⁢
                    
                      
                    
                    ⁢
                    a
                  
                
                /
                
                  T
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    g
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          e
        
        )
      
    
  




where g=g/l, ā=a/l, and μ=μa/μg.

In our simulations, we used different window sizes, depending on the pulse evolution, and we verified that the pulse intensities are always zero at the window boundaries. We spaced our node points 1-5 fs apart, and we chose a step size between 1 and 10 μm, depending on the material parameters in the simulation. In each simulation these values were constant and we checked that reducing these values made no visible difference in our plotted results. Finally, we verified by extending the propagation length that we were following the pulses over a sufficiently long length to reliably determine whether the pulses are stable or not. In FIG. 6, we show the limits of g and ā for stable operation with different values of T1/T2 when T1g=T1a≡T1 and T2g=T2a≡T. We begin by assuming that a hyperbolic-secant-shaped pulse having an area of π is injected into the QCL. Before the pulse is injected, the population is completely in the upper lasing state in the gain medium, i.e., Δg0=1.0 and is completely in the ground state in the absorbing medium, i.e., Δa0=−1.0. In FIG. 6, the black solid line at the bottom defines the lower limits of ā for any T1/T2. The solid curves on the top are the loss-limited boundaries for different values of T1/T2. The injected pulses are only stable when the gain and absorption parameters are set between these two boundary limits. Stable pulses propagate in the laser cavity with no change in shape and energy. We show the pulse evolution in the stable regime and unstable regimes in FIG. 7. FIG. 7a shows stable pulse evolution when g=4.0 and ā=3.5. The laser becomes unstable when operated with ā smaller than the lower limits given in FIG. 6 due to the growth of continuous waves. In this case, the net gain of the laser becomes positive, i.e., g−ā−1>0, and multiple pulses may form in the cavity, leading to the generation of multiple pulses in our simulations. We give an example of this behavior in FIG. 7b. In this case, we set g=4.0 and ā=1.0; the laser becomes unstable when z=20 and the laser cavity develops more than one pulse. With a greater than the upper limit, the gain medium cannot compensate for absorption and the linear loss. As a result, pulses damp. In FIG. 7c, which exhibits this behavior, we have set g=4.0 and ā=7.8. The upper limit for ā decreases when T1/T2 decreases as shown in FIG. 6, because the damping increases. We also show contours of the stable normalized pulse duration, τ=τFWHM/1.763T2, with dashed lines in FIG. 6. Pulse durations are on the order of T2 when g≧2.5 and ā≧2.0. The pulse durations can be made arbitrarily short by increasing g and ā. However, g and ā are proportional to the current, so that the current must be increased. If T1/T2 decreases, then g must increase if ā is constant in order to maintain ti at a constant value.
As we discussed, in Sec. III, in a practical QCL design, we must have T1g>T1a. For generality, we consider here the stability limits as T1a/T1g varies between 0.5 and 2.0. FIG. 8 shows the stability limits of g and ā as T1a/T1g is varied. The solid black line at the bottom is the lower limit of ā and remains the same for any T1a/T1g. However, the upper limit of ā decreases when T1a/T1g decreases. The analytical solution of the Maxwell-Bloch equations given in Eq. (5) assumes that the absorbing medium has a dipole moment twice that of the gain medium, i.e., μa=2μg. The condition μa=2μg will not be exactly satisfied due to design constraints and growth limitations. The QCL gain is determined by μg and T1g. To produce large gain, it is preferable that μg is large. In a vertical transition QCL, the dipole moment is generally 2 nm. In diagonal-transition QCLs, the dipole moment is ≧1.4 nm, which is smaller. Despite the smaller value of μg with diagonal transitions, we must have μa/e≧2.8 to satisfy the condition μa=2μg. Therefore, it is useful in practice if SIT mode locking is possible when μa<2μg. We determine the stability limits of μ=μa/μg for stable operation. FIG. 9 shows the lower and upper stability limits of μ vs g as ā varies. The solid lines in FIG. 9 indicate the lower limits for μ while the dashed lines indicate the upper limits. The two ends of each of the lines are at the stability boundaries for g at each particular ā. In each of the cases, the minimum value of μ is approximately 2 when g is near its minimum, below which an input pulse attenuates. As g increases toward the limit at which continuous waves become unstable, the minimum value of μ required for stable operation decreases significantly. Pulses are stable with μ ̃1.2 when g=3.5, 4.0, and 4.5 with ā=2.5, 3.0, and 3.5, respectively, with g just below the stability limit for generating continuous waves. However, the stable pulse duration increases significantly as μ decreases. When μ is below the solid lines in FIG. 9, pulses attenuate. The minimum value of μ required for stable operation increases as ā increases for any fixed g. We have found no hard upper limit to stability as μ increases, although the pulses are increasingly distorted. The dashed lines in FIG. 9 indicate the values of μ at which the pulses become double peaked. We have derived an energy balance equation in Eq. (16) that defines the limits of the input energy for stable operation. However, Eq. (16) assumes that input pulse has a duration, τi=τFWHM/1.763T2T1, so that the effects of a finite coherence time T2 and damping due to finite T1 may be ignored. If this condition is not satisfied, then Eq. (16) is no longer valid. From a practical standpoint, an input pulse having a duration on the order of T2 or longer than T2 is advantageous. We have calculated the dependence of the minimum and maximum input energy on the input pulse duration for two different combinations of gain and loss. We show the results in FIG. 10. The input pulse duration is normalized to T2g=T2a≡T2 and is plotted on a logarithmic scale. The value of T1/T2 has been set to 10. When τi/T2=0.1, we find that the minimum normalized energy θ=(μg/ h)∫−∞∞Edt=∫−∞∞Ēd t that is required for stable operation is 0.30π when g=3.5 and ā=3.0. However, as we increase τi/T2, the minimum normalized pulse energy that is required for stable operation increases significantly due to the pulse's decorrelation over its duration. It increases to 0.42π when τi/T2=1, 1.31π when τi/T2=10, and 9.59π when τi/T2=100. Pulses are stable for an input energy of at least 20π when τi/T2≦3. We find that pulses split into multiple pulses when the input pulse energy is ≧2π. However, at the stable pulse duration τ ̃0.5 for the parameters g=3.5, ā=3.0, only one pulse is stable, and the others damp even with an initial normalized energy of 20π when τi/T2<4. When τi/T2<4, continuous waves become unstable. We find that multiple pulses are generated when the input energy is ≧3π with τi/T2=10. The upper stability limit for the input energy decreases as τi/T2 increases when τi/T2≦10. However, beyond that point, the upper stability limit increases with τi/T2 as damping due to T1 comes into effect. The stability limits for the input normalized energy when g=3.5 and ā=3.5 show a similar trend with the exception that both the stability limits are shifted upward due to an increase in absorption. We simulated a number of cases in which we investigated the effect of detuning the absorbing medium from the gain medium and the carrier frequency of the light. Setting T1/T2=10, Δg0=1.0, and Δa0=−1.0, we found that stable operation can be obtained with a detuning ΔωT2≦0.53 when g=3.5, ā=3.5. Stable operation can be obtained with ΔωT2≦0.36 when g=3.5, ā=3.0, and with ΔωT2≦0.15 when g=3.5, ā=2.5.
Accordingly, in the context of the above teaching, it has been shown that by combining absorbing with gain periods in a QCL, one can create nearly ideal conditions to observe SIT mode locking and thereby obtain pulses that are less than 100 fs long from a mid-infrared laser. Here we also show that this analysis can be extended using detailed computational studies of the Maxwell-Bloch equations, to illustrate the stability of the solutions as the equation parameters vary. The solutions demonstrate the robustness of the SIT mode-locking technique and that QCLs can be mode-locked using the SIT effect within practically achievable parameter regimes.
Any references recited or referred to herein are incorporated herein in their entirety, particularly as they relate to teaching the level of ordinary skill in this art and for any disclosure necessary for the commoner understanding of the subject matter of the claimed invention. It will be clear to a person of ordinary skill in the art that the above embodiments may be altered or that insubstantial changes may be made without departing from the scope of the invention. Accordingly, the scope of the invention is determined by the scope of the following claims and their equitable equivalents.

## Claims

1. A self-induced transparency mode-locked quantum cascade laser, comprising:
(i) an active section comprising a plurality of quantum well layers deposited in alternating layers on a plurality of quantum barrier layers and form a sequence of alternating gain and absorbing periods, said alternating gain and absorbing periods interleaved along the growth axis of the active section, wherein the absorbing periods have a dipole moment of about twice that of the gain periods;
(ii) an optical cavity that houses the active section and permits amplified light to escape;
(iii) an externally supplied seed pulse; and
(iv) current injectors structured and arranged to apply an electric control field to the active section.

2. The laser of claim 1, further comprising wherein the quantum barrier layers are made using Indium-Aluminum-Arsenide, and the quantum well layers are made using Indium-Gallium-Arsenide.

3. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser.

4. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser having a wavelength of between about 8 micrometers and about 12 micrometers.

5. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide coherence times  gain recovery times  round trip times.

6. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide coherence times of at least 100 femtoseconds, gain recovery times of at least 1 picosecond, and round trip times of at least 50 picoseconds.

7. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide a heterostructure having a ratio of 4 gain periods:1 absorbing period.

8. The laser of claim 1, further comprising wherein the input pulse is a it pulse in the gain medium.

9. The laser of claim 1, further comprising wherein each gain period and each absorbing period have over 16 quantum layers.

10. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to suppress continuous waves and eliminate spatial hole burning.

11. The laser of claim 1, further comprising wherein the gain and absorbing periods are designed to provide a mid-IR wavelength laser having pulse length less than 100 femtoseconds long.