Activities
Derivative of Fermi-Dirac function Show that the magnitude of the slope of the Fermi-Direc function \(f\) evaluated at the Fermi level \(\varepsilon =\mu\) is inversely proportional to its temperature. This means that at lower temperatures the Fermi-Dirac function becomes dramatically steeper.
These lecture notes from week 7 of https://paradigms.oregonstate.edu/courses/ph441 apply the grand canonical ensemble to fermion and bosons ideal gasses. They include a few small group activities.
For electrons with an energy \(\varepsilon\gg mc^2\), where \(m\) is the mass of the electron, the energy is given by \(\varepsilon\approx pc\) where \(p\) is the momentum. For electrons in a cube of volume \(V=L^3\) the momentum takes the same values as for a non-relativistic particle in a box.
Show that in this extreme relativistic limit the Fermi energy of a gas of \(N\) electrons is given by \begin{align} \varepsilon_F &= \hbar\pi c\left(\frac{3n}{\pi}\right)^{\frac13} \end{align} where \(n\equiv \frac{N}{V}\) is the number density.
Show that the total energy of the ground state of the gas is \begin{align} U_0 &= \frac34 N\varepsilon_F \end{align}
Show that a Fermi electron gas in the ground state exerts a pressure \begin{align} p = \frac{\left(3\pi^2\right)^{\frac23}}{5} \frac{\hbar^2}{m}\left(\frac{N}{V}\right)^{\frac53} \end{align} In a uniform decrease of the volume of a cube every orbital has its energy raised: The energy of each orbital is proportional to \(\frac1{L^2}\) or to \(\frac1{V^{\frac23}}\).
Find an expression for the entropy of a Fermi electron gas in the region \(kT\ll \varepsilon_F\). Notice that \(S\rightarrow 0\) as \(T\rightarrow 0\).
Problem
(K&K 7.11) Show for a single orbital of a fermion system that \begin{align} \left<(\Delta N)^2\right> = \left<N\right>(1+\left<N\right>) \end{align} if \(\left<N\right>\) is the average number of fermions in that orbital. Notice that the fluctuation vanishes for orbitals with energies far enough from the chemical potential \(\mu\) so that \(\left<N\right>=1\) or \(\left<N\right>=0\).
Problem
The following two problems ask you to make Fermi estimates. In a good Fermi estimate, you start from basic scientific facts you already know or quantities that you can reasonably estimate based on your life experiences and then reason your way to estimate a quantity that you would not be able guess. You may look up useful conversion factors or constants. Use words, pictures, and equations to explain your reasoning:
- Imagine that you send a pea-sized bead of silver through a Stern-Gerlach device oriented to measure the z-component of intrinsic spin. Estimate the total z-component of the intrinsic spin of the ball you would measure in the HIGHLY improbable case that every atom is spin up.
- Protons, neutrons, and electrons are all spin-1/2 particles. Give a (very crude) order of magnitude estimate of the number of these particles in your body.
Let us imagine a new mechanics in which the allowed occupancies of an orbital are 0, 1, and 2. The values of the energy associated with these occupancies are assumed to be \(0\), \(\varepsilon\), and \(2\varepsilon\), respectively.
Derive an expression for the ensemble average occupancy \(\langle N\rangle\), when the system composed of this orbital is in thermal and diffusive contact with a resevoir at temperature \(T\) and chemical potential \(\mu\).
Return now to the usual quantum mechanics, and derive an expression for the ensemble average occupancy of an energy level which is doubly degenerate; that is, two orbitals have the identical energy \(\varepsilon\). If both orbitals are occupied the toal energy is \(2\varepsilon\). How does this differ from part (a)?
In this small group activity, students multiply a general 3x3 matrix with standard basis row/column vectors to pick out individual matrix elements. Students generate the expressions for the matrix elements in bra/ket notation.
Problem
- \(\left\langle {\Psi}\middle|{\Psi}\right\rangle =1\) Identify and discuss the dimensions of \(\left|{\Psi}\right\rangle \).
- For a spin \(\frac{1}{2}\) system, \(\left\langle {\Psi}\middle|{+}\right\rangle \left\langle {+}\middle|{\Psi}\right\rangle + \left\langle {\Psi}\middle|{-}\right\rangle \left\langle {-}\middle|{\Psi}\right\rangle =1\). Identify and discuss the dimensions of \(\left|{+}\right\rangle \) and \(\left|{-}\right\rangle \).
- In the position basis \(\int \left\langle {\Psi}\middle|{x}\right\rangle \left\langle {x}\middle|{\Psi}\right\rangle dx = 1\). Identify and discuss the dimesions of \(\left|{x}\right\rangle \).
Given the polar basis kets written as a superposition of Cartesian kets \begin{eqnarray*} \left|{\hat{s}}\right\rangle &=& \cos\phi \left|{\hat{x}}\right\rangle + \sin\phi \left|{\hat{y}}\right\rangle \\ \left|{\hat{\phi}}\right\rangle &=& -\sin\phi \left|{\hat{x}}\right\rangle + \cos\phi \left|{\hat{y}}\right\rangle \end{eqnarray*}
Find the following quantities: \[\left\langle {\hat{x}}\middle|{\hat{s}}\right\rangle ,\quad \left\langle {\hat{y}}\middle|{{\hat{s}}}\right\rangle ,\quad \left\langle {\hat{x}}\middle|{\hat{\phi}}\right\rangle ,\quad \left\langle {\hat{y}}\middle|{\hat{\phi}}\right\rangle \]
- Given a vector written in the polar basis \[\left|{\vec{v}}\right\rangle = a\left|{\hat{s}}\right\rangle + b\left|{\hat{\phi}}\right\rangle \] where \(a\) and \(b\) are known. Find coefficients \(c\) and \(d\) such that \[\left|{\vec{v}}\right\rangle = c\left|{\hat{x}}\right\rangle + d\left|{\hat{y}}\right\rangle \] Do this by using the completeness relation: \[\left|{\hat{x}}\right\rangle \left\langle {\hat{x}}\right| + \left|{\hat{y}}\right\rangle \left\langle {\hat{y}}\right| = 1\]
- Using a completeness relation, change the basis of the spin-1/2 state \[\left|{\Psi}\right\rangle = g\left|{+}\right\rangle + h\left|{-}\right\rangle \] into the \(S_y\) basis. In otherwords, find \(j\) and \(k\) such that \[\left|{\Psi}\right\rangle = j\left|{+}\right\rangle _y + k\left|{-}\right\rangle _y\]
For this problem, use the vectors \(|a\rangle = 4 |1\rangle - 3 |2\rangle\) and \(|b\rangle = -i |1\rangle + |2\rangle\).
- Find \(\langle a | b \rangle\) and \(\langle b | a \rangle\). Discuss how these two inner products are related to each other.
- For \(\hat{Q}\doteq \begin{pmatrix} 2 & i \\ -i & -2 \end{pmatrix} \), calculate \(\langle1|\hat{Q}|2\rangle\), \(\langle2|\hat{Q}|1\rangle\), \(\langle a|\hat{Q}| b \rangle\) and \(\langle b|\hat{Q}|a \rangle\).
- What kind of mathematical object is \(|a\rangle\langle b|\)? What is the result if you multiply a ket (for example, \(| a\rangle\) or \(|1\rangle\)) by this expression? What if you multiply this expression by a bra?
This activity acts as a reintroduction to doing quantum calculations while also introducing the matrix representation on the ring, allowing students to discover how to index and form a column vector representing the given quantum state. In addition, this activity introduces degenerate measurements on the quantum ring and examines the state after measuring both degenerate and non-degenerate eigenvalues for the state.
One way to write volume charge densities without using piecewise functions is to use step \((\Theta)\) or \(\delta\) functions. Consider a spherical shell with charge density \[\rho (\vec{r})=\alpha3e^{(k r)^3} \]
between the inner radius \(a\) and the outer radius \(b\). The charge density is zero everywhere else.
- (2 pts) What are the dimensions of the constants \(\alpha\) and \(k\)?
- (2 pts) By hand, sketch a graph the charge density as a function of \(r\) for \(\alpha > 0\) and \(k>0\) .
- (2 pts) Use step functions to write this charge density as a single function valid everywhere in space.
Students compute inner products to expand a wave function in a sinusoidal basis set. This activity introduces the inner product for wave functions, and the idea of approximating a wave function using a finite set of basis functions.
Students calculate the Fourier transform of the Dirac delta function.
Remember that the delta function is defined so that \[ \delta(x-a)= \begin{cases} 0, &x\ne a\\ \infty, & x=a \end{cases} \]
Also: \[\int_{-\infty}^{\infty} \delta(x-a)\, dx =1\].
- Find a set of functions that approximate the delta function \(\delta(x-a)\) with a sequence of isosceles triangles \(\delta_{\epsilon}(x-a)\), centered at \(a\), that get narrower and taller as the parameter \(\epsilon\) approaches zero.
- Using the test function \(f(x)=3x^2\), find the value of \[\int_{-\infty}^{\infty} f(x)\delta_{\epsilon}(x-a)\, dx\] Then, show that the integral approaches \(f(a)\) in the limit that \(\epsilon \rightarrow 0\).
Consider the normalized wavefunction \(\Phi\left(\phi\right)\) for a quantum mechanical particle of mass \(\mu\) constrained to move on a circle of radius \(r_0\), given by: \begin{equation} \Phi\left(\phi\right)= \frac{N}{2+\cos(3\phi)} \end{equation} where \(N\) is the normalization constant.
Find \(N\).
- Plot this wave function.
- Plot the probability density.
- Find the probability that if you measured \(L_z\) you would get \(3\hbar\).
- What is the expectation value of \(L_z\) in this state?
Problem
You have a charge distribution on the \(x\)-axis composed of two point charges: one with charge \(+3q\) located at \(x=-d\) and the other with charge \(-q\) located at \(x=+d\).
- Sketch the charge distribution.
- Write an expression for the volume charge density \(\rho (\vec{r})\) everywhere in space.
Suppose you have a definite function \(f(x)\) in mind and you already know its Fourier transform, i.e. you know how to do the integral \begin{equation} \tilde{f}(k)=\frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty}e^{-ikx}\, f(x)\, dx \end{equation} Find the Fourier transform of the shifted function \(f(x-x_0)\).
Instructor's Guide
Introduction
Students will need a short lecture giving the definition of the Fourier Transform \begin{equation} {\cal{F}}(f) =\tilde{f} (k)= \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} e^{-ikx}\, f(x)\, dx \end{equation}
Student Conversations
This example will feel very abstract to some students. It may be difficult for them to understand that the conditions of the problem state that the know both \(f(x)\) and \(\tilde{f}(k)\). This problem is about changing \(f\) slightly (by shifting its argument by \(x_0\)) and then asking how \(\tilde{f}\) changes, in response.Wrap-up
The result from this calculation underlies why it is possible to factor out the time dependence in the Fourier transform of a plane wave, Fourier Transform of a Plane Wave. Even though the problem is somewhat abstract, it is super important in applications for this reason.
Students compute probabilities and averages given a probability density in one dimension. This activity serves as a soft introduction to the particle in a box, introducing all the concepts that are needed.
Students learn about the geometric meaning of the amplitude and period parameters in the sine function. They also practice sketching the sum of two functions by hand.
In this small group activity, students calculate a (linear) function to represent the charge density on a one-dimensional rod from a description of the charge density in words.
This small group activity introduces students to constrained optimization problems. Students work in small groups to optimize a simple function on a given region. The whole class wrap-up discussion emphasizes the importance of the boundary.
In this small group activity, students determine various aspects of local points on an elliptic hill which is a function of two variables. The gradient is emphasized as a local quantity which points in the direction of greatest change at a point in the scalar field.
In this small group activity, students determine various aspects of local points on an elliptic hill which is a function of two variables. The gradient is emphasized as a local quantity which points in the direction of greatest change at a point in the scalar field.
The students are shown the graph of a function that is a superposition of three harmonic functions and asked to guess the harmonic terms of the Fourier series. Students then use prewritten Sage code to verify the coefficients from their guess. The program allows the students to enter functions of their own choice as well as the one that is preset.
Students are placed into small groups and asked to calculate the total differential of a function of two variables, each of which is in turn expressed in terms of two other variables.
Students work out heat and work for rectangular paths on \(pV\) and \(TS\) plots. This gives with computing heat and work, applying the First Law, and recognizing that internal energy is a state function, which cannot change after a cyclic process.
The students will set up a Styrofoam cup with heating element and a thermometer in it. They will measure the temperature as a function of time, and thus the energy transferred from the power supply, from which they compute changes in entropy.
Students use a pre-written Mathematica notebook or a Geogebra applet to explore how the shape of the effective potential function changes as the various parameters (angular momentum, force constant, reduced mass) are varied.
Students use prepared Sage code or a prepared Mathematica notebook to plot \(\sin\theta\) simultaneously with several terms of a power series expansion to judge how well the approximation fits. Students can alter the worksheet to change the number of terms in the expansion and even to change the function that is being considered. Students should have already calculated the coefficients for the power series expansion in a previous activity, Calculating Coefficients for a Power Series.
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.
Students calculate probabilities for energy, angular momentum, and position as a function of time for an initial state that is a linear combination of energy/angular momentum eigenstates for a particle confined to a ring written in bra-ket notation. This activity helps students build an understanding of when they can expect a quantity to depend on time and to give them more practice moving between representations.
Students, working in pairs, use their left arms to represent each component in a two-state quantum spin 1/2 system. Reinforces the idea that quantum states are complex valued vectors. Students make connections between Dirac, matrix, and Arms representation.
In this sequence of small whiteboard questions, students are shown the contour graph of a function of two variables and asked to find the derivative. They discover that, without a function to differentiate, they must instead think of the derivative as a ratio of small changes. This requires them to pick two nearby points. Which two?
Students calculate probabilities for a particle on a ring using three different notations: Dirac bra-ket, matrix, and wave function. After calculating the angular momentum and energy measurement probabilities, students compare their calculation methods for notation.