Skip to content

Stimulated emission is a downward matter transition induced by radiation already occupying the receiving field mode. For an upper state ∣e⟩|e\rangle, a lower state ∣g⟩|g\rangle, and one photon mode μ\mu,

∣e;nμ⟩⟶∣g;nμ+1⟩.|e;n_\mu\rangle \longrightarrow |g;n_\mu+1\rangle.

The matter system loses energy ℏω0\hbar\omega_0, while the selected mode gains one photon. The transition amplitude contains

nμ+1.\sqrt{n_\mu+1}.

Its squared magnitude therefore contains nμ+1n_\mu+1: the term proportional to nμn_\mu is stimulated emission, and the unity is the vacuum contribution associated with spontaneous emission.

The same interaction also causes absorption. Net amplification therefore depends on level populations, line strengths, and propagation losses. A laser additionally needs feedback and mode selection.

Required background. Two-Level Atom supplies the state reduction; Spontaneous Emission explains vacuum decay; Transition Rates in Light–Matter Interaction derives the mode normalization and the nn versus n+1n+1 factors.

Use two matter states with

Ee−Eg=ℏω0>0.E_e-E_g = \hbar\omega_0 > 0.

Throughout this page:

  • nμn_\mu is the photon occupation of one specified field mode μ\mu;
  • N1\mathcal N_1 and N2\mathcal N_2 are number densities in lower and upper matter levels;
  • g1g_1 and g2g_2 are the corresponding level degeneracies;
  • deg=⟨e∣d∣g⟩\mathbf d_{eg}=\langle e|\mathbf d|g\rangle is the transition dipole;
  • Δ=ω0−ω\Delta=\omega_0-\omega is the atom-minus-field detuning;
  • II is cycle-averaged intensity and Φ=I/(ℏω)\Phi=I/(\hbar\omega) is photon flux;
  • σa\sigma_a and σe\sigma_e are absorption and stimulated-emission cross sections for a declared polarization and propagation mode;
  • g(ω)g(\omega) is a power-gain coefficient with dimensions of inverse length, not a dimensionless gain;
  • G=Iout/IinG=I_{\mathrm{out}}/I_{\mathrm{in}} is a dimensionless power gain;
  • αbg\alpha_{\mathrm{bg}} is a distributed background-loss coefficient.

The symbol gg is also widely used for an atom–mode coupling constant and for a level degeneracy. Here the coupling is written gμg_\mu, degeneracies carry subscripts g1,g2g_1,g_2, and propagation gain is written g(ω)g(\omega).

For one near-resonant field mode, the electric-dipole interaction can be reduced to

HI=ℏgμσ+aμ+ℏgμ∗σ−aμ†,H_I = \hbar g_\mu \sigma_+a_\mu + \hbar g_\mu^* \sigma_-a_\mu^\dagger,

where

σ+=∣e⟩⟨g∣,σ−=∣g⟩⟨e∣.\sigma_+ = |e\rangle\langle g|, \qquad \sigma_- = |g\rangle\langle e|.

The term σ+aμ\sigma_+a_\mu absorbs a photon while exciting the matter system. The Hermitian-conjugate term σ−aμ†\sigma_-a_\mu^\dagger lowers the matter system while creating a photon.

For absorption,

⟨e;nμ−1∣HI∣g;nμ⟩=ℏgμnμ.\langle e;n_\mu-1| H_I |g;n_\mu\rangle = \hbar g_\mu \sqrt{n_\mu}.

For downward emission,

⟨g;nμ+1∣HI∣e;nμ⟩=ℏgμ∗nμ+1.\langle g;n_\mu+1| H_I |e;n_\mu\rangle = \hbar g_\mu^* \sqrt{n_\mu+1}.

The corresponding squared matrix elements are proportional to

nμandnμ+1.n_\mu \quad\text{and}\quad n_\mu+1.

No additional interaction among photons is required. The enhancement follows from the normalization of bosonic occupation states.

The creation operator obeys

aμ†∣nμ⟩=nμ+1∣nμ+1⟩.a_\mu^\dagger|n_\mu\rangle = \sqrt{n_\mu+1} |n_\mu+1\rangle.

Equivalently,

∥aμ†∣nμ⟩∥2=⟨nμ∣aμaμ†∣nμ⟩=⟨nμ∣(aμ†aμ+1)∣nμ⟩=nμ+1.\begin{aligned} \left\| a_\mu^\dagger|n_\mu\rangle \right\|^2 &= \langle n_\mu| a_\mu a_\mu^\dagger |n_\mu\rangle \\ &= \langle n_\mu| \left( a_\mu^\dagger a_\mu+1 \right) |n_\mu\rangle \\ &= n_\mu+1. \end{aligned}

The commutator

[aμ,aμ†]=1[a_\mu,a_\mu^\dagger] = 1

provides the vacuum unity. The number operator aμ†aμa_\mu^\dagger a_\mu provides the occupation-dependent part.

This is final-state Bose enhancement: a transition into a bosonic mode is more likely when that mode is already occupied. It is a consequence of indistinguishability and symmetrization, not a force exerted by one photon on another.

The label μ\mu may include:

  • spatial mode or wavevector;
  • polarization;
  • frequency or wavepacket spectrum;
  • transverse profile;
  • propagation direction;
  • cavity or waveguide boundary conditions.

Saying that stimulated emission enters “the same mode” means that the newly created excitation occupies the same complete mode label used by aμ†a_\mu^\dagger. Frequency alone is not enough.

For a traveling plane-wave mode, this selection also fixes the added photon’s momentum and polarization. Momentum conservation then appears in the recoil of the emitter or of the larger medium that contains it.

The algebraic split

nμ+1=nμ⏟stimulated+1⏟vacuumn_\mu+1 = \underbrace{n_\mu}_{\text{stimulated}} + \underbrace{1}_{\text{vacuum}}

is exact for the mode-resolved squared matrix element. Its physical interpretation is:

  • Absorption: begins with nμn_\mu photons, carries the factor nμn_\mu, and vanishes when that mode is empty.
  • Stimulated emission: is the part of downward emission proportional to the initial occupation nμn_\mu and therefore requires incident radiation in the receiving mode.
  • Vacuum contribution: supplies the remaining factor 11 and survives when nμ=0n_\mu=0.

For fixed coupling and a common final-state linewidth, one may write schematically

Γμ,↓=γμ(nμ+1),\Gamma_{\mu,\downarrow} = \gamma_\mu \left( n_\mu+1 \right),

so that

Γμ,stim=γμnμ,Γμ,sp=γμ.\Gamma_{\mu,\mathrm{stim}} = \gamma_\mu n_\mu, \qquad \Gamma_{\mu,\mathrm{sp}} = \gamma_\mu.

Then

Γμ,stimΓμ,sp=nμ.\frac{ \Gamma_{\mu,\mathrm{stim}} }{ \Gamma_{\mu,\mathrm{sp}} } = n_\mu.

This ratio is mode resolved. It does not justify comparing the occupation of one driven mode with the total spontaneous rate into every vacuum mode. Free-space spontaneous emission sums a small contribution over a continuum of directions, polarizations, and frequencies. A narrow laser beam may put a large occupation into only one family of modes.

A discrete lossless mode is not automatically a rate

Section titled “A discrete lossless mode is not automatically a rate”

The matrix elements above are exact, but a constant irreversible rate requires additional conditions: a continuum, finite linewidth, dephasing, loss, coarse graining, or another mechanism that destroys coherent recurrences.

One isolated emitter resonantly coupled to one ideal mode can exchange an excitation reversibly. In that regime, the n+1\sqrt{n+1} factor sets a Rabi frequency rather than a golden-rule decay rate. The irreversible rate language is recovered only in a suitable weak-coupling or coarse-grained limit.

Stimulated emission appears in several formalisms. They answer different questions and should not be mixed without checking their assumptions.

The Fock-state transition

∣e;n⟩⟶∣g;n+1⟩|e;n\rangle \longrightarrow |g;n+1\rangle

is best for:

  • identifying the n+1n+1 operator factor;
  • separating vacuum and occupation-dependent contributions;
  • counting photons in a specified mode;
  • analyzing few-photon or cavity processes.

A number state has no definite optical phase. Therefore the statement “stimulated emission has the phase of the incident photon” is not meaningful for a pure Fock state.

A coherent state satisfies

a∣α⟩=α∣α⟩,n‾=∣α∣2.a|\alpha\rangle = \alpha|\alpha\rangle, \qquad \overline n = |\alpha|^2.

For large n‾\overline n, replacing the field operator by its coherent amplitude leads to a prescribed classical drive. The induced matter polarization can radiate phase coherently into that drive, increasing the field amplitude in a phase-matched mode.

For a weak, short interaction followed by detection of the atom in its lower state, the leading-order conditional field operation is

a†∣α⟩,a^\dagger|\alpha\rangle,

which is a photon-added coherent state and is not exactly proportional to ∣α⟩|\alpha\rangle. It is not the exact operation for an arbitrary interaction time: resonant Jaynes–Cummings evolution weights each input coefficient cnc_n by sin⁡(∣gμ∣tn+1)\sin(|g_\mu|t\sqrt{n+1}) before producing ∣n+1⟩|n+1\rangle. The creation-operator rule follows when ∣gμ∣tn+1≪1|g_\mu|t\sqrt{n+1}\ll1 over the occupied distribution. A classical drive further assumes negligible depletion and atom–field correlations on the timescale of interest.

If field phase memory and optical coherence can be neglected, populations may obey Einstein or cross-section rate equations. This is appropriate for:

  • broadband or incoherent radiation;
  • weak probes in media with rapid dephasing;
  • population kinetics averaged over many events;
  • small-signal amplifier calculations.

It is not appropriate for a short coherent pulse that produces resolved Rabi oscillations. Coherent evolution may later reduce to a rate equation after adiabatic elimination of the optical coherence, but that reduction must be stated.

One process, different retained information

Section titled “One process, different retained information”
Description Retains Common output
Quantized mode photon number and atom–field amplitudes n+1n+1, few-photon dynamics
Semiclassical coherent field optical phase and matter coherence Rabi frequency, susceptibility
Einstein kinetics populations and spectral energy density upward and downward rates
Propagation model intensity, cross sections, spatial loss gain per unit length

These descriptions agree in their common weak-coupling limits. They are not interchangeable outside those limits.

For lower and upper level populations N1\mathcal N_1 and N2\mathcal N_2 in an isotropic unpolarized radiation field, with equal populations within each degenerate level, the sharp-line Einstein model gives

Rabs=N1B12(ω)uω(ω0),\mathcal R_{\mathrm{abs}} = \mathcal N_1 B_{12}^{(\omega)} u_\omega(\omega_0), Rstim=N2B21(ω)uω(ω0),\mathcal R_{\mathrm{stim}} = \mathcal N_2 B_{21}^{(\omega)} u_\omega(\omega_0),

and

Rsp=N2A21.\mathcal R_{\mathrm{sp}} = \mathcal N_2 A_{21}.

Here uωu_\omega is total isotropic radiation energy density per unit angular frequency. The superscript is essential because changing the spectral variable changes the numerical BB coefficient.

The net stimulated production of radiation quanta is

Rnet=Rstim−Rabs=[N2B21(ω)−N1B12(ω)]uω.\begin{aligned} \mathcal R_{\mathrm{net}} &= \mathcal R_{\mathrm{stim}} - \mathcal R_{\mathrm{abs}} \\ &= \left[ \mathcal N_2 B_{21}^{(\omega)} - \mathcal N_1 B_{12}^{(\omega)} \right] u_\omega. \end{aligned}

Thermal detailed balance gives

g1B12(ω)=g2B21(ω).g_1B_{12}^{(\omega)} = g_2B_{21}^{(\omega)}.

Therefore,

Rnet=B21(ω)uω×(N2−g2g1N1).\mathcal R_{\mathrm{net}} = B_{21}^{(\omega)} u_\omega \times \left( \mathcal N_2 - \frac{g_2}{g_1} \mathcal N_1 \right).

Ordinary two-level gain requires

N2g2>N1g1.\frac{\mathcal N_2}{g_2} > \frac{\mathcal N_1}{g_1}.

This is population inversion per substate. The shortcut N2>N1\mathcal N_2>\mathcal N_1 is valid only when g1=g2g_1=g_2 and the two line profiles are reciprocal.

For a thermal free-space field and matching Einstein conventions,

RstimRsp=B21(ω)uωA21=n‾(ω0,T),\frac{ \mathcal R_{\mathrm{stim}} }{ \mathcal R_{\mathrm{sp}} } = \frac{ B_{21}^{(\omega)}u_\omega }{ A_{21} } = \overline n(\omega_0,T),

where

n‾(ω,T)=1eℏω/(kBT)−1.\overline n(\omega,T) = \frac{1}{ e^{\hbar\omega/(k_{\mathrm B}T)}-1 }.

Thermal stimulation is consequently negligible for most optical transitions at room temperature but can dominate at microwave frequencies. Absorption from the same thermal field must be included at the same time.

Einstein coefficients are natural for isotropic radiation. A directed probe beam is usually described by cross sections.

For a monochromatic traveling wave with photon flux

Φ=Iℏω,\Phi = \frac{I}{\hbar\omega},

define the weak-field rates per particle by

Wabs=σa(ω)Φ,W_{\mathrm{abs}} = \sigma_a(\omega)\Phi, Wstim=σe(ω)Φ.W_{\mathrm{stim}} = \sigma_e(\omega)\Phi.

Both σa\sigma_a and σe\sigma_e have dimensions of area. They include the declared polarization, line shape, and averaging over magnetic or molecular substates.

For resolved reciprocal substates,

σa(ω)=σe(ω)\sigma_a(\omega) = \sigma_e(\omega)

when the same field mode and normalization are used. For level-averaged cross sections with reciprocal profiles,

g1σa(ω)=g2σe(ω).g_1\sigma_a(\omega) = g_2\sigma_e(\omega).

Broad vibronic manifolds need more care: emission and absorption profiles can differ because the initial sublevel distributions differ. Relations such as the McCumber relation require thermalized manifolds and additional spectroscopic assumptions; they are not universal identities between arbitrary measured spectra.

For a nondegenerate electric-dipole transition driven by a plane wave, define

Ω=∣ϵ⋅deg∣E0ℏ,\Omega = \frac{ \left| \boldsymbol\epsilon\mathbin{\cdot}\mathbf d_{eg} \right| E_0 }{ \hbar },

with

I=12ϵ0cE02.I = \frac12\epsilon_0cE_0^2.

If the optical coherence decays at rate γ2\gamma_2, setting the weak-field coherence to its stationary value gives the stimulated rate coefficient

W(Δ)=Ω22γ2Δ2+γ22.W(\Delta) = \frac{\Omega^2}{2} \frac{\gamma_2}{ \Delta^2+\gamma_2^2 }.

Using this as a time-local population rate additionally requires populations and the field envelope to vary slowly compared with coherence relaxation. It does not describe the initial coherent transient. Dividing by photon flux yields

σ(Δ)=ω∣ϵ⋅deg∣2ϵ0cℏ×γ2Δ2+γ22.\sigma(\Delta) = \frac{ \omega \left| \boldsymbol\epsilon\mathbin{\cdot}\mathbf d_{eg} \right|^2 }{ \epsilon_0c\hbar } \times \frac{\gamma_2}{ \Delta^2+\gamma_2^2 }.

This formula is state and polarization specific. For a closed, lifetime-limited, optimally polarized transition,

γ2=Γ2,\gamma_2 = \frac{\Gamma}{2},

and the resonant value becomes

σ0=3λ022π.\sigma_0 = \frac{3\lambda_0^2}{2\pi}.

Branching, Clebsch–Gordan coefficients, polarization mismatch, pure dephasing, Doppler averaging, and inhomogeneous broadening reduce or redistribute the effective peak cross section. The compact resonant formula must not be applied to an arbitrary multilevel line without those factors.

Consider a narrow probe propagating along zz. In the weak-probe, slowly-varying limit,

dI(ω,z)dz=g0(ω)I(ω,z),\frac{dI(\omega,z)}{dz} = g_0(\omega)I(\omega,z),

where the small-signal net gain coefficient is

g0(ω)=N2σe(ω)−N1σa(ω)−αbg(ω).g_0(\omega) ={} \mathcal N_2\sigma_e(\omega) - \mathcal N_1\sigma_a(\omega) - \alpha_{\mathrm{bg}}(\omega).

Each term has dimensions of inverse length. The signs have direct meanings:

  • N2σe\mathcal N_2\sigma_e adds probe photons by stimulated emission;
  • N1σa\mathcal N_1\sigma_a removes probe photons by absorption;
  • αbg\alpha_{\mathrm{bg}} includes scattering, parasitic absorption, and other distributed losses.

If the populations and cross sections are constant over a length LL,

I(L)=I(0)eg0L.I(L) = I(0)e^{g_0L}.

The power gain is

G=I(L)I(0)=eg0L,G = \frac{I(L)}{I(0)} = e^{g_0L},

and its decibel value is

GdB=10log⁡10G.G_{\mathrm{dB}} = 10\log_{10}G.

A positive stimulated-emission term does not guarantee net amplification. The net coefficient must satisfy

g0(ω)>0.g_0(\omega)>0.

With the propagation convention

dIdz=−αI,\frac{dI}{dz} = -\alpha I,

an amplifying medium has

α(ω)<0,g(ω)=−α(ω)>0.\alpha(\omega)<0, \qquad g(\omega)=-\alpha(\omega)>0.

In linear-response language, inversion reverses the sign of the dissipative part of the resonant susceptibility. The pumped medium is active, so the equilibrium passivity constraint does not apply: the amplified field draws energy from stored excitation supplied by the pump.

The sign and size of g0g_0 depend on frequency through:

  • homogeneous natural, collisional, and power-broadened profiles;
  • Doppler or static inhomogeneous distributions;
  • upper- and lower-manifold populations;
  • polarization-dependent line strengths;
  • reabsorption and background loss;
  • the pump-induced spatial distribution of inversion.

An integrated line strength constrains area under a cross-section profile, whereas peak gain also depends on linewidth. Broadening can lower a peak without changing the relevant integrated strength.

Spontaneous emission into the observed spatial and spectral mode acts as a source term:

dIdz=g0I+jsp.\frac{dI}{dz} = g_0I+j_{\mathrm{sp}}.

For constant g0≠0g_0\ne0 and jspj_{\mathrm{sp}},

I(L)=I(0)eg0L+jspg0(eg0L−1).I(L) ={} I(0)e^{g_0L} + \frac{j_{\mathrm{sp}}}{g_0} \left( e^{g_0L}-1 \right).

The first term is amplified input. The second is amplified spontaneous emission. Its directional and spectral narrowing can resemble laser output, but it does not by itself establish resonator-defined oscillation or laser-like phase coherence.

Small-signal gain assumes that the probe does not appreciably change the populations. A strong field stimulates enough transitions to deplete the inversion, so the gain becomes intensity dependent.

Take equal degeneracies and equal cross sections,

σa=σe=σ.\sigma_a = \sigma_e = \sigma.

Let

D=N2−N1D = \mathcal N_2-\mathcal N_1

be the inversion density. Suppose pumping and relaxation restore a no-probe inversion D0D_0 on a time scale T1T_1. A stimulated transition moves one particle from one level to the other, changing DD by two. With

W=σIℏω,W = \frac{\sigma I}{\hbar\omega},

the rate equation is

dDdt=D0−DT1−2WD.\frac{dD}{dt} = \frac{D_0-D}{T_1} - 2WD.

The steady inversion is

Dss=D01+I/Isat,D_{\mathrm{ss}} = \frac{D_0}{ 1+I/I_{\mathrm{sat}} },

with

Isat=ℏω2σT1.I_{\mathrm{sat}} = \frac{ \hbar\omega }{ 2\sigma T_1 }.

If

ggross,0=σD0,g_{\mathrm{gross},0} = \sigma D_0,

then

gnet(I)=ggross,01+I/Isat−αbg.g_{\mathrm{net}}(I) = \frac{ g_{\mathrm{gross},0} }{ 1+I/I_{\mathrm{sat}} } - \alpha_{\mathrm{bg}}.

This is a model, not a universal law. The factor of two follows from a closed two-level inversion definition. A rapidly emptied lower laser level, multilevel branching, coherent saturation, pulsed extraction, spatial hole burning, or pump depletion changes the saturation scale and sometimes the functional form.

In a homogeneously broadened ensemble, all emitters within the line share the same transition frequency to the accuracy of the model. A strong narrow field can therefore reduce gain across the homogeneous profile.

In an inhomogeneously broadened ensemble, a narrow field first depletes the resonant subensemble. This can burn a spectral hole while off-resonant subensembles retain inversion. Collisions, spectral diffusion, and velocity changing processes can refill the hole.

Gain saturation is reversible population redistribution caused by the signal. It should be distinguished from:

  • irreversible photochemical or dielectric damage;
  • thermal lensing;
  • excited-state absorption;
  • multiphoton ionization;
  • pump depletion;
  • nonlinear refractive effects.

Several can occur together in a real amplifier, but they have different scalings and recovery times.

An ordinary gain medium needs:

  1. a pump that supplies free energy;
  2. a pathway that populates the upper laser level;
  3. sufficiently slow loss from that upper level;
  4. sufficiently rapid removal from, or weak occupation of, the lower laser level;
  5. a transition cross section large enough to overcome propagation losses.

The pump can be optical, electrical, chemical, collisional, or another nonequilibrium source. Stimulated emission transfers stored excitation to the signal; it does not supply the energy being amplified.

Why resonant two-level pumping does not create inversion

Section titled “Why resonant two-level pumping does not create inversion”

Let one incoherent resonant field drive both directions at rate WW, and let the upper state decay at rate Γ\Gamma. For a closed two-level system,

dN2dt=WN1−WN2−ΓN2.\frac{d\mathcal N_2}{dt} = W\mathcal N_1 - W\mathcal N_2 - \Gamma\mathcal N_2.

At steady state,

N2N1=WW+Γ<1.\frac{\mathcal N_2}{\mathcal N_1} = \frac{W}{W+\Gamma} < 1.

Increasing the resonant pump only approaches equal populations. It does not invert the same isolated transition. Three- and four-level schemes evade this ceiling by pumping through additional states and using asymmetric relaxation pathways.

Population inversion is required for gain in an ordinary reciprocal two-level medium whose response is described only by populations. It is not a theorem covering every driven multilevel quantum system.

Coherences can make distinct absorption amplitudes interfere destructively while leaving an emissive pathway. Such gain without inversion has been demonstrated in coherently driven multilevel atoms. This is a specialized interference effect, not evidence that the two-level balance condition was wrong.

Stimulated emission is the microscopic amplification process in a conventional laser, but it is only one part of laser operation.

Amplifier, amplified spontaneous emission, and laser

Section titled “Amplifier, amplified spontaneous emission, and laser”
  • Traveling-wave amplifier: a pumped gain medium gives one-pass gain to an externally supplied signal without requiring feedback.
  • Amplified spontaneous emission: spontaneous photons seed the pumped medium and are amplified without a required resonator or external input.
  • Laser oscillator: feedback selects a self-consistent mode whose round-trip small-signal gain reaches threshold; no external seed is required.

A pumped single-pass medium can show stimulated emission and large gain without being a laser. Conversely, a resonator without net gain only stores and filters light; it does not sustain oscillation.

Consider a linear cavity uniformly filled with gain over length LL, mirror power reflectivities R1R_1 and R2R_2, gross material gain coefficient gmg_{\mathrm m}, and distributed internal loss αi\alpha_{\mathrm i}. One round trip multiplies the intracavity intensity by

Grt=R1R2exp⁡[2(gm−αi)L].\mathcal G_{\mathrm{rt}} = R_1R_2 \exp \left[ 2 \left( g_{\mathrm m}-\alpha_{\mathrm i} \right) L \right].

The small-signal threshold is

Grt=1,\mathcal G_{\mathrm{rt}} = 1,

or

gm,th=αi+12Lln⁡(1R1R2).g_{\mathrm{m,th}} ={} \alpha_{\mathrm i} + \frac{1}{2L} \ln \left( \frac{1}{R_1R_2} \right).

Below threshold, a fluctuation decays after repeated round trips. Above threshold, it initially grows. In the classical large-photon-number limit, saturation and pump depletion reduce the gain until it balances loss. With a finite spontaneous source, stationary total photon balance includes that source: the stimulated gain can remain slightly below the loss. High-spontaneous-emission-fraction lasers need this qualification when defining a threshold.

This formula assumes uniform traveling-wave power gain on each pass and neglects standing-wave structure, spatial hole burning, mode competition, dispersion, and time dependence. Ring cavities and pulsed lasers require appropriately modified round-trip maps.

Stimulated emission preferentially feeds whatever field modes are occupied, but occupation alone does not choose a unique mode. Selection also depends on:

  • cavity resonances and diffraction loss;
  • gain bandwidth and dispersion;
  • spatial overlap with the pumped region;
  • polarization-dependent gain and loss;
  • nonlinear saturation and mode competition;
  • spontaneous-emission noise and technical perturbations.

The coherent oscillator phase is a collective dynamical property of field, medium, feedback, and noise. It is not attached to an individual photon before the process.

The phrase “the emitted photon is identical to the stimulating photon” is a useful mode-label shorthand, but it can be misleading.

No individual photon is selected for copying

Section titled “No individual photon is selected for copying”

Photons in one mode are indistinguishable. The field changes from an nn-excitation state to an (n+1)(n+1)-excitation state; quantum mechanics does not label one photon as the original and another as its copy.

For a coherent input,

a†∣α⟩a^\dagger|\alpha\rangle

is not the original coherent state with a separately cloned photon appended. Likewise, stimulated emission cannot clone an arbitrary unknown photonic quantum state.

A phase-preserving amplifier must add noise

Section titled “A phase-preserving amplifier must add noise”

An idealized linear phase-preserving power amplifier with G≥1G\ge1, independent input and auxiliary modes, and no input–auxiliary correlations can be represented as

aout=G ain+eiϕG−1 bin†.a_{\mathrm{out}} = \sqrt G\,a_{\mathrm{in}} + e^{i\phi} \sqrt{G-1}\, b_{\mathrm{in}}^\dagger.

The auxiliary mode is necessary because

[aout,aout†]=G−(G−1)=1.\begin{aligned} [a_{\mathrm{out}},a_{\mathrm{out}}^\dagger] &= G - (G-1) \\ &= 1. \end{aligned}

Without the second term, amplification would multiply the canonical commutator by GG. Even when the auxiliary input is vacuum, its fluctuations appear as added amplifier noise. Spontaneous emission into the amplified mode is one physical manifestation of this unavoidable noise in an inverted-medium amplifier.

Experimental Signatures and Forward Models

Section titled “Experimental Signatures and Forward Models”

A weak probe can measure small-signal gain through

Trel(ω)=Iout(pump on)Iout(reference).T_{\mathrm{rel}}(\omega) = \frac{ I_{\mathrm{out}}^{(\mathrm{pump\ on})} }{ I_{\mathrm{out}}^{(\mathrm{reference})} }.

The reference must be declared. Pump-off transmission may include ground-state absorption that disappears when the pump redistributes population. A measured increase can therefore contain both reduced absorption and true stimulated emission. This relative transmission TrelT_{\rm rel} is not automatically the absolute gain G=Iout/IinG=I_{\rm out}/I_{\rm in}. Subtract ASE/background and calibrate the incident seed to establish net amplification of that seed.

In the linear regime,

Iout≈GIin+IASE.I_{\mathrm{out}} \approx GI_{\mathrm{in}} + I_{\mathrm{ASE}}.

Varying IinI_{\mathrm{in}} separates a seed-proportional term from the zero-seed intercept. At higher intensity, gain saturation bends the seed-proportional response.

Lifetime shortening under a resonant field

Section titled “Lifetime shortening under a resonant field”

In an incoherent or coherence-eliminated rate regime, if lower-state reabsorption and repumping are negligible, an applied resonant mode adds a stimulated rate:

Γobs=Γtot+Wstim.\Gamma_{\mathrm{obs}} = \Gamma_{\mathrm{tot}} + W_{\mathrm{stim}}.

An intensity-dependent lifetime can therefore measure a stimulated-emission cross section. In a dense or multilevel medium, radiation trapping, pump redistribution, cooperative emission, and excited-state absorption must be excluded before using this simple relation.

Stimulated output should appear in the phase-matched spatial and polarization mode of the seed, whereas uncollected free-space fluorescence occupies many modes. Finite numerical aperture, birefringence, recoil, collisions, and mode conversion blur this ideal distinction.

Starting from the ladder-operator algebra, find the ratio of:

  1. the total downward squared matrix element for n=24n=24 to that for vacuum;
  2. the stimulated part to the spontaneous part in the same mode;
  3. the absorption squared matrix element to the total downward one at n=24n=24.
Solution

The downward matrix element contains

n+1,\sqrt{n+1},

so its squared magnitude is proportional to n+1n+1. Relative to vacuum,

n+10+1=25.\frac{n+1}{0+1} = 25.

The stimulated and spontaneous pieces are proportional to nn and 11, respectively, so

ΓstimΓsp=n=24.\frac{ \Gamma_{\mathrm{stim}} }{ \Gamma_{\mathrm{sp}} } = n = 24.

Absorption is proportional to nn, while total downward emission is proportional to n+1n+1. Thus

∣Mabs∣2∣Mdown∣2=2425.\frac{ |M_{\mathrm{abs}}|^2 }{ |M_{\mathrm{down}}|^2 } = \frac{24}{25}.

These comparisons are for the same mode, matter matrix element, and line-broadening convention.

For a thermal mode, stimulated and spontaneous downward rates are equal when

n‾(ω,T)=1.\overline n(\omega,T)=1.

Derive the crossover temperature and evaluate it for:

  1. ν=10.0 GHz\nu=10.0\ \mathrm{GHz};
  2. ν=500 THz\nu=500\ \mathrm{THz}.
Solution

Set

1ehν/(kBT)−1=1.\frac{1}{ e^{h\nu/(k_{\mathrm B}T)}-1 } = 1.

Then

ehν/(kBT)=2,e^{h\nu/(k_{\mathrm B}T)} = 2,

so

T×=hνkBln⁡2.T_\times = \frac{ h\nu }{ k_{\mathrm B}\ln2 }.

For 10.0 GHz10.0\ \mathrm{GHz},

hνkB≃0.480 K,\frac{h\nu}{k_{\mathrm B}} \simeq 0.480\ \mathrm K,

and therefore

T×≃0.692 K.T_\times \simeq 0.692\ \mathrm K.

For 500 THz500\ \mathrm{THz},

hνkB≃2.40×104 K,\frac{h\nu}{k_{\mathrm B}} \simeq 2.40\times10^4\ \mathrm K,

giving

T×≃3.46×104 K.T_\times \simeq 3.46\times10^4\ \mathrm K.

Thermal stimulation is therefore easy to encounter for microwave transitions and negligible for a visible-frequency transition at ordinary laboratory temperatures.

A transition connects levels with

g1=2,g2=4.g_1=2, \qquad g_2=4.

The number densities are

N1=8.0×1018 m−3,\mathcal N_1 = 8.0\times10^{18}\ \mathrm{m^{-3}}, N2=1.2×1019 m−3.\mathcal N_2 = 1.2\times10^{19}\ \mathrm{m^{-3}}.

Assume reciprocal level-averaged line profiles and no background loss. Does the medium amplify even though N2>N1\mathcal N_2>\mathcal N_1? Find the minimum upper-level density for gain.

Solution

The relevant populations are per substate:

N1g1=4.0×1018 m−3,\frac{\mathcal N_1}{g_1} = 4.0\times10^{18}\ \mathrm{m^{-3}}, N2g2=3.0×1018 m−3.\frac{\mathcal N_2}{g_2} = 3.0\times10^{18}\ \mathrm{m^{-3}}.

The upper population per substate is smaller, so absorption wins. The fact that the total upper population is larger is insufficient.

Threshold occurs at

N2g2=N1g1.\frac{\mathcal N_2}{g_2} = \frac{\mathcal N_1}{g_1}.

Therefore,

N2,tr=g2g1N1=1.6×1019 m−3.\mathcal N_{2,\mathrm{tr}} = \frac{g_2}{g_1} \mathcal N_1 = 1.6\times10^{19}\ \mathrm{m^{-3}}.

Strict gain requires a value above this transparency density.

At one frequency, a medium has

N2=2.0×1021 m−3,N1=4.0×1020 m−3,\begin{aligned} \mathcal N_2&=2.0\times10^{21}\ \mathrm{m^{-3}}, \\ \mathcal N_1&=4.0\times10^{20}\ \mathrm{m^{-3}}, \end{aligned}

with

σe=2.0×10−20 m2,σa=1.5×10−20 m2,\begin{aligned} \sigma_e&=2.0\times10^{-20}\ \mathrm{m^2}, \\ \sigma_a&=1.5\times10^{-20}\ \mathrm{m^2}, \end{aligned}

and

αbg=1.0 m−1.\alpha_{\mathrm{bg}} = 1.0\ \mathrm{m^{-1}}.

Find the net gain coefficient, the power gain through L=5.0 cmL=5.0\ \mathrm{cm}, and the gain in decibels.

Solution

The three propagation contributions are

N2σe=40 m−1,\mathcal N_2\sigma_e = 40\ \mathrm{m^{-1}}, N1σa=6.0 m−1,\mathcal N_1\sigma_a = 6.0\ \mathrm{m^{-1}},

and

αbg=1.0 m−1.\alpha_{\mathrm{bg}} = 1.0\ \mathrm{m^{-1}}.

Hence

g0=40−6−1=33 m−1.g_0 = 40-6-1 = 33\ \mathrm{m^{-1}}.

The dimensionless gain exponent is

g0L=33(0.050)=1.65.g_0L = 33(0.050) = 1.65.

Therefore,

G=e1.65≃5.21,G = e^{1.65} \simeq 5.21,

and

GdB=10log⁡10(5.21)≃7.17 dB.\begin{aligned} G_{\mathrm{dB}} &= 10\log_{10}(5.21) \\ &\simeq 7.17\ \mathrm{dB}. \end{aligned}

5. Propagation through a saturated amplifier

Section titled “5. Propagation through a saturated amplifier”

Set background loss to zero, and let a uniform amplifier obey g(I)=ggross,0/(1+I/Isat)g(I)=g_{\rm gross,0}/(1+I/I_{\rm sat}). Relate its input and output intensities over length LL. What replaces exponential growth at high input?

Solution

Separate dI/dz=ggross,0I/(1+I/Isat)dI/dz=g_{\rm gross,0}I/(1+I/I_{\rm sat}):

ln⁡IoutIin+Iout−IinIsat=ggross,0L.\ln\frac{I_{\rm out}}{I_{\rm in}}+ \frac{I_{\rm out}-I_{\rm in}}{I_{\rm sat}}=g_{\rm gross,0}L.

When both intensities are small compared with IsatI_{\rm sat}, the second term is negligible and exponential gain is recovered. At high input, dI/dz≃ggross,0IsatdI/dz\simeq g_{\rm gross,0}I_{\rm sat}, so the added intensity tends to ggross,0IsatLg_{\rm gross,0}I_{\rm sat}L rather than growing in proportion to the input. The pump and recovery model fixes the available power per length. This assumes a steady homogeneous inversion and negligible pump depletion.

A linear cavity has

L=0.300 m,L=0.300\ \mathrm m, R1=0.990,R2=0.900,R_1=0.990, \qquad R_2=0.900,

and

αi=0.050 m−1.\alpha_{\mathrm i} = 0.050\ \mathrm{m^{-1}}.

Find the threshold gross material gain coefficient.

Solution

Use

gm,th=αi+12Lln⁡(1R1R2).g_{\mathrm{m,th}} = \alpha_{\mathrm i} + \frac{1}{2L} \ln \left( \frac{1}{R_1R_2} \right).

The mirror product is

R1R2=0.891.R_1R_2 = 0.891.

Therefore,

gm,th=0.050+10.600ln⁡(10.891)≃0.242 m−1.\begin{aligned} g_{\mathrm{m,th}} &= 0.050 + \frac{1}{0.600} \ln \left( \frac{1}{0.891} \right) \\ &\simeq 0.242\ \mathrm{m^{-1}}. \end{aligned}

This is a power-gain coefficient. Using amplitude reflectivities without changing the exponent would introduce a factor-of-two error.

Normalize the photon-added coherent state

∣ψ⟩∝a†∣α⟩|\psi\rangle \propto a^\dagger|\alpha\rangle

and find its fidelity with ∣α⟩|\alpha\rangle. Interpret the limits ∣α∣→0|\alpha|\to0 and ∣α∣→∞|\alpha|\to\infty.

Solution

Its norm is

⟨α∣aa†∣α⟩=⟨α∣(a†a+1)∣α⟩=∣α∣2+1.\begin{aligned} \langle\alpha| aa^\dagger |\alpha\rangle &= \langle\alpha| \left( a^\dagger a+1 \right) |\alpha\rangle \\ &= |\alpha|^2+1. \end{aligned}

Thus

∣ψ⟩=a†∣α⟩∣α∣2+1.|\psi\rangle = \frac{ a^\dagger|\alpha\rangle }{ \sqrt{|\alpha|^2+1} }.

The overlap is

⟨α∣ψ⟩=α∗∣α∣2+1,\langle\alpha|\psi\rangle = \frac{\alpha^*}{ \sqrt{|\alpha|^2+1} },

so the fidelity is

F=∣α∣2∣α∣2+1.F = \frac{|\alpha|^2}{ |\alpha|^2+1 }.

For ∣α∣→0|\alpha|\to0, the operation creates a one-photon state orthogonal to vacuum, so F→0F\to0. For a highly occupied coherent mode, F→1F\to1: one added quantum causes a small relative disturbance. It is still not an exact cloning operation at finite occupation.

8. The added noise of a quantum-limited amplifier

Section titled “8. The added noise of a quantum-limited amplifier”

Use the two-mode amplifier map above with a coherent signal input and an independent vacuum auxiliary input. For X=(a+a†)/2X=(a+a^\dagger)/\sqrt2, find the output mean photon number and quadrature variance. Refer the added quadrature noise to the input by dividing by the power gain.

Solution

The vacuum auxiliary contributes ⟨bb†⟩=1\langle bb^\dagger\rangle=1, giving

⟨nout⟩=G∣α∣2+(G−1).\langle n_{\rm out}\rangle=G|\alpha|^2+(G-1).

Each independent input quadrature has variance 1/21/2, so Var⁡Xout=G/2+(G−1)/2\operatorname{Var}X_{\rm out}=G/2+(G-1)/2. Relative to amplification of the input variance, the added variance is (G−1)/2(G-1)/2 at the output and (1−1/G)/2(1-1/G)/2 referred to the input. It approaches a half quantum at high gain (Caves, 1982). The photon-number offset G−1G-1 and quadrature added noise are different normalizations of the same auxiliary-vacuum contribution.