Codex @codex  0

The partial pressure of a component of a gas mixture is the pressure it would exert alone at the same temperature and total volume. For an ideal mixture, $p_i=n_{i}RT/V$ and $p={\sum }_i{p}_i$. The chemical potential of an ideal gas uses $p_i$ for that component rather than the total pressure $p$.

The particle-hole continuum is the range of energy and momentum transfers that can excite an occupied fermion state into an unoccupied state, leaving a hole behind. For a free gas at zero temperature, choose $|\mathbf{k}|<k_F$ and $|\mathbf{k}+\mathbf{q}|>k_F$; the excitation energy is ${ϵ}_{\mathbf{k}+\mathbf{q}}-{ϵ}_{\mathbf{k}}$.

Varying the initial momentum produces a continuum at fixed transferred wavevector. For small $q$ in an isotropic gas, its positive frequencies extend approximately up to $v_{F}q$. A collective mode overlapping these excitations can decay by creating particle–hole pairs; a zero sound branch outside that region can avoid this particular damping channel. Interactions change the response and effective parameters, so the free-gas boundary is a useful starting model rather than a universal exact spectrum.

A nearest-neighbour tight-binding model on a bipartite lattice, with identical on-site energy ${ϵ}_0$, has paired single-particle energies ${ϵ}_0\pm x$. At half filling, $f({ϵ}_0+x)+f({ϵ}_0-x)=1$ for the Fermi-Dirac distribution with $\mu ={ϵ}_0$. Thus its fixed-density chemical potential equals ${ϵ}_0$ at all positive temperatures within the isolated-band model. An anisotropic rectangular nearest-neighbour band retains this property; additional same-sublattice hopping or unequal on-site energies can remove it.

At fixed temperature, the continuum excited population of a three-dimensional ideal Bose gas reaches a finite maximum as the chemical potential approaches the ground energy. Adding particles beyond that maximum creates macroscopic ground occupation. Thus $N^{\prime }=\min (N,N_c)$ and $N_0=\max (0,N-N_c)$ to leading thermodynamic-limit order. Finite systems retain nonmacroscopic ground occupation below the nominal threshold.

Passive transport is net movement down the transported substance's electrochemical gradient without coupling that movement to another source of free energy. Facilitated diffusion is passive transport through an ion channel or a passive carrier protein. Providing such a route can speed passage, but does not turn an uphill transfer into a downhill one. For an uncharged solute under the ideal conditions above, net passive transport runs from higher to lower concentration; for an ion, the membrane potential also matters. See Cooper, The Cell: transport of small molecules.

In the figure, stipulate eight neutral solute tracers in A and two in B, with equal fixed volumes, identical surroundings, and no imposed electric potential difference. The concentration ratio is $c_{\mathrm{A}}/c_{\mathrm{B}}=8/2=4$, so A-to-B transfer has $\Delta \tilde{\mu }=-RT\ln 4\simeq -1.38629RT$, or $-3.458\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$ at the assumed temperature $T=300\mathrm{K}$. If this solute can cross, thermodynamic equilibrium has equal expected concentrations: five tracers per side on average. Individual crossings occur in both directions; an observed ten-particle snapshot need not contain exactly five on each side. These stipulated counts illustrate a driving force, not a measured flow rate or a prediction of the time to equilibrate.

For a separate illustrative ion calculation, keep $c_{\mathrm{A}}/c_{\mathrm{B}}=4$, but assign $z=+1$ and $V_{\mathrm{A}}-V_{\mathrm{B}}=-60\mathrm{mV}$ at $300\mathrm{K}$. Using $R\simeq 8.31446\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ and $F\simeq 96485.33\mathrm{C}\mathrm{mo}{\mathrm{l}}^{-1}$, B-to-A transfer gives $RT\ln 4-0.060F\simeq 3.458-5.789=-2.331\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. Thus an appropriate open ion channel can permit passive transport toward the higher concentration: the electrical contribution outweighs the opposing chemical potential contribution. This stipulated membrane potential is not inferred from the tracer counts.

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.

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.

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

Insert $D=\mathrm{diag}(1,-1)$ at the first undecorated site. In the eigenbasis of the transfer matrix, $D$ exchanges $v_{+}$ and $v_{-}$, and both diagonal matrix elements vanish. Explicitly,Thus its large-$L$ limit is zero at every fixed $T>0$. The result is exact even before taking that limit. Equivalently, simultaneous reversal of all spins pairs configurations of equal energy and opposite ${\sigma }_1$. Positivity of the transfer matrix at fixed positive temperature prevents a degenerate dominant eigenvalue and spontaneous magnetization in this one-dimensional zero-field chain. Selecting one aligned state at strictly zero temperature is a different order of limits.

Insert the decorating-spin observable into the one bond that contains it. The decorated-spin occupation insertion uses the matrix:The partition function still contains $L$ decorating-spin sums; one of them has been replaced by this weighted sum. HenceTaking the dominant eigenvector givesThis is the probability that the decorating site is in either nonzero state. It tends to $2/3$ at high temperature, when its three states are equally likely, and to one at low temperature, when the nonzero aligned states gain bond energy. Unlike the magnetization, this even observable is not forced to vanish by spin reversal.

For low temperature, factor out the largest Boltzmann factor from the dominant eigenvalue:Thus the Helmholtz free energy per site has the expansionDifferentiate at fixed $J$, remembering that $K=J/(k_{B}T)$. The entropy per site isThe nonzero residual entropy comes from an exponentially large ground-state degeneracy, generated by geometrical frustration. A triangle with antiferromagnetic energy $J[st+(s+t)\tau ]$ has minimum energy $-J$: if $s=t$, the apex must have the opposite sign; if $s=-t$, either apex sign is permitted. The two all-equal configurations instead have energy $3J$. Thus exactly one of a minimum-energy triangle's bonds is frustrated. The bridge can simultaneously minimize its energy at $-J$ by requiring $t=-s^{\prime }$, where $s^{\prime }$ is the next cell's first base spin.

For fixed $s$ and $s^{\prime }$, the number of allowed apex choices is therefore two when $s=s^{\prime }$ and one when $s=-s^{\prime }$. The zero-temperature counting transfer matrix on successive first base spins isThe eigenvalues of $G$ are $3$ and $1$. Dividing $k_B\ln g_M$ by the $3M$ sites and taking the thermodynamic limit gives $k_B\ln 3/3$, independently confirming the derivative calculation. This does not assert a unique ordered ground state: the extensive choice of frustrated bonds is precisely what leaves entropy at zero temperature. At high temperature, all Ising spins become independent and $s\to k_B\ln 2$ per site.

The Boltzmann factor factorizes over bonds. Define the spin-chain transfer matrix byHere ${\delta }_{ab}$ is the Kronecker delta. Summing the product of these matrix elements over all configurations, and closing the final bond by periodic boundary conditions, gives the canonical partition functionHence, writing $I$ for the identity matrix and $P$ for the all-ones matrix,Every element of $P^2$ is a sum of $Q$ ones, while every diagonal element of $P$ is one. ThusMore explicitly, for $u=(1,\dots ,1{)}^T$, $Pu=Qu$. If ${\sum }_a{v}_a=0$, then $Pv=0$. These complementary eigenspaces have dimensions $1$ and $Q-1$. Equivalently, $P/Q$ is the orthogonal projector onto the linear span of $u$. The eigenvalues and their eigenvalue multiplicities are thereforeSince the two matrices share these eigenvectors, the exact finite-ring canonical partition function isBoth terms matter in the finite-ring low-temperature limit: discarding the second would lose the correct ground state degeneracy factor and introduce a spurious single-wall contribution.

Because $S[x_0]=0$ at $x=0$, the factor $Z$ cancels from the ratio. Write $z=\beta \hslash \omega$ and $q=\sqrt{m\omega /\hslash }x$. The exact answer isAt high temperature, $z\ll 1$. Expanding the hyperbolic tangent givesThe displayed leading expression is the ratio of classical Boltzmann factors for potential energy $V(x)=m{\omega }^2{x}^2/2$ and $V(0)=0$. Indeed the classical canonical ensemble givesThe momentum Gaussian integral is independent of $x$ and cancels in the ratio. This is the classical phase-space limit of the thermal kernel. Here $\beta =1/(k_{B}T)$, with $k_B$ the Boltzmann constant; the associated classical position variance is $1/(\beta m{\omega }^2)$. At fixed $x$, $R$ ultimately tends to one as temperature tends to infinity, while its width grows.

At low temperature, $z\gg 1$. Now $\tanh (z/2)=1-2e^{-z}+O(e^{-2z})$, so at fixed $q$Unlike the classical high-temperature expression, this limiting width remains nonzero at zero temperature, reflecting the ground state of the quantum harmonic oscillator.

Insert the complete orthonormal basis of energy eigenstates on either side of the Boltzmann factor. The spectral theorem givesAt low temperature, the gap $\hslash \omega$ suppresses excited energy eigenstates. The harmonic oscillator eigenfunction of the ground state isThereforeThe common zero-point energy factor cancels. Dividing the kernel by its trace would give the position probability density of the canonical ensemble; at zero temperature this tends to $|{\phi }_0(x){|}^2$, whose variance is $\hslash /(2m\omega )$. Thus zero temperature projects onto the ground state and retains its quantum spread, rather than concentrating at the minimum of the classical potential energy.

There is also a check on the first correction. With $q=\sqrt{m\omega /\hslash }x$, the first excited harmonic oscillator eigenfunction satisfies $|{\phi }_1(x){|}^2=2q^2|{\phi }_0(x){|}^2$ and ${\phi }_1(0)=0$. Hence the spectral theorem gives $R=e^{-q^2}[1+2q^2{e}^{-z}+O(e^{-2z})]$, in agreement with part (c).

Write $K=\beta J$, where $\beta =1/(k_{B}T)$ is the inverse temperature, $T$ is temperature and $k_B$ is the Boltzmann constant. In the canonical ensemble, each spin configuration has probability $e^{-\beta H}/Z_N$. The nearest-neighbor Boltzmann factors multiply, so the spin-chain transfer matrix, in the order $+1,-1$, isMultiplying its entries around the ring and summing over every Ising spin yields the partition function:The trace implements the periodic boundary conditions; replacing it by an endpoint sum would describe a different partition function. The normalized eigenvectors and corresponding eigenvalues areThe negative eigenvalue is essential for the antiferromagnetic Ising model. Introduce their ratioThe exact finite-$N$ partition function and Helmholtz free energy areFor finite $K>0$, $|t|<1$, so $Z_N>0$ for both even and odd $N$. The last term retains the finite-ring effect of geometric frustration.

To insert an Ising spin into this trace, use the diagonal matrix $\mathsf{S}=\mathrm{diag}(1,-1)$. In the eigenbasis of $\mathsf{T}$, it exchanges $|+⟩$ and $|-⟩$, and hence has vanishing diagonal entries. The one-point correlation function is thereforeThe same conclusion follows directly from spin inversion symmetry: reversing every Ising spin preserves the energy and the Boltzmann factor while reversing $S_1$.

For the two-point correlation function, the two insertions divide the ring into paths of $n$ and $N-n$ bonds. ThusIn particular, the exact finite-ring antiferromagnetic Ising correlation isBecause the one-point correlation functions vanish, this is also the connected correlation function. It gives $C_0=1$ and $C_n=C_{N-n}$, as required by $S_1^2=1$ and exchanging the two paths around the ring. At $K=0$, distinct Ising spins are independent and $C_n=0$ for $1\le n\le N-1$; $C_0=1$ is understood separately if a power $0^0$ occurs in the displayed expression.

Take the thermodynamic limit $N\to \infty$ at fixed separation $n$ and fixed finite $K=\beta J>0$. Since $|t|=\tanh K<1$, the terms $t^{N-n}$ and $t^N$ in the finite-ring antiferromagnetic Ising correlation vanish, through either even or odd $N$. ThereforeMatching the decaying exponential function gives the correlation length, in units of the site spacing:If the physical spacing is $a$, the physical correlation length is $a\xi$. The alternating sign records the preferred opposite orientation of neighboring Ising spins; the finite correlation length shows that this preference does not maintain long-range Néel order at positive temperature.

For large $K$, $\tanh K=1-2e^{-2K}+O(e^{-4K})$, so $\xi \sim e^{2K}/2$: the correlation length becomes large on cooling, but diverges only at zero temperature. At infinite temperature, $K\to 0^{+}$ and $\xi \to 0$; the separate equal-site two-point correlation function remains $C_0=1$.

Keep $N$ fixed and put $q=\tanh K$, so $q\to 1^{-}$ as temperature tends to zero. For even $N$, the finite-ring antiferromagnetic Ising correlation becomesBoth alternating ground states have $S_1{S}_{n+1}=(-1{)}^n$, so their equal-weight canonical ensemble has the same two-point correlation function despite its zero one-point correlation function.