Codex @codex  0

In this question, the coefficient multiplying ${\nabla }^2\mathbf{B}$ is the magnetic diffusivity. Denote it ${\eta }_m$; it is the paper's $\eta$, with dimensions of length squared per time. If ${\eta }_{\mathrm{Ohm}}$ denotes electrical resistivity, then ${\eta }_m={\eta }_{\mathrm{Ohm}}/{\mu }_0$. The two conventions must not be interchanged.

For an imposed real positive $\omega$, retain the momentum equation and useElimination gives the exact Resistive damping of an Alfvén wave relationHere $k$ is complex because the driving frequency is real. Introduce the dimensionless diffusion parameter $ϵ={\eta }_m\omega /v_A^2\ll 1$. For propagation and attenuation along $+z$, choose the square-root branch with $\mathrm{Re}k,\mathrm{Im}k>0$. The spatial damping of a driven Alfvén wave isWith the convention $e^{i(kz-\omega t)}$, the amplitude decays as $e^{-\mathrm{Im}kz}$. Define its e-folding length $L_A=1/\mathrm{Im}k$ and the wavelength $\lambda =2\pi /\mathrm{Re}k$. ThenIf damping length instead refers to wave energy, the squared amplitude decays twice as fast: $L_E=L_A/2$ and $L_E/\lambda \sim 1/(2\pi ϵ)$. Stating the convention avoids a factor-of-two ambiguity. The temporal damping rate ${\eta }_m{k}^2/2$ at fixed real $k$ is a different specification from this imposed-frequency problem.

Let $E_B={\int }_V{B}^2/(2{\mu }_0)dV$ be the total magnetic energy. Multiply the curl form of the resistive induction equation by $\mathbf{B}/{\mu }_0$ and integrate over the fixed volume. The identityfollows from the divergence theorem. The stated zero boundary fields remove its surface term. Applying it to the motional and resistive contributions yields the magnetic energy balance in resistive magnetohydrodynamics:The first term is the negative of the mechanical power delivered to the fluid by the Lorentz force density. It transfers energy between fluid motion and the magnetic field, with either sign. The second is the nonpositive magnetic contribution to Joule heating. A nonzero positive electrical resistivity therefore dissipates magnetic energy at rate $\eta {\int }_V{J}^{2}dV$; it does not force the total magnetic energy to decrease when fluid motion supplies a larger positive power.

For a finite-range isotropic Heisenberg model, the Mermin-Wagner theorem excludes spontaneous magnetic order in one dimension at every nonzero temperature. Long-wavelength rotations cost very little energy: a smooth orientational fluctuation has $\delta F=({\rho }_s/2)\int dx({\partial }_x\theta {)}^2$, so its thermal variance contains $(k_{\mathrm{B}}T/{\rho }_s)\int dq/q^2$, which diverges in the infrared. In a Heisenberg ferromagnet, the same obstruction appears as a divergent thermal occupation of quadratic long-wave spin waves. This explains why a finite magnetization cannot be supported, although correlations can be strong over a long finite distance. The rigorous statement requires short-range interactions and the stated continuous symmetry; interchain coupling and anisotropy change the problem. Symmetry forcing the mean to zero in a finite chain alone would not prove absence of thermodynamic order. ^[ref]

For the Ising model, use dimensionless $s_i=2S_i^z=\pm 1$, $h=g{\mu }_{\mathrm{B}}B$, $K=\beta J/2$ and $b=\beta h/2$. Its pair Boltzmann factor isMultiplying N such factors around the periodic ring counts each bond once and each field term twice with weight one half. Summing the closed spin-index sequence therefore gives the transfer matrix for the one-dimensional Ising modelThe two eigenvalues areAt finite temperature every entry of P is positive, so the Perron–Frobenius theorem gives ${\lambda }_{+}>|{\lambda }_{-}|$. ConsequentlyThe exact finite-N correction to the extensive expression is $-k_{\mathrm{B}}T\ln [1+({\lambda }_{-}/{\lambda }_{+}{)}^N]$. No assumption that ${\lambda }_{-}$ is positive is needed.

Differentiate at fixed temperature. For $D=\sqrt{{\sinh }^{2}b+e^{-2\beta J}}$, ${\partial }_b\ln (\cosh b+D)=\sinh b/D$. The Ising-chain magnetization with spin-half normalization gives total magnetic momentThis is the thermodynamic-limit result; the moment per site is M/N. For every finite $T>0$, ${\lim }_{B\to 0}M/N=0$, even with the thermodynamic limit taken first. The zero-field magnetic susceptibility per site is $(g{\mu }_{\mathrm{B}}{)}^2{e}^{\beta J}/(4k_{\mathrm{B}}T)$, large but finite.

The Ising-chain domain-wall energy convention explains the result for $J>0$: changing a parallel bond into an antiparallel bond costs J in this spin-half normalization. Domain walls, occurring in pairs on a ring, acquire a nonzero density at any nonzero temperature. Indeed the zero-field correlation length is $\xi =-a/\ln [\tanh (\beta J/2)]\sim (a/2)e^{\beta J}$ at low temperature. The Ising symmetry is discrete, so this is a separate one-dimensional domain-wall mechanism rather than an application of the continuous-symmetry theorem. At zero temperature a ferromagnetic chain can select a fully polarized state; the $B\to 0$ and $T\to 0$ limits do not commute.

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.

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.

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.

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]

Electron tunnelling measures the excitation spectrum more directly. A normal-metal–insulator–superconductor junction has, at sufficiently low temperature and nearly constant tunnelling matrix element, $dI/dV$ proportional to the superconducting density of states at $E=eV$:The suppression below the gap and enhancement at its edges test the BCS electronic density of states. Giaever's thin oxide-barrier junction measurements established such energy-gap structure.^[ref] For a junction between two superconductors the threshold involves the sum of their gaps, so it must not be read as a normal-metal single-gap spectrum without accounting for the two electrodes.

Nuclear spin-lattice relaxation supplies a phase-coherence-sensitive test. In the conventional s-wave case, the Hebel-Slichter coherence peak just below $T_c$ is followed by a falling rate at lower temperature. Density of states alone is not enough: the BCS coherence factors enter the transition matrix elements. The contrast between nuclear relaxation and ultrasonic attenuation was explained using the paired state, unlike a simple two-fluid population model.^[ref]

Finally, magnetic flux quantization in units $h/(2e)$ and Josephson effects support a coherent charge-$2e$ condensate. Those observations establish pairing and phase coherence more broadly; their existence alone does not specify a phonon mechanism or an isotropic weak-coupling gap. The converging isotope, spectroscopic, thermodynamic and coherence-factor tests support conventional BCS superconductivity, with explicit limits to its simplest universal ratios.

The first term assigns energy ${ϵ}_i$ to occupying the localized orbital at site $i$. The hopping term moves an electron between neighbouring sites, producing kinetic delocalization and band dispersion; Hermitian conjugate hops make the Hamiltonian Hermitian. The final term costs $U>0$ only when both spin states occupy the same local orbital. This is a Hubbard model with an on-site Coulomb penalty, not a classical interaction between magnetic dipoles.

At one electron per site and large $U$, a real hop is energetically costly, but a virtual hop into a doubly occupied intermediate state is possible. Opposite-spin states can exploit these virtual processes and lower their energy, while parallel spins cannot occupy the same one-orbital site. The resulting singlet–triplet splitting produces antiferromagnetic Hubbard kinetic exchange.

There is a terminology qualification. The paper calls this overlap-mediated process “direct exchange” because the hopping is directly between neighbouring magnetic orbitals. In the stricter distinction between direct Coulomb exchange versus kinetic exchange, the displayed Hamiltonian omits an explicit intersite Coulomb exchange integral. The antiferromagnetic $t^2/U$ mechanism derived below is kinetic exchange; a genuine intersite direct exchange term is separate and can compete with it. This distinction also separates the model from ligand-mediated superexchange in part f.

Use $|j⟩$ for the six printed basis vectors. The three singly occupied triplets areEach has energy zero: either Pauli exclusion prevents hopping to double occupancy or the two contributions cancel. DefineThe state $D_{-}$ has energy $U$. In the $S,D_{+}$ subspace the Hamiltonian is $\left(\begin{array}{cc}0& -2t\\ -2t& U\end{array}\right)$, withFor $U\gg |t|>0$, the fourth low-energy eigenstate is the normalized singletwhere $\mathcal{N}=[1+(E_{-}/(2t){)}^2{]}^{-1/2}$. Its energy isThe four lowest energies are $-4t^2/U,0,0,0$ to leading order, with the singlet and triplets just given. The upper pair lie near $U$ and $U+4t^2/U$.

The low-energy Hamiltonian can be writtenusing dimensionless spin-half operators. It places the singlet below the triplet by $J$, hence favors antiferromagnetic alignment. A lattice of such bonds can support antiferromagnetism, although its ordering pattern also depends on geometry and competing couplings. The physical energy gain comes from virtual charge fluctuations, not from a permanent double occupancy in the large-$U$ limit.

Use an oxygen-mediated superexchange path between two magnetic metal orbitals. A minimal one-active-orbital three-site charge-transfer superexchange Hamiltonian is clearest in a hole representation:Here $d_i$ describes a metal hole, $p$ a ligand hole, and $\Delta >0$ is the energy needed to create the virtual ligand hole. The low-energy configuration has one hole on each metal orbital and no ligand hole, corresponding to a filled bridging ligand orbital in the electron picture. $U_d$ and $U_p$ penalize double hole occupancy on the corresponding sites.

Virtual metal–ligand hopping transfers charge through oxygen without producing a real mobile carrier. Eliminating the ligand gives an effective metal–metal hopping $t_{dd}^{\mathrm{eff}}\sim t_{pd}^2/\Delta$. Combining it with the large metal repulsion givesA ligand double-occupancy channel adds a contribution of the same sign for the simple half-filled straight-bond model, often organized as $J\sim (4t_{pd}^4/{\Delta }^2)[1/U_d+2/(2\Delta +U_p)]$. The precise coefficients depend on the orbitals and hopping conventions; the essential mechanism is virtual charge transfer plus spin-selective access to intermediate states. The effective spin interaction favors antiparallel alignment while the material can remain insulating.

For MnO's rock salt structure, the strong approximately $180^{\circ }$ Mn–O–Mn path joins next-nearest Mn sites. It is not the nearest fcc Mn–Mn separation used in the direct-overlap estimate of part e. The oxygen-mediated path avoids relying on the tiny direct overlap of the two distant 3d tails. A realistic model includes several Mn 3d orbitals, Hund coupling and the high-spin Mn ion; the single-orbital Hamiltonian above isolates the antiferromagnetic exchange channel. Approximately $90^{\circ }$ paths and different orbital occupancies need not have the same sign.

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.

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.

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 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.