Forest FermiGoldenRule Notes

Fermi's Golden Rule

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]

Transition Amplitude

[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,\vec{p}_2) = (|\vec{p}_2|,\vec{p}_2) = (p_2,\vec{p}_2)[/math]

[math]p_3^\mu = (E_3,\vec{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] :

[math]\int d \Gamma = \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 = f( \vec{p}_2) \ne f( \vec{p}_2,\vec{p}_3)[/math]( the transition does not depend on the momentum vector directions)

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{1}{2} \;\; : \delta (kx) = \frac{\delta(x)}{k} \Rightarrow \delta (m_{\pi} -2p_2)=\frac{\delta (\frac{m_{\pi}}{2} -p_2)}{2} [/math]

[math] \Gamma= \frac{|M|^2}{(32 \pi) \hbar m_{\pi}} [/math] with the additional conditions that [math]p_2 = \frac{m_{\pi}}{2}[/math]and [math]\vec{p}_2 = - \vec{p}_3[/math] which must be applied when evaluating[math] |M|^2[/math]

Units check

[math]\frac{1}{\hbar m_{\pi}} = \left ( 6.582 \times 10^{-10} MeV \cdot s\right ) (150 MeV) = 9.97 \times 10^{-8} MeV^2 \cdot s[/math]

[math]\left [ |M|^2\right ] = MeV^2 [/math]= energy eigenvalues squared

[math] \left [ \Gamma \right ] = \frac{MeV^2 }{MeV^2 \cdot s} = \frac{1}{s} =[/math] transition probability per unit time

Caveat

Sometime [math]|M|[/math] will depend on [math]\vec{p}_i[/math] in which case the integral must be done after evaluating the matrix element amplitude. An example of this is when the transition is spin dependent [math](\vec{p} \cdot \vec{s})[/math] as in the hyperfine interaction or polarization based transitions.

Example: Two -Body decay (fission fragments)

Now consider a more general case in whigh the decay daughters have mass:

[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]

We don't know if the two daughter particles are identical so leave [math]S[/math] as a funtion.

[math]S = ?[/math]

We can still do the calculation in the mother particles rest frame.

[math]p_1^\mu = (E_1,p_1) = (E_1,0) = (m_{\pi^0},0)[/math]

Now that the daughters have mass we need:

[math]p_2^\mu = (E_2,\vec{p}_2) \; E_2 = \sqrt{p_2^2 + m_2^2}[/math]

[math]p_3^\mu = (E_3,\vec{p}_3) \; E_3 = \sqrt{p_3^2 + m_3^2}[/math]

[math]\int d \Gamma = \int \frac{S|M|^2}{8 (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]

Recast the delta function:

[math]\delta^4(p_1^\mu-p_2^\mu-p_3^\mu) = \delta(m_{\pi} -E_2-E_3) \delta(0-\vec{p}_2 - \vec{p}_3)[/math]

Upon integrating over [math]d^3 \vec{p}_3[/math]

[math]\delta(0-\vec{p}_2 - \vec{p}_3) \Rightarrow |\vec{p}_2| = | \vec{p}_3| \Rightarrow E_3 = \sqrt{p_3^2 + m_3^2} = E_3 = \sqrt{p_2^2 + m_3^2}[/math]

[math]\delta^4(p_1^\mu-p_2^\mu-p_3^\mu) = \delta(m_{\pi} -\sqrt{p_2^2 + m_2^2} - \sqrt{p_2^2 + m_3^2}) [/math]

After integrating over [math]d^3 \vec{p}_3[/math] the delta function gives you

[math]\int d \Gamma = \int \frac{S|M|^2}{8 (2 \pi)^2 \hbar m_{\pi} }\left [ \left ( \frac{d^3\vec{p}_2}{E_2 \sqrt{p_2^2 + m_3^2}} \right ) \right ] \delta(m_{\pi} -\sqrt{p_2^2 + m_2^2} - \sqrt{p_2^2 + m_3^2}) [/math]