2.5: Fermi's Golden Rule in a 1D Resonance Model
Fermi’s golden rule can be used to study a wide variety of quantum systems, and we will see examples of its usefulness in subsequent chapters. In this section, we will apply it to the simple 1D model from Section 2.2, and evaluate how well it works. The model potential is
\[\begin{align}\begin{aligned}&V(x) = V_0(x) + V_1(x), \;\;\;\mathrm{where} \\ &V_0(x) = -U \, \Theta(a-|x|), \;\; V_1(x) = V_b\, \Theta(b-|x|), \;\;\; 0<U<V_b.\end{aligned}\end{align}\]
The potential well \(V_0\) supports one or more bound states. For simplicity, we focus on the ground state, whose energy \(E_0\) lies in the range \([-U,0]\) . Once \(V_1\) is introduced, this will turn into the lowest quasi-bound state, with resonance energy \(E_{\mathrm{res}} \approx E_0 + V_b > 0\) .
When the potential consists only of \(V_0\) , the ground state wavefunction can be obtained by solving the Schrödinger wave equation for the boundary condition \(\varphi(x) \rightarrow 0\) as \(|x|\rightarrow\infty\) . It has the form
\[\varphi(x) = \begin{cases}\mathcal{A}\,\cos(qx), & |x| < a \\ \mathcal{B} \, \exp\left(-\eta|x|\right), & |x| \ge a,\end{cases}\]
where \(\mathcal{A}\) and \(\mathcal{B}\) are constants to be determined, and
\[q = \sqrt{\frac{2m}{\hbar^2}(E_0+U)}, \;\; \eta = \sqrt{\frac{2m}{\hbar^2}|E_0|}.\]
By matching both \(\varphi(x)\) and \(d\varphi/dx\) across the \(x=a\) interface, we can derive a transcendental equation for \(E_0\) , which can be solved numerically. Once we know \(E_0\) , we also know \(q\) and \(\eta\) . We can then relate \(\mathcal{A}\) and \(\mathcal{B}\) by matching \(\varphi(x)\) across the \(x=a\) interface:
\[\mathcal{A} = \mathcal{B}\; \frac{\exp(-\eta a)}{\cos(qa)}.\]
Moreover, the wavefunction must be normalized to unity:
\[\begin{align} \begin{aligned}1 &= \mathcal{A}^2 \int_{-a}^a \cos^2(qx) + 2\mathcal{B}^2 \int_a^\infty dx\; \exp(-2\eta x) \\ &= \mathcal{A}^2 \left(a+\frac{\sin(2qa)}{2q}\right) + \mathcal{B}^2\, \frac{\exp(-2\eta a)}{\eta}.\end{aligned}\end{align}\]
Putting the last two equations together yields
\[\mathcal{B}^2 = \frac{\exp\big(2\eta a\big)}{a} \left[\frac{1+\sin(2qa)/2qa}{\cos^2(qa)}+\frac{1}{\eta a}\right]^{-1}.\]
We now wish to compute the transition amplitude
\[\langle\psi_k|\hat{V}_1|\varphi\rangle \,=\, \int_{-\infty}^\infty dx \, \psi_k^*(x)\, V_1(x)\, \varphi(x) \,=\, V_b \int_{-b}^b dx \; \psi_k^*(x) \,\varphi(x),\]
where \(\psi_k(x)\) is the wavefunction of a free state for the system with potential \(V_0(x)\) , which is labeled by some continuous index \(k\) . We can use a trick to avoid calculating the exact form of \(\psi_k(x)\) . Because \(\psi_k(x)\) and \(\varphi(x)\) must be orthogonal,
\[\int_{-b}^b dx \; \psi_k^*(x) \,\varphi(x) = - \int_{|x| > a} dx \; \psi_k^*(x) \,\varphi(x).\]
When evaluating the free states outside the potential well, we can approximate them as the free states of the particle without the potential well, i.e., simple plane waves:
\[\psi_k(x) = \frac{1}{\sqrt{2\pi}}\, e^{ikx}, \;\;\;k\in\mathbb{R}, \;\;|x|>a.\]
We likewise approximate their energies by \(E_k \approx \hbar^2k^2/2m\) . Plugging \(\psi_k(x)\) into the formula for the transition amplitude, and performing the necessary integrals, yields
\[\langle\psi_k|\hat{V}_1|\varphi\rangle \,=\, - \mathcal{B} \, \sqrt{\frac{2}{\pi}} \, \frac{\hbar^2}{2m} \; \big[\eta\cos(kb) - k\sin(kb)\big] \, \exp(-\eta b).\]
Within Fermi’s golden rule, we must use a value of \(k\) such that
\[E_k = E_{\mathrm{res}} \approx E_0 + V_b \;\;\;\Rightarrow \;\;\; k \approx \sqrt{\frac{2m}{\hbar^2}(E_0+V_b)}.\]
The next thing that we need to calculate is the density of free states. This can be done by once again taking \(E_k \approx \hbar^2k^2/2m\) , and performing a change of variables:
\[\begin{align} \begin{aligned}\mathcal{D}(E) &= \int_{-\infty}^\infty dk \; \delta\left(E-\frac{\hbar^2k^2}{2m}\right) \\ &= 2 \cdot \int_0^\infty dE' \, \frac{dk}{dE'} \, \delta(E-E') \\ &= \sqrt{\frac{2m}{\hbar^2 E}}. \end{aligned}\end{align}\]
Note the factor of 2 on the second line; it is there because there is both a positive and a negative value of \(k\) for each \(E\) .
Having obtained expressions for the transition amplitude and the density of states, we just have to plug them into Fermi’s golden rule to obtain the decay rate \(\kappa\) . The figure below shows how \(\kappa\) varies with the barrier thickness \(b-a\) , with all other model parameters fixed ( \(U = 20, V_b = 30, a = 1\) ):
For comparison, the figure also plots the width of the resonant scattering peak, which according to our preceding discussion is supposed to be equal to \(\kappa\) . To obtain these values, we use the transfer matrix method to compute the \(E\) -dependent transmittance for particles incident on one side of the potential (see Appendix B). The transmittance is peaked near \(E_{\mathrm{res}}\) , and the full-width at half-maximum (see Section 2.3) is estimated numerically. The results are shown as blue dots in the figure.
We see that Fermi’s golden rule generally agrees well with the results from the transfer matrix method. There is some discrepancy, notably when the barrier is thin. Remember that Fermi’s golden rule relies on the approximation that the quasi-bound state is strongly confined (i.e., having \(\langle \varphi|G|\varphi\rangle\) be large in the derivations of Sections 2.3 and 2.4). Moreover, we have approximated the free states \(|\psi_k\rangle\) using the plane wave states in the absence of any potential, which is an assumption commonly employed when using Fermi’s golden rule.