Skip to content

Einstein Coefficients

Einstein coefficients package three weak-field radiative processes between a lower level 11 and an upper level 22:

R1→2=B12(ω)uω(ω0),R2→1=A21+B21(ω)uω(ω0).\begin{aligned} R_{1\rightarrow2} &= B_{12}^{(\omega)} u_\omega(\omega_0), \\ R_{2\rightarrow1} &= A_{21} + B_{21}^{(\omega)} u_\omega(\omega_0). \end{aligned}

Here

ℏω0=E2−E1>0.\hbar\omega_0=E_2-E_1>0.

The three terms describe:

  • absorption, proportional to B12(ω)B_{12}^{(\omega)};
  • stimulated emission, proportional to B21(ω)B_{21}^{(\omega)};
  • spontaneous emission, proportional to A21A_{21}.

The superscript (ω)(\omega) is not decorative. It declares that uω(ω)u_\omega(\omega) is the total isotropic radiation energy density per unit angular frequency, summed over propagation directions and polarizations. A different spectral variable changes the numerical BB coefficient.

The AA coefficient is a rate. A BB coefficient is not: it becomes a rate only after multiplication by the matching spectral energy density. Neither coefficient is a finite-time Born-rule probability.

This page owns the spectroscopy-facing definitions of the Einstein AA and BB coefficients, the detailed-balance derivation of their relations, their connection to Planck radiation, and the convention ledger needed to compare tabulated or calculated values.

Several nearby pages keep narrower canonical responsibilities:

The Einstein picture is a rate model for weak, incoherent radiation and population kinetics. It is not a replacement for coherent amplitudes, strong-drive dynamics, or environment-specific quantum electrodynamics.

Unless stated otherwise:

  • E2>E1E_2>E_1 and ω0=(E2−E1)/ℏ\omega_0=(E_2-E_1)/\hbar;
  • g1g_1 and g2g_2 are the degeneracies of the lower and upper levels;
  • N1N_1 and N2N_2 are populations or number densities of those levels;
  • populations within each degenerate level are uniformly distributed;
  • uω(ω)u_\omega(\omega) is total isotropic vacuum radiation energy per volume per unit angular frequency;
  • all propagation directions and both transverse polarizations are included in uωu_\omega;
  • the line is sufficiently narrow that the radiation density can be evaluated at ω0\omega_0;
  • the radiation is weak enough for rates to remain linear in spectral energy density;
  • coefficients refer to one declared radiative channel or to an explicitly stated sum of channels.

The simple Einstein relations concern sharp levels in thermal equilibrium. Broad molecular bands, overlapping lines, anisotropic reservoirs, polarization-selected transitions, and structured photonic environments require more detailed spectral or mode-resolved descriptions.

The coefficient A21A_{21} is the spontaneous transition rate per system prepared in level 22 for the specified downward radiative channel:

A21=dP2→1dt∣sp,[A21]=s−1.A_{21} = \left. \frac{dP_{2\rightarrow1}}{dt} \right|_{\mathrm{sp}}, \qquad [A_{21}]=\mathrm{s}^{-1}.

For a population N2N_2, the expected number of photons emitted into that channel per unit time is

N˙γ,21=N2A21\dot N_{\gamma,21} = N_2A_{21}

when reabsorption, stimulated processes, and detector effects are excluded. The corresponding emitted power is

P21=N2ℏω0A21.\mathcal P_{21} = N_2\hbar\omega_0A_{21}.

Atomic databases often call A21A_{21} a transition probability, but its units show that it is a probability rate, not a dimensionless probability.

If upper level 22 can decay radiatively to several lower levels ℓ\ell, its total radiative width in inverse-time units is

Γrad,2=∑ℓA2ℓ.\Gamma_{\mathrm{rad},2} = \sum_\ell A_{2\ell}.

The radiative lifetime is

τrad,2=1Γrad,2.\tau_{\mathrm{rad},2} = \frac{1}{\Gamma_{\mathrm{rad},2}}.

The radiative branching fraction into channel 2→12\rightarrow1 is

b21(rad)=A21∑ℓA2ℓ.b_{21}^{(\mathrm{rad})} = \frac{A_{21}} {\sum_\ell A_{2\ell}}.

If nonradiative, collisional, predissociative, or quenching channels are present,

Γtot,2=Γrad,2+Γnr,2+Γcoll,2+⋯ ,\Gamma_{\mathrm{tot},2} = \Gamma_{\mathrm{rad},2} + \Gamma_{\mathrm{nr},2} + \Gamma_{\mathrm{coll},2} + \cdots,

and

τobs,2=1Γtot,2.\tau_{\mathrm{obs},2} = \frac{1}{\Gamma_{\mathrm{tot},2}}.

The radiative quantum yield is

Φrad=Γrad,2Γtot,2.\Phi_{\mathrm{rad}} = \frac{\Gamma_{\mathrm{rad},2}} {\Gamma_{\mathrm{tot},2}}.

Thus one AA coefficient, a radiative lifetime, an observed lifetime, and a fluorescence quantum yield are four different quantities.

The coefficient B12(ω)B_{12}^{(\omega)} converts spectral energy density into an absorption rate:

R1→2(abs)=B12(ω)uω(ω0).R_{1\rightarrow2}^{(\mathrm{abs})} = B_{12}^{(\omega)} u_\omega(\omega_0).

Similarly,

R2→1(stim)=B21(ω)uω(ω0).R_{2\rightarrow1}^{(\mathrm{stim})} = B_{21}^{(\omega)} u_\omega(\omega_0).

In SI units,

[uω]=J m−3(rad s−1)−1,\left[ u_\omega \right] = \mathrm{J\,m^{-3}} \left(\mathrm{rad\,s^{-1}}\right)^{-1},

so, treating radians as dimensionless,

[B(ω)]=m3 J−1 s−2.\left[ B^{(\omega)} \right] = \mathrm{m^3\,J^{-1}\,s^{-2}}.

Only the products B12(ω)uωB_{12}^{(\omega)}u_\omega and B21(ω)uωB_{21}^{(\omega)}u_\omega have dimensions of inverse time.

For a normalized absorption profile,

∫0∞ϕ12(ω)(ω) dω=1,\int_0^\infty \phi_{12}^{(\omega)}(\omega)\,d\omega =1,

the weak-field rate is more generally

R1→2(abs)=B12(ω)∫0∞uω(ω)ϕ12(ω)(ω) dω.R_{1\rightarrow2}^{(\mathrm{abs})} = B_{12}^{(\omega)} \int_0^\infty u_\omega(\omega) \phi_{12}^{(\omega)}(\omega)\,d\omega.

If uωu_\omega varies negligibly across the line, this reduces to

R1→2(abs)≃B12(ω)uω(ω0).R_{1\rightarrow2}^{(\mathrm{abs})} \simeq B_{12}^{(\omega)} u_\omega(\omega_0).

For spontaneous emission, a normalized profile can distribute the total coefficient:

dA21dω=A21ϕ21(ω)(ω),∫0∞dA21dω dω=A21.\begin{aligned} \frac{dA_{21}}{d\omega} &= A_{21} \phi_{21}^{(\omega)}(\omega), \\ \int_0^\infty \frac{dA_{21}}{d\omega}\,d\omega &= A_{21}. \end{aligned}

The integrated coefficient and the spectral line shape are separate inputs. For broad vibronic bands, absorption and emission profiles need not be identical, and the sharp-line detailed-balance derivation below should not be applied point by point without additional assumptions.

Ordinary frequency versus angular frequency

Section titled “Ordinary frequency versus angular frequency”

Let uν(ν)u_\nu(\nu) be energy density per unit ordinary frequency. Conservation of energy in a spectral interval requires

uν(ν) dν=uω(ω) dω.u_\nu(\nu)\,d\nu = u_\omega(\omega)\,d\omega.

Because

ω=2πν,dω=2π dν,\omega=2\pi\nu, \qquad d\omega=2\pi\,d\nu,

the densities satisfy

uν(ν)=2πuω(ω).u_\nu(\nu) = 2\pi u_\omega(\omega).

The same physical rate can be written

Bij(ν)uν=Bij(ω)uω,B_{ij}^{(\nu)}u_\nu = B_{ij}^{(\omega)}u_\omega,

so

Bij(ν)=12πBij(ω).B_{ij}^{(\nu)} = \frac{1}{2\pi} B_{ij}^{(\omega)}.

An unlabeled BijB_{ij} is therefore incomplete. The factors of 2π2\pi cannot be repaired by dimensional analysis because the two BB conventions have the same SI dimensions when radians are treated as dimensionless.

Einstein’s thermal argument uses isotropic spectral energy density. For isotropic specific intensity IνI_\nu,

uν=4πcIν.u_\nu = \frac{4\pi}{c} I_\nu.

A collimated polarized laser beam is not isotropic. Its irradiance, polarization, propagation direction, and coherence must be converted to the field quantity appropriate to the transition. In coherent spectroscopy, a Rabi frequency or optical Bloch description is usually clearer than inserting beam intensity into BuB u as though it were an isotropic blackbody density.

At temperature TT, thermal equilibrium gives

N2N1=g2g1exp⁡(−ℏω0kBT).\frac{N_2}{N_1} = \frac{g_2}{g_1} \exp\left( -\frac{\hbar\omega_0}{k_{\mathrm B}T} \right).

Define

β=1kBT.\beta = \frac{1}{k_{\mathrm B}T}.

The ratio can then be written

N2N1=g2g1e−βℏω0.\frac{N_2}{N_1} = \frac{g_2}{g_1} e^{-\beta\hbar\omega_0}.

The degeneracies matter because thermal equilibrium populates individual microstates according to Boltzmann weights. A level with more substates has a larger total population at the same energy per substate.

At equilibrium, the total upward transition count equals the total downward count:

N1B12(ω)uω=N2[A21+B21(ω)uω].N_1B_{12}^{(\omega)}u_\omega = N_2 \left[ A_{21} + B_{21}^{(\omega)}u_\omega \right].

Solving for the equilibrium radiation density gives

uω=A21N1N2B12(ω)−B21(ω).u_\omega = \frac{A_{21}} { \dfrac{N_1}{N_2} B_{12}^{(\omega)} - B_{21}^{(\omega)} }.

Using the Boltzmann population ratio,

uω=A21g1g2B12(ω)eβℏω0−B21(ω).u_\omega = \frac{A_{21}} { \dfrac{g_1}{g_2} B_{12}^{(\omega)} e^{\beta\hbar\omega_0} - B_{21}^{(\omega)} }.

Two energy levels with absorption, stimulated emission, spontaneous emission, and thermal-balance relations

Einstein’s three-process rate model for two sharp levels in an isotropic radiation field. Thermal equilibrium fixes the degeneracy relation between the two stimulated coefficients and the ratio of spontaneous to stimulated emission coefficients.

The blackbody energy density per unit angular frequency is

uω(ω,T)=ℏω3π2c31eβℏω−1.u_\omega(\omega,T) = \frac{\hbar\omega^3} {\pi^2c^3} \frac{1} {e^{\beta\hbar\omega}-1}.

For the rate-balance expression to have this form at every temperature, the coefficients must obey

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

and

A21B21(ω)=ℏω03π2c3.\frac{A_{21}} {B_{21}^{(\omega)}} = \frac{\hbar\omega_0^3} {\pi^2c^3}.

In ordinary-frequency notation, the same relations are

g1B12(ν)=g2B21(ν)g_1B_{12}^{(\nu)} = g_2B_{21}^{(\nu)}

and

A21B21(ν)=8πhν03c3.\frac{A_{21}} {B_{21}^{(\nu)}} = \frac{8\pi h\nu_0^3} {c^3}.

The two A/BA/B formulas are physically equivalent. Their different numerical prefactors follow entirely from the definitions of spectral density and BB.

What equilibrium does and does not determine

Section titled “What equilibrium does and does not determine”

Thermodynamic balance determines:

  • the ratio of absorption to stimulated-emission coefficients after degeneracies are included;
  • the ratio of spontaneous to stimulated-emission coefficients after the radiation-density convention is fixed.

It does not determine the absolute magnitude of A21A_{21} or B21(ω)B_{21}^{(\omega)}. That magnitude depends on the transition matrix element. Modern quantum electrodynamics supplies it microscopically.

The first relation is sometimes written B12=B21B_{12}=B_{21}. That is correct for two nondegenerate states or for matched state-to-state coefficients. For level-averaged coefficients,

B12(ω)=g2g1B21(ω).B_{12}^{(\omega)} = \frac{g_2}{g_1} B_{21}^{(\omega)}.

Ignoring this distinction is a common source of factors 2J+12J+1.

For resolved substates connected by the same polarization component, the Hermiticity of the interaction gives equal squared forward and reverse matrix elements. Level-averaged absorption sums over final substates while stimulated emission averages over a different initial manifold. Counting those substates produces

g1B12=g2B21.g_1B_{12}=g_2B_{21}.

The degeneracy relation is therefore consistent with microscopic reversibility; it does not assert that the two level-averaged numbers must be equal.

The thermal mean occupation of a photon mode is

n‾(ω,T)=1eβℏω−1.\overline n(\omega,T) = \frac{1} {e^{\beta\hbar\omega}-1}.

Planck’s angular-frequency energy density can be factored as

uω(ω,T)=u0(ω)(ω)n‾(ω,T),u_\omega(\omega,T) = u_0^{(\omega)}(\omega) \overline n(\omega,T),

where

u0(ω)(ω)=ℏω3π2c3.u_0^{(\omega)}(\omega) = \frac{\hbar\omega^3} {\pi^2c^3}.

This is a convenient free-space mode-density scale, not the divergent zero-point energy density of the electromagnetic vacuum. At the transition frequency it is precisely the ratio

u0(ω)(ω0)=A21B21(ω).u_0^{(\omega)}(\omega_0) = \frac{A_{21}} {B_{21}^{(\omega)}}.

It converts the spontaneous rate to the stimulated coefficient.

Using the Einstein relations,

B21(ω)uω(ω0,T)=A21n‾(ω0,T).B_{21}^{(\omega)} u_\omega(\omega_0,T) = A_{21} \overline n(\omega_0,T).

Thus the downward rate per system in level 22 is

R2→1(th)=A21[n‾(ω0,T)+1].R_{2\rightarrow1}^{(\mathrm{th})} = A_{21} \left[ \overline n(\omega_0,T)+1 \right].

The upward rate per system in level 11 is

R1→2(th)=g2g1A21n‾(ω0,T).R_{1\rightarrow2}^{(\mathrm{th})} = \frac{g_2}{g_1} A_{21} \overline n(\omega_0,T).

For equal degeneracies, these reduce to the familiar quantum-optical expressions proportional to n‾\overline n and n‾+1\overline n+1.

The identity

n‾+1=eβℏω0n‾\overline n+1 = e^{\beta\hbar\omega_0} \overline n

ensures that the Boltzmann population ratio makes total upward and downward event counts equal.

When

ℏω0≫kBT,\hbar\omega_0\gg k_{\mathrm B}T,

the occupation is

n‾≃e−βℏω0≪1.\overline n \simeq e^{-\beta\hbar\omega_0}\ll1.

Stimulated processes from a thermal field are then negligible compared with spontaneous emission. This is typical for visible transitions at room temperature.

When

ℏω0≪kBT,\hbar\omega_0\ll k_{\mathrm B}T,

the occupation is

n‾≃kBTℏω0≫1.\overline n \simeq \frac{k_{\mathrm B}T} {\hbar\omega_0}\gg1.

Stimulated emission and absorption dominate over the extra spontaneous +1+1. Microwave and radio-frequency transitions can therefore be strongly affected by thermal photons even at ordinary laboratory temperatures.

The stimulated and spontaneous downward contributions are equal when

n‾=1,\overline n=1,

or

T=ℏω0kBln⁡2.T = \frac{\hbar\omega_0} {k_{\mathrm B}\ln 2}.

Einstein introduced the three-process rate structure in 1916–1917 and used thermal equilibrium with the known radiation law to constrain the coefficients. Stimulated emission emerged as a necessary process: absorption and spontaneous emission alone cannot reproduce Planck’s denominator.

The argument was extraordinarily predictive, but it was not yet a modern quantized-field calculation of A21A_{21}. Dirac’s 1927 radiation theory and later quantum electrodynamics supplied a microscopic account of emission and absorption amplitudes. The historical development should not be compressed into the claim that equilibrium thermodynamics alone calculated every coefficient.

For a nondegenerate electric-dipole transition in vacuum, define

d21=⟨2∣D^∣1⟩.\mathbf d_{21} = \langle2| \widehat{\mathbf D} |1\rangle.

The free-space spontaneous-emission coefficient is

A21=ω033πϵ0ℏc3∣d21∣2.A_{21} = \frac{\omega_0^3} {3\pi\epsilon_0\hbar c^3} \left| \mathbf d_{21} \right|^2.

Combining this with the angular-frequency Einstein relation gives

B21(ω)=π3ϵ0ℏ2∣d21∣2.B_{21}^{(\omega)} = \frac{\pi} {3\epsilon_0\hbar^2} \left| \mathbf d_{21} \right|^2.

For two nondegenerate states,

B12(ω)=B21(ω).B_{12}^{(\omega)} = B_{21}^{(\omega)}.

In ordinary-frequency notation,

B21(ν)=16ϵ0ℏ2∣d21∣2.B_{21}^{(\nu)} = \frac{1}{6\epsilon_0\hbar^2} \left| \mathbf d_{21} \right|^2.

These equations display the convention dependence directly:

B21(ω)=2πB21(ν).B_{21}^{(\omega)} = 2\pi B_{21}^{(\nu)}.

For an E1 transition between levels with reduced line strength

S21=∣⟨γ2J2∥D^(1)∥γ1J1⟩∣2,S_{21} = \left| \left\langle \gamma_2J_2 \left\| \widehat D^{(1)} \right\| \gamma_1J_1 \right\rangle \right|^2,

the level-averaged coefficients are

A21=ω033πϵ0ℏc3g2S21,B21(ω)=π3ϵ0ℏ2g2S21,B12(ω)=π3ϵ0ℏ2g1S21.\begin{aligned} A_{21} &= \frac{\omega_0^3} {3\pi\epsilon_0\hbar c^3g_2} S_{21}, \\ B_{21}^{(\omega)} &= \frac{\pi} {3\epsilon_0\hbar^2g_2} S_{21}, \\ B_{12}^{(\omega)} &= \frac{\pi} {3\epsilon_0\hbar^2g_1} S_{21}. \end{aligned}

The first coefficient averages over upper substates; the absorption coefficient averages over lower substates. The common reduced line strength makes the degeneracy relation manifest.

For an E1 absorption oscillator strength f12f_{12},

A21=e2ω022πϵ0mec3g1g2f12.A_{21} = \frac{e^2\omega_0^2} {2\pi\epsilon_0m_ec^3} \frac{g_1}{g_2} f_{12}.

The stimulated coefficients are

B12(ω)=πe22ϵ0meℏω0f12,B21(ω)=g1g2B12(ω).\begin{aligned} B_{12}^{(\omega)} &= \frac{\pi e^2} {2\epsilon_0m_e\hbar\omega_0} f_{12}, \\ B_{21}^{(\omega)} &= \frac{g_1}{g_2} B_{12}^{(\omega)}. \end{aligned}

The complete convention and line-strength ledger belongs to Oscillator Strengths.

For a fixed transition dipole,

A21∝ω03,B21(ω)∝∣d21∣2.A_{21}\propto\omega_0^3, \qquad B_{21}^{(\omega)} \propto \left|\mathbf d_{21}\right|^2.

The spontaneous rate carries the free-space photon mode-density factor ω02\omega_0^2 and one additional frequency from field normalization. The angular-frequency BB coefficient does not carry that ω03\omega_0^3 factor. The relation

A21B21(ω)∝ω03\frac{A_{21}} {B_{21}^{(\omega)}} \propto \omega_0^3

therefore reflects the available free-space radiation modes, not a temperature dependence of either intrinsic coefficient.

Einstein AA coefficients are also tabulated for M1, E2, and higher-multipole transitions, but their microscopic line-strength formulas differ from the E1 expressions above. The equilibrium A/BA/B relation for a sharp one-photon transition still relies on the declared radiation-mode convention.

In a cavity, waveguide, dielectric, near a surface, or in another structured photonic environment, spontaneous emission depends on the local electromagnetic density of states and field mode functions. The free-space formula can be enhanced, suppressed, or made anisotropic. One should then use an environment-specific mode sum or electromagnetic Green tensor rather than silently reusing the vacuum A21A_{21}. See Spontaneous Emission for the Green-tensor and Purcell-factor preview.

In a radiation field, the net stimulated downward event rate per volume is

Rstim,net=N2B21(ω)uω−N1B12(ω)uω=B21(ω)uω(N2−g2g1N1).\begin{aligned} \mathcal R_{\mathrm{stim,net}} &= N_2B_{21}^{(\omega)}u_\omega - N_1B_{12}^{(\omega)}u_\omega \\ &= B_{21}^{(\omega)}u_\omega \left( N_2 - \frac{g_2}{g_1}N_1 \right). \end{aligned}

Net gain requires

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

This is population inversion per substate. The simpler condition N2>N1N_2>N_1 is sufficient only when g1=g2g_1=g_2.

A thermal two-level ensemble cannot satisfy this inequality because

N2/g2N1/g1=e−βℏω0<1.\frac{N_2/g_2}{N_1/g_1} = e^{-\beta\hbar\omega_0} <1.

Pumping and additional levels are therefore needed to maintain inversion in a laser medium.

Spontaneous emission:

  • depletes the upper laser level;
  • seeds photons with random emission times, directions, phases, and polarizations;
  • contributes quantum noise;
  • competes with nonradiative decay and cavity loss;
  • supplies fluorescence even below threshold.

Stimulated emission into a selected mode copies that mode’s frequency, polarization, propagation direction, and phase relation. The Einstein rate picture captures population transfer and small-signal gain, but a coherent laser field and its phase dynamics require Maxwell–Bloch, master-equation, or quantum-optical methods.

The Einstein equations assume that optical coherence can be neglected and that transition events can be described by approximately constant rates. They become inadequate when:

  • a coherent pulse drives resolved Rabi oscillations;
  • the field is strong enough to saturate or power broaden the transition;
  • coherence between degenerate substates matters;
  • one cavity mode exchanges excitations reversibly with the emitter;
  • reservoir memory or strong light–matter coupling invalidates Markovian decay;
  • interference among multiple pathways changes the transition amplitude.

The appropriate coherent two-level treatment is developed in Optical Bloch Equations.

Spontaneous emission is not classical driving

Section titled “Spontaneous emission is not classical driving”

A prescribed classical field can produce absorption and stimulated emission, but it does not derive spontaneous emission from an initially empty field. In the quantized-field description, downward coupling to a photon mode has the matrix-element factor

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

while absorption has

n.\sqrt n.

The surviving unity at n=0n=0, summed over vacuum modes, produces the spontaneous rate. The full derivation and free-space dipole radiation pattern belong to Transition Rates in Light–Matter Interaction.

B(ν)B^{(\nu)} and B(ω)B^{(\omega)} differ by 2π2\pi. Additional conventions may use energy per solid angle, per polarization, or per wavenumber. State the field density explicitly.

BB requires multiplication by a compatible spectral energy density. Its units are not inverse seconds.

Setting B12 equal to B21 for degenerate levels

Section titled “Setting B12 equal to B21 for degenerate levels”

State-to-state coefficients can be equal, but level-averaged coefficients obey g1B12=g2B21g_1B_{12}=g_2B_{21}.

A21A_{21} is a rate. Over a sufficiently short rate interval,

P2→1(t)≃A21t,P_{2\rightarrow1}(t) \simeq A_{21}t,

but the exact finite-time probability and depletion law require a dynamical model.

Confusing radiative and observed lifetimes

Section titled “Confusing radiative and observed lifetimes”

The inverse observed lifetime contains all decay channels. It equals ∑ℓA2ℓ\sum_\ell A_{2\ell} only when nonradiative and collisional losses are negligible.

Inserting beam irradiance into an isotropic formula

Section titled “Inserting beam irradiance into an isotropic formula”

Irradiance, radiance, energy density, and spectral energy density are different quantities. Directional polarization coupling cannot be recovered from total isotropic uωu_\omega after the information has been discarded.

Assuming equilibrium fixes absolute strengths

Section titled “Assuming equilibrium fixes absolute strengths”

Detailed balance fixes coefficient ratios. A matrix element or measured quantity is needed to determine their absolute scale.

Using free-space A in a structured reservoir

Section titled “Using free-space A in a structured reservoir”

Cavities, interfaces, and photonic materials modify the available modes. Their emission rates are not universal properties of the isolated emitter alone.

Applying rate equations to coherent strong driving

Section titled “Applying rate equations to coherent strong driving”

When phase coherence, Rabi cycling, or saturation matters, use amplitudes, density matrices, or optical Bloch equations.

When using Einstein coefficients:

  1. identify the lower and upper levels and their degeneracies;
  2. declare whether the coefficients are state-resolved, level-averaged, or multiplet-averaged;
  3. specify ordinary frequency, angular frequency, wavelength, wavenumber, or another spectral coordinate;
  4. define whether the field density includes all directions and polarizations;
  5. normalize any absorption and emission profiles;
  6. distinguish one AA channel from the sum controlling a lifetime;
  7. include nonradiative, collisional, transfer, and reabsorption processes when comparing with measured kinetics;
  8. use dipole moments, line strengths, oscillator strengths, and database values only after matching degeneracy and unit conventions;
  9. replace the rate model by coherent or environment-specific dynamics when its assumptions fail.
  • A21A_{21} is a spontaneous-emission rate; BijB_{ij} becomes a stimulated rate only after multiplication by a matching spectral energy density.
  • The numerical value of a BB coefficient depends on the spectral-density convention.
  • Level-averaged coefficients satisfy g1B12=g2B21g_1B_{12}=g_2B_{21}.
  • In the angular-frequency convention, A21/B21(ω)=ℏω03/(π2c3)A_{21}/B_{21}^{(\omega)}=\hbar\omega_0^3/(\pi^2c^3).
  • The same relations reproduce Planck’s law in thermal equilibrium but do not determine the absolute transition strength.
  • Free-space E1 coefficients are fixed by one transition dipole or reduced line strength.
  • Radiative lifetime requires a sum over all spontaneous channels; measured lifetime can include additional losses.
  • Laser gain requires inversion per degenerate substate, not merely more particles in the upper level.
  • Strong coherent driving and structured reservoirs lie beyond the simple Einstein rate model.

Exercise 1: Convert B between spectral conventions

Section titled “Exercise 1: Convert B between spectral conventions”

Show that

B21(ν)=B21(ω)2πB_{21}^{(\nu)} = \frac{B_{21}^{(\omega)}}{2\pi}

when uνu_\nu and uωu_\omega describe the same radiation field. Verify that the stimulated rate is unchanged.

Solution

The spectral energy in a physical interval is invariant:

uν dν=uω dω.u_\nu\,d\nu = u_\omega\,d\omega.

Since dω=2π dνd\omega=2\pi\,d\nu,

uν=2πuω.u_\nu=2\pi u_\omega.

Rate invariance requires

B21(ν)uν=B21(ω)uω.B_{21}^{(\nu)}u_\nu = B_{21}^{(\omega)}u_\omega.

Substitution gives

B21(ν)(2πuω)=B21(ω)uω,B_{21}^{(\nu)} (2\pi u_\omega) = B_{21}^{(\omega)}u_\omega,

and therefore

B21(ν)=B21(ω)2π.B_{21}^{(\nu)} = \frac{B_{21}^{(\omega)}}{2\pi}.

The coefficient and density change inversely, leaving their product unchanged.

Starting from thermal population balance and Planck’s angular-frequency energy density, derive the two Einstein relations for levels with degeneracies g1g_1 and g2g_2.

Solution

Equilibrium requires

N1B12(ω)uω=N2(A21+B21(ω)uω).N_1B_{12}^{(\omega)}u_\omega = N_2 \left( A_{21} + B_{21}^{(\omega)}u_\omega \right).

With

N2N1=g2g1e−βℏω0,\frac{N_2}{N_1} = \frac{g_2}{g_1} e^{-\beta\hbar\omega_0},

solving for uωu_\omega gives

uω=A21g1g2B12(ω)eβℏω0−B21(ω).u_\omega = \frac{A_{21}} { \dfrac{g_1}{g_2} B_{12}^{(\omega)} e^{\beta\hbar\omega_0} - B_{21}^{(\omega)} }.

To match a denominator proportional to eβℏω0−1e^{\beta\hbar\omega_0}-1, require

g1g2B12(ω)=B21(ω).\frac{g_1}{g_2} B_{12}^{(\omega)} = B_{21}^{(\omega)}.

Hence

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

The remaining prefactor must match Planck’s law:

A21B21(ω)=ℏω03π2c3.\frac{A_{21}} {B_{21}^{(\omega)}} = \frac{\hbar\omega_0^3} {\pi^2c^3}.

An E1 line connects

J1=12,J2=32.J_1=\frac{1}{2}, \qquad J_2=\frac{3}{2}.

Find B12(ω)/B21(ω)B_{12}^{(\omega)}/B_{21}^{(\omega)} for level-averaged coefficients.

Solution

The statistical weights are

g1=2J1+1=2,g2=2J2+1=4.\begin{aligned} g_1&=2J_1+1=2, \\ g_2&=2J_2+1=4. \end{aligned}

Detailed balance gives

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

Therefore

B12(ω)B21(ω)=g2g1=2.\frac{B_{12}^{(\omega)}} {B_{21}^{(\omega)}} = \frac{g_2}{g_1} =2.

This factor does not violate microscopic reversibility. It arises because the absorption coefficient averages over two lower substates while the stimulated-emission coefficient averages over four upper substates.

Exercise 4: Thermal photons at optical and microwave frequencies

Section titled “Exercise 4: Thermal photons at optical and microwave frequencies”

Evaluate the thermal occupation

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

at T=300 KT=300\ \mathrm K for:

  1. ν=5.00×1014 Hz\nu=5.00\times10^{14}\ \mathrm{Hz};
  2. ν=1.00×1010 Hz\nu=1.00\times10^{10}\ \mathrm{Hz}.

Interpret the ratio of stimulated to spontaneous downward rates.

Solution

For the optical transition,

hνkBT≃79.99,\frac{h\nu}{k_{\mathrm B}T} \simeq 79.99,

so

n‾opt≃1.83×10−35.\overline n_{\mathrm{opt}} \simeq 1.83\times10^{-35}.

Thermal stimulated emission is negligible compared with spontaneous emission.

For the microwave transition,

hνkBT≃1.600×10−3,\frac{h\nu}{k_{\mathrm B}T} \simeq 1.600\times10^{-3},

so

n‾mw≃6.25×102.\overline n_{\mathrm{mw}} \simeq 6.25\times10^2.

Because

RstimRsp=n‾,\frac{R_{\mathrm{stim}}}{R_{\mathrm{sp}}} = \overline n,

thermal stimulated emission dominates for the microwave transition. The same thermal field also drives strong absorption.

Exercise 5: Lifetime, branching, and quantum yield

Section titled “Exercise 5: Lifetime, branching, and quantum yield”

An upper level has two radiative channels,

A21=2.0×107 s−1,A20=6.0×107 s−1,\begin{aligned} A_{21}&=2.0\times10^7\ \mathrm{s}^{-1}, \\ A_{20}&=6.0\times10^7\ \mathrm{s}^{-1}, \end{aligned}

and a nonradiative decay rate

Γnr=2.0×107 s−1.\Gamma_{\mathrm{nr}} = 2.0\times10^7\ \mathrm{s}^{-1}.

Find the radiative lifetime, observed lifetime, radiative quantum yield, and the total-event branching probabilities.

Solution

The radiative rate is

Γrad=A21+A20=8.0×107 s−1,\Gamma_{\mathrm{rad}} = A_{21}+A_{20} = 8.0\times10^7\ \mathrm{s}^{-1},

so

τrad=12.5 ns.\tau_{\mathrm{rad}} = 12.5\ \mathrm{ns}.

The total rate is

Γtot=1.0×108 s−1,\Gamma_{\mathrm{tot}} = 1.0\times10^8\ \mathrm{s}^{-1},

giving

τobs=10.0 ns.\tau_{\mathrm{obs}} = 10.0\ \mathrm{ns}.

The radiative quantum yield is

Φrad=8.0×1071.0×108=0.80.\Phi_{\mathrm{rad}} = \frac{8.0\times10^7}{1.0\times10^8} =0.80.

The probabilities that a decay event uses the three channels are

p21=0.20,p20=0.60,pnr=0.20.\begin{aligned} p_{21}&=0.20, \\ p_{20}&=0.60, \\ p_{\mathrm{nr}}&=0.20. \end{aligned}

Conditioned on a radiative event, the two radiative branching fractions are 0.250.25 and 0.750.75.

Exercise 6: A coefficient from a transition dipole

Section titled “Exercise 6: A coefficient from a transition dipole”

A nondegenerate E1 transition in vacuum has

λ=500 nm,∣d21∣=ea0.\lambda=500\ \mathrm{nm}, \qquad \left|\mathbf d_{21}\right| = ea_0.

Estimate A21A_{21} and the single-channel radiative lifetime. Use ea0=8.478×10−30 C mea_0=8.478\times10^{-30}\ \mathrm{C\,m}.

Solution

The angular frequency is

ω0=2πcλ≃3.767×1015 s−1.\omega_0 = \frac{2\pi c}{\lambda} \simeq 3.767\times10^{15}\ \mathrm{s}^{-1}.

Then

A21=ω03∣d21∣23πϵ0ℏc3≃1.62×107 s−1.A_{21} = \frac{\omega_0^3|\mathbf d_{21}|^2} {3\pi\epsilon_0\hbar c^3} \simeq 1.62\times10^7\ \mathrm{s}^{-1}.

If this is the only decay channel,

τrad=1A21≃6.17×10−8 s=61.7 ns.\tau_{\mathrm{rad}} = \frac{1}{A_{21}} \simeq 6.17\times10^{-8}\ \mathrm s = 61.7\ \mathrm{ns}.

A transition has g1=2g_1=2 and g2=6g_2=6. Determine whether each population pair provides net stimulated gain:

N1N2A3.0×10155.0×1015B3.0×10151.2×1016\begin{array}{c|cc} &N_1&N_2\\ \hline \mathrm{A} &3.0\times10^{15} &5.0\times10^{15} \\ \mathrm{B} &3.0\times10^{15} &1.2\times10^{16} \end{array}
Solution

The inversion criterion is

N2g2>N1g1,\frac{N_2}{g_2} > \frac{N_1}{g_1},

or equivalently

N2>g2g1N1=3N1.N_2 > \frac{g_2}{g_1}N_1 = 3N_1.

For case A, gain would require N2>9.0×1015N_2>9.0\times10^{15}, but N2=5.0×1015N_2=5.0\times10^{15}. The medium absorbs.

For case B, N2=1.2×1016>9.0×1015N_2=1.2\times10^{16}>9.0\times10^{15}. The stimulated contribution gives net gain.

Exercise 8: Audit a laser-beam calculation

Section titled “Exercise 8: Audit a laser-beam calculation”

A calculation takes a tabulated B12(ω)B_{12}^{(\omega)}, multiplies it directly by a laser irradiance in W m−2\mathrm{W\,m^{-2}}, and calls the result an absorption probability. Identify the problems and state what is needed.

Solution

There are several independent errors:

  1. B12(ω)B_{12}^{(\omega)} multiplies spectral energy density per angular frequency, not irradiance.
  2. The tabulated coefficient may already average over directions, polarizations, and magnetic substates that do not match the laser preparation.
  3. A monochromatic coherent beam requires a linewidth or spectral profile; total irradiance alone does not specify uω(ω0)u_\omega(\omega_0).
  4. The product BuωB u_\omega is a rate, not a dimensionless probability.
  5. If coherent evolution or saturation matters, the Einstein rate model is insufficient.

One must specify the beam spectrum, polarization, propagation direction, field normalization, transition profile, level degeneracies, and interaction time. In the weak incoherent limit, convert the beam to a compatible spectral photon flux or energy density and calculate a rate. For coherent narrowband driving, use the transition dipole to form a Rabi frequency and solve the appropriate density-matrix dynamics.