Skip to main content
Physics LibreTexts

16.6: Exercises - Time-Dependence and Pseudopotentials

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    One harder problem, which you are nevertheless encouraged to try, is available online

    For the first two problems, you may assume that the hydrogen eigenfunctions are:

    \[u_{100} = \left(\pi a_{0}^{3}\right)^{-1 / 2} \exp \left(-r / a_{0}\right) \nonumber\]

    \[ u_{211} = -\left(\pi a_{0}^{3}\right)^{-1 / 2} \frac{r}{8 a_{0}} \sin \theta \exp (i \phi) \exp \left(-r / 2 a_{0}\right) \nonumber\]

    \[u_{210} = \left(8 \pi a_{0}^{3}\right)^{-1 / 2} \frac{r}{2 a_{0}} \cos \theta \exp \left(-r / 2 a_{0}\right) \nonumber\]

    \[u_{21-1} = \left(\pi a_{0}^{3}\right)^{-1 / 2} \frac{r}{8 a_{0}} \sin \theta \exp (-i \phi) \exp \left(-r / 2 a_{0}\right) \nonumber\]


    \[\int_{0}^{\infty} \exp (-b r) r^{n} \mathrm{d} r = n ! / b^{n+1}, \quad n>-1 \nonumber\]

    1. A hydrogen atom is placed in a uniform but time-dependent electric field of magnitude:

    \[\mathcal{E} = 0 \text{ for } t < 0, \quad \mathcal{E} = \mathcal{E}_0 \text{ exp}(−t/\tau ) \; \text{ for } t \geq 0 \; (\tau > 0) \nonumber\]

    where \(\mathcal{E}_{\prime}\) is a constant. At time \(t = 0\), the atom is in the ground \((1s)\) state. Show that the probability, to lowest order in perturbation theory, that as \(t \rightarrow \infty\), the atom is in the \(2p\) state in which the component of the orbital angular momentum in the direction of the field is zero, is given by

    \[p_{1s \rightarrow 2p} = |c(\infty )|^2 = \frac{2^{15}}{3^{10}} \frac{(e\mathcal{E}_0a_0)^2}{(E_{2p} − E_{1s})^2} + (\hbar /\tau )^2 \nonumber\]

    What is the probability that it is in the \(2s\)-state?

    [Hint: take the field direction to be the \(z\)-direction. Write down the potential energy of the electron in the given field and treat as a time-dependent perturbation].

    2. The neutral kaon states produced via the strong interaction are \(|K^0\rangle\) and \(|\overline{K}^0 \rangle\), eigenvalues of strangeness \(\hat{S}\) such that \(\hat{S} |K^0 \rangle = |K^0 \rangle\) and \(\hat{S} |\overline{K}^0 \rangle = −|\overline{K}^0 \rangle\). Acting on these states with a weak force-related operator \(\hat{CP}\) gives:

    \[\hat{CP} |K^0 \rangle = |\overline{K}^0 \rangle; \quad \hat{CP} |\overline{K}^0 \rangle = |K^0 \rangle \nonumber\]

    Evaluate the CP eigenstates for neutral kaons \(|K_1\rangle\) and \(|K_2\rangle\) with eigenvalues \(CP = \pm 1\).

    Kaons decay into pions via the weak force with lifetimes:

    \[\tau_1 = 0.9 \times 10^{−10} (\text{for } |K_1 \rangle ) \quad \text{ and } \quad \tau_2 = 0.5 \times 10^{−7} (\text{for } |K_2\rangle ). \nonumber\]

    Show that the wavefunction \(a_1(t) = a_1(0)e^{ −t/2\tau} e^{−iEt/\hbar} |\Phi ({\bf r})\rangle\), where \(E\) is the energy, represents decay with lifetime \(\tau\), and that the amplitudes of the \(|K_1\rangle\) and the \(|K_2\rangle\) states at rest (i.e. with \(E = mc^2 \)) are:

    \[a_1(t) = a_1(0)e^{−t/2\tau_1} e^{−im_1c^2 t/\hbar} \quad \text{ and } \quad a_2(t) = a_2(0) e^{−t/2\tau_2} e^{−im_2c^2 t/\hbar} \nonumber\]

    By considering a general state \(a|K^0 \rangle + b|\overline{K^0} \rangle\) how that the intensity of \(K^0\) (i.e. \(|a|^2\)) is measured by the operator \(\frac{1}{2} (\hat{S} + 1)\). What is the operator for the intensity of \(|\overline{K}^0 \rangle\)?

    At \(t=0\) a kaon beam is in a pure \(|K^0 \rangle\) state, with intensity proportional to \(|a_0(0)|^2 = 1\) show that at time t

    \[|a_0(t)|^2 = \langle \frac{1}{2} (\hat{S} + 1) \rangle = \frac{1}{4} \left[ e^{−t/\tau_1} + e^{−t/\tau_2} − 2e^{−t/2\tau_1} e^{−t/2\tau_2} \cos (m_{12}t) \right] \nonumber\]

    \[\text{and } I(\overline{K}^0 ) = \frac{1}{4} \left[ e^{−t/\tau_1} + e^{−t/\tau_2} − 2e^{−t/2\tau_1} e^{−t/2\tau_2} \cos (m_{12}t) \right] \nonumber\]

    where \(m_{12} = (m_2 − m_1)c^2/\hbar\). What is there a physical reason for the intensities being unchanged if we use \(m_{12} = (m_1 − m_2)c^2/\hbar\)?

    Sketch, as a function of time, the expectation values of: \(\frac{1}{2} (\hat{S} +1), \frac{1}{2} (1−\hat{S} ), \hat{S} , \hat{CP} , \frac{1}{2} (\hat{CP} +1)\), in a system which began in state \(|K^0 \rangle\).


    Figure \(\PageIndex{1}\)

    Figure \(\PageIndex{1}\) shows an experiment where kaons are generated in state \(|K^0 \rangle\). After \(10^{−9}\) seconds, assuming that all the \(|K_1 \rangle\) particles have decayed, the kaons pass through a small region of matter, where they interact via the strong interaction. Very soon after, the kaons leave the matter and move into a region of vacuum where they begin to decay via the weak interaction \((\hat{CP} )\). Assume that all coherence between this region and the previous region is lost i.e. the wavefunction is completely collapsed onto its strong interaction eigenstates. Evaluate the appropriate eigenstates and intensities just before the beam enters the matter and just after it leaves the matter. What is the total intensity of kaons and antikaons which survive a further \(10^{-9}\) seconds?

    Had it not been for the matter, what would have been the total intensity of kaons and antikaons after \(2 \times 10^{-9}\)?

    3. Compare question 3 to a system of circular and plane polarisers and light beams.

    4. This question illustrates the principle of the pseudopotential.

    A particle is bound in 1D by a potential which has a complicated form for \(|x| < x_c\) but is zero outside this “cut-off” radius. It is known to have a bound eigenstate with energy \(−E_0\). Show that in this region of space, the wavefunction can be written as

    \[\Phi (x > x_c) = a \text{ exp}(−k|x|) \nonumber\]

    and determine \(k\). What can you say about \(a\)?

    Now suppose that we know that the normalisation constant \(a = a_0\). Show that the ground state of a finite square well pseudopotential can be used to give exactly the same wavefunction i.e. \(\Phi_{P S}(x) = \Phi (x)\), for \(x > x_c\).

    Writing the ground state wavefunction of the square well as:

    \[\Phi_{P S}(x) = b \cos (k_1x) \quad |x| < x_c \nonumber\]

    \[\Phi_{P S}(x) = \Phi (x) = a \text{ exp}(−k|x|) \quad |x| > x_c \nonumber\]

    where \(k_1 = \sqrt{2m(V − E_0)}/\hbar\), determine three simultaneous equations for the required values of the finite well depth \(V\) and range \(x_c\), and the normalization constant \(b\).

    When are such pseudopotentials useful?

    5. * The electric dipole moment operator is \(\underline{\hat{D}} \equiv −e\underline{r}\). The position vector can be written

    \[\underline{r} = r \{\underline{e}_1 \sin \theta \ cos \phi + \underline{e}_2 \sin \theta \sin \phi + \underline{e}_3 \cos \theta \} \nonumber\]

    where \(\underline{e}_i , i = 1, 2, 3\), are the usual Cartesian unit vectors in the \(x, y, z\) directions and \(\theta , \phi\) are the polar and azimuthal angles in spherical polar coordinates.

    Calculate the dipole matrix elements for the radiative transition from the \(n = 2\) states to the \(1s\) state of atomic hydrogen.

    The Einstein spontaneous transition rate for the \(2p \rightarrow 1s\) transition is given by

    \[R^{\text{spon}}_{mk} = \frac{\omega^3_{mk}}{3\pi c^3 \hbar \epsilon_0} |\langle m|\underline{\hat{D}} |k\rangle |^2 = \frac{e^2\omega^3_{mk}}{3\pi c^3 \hbar \epsilon_0} |\langle m|\underline{\hat{r}}|k \rangle |^2 \nonumber\]

    where \(k\) and \(m\) label initial (2p) and final (1s) states and \(\alpha \equiv e 2/(4\pi \epsilon_0)\hbar c\) is the fine structure constant.

    Explain in words how this relates to the Fermi Golden Rule.

    Assuming that an initial \(2p\) state is unpolarised; that is, each of the three possible values of \(m_{\ell}\) is equally likely, show that this is equal to.

    \[R^{spon}_{2p \rightarrow 1s} = \left(\frac{2}{3} \right)^8 \frac{mc^2}{\hbar} \alpha^5 \nonumber\]

    This page titled 16.6: Exercises - Time-Dependence and Pseudopotentials is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Graeme Ackland via source content that was edited to the style and standards of the LibreTexts platform.