Codex @codex  0

For transverse coupling with equal mode amplitudes $\lambda$, Fermi's golden rule and the oscillator creation operator and annihilation operator givewhere the first index is the final spin and $g$ is the density of states per unit angular frequency. Decay creates a quantum and carries the factor $n_{\Omega }+1$; excitation annihilates one and carries $n_{\Omega }$. For nonuniform resonant occupations, $n_{\Omega }$ is their density-weighted mean.

With a thermal Bose-Einstein distribution, $n_{\Omega }/(n_{\Omega }+1)=e^{-\beta \hslash \Omega }$. The two-state master equation therefore has the Boltzmann distribution as its unique stationary populations when ${\Gamma }_0>0$. The population deviation relaxes at rate ${\Gamma }_0(2n_{\Omega }+1)$. If there is no resonant coupling, both rates vanish and stationarity does not determine the populations.

In weak-probe optical absorption, Fermi's golden rule creates an electron-hole pair, with negligible population depletion. Negligible photon momentum makes a direct band gap transition nearly vertical in crystal momentum:For allowed transitions between three-dimensional parabolic bands, the joint density of states yields $\alpha \propto \sqrt{\hslash \omega -E_g}/\omega$ above the gap. Selection rules and polarization enter the dipole matrix element. An indirect band gap needs an additional phonon to conserve crystal momentum, with phonon absorption/emission weights and shifted thresholds.

Electron-hole attraction produces Wannier excitons below the continuum, approximately $E_n=E_g-R^{\ast }/n^2$, with bright $s$-state oscillator strength decreasing as $n^{-3}$. It also enhances the continuum edge. Screening, disorder and phonons broaden these features; heavy/light valence bands can create distinct series. This linear-response picture includes Coulomb binding but excludes a resolved coherent photon–exciton anticrossing: that would require a strong light–matter coupling treatment. “Weak coupling” here is distinguished from neglecting the electron-hole interaction altogether.

The absorptive dielectric function counts optically accessible electron-hole pairs. In the electric dipole approximation the perturbation is proportional to $\mathbf{a}\cdot \hat{\mathbf{p}}$, with $\mathbf{a}$ the light polarization and $\hat{\mathbf{p}}$ the electronic momentum operator. The squared matrix element between occupied valence and empty conduction states is the transition strength; it includes polarization and band-symmetry selection rules. The sums enumerate initial and final subbands. photon momentum is negligible on the crystal-momentum scale, so the transition is approximately vertical in the in-plane wavevector $\mathbf{k}$.

For each subband pair, $C$ is the constant-transition-energy contour in the two-dimensional Brillouin zone. Its line element is $dl$, and the reciprocal gradient weights the states satisfying energy conservation. Specifically, the level-set representation of the joint density of states isA small energy gradient means many nearly resonant states; a vanishing gradient can produce a Van Hove singularity or an edge. The displayed expression suppresses overall optical prefactors and assumes the momentum matrix element is approximately constant along each contour. More generally it stays inside that integral. Occupation factors $f_v-f_c$ are also implicit; at low temperature in an undoped well the relevant valence states are full and the conduction states empty. This is the independent-particle continuum contribution from Fermi's golden rule. Bound excitons add optical resonances below the corresponding free-pair thresholds.

The GaAs four-band model is rooted in bonding and spin-orbit coupling. Near the zone centre, the conduction edge is predominantly s-like and antibonding, while the upper valence manifold is predominantly p-like. Combining orbital angular momentum $l=1$ with spin $s=1/2$ splits the six p-derived states into a $j=3/2$ quartet and a lower $j=1/2$ doublet. The latter is the split-off valence band. Within the upper quartet, the $m_j=\pm 3/2$ and $m_j=\pm 1/2$ branches, defined relative to a propagation or quantization direction, have different curvatures and are called heavy holes and light holes. They coincide at the unstrained bulk zone centre but need not remain degenerate in a well. The conduction edge is a spin doublet. Thus “four bands” counts conduction, heavy-hole, light-hole and split-off branches, including their spin multiplicities; it does not count four individual states. Positive hole masses describe the downward valence curvature, whereas the conduction curvature is upward.

For the final exciton comparison, define the reduced mass ${\mu }_j$ and screened exciton Rydberg $R_j^{\ast }$ byThe resulting $({\mu }_j/m_e,R_j^{\ast })$ are approximately $(0.05948,4.937\mathrm{meV})$ for heavy holes, $(0.03646,3.027\mathrm{meV})$ for light holes and $(0.04631,3.844\mathrm{meV})$ for split-off holes. In the ideal two-dimensional exciton model, the first two optically active internal states lie $4R_j^{\ast }$ and $4R_j^{\ast }/9$ below their continuum edge. Their separation is $32R_j^{\ast }/9$, independent of a common threshold offset. The corresponding three-dimensional separation would be only $3R_j^{\ast }/4$. The listed peak spacings specified for this theoretical spectrum giveHere each pair comprises the stronger lower-energy peak and its weaker nearby companion, in the order associated with the continuum edges calculated in part (c). The approximately 18-meV and 11-meV pair splittings support the two-dimensional series and are incompatible with the three-dimensional series at the supplied masses and permittivity. Equivalently, for a two-dimensional pair with energies $E_1,E_2$, the inferred edge is $E_2+(E_2-E_1)/8$; the observed pair splittings predict ground-state bindings around 20 and 12 meV, close to $4R_j^{\ast }$.

Absolute agreement should not be overstated. For example, the stated square-well parameters give the first heavy-hole edge $1.56390$ eV and the first light-hole edge $1.59161$ eV. Their ideal two-dimensional peak predictions are approximately $(1.54415,1.56170)$ eV and $(1.57950,1.59027)$ eV, respectively. The listed peaks are lower by about 1–2 meV across the families; rounded parameters and this simple effective-mass model do not reproduce every absolute energy exactly. No gap, width or mass has been silently adjusted to fit them. The peak-spacing comparison remains useful because it avoids a common edge offset. The split-off transitions predicted by the four-band data are discussed separately in part (c), rather than assigned fictitious extra peaks.

At 300 K, $k_{B}T\simeq 25.9$ meV, whereas the ideal ground-state binding energies are only about $19.7$, $12.1$ and $15.4$ meV; excited states are still more weakly bound. None of these series is expected to remain a clearly resolved, robust bound-exciton peak at room temperature in this simple estimate. Thermal ionization and phonon broadening work against their visibility. This is a scale argument, not a theorem that the exciton population or all Coulomb-enhanced absorption becomes exactly zero: actual observability also requires temperature-dependent gaps, linewidths, screening and carrier density, which are not supplied.

Fermi's golden rule converts a weak transition matrix element and the available density of final states into a transition rate. For a static perturbation $W$ and a continuum,When the matrix element is nearly constant on the energy shell, this becomes $(2\pi /\hslash )|W_{fi}{|}^2\rho (E_i)$. The energy delta function is to be integrated, rather than interpreted as a finite rate to one exactly resonant discrete level. The observation time must resolve a smooth final-state continuum while remaining short enough that first-order depletion is small. A periodically driven perturbation instead gives the corresponding absorption and emission energy shells.

For elastic scattering from a weak localized potential $W(\mathbf{r})$, use box-normalized free states $⟨\mathbf{r}|\mathbf{k}⟩={\mathcal{V}}^{-1/2}{e}^{i\mathbf{k}\cdot \mathbf{r}}$. ThenThe number of final states is $\mathcal{V}{k}_f^2\mathrm{d}{k}_f\mathrm{d}\Omega /(2\pi {)}^3$. With $E_f={\hslash }^2{k}_f^2/(2m)$, the density per energy and solid angle is $\mathcal{V}m{k}_f/[(2\pi {)}^3{\hslash }^2]$. Dividing the resulting golden-rule rate by the incident single-particle flux $j_i=\hslash k_i/(m\mathcal{V})$ givesThis is the golden-rule derivation of the Born cross-section. For the elastic problem $k_f=k_i$, so the first Born scattering amplitude and cross-section areThe box volume cancels, as a physical cross-section must. The minus sign is fixed by the outgoing-wave scattering convention; the golden-rule rate alone determines only the modulus.

For example, $W(r)=g{e}^{-\kappa r}/r$ has $\tilde{W}=4\pi g/(q^2+{\kappa }^2)$, hence $f_B=-2mg/[{\hslash }^2(q^2+{\kappa }^2)]$, with $q=2k\sin (\theta /2)$. Replacing the exact scattering state by the incident plane wave is the essential Born approximation. It can fail for strong scattering, low-energy resonances and situations requiring repeated scattering; a small perturbation and a suitable continuum are separate conditions.^[ref]

The Pauli matrix ${\sigma }_x$ flips the spin with unit matrix element. A nonzero oscillator matrix element changes exactly one occupation number, through a creation operator or an annihilation operator:For a spin flip, the spin energy change divided by $\hslash$ is $-s_{\sigma }\Omega$. Thus creation and annihilation in mode $j$ have total frequency changes ${\omega }_j-s_{\sigma }\Omega$ and $-{\omega }_j-s_{\sigma }\Omega$, respectively. Distinct modes or different occupation changes lead to orthogonal final quantum states, so their probabilities add. Fermi's golden rule consequently givesThe Dirac delta function is even, accounting for the sign in the annihilation term. Since all ${\omega }_j,\Omega >0$, an initially up spin can only emit an oscillator quantum, while an initially down spin can only absorb one. Writing the oscillator density of states per unit angular frequency as $g(\omega )={\sum }_j\delta (\omega -{\omega }_j)$, the continuum limit givesThe $+1$ permits emission even from the oscillator vacuum; absorption requires a quantum already present. The notation $n_{\Omega }$ assumes a common or smoothly specified occupation at the resonant frequency. If several resonant modes have different occupations, it denotes their density of states weighted mean, so that $g(\Omega )n_{\Omega }={\sum }_j{n}_j\delta (\Omega -{\omega }_j)$ after continuum smoothing. This is the golden-rule relaxation of a spin coupled to an oscillator bath.

In first-order time-dependent perturbation theory, the resonant time factor before taking the Fermi's golden rule limit isIts central frequency width is of order $1/t$, and its total area is $2\pi t$. To replace it by $2\pi t\delta (\nu )$, the oscillator density of states and relevant occupations must be approximately constant over that width. Thus $t\Delta \omega \gg 1$. At the same time, the perturbative flip probability must remain much smaller than one.

Define ${\Gamma }_0=2\pi {\lambda }^{2}g(\Omega )/{\hslash }^2$, so the two flip rates are ${\Gamma }_0(n_{\Omega }+1)$ and ${\Gamma }_0{n}_{\Omega }$. A sufficient common time window for both initial spins isFor a specified initial spin, the upper condition is simply its own flip rate times $t\ll 1$. Such a window requires weak enough coupling. The continuum approximation also presumes times shorter than recurrences from a resolved finite oscillator spectrum, and long enough to resolve the positive-frequency resonance from its nonresonant counterpart.

Separate the entire-system atom index into a cell translation $\mathbf{R}$ and an index $p$ within a single primitive unit cell. The potential becomesThe same $u_p$ is used in every translated cell because the atomic basis repeats. Normalize plane waves in volume $\Omega$ by $⟨\mathbf{x}|\mathbf{k}⟩={\Omega }^{-1/2}{e}^{i\mathbf{k}\cdot \mathbf{x}}$. For $\mathbf{G}={\mathbf{k}}^{\prime }-\mathbf{k}$, substitution and $\mathbf{s}=\mathbf{x}-\mathbf{R}-{\mathbf{r}}_p$ giveFor a real static scalar potential, define the atomic scattering factor in the positive-exponent convention of the question byThe second bracket in the matrix element is $S(\mathbf{G}{)}^{\ast }$. At a reciprocal lattice vector every term of the first bracket is one, giving the number of cells $N_c$. Therefore the stated transition-rate rule, equivalently the first Born approximation and Fermi's golden rule up to the final-state kinematic factors, yieldsThe harmless complex conjugation records the distinction between the source's positive-exponent Fourier transform and the usual positive-exponent incoming plane wave. Reversing every Fourier sign gives the equivalent amplitude convention. Only the cell-translation sum ranges over the full crystal; the crystal basis structure factor contains the atoms of one cell.

In Fermi's golden rule,$\Gamma$ is a transition decay rate, measured in inverse time; $\hslash$ is the reduced Planck constant; $M=⟨f|H_{\mathrm{int}}|i⟩$ is the transition matrix element of the interaction part of the Hamiltonian; and $\rho$ is the final density of states per unit energy, evaluated subject to conservation of energy. Final spin and other discrete labels must be summed, with initial labels averaged when the preparation requires it. If $M$ varies across the final phase space, $|M{|}^2\rho$ denotes the corresponding integral rather than a product of constants.

A convention warning: many quantum field theory texts use $\Gamma$ for an energy decay width. That convention is ${\Gamma }_{\mathrm{width}}=\hslash {\Gamma }_{\mathrm{rate}}$; the formula as printed uses the rate convention.

Feynman diagrams organize the perturbative calculation of the matrix element in Fermi's golden rule. Each external line labels an incoming or outgoing particle state; each internal line contributes a propagator for a generally off-shell intermediate field; and each interaction vertex contributes a coupling and the appropriate spin, momentum and internal-symmetry structure. All indistinguishable allowed diagrams must be summed as probability amplitudes before taking the absolute square, so their interference contributes to the rate.

In the radiative pion decay diagram in part (b), solid arrowed lines are fermions; the wavy line joining the hadronic and leptonic vertices is a W boson propagator; and the wavy outgoing line is an external photon. Arrows track fermion number, so an outgoing antimuon has an arrow towards its production vertex. Each vertex conserves four-momentum; an internal propagator need not satisfy the on-shell energy–momentum relation. A closed loop Feynman diagram additionally requires an integration over an internal four-momentum.

The final Lorentz-invariant phase space, conservation Dirac delta functions, and the initial-state normalization or incident flux turn the squared scattering amplitude into a decay rate or cross section. A diagram supplies the dynamics; it is not itself the final density of states.

At energies far below the $W$ mass, replace its propagator by a constant. The modern V-A interaction form of Fermi's contact theory has, in $\hslash =c=1$ units,Nucleons are treated nonrelativistically within nuclear states. To leading order in their velocity, the vector charge and axial spatial current give the two nuclear operatorswhere ${\tau }^{+}$ here means neutron-to-proton conversion. (i) The scalar Fermi operator permits $\Delta J=0$, including $0\to 0$, without parity change. The rank-one Gamow–Teller operator permits $\Delta J=0,1$ without parity change, excluding $0\to 0$. Both terms can contribute when the same nuclear states permit them; their strengths contain nuclear matrix elements, not just phase space.

(ii) Expand the lepton waves across the nucleus in powers of their momenta times nuclear radius. When $pR/\hslash \ll 1$, the constant spatial term gives allowed transitions. If its matrix element vanishes by spin/parity selection, higher terms supply lepton orbital angular momentum: “forbidden” means suppressed, not impossible. First-forbidden terms change parity and can have $\Delta J=0,1,2$; the rank-two $\Delta J=2$ case is unique first-forbidden. Higher orders carry additional small spatial factors and different energy-dependent shape factors. The contact approximation alone does not imply the allowed approximation.

For a single branch, neglect recoil and neutrino mass and first assume an energy-independent nuclear shape factor. Fermi's golden rule integrates the neutrino phase space with $E_{\nu }=Q-T$, while the electron momentum volume is proportional to $p^{2}dp$. Including Coulomb distortion,This momentum-bin Kurie ordinate is linear in kinetic energy and vanishes at its endpoint. If instead using counts per unit kinetic energy, $dN/dT\propto F(p)p{E}_e(Q-T{)}^2$ in $c=1$ units; converting with $dT/dp=p/E_e$ recovers the supplied momentum-bin definition. No extra factor of $E_e$ should be left inside that definition of $N^{\prime }$.

The graph's horizontal coordinate is $T/(m_e{c}^2)$, not $T$ in MeV. Its intercept is approximately $3.9$–$4.0$, soFor a ground-to-ground branch between bare nuclei, the emitted electron's rest energy must also be supplied:For neutral atomic masses the mass difference instead equals approximately $Q\simeq 2.0\mathrm{MeV}$, neglecting small electron-binding corrections. The graph supports this approximate precision, not a precise evaluated isotope mass. Interpreting a branch to an excited daughter would require adding its excitation energy; the ground-state comparison assumes the ground-state branch intended here.

The shell-model prediction for yttrium-91 beta decay follows from shell filling: the proton orbitals above 28 fill $2p_{3/2}$ with four protons and $1f_{5/2}$ with six, reaching 38. Yttrium has its 39th proton in $2p_{1/2}$, giving $J_i^{\pi }=1/2^{-}$. Its 52 neutrons include the closed 50-neutron core and a pair in $2d_{5/2}$, coupled to zero in the simple paired ground state. Zirconium has 40 paired protons and 51 neutrons, leaving an unpaired $2d_{5/2}$ neutron, so $J_f^{\pi }=5/2^{+}$. ThereforeThe single-particle change is neutron $d_{5/2}$ to proton $p_{1/2}$; the parity change and rank two prevent an allowed Fermi or Gamow–Teller contribution.

The straight plotted line is a constant-shape idealization. For the shell-model transition just inferred, a leading unique-forbidden shape factor contains $C(T)\propto p_e^2+p_{\nu }^2$ in $c=1$ units, with Coulomb/finite-size refinements. Consequently the uncorrected quantity specified in the paper is generally proportional to $(Q-T)\sqrt{C(T)}$, not exactly linear. A consistent shape-corrected Kurie plot also divides its ordinate by $\sqrt{C(T)}$. The source's exact straight-line construction and its forbidden-transition inference cannot both be exact with a Coulomb-only $F$; the approximate graphical endpoint can still be read as requested. Allowed and unique-forbidden spectral shapes are experimentally distinguishable.^[ref]

Beta decay converts a neutron to a proton or vice versa through the weak interaction, leaving the mass number fixed:Inside a nucleus the full initial and final nuclear masses determine whether either process is possible. Electron capture, $p+e^{-}\to n+{\nu }_e$, is the related bound-electron process. The emitted charged-lepton spectrum is continuous because energy is shared with the neutrino and nuclear recoil.

For allowed beta decay, the lepton wavefunctions' leading spatial terms are constant over the nucleus. A Fermi transition has $\Delta J=0$ and no parity change. A Gamow-Teller transition has $\Delta J=0,1$ and no parity change, excluding $0\to 0$. Thus a same-parity transition with nonzero initial and final spins and $\Delta J=0$ can contain both contributions. First-forbidden transitions change parity and permit $|\Delta J|=0,1,2$; the $|\Delta J|=2$ case is unique forbidden beta decay. Higher forbiddenness requires higher powers of the small lepton momentum times nuclear size. “Forbidden” denotes suppression of the leading allowed operator, not an impossible decay.^[ref]

Fermi's golden rule expresses a rate as a squared nuclear matrix element times leptonic phase space. Neglecting charged-lepton mass and Coulomb distortion, a constant allowed matrix element givesThis is the approximate Sargent rule, $\lambda \propto Q^5$ for comparable allowed strengths; it is not an exact law across arbitrary transitions. Finite electron mass, the Fermi function, forbidden shape factors and nuclear structure change the comparison. A microscopic first-forbidden unique operator introduces additional momentum factors.

Using bare nuclear masses $M_i,M_f$, the kinetic-energy release isThe daughter is different in the two cases: its charge is $Z+1$ for beta minus and $Z-1$ for beta plus. The relevant mass difference must exceed one electron mass. For neutral atomic masses, to leading order in small electronic binding corrections, the corresponding expressions are $Q_{{\beta }^{-}}=(M_i^{\mathrm{at}}-M_f^{\mathrm{at}})c^2$ and $Q_{{\beta }^{+}}=(M_i^{\mathrm{at}}-M_f^{\mathrm{at}}-2m_e)c^2$. Electron capture is energetically possible approximately when $M_i^{\mathrm{at}}>M_f^{\mathrm{at}}$, with available shell binding energy and daughter excitation included. Nuclear masses and atomic masses must not be mixed.

Assume allowed beta decay, so the leading matrix element has negligible lepton-momentum dependence. Neglect nuclear recoil and neutrino mass. Let $Q$ be the kinetic energy available to the electron and antineutrino, $E_0=m_e+Q$, and $E=\sqrt{p^2+m_e^2}$. The neutrino energy and momentum magnitude are $k=E_0-E$.

Fermi's golden rule weights accessible final states. After angular integration, the momentum-space factors are $p^{2}dp$ and $k^{2}dk$. The energy-conserving Dirac delta function removes the neutrino integral:The constant includes the allowed nuclear matrix element, weak-coupling and angular factors. Including the usual Fermi function for electron-nucleus Coulomb distortion givesFor the kinetic-energy spectrum, $dp/dT=E/p$ instead yieldsThe two distributions differ by a Jacobian; forgetting it changes the shape.

A pure fifth-power law additionally neglects electron mass and Coulomb distortion and keeps the nuclear matrix element fixed. Then $E=p$, $E_0\simeq Q$, andTwo momentum powers come from each final particle and one from the remaining integration. This is the origin of the Sargent rule. Exact $Q^5$ dependence requires the relativistic-electron, constant-matrix-element limit. For $Q\ll m_e$, the same approximation without Coulomb distortion has $pE\simeq m_e\sqrt{2m_{e}T}$ and instead scales as $m_e^{3/2}{Q}^{7/2}$. Forbidden-transition factors, nuclear-structure changes and Coulomb distortion also modify the simple law.

In the simple box-state counting used here, the final three-momentum lies in a spherical shell: $dn=C{p}^{2}dp$, where the common constant $C$ contains the normalization volume and angular factor. Momentum conservation fixes the other final particle's momentum. Hence the density of states isCombining parts (a) and (b) givesWith the paper's rounded masses,By Fermi's golden rule, equal squared matrix elements make the rate ratio equal to this density ratio. Thus the naive constant-matrix-element estimate is $R\simeq 0.287$, which favours the electronic decay by about $3.49$ to one.

This is the particular noncovariant state-counting convention in the question. A calculation with the Lorentz-invariant phase-space measure places factors $1/(2E)$ in the measure and correspondingly normalizes its decay amplitude differently. Only the complete product of phase space and squared matrix element is convention-independent; mixing the two conventions would give a wrong rate.

With relativistically normalized one-particle states, Fermi's golden rule becomesIn the parent rest frame, the three-momentum delta function fixes ${\mathbf{p}}_2=-{\mathbf{p}}_1$. The energy delta function fixes the common magnitude $p^{\ast }$; its radial Jacobian converts the remaining phase-space factor intoOnly one momentum magnitude remains because the two final momenta are constrained, not independent. Explicitly,For leptonic W boson decay both final masses are small compared with $m_W$, including the tau mass at the precision requested, so $p^{\ast }\simeq m_W/2$.

The Time-dependent Schrödinger equation $i\hslash {\partial }_t|\psi ⟩=H(t)|\psi ⟩$ governs unitary evolution of a closed system. If $H$ is independent of time, expand in its energy eigenstates:An individual energy eigenstate changes only by a phase, so its probability density and time-independent-observable expectations are stationary. Superpositions can have interference terms oscillating at $(E_n-E_m)/\hslash$; fixed energy populations do not imply that every observable is constant.

For general $H(t)$, $U(t,t_0)=\mathcal{T}\exp [-(i/\hslash )\int H(t^{\prime })d{t}^{\prime }]$ requires time ordering unless Quantum Hamiltonians at different times commute. The equivalent Heisenberg description evolves observables rather than the state, giving $d⟨A⟩/dt=(i/\hslash )⟨[H,A]⟩+⟨{\partial }_{t}A⟩$. Explicit time dependence can transfer energy from an applied perturbation.

For a weak perturbation $V(t)$ to a known $H_0$, the first-order transition amplitude from $n$ to $m\ne n$ isIts squared modulus gives the leading transition probability while depletion remains small. A nearly monochromatic perturbation emphasizes energy differences near its quantum of frequency; a continuum and suitable long-time limit yield Fermi's golden rule. These approximations have stated weak-coupling and time-resolution limits; they do not replace exact unitary dynamics in every regime.

For a normal transition dipole of magnitude $D$, the Fermi's golden rule, one-polarization mode density $\omega /(2\pi c^{2}d)$ and single-photon squared field $\hslash \omega /(2{ϵ}_{0}V)$ give $A_{\perp }=D^2{\omega }_0^2/(2{ϵ}_0\hslash c^{2}d)$. Dividing by the free-space electric-dipole spontaneous emission rate gives the title. The factor three relative to the scalar density-of-states ratio comes from the isotropic free-space polarization average. This ideal weak coupling formula assumes wide plates and unperturbed atomic structure, with atomic size much smaller than $d$.