Fermi's Golden rule is used to calculate the probability (per unit time) of a quantum mechanical transition between two quantum states. Although Fermi first coined the term "Golden Rule", Dirac developed most of the machinery.
The first part of the Golden rule is the transition matrix element (or "Amplitude") [math]M_{i,f}[/math]
- [math]| M_{i,f}| ^2 \equiv \int \psi_i^{*}(\vec{r}) H_{int} \psi_f(\vec{r}) dr^3[/math]
where
- [math]\psi_i[/math] = initial quantum state of the system which is an eigenstate of the time independent ("steady state") Hamiltonian ([math]H_0[/math])
- [math]\psi_f[/math] = final quantum state of system after a transition
- [math]H_{int}[/math] = the part of the total Hamiltonian ([math]H_{tot}[/math]) which describes the interaction responsible for the transition.
- [math]H_0[/math] = Unperturbed ("steady state") Hamiltonian
- [math]H_{tot} = H_0 + H_{int}[/math] = total Hamiltonian describing the quantum mechanical system
- [math]dr^2[/math] integration over all space
The off diagonal elements of the [math]M_{i,f}[/math] matrix tell you the transition probablility.
Single Particle decay
Consider the case when a single particle decays into multiple fragments (several other particles)
- [math]d \Gamma = \hbar W = | M |^2 \frac{S}{2 \hbar m_1} \left [ \left ( \frac{d^3\vec{p}_2}{(2 \pi)^3 2E_2} \right )\left ( \frac{d^3\vec{p}_3}{(2 \pi)^3 2E_3} \right )\left ( \frac{d^3\vec{p}_4}{(2 \pi)^3 2E_4} \right ) \cdots \left ( \frac{d^3\vec{p}_N}{(2 \pi)^3 2E_N} \right ) \right ] \times (2 \pi)^4 \delta^4(p_1^\mu-p_2^\mu-p_3^\mu- \cdots p_n^\mu)[/math]
where
- [math]W[/math] = probability per second that the particle will decay
- [math]S[/math] = a symmetry factor of [math]\left( \frac{1}{j!}\right )[/math] for every group of [math]j[/math] identical particles in the final state
- [math]p_i^\mu = (E_i, \vec{p}_i) [/math]= 4-momentum of the [math]i^{th}[/math] particle. ;[math] p_1 = (E_1,0)[/math]
- [math]\delta^4(p_1-p_2-p_3- \cdots p_n)[/math] = conservation of 4-momentum
- Note
- [math]\int \frac{d^3 \vec{p}_i}{2E_i} = \int d^4p_i \delta(p_i^2-m_i^2)[/math] = invariant under Lorentz transformations
Example: Pi-zero ([math]\pi^0[/math]) decay
- [math]\pi^0 \rightarrow \gamma + \gamma[/math]
we are interested in calculating
- [math]d \Gamma = \hbar W = | M |^2 \frac{S}{2 \hbar m_1} \left [ \left ( \frac{d^3\vec{p}_2}{(2 \pi)^3 2E_2} \right )\left ( \frac{d^3\vec{p}_3}{(2 \pi)^3 2E_3} \right ) \right ] (2 \pi)^4 \delta^4(p_1^\mu-p_2^\mu-p_3^\mu)[/math]
Consider the decay of a neutral pion ([math]\pi^0[/math]) into two photons ([math]\gamma[/math]).
The two gammas are identical particles so
- [math]S = \frac{1}{2!} = \frac{1}{1 \cdot 2} = \frac{1}{2}[/math]
Since the pion is initially at rest (or we can go to its rest fram and then Lorentz boost to back to the lab frame)
- [math]p_1^\mu = (E_1,p_1) = (E_1,0) = (m_{\pi^0},0)[/math]
Because photons have no mass, [math]m_2=m_3 = 0[/math] :
- [math]p_2^\mu = (E_2,p_2) = (|\vec{p}_2|,\vec{p}_2) = (p_2,\vec{p}_2)[/math]
- [math]p_3^\mu = (E_3,p_3) = (|\vec{p}_3|,\vec{p}_3)= (p_3,\vec{p}_3)[/math]
- [math]\delta^4(p_1^\mu-p_2^\mu-p_3^\mu) = \delta(m_{\pi} -p_2-p_3) \delta(0-\vec{p}_2 - \vec{p}_3)[/math]
- [math]\int d \Gamma = \int \frac{|M|^2}{16 (2 \pi)^6 \hbar m_{\pi} }\left [ \left ( \frac{d^3\vec{p}_2}{E_2} \right )\left ( \frac{d^3\vec{p}_3}{E_3} \right ) \right ] (2 \pi)^4 \delta(m_{\pi} -p_2-p_3) \delta(0-\vec{p}_2 - \vec{p}_3)[/math]
- [math]= \int \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi} }\left [ \left ( \frac{d^3\vec{p}_2}{E_2} \right )\left ( \frac{d^3\vec{p}_3}{E_3} \right ) \right ] \delta(m_{\pi} -p_2-p_3) \delta(0-\vec{p}_2 - \vec{p}_3)[/math]
- [math]= \int \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi}p_2 p_3 }\left [ \left ( d^3\vec{p}_2 \right )\left ( d^3\vec{p}_3 \right ) \right ] \delta(m_{\pi} -p_2-p_3) \delta(0-\vec{p}_2 - \vec{p}_3)[/math]
Integrating over[math] d^3\vec{p}_3 \Rightarrow \vec{p}_2 = - \vec{p}_3 , | p_2| = |p_3|[/math] :
- \int d \Gamma =[math]= \int \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi} p_2^2}\left [ \left ( d^3\vec{p}_2 \right ) \right ] \delta(m_{\pi} -2p_2) [/math]
- = [math]= \int \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi} p_2^2}\left [ \left ( |\vec{p}_2|^2 d|\vec{p}_2| \sin(\theta) d \theta d \phi \right ) \right ] \delta(m_{\pi} -2p_2)[/math]
If [math]|M|^2 \ne f( \vec{p}_2)[/math]
then
- [math]\int d \Gamma = \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi}}\int \frac{1}{ p_2^2}\left [ \left ( |\vec{p}_2|^2 d|\vec{p}_2| \sin(\theta) d \theta d \phi \right ) \right ] \delta(m_{\pi} -2p_2)[/math]
- [math]= \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi}}\int \delta(m_{\pi} -2p_2) d|\vec{p}_2| \sin(\theta) d \theta d \phi [/math]
- [math]= \frac{|M|^2}{(8 \pi)^2 \hbar m_{\pi}} (4 \pi) \int \delta(m_{\pi} -2p_2) d|\vec{p}_2| [/math]
- [math]= \frac{|M|^2}{(16 \pi) \hbar m_{\pi}} \int \delta(m_{\pi} -2p_2) d|\vec{p}_2| [/math]
- [math]= \frac{|M|^2}{(16 \pi) \hbar m_{\pi}} \frac{m_{\pi}}{2} [/math]
- [math]= \frac{|M|^2}{(32 \pi) \hbar} [/math]