Skip to content

Density of States in Transition Rates

In a transition rate, the density of states is the measure of how many distinct final states are available near the energy selected by the drive or interaction. It converts a probability for one normalized final state into a rate summed over a continuum. The result is trustworthy only when the state normalization, channel labels, degeneracy factors, and energy-shell Jacobian are all compatible with the matrix element.

This page owns that transition-rate bookkeeping. The elementary box-counting derivation lives in Density of States: First Encounter. The resolvent identity, local density of states, and spectral broadening live in Green Functions and Density of States. Compact state-counting formulas live in the Density of States formula card.

Let ∣i⟩\lvert i\rangle be the prepared state and let ∣α⟩\lvert\alpha\rangle label exact final states with energies EαE_\alpha. A standard golden-rule expression is

Γi=2πℏ∑α∣Wαi∣2δ(Eα−E∗),\Gamma_i = \frac{2\pi}{\hbar} \sum_\alpha \lvert W_{\alpha i}\rvert^2 \delta(E_\alpha-E_*),

where

Wαi=⟨α∣W∣i⟩.W_{\alpha i} = \langle\alpha\rvert W\lvert i\rangle.

The selected energy E∗E_* depends on the process. For a static perturbation,

E∗=Ei.E_*=E_i.

For absorption from a monochromatic drive,

E∗=Ei+ℏω,E_*=E_i+\hbar\omega,

while stimulated emission selects Ei−ℏωE_i-\hbar\omega. The delta function is the long-time energy-selection limit; at finite time it is replaced by a kernel of width of order ℏ/T\hbar/T.

The exact density of final states is the distribution

ρ(E)=∑αδ(E−Eα).\rho(E) = \sum_\alpha \delta(E-E_\alpha).

It has units of inverse energy when it counts a dimensionless number of states. A density per volume has units

[ρ/V]=(energy×volume)−1.[\rho/\mathcal V] = (\text{energy}\times\text{volume})^{-1}.

Density of states is not an occupation probability. It says what states exist, not which of them are populated.

The Weighted Density Is the Physical Object

Section titled “The Weighted Density Is the Physical Object”

The rate does not usually depend on ρ(E)\rho(E) alone. Define the matrix-element-weighted density

Si(E)≡∑α∣Wαi∣2δ(E−Eα).\mathcal S_i(E) \equiv \sum_\alpha \lvert W_{\alpha i}\rvert^2 \delta(E-E_\alpha).

Then

Γi=2πℏSi(E∗).\Gamma_i = \frac{2\pi}{\hbar} \mathcal S_i(E_*).

Only when the matrix element varies slowly among the states in the relevant energy window may one factorize

Si(E)≃∣W(E)∣2‾ρ(E).\mathcal S_i(E) \simeq \overline{ \lvert W(E)\rvert^2 } \rho(E).

This distinction prevents several common errors:

  • a large density of states produces no transition if symmetry makes every relevant matrix element vanish;
  • two channels with the same density of states can have very different rates;
  • a threshold singularity in ρ(E)\rho(E) can be softened or removed by a vanishing matrix element;
  • pulling ∣W∣2\lvert W\rvert^2 outside an integral is an approximation, not an identity.

The dimensions also check the formula. Since

[Si]=energy,[\mathcal S_i] = \text{energy},

the factor 2π/ℏ2\pi/\hbar converts it to inverse time.

Suppose final states are organized into channels cc and continuous coordinates q\mathbf q. A channel may encode spin, polarization, angular momentum, particle species, propagation direction, band index, or any other discrete label. Write the completeness measure schematically as

∑α⟶∑c∫dμc(q).\sum_\alpha \longrightarrow \sum_c \int d\mu_c(\mathbf q).

The channel density of states is

ρc(E)=∫dμc(q) δ(E−Ec(q)),\rho_c(E) = \int d\mu_c(\mathbf q)\, \delta\big(E-E_c(\mathbf q)\big),

and the total density is

ρ(E)=∑cρc(E).\rho(E) = \sum_c\rho_c(E).

The rate is

Γi=2πℏ∑c∫dμc(q) ×∣Wci(q)∣2δ(Ec(q)−E∗).\begin{aligned} \Gamma_i &= \frac{2\pi}{\hbar} \sum_c \int d\mu_c(\mathbf q)\, \\ &\quad\times \lvert W_{ci}(\mathbf q)\rvert^2 \delta\big(E_c(\mathbf q)-E_*\big). \end{aligned}

This form is safer than writing ∣W∣2ρ\lvert W\rvert^2\rho prematurely because the integration measure and matrix element remain visibly paired.

For coordinates q∈Rn\mathbf q\in\mathbb R^n with

dμc(q)=wc(q)dnq,d\mu_c(\mathbf q) = w_c(\mathbf q)d^nq,

the delta function restricts the integral to the constant-energy surface

Σc,E={q:Ec(q)=E}.\Sigma_{c,E} = \left\{ \mathbf q: E_c(\mathbf q)=E \right\}.

Away from critical points where ∇qEc=0\nabla_{\mathbf q}E_c=0, the coarea formula gives

ρc(E)=∫Σc,Ewc(q)dΣ∥∇qEc∥.\rho_c(E) = \int_{\Sigma_{c,E}} \frac{ w_c(\mathbf q)d\Sigma }{ \lVert\nabla_{\mathbf q}E_c\rVert }.

The factor

1∥∇qEc∥\frac{1}{ \lVert\nabla_{\mathbf q}E_c\rVert }

is the Jacobian that converts a coordinate-space shell thickness into an energy interval. Large constant-energy surfaces increase the density of states; a small energy gradient also increases it because many states fit into a narrow energy window.

At a band extremum, saddle point, flat band, or threshold, the gradient can vanish. The resulting density-of-states singularity may be physical, integrable, regularized by finite size or broadening, or modified by interactions. One should return to the original integral rather than substitute a formally divergent shell expression without analysis.

Continuum rates are often regulated by a periodic box of volume V=Ld\mathcal V=L^d. The normalized plane waves are

⟨x∣k⟩V=1Veik⋅x,\langle\mathbf x\vert\mathbf k\rangle_{\mathcal V} = \frac{1}{\sqrt{\mathcal V}} e^{i\mathbf k\cdot\mathbf x},

with

k=2πLn,n∈Zd.\mathbf k = \frac{2\pi}{L} \mathbf n, \qquad \mathbf n\in\mathbb Z^d.

The large-box sum becomes

∑k⟶V(2π)d∫ddk.\sum_{\mathbf k} \longrightarrow \frac{\mathcal V}{(2\pi)^d} \int d^dk.

For a localized initial state and localized interaction WW, a typical box-normalized matrix element scales as

Wki(V)=1VMi(k),W_{\mathbf k i}^{(\mathcal V)} = \frac{1}{\sqrt{\mathcal V}} \mathcal M_i(\mathbf k),

where

Mi(k)=∫ddx e−ik⋅xW(x)ψi(x)\mathcal M_i(\mathbf k) = \int d^dx\, e^{-i\mathbf k\cdot\mathbf x} W(\mathbf x)\psi_i(\mathbf x)

in this convention. Consequently,

∣Wki(V)∣2=1V∣Mi(k)∣2.\left\lvert W_{\mathbf k i}^{(\mathcal V)} \right\rvert^2 = \frac{1}{\mathcal V} \lvert\mathcal M_i(\mathbf k)\rvert^2.

Substitution into the golden-rule sum gives

Γi=2πℏV(2π)d∫ddk ∣Mi(k)∣2V×δ(E(k)−E∗)=2πℏ∫ddk(2π)d ∣Mi(k)∣2δ(E(k)−E∗).\begin{aligned} \Gamma_i &= \frac{2\pi}{\hbar} \frac{\mathcal V}{(2\pi)^d} \int d^dk\, \frac{ \lvert\mathcal M_i(\mathbf k)\rvert^2 }{\mathcal V} \\ &\qquad\times \delta\big(E(\mathbf k)-E_*\big) \\ &= \frac{2\pi}{\hbar} \int\frac{d^dk}{(2\pi)^d}\, \lvert\mathcal M_i(\mathbf k)\rvert^2 \delta\big(E(\mathbf k)-E_*\big). \end{aligned}

The artificial box volume cancels. If a purported rate for one localized emitter grows with the normalization volume, the continuum normalization has almost certainly been mixed inconsistently.

Cancellation of box volume between continuum state counting and normalized matrix elements

Box normalization makes the number of continuum states grow as V\mathcal V, while each squared matrix element into a plane-wave state scales as V−1\mathcal V^{-1}. Their product has a finite continuum limit. The cancellation must occur before interpreting the result as a physical rate.

One may bypass the box by choosing

⟨x∣k⟩=eik⋅x,\langle\mathbf x\vert\mathbf k\rangle = e^{i\mathbf k\cdot\mathbf x},

so that

⟨k∣k′⟩=(2π)dδ(d)(k−k′)\langle\mathbf k\vert\mathbf k'\rangle = (2\pi)^d \delta^{(d)}(\mathbf k-\mathbf k')

and

I=∫ddk(2π)d∣k⟩⟨k∣.I = \int\frac{d^dk}{(2\pi)^d} \lvert\mathbf k\rangle \langle\mathbf k\rvert.

The continuum rate then has the volume-free form obtained above. Other conventions place factors of (2π)d/2(2\pi)^{d/2} in the wavefunction and use δ(d)(k−k′)\delta^{(d)}(\mathbf k-\mathbf k') without the prefactor. Either convention is valid if completeness, matrix elements, and density-of-states measures are transformed together.

Sometimes the density-of-states Jacobian is absorbed into the state normalization. Let qq be a one-dimensional continuum coordinate within channel cc, with

⟨q,c∣q′,c′⟩=δcc′δ(q−q′).\langle q,c\vert q',c'\rangle = \delta_{cc'}\delta(q-q').

Assume Ec(q)E_c(q) is monotonic on that branch. Define

∣E,c⟩=∣q(E),c⟩∣dEc/dq∣.\lvert E,c\rangle = \frac{ \lvert q(E),c\rangle }{ \sqrt{ \left\lvert dE_c/dq\right\rvert } }.

Then

⟨E,c∣E′,c′⟩=δcc′δ(E−E′)\langle E,c\vert E',c'\rangle = \delta_{cc'}\delta(E-E')

and

Icont=∑c∫dE ∣E,c⟩⟨E,c∣.I_{\mathrm{cont}} = \sum_c \int dE\, \lvert E,c\rangle \langle E,c\rvert.

In this convention, the golden-rule rate is

Γi=2πℏ∑c∣⟨E∗,c∣W∣i⟩∣2.\Gamma_i = \frac{2\pi}{\hbar} \sum_c \left\lvert \langle E_*,c\rvert W\lvert i\rangle \right\rvert^2.

No explicit ρ(E∗)\rho(E_*) appears because the Jacobian is already inside ∣E,c⟩\lvert E,c\rangle. The energy-normalized matrix element has different units from a box-normalized matrix element. Inserting an additional density of states would double count the same Jacobian.

If a given energy has several roots qj(E)q_j(E), each monotonic branch must be a separate channel or root label. In one dimension, for example, right- and left-moving states at the same energy account for the familiar factor of two unless the physical boundary condition selects only one direction.

ConventionInner productCompleteness measureWhere the density factor lives
Finite boxδkk′\delta_{\mathbf k\mathbf k'}∑k\sum_{\mathbf k}in the increasing number of discrete states
Momentum delta(2π)dδ(d)(k−k′)(2\pi)^d\delta^{(d)}(\mathbf k-\mathbf k')∫ddk/(2π)d\int d^dk/(2\pi)^din the momentum measure and energy-shell Jacobian
Energy deltaδcc′δ(E−E′)\delta_{cc'}\delta(E-E')∑c∫dE\sum_c\int dEabsorbed into the definition of ∣E,c⟩\lvert E,c\rangle

These are different coordinate descriptions of the same continuum. A calculation may change convention, but it must transform the states, completeness relation, and matrix element as a single package.

Free-Particle Examples in One, Two, and Three Dimensions

Section titled “Free-Particle Examples in One, Two, and Three Dimensions”

For a spinless nonrelativistic particle with threshold energy EthE_{\mathrm{th}}, define

ε≡E−Eth,ε>0,\varepsilon \equiv E-E_{\mathrm{th}}, \qquad \varepsilon\gt0,

and dispersion

ε=ℏ2k22m.\varepsilon = \frac{\hbar^2k^2}{2m}.

Per box volume Vd\mathcal V_d, the density of states in dd spatial dimensions is

ρd(ε)Vd=gintSd−1(2π)dmℏ2(2mεℏ)d−2,\frac{\rho_d(\varepsilon)}{\mathcal V_d} = g_{\mathrm{int}} \frac{S_{d-1}}{(2\pi)^d} \frac{m}{\hbar^2} \left( \frac{\sqrt{2m\varepsilon}}{\hbar} \right)^{d-2},

where gintg_{\mathrm{int}} counts included internal degeneracies and

Sd−1=2πd/2Γ(d/2)S_{d-1} = \frac{2\pi^{d/2}}{ \Gamma(d/2) }

is the area of the unit sphere in dd dimensions. A factor Θ(ε)\Theta(\varepsilon) is understood.

The three most common cases are

dρd(ε)/Vd1gintπℏm2ε2gintm2πℏ23gint4π2(2mℏ2)3/2ε\begin{array}{c|c} d & \rho_d(\varepsilon)/\mathcal V_d \\ \hline 1 & \displaystyle \frac{g_{\mathrm{int}}}{\pi\hbar} \sqrt{\frac{m}{2\varepsilon}} \\[0.9em] 2 & \displaystyle \frac{g_{\mathrm{int}}m}{2\pi\hbar^2} \\[0.9em] 3 & \displaystyle \frac{g_{\mathrm{int}}}{4\pi^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \sqrt{\varepsilon} \end{array}

Thus

ρd(ε)∝εd/2−1.\rho_d(\varepsilon) \propto \varepsilon^{d/2-1}.
  • In one dimension the state density diverges as ε−1/2\varepsilon^{-1/2} at threshold.
  • In two dimensions it approaches a constant step.
  • In three dimensions it vanishes as ε\sqrt{\varepsilon}.

These are density-of-states statements, not complete rate laws. The matrix element, channel boundary conditions, and any incident-flux factor can introduce additional powers of momentum.

Worked Energy-Shell Integral in Three Dimensions

Section titled “Worked Energy-Shell Integral in Three Dimensions”

Consider the volume-independent continuum expression

Γi=2πℏ∫d3k(2π)3 ×∣Mi(k)∣2δ(ℏ2k22m−ε∗).\begin{aligned} \Gamma_i &= \frac{2\pi}{\hbar} \int\frac{d^3k}{(2\pi)^3}\, \\ &\quad\times \lvert\mathcal M_i(\mathbf k)\rvert^2 \delta\left( \frac{\hbar^2k^2}{2m}-\varepsilon_* \right). \end{aligned}

Let

k∗=2mε∗ℏ.k_* = \frac{\sqrt{2m\varepsilon_*}}{\hbar}.

For k≥0k\ge0,

δ(ℏ2k22m−ε∗)=mℏ2k∗δ(k−k∗).\delta\left( \frac{\hbar^2k^2}{2m}-\varepsilon_* \right) = \frac{m}{\hbar^2k_*} \delta(k-k_*).

Using d3k=k2dk dΩkd^3k=k^2dk\,d\Omega_{\mathbf k} gives

Γi=2πℏmk∗(2π)3ℏ2∫dΩk ∣Mi(k∗,Ωk)∣2.\Gamma_i = \frac{2\pi}{\hbar} \frac{mk_*}{(2\pi)^3\hbar^2} \int d\Omega_{\mathbf k}\, \left\lvert \mathcal M_i(k_*,\Omega_{\mathbf k}) \right\rvert^2.

If the on-shell matrix element is isotropic,

Γi=2πℏ∣Mi(k∗)∣2mk∗2π2ℏ2.\Gamma_i = \frac{2\pi}{\hbar} \lvert\mathcal M_i(k_*)\rvert^2 \frac{mk_*}{2\pi^2\hbar^2}.

The last factor is precisely the spinless three-dimensional density of states per volume in the momentum-delta convention appropriate to Mi\mathcal M_i.

If the on-shell matrix element approaches a nonzero constant, the free-particle threshold scaling follows the density of states:

Γ(ε)∝εd/2−1.\Gamma(\varepsilon) \propto \varepsilon^{d/2-1}.

But angular momentum and symmetry often force the matrix element to vanish. If a partial-wave amplitude behaves as

Mℓ(k)∝kℓ,\mathcal M_\ell(k) \propto k^\ell,

then in three dimensions

Γℓ∝k2ℓ+1∝εℓ+1/2,\Gamma_\ell \propto k^{2\ell+1} \propto \varepsilon^{\ell+1/2},

before long-range forces or other threshold effects are included. The density of states supplies only one part of the threshold law.

Selection Rules in Transition Rates explains how symmetry determines the zeros and allowed tensor structure of the matrix element.

A finite observation time TT resolves an energy window of order

ΔET∼ℏT.\Delta E_T \sim \frac{\hbar}{T}.

Let δElevel\delta E_{\mathrm{level}} be the typical spacing of final levels in a large but finite box. A smooth continuum approximation requires many levels inside the finite-time window:

δElevel≪ℏT.\delta E_{\mathrm{level}} \ll \frac{\hbar}{T}.

Since

δElevel∼1ρ(E∗),\delta E_{\mathrm{level}} \sim \frac{1}{\rho(E_*)},

the number of sampled levels is approximately

Nwin∼ρ(E∗)ℏT.N_{\mathrm{win}} \sim \rho(E_*) \frac{\hbar}{T}.

One wants

Nwin≫1.N_{\mathrm{win}} \gg 1.

At the same time, the weighted density should vary little across that window. If EvarE_{\mathrm{var}} is its characteristic variation scale, require

ℏT≪Evar.\frac{\hbar}{T} \ll E_{\mathrm{var}}.

Together,

δElevel≪ℏT≪Evar.\delta E_{\mathrm{level}} \ll \frac{\hbar}{T} \ll E_{\mathrm{var}}.

The rate description also requires weak depletion,

ΓiT≪1,\Gamma_iT \ll 1,

unless the exponential decay or master-equation resummation is performed. These inequalities define a time window; a continuum rate is not an exact all-time law for a finite closed system.

A reliable channel sum states what is included in cc and in ρc\rho_c:

Γi=∑cΓi→c.\Gamma_i = \sum_c\Gamma_{i\to c}.

Typical discrete factors include

  • spin or helicity of the final particles;
  • photon polarization;
  • valleys, bands, subbands, or transverse modes;
  • magnetic quantum numbers;
  • right- and left-moving branches;
  • particle species and reaction channels.

Initial-state quantum numbers are usually averaged only when the preparation is an incoherent mixture or the experiment does not resolve them. Final-state quantum numbers are summed when the measurement is inclusive over them. Identical final particles may require a symmetry factor rather than a naive label sum.

A degeneracy factor must appear exactly once. If gintg_{\mathrm{int}} is already included in ρ(E)\rho(E), summing the same internal states again in ∑c\sum_c double counts them.

In nonrelativistic many-particle problems, each final momentum contributes a measure such as

V d3p(2πℏ)3,\frac{\mathcal V\,d^3p}{(2\pi\hbar)^3},

and conservation delta functions reduce the number of independent integrations. Relativistic QFT packages the corresponding final-state density into Lorentz-invariant phase space.

In natural units ℏ=c=1\hbar=c=1,

dΦn(P)=(2π)4δ(4)(P−∑r=1npr)×∏r=1nd3pr(2π)32Er.\begin{aligned} d\Phi_n(P) &= (2\pi)^4 \delta^{(4)}\left( P-\sum_{r=1}^{n}p_r \right) \\ &\quad\times \prod_{r=1}^{n} \frac{d^3p_r}{ (2\pi)^3 2E_r }. \end{aligned}

For a one-particle initial state of energy EiE_i, a decay-rate element has the schematic form

dΓ=12Ei1S∣M∣2‾dΦn,d\Gamma = \frac{1}{2E_i} \frac{1}{S} \overline{ \lvert\mathcal M\rvert^2 } d\Phi_n,

where SS accounts for identical-particle overcounting and the bar denotes any required spin sums or averages. Cross sections include an additional incident-flux factor.

The analogy with the nonrelativistic golden rule is structural:

rate=matrix element squared×final-state measure×conservation constraints.\begin{aligned} \text{rate} &= \text{matrix element squared} \\ &\quad\times \text{final-state measure} \\ &\quad\times \text{conservation constraints}. \end{aligned}

Relativistic normalization changes every factor, so one should not transplant a nonrelativistic density of states into a QFT formula. See Scattering: QFT Bridge Reference for the broader dictionary.

  1. List complete final-state labels. Separate continuous coordinates from discrete channels.
  2. Write the normalization and completeness relation. Do this before evaluating matrix elements.
  3. Keep the weighted sum intact. Start from ∑α∣Wαi∣2δ(Eα−E∗)\sum_\alpha\lvert W_{\alpha i}\rvert^2\delta(E_\alpha-E_*).
  4. Convert sums to integrals. Include every factor of volume, 2π2\pi, ℏ\hbar, and degeneracy.
  5. Apply the delta-function Jacobian. Retain all roots and branches of the energy constraint.
  6. Check regulator cancellation. A physical localized rate should not depend on an artificial box size.
  7. Test the continuum time window. Require many levels in ℏ/T\hbar/T and slow variation across that window.
  8. State what is summed or averaged. Include initial mixtures, final inclusiveness, and identical-particle factors.
  9. Check dimensions. The final Γ\Gamma must have units of inverse time.
  • Treating density of states as occupation. Availability and population are different objects.
  • Multiplying by ρ(E)\rho(E) after using energy-normalized states. That double counts the Jacobian.
  • Mixing box-normalized and delta-normalized matrix elements. Their dimensions and volume scaling differ.
  • Dropping roots of the energy equation. In one dimension, kk and −k-k are distinct unless a channel condition removes one.
  • Counting a degeneracy twice. Decide whether it belongs in gintg_{\mathrm{int}}, the channel sum, or the matrix element basis.
  • Pulling a rapidly varying matrix element off shell. The factorized form ∣W∣2ρ\lvert W\rvert^2\rho can fail near thresholds and resonances.
  • Using the shell Jacobian at a critical point without care. Return to the unreduced integral when ∇E=0\nabla E=0.
  • Taking the continuum limit at arbitrarily late times. A finite closed spectrum eventually resolves levels and can show recurrences.
  • Using state density as a scattering cross section. Cross sections also require incident flux and the correct asymptotic normalization.

Two final channels have the same density of states ρ(E∗)\rho(E_*). Their on-shell matrix elements are W1i=wW_{1i}=w and W2i=0W_{2i}=0. Compute the total rate and explain why 2ρ(E∗)∣w∣22\rho(E_*)\lvert w\rvert^2 is incorrect.

Solution

The channel sum gives

Γi=2πℏρ(E∗)(∣w∣2+0)=2πℏρ(E∗)∣w∣2.\begin{aligned} \Gamma_i &= \frac{2\pi}{\hbar} \rho(E_*) \left( \lvert w\rvert^2+0 \right) \\ &= \frac{2\pi}{\hbar} \rho(E_*)\lvert w\rvert^2. \end{aligned}

The second channel contributes states but no transition strength. Replacing both channel matrix elements by ww ignores the selection-rule zero. The physical object is the weighted density, not the unweighted count alone.

In dd dimensions, suppose

Wki(V)=V−1/2Mi(k).W_{\mathbf k i}^{(\mathcal V)} = \mathcal V^{-1/2}\mathcal M_i(\mathbf k).

Show explicitly that the golden-rule rate is independent of V\mathcal V in the continuum limit.

Solution

Use

∑k⟶V(2π)d∫ddk\sum_{\mathbf k} \longrightarrow \frac{\mathcal V}{(2\pi)^d} \int d^dk

and

∣Wki(V)∣2=1V∣Mi(k)∣2.\left\lvert W_{\mathbf k i}^{(\mathcal V)} \right\rvert^2 = \frac{1}{\mathcal V} \lvert\mathcal M_i(\mathbf k)\rvert^2.

Then

Γi=2πℏV(2π)d∫ddk ×∣Mi(k)∣2Vδ(Ek−E∗)=2πℏ∫ddk(2π)d ∣Mi(k)∣2δ(Ek−E∗).\begin{aligned} \Gamma_i &= \frac{2\pi}{\hbar} \frac{\mathcal V}{(2\pi)^d} \int d^dk\, \\ &\quad\times \frac{\lvert\mathcal M_i(\mathbf k)\rvert^2}{\mathcal V} \delta(E_{\mathbf k}-E_*) \\ &= \frac{2\pi}{\hbar} \int\frac{d^dk}{(2\pi)^d}\, \lvert\mathcal M_i(\mathbf k)\rvert^2 \delta(E_{\mathbf k}-E_*). \end{aligned}

The volume from state counting cancels the inverse volume from the squared normalized matrix element.

For ε=ℏ2k2/(2m)\varepsilon=\hbar^2k^2/(2m), derive the scaling of the free-particle density of states in dd dimensions and specialize to d=1,2,3d=1,2,3.

Solution

The number of states in a shell is proportional to

kd−1dk.k^{d-1}dk.

Since

dεdk=ℏ2km,\frac{d\varepsilon}{dk} = \frac{\hbar^2k}{m},

one has

ρd(ε)∝kd−1dε/dk∝kd−2.\rho_d(\varepsilon) \propto \frac{k^{d-1}}{d\varepsilon/dk} \propto k^{d-2}.

Using k∝εk\propto\sqrt\varepsilon gives

ρd(ε)∝εd/2−1.\rho_d(\varepsilon) \propto \varepsilon^{d/2-1}.

Therefore

ρ1∝ε−1/2,ρ2∝ε0,ρ3∝ε1/2.\rho_1\propto\varepsilon^{-1/2}, \qquad \rho_2\propto\varepsilon^0, \qquad \rho_3\propto\varepsilon^{1/2}.

Assume ⟨q∣q′⟩=δ(q−q′)\langle q\vert q'\rangle=\delta(q-q') and a monotonic dispersion E(q)E(q). Verify that

∣E⟩=∣q(E)⟩∣dE/dq∣\lvert E\rangle = \frac{\lvert q(E)\rangle}{ \sqrt{\lvert dE/dq\rvert} }

satisfies ⟨E∣E′⟩=δ(E−E′)\langle E\vert E'\rangle=\delta(E-E').

Solution

The coordinate delta function transforms as

δ(q(E)−q(E′))=∣dEdq∣δ(E−E′).\delta\big(q(E)-q(E')\big) = \left\lvert \frac{dE}{dq} \right\rvert \delta(E-E').

Therefore

⟨E∣E′⟩=δ(q(E)−q(E′))∣dE/dq∣E∣dE/dq∣E′=δ(E−E′).\begin{aligned} \langle E\vert E'\rangle &= \frac{ \delta\big(q(E)-q(E')\big) }{ \sqrt{\lvert dE/dq\rvert_E} \sqrt{\lvert dE/dq\rvert_{E'}} } \\ &= \delta(E-E'). \end{aligned}

On the support of the delta function, the two Jacobians are equal. If the dispersion has several roots at the same energy, each branch needs a separate channel label.

Near E∗E_*, let the mean final-level spacing be δE\delta E, let the weighted density vary on scale EvarE_{\mathrm{var}}, and let the predicted rate be Γ\Gamma. State inequalities on TT for a continuum golden-rule calculation without substantial depletion.

Solution

Many levels must lie inside the finite-time energy window:

δE≪ℏT.\delta E \ll \frac{\hbar}{T}.

The weighted density must be nearly constant across that window:

ℏT≪Evar.\frac{\hbar}{T} \ll E_{\mathrm{var}}.

Finally, first-order depletion must remain small:

ΓT≪1.\Gamma T \ll 1.

Thus a useful interval requires

ℏEvar≪T≪min⁡(ℏδE,1Γ),\frac{\hbar}{E_{\mathrm{var}}} \ll T \ll \min\left( \frac{\hbar}{\delta E}, \frac{1}{\Gamma} \right),

up to numerical factors set by the precise window convention. If no such interval exists, a constant continuum rate is not controlled.

For a decay into two identical spinless particles, list the ingredients beyond ∣M∣2\lvert\mathcal M\rvert^2 that are required in a relativistic rate.

Solution

The rate requires

  • the invariant one-particle measures d3pr/[(2π)32Er]d^3p_r/[(2\pi)^3 2E_r];
  • a four-momentum delta function enforcing total energy and momentum conservation;
  • the initial normalization factor 1/(2Ei)1/(2E_i);
  • a symmetry factor 1/2!1/2! for two identical final particles when phase space labels would otherwise count each physical state twice;
  • any required sums or averages over unresolved quantum numbers.

For a cross section rather than a decay rate, an incident-flux factor is also required. All factors must use one consistent relativistic state-normalization convention.

  • E. Fermi, Nuclear Physics, University of Chicago Press, 1950.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976.
  • M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Addison-Wesley, 1995.