Researchers have formulated an advanced theoretical framework coupling microscopic density matrix transport calculations with effective semiconductor Maxwell–Bloch equations to simulate the nonlinear dynamics of a fabricated terahertz quantum cascade laser under external optical feedback. By extracting drive-current-dependent physical parameters directly from the quantum heterostructure layout, the study establishes an end-to-end modeling pipeline connecting quantum-level carrier scattering to macroscopic laser behavior at threshold and maximum operating currents. Numerical evaluations demonstrate that external optical reinjection reorganizes the laser emission, inducing transitions from a single-mode state to structured multimode regimes while simultaneously lowering the threshold current required for lasing. The investigated bound-to-continuum gallium arsenide device exhibits distinct dynamical states dictated by the distance to an external reflector, producing fundamental frequency combs, second harmonic combs, or clustered spectra. This predictive approach provides quantitative insights into self-mixing and mode competition, enabling the optimization of semiconductor laser designs for reliable high-frequency sensing and optical communication.
Quantum cascade lasers emitting in the terahertz spectrum operate within the traditionally challenging far-infrared spectral range, generating coherent radiation between 1.2 and 6.03 terahertz in pulsed mode and between 1.2 and 5.71 terahertz in continuous-wave operation. These versatile semiconductor sources support applications in medical diagnostics, chemical detection, and non-destructive imaging, achieving continuous-wave operation up to 129 kelvin and pulsed functioning up to 261 kelvin. Demonstrating experimental performance limits with peak pulsed powers reaching 2.4 watts and continuous-wave powers reaching 312 milliwatts, these lasers frequently experience external optical feedback during practical sensing. When emitted light reflects off an external target and re-enters the laser cavity, delayed reinjection alters the internal electric field and carrier densities through self-mixing. Depending on external reflection conditions, this self-mixing interaction can dramatically perturb device stability or yield precise interferometric information, making the accurate simulation of optical feedback dynamics crucial for next-generation terahertz instrumentation.
Analyzing semiconductor laser dynamics has relied extensively on the coupled Maxwell–Bloch equations since the 1980s, tracing atomic state populations, electric fields, and macroscopic medium polarizations through density matrix formalism. Although indispensable for modeling nonlinear light–matter interactions, applying Maxwell–Bloch equations to quantum cascade devices involves notable numerical complexity and analytical constraints. Conventional formulations incorporate relaxation and dephasing through fixed phenomenological time constants that fail to reflect complex microscopic scattering and carrier transport across varied operational biases. Alternative frameworks, such as cavity mean-field theory, approximate the intracavity optical field by averaging across round trips, effectively describing global phase locking and frequency combs but omitting longitudinal spatial hole burning and delayed reflection boundaries. Master-equation techniques capture dispersion and spatial hole burning yet retain simplifying approximations. By contrast, effective semiconductor Maxwell–Bloch equations maintain full spatial and traveling-wave field resolution, offering a rigorous platform for describing local microscopic interactions under external optical feedback.
To enhance physical realism, the researchers modified the effective semiconductor Maxwell–Bloch framework by incorporating bias-dependent input parameters computed via microscopic density matrix transport simulations. Rather than adopting static empirical estimates, the model directly evaluates material gain, electron volume concentration, non-radiative carrier lifetime, and polarization decoherence time as functions of operating current. The transport calculations account for electron interactions with longitudinal optical phonons, acoustic phonons, alloy disorder, interface roughness, ionized impurities, and electron–electron scattering evaluated through Fermi’s golden rule. These microscopic properties inform a phenomenological optical susceptibility that accounts for asymmetric, frequency-dependent gain, refractive index dynamics, and carrier density, while multiple-scale expansion describes spatial hole burning and standing wave formation. Differential gain is calculated directly from peak material gain, electron concentration, and polarization lifetime. Incorporating reflected light through time-delayed facet boundary conditions establishes the first comprehensive simulation chain predicting dynamic behavior under external optical feedback directly from the heterostructure specifications of a fabricated laser.
The integrated simulation approach was implemented for a specific bound-to-continuum gallium arsenide and aluminum gallium arsenide terahertz quantum cascade laser engineered for 2.0-terahertz emission. The semiconductor active region features an alternating layer structure starting with an injection barrier: 5.0, 14.4, 1.0, 11.8, 1.0, 14.4, 2.4, 14.4, 2.4, 13.2, 3.0, 12.4, 3.2, 12.0, 4.4, and 12.6 nanometers. Barriers comprise aluminum gallium arsenide with ten percent aluminum fraction, while active wells are doped to 1.3 times ten to the sixteenth electrons per cubic centimeter. In this architecture, electron extraction from the lower laser level to the injection level occurs through diagonal transitions within an energy miniband, with the injection state strongly coupled to the upper laser level of the adjacent period. Density matrix calculations established a free-running lasing threshold current of 430 milliamperes, alongside a negative differential resistance current of 580 milliamperes that marks the upper operational boundary of the device.
Operating under free-running conditions at the 430-milliampere threshold current without external reflections, numerical discretization via the time-domain traveling-wave finite difference method reveals pure single-mode continuous-wave operation. Due to rapid depopulation of the lower laser level, the carrier density tracks the upper laser level population, and the active medium oscillates collectively without forming traveling-wave modulations. The emitted optical spectrum exhibits an isolated single mode at 2.2 terahertz, coinciding with the cold cavity mode nearest the gain peak. The simulated output power waveform reaches a stable steady state in less than 10 nanoseconds, producing an average power output of 0.28 milliwatts. On sub-picosecond scales, the instantaneous output power oscillates at 4.4 terahertz, twice the electric field frequency, exhibiting consecutive intensity peaks separated by 0.227 picoseconds. This initial baseline confirms that the 2.2-millimeter Fabry–Pérot cavity behaves as a pristine, low-intensity single-frequency source when operating just above lasing threshold in the absence of feedback.
Introducing optical feedback from an external reflector with a reflectivity of 0.3 dramatically alters near-threshold dynamics, driving adjacent Fabry–Pérot modes into lasing and initiating multimode emission. For an internal cavity length of 2.2 millimeters and refractive index of 3.6, the internal free spectral range is roughly 18.9 gigahertz, defining a round-trip time of 52.8 picoseconds. Adjusting the external cavity length to 8.0 millimeters satisfies the fundamental comb condition, transitioning the emission into a fundamental frequency comb spaced by 18.9 gigahertz. Stabilizing after roughly 50 nanoseconds, the laser emits an increased average power of 0.68 milliwatts, displaying repetitive intensity peaks. Alternatively, configuring a 4.0-millimeter external cavity satisfies the second-harmonic comb condition, amplifying alternate modes while suppressing intermediate ones. This generates a second harmonic comb with 37.8-gigahertz spacing, reaching steady state in approximately 10 nanoseconds with an average power of 1.31 milliwatts.
A crucial physical effect demonstrated and quantified by the model is the feedback-induced reduction of the laser threshold current. Coherent reinjection effectively supplements the circulating intracavity field, lowering net resonator losses and reducing the optical gain required to initiate lasing. To determine the new operating threshold under feedback, the researchers simulated the 4.0-millimeter external cavity while progressively lowering the injection current and updating transport parameters accordingly. The simulations revealed that external optical feedback reduced the effective threshold current from 430 milliamperes down to 417 milliamperes. At this reduced operating current of 417 milliamperes, the device returns to a single-mode regime near threshold, with the optical power settling after 100 nanoseconds to an average value of 0.263 milliwatts. While conventional Maxwell–Bloch simulations with static parameters also demonstrate qualitative threshold reductions, incorporating density matrix transport physics establishes accurate quantitative values for device-specific performance under optical feedback.
Extending the external reflector distance to 20 centimeters introduces a substantially longer propagation delay, creating a mixed multimode regime characterized by clustered mode dynamics. Under these conditions, the external cavity imposes a distinct resonance structure with an external free spectral range of approximately 0.75 gigahertz. At the 430-milliampere threshold current, the output power waveform requires roughly 300 nanoseconds to achieve steady state due to prolonged feedback round-trip times, delivering an average emitted power of 1.27 milliwatts. Rather than sustaining isolated lines, each longitudinal mode of the internal resonator, spaced at the 18.9-gigahertz internal free spectral range, splits into a discrete cluster of closely spaced sub-modes separated by 0.75 gigahertz. This spectral clustering reflects direct competition between the internal cavity and the external resonator, generating disconnected mode groups that prevent the formation of a continuous optical frequency comb while displaying complex interference patterns across multiple internal cavity round trips.
Operating at the maximum operating current of 580 milliamperes, corresponding to the negative differential resistance point, the laser shifts into an inherently multimode state with 14 modes within a 15-decibel interval centered at 2.15 terahertz, producing 1.8 milliwatts of average power. Introducing optical feedback alters this multimode emission: a 4.0-millimeter cavity enforces a second harmonic comb of 8 modes spaced by 37.8 gigahertz with 1.3 milliwatts average power, while a 20-centimeter cavity creates mode clusters spaced by 18.9 gigahertz with internal 0.75-gigahertz sub-mode spacing, emitting 1.04 milliwatts. The research team evaluated feedback-based displacement sensing by varying external cavity length in 0.5-micrometer steps between 0.1 and 0.1002 meters. Near threshold at 430 milliamperes, interference produces sharp, stable periodic fringes with a 68-micrometer period corresponding to half the 2.2-terahertz wavelength, providing high measurement sensitivity. Conversely, at 580 milliamperes, chaotic multimode interference degrades signal clarity, demonstrating that near-threshold, low-mode-count operation is essential for reliable terahertz sensing.
Subject of Research: Terahertz quantum cascade laser dynamics under external optical feedback
Article Title: Analysis of self-mixing in THz quantum cascade lasers under varying mode dynamics and external optical feedback
Article References: Stanojević, N., Vuković, N., Ignjatović, M., Basta, N., & Radovanović, J. (2026). Analysis of self-mixing in THz quantum cascade lasers under varying mode dynamics and external optical feedback. Applied Physics B, 132(6), Article 76. https://doi.org/10.1007/s00340-026-08682-5
Image Credits: AI Generated
DOI: 10.1007/s00340-026-08682-5
Keywords: Quantum cascade lasers, Terahertz radiation, Optical feedback, Self-mixing interferometry, Maxwell-Bloch equations, Frequency combs, Laser dynamics, Semiconductor physics, Analysis, self-mixing, quantum, cascade
News Source: Katie Riggs. (November 12, 2026). Dynamic Modeling Uncovers Feedback Effects in Terahertz Quantum Cascade Lasers. Scienmag.
