Past exam of the advanced Quantum Condensed Matter Physics course of the University of Cambridge

2016 paper: Public original PDF, physical page 783. Student-answer availability for this paper remains unverified.

Original solutions to the released 2016 1/AQC paper. ^[ref]

The absorptive dielectric function counts optically accessible electron-hole pairs. In the electric dipole approximation the perturbation is proportional to $\mathbf{a}\cdot \hat{\mathbf{p}}$, with $\mathbf{a}$ the light polarization and $\hat{\mathbf{p}}$ the electronic momentum operator. The squared matrix element between occupied valence and empty conduction states is the transition strength; it includes polarization and band-symmetry selection rules. The sums enumerate initial and final subbands. photon momentum is negligible on the crystal-momentum scale, so the transition is approximately vertical in the in-plane wavevector $\mathbf{k}$.

For each subband pair, $C$ is the constant-transition-energy contour in the two-dimensional Brillouin zone. Its line element is $dl$, and the reciprocal gradient weights the states satisfying energy conservation. Specifically, the level-set representation of the joint density of states isA small energy gradient means many nearly resonant states; a vanishing gradient can produce a Van Hove singularity or an edge. The displayed expression suppresses overall optical prefactors and assumes the momentum matrix element is approximately constant along each contour. More generally it stays inside that integral. Occupation factors $f_v-f_c$ are also implicit; at low temperature in an undoped well the relevant valence states are full and the conduction states empty. This is the independent-particle continuum contribution from Fermi's golden rule. Bound excitons add optical resonances below the corresponding free-pair thresholds.

The GaAs four-band model is rooted in bonding and spin-orbit coupling. Near the zone centre, the conduction edge is predominantly s-like and antibonding, while the upper valence manifold is predominantly p-like. Combining orbital angular momentum $l=1$ with spin $s=1/2$ splits the six p-derived states into a $j=3/2$ quartet and a lower $j=1/2$ doublet. The latter is the split-off valence band. Within the upper quartet, the $m_j=\pm 3/2$ and $m_j=\pm 1/2$ branches, defined relative to a propagation or quantization direction, have different curvatures and are called heavy holes and light holes. They coincide at the unstrained bulk zone centre but need not remain degenerate in a well. The conduction edge is a spin doublet. Thus “four bands” counts conduction, heavy-hole, light-hole and split-off branches, including their spin multiplicities; it does not count four individual states. Positive hole masses describe the downward valence curvature, whereas the conduction curvature is upward.

For the final exciton comparison, define the reduced mass ${\mu }_j$ and screened exciton Rydberg $R_j^{\ast }$ byThe resulting $({\mu }_j/m_e,R_j^{\ast })$ are approximately $(0.05948,4.937\mathrm{meV})$ for heavy holes, $(0.03646,3.027\mathrm{meV})$ for light holes and $(0.04631,3.844\mathrm{meV})$ for split-off holes. In the ideal two-dimensional exciton model, the first two optically active internal states lie $4R_j^{\ast }$ and $4R_j^{\ast }/9$ below their continuum edge. Their separation is $32R_j^{\ast }/9$, independent of a common threshold offset. The corresponding three-dimensional separation would be only $3R_j^{\ast }/4$. The listed peak spacings specified for this theoretical spectrum giveHere each pair comprises the stronger lower-energy peak and its weaker nearby companion, in the order associated with the continuum edges calculated in part (c). The approximately 18-meV and 11-meV pair splittings support the two-dimensional series and are incompatible with the three-dimensional series at the supplied masses and permittivity. Equivalently, for a two-dimensional pair with energies $E_1,E_2$, the inferred edge is $E_2+(E_2-E_1)/8$; the observed pair splittings predict ground-state bindings around 20 and 12 meV, close to $4R_j^{\ast }$.

Absolute agreement should not be overstated. For example, the stated square-well parameters give the first heavy-hole edge $1.56390$ eV and the first light-hole edge $1.59161$ eV. Their ideal two-dimensional peak predictions are approximately $(1.54415,1.56170)$ eV and $(1.57950,1.59027)$ eV, respectively. The listed peaks are lower by about 1–2 meV across the families; rounded parameters and this simple effective-mass model do not reproduce every absolute energy exactly. No gap, width or mass has been silently adjusted to fit them. The peak-spacing comparison remains useful because it avoids a common edge offset. The split-off transitions predicted by the four-band data are discussed separately in part (c), rather than assigned fictitious extra peaks.

At 300 K, $k_{B}T\simeq 25.9$ meV, whereas the ideal ground-state binding energies are only about $19.7$, $12.1$ and $15.4$ meV; excited states are still more weakly bound. None of these series is expected to remain a clearly resolved, robust bound-exciton peak at room temperature in this simple estimate. Thermal ionization and phonon broadening work against their visibility. This is a scale argument, not a theorem that the exciton population or all Coulomb-enhanced absorption becomes exactly zero: actual observability also requires temperature-dependent gaps, linewidths, screening and carrier density, which are not supplied.

Choose the growth direction as $z$ and the well as $0<z<L$. In the envelope function approximation, write a subband state as $A^{-1/2}{e}^{i{\mathbf{k}}_{\parallel }\cdot {\mathbf{r}}_{\parallel }}{F}_n(z)u_{j0}(\mathbf{r})$, where $u_{j0}$ is the band-edge cell-periodic Bloch factor. The infinite square well hasThe envelope has the same shape for electron and hole states: changing their effective masses changes the eigenenergies, but not these infinite-well eigenfunctions. With in-plane polarization, the leading interband momentum matrix element acts on the Bloch factors, while the slow envelopes giveorthonormality therefore yields the infinite-well interband selection rule $n^{\prime }=n$. A nonzero band-level matrix element is still required; the overlap rule does not replace polarization selection rules. Finite barriers can make electron and hole envelopes different, and interface asymmetry or band mixing can relax this ideal same-index rule.

For a fixed allowed band/subband pair, parabolic dispersions give $E_c(\mathbf{k})-E_v(\mathbf{k})=E_{j,n}^{\mathrm{edge}}+{\hslash }^2{k}_{\parallel }^2/(2{\mu }_j)$. If the photon energy exceeds the edge, the contour is a circle of radius $k_{\parallel }$ and the gradient magnitude is ${\hslash }^2{k}_{\parallel }/{\mu }_j$. ThusBelow threshold the contour is absent. With the conventional $d^{2}k/(2\pi {)}^2$ normalization, the two-dimensional joint density of states per area and independent channel isWith a locally energy-independent optical matrix element, each allowed pair therefore adds a constant step to the continuum dielectric function in the approximation of the supplied equation. Different hole families need not have equal step heights: reduced masses, spin multiplicities and Bloch-state matrix elements differ. Frequency-dependent overall optical prefactors, nonparabolicity, broadening and excitonic Coulomb enhancement can modify the literal step shape outside this approximation.

electron confinement raises the conduction edge and hole confinement lowers the corresponding valence edge, so their positive confinement energies add. With $L=120\mathring{\mathrm{A}}=12$ nm, the common bare-mass scale isWriting ${\eta }_e=0.067$ and ${\eta }_j=m_{h,j}^{\ast }/m_e$, the quantum-well interband thresholds areUse $E_{g,\mathrm{hh}}=E_{g,\mathrm{lh}}=1.52$ eV and $E_{g,\mathrm{so}}=1.86$ eV. The thresholds in the requested window, using the printed masses without fitting, areThe next heavy-hole and split-off thresholds are about $2.2224$ and $2.3674$ eV, outside the interval. Each threshold marks a continuum onset or shoulder, not the stronger bound-exciton peak below it. The six heavy- and light-hole onsets follow the six paired peak groups in the supplied diagram. Their first members lie roughly 20 and 12 meV below their respective onsets, as the common solution above explains.

A qualification is necessary for the split-off branch. Its p-like $j=1/2$ bloch states have nonzero $p_x$ matrix elements to the s-like conduction states, so in-plane polarization does not generally forbid them. For example, a $j=1/2$ state has a nonzero $p_x\pm i{p}_y$ orbital component with the appropriate spin, which couples to an s-like state through $p_x$. The same-envelope rule therefore allows the two split-off thresholds in the table. The first nearly coincides with the third heavy-hole edge, while the second would give an additional onset near 2.086 eV. The supplied figure and twelve peak energies appear to retain only the heavy- and light-hole series; they do not specify a split-off-series oscillator strength or identify its extra features. Omitting the split-off transitions would require an additional approximation beyond the four-band information and polarization given. Their expected ideal two-dimensional first-state peaks would be near 1.9010 and 2.0701 eV; these are predictions of the stated model, not additional listed data.

Start from electrons moving in the periodic equilibrium ionic potential. Displace ion $l$ by ${\mathbf{u}}_l$ and expand its contribution as $V(\mathbf{r}-{\mathbf{R}}_l-{\mathbf{u}}_l)\simeq V(\mathbf{r}-{\mathbf{R}}_l)-{\mathbf{u}}_l\cdot \nabla V(\mathbf{r}-{\mathbf{R}}_l)$. Quantizing the harmonic normal modes expresses ${\mathbf{u}}_l$ as a sum of creation and annihilation operators with displacement scale $\sqrt{\hslash /(2{\rho }_{\mathrm{ion}}V{\omega }_{\mathbf{q}})}$ for a single simple branch. Expanding the electronic field in bloch states then gives a band-projected Fröhlich HamiltonianHere $c,c^{†}$ obey canonical anticommutation relations and remove/create an electron with crystal momentum $\mathbf{k}$ and spin $\sigma$; $b,b^{†}$ are boson operators for lattice modes. The matrix element $M_{\mathbf{q}}$ contains the derivative of the ionic potential, phonon polarization, displacement normalization and electronic wavefunction overlaps. With phases chosen appropriately, $M_{-\mathbf{q}}=M_{\mathbf{q}}^{\ast }$ ensures a Hermitian Hamiltonian. The displayed one-band, one-branch notation suppresses additional branch and band labels. “Fröhlich Hamiltonian” here means the linear electron-phonon model used in the question; it is not a requirement to use specifically the polar-optical-phonon coupling of a polar semiconductor. A harmonic lattice, small displacements and a weak enough coupling to organize perturbation theory are the underlying approximations.

The Fröhlich canonical transformation removes real linear phonon operators from the electron-only effective problem to first order. Choose an anti-Hermitian $S$ with $[S,H_0]=-V_{e\text{-}ph}$. For $A_{\mathbf{kq}\sigma }=c_{\mathbf{k}+\mathbf{q},\sigma }^{†}{c}_{\mathbf{k}\sigma }$ and ${\Delta }_{\mathbf{kq}}={ϵ}_{\mathbf{k}+\mathbf{q}}-{ϵ}_{\mathbf{k}}$, one possible generator isIts Hermitian conjugate is its negative after relabelling momenta. The Baker--Campbell--Hausdorff formula givesThe commutator describes virtual emission and absorption of a phonon. normal ordering yields an electron-electron term as well as one-electron energy renormalizations and phonon terms. Retain the former, absorb single-particle shifts into ${ϵ}_{\mathbf{k}}$ and eliminate the explicit phonon sector. For elastic pair scattering the two virtual time orderings have the combined amplitudeUsing the usual factor one-half to count unordered pairs gives the coefficient in the printed Hamiltonian, with $\Delta ={ϵ}_{\mathbf{k}}-{ϵ}_{\mathbf{k}-\mathbf{q}}$. The two outgoing momenta are $\mathbf{k}-\mathbf{q}$ and ${\mathbf{k}}^{\prime }+\mathbf{q}$, so total momentum and the spins are conserved. A fully off-shell effective operator can instead use symmetrized denominators for both external legs; the displayed expression is the usual elastic-scattering reduction. Eliminating nearly resonant real phonons requires more care than this second-order construction.

For $|\Delta |<{\Omega }_{\mathbf{q}}$ the denominator is negative: virtual phonon exchange is attractive in this low-energy channel. Physically, an electron changes the ionic density and a second electron can gain energy by responding to the delayed distortion. This is a retarded interaction, not an unscreened instantaneous attraction at every energy transfer. At large transfer the sign changes, and screened Coulomb repulsion must also be included when deciding the net low-energy pairing strength. The connection between lattice dynamics and superconductivity is supported by the isotope effect, but the mathematical attraction alone does not determine a measured transition temperature.^[ref]

For the pair calculation, make the model restrictions explicit. The filled Fermi sea is inert, and $|0⟩$ denotes that sea with no added particles, not an empty electronic vacuum. Project the interaction onto the previously empty shell $0\le {\xi }_{\mathbf{k}}\le D$, where ${\xi }_{\mathbf{k}}={\hslash }^2{k}^2/(2m)-E_F$ and $D=\hslash {\omega }_D$. Measure the added-pair energy relative to $2E_F$; this is the negative $E$ in the final formula. If absolute one-electron energies are retained in the preceding eigenvalue equation, its eigenvalue must instead be $E_{\mathrm{abs}}=2E_F+E$. A shell-restricted state would not in general be an exact eigenstate of the completely unprojected, all-momentum contact Hamiltonian.

Write $g={\gamma }^2>0$ and $|\mathbf{k}⟩=c_{-\mathbf{k}\downarrow }^{†}{c}_{\mathbf{k}\uparrow }^{†}|0⟩$. Within the zero-total-momentum pair subspace the two opposite-spin terms cancel the interaction's factor one-half, while same-spin s-wave terms cancel by fermionic antisymmetry. Thus the projected Hamiltonian acts asFor $|\Psi ⟩={\sum }_{⟨\mathbf{k}⟩}{\phi }_{\mathbf{k}}|\mathbf{k}⟩$ and $C={\sum }_{⟨\mathbf{k}⟩}{\phi }_{\mathbf{k}}$, matching each basis coefficient givesThis is the stated coefficient formula in the consistent Fermi-level energy convention. Its even dependence on $\mathbf{k}$ produces an even orbital state with a spin singlet. The constant $C$ is fixed up to an overall phase by normalization; for a bound state, $|C{|}^{-2}=g^2{\sum }_{⟨\mathbf{k}⟩}(E-2{\xi }_{\mathbf{k}}{)}^{-2}$. It cannot be set to zero for the nontrivial bound solution.

Put $E=-B$ with $B>0$ and impose the sum defining $C$:Here $\rho (E_F)$ counts orbital states for one spin species, as appropriate to the unrestricted $\mathbf{k}$ sum in this pair basis, and is approximately constant over the thin shell. If a two-spin total density of states is used, it must be divided by two in this integral. State-counting and the volume convention in $g$ and $\rho$ must be consistent. With $\xi =g\rho (E_F)$ as the paper's dimensionless coupling, inversion yieldsThe first approximation is the constant-density-of-states, restricted-shell model; the final exponential additionally assumes weak attraction. The binding energy is nonanalytic at zero coupling, so a finite-order power expansion in $g$ cannot produce it.

Cooper's problem asks whether two electrons added above an inert Fermi sea bind despite Pauli exclusion of lower states. An arbitrarily weak net attraction binds a zero-momentum singlet below the free-pair threshold, because the near-Fermi-surface integral diverges logarithmically as $B\downarrow 0$. Its wavefunction draws on many time-reversed momentum states and becomes dominated by a narrow region near the Fermi surface at weak coupling; the corresponding pair is spatially extended. This establishes a Cooper instability of the unpaired reference state.^[ref] The full superconducting ground state needs a coherent many-electron treatment in which overlapping pairs and the reorganized sea are included. BCS theory supplies that step; the two-electron binding scale alone is neither a calculation of the many-body gap nor a derivation of zero resistance.^[ref]

A Bravais lattice is the translation set $\Lambda =\{n_1{\mathbf{a}}_1+n_2{\mathbf{a}}_2+n_3{\mathbf{a}}_3:n_i\in \mathbb{Z}\}$ generated by primitive lattice translation vectors. Translating by any member of $\Lambda$ leaves its surroundings equivalent. A crystal structure requires both this lattice and an atomic basis attached to every lattice point. The basis can lower the symmetry of the bare point lattice; a familiar-looking geometrical arrangement of atoms is not necessarily itself a Bravais lattice. For example, the diamond cubic structure is a face-centred cubic lattice with an identical-atom basis at fractional positions $(0,0,0)$ and $(1/4,1/4,1/4)$ in conventional cubic coordinates. The displacement between those two basis atoms is not itself a lattice translation.

The seven three-dimensional crystal systems are triclinic, monoclinic, orthorhombic, tetragonal, trigonal, hexagonal and cubic; the distinct translation types give 14 Bravais lattices. A conventional cell can contain several lattice points, while a primitive unit cell contains one. Simple cubic primitive lattice translation vectors can be chosen along the cube edges; for face-centred cubic they can be $(a/2)(0,1,1)$, $(a/2)(1,0,1)$ and $(a/2)(1,1,0)$. Their primitive volume is $a^3/4$, explaining the four lattice points per conventional face-centred cell.

A crystallographic point group comprises the orthogonal parts of the crystal's spatial symmetry operations. It constrains possible rotations, reflections, inversion and their combinations, while forgetting translations. Lattice preservation restricts proper rotational orders in three dimensions to $1,2,3,4,6$; a fivefold axis cannot be a symmetry of a periodic three-dimensional lattice. One quick argument restricts a rotation's integral lattice representation through $1+2\cos \theta$, an integer trace, producing these orders. There are 32 crystallographic point-group types. Examples are the cubic full octahedral group $O_h$ with Hermann–Mauguin symbol $m\bar{3}m$, the tetrahedral class $T_d$ with symbol $\bar{4}3m$, and the polar tetragonal class $C_{4v}$ with symbol $4mm$.

The space group retains the full isometry. In Seitz notation $\{R\mid \mathbf{t}\}$ acts by $\mathbf{x}\mapsto R\mathbf{x}+\mathbf{t}$, and multiplication isIts translation subgroup is normal and the quotient is the crystallographic point group. A screw axis combines rotation with translation along the axis; for example, applying a $2_1$ screw twice gives one whole lattice translation. A glide plane combines reflection with a translation parallel to the plane; applying the glide twice likewise gives a translation. These operations can preserve the complete crystal without their rotation or reflection alone preserving it.

A symmorphic space group admits an origin and representatives for which the point-group operations have no translation component. A nonsymmorphic space group has essential fractional translations that cannot all be removed by moving the origin. There are 230 three-dimensional space-group types, distinguished more finely than their point groups. These counts refer to ordinary, nonmagnetic periodic space groups; magnetic and aperiodic generalizations are different classifications.^[ref]

For concrete comparisons, a monatomic simple-cubic crystal has space group $Pm\bar{3}m$ and point group $m\bar{3}m$. Rock salt has a face-centred cubic lattice, a two-species basis and symmorphic space group $Fm\bar{3}m$. Diamond has the same cubic point group and translation-lattice type, but its basis leads to nonsymmorphic $Fd\bar{3}m$. Replacing the two diamond basis sites by distinct species gives the zinc blende structure, exemplified by GaAs, with $F\bar{4}3m$ and point group $\bar{4}3m$; inversion that would interchange the species is no longer a symmetry. The lattice specifies translations, the point group specifies the linear symmetry classes, and the space group specifies their complete spatial implementation. These examples show why none alone uniquely describes an atomic structure.

A quantum dot confines electronic motion in all three directions. A semiconductor dot can be thin enough that only one transverse subband is occupied, with discrete lateral states controlled by gates. In a rectangular hard-wall model,A smooth circular lateral confinement is often better approximated by $m^{\ast }{\omega }_0^2{r}^2/2$, giving ${ϵ}_{n_x{n}_y}=\hslash {\omega }_0(n_x+n_y+1)$ above the transverse ground energy. Shell $q=n_x+n_y$ contains $q+1$ orbitals, twice that many electrons before interactions split the spin degeneracy. Hence closed-shell electron numbers are $2,6,12,20,\dots$. These are single-particle quantum numbers; an interacting dot's states are many-electron eigenstates, not simply one orbital occupied by a whole charged island.

A perpendicular magnetic field changes the harmonic spectrum to the Fock-Darwin spectrum. For an electron of charge $-e$, take ${\omega }_c=eB/m^{\ast }>0$ and $L_z=\ell \hslash$ in the usual symmetric gauge. ThenReversing the definition of the angular label changes the displayed sign of its orbital term, not the set of energies. The Zeeman term lifts spin degeneracy, while exchange can favor parallel-spin filling in partially filled shells. These effects are distinguished by how peaks evolve with field rather than by identifying every irregular peak spacing with one-particle confinement.

For electrostatics, let $C_{\Sigma }=C_s+C_d+C_g$ and $Q_g=C_s{V}_s+C_d{V}_d+C_g{V}_g+Q_0$, where $Q_0$ represents a fixed offset charge in the chosen convention. The constant interaction model givesThe charging energy is conventionally $E_C=e^2/(2C_{\Sigma })$, whereas adjacent charge-state potentials differ byThis addition energy separates a capacitive contribution from orbital spacing. Exchange, nonuniform screening and correlations are corrections to the constant interaction model; a change in dot shape or effective capacitance can also alter spacings.

Connect the island weakly to source and drain reservoirs through tunnel barriers. sequential tunnelling is suppressed in a charge-stable region when neither adding nor removing an electron is energetically accessible. Sweeping the gate brings a charge electrochemical potential through the reservoir Fermi level, giving a Coulomb-blockade conductance peak. The gate lever arm is ${\alpha }_g=C_g/C_{\Sigma }$, so $\Delta V_g=E_{\mathrm{add}}/(e{\alpha }_g)$ locally if the capacitances and level spectrum do not vary appreciably. At finite source-drain bias, alignment with either reservoir defines the boundaries of a Coulomb diamond. Its slopes constrain capacitance ratios, and lines parallel to its edges can reveal excited-state transport. The conversion from a diamond's bias height to an energy depends on bias convention; a full-height/half-height factor must not be assumed without declaring the circuit.

A clear charging spectrum requires $k_{B}T$ and tunnel broadening $\hslash \Gamma$ smaller than the relevant addition energy; resolving individual orbital levels additionally requires them smaller than the orbital spacing. Tunnel resistances large compared with $h/e^2$ suppress strong charge fluctuations in the usual weak-junction regime. Blockade is not an assertion of mathematically zero current: higher-order cotunnelling can remain, and strong coupling or low-temperature spin correlations can invalidate the sequential model. An experimental check therefore combines gate sweeps with bias spectroscopy, temperature dependence and barrier tuning to separate level structure from broadening or uncontrolled charge traps.

In a primary vertical-dot experiment, gate-controlled current oscillations were converted into addition energies using finite-bias diamond spectroscopy; varying magnetic field then followed shell and spin filling. Reproducing enhanced addition energies at electron numbers 2 and 6 in several devices strengthened the shell assignment. The electron temperature was independently estimated above the refrigerator base temperature, an essential distinction when judging spectral resolution. Field-dependent reordering also revealed exchange effects beyond a purely spin-degenerate orbital model.^[ref]

Use a three-dimensional electron gas with a uniform positive background, the jellium model. The Hartree-Fock method minimizes the expectation of the interacting Hamiltonian over normalized Slater determinants. Antisymmetry is exact within that trial class, but correlations requiring a superposition of determinants are omitted. For orthonormal occupied spin orbitals ${\varphi }_i$, the energy iswhere $J_{ij}$ is the direct Coulomb integral and $K_{ij}$ the exchange integral; the spin overlap makes exchange vanish between opposite-spin orbitals. The self-interaction terms $J_{ii}$ and $K_{ii}$ cancel. Variation gives a self-consistent Fock operator $\hat{f}=h_0+{\sum }_j({\hat{J}}_j-{\hat{K}}_j)$; exchange is nonlocal. Its eigenfunctions define occupied orbitals of a stationary determinant, and the determinant must be compared with competing self-consistent solutions to select a variational ground state.

The eigenvalues are orbital Lagrange multipliers. Their sum counts each interaction twice:Frozen-orbital removal energies may be estimated using Koopmans' theorem, but measured quasiparticle energies can contain relaxation and correlations absent from this approximation. Unoccupied orbital energies are not automatically an exact excitation spectrum. In particular, unscreened exchange in a uniform three-dimensional gas produces an unphysical logarithmic slope singularity at its Fermi surface; a one-determinant mean field does not reproduce the interacting Fermi liquid's screened response.

For a spatially uniform plane-wave determinant, electron-background and background-background terms cancel the uniform direct Coulomb contribution. Retain the background when making this cancellation: an isolated infinite electron charge distribution would have a divergent direct energy. The surviving contributions per particle are kinetic energy and negative same-spin exchange energy. Let $(3/(4\pi n){)}^{1/3}=r_s{a}_B^{\ast }$ define the Dimensionless Wigner-Seitz radius, with $a_B^{\ast }=4\pi \epsilon {\hslash }^2/(m^{\ast }{e}^2)$, $E_H^{\ast }={\hslash }^2/(m^{\ast }{a}_B^{\ast 2})$ and $\zeta =(n_{\uparrow }-n_{\downarrow })/n$. ThenFor the unpolarized state this is approximately $1.105/r_s^2-0.4582/r_s$. The kinetic contribution penalizes spin polarization because it enlarges an occupied Fermi sphere; exchange favors polarization because equal-spin occupation gains more exchange energy. As density falls, the $r_s^{-1}$ exchange scale overtakes the $r_s^{-2}$ kinetic scale. The uniform fully polarized and unpolarized energies cross atThis is a variational comparison within the homogeneous Hartree-Fock electron gas, not an exact physical ferromagnetic phase boundary. It demonstrates how Pauli kinetic energy and exchange compete and why the allowed trial states matter.

Allowing spatial symmetry breaking changes the result: the Overhauser instability lowers the paramagnetic determinant's Hartree-Fock energy through a spin density wave, even in the high-density limit. A calculation restricted to uniform Fermi spheres therefore does not give the unrestricted Hartree-Fock ground state.^[ref] At still lower density, localization and a Wigner crystal compete with fluid states. electron-gas correlation energy beyond a single determinant is important in assessing their competition. Quantum Monte Carlo provides a concrete route to comparing correlated fluid and crystal energies; the original stochastic calculation demonstrates why exchange-only energy differences are insufficient.^[ref]

Hartree-Fock is a useful variational baseline with exact exchange, not a controlled complete low-density phase diagram. At high density the kinetic scale dominates, but even a small symmetry-breaking energy gain can invalidate a restricted determinant as the exact Hartree-Fock minimizer. At low density neglected correlations, screening and localization can change both transition locations and the ordering of candidate phases. A numerical phase boundary should therefore state dimension, spin constraints, background model and whether correlation and spatially modulated states were included.

The Hartree-Fock method minimizes the expectation of an interacting-electron Hamiltonian over normalized Slater determinants of orthonormal one-particle spin orbitals. It is a variational, self-consistent mean-field approximation. The determinant enforces fermionic antisymmetry exactly within that class, while each orbital responds to the mean charge density and a nonlocal exchange operator generated by the occupied orbitals.

For an electron gas, this accounts for kinetic pressure, the average electrostatic field, same-spin exchange and the associated exchange hole, and possible spin polarization. It does not account for general dynamical electron correlation, such as the correlated avoidance of opposite-spin electrons, and cannot capture all screening or collective fluctuations merely by optimizing one determinant. Uniform plane waves are an additional restriction, not the definition of Hartree-Fock itself.

In the energy functional, the one-body term is the sum of electron kinetic energy and ionic-potential energy. The direct two-body term is the Hartree energy: the classical interaction of the mean electron density with itself. Its factor $1/2$ avoids counting every pair twice. The negative Fock term comes from exchanging the two indistinguishable electron labels in the antisymmetric state. For a spin-independent interaction, spin orthogonality makes this exchange act only between equal spin components. For an orbital interacting with itself, direct and exchange contributions cancel exactly; the exchange term removes the spurious one-electron self-interaction. Ionic background terms and their prescribed cancellation must be included consistently with the direct term.

Use normalized plane-wave spin orbitals $V^{-1/2}{e}^{i\mathbf{k}\cdot \mathbf{r}}{\chi }_{\sigma }$ and fill a separate Fermi sphere for each spin species. Each spin state counts once, soIts kinetic energy per volume isThe spin-independent contact interaction has uncancelled direct energy $E_H/V=\gamma n^2/2$, with $n=n_{\uparrow }+n_{\downarrow }$. By the imposed background assumption this direct contribution cancels. For the Fock term the contact delta function sets the two spatial positions equal; its spin factors then select equal spins. The local one-body density matrix for each spin has diagonal value $n_{\sigma }$, givingCombining these results gives the contact-interaction Hartree-Fock spin energyHere $\gamma$ has units of energy times volume. The result is the model with the stated background subtraction; an isolated contact-interacting gas without that subtraction retains its direct contribution.

For the unpolarized gas write $x=n/2$. The positive direct Hartree energy before background cancellation and the kinetic energy areTheir ratio is thereforeThis compares the separate direct and kinetic terms. The net Hartree-plus-background contribution after the stipulated cancellation is zero, so it would be incorrect to interpret $r_s$ as that residual contribution divided by $T$. The exchange-energy magnitude in the unpolarized state is $E_H/2$, not $E_H$.

At fixed total density, the unpolarized and fully polarized candidate energies arePolarization increases the kinetic energy because one spin species must occupy a larger Fermi sphere, but gains negative exchange energy. The fully polarized candidate wins this endpoint comparison ifUsing the defined $r_s$ givesThus the specified two-state comparison indeed predicts a polarized state below the unpolarized state's energy for $r_s>1.175$.

The endpoint and continuous polarization thresholds of a contact electron gas are different. To check the actual onset within the same uniform Hartree-Fock model, allow $P=(n_{\uparrow }-n_{\downarrow })/n$ to vary continuously. In units $2A(n/2{)}^{5/3}$,The unpolarized state loses local stability already at $r_s=10/9\simeq 1.111$. Full polarization becomes the minimizing endpoint only when $f^{\prime }(1)\le 0$, or $r_s\ge (5/6)2^{2/3}\simeq 1.323$. Between those values the optimum is partly polarized. Consequently $1.175$ is the requested endpoint crossing, not the exact onset of spontaneous polarization after all spin populations are allowed.

The central limitation is that a Hartree-Fock electron gas includes exchange but omits the correlation energy that is not representable by one determinant. Opposite-spin electrons can lower their repulsion by avoiding one another in a correlated unpolarized state. Hartree-Fock misses that avoidance while retaining the exchange advantage of aligned spins, so it need not place a magnetic instability at a reliable interaction strength. The relevant energy differences between competing spin states can be much smaller than their total energies. Correlation calculations of the electron gas explicitly address this competition.^[ref]

There are further model restrictions. A constant contact coupling discards the momentum dependence of the Coulomb interaction and the resulting nontrivial screening. Real materials have a lattice, band structure and material-dependent effective interactions. Spatial magnetic textures or other competing states are excluded by a comparison restricted to uniform Fermi spheres. Even before those changes, part c's two endpoints do not find the continuous polarization onset.

The symbol $r_s$ also needs care. Here it is an energy ratio proportional to $\gamma n^{1/3}$. The usual Dimensionless Wigner-Seitz radius of a three-dimensional Coulomb electron gas is $r_s^{\mathrm{WS}}=[3/(4\pi n){]}^{1/3}/a_B$, proportional to $n^{-1/3}$. They are not the same density parameter. A quoted experimental value of the latter cannot be compared numerically to $1.175$ without a model-specific mapping. The contact interaction, missing correlations, restricted variational class and parameter convention prevent the number $1.175$ from being a quantitative experimental prediction. The paper supplies no particular experimental threshold, so none is invented here.

Landau quasiparticles describe the long-lived low-energy excitations of an interacting normal Fermi liquid. Starting from a free Fermi gas, imagine switching on interactions without changing the phase or destroying the one-to-one labelling of low-energy states. The resulting excitation carries the quantum numbers of a dressed particle but has a renormalized energy and interactions with other excitations. Its microscopic wavefunction need not resemble a bare electron state. The occupied momenta define a Fermi surface, and the quasiparticle energy measured from the chemical potential vanishes on that surface.

The Landau quasiparticle energy functional expresses the change near equilibrium in terms of occupation changes:Differentiation gives the energy of one quasiparticle in the distorted distribution. This is a phenomenological expansion around the true interacting ground state, not the assertion that its total energy equals that of independent particles. For an isotropic system, spin-symmetric and spin-antisymmetric components of $f$ are expanded in Legendre polynomials. Multiplication by the appropriate Fermi-surface density of states defines dimensionless Landau parameters $F_l^s,F_l^a$. The original formulation relates a change in the distribution to the energy shift.^[ref]

Pauli blocking makes the low-energy excitations sharp. The phase space for a quasiparticle to scatter while creating another particle–hole pair shrinks as the Fermi surface is approached, giving the regular three-dimensional low-energy Fermi-liquid quasiparticle lifetimeThus at zero temperature its width is small compared with its excitation energy. The single-particle spectral function has a pole with quasiparticle residue $Z$, generally less than one, and an incoherent remainder. In a conventional Fermi liquid the occupation distribution has a discontinuity of height $Z$ at the Fermi surface, rather than a unit free-particle step.

For an isotropic three-dimensional system, define $N^{\ast }(0)$ to include both spins. The quasiparticle density of states determines the low-temperature electronic heat capacity, $C_V/V=({\pi }^2/3)k_B^2{N}^{\ast }(0)T$. In the usual Landau-parameter normalization, the spin response is modified by $1/(1+F_0^a)$ and the compressibility by $1/(1+F_0^s)$. In a Galilean-invariant continuum, $m^{\ast }/m=1+F_1^s/3$; that particular relation must not be imposed unchanged on a crystal band. A denominator approaching zero signals a possible instability of the assumed normal state.

The same interaction functional and a kinetic equation describe collective zero sound when the excitation frequency exceeds the collision rate. Ordinary first sound is the collision-dominated hydrodynamic limit. The theory organizes thermodynamics, response and transport through a small set of measured effective parameters while explaining why interactions coexist with long-lived fermionic excitations. It is not applicable without qualification in phases where the Fermi-liquid quasiparticle pole disappears or a gap opens over the relevant Fermi surface.

A Wannier-Mott exciton is a bound conduction-band electron and valence-band hole whose envelope extends over many lattice cells. That large radius permits an effective mass description and dielectric screening of the electron–hole attraction. With electron and hole masses $m_e^{\ast },m_h^{\ast }$, define reduced mass $\mu =m_e^{\ast }{m}_h^{\ast }/(m_e^{\ast }+m_h^{\ast })$ and center-of-mass mass $M=m_e^{\ast }+m_h^{\ast }$. The bulk relative-motion Hamiltonian is hydrogenic with an appropriate relative permittivity ${ϵ}_r$:A semiconductor quantum well confines the electron and hole perpendicular to the plane while leaving approximately free in-plane center-of-mass motion. If the well width is small compared with the bulk exciton Bohr radius, their relative motion becomes quasi-two-dimensional. The ideal zero-width, same-dielectric limit isIts ground-state binding is $4R_X$, compared with $R_X$ in three dimensions, and the ground envelope decays on the scale $a_X/2$. Stronger binding follows because confinement prevents the pair from separating in the third direction.

For a finite-width quantum-well exciton, average the Coulomb attraction over the normalized electron and hole confinement envelopes:This smooths the strict planar singularity and changes the binding. The ideal factor four is a useful reference, not a universal value for every real well. Finite barriers, dielectric mismatch, several confinement subbands and nonparabolic valence bands can require a coupled treatment.

Optical absorption creates a nearly zero in-plane center-of-mass wavevector exciton. For the lowest confinement subbands, its line lies approximately atThe discrete bound-state peaks occur below the free electron–hole continuum. Dipole-allowed envelope states have nonzero relative wavefunction at zero separation; the absorption strength also contains the confinement overlap $|\int {\chi }_e^{\ast }(z){\chi }_h(z)dz{|}^2$. Confinement can increase the overlap and the binding, making excitonic optical structure more prominent. A perpendicular electric field shifts and separates the confinement envelopes, changing the transition energy and usually reducing their optical overlap: this is the quantum-confined Stark effect.

Quantum wells tune both the continuum threshold and the bound-state energy; an optical line shift is not automatically a change in binding alone. The hydrogenic envelope model requires a Wannier-sized state, rather than an exciton localized to one lattice cell.

Experimental tests of BCS theory compare several properties of the same superconducting state. In the conventional weak-coupling setting, phonon-mediated attraction pairs opposite momenta and spins, and the BCS quasiparticle spectrum is $E_k=\sqrt{{\xi }_k^2+{\Delta }^2(T)}$. Agreement with a single energy scale and its temperature dependence is more discriminating than observing zero resistance alone.

The superconducting isotope effect tests the lattice's role. Replacing one isotope by another changes ionic masses while leaving the chemical valence nearly unchanged. Phonon frequencies scale approximately as $M^{-1/2}$, so a weak-coupling phonon mechanism predicts $T_c\propto M^{-\alpha }$ with $\alpha \simeq 1/2$ when other parameters stay fixed. The shift was measured in mercury before the microscopic pairing theory was developed.^[ref] Deviations from $1/2$ can result from Coulomb, anharmonic and strong-coupling effects; the isotope observation supports lattice-mediated pairing without uniquely proving every assumption of the simplest model.

Thermodynamic measurements test the gapped excitations and condensation energy. An isotropic fully gapped superconductor has exponentially small low-temperature electronic heat capacity, rather than the normal linear law. At the transition the weak-coupling BCS thermodynamic ratios areHere $C_n$ means the electronic contribution after subtracting the lattice background. The gap closes at $T_c$, and the same framework predicts the temperature dependence of penetration depth and thermodynamic critical field. These quantities were quantitative comparisons in the original theory.^[ref]