Students calculate the expectation value of energy and angular momentum as a function of time for an initial state for a particle on a ring. This state is a linear combination of energy/angular momentum eigenstates written in bra-ket notation.

Expectation Values and Uncertainty

You have a system that consists of quantum particles with spin. On this system, you will perform a Stern-Gerlach experiment with an analyzer oriented in the \(z\)-direction.

Consider one of the different initial spin states described below:

A spin 1/2 particle described by:

1. \(\left|{+}\right\rangle \)
2. \(\frac{i}{2}\left|{+}\right\rangle -\frac{\sqrt{3}}{2}\left|{-}\right\rangle \)

3. \(\left|{+}\right\rangle _x\)

   A spin 1 particle described by:

4. \(\left|{0}\right\rangle \)
5. \(\left|{-1}\right\rangle _x\)
6. \(\frac{2}{3}\left|{1}\right\rangle +\frac{i}{3}\left|{0}\right\rangle -\frac{2}{3}\left|{-1}\right\rangle \)

• List the possible values of spin you could measure and determine the probability associated with each value of the z-component of spin.

  
  

• Plot a histogram of the probabilities.

  
  

• Find the expectation value of the z-component of spin.

  
  

• Find the uncertainty of the z-component of spin.

Introduction

I like to break this activity into two parts:

(1) Calculating expectation values and relating them to the associated distributions of the probabilities of results, and

(2) Calculating the quantum uncertainty of the state and relating the uncertainty to distributions of the probabilities of results.

Therefore, I have my students do the first part of the activity before I introduce quantum uncertainty.

I introduce the activity by reminding students about two ways of calculating the expectation value. Given a quantum state \(\left|{\psi}\right\rangle \), for a measurement of an observable represented by an operator \(\hat{A}\) with eigenstates \(\left|{a_i}\right\rangle \) and eigenvalues \(a_i\): \begin{align*} \langle \hat{A} \rangle &= \sum_{i} a_i\mathcal{P}(a_i) \\ &= \left\langle {\psi}\right|\hat{A}\left|{\psi}\right\rangle \end{align*}

After the students calculate expectation values and we have a whole class discussion about 1 of the examples, then I do a lecture introducing the quantity of quantum uncertainty (relating it to the standard deviation of the distribution of probabilities by spin component value) and deriving the simplified equation:

\[\Delta A = \sqrt{\langle A^2 \rangle - \langle A \rangle ^2}\]

Student Conversations

1. One could have each group report out, or the instructor could discuss a few key examples.

   For expectation value, I like to talk about Case 2: \(\frac{i}{2}\left|{+}\right\rangle -\frac{\sqrt{3}}{2}\left|{-}\right\rangle \), where the probabilities of the two outcomes are not equal to show how the weighting plays out. Also, the expectation value is not a possible measurement value, and I like to talk about that. “Expectation” value is a misleading name for this quantity - it characterizes the distribution and is not necessarily a result of an individual measurement.

   I also like to discuss an example like Case 5: \(\left|{1}\right\rangle _x\) where the distribution is symmetric around \(0\hbar\).

2. I think it's important to encourage students to calculate expectation values both ways (with probabilities and as a bracket with matrix notation) while the teaching team is available to help them.

3. For quantum uncertainty, I like to talk about an example like Case 3: \(\left|{+}\right\rangle _x\) where all the individual measurements are the same ”distance” away from the expectation value as a sensemaking exercise to connect to a conceptual interpretation of physics.

   I also like to discuss an example like Case 5: \(\left|{-1}\right\rangle _x\), where the fact that we're taking an rms average is apparent: half the measurements are \(\hbar\) away from the expectation value and the other half are \(0\hbar\) away, but the uncertainty is \(\hbar/\sqrt{2}\).

This activity reinforces the strategies students have been practicing on each system by letting them create their own matrix operators and columns on the hydrogen atom and do some calculations with them.

Students are asked to find eigenvalues, probabilities, and expectation values for \(H\), \(L^2\), and \(L_z\) for a superposition of \(\vert n \ell m \rangle\) states. This can be done on small whiteboards or with the students working in groups on large whiteboards.

Students then work together in small groups to find the matrices that correspond to \(H\), \(L^2\), and \(L_z\) and to redo \(\langle E\rangle\) in matrix notation.

You have a system that consists of two identical (fair) six-sided dice. Imagine that you will perform an experiment where you roll the pair of dice together and record the observable: the norm of the difference between the values displayed by the two dice.

1. What are the possible results of the observable for each roll?

   
   

2. What is the theoretical probability of measuring each of those results? Assume the results are fair.

   Plot a probability histogram. Use your histogram to make a guess about where the average value is and the standard deviation.

   
   

3. Use your theoretical probabilities to determine a theoretical average value of the observable (the expectation value)? Indicate the expectation value on your histogram.

   
   

4. Use your theoretical probabilities to determine the standard deviation (the uncertainty) of the distribution of possible results. Indicate the uncertainty on your histogram.

   
   

5. Challenge: Use

   1. Dirac bra-ket notation
   2. matrices

   to represent:

   • the possible states of the dice after a measurement is made;

   • the state of the dice when you're shaking them up in your hand;

   • an operator that represents the norm of the difference of the dice.

Consider a two-state quantum system with a Hamiltonian \begin{equation} \hat{H}\doteq \begin{pmatrix} E_1&0\\ 0&E_2 \end{pmatrix} \end{equation} Another physical observable \(M\) is described by the operator \begin{equation} \hat{M}\doteq \begin{pmatrix} 0&c\\ c&0 \end{pmatrix} \end{equation} where \(c\) is real and positive. Let the initial state of the system be \(\left|{\psi(0)}\right\rangle =\left|{m_1}\right\rangle \), where \(\left|{m_1}\right\rangle \) is the eigenstate corresponding to the larger of the two possible eigenvalues of \(\hat{M}\). What is the expectation value of \(M\) as a function of time? What is the frequency of oscillation of the expectation value of \(M\)?

This activity lets students explore translating a wavefunction that isn't obviously made up of eigenstates at first glance into ket and matrix form. Then students explore wave functions, probabilities in a region, expectation values, and what wavefunctions can tell you about measurements of \(L_z\).

Students work in small groups to use the superposition principle \[\vec{A}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert}\, d\tau^{\prime}\] to find an integral expression for the magnetic vector potential, \(\vec{A}(\vec{r})\), due to a spinning ring of charge.

In an optional extension, students find a series expansion for \(\vec{A}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students work in small groups to use the Biot-Savart law \[\vec{B}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})\times \left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the magnetic field, \(\vec{B}(\vec{r})\), due to a spinning ring of charge.

In an optional extension, students find a series expansion for \(\vec{B}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students work in small groups to use Coulomb's Law \[\vec{E}(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})\left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the electric field, \(\vec{E}(\vec{r})\), everywhere in space, due to a ring of charge.

In an optional extension, students find a series expansion for \(\vec{E}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students work in small groups to use the superposition principle \[V(\vec{r}) = \frac{1}{4\pi\epsilon_0}\sum_i \frac{q_i}{\vert\vec{r}-\vec{r}_i\vert}\] to find the electrostatic potential \(V\) everywhere in space due to a pair of charges (either identical charges or a dipole). Different groups are assigned different arrangements of charges and different regions of space to consider: either on the axis of the charges or in the plane equidistant from the two charges, for either small or large values of the relevant geometric variable. Each group is asked to find a power series expansion for the electrostatic potential, valid in their group's assigned region of space. The whole class wrap-up discussion then compares and contrasts the results and discuss the symmetries of the two cases.

Students work in small groups to use the superposition principle \[V(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert} \, d\tau^{\prime}\] to find an integral expression for the electrostatic potential, \(V(\vec{r})\), everywhere in space, due to a ring of charge.

In an optional extension, students find a series expansion for \(V(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students use Mathematica to visualize the probability density distribution for the hydrogen atom orbitals with the option to vary the values of \(n\), \(\ell\), and \(m\).

Students draw the 3D graphs of equations using three variables. They make choices for drawing a stack of curves in parallel planes and a curve in a perpendicular plane (e.g. substituting in values for \(x\), \(y\), or \(z\). )

• The gradient is perpendicular to the level curves.
• The gradient is a local quantity, i.e. it only depends on the values of the function at infinitesimally nearby points.
• Although students learn to chant that "the gradient points uphill," the gradient does not point to the top of the hill.
• The gradient path is not the shortest path between two points.