Codex @codex  0

Exchange-split bands in a ferromagnet generally give a different spin-resolved density of states for each spin branch, $D_{\uparrow }(E_F)$ and $D_{\downarrow }(E_F)$. If spin is approximately conserved across a multilayer, its two spin channels conduct in parallel. A spin valve combines a comparatively fixed reference-layer magnetization with a soft free layer that rotates in a small applied field. Its readout depends on the relative angle, not simply the absolute spin population. Biasing the layers near a perpendicular alignment gives a useful linear response instead of operating at the flat parallel/antiparallel extrema.

In giant magnetoresistance, a metallic spacer separates ferromagnetic layers. Spin-dependent bulk and interface scattering make their channel resistances unequal. The density of available final $d$ states influences scattering of more mobile states; however, a large total density of states does not by itself guarantee a large conductivity. The scattering-weighted spin-channel asymmetry also involves velocity and scattering matrix elements. In the elementary identical-layer two-current resistor model, let each magnetic layer have resistances $a,b$ in the two spin channels and neglect spacer resistance. Parallel moments give channels $2a$ and $2b$, while antiparallel moments give $a+b$ in each channel. ThereforeOne favorable spin channel traverses both aligned layers efficiently; reversal removes that common easy path. A conducting spacer adds resistance, and spin mixing weakens the separate-channel approximation. The channel asymmetry is material/interface dependent, so positive magnetoresistance is not a theorem for every pair of different ferromagnets.^[ref]

In tunnelling magnetoresistance, the spacer is an ultrathin insulator and current flows across a magnetic tunnel junction. In the simple spin-conserving density-of-states tunnelling magnetoresistance model, assume spin-independent matrix elements and low voltage. Available initial and final states giveWriting $P_i=(D_{i\uparrow }-D_{i\downarrow })/(D_{i\uparrow }+D_{i\downarrow })$ and $G_0=C{D}_1{D}_2/2$, these become $G_P=G_0(1+P_1{P}_2)$ and $G_{AP}=G_0(1-P_1{P}_2)$. ThusThe effective tunnelling polarization can depend on the barrier and orbital symmetry, not just a bulk density-of-states sketch. Crystalline MgO junctions exploit symmetry-selective transmission, so the elementary product model is not a quantitative universal formula for them.^[ref]

Both devices convert a free-layer rotation into an electrical change. A useful approximate GMR interpolation is $R(\theta )=R_P+\Delta R(1-\cos \theta )/2$; the elementary tunnelling model instead gives $G(\theta )=G_0(1+P_1{P}_2\cos \theta )$. The spin-valve small-field sensitivity is a chain rule,magnetic anisotropy, biasing and pinning set the rotation susceptibility. Conducting-spacer GMR usually offers low resistance and supports current-in-plane or current-perpendicular geometries. TMR uses perpendicular tunnelling, often offers a larger fractional signal and a larger resistance-area product, but demands uniform, stable barriers and has bias dependence and breakdown limits. Neither is automatically the better small-field sensor: heating, Johnson and low-frequency noise, domain switching, hysteresis and stability set the detectable field. Reverse sweep direction, vary current and temperature, and independently test reference-layer pinning to distinguish a reproducible magnetic signal from drift or heating.

Oceanic magnetic anomalies. Newly formed oceanic crust contains magnetic minerals. As it cools through their magnetic blocking temperatures, it acquires thermoremanent magnetization recording the ambient field polarity. Continuing seafloor spreading carries that record away from the ridge. Geomagnetic reversals therefore generate approximately paired marine magnetic stripes parallel to the ridge, with older crust farther from its axis. A magnetometer measures the superposed field above this magnetized material; subtracting a reference field reveals magnetic anomalies. The anomaly sign and shape depend on magnetization direction, latitude, source depth and geometry, so a positive measured anomaly is not universally identical to normal rock polarity.^[ref]

Matching a sequence of stripes to an independently dated reversal chronology gives plate ages and spreading histories. A stripe of width $w$ formed over interval $\Delta t$ implies half-spreading rate $v=w/\Delta t$ under steady symmetric spreading; full separation speed is $2v$. Unequal widths on opposite flanks can reflect asymmetric spreading or later deformation. Alteration, remagnetization and overlapping magnetic sources complicate inversion. The sequence correlation is stronger evidence than an isolated stripe or a field-strength peak.

Rayleigh-Bénard convection. Consider a horizontal fluid layer of thickness $d$, hotter below by $\Delta T>0$. The conductive state is motionless with a linear temperature gradient. A rising warm parcel is buoyant, but viscosity resists its motion and thermal diffusion erases its temperature anomaly. In the Boussinesq approximation, density variations are retained in buoyancy while the flow is otherwise incompressible. The competing effects are measured by the Rayleigh number and Prandtl number,For stress-free impermeable, fixed-temperature boundaries, linear stationary disturbances with horizontal wavenumber $q/d$ giveMinimizing gives $q_c=\pi /\sqrt{2}$ and the critical Rayleigh number for an isothermal layer, ${\mathrm{Ra}}_c=27{\pi }^4/4\simeq 657.5$. Two rigid no-slip, fixed-temperature boundaries instead give approximately $1708$; the critical value is not independent of boundary conditions.^[ref]

Above onset, buoyancy can sustain overturning rolls and heat transfer exceeds pure conduction, with Nusselt number $\mathrm{Nu}>1$. Nonlinear advection eventually limits the disturbance; larger driving can produce time-dependent or turbulent states. The ideal Boussinesq problem does not require hexagonal cells: planform selection depends on boundaries and non-Boussinesq effects. Convection driven by surface tension at a free surface is a distinct mechanism. An adverse temperature gradient supplies driving, but overcoming diffusion and viscosity requires a finite threshold.

Assume an Earth-like reversing planetary magnetic field and magnetizable crust that acquires thermoremanent magnetization on cooling. Alternating polarity intervals then produce marine magnetic stripes parallel to each spreading ridge, with matching age/polarity sequences on its two sides. AB has half-spreading rate 10 mm/year, whereas AC has half-spreading rate 24 mm/year. Thus the same polarity interval $\Delta t$ produces widths $10\Delta t$ and $24\Delta t$ mm when time is in years: the AC stripes are 2.4 times wider. Age increases away from each ridge. The transform produces no new crust; bands end or are displaced at it, rather than forming a third symmetric stripe sequence centred on BC. This is the remanent-magnetization/spreading mechanism developed from marine magnetic observations.^[ref]

A more precise local sketch can be made in the northern junction frame. Let $x,y$ be distances in kilometres along the fixed tangent axes and $t$ the crustal age in Myr; numerically, mm/year equals km/Myr. Relative to the moving junction, the velocities are ${\mathbf{w}}_A=(24,10)$, ${\mathbf{w}}_B=(24,-10)$ and ${\mathbf{w}}_C=(-24,10)$. In A the steady age field is $t_A=\min (x/24,y/10)$, since backwards advection first reaches either the AC or AB ridge. Equal-age bands therefore join as right-angle corners within A, with the corners migrating along the ray $(24,10)$. In the B wedge, $t_B=-y/10$ and the stripes meet the transform. On the C flank directly supplied from AC, $t_C=-x/24$; other parts of C carry older crust whose inherited pattern requires earlier boundary history.

The sketch displays magnetization polarity, not a computed total-field anomaly amplitude. Actual anomalies depend on the planetary field direction, source depth, magnetization strength and observation geometry. The paper supplies no field history, so a numerical reversal chronology or an absolute positive/negative anomaly sequence cannot be inferred. Without field reversals, ridge accretion alone does not imply alternating polarity stripes.

The de Haas-van Alphen effect is an oscillation of equilibrium magnetization, and hence magnetic susceptibility, as Landau levels pass through the Fermi energy when the field changes. For a fixed slice perpendicular to the field, Onsager quantization givesAt $E_F$, successive crossings change $n$ by one, soHere $F=1/\Delta (1/B)$ is the oscillation frequency in inverse field. In three dimensions, integrating the oscillatory response over the momentum parallel to the field causes most sections to cancel by phase averaging. The leading surviving contributions come from stationary, hence extremal, cross-sectional areas. Low temperature and sufficiently weak scattering are needed to resolve the levels. Measuring several field directions therefore probes the shape of the Fermi surface rather than only a single scalar electron density.

For the cylindrical $\gamma$ sheet, $A_k=\pi k_F^2$. The measured frequency givesThere are $4/3$ electrons per site in this band, with two spin states per wavevector. Thus the occupied disk takes fraction $2/3$ of the square Brillouin zone:This uses the paper's circular, electron-occupied disk approximation and equal band populations.

In the tight-binding approximation, hopping to the two neighbors at $\pm a\hat{x}$ adds phase factors $e^{i{k}_{x}a}$ and $e^{-i{k}_{x}a}$. Their sum is $2\cos (k_{x}a)$; absorb twice the hopping matrix element into $t$. With no $y$ or $z$ hopping this yields $E=E_0+t\cos (k_{x}a)$. For $t<0$, occupancy is centered on $k_x=0$. The filling fraction of an independent cosine band is $2/3$, requiring the occupied strip $|k_x|<2\pi /(3a)$, so sheet $A$ consists of the two planes $k_x=\pm 2\pi /(3a)$; sheet $B$ similarly has $k_y=\pm 2\pi /(3a)$. Each extends along $k_z$, with no dispersion in that direction. The $\gamma$ sheet is a cylinder along $k_z$, of radius $k_F$.

Numerically the planar offsets are $5.35\times 10^9{\mathrm{m}}^{-1}$, the zone boundary is $\pi /a=8.03\times 10^9{\mathrm{m}}^{-1}$, and the cylinder radius is $k_{F}a=2.894$ in dimensionless coordinates. If $t>0$ instead, the band minimum is at the zone edge; occupying its lowest two-thirds puts the planes at $k_x=\pm \pi /(3a)$ or $k_y=\pm \pi /(3a)$, with the central third unoccupied. The sign of the hopping must therefore accompany an absolute placement of these open sheets. The following sketch uses $t<0$ and the stipulated independent-band model, before any avoided crossings caused by additional interorbital hybridization.

Different measurements reveal different aspects of electronic band structure; a convincing determination combines their constraints.

Angle-resolved photoemission. In angle-resolved photoemission spectroscopy, photons of known energy eject occupied electrons. Measuring their kinetic energy and emission angle reconstructs binding energy and the component of crystal momentum parallel to the surface:Here $\Phi$ is the sample work function when kinetic energy is referenced to the nearby vacuum; an analyzer energy scale requires its own work-function calibration. Momentum is conserved modulo a surface reciprocal lattice vector, neglecting photon momentum. Tracking the peaks maps occupied dispersion relations and the Fermi surface. The measured intensity also contains the matrix element, the Fermi-Dirac distribution and the spectral function, so an absent bright feature does not establish an absent band. For a three-dimensional solid, the perpendicular momentum requires a final-state model. Inverse photoemission spectroscopy probes unoccupied states by injecting electrons and detecting emitted photons.

Optical spectroscopy. Photon absorption and ellipsometry measure interband transitions, optical selection rules and the complex dielectric response. Because optical photon momentum is small on the scale of the Brillouin zone, direct transitions approximately preserve crystal momentum. Their onset constrains a direct band gap; a threshold involving a phonon can reveal an indirect band gap. Peak positions constrain energy differences between bands, but the response also depends on the joint density of states and optical matrix elements, so it does not uniquely determine every $E_n(\mathbf{k})$.

Quantum oscillations and cyclotron resonance. In clean samples at low temperature, the de Haas-van Alphen effect measures oscillations in magnetization and the Shubnikov-de Haas effect measures oscillations in electrical resistivity as the inverse magnetic field is varied. Their period in inverse magnetic field is $1/F$, where Onsager quantization relates $F$ to an extremal Fermi-surface area perpendicular to the field:Rotating the magnetic field maps the geometry of extremal orbits. The thermal damping determines the cyclotron effective mass,with the appropriate magnitude for a hole orbit. Cyclotron resonance supplies the corresponding dynamical mass through ${\omega }_c=eB/m_c$. These are detailed probes near the Fermi level, rather than measurements of all occupied bands. Diffraction supplies the reciprocal lattice needed to label the measured momenta, but diffraction by itself is not a measurement of electronic dispersion.

In a weak applied field, magnetic susceptibility relates magnetization to magnetic field intensity, $M=\chi H$. Its sign, temperature dependence and nonlinear field dependence distinguish several mechanisms.

Diamagnetism is the orbital response opposing an applied magnetic field, with $\chi <0$. Closed-shell atoms and filled bands generally give a weak, nearly temperature-independent response. In a simple free electron metal, Landau diamagnetism contributes ${\chi }_L=-{\chi }_P/3$; the factor is specific to the isotropic parabolic-band model and is not universal in crystals.

Paramagnetism gives $\chi >0$. Localized magnetic moments align statistically with the field against thermal disorder. For independent moments of angular momentum $J$,which is the Curie law. Exchange interactions often lead approximately to the Curie–Weiss law, $\chi =C/(T-\Theta )$, sufficiently far above any ordering temperature. The sign of $\Theta$ suggests predominantly parallel or antiparallel exchange, though it is not a complete classification. Conduction electrons instead exhibit Pauli paramagnetism:for $g\simeq 2$ and a spin-summed density of states per unit volume. Only states close to the Fermi energy change occupation, making this response nearly temperature independent when $k_{B}T\ll E_F$.

Ferromagnetism has spontaneous parallel magnetic order below a Curie temperature. Exchange interaction, not simply classical dipole attraction, provides the dominant alignment mechanism. In the convention $H=-J_{\mathrm{ex}}{\sum }_{⟨ij⟩}{\mathbf{S}}_i\cdot {\mathbf{S}}_j$, positive $J_{\mathrm{ex}}$ favours alignment. Magnetic domains can make a ferromagnet's total magnetization vanish before an external field aligns its domains. Domain motion and anisotropy produce saturation, hysteresis, remanent magnetization and a reverse field measured by magnetic coercivity.

Antiferromagnetism has oppositely ordered sublattices below a Néel temperature; equal opposing moments give zero net magnetization while retaining a nonzero staggered order parameter. A susceptibility anomaly and magnetic neutron diffraction reveal this order even when bulk magnetization vanishes. Ferrimagnetism also has opposing sublattices but unequal contributions, giving a nonzero spontaneous net moment. Above their ordering temperatures both can become paramagnetic. Measurements with a magnetometer as a function of field and temperature, supplemented by magnetic diffraction, therefore distinguishes a compensated ordered material from an ordinary paramagnet.

A magnetic field quantizes a charged carrier's cyclotron motion into Landau levels. For an isotropic parabolic band, neglecting spin splitting,The degeneracy of each orbital level per unit perpendicular area is $eB/h$ per spin. Increasing the field changes both the level spacing and the number of states in each level. Successive levels cross the Fermi energy, producing an oscillatory density of states and hence quantum oscillations in thermodynamic and transport quantities.

For a general closed semiclassical orbit, Onsager quantization relates its wavevector-space area to field:Setting this area to an extremal Fermi-surface cross-section gives successive inverse-field crossings separated byA typical oscillatory contribution is proportional to $\cos (2\pi F/B+\phi )$. The oscillations are periodic in $1/B$, not in $B$. In three dimensions, extremal cross-sections dominate because contributions from other longitudinal momenta cancel by phase averaging. Rotating the field therefore maps the shape of the Fermi surface.

Oscillations in magnetization are the de Haas-van Alphen effect; oscillations in longitudinal electrical resistivity are the Shubnikov-de Haas effect. They are most readily observed in clean samples at low temperature. The thermal damping factor isand a common first-harmonic disorder factor is $R_D=\exp [-\pi /({\omega }_c{\tau }_q)]$. Thus the Landau spacing must be sufficiently large compared with thermal broadening and disorder broadening, approximately $k_{B}T\lesssim \hslash {\omega }_c$ and ${\omega }_c{\tau }_q\gtrsim 1$, with the displayed factors giving the quantitative suppression.

The temperature dependence measures the cyclotron effective mass, and the field-dependent damping probes the quantum lifetime ${\tau }_q$. That lifetime need not equal the momentum relaxation time controlling electrical conductivity: small-angle scattering can broaden Landau levels while relaxing electric current inefficiently. Spin splitting can generate further modulation or resolved spin structure. Oscillations are not the same as the ordinary smooth magnetoresistance or cyclotron resonance; they probe quantized orbital state counting near the Fermi energy.

Specify the spin counting first. Let $g_s=g_s(E_F)$ denote the unpolarized single-spin density of states per atom per energy, and put $\delta n=n_{\uparrow }-n_{\downarrow }$. At fixed total occupation $n$, the two Hubbard model mean-field energies can be writtenThe common shift $Un/2$ can be absorbed into the chemical potential. Lowering the up-spin band by $\Delta$ adds $g_s\Delta$ carriers, while raising the down-spin band removes the same number. ThusIf $n_a$ is the atom number density, the magnetization is $M=n_a{\mu }_B\delta n$. The Stoner-enhanced magnetic susceptibility in SI volume units is thereforeHere ${\chi }_P$ is the noninteracting Pauli paramagnetism result. If the paper's $g(E_F)$ counts both spins, $g=2g_s$ and the same result isA susceptibility per atom omits $n_a$. Stating the spin and volume conventions avoids an otherwise spurious factor of two or a dimensional mismatch.

The Stoner criterion is $U{g}_s(E_F)>1$, equivalently $Ug(E_F)>2$ with spin-summed $g$. To see the instability, transfer $\eta$ electrons per atom from one spin band to the other. The band-filling cost is ${\eta }^2/g_s$, whereas the interaction energy changes by $-U{\eta }^2$. Their sum $(1/g_s-U){\eta }^2$ becomes negative at the threshold. The unpolarized state then ceases to be a local minimum, and nonlinear band filling determines the spontaneous ferromagnetism. The negative continuation of the displayed linear susceptibility beyond the threshold is not a stable paramagnetic response.

Energy and momentum measurement. In angle-resolved photoemission spectroscopy, monochromatic photons illuminate a clean surface in vacuum, and an electron energy analyzer measures the outgoing kinetic energies and emission directions. Measuring both gives the occupied electronic band structure, rather than only a momentum-integrated distribution of removal energies.

Choose a common energy reference for initial and final states. Let $E_i$ be the initial electron energy, $E_F$ the Fermi energy, and $E_{\mathrm{vac}}=E_F+\varphi$ the vacuum level immediately outside the sample. The binding energy is $E_B=E_F-E_i\ge 0$. For a direct emission event without additional losses, conservation of energy givesHere the work function $\varphi$ is the energy from the Fermi level to the local vacuum. The initial and final energies are total single-electron energies on the same scale, while $E_{\mathrm{kin}}$ is measured above vacuum. The diagram illustrates each interval independently.

Neglecting ultraviolet photon momentum, photoemission momentum conversion giveswhere $\theta$ is measured from the surface normal. Parallel crystal momentum is conserved modulo a surface reciprocal lattice vector; the surface breaks perpendicular translation symmetry. A three-dimensional sample therefore generally needs a final-state model to infer perpendicular momentum. The stated two-dimensional sheet avoids that ambiguity. The kinetic energy used in this formula is the sample-vacuum value; if the analyzer and sample have different work functions, their contact potential must be included when converting the analyzer's reading, as in photoelectron binding energy reference.

Spectral information. Define $ϵ=E_i-E_F=-E_B$ and the electron spectral function by $A(\mathbf{k},ϵ)=-\mathrm{Im}{G}^R(\mathbf{k},ϵ)/\pi$. In the sudden approximation, the basic photoemission intensity and spectral function relation isbefore instrumental broadening and background are included. Here $f$ is the Fermi-Dirac distribution, and $\mathcal{M}$ is the transition matrix element. The electron removal spectral function contains the occupied part of the spectrum. In a noninteracting band it has a delta peak at the band energy; interactions shift and broaden peaks and can add incoherent weight. Thus peak trajectories reveal dispersion, while line shapes can constrain lifetimes after resolution effects are accounted for. Dividing by known matrix and occupation factors can estimate $A$ where those factors are appreciable; intensity alone is not identical to $A$, and a missing peak can reflect a matrix-element suppression. States well above $E_F$ are not accessible to ordinary low-temperature photoemission.

Two independent band-structure probes. First, the de Haas-van Alphen effect measures oscillatory magnetization as Landau levels cross the Fermi level when a magnetic field changes. The inverse-field frequency iswhere $S_{\mathrm{ext}}$ is an extremal Fermi surface cross-sectional area perpendicular to the field. Changing field orientation maps its shape. This reveals the momentum-space geometry of occupied bands; the related Shubnikov-de Haas effect measures oscillatory resistance.

Second, absorption or reflectance spectroscopy detects interband transitions. Photon energies at absorption edges and peaks constrain the band gap and interband energy differences. For a direct band gap, negligible photon momentum selects approximately vertical transitions in the band diagram, weighted by the joint density of states and transition matrix elements. An indirect transition needs a phonon to carry the extra crystal momentum. These measurements provide evidence beyond a free-electron parabola, although an optical spectrum alone does not reconstruct every band's dispersion.

Reading the graphite spectra. For each direction, follow a moving peak through successive angles and convert each peak's own kinetic energy to $k_{\parallel }$; angle is not itself proportional to momentum over this wide range. The prominent dispersing valence feature starts at roughly $7$–$9\mathrm{eV}$ binding energy near $\Gamma$ and moves towards lower binding energy as $k_{\parallel }$ increases. Along $\Gamma K$ it approaches the Fermi level at the corner, then turns back to larger binding energy on continuing beyond the corner. Along $\Gamma M$ it reaches a finite minimum binding energy of roughly $3$–$4\mathrm{eV}$ near the edge midpoint and turns back; it does not reach the Fermi level there. Additional deeper peaks describe other occupied valence branches, with visibility varying across the spectra.

A useful qualitative description is the upper occupied $\pi$ band of the nearest-neighbour model of graphene: $E_{\pi }=-t|f(\mathbf{k})|$. It rises from $-3t$ at $\Gamma$ to zero at $K$, but only to $-t$ at $M$. The schematic curve in the figure captures those two different behaviours. It is not a full multi-orbital assignment or a quantitative fit to graphite, whose deeper valence features are also visible in the supplied spectra.

Angle needed to reach the corner. The corner's parallel wavevector isUse the near-zero binding energy of the upper valence band at that point. With the sample work function interpreted as the stated $3.8\mathrm{eV}$ separation below vacuum,HenceThe energy and angle readings in the PDF are approximate. If a peak at the corner were instead assigned a small binding energy $E_B$, the estimate becomesFor $E_B=0$–$1\mathrm{eV}$ this gives about $57$–$60^{\circ }$, consistent with the near-Fermi feature in that angular range. Substituting the full photon energy for the kinetic energy would omit the work-function cost and underestimate the required angle.

As $B$ changes, each quantized orbit area grows linearly with field while the Fermi surface cross-section stays approximately fixed. Landau energies repeatedly cross the Fermi energy, making the density of states there oscillate. If the relevant extremal orbit has area $S_F$, its crossings satisfyThus the oscillations are periodic in inverse field rather than field. In three dimensions integration over longitudinal momentum makes extremal cross-sections dominate.

One observable is the de Haas-van Alphen effect: the field-dependent occupation changes the electronic free energy, so its field derivative, the magnetization, oscillates. The corresponding oscillations of longitudinal electrical resistance are the Shubnikov-de Haas effect. Both probe Fermi-surface geometry.

The large effective field primarily represents exchange interactions, which arise from Coulomb interactions, electronic hopping and the Pauli exclusion principle. It is not usually the much smaller classical magnetic dipole field of neighbouring ions. Direct exchange involves overlapping magnetic-ion orbitals. Superexchange instead uses virtual charge transfer through a ligand while real charge transport remains gapped, allowing substantial magnetic order in an electrical insulator.

Use dimensionless spin operators ${\mathbf{S}}_i$ with eigenvalue $S=1/2$, and count each bond once. One consistent Heisenberg model convention isThen $J_{ij}<0$ favours a Heisenberg ferromagnet, while $J_{ij}>0$ favours a Heisenberg antiferromagnet. Replacing neighbouring spins by their mean values gives the effective fieldFor $z$ equivalent ferromagnetic neighbours with $J_{ij}=-J_F<0$ and uniform order, $M=n{g}_J{\mu }_B⟨S_z⟩$, soThe last equality is the same molecular-field result derived above, not an exact transition temperature for every lattice.

Microscopically, a half-filled orbital pair with virtual hopping $t_{ij}$ and a large on-site repulsion $U$ permits an energy lowering of the singlet but not the parallel-spin triplet, yielding antiferromagnetic $J\simeq 4|t_{ij}{|}^2/U$ in this convention. Ligand-mediated hopping is the analogous mechanism in many ionic insulators. Different orbital occupancies and nearly orthogonal ligand paths can suppress this antiferromagnetic term and permit ferromagnetic exchange to dominate. Thus the sign depends on electronic structure and bonding geometry, not simply whether the sample conducts.

For antiferromagnetic coupling on a bipartite lattice, the field on each sublattice is proportional to the opposite sublattice magnetization, with opposite alignment preferred. Two coupled self-consistency equations then permit $M_A=-M_B\ne 0$ and vanishing total magnetization: the order parameter is staggered. The single positive-$\lambda$ uniform equation describes ferromagnetism and must be extended to these sublattices to describe an antiferromagnet.

The spin commutation relations give $[S_{\mathrm{tot}}^z,S_n^{\pm }]=\pm S_n^{\pm }$. Using the product rule for a commutator,and the conjugate term also vanishes. The longitudinal field commutes with every $S_n^z$, so $[H,S_{\mathrm{tot}}^z]=0$. The Hamiltonian is invariant under the global spin rotation $U(\theta )=e^{-i\theta S_{\mathrm{tot}}^z}$: its continuous symmetry group is $U(1)$. Simultaneous energy and total-magnetization eigenstates can be chosen, and evolution preserves each fixed-magnetization sector. Local magnetization can still flow between sites; generally $[H,S_n^z]\ne 0$. In the fermion representation the conserved number will be $N_f=N/2-S_{\mathrm{tot}}^z$.

Because $S_{\mathrm{tot}}^z=N/2-N_f$, the exact finite-size answer is $m_N=1/2-N_f^{\mathrm{gs}}/N$, with $N_f^{\mathrm{gs}}$ determined by the parity-consistent filling above. For the thermodynamic magnetization of the quantum XX chain, the occupied fraction is the length of the negative-energy interval divided by the Brillouin-zone length:ThereforeThe zero-field limit has zero magnetization and the high-field limit has all spins up. The finite-ring curve consists of magnetization steps approaching this continuous limit. Strictly above $h_c$ the fully polarized ground state is unique; exactly at $h_c$ a zero-energy spin-flip state can be degenerate with it. This finite-size qualification does not alter the saturation field or limiting magnetization density.

The Holstein–Primakoff transformation converts a spin algebra into bosonic operators while retaining the finite physical spin Hilbert space. With $[a_n,a_m^{†}]={\delta }_{nm}$, writeThe Holstein–Primakoff occupation constraint is $0\le n_n\le 2S$. Acting on $|n⟩$, the raising matrix element is $\sqrt{n(2S-n+1)}$ and the lowering element is $\sqrt{(n+1)(2S-n)}$. Their squared difference is $2(S-n)$, reproducing $[S^{+},S^{-}]=2S^z$. Operator ordering inside the square roots is essential. For low boson occupation relative to $S$, expand them to obtain $S^{+}\simeq \sqrt{2S}a$ and $S^{-}\simeq \sqrt{2S}{a}^{†}$. Keeping quadratic terms gives the linear spin-wave approximation; higher terms describe magnon interactions. Large $S$ helps control the expansion, but a low-dimensional infrared divergence can invalidate it even at large finite $S$.

For the ferromagnet choose $H_F=-J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, $J>0$, with unit lattice spacing and $\hslash =1$. Expand around the all-up state. A longitudinal bond contributes $-J{S}^2+JS(n_n+n_{n+1})$, and the transverse exchange contributes $-JS(a_n^{†}{a}_{n+1}+a_{n+1}^{†}{a}_n)$. After a Fourier transform,These are ferromagnetic magnons. The polarized vacuum is an exact ground state of the full ferromagnetic Hamiltonian: every bond is in its maximum total-spin sector, minimizing that bond energy. It therefore has no quantum zero-point reduction of its magnetization. The one-magnon dispersion is also exact, although interacting many-magnon states need not be free bosons. Rotating the polarization generates the maximum-total-spin ground-state multiplet. Selecting a definite direction in a thermodynamic symmetry-breaking description breaks $SU(2)$ to $U(1)$; the two transverse spin components form a canonical pair, producing a single quadratic type-B Goldstone boson. A small aligning field gaps it by the Zeeman energy.

At nonzero temperature the one-dimensional ferromagnet is different: its thermal magnon number contains $\int dk[e^{\beta JS{k}^2}-1{]}^{-1}$, with infrared behaviour $\int dkT/k^2$. The divergent magnetization reduction signals the absence of finite-temperature long-range order in an infinite short-range isotropic chain, consistent with the Mermin-Wagner theorem. It does not disallow the polarized zero-temperature ground states just constructed.^[ref]

For the antiferromagnet use $H_{AF}=+J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, still $J>0$, on an even bipartite chain. The classical reference state alternates up and down. Rotate one sublattice through $\pi$ so that both local quantization axes point along their reference spins. In unrotated physical components, to leading order $S_A^z=S-a^{†}a$, $S_B^z=-S+b^{†}b$, $S_A^{+}\simeq \sqrt{2S}a$ and $S_B^{+}\simeq \sqrt{2S}{b}^{†}$. Choosing Fourier phases using the actual site positions givesThe reduced Brillouin zone has $M=N/2$ momenta and length $\pi$. The pair creation and annihilation terms express the fact that an antiparallel reference product state is not an eigenstate. A bosonic Bogoliubov transformation with $u_k=\cosh {\theta }_k$, $v_k=\sinh {\theta }_k$, $\tanh 2{\theta }_k=\cos k$ diagonalizes each two-mode block. With the sign convention $a_k=u_k{\alpha }_k-v_k{\beta }_{-k}^{†}$ and $b_{-k}=u_k{\beta }_{-k}-v_k{\alpha }_k^{†}$, the result isThus the harmonic antiferromagnetic modes have linear, rather than quadratic, long-wave dispersion. In an unfolded Brillouin-zone description the same result is represented by zeros at both zero and the antiferromagnetic wavevector; these are not extra independent degrees of freedom beyond the two-sublattice accounting.

The Bogoliubov vacuum has no $\alpha$ or $\beta$ excitations but contains original spin bosons, with $v_k^2=\frac{1}{2}(|\sin k{|}^{-1}-1)$. Its staggered magnetization is $m_s=S-M^{-1}{\sum }_{k\in \mathrm{RBZ}}{v}_k^2$. This quantum depletion of Néel order diverges logarithmically at small $k$ in one dimension. Linear spin-wave theory therefore does not establish a stable ordered antiferromagnetic ground state for this chain. Its nominal gapless dispersion is a result of expanding around an invalid ordered reference, not a general proof of the exact excitation spectrum. For the nearest-neighbour isotropic chain, integer-spin systems have a Haldane gap, whereas half-integer-spin systems are gapless; the spin-one-half chain in particular has fractional spin excitations rather than a simple dilute gas of ordered-state magnons. This distinction requires physics beyond the harmonic expansion, including the topological term in the continuum theory.^[ref]

The generic instrument uses a magnetically responsive bead, a controllable external field and optical tracking of bead position or orientation. A field gradient supplies translational force; rotating an appropriate anisotropic field or bead moment can supply torque. For a permanent dipole,For an isotropic unsaturated induced moment with effective susceptibility contrast $\Delta \chi$ and bead volume $V$, the induced magnetic bead force instead follows fromThe factor one-half accounts for inducing the moment. A spatially uniform field alone supplies no translational force; an isotropic induced moment parallel to the field supplies no orienting torque either. Saturation or magnetic anisotropy requires the corresponding magnetization model.

Position tracking needs a length calibration, and force calibration needs measured field gradients, a known reference load, or calibrated thermal fluctuations. For a generic tensioned tether of extension $\ell$, small transverse displacement $x$ costs energy $F{x}^2/(2\ell )$, so $F=k_{B}T\ell /⟨x^2⟩$ under equilibrium, small-angle and single-mode assumptions. Attachment geometry, bead rotation, wall drag and finite acquisition bandwidth can affect this calibration. These equations describe generic magnetic tweezers mechanics, without assigning a biological response.

TODO: biological single-molecule applications are user-deferred; no biological systems or protocols have been researched or solved.

The expansion is a local phenomenological description in a small scalar order parameter near a transition. Its coefficients summarize microscopic physics and generally depend smoothly on temperature. At zero field, spin-reversal Z2 symmetry requires even powers of magnetization; the conjugate field supplies the term $-hm$. The additive $F_0$ does not affect minimization with respect to $m$.

For the ordinary continuous-transition quartic model, $b>0$ ensures stability, and the quadratic coefficient changes sign, $a=a_{0}t+O(t^2)$ with $a_0>0$. If $b$ vanishes or is negative, stabilizing higher powers are required and the phase structure can change. Treating one uniform magnetization neglects spatial gradients, fluctuations and domains; the resulting exponents are predictions of mean-field theory.

Differentiate the equation of state at fixed temperature:Above the transition, zero-field equilibrium has $m=0$. Below it the stable pure phases have $m^2=-a/b$. ThusThis is the mean-field magnetic-susceptibility amplitude ratio. Below the transition, take a one-sided response in a selected ordered phase: the globally minimizing magnetization switches discontinuously between opposite phases when $h$ crosses zero. Its two-sided derivative at precisely zero field is therefore not the finite pure-phase magnetic susceptibility used in the ratio.

Assume $b>0$ and temperature-independent $d$ for the local model. The zero-field stationary points obeyFor $d\ne 0$, the favoured ordered sign is opposite to $d$. The nonzero extrema first appear at $a=d^2/(4b)$; their existence is not yet a coexistence condition. To find the equilibrium transition, require both $F^{\prime }(m_{\ast })=0$ and $F(m_{\ast })=F(0)$. Substituting $a=-d{m}_{\ast }-b{m}_{\ast }^2$ into the free energy givesHenceAt coexistence the entire potential factorizes:It has equal minima at zero and $m_{\ast }$ separated by a barrier. Every fixed nonzero cubic coefficient produces a first-order transition in this quartic model, with a finite magnetization jump before the quadratic coefficient reaches zero.

For $d^2/(4b)>a>a_{\ast }$, the ordered minimum is metastable and zero is globally stable. For $a_{\ast }>a>0$, the ordered state is globally stable and zero is metastable. At $a=0$ the zero minimum loses local stability. Thus the spinodals are $a=d^2/(4b)$ and $a=0$, distinct from $a_{\ast }$. For $a<0$, the opposite-sign nonzero state can also be a local minimum, but the cubic term favours the state of sign $-d$.

If $a=a_0(T-T_c)/T_c$, the transition is shifted to $T_{\ast }=T_c[1+2d^2/(9b{a}_0)]$. With $b,d$ constant and $a^{\prime }(T)>0$, the entropy jump from the disordered to ordered phase is $-a^{\prime }(T_{\ast })m_{\ast }^2/2$, giving a latent heat of magnitude $T_{\ast }{a}^{\prime }(T_{\ast })m_{\ast }^2/2$. As $d\to 0$, the jump and shift vanish and the original continuous transition is recovered.

A genuine zero-field uniaxial ferromagnet with exact spin-reversal Z2 symmetry forbids this cubic term. The modified model therefore presumes that symmetry has been changed; it must not be presented as a symmetry-allowed correction for the original ferromagnet.

Near zero the Langevin function has expansionFor positive $x$, $\mathcal{L}(x)$ rises toward one and its ratio $\mathcal{L}(x)/x$ decreases from $1/3$ to zero. Therefore a line $\tau x$ has only the zero intersection for $\tau \ge 1/3$, and also a nonzero positive intersection for $0<\tau <1/3$. In a fixed signed component there is a corresponding negative solution; in the isotropic vector problem the nonzero solutions have all spin-space orientations.

The mean-field transition temperature isAt the critical temperature the line is tangent at the origin. Just below it, $S\simeq \sqrt{(5/3)(1-T/T_c^{\mathrm{MF}})}$, the usual mean-field square-root onset. The ordered minimum is stable below this temperature while $S=0$ becomes unstable.

This is a mean-field prediction, not a genuine finite-temperature ordering transition for the ideal two-dimensional short-range isotropic model. Long-wavelength spin fluctuations invalidate the predicted nonzero thermodynamic magnetization at every $T>0$. The relevant obstruction is the Mermin-Wagner theorem for continuous spin symmetry; anisotropy, long-range interactions or a finite sample change the assumptions. Its original paper establishes the absence of finite-temperature magnetic ordering for one- and two-dimensional finite-range isotropic Heisenberg models.^[ref]

For uniform unconstrained boundary conditions, the positive gradient coefficient favors a spatially constant minimizer. Put $u=|\mathbf{m}{|}^2\ge 0$. The remaining free-energy difference is $a(T-T_c)u+b{u}^2$. Its minimum isBelow $T_c$ the amplitude is fixed but its direction is arbitrary; the set of uniform minima is a two-sphere. The free-energy minimum is $f_0-a^2(T_c-T{)}^2/(4b)$ below $T_c$ and $f_0$ above. The order parameter vanishes continuously, so this is a continuous, second-order transition within Landau theory.

The free energy has internal $O(3)$ spin-rotation/reflection symmetry; a chosen nonzero vector leaves an $O(2)$ subgroup invariant. Restricting to proper rotations, this is $SO(3)\to SO(2)$ spontaneous symmetry breaking. There are two transverse orientation directions and a massive amplitude direction below the transition. Spatially varying directions cost gradient energy, and imposed boundaries or defects can prevent a uniform minimizer. In strictly two-dimensional equilibrium with short-range interactions, the fluctuation caveat of part (d) still applies: minimizing this polynomial functional at mean-field level does not overturn the absence of finite-temperature spontaneous magnetization.