The Aharonov-Bohm effect

A charged particle confined to a 1D ring of radius \(R\) threads a magnetic flux \(\Phi\) through its center; even though the field \(B=0\) everywhere the particle can actually be, the eigenspectrum

\[E_n(\Phi) = \frac{\hbar^2}{2mR^2}\left(n - \frac{\Phi}{\Phi_0}\right)^2, \qquad \Phi_0 = \frac{2\pi\hbar}{q},\]

still shifts periodically with \(\Phi/\Phi_0\) – the Aharonov-Bohm effect, a purely topological/boundary-condition consequence of the vector potential.

import matplotlib.pyplot as plt
import numpy as np

from physicskit.quantum.chapters.entanglement import AharonovBohmRing

The ring’s flux-periodic energy spectrum

The ring’s energy levels \(E_n(\Phi)\) for angular-momentum quantum numbers \(n=-2,\dots,2\), plotted against the enclosed flux in units of the flux quantum \(\Phi_0\) – each parabola is centered on the \(n\) that minimizes the shifted quantum number, giving the sawtooth-like periodic ground-state energy.

ring = AharonovBohmRing(R=1.0)
Phi_over_Phi0 = np.linspace(-2, 2, 400)

fig, ax2 = plt.subplots(figsize=(7, 4.5))
for n in range(-2, 3):
    E_n = [ring.energy(n, f * ring.flux_quantum) for f in Phi_over_Phi0]
    ax2.plot(Phi_over_Phi0, E_n, label=f"n={n}")
ax2.set_xlabel(r"$\Phi / \Phi_0$")
ax2.set_ylabel("E_n(Phi)")
ax2.set_title("Aharonov-Bohm ring: flux-periodic energy spectrum")
ax2.legend(fontsize=7, ncol=2)
fig.tight_layout()

print(f"ground-state energy at Phi=0: {ring.energy(0, 0.0):.6f}")
print(f"ground-state energy at Phi=Phi0/2: {min(ring.energy(n, 0.5 * ring.flux_quantum) for n in range(-2, 3)):.6f}")

ground-state energy at Phi=0: 0.000000
ground-state energy at Phi=Phi0/2: 0.125000