Skip to content

Einstein Coefficient Reference

Einstein coefficients describe absorption, stimulated emission, and spontaneous emission between two declared levels in a weak-field rate model. For a lower level 11, an upper level 22, and a narrow transition at ω0=(E2−E1)/ℏ>0\omega_0=(E_2-E_1)/\hbar>0,

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

Here uωu_\omega is the total isotropic radiation energy density per unit angular frequency. The superscript on BB is essential: a BB coefficient has a numerical value only relative to the spectral density, directional average, polarization sum, and line-profile convention that it multiplies.

This page is a lookup and audit sheet. The detailed-balance derivation and physical discussion live in Einstein Coefficients.

This page owns:

  • operational definitions and units of A21A_{21}, B12(ω)B_{12}^{(\omega)}, and B21(ω)B_{21}^{(\omega)};
  • the degeneracy and Planck-density relations among the three coefficients;
  • conversion among angular-frequency, ordinary-frequency, wavenumber, and wavelength BB conventions;
  • radiative widths, lifetimes, branching fractions, observed lifetimes, and quantum yields;
  • E1 conversion to transition dipoles, reduced line strengths, and oscillator strengths;
  • compact atomic and molecular database checks; and
  • applicability tests for rate equations.

It does not own:

  • the historical or modern derivation of Planck’s radiation law;
  • mode-quantization derivations of the nn and n+1n+1 factors;
  • Wigner–Weisskopf decay and photon wavepackets;
  • optical Bloch dynamics under coherent driving;
  • detailed line-shape theory;
  • Purcell enhancement and Green-tensor calculations; or
  • the derivation of forbidden-multipole line-strength formulas.

The reference formulas are valid only with the conventions stated beside them.

Unless a section says otherwise:

  • E2>E1E_2>E_1 and ω0=(E2−E1)/ℏ\omega_0=(E_2-E_1)/\hbar;
  • g1g_1 and g2g_2 are the total degeneracies of the lower and upper objects;
  • populations are uniform within each degenerate level;
  • A21A_{21} describes one declared upper-to-lower radiative channel;
  • uω(ω)u_\omega(\omega) is total isotropic vacuum radiation energy per volume per angular-frequency interval, summed over directions and both transverse polarizations;
  • line profiles are normalized to unit area in the variable printed in their subscript;
  • the narrow-line approximation is used only when the radiation density is nearly constant across the profile; and
  • microscopic matrix-element formulas refer to an electric-dipole, or E1, transition in free space.

The three coefficients are intrinsic to the declared transition model. A rate additionally requires a radiation field or a prepared upper-state population.

QuantityOperational definitionSI unitsCommon database meaning
A21A_{21}spontaneous downward rate per system in level 22s−1\mathrm{s^{-1}}emission transition probability per unit time
B12(ω)B_{12}^{(\omega)}absorption rate divided by matching uωu_\omegam3 J−1 s−2\mathrm{m^3\,J^{-1}\,s^{-2}}usually derived rather than tabulated
B21(ω)B_{21}^{(\omega)}stimulated-emission rate divided by matching uωu_\omegam3 J−1 s−2\mathrm{m^3\,J^{-1}\,s^{-2}}convention dependent
g2A21g_2A_{21}upper-weighted spontaneous rates−1\mathrm{s^{-1}}a distinct column from A21A_{21}
∑ℓA2ℓ\sum_\ell A_{2\ell}total radiative decay rate out of level 22s−1\mathrm{s^{-1}}inverse radiative lifetime

For a system prepared in level 22,

A21=dP2→1dt∣sp.A_{21} = \left. \frac{dP_{2\rightarrow1}}{dt} \right|_{\mathrm{sp}}.

For N2N_2 independent emitters, and before collection or detector losses,

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

and the emitted optical power in that channel is

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

Atomic databases often call A21A_{21} a transition probability. Its units show that it is a probability rate, not a dimensionless finite-time probability.

The absorption and stimulated-emission rates per system are

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

With

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

and radians treated as dimensionless,

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

A BB coefficient alone is not a rate. It becomes one only after multiplication by the matching spectral energy density.

Let the absorption profile be normalized by

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

The weak-field rate is

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.

Only when uωu_\omega changes negligibly across the line may one write

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

The spontaneous coefficient may likewise be distributed over a normalized emission profile:

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

The integrated AA and the line profile are separate objects. Broad vibronic bands, overlapping transitions, or rapidly varying radiation fields require the integral rather than a line-center shortcut.

Spectral energy is invariant under a coordinate change:

ux(x) dx=uy(y) dy.u_x(x)\,dx = u_y(y)\,dy.

The numerical BB must change inversely so that the physical rate remains invariant:

B(x)ux=B(y)uy.B^{(x)}u_x = B^{(y)}u_y.

Because

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

the coefficients satisfy

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

The two conventions have the same reduced SI dimensions when radians are dimensionless, so dimensional analysis does not reveal a missing 2π2\pi.

For ν~=ν/c\widetilde\nu=\nu/c measured in m−1\mathrm{m^{-1}},

uν~=c uν,Bij(ν~)=1cBij(ν).u_{\widetilde\nu} = c\,u_\nu, \qquad B_{ij}^{(\widetilde\nu)} = \frac{1}{c} B_{ij}^{(\nu)}.

If wavenumber is measured in cm−1\mathrm{cm^{-1}}, then

Bij(cm−1)=1100cBij(ν).B_{ij}^{(\mathrm{cm}^{-1})} = \frac{1}{100c} B_{ij}^{(\nu)}.

The unit attached to the wavenumber variable is part of the coefficient definition.

Since ν=c/λ\nu=c/\lambda,

uλ(λ)=cλ2uν(ν).u_\lambda(\lambda) = \frac{c}{\lambda^2} u_\nu(\nu).

At a narrow line centered on λ0\lambda_0,

Bij(λ)=λ02cBij(ν).B_{ij}^{(\lambda)} = \frac{\lambda_0^2}{c} B_{ij}^{(\nu)}.

For a broad band, the Jacobian varies across the profile. A single line-center factor is then insufficient.

Density used in the rateCoordinate unitsCoefficient relative to B(ω)B^{(\omega)}
uωu_\omegarad s−1\mathrm{rad\,s^{-1}}B(ω)B^{(\omega)}
uνu_\nuHz\mathrm{Hz}B(ν)=B(ω)/(2π)B^{(\nu)}=B^{(\omega)}/(2\pi)
uν~u_{\widetilde\nu}m−1\mathrm{m^{-1}}B(ν~)=B(ω)/(2πc)B^{(\widetilde\nu)}=B^{(\omega)}/(2\pi c)
uν~u_{\widetilde\nu}cm−1\mathrm{cm^{-1}}B(cm−1)=B(ω)/(200πc)B^{(\mathrm{cm}^{-1})}=B^{(\omega)}/(200\pi c)
uλu_\lambda near λ0\lambda_0m\mathrm mB(λ)=λ02B(ω)/(2πc)B^{(\lambda)}=\lambda_0^2B^{(\omega)}/(2\pi c)

For sharp levels in free space and the conventions above, thermal detailed balance gives

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

Thus

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

The two level-averaged BB coefficients are equal only when the degeneracies are equal. State-to-state coefficients for matched resolved substates can be equal even when the level-averaged coefficients are not.

The spontaneous-to-stimulated relation is

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

In ordinary-frequency notation,

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

These are the same physical relation expressed with different spectral densities. Thermodynamic balance fixes ratios among coefficients; it does not determine their absolute scale. The transition matrix element supplies that scale.

The mean occupation of a thermal photon mode is

n‾(ω0,T)=1exp⁡ ⁣(ℏω0kBT)−1.\overline n(\omega_0,T) = \frac{1} { \exp\!\left( \dfrac{\hbar\omega_0}{k_{\mathrm B}T} \right)-1 }.

Planck’s energy density per angular frequency is

uω(ω0,T)=ℏω03π2c3n‾.u_\omega(\omega_0,T) = \frac{\hbar\omega_0^3} {\pi^2c^3} \overline n.

Using the Einstein relations gives the compact rates

R2→1(th)=A21(n‾+1),R_{2\rightarrow1}^{(\mathrm{th})} = A_{21}(\overline n+1), R1→2(th)=g2g1A21n‾.R_{1\rightarrow2}^{(\mathrm{th})} = \frac{g_2}{g_1} A_{21}\overline n.

The +1+1 is the spontaneous contribution. The thermal population ratio

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

ensures

N1R1→2(th)=N2R2→1(th).N_1R_{1\rightarrow2}^{(\mathrm{th})} = N_2R_{2\rightarrow1}^{(\mathrm{th})}.

For optical transitions at room temperature, n‾\overline n is usually negligible. For microwave and radio-frequency transitions it can be much larger than one.

If level 22 has several radiative channels,

Γ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}}.

One channel coefficient equals the inverse radiative lifetime only when it is the sole radiative decay channel.

The radiative branching fraction into channel 2→ℓ2\rightarrow\ell is

b2ℓ(rad)=A2ℓ∑jA2j,∑ℓb2ℓ(rad)=1.b_{2\ell}^{(\mathrm{rad})} = \frac{A_{2\ell}} {\sum_jA_{2j}}, \qquad \sum_\ell b_{2\ell}^{(\mathrm{rad})} = 1.

Conversely, if a radiative lifetime and complete radiative branching fractions are known,

A2ℓ=b2ℓ(rad)τrad,2.A_{2\ell} = \frac{ b_{2\ell}^{(\mathrm{rad})} }{ \tau_{\mathrm{rad},2} }.

With additional first-order channels,

Γ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=τobs,2τrad,2.\Phi_{\mathrm{rad}} = \frac{\Gamma_{\mathrm{rad},2}} {\Gamma_{\mathrm{tot},2}} = \frac{\tau_{\mathrm{obs},2}} {\tau_{\mathrm{rad},2}}.

The last equality assumes the same state and first-order kinetics. In multiexponential, transfer-coupled, diffusion-limited, or time-dependent environments, no single reciprocal rate need describe the measured trace.

For a stable lower state, no pure dephasing, and exponential upper-state decay,

Γ1=∑ℓA2ℓ,Γ2=Γ12.\Gamma_1 = \sum_\ell A_{2\ell}, \qquad \Gamma_2 = \frac{\Gamma_1}{2}.

The Lorentzian full width at half maximum is then

ΔωFWHM=Γ1,\Delta\omega_{\mathrm{FWHM}} = \Gamma_1,

or

ΔνFWHM=Γ12π.\Delta\nu_{\mathrm{FWHM}} = \frac{\Gamma_1}{2\pi}.

More generally,

Γ2=Γupper+Γlower2+γϕ,\Gamma_2 = \frac{\Gamma_{\mathrm{upper}} +\Gamma_{\mathrm{lower}}}{2} + \gamma_\phi,

so

ΔωFWHM=2Γ2.\Delta\omega_{\mathrm{FWHM}} = 2\Gamma_2.

The natural-width conversion therefore uses the sum of all relevant AA channels and any lower-state width or pure dephasing. The complete width ledger belongs to Line Shape Reference.

For

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

the free-space E1 coefficients are

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, B21(ω)=π3ϵ0ℏ2∣d21∣2.B_{21}^{(\omega)} = \frac{\pi} {3\epsilon_0\hbar^2} \left|\mathbf d_{21}\right|^2.

For nondegenerate endpoints,

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

Degenerate levels and reduced line strength

Section titled “Degenerate levels and reduced line strength”

For a reduced E1 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 spontaneous and stimulated-emission coefficients average over upper substates; the absorption coefficient averages over lower substates.

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 angular-frequency stimulated coefficients are

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

For vacuum wavelength in ångströms and A21A_{21} in s−1\mathrm{s^{-1}},

f12=1.49919×10−16g2g1λA˚ 2A21.f_{12} = 1.49919\times10^{-16} \frac{g_2}{g_1} \lambda_{\mathring{\mathrm A}}^{\,2} A_{21}.

For S21S_{21} in e2a02e^2a_0^2,

A21=2.02613×1018g2λA˚ 3S21.A_{21} = \frac{2.02613\times10^{18}} {g_2\lambda_{\mathring{\mathrm A}}^{\,3}} S_{21}.

The complete ff, gfgf, log⁡(gf)\log(gf), and SS ledger is in Oscillator Strength Reference.

For a fixed E1 line strength,

A21∝ω03,B21(ω)∝S21.A_{21}\propto\omega_0^3, \qquad B_{21}^{(\omega)} \propto S_{21}.

For a fixed oscillator strength, however,

A21∝ω02,B12(ω)∝ω0−1.A_{21}\propto\omega_0^2, \qquad B_{12}^{(\omega)} \propto\omega_0^{-1}.

Always state which structural quantity is being held fixed before quoting a frequency scaling.

Einstein AA values are also used for M1, E2, and higher-multipole transitions. Their microscopic operators, dimensions, angular factors, and frequency powers differ from the E1 formulas above. The transition type must be read before converting a database line strength.

In a cavity, waveguide, dielectric, photonic crystal, near a surface, or in another structured reservoir, spontaneous emission depends on the local electromagnetic density of states and mode functions. A vacuum AA can be enhanced, suppressed, redirected, or split among channels. Use a mode-resolved calculation or electromagnetic Green tensor rather than silently applying the free-space formula.

The equilibrium relations among upward and downward processes also assume a specified reservoir and consistent mode normalization. They are not a universal recipe for replacing an environment-dependent decay rate while leaving every other coefficient unchanged.

Atomic line lists commonly use lower level ii and upper level kk:

i→kfor absorption,k→ifor emission.i\rightarrow k \quad\text{for absorption}, \qquad k\rightarrow i \quad\text{for emission}.
ColumnMeaningRequired check
AkiA_{ki}spontaneous upper-to-lower rates−1\mathrm{s^{-1}} or units of 108 s−110^8\,\mathrm{s^{-1}}
gkAkig_kA_{ki}upper-weighted rateuse gk=2Jk+1g_k=2J_k+1 for a fine-structure level
accuracy codeassessed uncertainty categoryretain the cited transition-probability reference
typeE1, M1, E2, mixed, induced, or otherchoose the matching line-strength formula
relative intensitysource-dependent observed strengthdo not equate with AkiA_{ki}

The lifetime of upper level kk requires all radiative channels:

τrad,k−1=∑iAki.\tau_{\mathrm{rad},k}^{-1} = \sum_i A_{ki}.

A search window, wavelength table, or detector band may omit important branches. Sum over the level’s complete transition set, not merely the lines visible in one query.

Molecular AjiA_{ji} values refer to fully or partially resolved rovibronic transitions. Statistical weights can include rotational, electronic, nuclear-spin, parity, and symmetry factors. The same convention must be used in the line list and its partition function.

HITRAN, for example, stores a temperature-independent AjiA_{ji} in s−1\mathrm{s^{-1}}, while its line intensity is temperature dependent. In HITRAN’s wavenumber and cgs-style convention,

Sij(HIT)(T)=IaAji8πcν~ij 2gjexp⁡(−c2Ei/T)Q(T)×[1−exp⁡(−c2ν~ij/T)],\begin{aligned} S_{ij}^{(\mathrm{HIT})}(T) &= I_{\mathrm a} \frac{A_{ji}} {8\pi c\widetilde\nu_{ij}^{\,2}} \frac{ g_j \exp(-c_2E_i/T) }{ Q(T) } \\ &\quad\times \left[ 1-\exp(-c_2\widetilde\nu_{ij}/T) \right], \end{aligned}

where:

  • Sij(HIT)S_{ij}^{(\mathrm{HIT})} has units cm−1/(molecule cm−2)\mathrm{cm^{-1}/(molecule\,cm^{-2})}, equivalently cm molecule−1\mathrm{cm\,molecule^{-1}};
  • ν~ij\widetilde\nu_{ij} and EiE_i are in cm−1\mathrm{cm^{-1}};
  • cc is in cm s−1\mathrm{cm\,s^{-1}};
  • c2=hc/kBc_2=hc/k_{\mathrm B} is in cm K\mathrm{cm\,K};
  • gjg_j is the upper-state statistical weight;
  • Q(T)Q(T) uses the matching weight convention; and
  • IaI_{\mathrm a} is the terrestrial isotopologue abundance factor used by the database.

The database line intensity Sij(HIT)S_{ij}^{(\mathrm{HIT})} is not the reduced matrix-element line strength S21S_{21} used in the E1 formulas. The shared letter SS does not make them the same physical quantity.

The Einstein rate model is appropriate when:

  • radiation is weak enough for induced rates to remain linear in spectral energy density;
  • phase coherence between the two levels is irrelevant or rapidly lost;
  • populations within a degenerate level follow the assumed distribution;
  • memory effects and strong system–reservoir coupling are negligible;
  • line overlap can be handled by normalized profile integrals; and
  • the reservoir convention matches the coefficient definition.

Use a coherent or mode-resolved description when:

  • a narrow laser establishes a definite phase and Rabi oscillations;
  • power broadening or saturation matters;
  • short pulses prepare coherent superpositions;
  • interference among pathways is observable;
  • a cavity or nanophotonic environment reshapes the mode density; or
  • non-Markovian decay invalidates a constant AA.

A collimated laser irradiance is not the isotropic uωu_\omega of Einstein’s thermal argument. Converting irradiance into a field amplitude and using a Rabi frequency is usually the cleanest route.

BB becomes a rate only after multiplication by its matching spectral energy density or after integration over a line profile.

B(ω)B^{(\omega)}, B(ν)B^{(\nu)}, B(ν~)B^{(\widetilde\nu)}, and B(λ)B^{(\lambda)} are numerically different. The omission is especially dangerous for ω\omega versus ν\nu because their reduced SI dimensions are the same.

For level-averaged coefficients,

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

Equality requires matched degeneracies or resolved state-to-state coefficients.

The radiative inverse lifetime is ∑ℓA2ℓ\sum_\ell A_{2\ell}. Nonradiative and collisional rates further shorten an observed lifetime.

Emission also depends on upper-state population, photon energy, reabsorption, collection geometry, branching, and detector response.

For a fixed free-space transition model, AA is intrinsic. Thermal populations and stimulated processes change line intensity and total transition rates, not the stored spontaneous coefficient.

M1, E2, mixed, and induced transitions require their own microscopic strength conventions.

Applying rate equations to coherent strong driving

Section titled “Applying rate equations to coherent strong driving”

Rabi dynamics, optical coherences, saturation, and dressed states require optical Bloch or fuller quantum-optical equations.

A reproducible statement should specify:

  1. lower and upper state or level labels, including all resolved quantum numbers;
  2. g1g_1 and g2g_2 and what each degeneracy counts;
  3. whether A21A_{21} is one channel or a sum;
  4. transition frequency and vacuum, air, or medium wavelength convention;
  5. whether BB multiplies uωu_\omega, uνu_\nu, wavenumber density, wavelength density, radiance, irradiance, or another field quantity;
  6. line-profile normalization and whether the narrow-line approximation was used;
  7. transition multipole and environment;
  8. radiative and nonradiative channels used in a lifetime;
  9. uncertainty, accuracy code, and bibliographic provenance; and
  10. temperature, population, partition-function, and isotopic-abundance conventions for molecular intensities.

Exercise 1: Lifetime, branching, and quantum yield

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

An upper level has radiative channels A21=3.0×107 s−1A_{21}=3.0\times10^7\,\mathrm{s^{-1}} and A20=1.0×107 s−1A_{20}=1.0\times10^7\,\mathrm{s^{-1}}, plus a nonradiative rate Γnr=1.0×107 s−1\Gamma_{\mathrm{nr}}=1.0\times10^7\,\mathrm{s^{-1}}. Find the radiative lifetime, radiative branching fraction into level 11, observed lifetime, and radiative quantum yield.

Solution

The radiative rate is

Γrad=3.0×107+1.0×107=4.0×107 s−1.\Gamma_{\mathrm{rad}} = 3.0\times10^7+1.0\times10^7 = 4.0\times10^7\ \mathrm{s^{-1}}.

Therefore

τrad=14.0×107=25 ns.\tau_{\mathrm{rad}} = \frac{1}{4.0\times10^7} = 25\ \mathrm{ns}.

The radiative branching fraction is

b21(rad)=3.04.0=0.75.b_{21}^{(\mathrm{rad})} = \frac{3.0}{4.0} = 0.75.

Including the nonradiative channel,

Γtot=5.0×107 s−1,τobs=20 ns.\Gamma_{\mathrm{tot}} = 5.0\times10^7\ \mathrm{s^{-1}}, \qquad \tau_{\mathrm{obs}} = 20\ \mathrm{ns}.

Finally,

Φrad=4.05.0=0.80.\Phi_{\mathrm{rad}} = \frac{4.0}{5.0} = 0.80.

The single channel A21A_{21} is therefore not the inverse lifetime in this example.

A level has radiative lifetime τrad=40.0 ns\tau_{\mathrm{rad}}=40.0\,\mathrm{ns} and complete radiative branching fractions 0.700.70, 0.200.20, and 0.100.10. Find the three AA coefficients.

Solution

The total radiative rate is

Γrad=140.0×10−9 s=2.50×107 s−1.\Gamma_{\mathrm{rad}} = \frac{1}{40.0\times10^{-9}\,\mathrm s} = 2.50\times10^7\ \mathrm{s^{-1}}.

Multiplying by each branching fraction gives

A1=1.75×107 s−1,A2=5.00×106 s−1,A3=2.50×106 s−1.\begin{aligned} A_1 &= 1.75\times10^7\ \mathrm{s^{-1}}, \\ A_2 &= 5.00\times10^6\ \mathrm{s^{-1}}, \\ A_3 &= 2.50\times10^6\ \mathrm{s^{-1}}. \end{aligned}

Their sum recovers 2.50×107 s−12.50\times10^7\,\mathrm{s^{-1}}.

Exercise 3: Degeneracy and spectral convention

Section titled “Exercise 3: Degeneracy and spectral convention”

Suppose g1=2g_1=2, g2=6g_2=6, and B21(ω)=4.00×1020B_{21}^{(\omega)}=4.00\times10^{20} in the matching SI units. Find B12(ω)B_{12}^{(\omega)} and B12(ν)B_{12}^{(\nu)}.

Solution

Detailed balance gives

B12(ω)=g2g1B21(ω)=3(4.00×1020)=1.20×1021.B_{12}^{(\omega)} = \frac{g_2}{g_1} B_{21}^{(\omega)} = 3(4.00\times10^{20}) = 1.20\times10^{21}.

Changing from angular to ordinary frequency,

B12(ν)=B12(ω)2π=1.91×1020.\begin{aligned} B_{12}^{(\nu)} &= \frac{B_{12}^{(\omega)}}{2\pi} \\ &= 1.91\times10^{20}. \end{aligned}

The degeneracy factor and the 2π2\pi Jacobian answer different questions and must both be applied.

An E1 transition has λ=5000.0 A˚\lambda=5000.0\,\mathring{\mathrm A}, f12=0.200f_{12}=0.200, g1=2g_1=2, and g2=4g_2=4. Find A21A_{21}. If it is the sole decay channel, find the radiative lifetime.

Solution

Rearrange the numerical conversion:

A21=f121.49919×10−16(g2/g1)λA˚ 2.A_{21} = \frac{f_{12}} { 1.49919\times10^{-16} (g_2/g_1) \lambda_{\mathring{\mathrm A}}^{\,2} }.

Substitution gives

A21=2.668×107 s−1.A_{21} = 2.668\times10^7\ \mathrm{s^{-1}}.

For a single radiative channel,

τrad=1A21=37.48 ns.\tau_{\mathrm{rad}} = \frac{1}{A_{21}} = 37.48\ \mathrm{ns}.

If other channels exist, this last step is invalid until their AA values are included.

For a transition with A21=1.00×107 s−1A_{21}=1.00\times10^7\,\mathrm{s^{-1}}, g2/g1=3g_2/g_1=3, and thermal occupation n‾=0.20\overline n=0.20, find the upward and downward rates per system. Verify event balance using the equilibrium population ratio.

Solution

The downward rate per upper-state system is

R2→1=A21(n‾+1)=1.20×107 s−1.R_{2\rightarrow1} = A_{21}(\overline n+1) = 1.20\times10^7\ \mathrm{s^{-1}}.

The upward rate per lower-state system is

R1→2=g2g1A21n‾=6.00×106 s−1.R_{1\rightarrow2} = \frac{g_2}{g_1} A_{21}\overline n = 6.00\times10^6\ \mathrm{s^{-1}}.

The Bose factor implies

e−ℏω0/(kBT)=n‾n‾+1=16.e^{-\hbar\omega_0/(k_{\mathrm B}T)} = \frac{\overline n}{\overline n+1} = \frac{1}{6}.

Hence

N2N1=3(16)=12.\frac{N_2}{N_1} = 3\left(\frac{1}{6}\right) = \frac{1}{2}.

Therefore

N1(6.00×106)=N2(1.20×107),N_1(6.00\times10^6) = N_2(1.20\times10^7),

as required.

Estimate the thermal photon occupation at ν0=10.0 GHz\nu_0=10.0\,\mathrm{GHz} and T=300 KT=300\,\mathrm K. Is spontaneous emission the dominant downward contribution?

Solution

The dimensionless ratio is

hν0kBT=1.600×10−3.\frac{h\nu_0}{k_{\mathrm B}T} = 1.600\times10^{-3}.

Thus

n‾=1exp⁡(1.600×10−3)−1=624.6.\overline n = \frac{1} {\exp(1.600\times10^{-3})-1} = 624.6.

The stimulated downward contribution is A21n‾A_{21}\overline n, whereas the spontaneous contribution is A21A_{21}. Thermal stimulated emission is therefore about 625625 times larger than the extra spontaneous term. Room-temperature microwave experiments must account for thermal photons.

An upper level decays through channels with A1=6.0×107 s−1A_1=6.0\times10^7\,\mathrm{s^{-1}} and A2=2.0×107 s−1A_2=2.0\times10^7\,\mathrm{s^{-1}}. The lower level is stable and pure dephasing is negligible. Find the lifetime and ordinary-frequency FWHM.

Solution

The total population-decay rate is

Γ1=A1+A2=8.0×107 s−1.\Gamma_1 = A_1+A_2 = 8.0\times10^7\ \mathrm{s^{-1}}.

Hence

τ=1Γ1=12.5 ns.\tau = \frac{1}{\Gamma_1} = 12.5\ \mathrm{ns}.

The natural Lorentzian width is

ΔνFWHM=Γ12π=1.273×107 Hz=12.73 MHz.\begin{aligned} \Delta\nu_{\mathrm{FWHM}} &= \frac{\Gamma_1}{2\pi} \\ &= 1.273\times10^7\ \mathrm{Hz} \\ &= 12.73\ \mathrm{MHz}. \end{aligned}

Using only A1A_1 would omit the coherence loss caused by the second decay channel.

A report says: “The HITRAN Einstein AA coefficient doubles from 200 K to 400 K, and the stimulated-emission rate is BIBI for laser irradiance II.” Identify the two convention errors.

Solution

First, a stored free-space AjiA_{ji} for a fixed molecular transition is temperature independent. The HITRAN line intensity changes with temperature because level populations, the partition function, stimulated-emission correction, and possibly isotopic conventions enter. That does not make AjiA_{ji} temperature dependent.

Second, the Einstein BB used here multiplies a specified isotropic spectral energy density, not an unlabeled irradiance. A directional coherent laser also selects propagation direction and polarization and can preserve optical phase. One must either convert the laser spectrum and geometry to the exact field convention used to define BB, or, more transparently, use the field amplitude, transition dipole, Rabi frequency, and optical Bloch equations.