Heavy fermion materials are correlated metals whose low-energy quasiparticles have exceptionally large effective masses. Partly localized -electron moments hybridize with broader conduction states. The Kondo effect and coherence below the Kondo-lattice coherence temperature produce a narrow low-energy band, so a small energy change corresponds to a large momentum-space state density.
For the electronic heat capacity of a heavy Fermi liquid, the resulting density of states gives a large low-temperature electronic specific-heat coefficient,
and often an enhanced magnetic susceptibility. Here includes both spins and must be expressed per sample, mole or volume consistently with . A coherent Fermi-liquid regime often has with an enhanced coefficient . Above coherence, scattering by local moments can instead give a Kondo-like resistivity upturn.
The competition between moment screening and intersite magnetic interactions can yield a heavy paramagnetic metal, magnetic order, quantum-critical behaviour or unconventional superconductivity. Ce- and U-based compounds supply familiar examples. “Heavy” describes the emergent electronic quasiparticles, not the ionic mass; the simple heavy-band picture need not survive at a quantum critical point.[ref]
Figure 1.
Hybridization flattens a conduction band near a localized level, while an itinerant exchange field splits spin bands
. Source. Details.
Assume three-dimensional parabolic band edges, with density-of-states masses and one valley per edge, each with spin degeneracy two. Their effective densities of states are
For the semiconductor mass-action law, integrating the square-root density of states against the Boltzmann tail gives and . Multiplication cancels the chemical potential, giving
If there are additional valley degeneracies, multiply by their product; anisotropic bands use the appropriate density-of-states masses. In the later equal-mass model, . No numerical mass or valley multiplicity has been invented.
When , the conduction-edge occupation is exactly , so the nondegenerate approximation for electrons fails at every positive temperature. Define the Normalized complete Fermi-Dirac integral
For the carrier densities with chemical potential pinned to a band edge, with the two specified equal masses, put . Direct Fermi-Dirac integration then gives
Thus exactly within the parabolic model. At low temperature, is thermally activated; at , from below. Both vanish at zero temperature and increase at fixed chemical potential. The other states that pin can supply the compensating charge; one must not additionally impose .
For the Pauli reduction of the pinned-edge mass-action product, the exact product relative to the classical law is
It approaches at low temperature and in the formal high-temperature limit. The low-temperature activation factor is unchanged but its prefactor is reduced; the classical mass-action formula overestimates the product. The high-temperature limit here refers to the given unbounded parabolic bands, not a claim that real semiconductor band structures remain unchanged at arbitrarily high temperature.[ref]
Finally, for a fixed positive , the zero-temperature Fermi sphere has radius , so
The low-temperature expansion at fixed gives the fixed-chemical-potential Sommerfeld density increase:
The monotonic density at fixed chemical potential in a square-root density of states initially rises quadratically from . For , its leading behaviour crosses to . The fixed chemical potential is what allows the particle number to change.
Figure 1.
Pinning the chemical potential to or above a band edge changes carrier populations and the mass-action product
. Source. Details.
In the three-dimensional Debye model, the lattice has normal modes for atoms. Replacing their dispersion by a linear acoustic dispersion up to a cutoff gives the Debye heat capacity
At low temperature the upper limit can be replaced by infinity. Since the integral then equals ,
At high temperature each harmonic mode contributes by the equipartition theorem, giving the Dulong-Petit law
For one mole replace by . The curve starts with zero slope, rises as , and levels off toward the classical limit. This is the lattice contribution; a metal's electronic heat capacity is a separate low-temperature term.
The Debye wavevector follows by counting three polarizations inside the cutoff sphere:
With the supplied speed of sound treated as the effective Debye speed, and
Here . Real crystals have different longitudinal and transverse sound speeds; the one-speed estimate uses the stated value, rather than inferring those unprovided branches.
Figure 1.
Debye lattice heat capacity and the squared magnitude of the Drude conductivity
. Details.
For nondegenerate electrons in a parabolic conduction band,
when the density-of-states effective mass and valley multiplicity are fixed. The intrinsic carrier extrapolation across a crossover must therefore specify how the activation energy changes with temperature.
The usual intended estimate treats the quoted as a constant intrinsic activation energy between the two temperatures, and identifies the crossover by . Then
For , and , the exponent is and the prefactor is . Hence
Equivalently, a fixed-gap intrinsic model with gives this result. A crossover is gradual: if denotes the net ionized donor density and , charge neutrality gives the actual electron population , rather than exactly .
There are two limits to interpreting the numerical estimate literally. First, the printed data specify only at , without specifying or the valence density of states. They do not mathematically determine : without the intended activation assumption the answer remains
Second, using the original PDF's and density implies , extraordinarily large for a conventional semiconductor. This confirms that the numerical data should not be treated as a realistic material determination. No unprinted decimal-point correction has been assumed.
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 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 by
A typical oscillatory contribution is proportional to . The oscillations are periodic in , not in . 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.
Figure 1.
Quantized orbital crossings and their periodicity in inverse magnetic field
. Details.
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 is
and a common first-harmonic disorder factor is . Thus the Landau spacing must be sufficiently large compared with thermal broadening and disorder broadening, approximately and , 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 . 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.
Relaxation-time model. The Drude model treats mobile charge carriers as independent particles accelerated by the Lorentz force between scattering events. Collisions randomize the drift momentum toward local equilibrium and are represented by a uniform, field-independent momentum relaxation time. For the present calculation the carrier density, isotropic mass and relaxation time are fixed, the fields are uniform and stationary, and quantum orbital quantization and intercarrier correlations are neglected. The steady drift velocity obeys
Here is the number of carriers of that species per unit volume, their signed charge, their inertial or isotropic band mass, the momentum relaxation time, and their single-species electrical conductivity at zero field.
One-carrier tensor. With and , solving for the electric field gives
Since and ,
Thus the resistivity tensor, defined by , is
Its off-diagonal signs follow from the stated axes and signed . For zero transverse electric current, is the Hall electrical resistivity; reversing charge reverses it. The single-species longitudinal electrical resistivity is independent of in this ideal isotropic model.
Adding electrons and electron holes. Let both species have density and charges , with . The stipulated equal zero-field electrical conductivity means equal positive charge carrier mobility , even if their masses and momentum relaxation times separately differ. Put and . Inverting each electrical resistivity matrix gives
Both species experience the same electric field, so their electric current densities add. Therefore it is their conductivity tensors, rather than their resistivity tensors, that must be summed. The Hall parts cancel and
This is compensated two-carrier magnetoresistance. Although the net Hall field vanishes, each carrier species is bent in its own direction; maintaining the total longitudinal electric current requires an increasing longitudinal field. Equal carrier numbers alone do not ensure zero Hall response if charge carrier mobilities differ; both equalities given here are used.
Large-field limit and quantum oscillations. Expanding the exact diagonal entry,
For ,
Here the last equality applies to either species with its own mass and transport time; the stipulated charge carrier mobilities give the same threshold. This is the classical strong-field condition : carriers turn substantially before their drift momentum relaxes.
Observable quantum oscillations additionally require sufficiently sharp Landau levels and a sufficiently cold Fermi-Dirac distribution, approximately
When , the disorder field scale is comparable to the onset scale for large classical magnetoresistance. The thermal condition can demand a much larger field, and when even the disorder criterion is stronger. The quantitative damping factors in Q2(c) show that these are crossover scales rather than sharp visibility boundaries. A large smooth magnetoresistance therefore does not itself demonstrate Landau level; it can exist while thermal smearing suppresses oscillations. Conversely, the asymptotic condition is stronger than merely reaching unity, so there is no universal ordering of the two experimental onset fields without temperatures and lifetime information.
Scaling function. The material's actual zero-field longitudinal electrical resistivity is . Dividing the exact answer by this value gives
The assumption that stays constant is essential to eliminating further temperature-dependent parameters. is a dimensionful scaling variable; is its dimensionless version. The corresponding fractional increase is .
Figure 1.
Original equal-mobility two-carrier curves collapse under zero-field resistivity scaling
. Details.
What the supplied data test. In the original PDF's left panel, the field dependence is strongly temperature dependent. In the right panel, plotting the fractional magnetoresistance against the squared rescaled field makes the displayed 125, 150, 175, 200 and 225 K curves collapse approximately onto one straight line, over the field intervals where they overlap. This supports the proposed one-parameter scaling and, on that range, the quadratic field law. The 100 K curve does not collapse: it bends away and has a different slope in the rescaled plot. The indicated scaling is therefore supported for the displayed 125–225 K data, not for the entire displayed temperature range. The graph alone does not locate a sharp transition between 100 and 125 K or establish behavior above 225 K.
This empirical collapse is Kohler's rule. It is consistent with a common temperature-dependent charge carrier mobility or relaxation scale and a fixed carrier population; it is not a unique proof of exact compensation or of one microscopic scattering mechanism. A change in relative scattering times, carrier population or Fermi-surface structure could spoil the collapse, and this figure alone does not select the cause of the 100 K departure.
There is a factor-of-two normalization convention in the printed text. The introduced for each species gives , whereas the final paragraph calls “.” If denotes the experimentally measured zero-field electrical resistivity, the consistent slope is
If one instead keeps the literal single-species definition , the same model has . The choice rescales and offsets the axes by temperature-independent constants, so it does not change which temperature curves collapse, but it must be resolved before inferring a carrier density from the slope.
Dielectric screening is the rearrangement of a medium's charges in response to an applied charge. In a metal, mobile electrons accumulate near a positive perturbation and reduce its long-range potential. Thomas-Fermi screening treats the electron gas locally as an equilibrium gas whose density responds to the local electrostatic potential, ignoring gradients of the density on the microscopic wavelength scale.
Take , background permittivity , and electrostatic potential . An electron's potential energy is , so the local kinetic Fermi energy becomes . For a three-dimensional spin-degenerate parabolic band at zero temperature,
where is the spin-summed density of states per volume per energy. The induced charge is . Combining it with Poisson equation gives
For a point external charge , spherical symmetry and decay at infinity select
Thus the bare potential is cut off exponentially. Fourier transformation makes the same result transparent: , so the static longitudinal dielectric factor is relative to the background. For a free electron band, and .
Figure 1.
Local Thomas-Fermi screening and optical/acoustic modes of an illustrative diatomic chain
. Details.
The local linear approximation requires a weak perturbation and slowly varying fields on the relevant microscopic scale. It omits nonlocal Friedel oscillations and atom-scale structure; a nonlinear local-density treatment is needed near a sufficiently strong charge. Finite temperature replaces by the appropriate compressibility derivative . The nondegenerate limit gives classical screening rather than the degenerate Thomas-Fermi length. A larger density of states means that a smaller potential can move more carriers, hence stronger screening.
A normal Fermi liquid has interacting fermions whose low-energy excitations can be followed continuously from those of a Fermi gas, provided no phase transition destroys that correspondence. A Landau quasiparticle is a particle together with the surrounding many-body response. It has the appropriate charge and spin, a renormalized dispersion and a long lifetime close to the Fermi surface. Interactions remain in both the quasiparticle parameters and the interactions between quasiparticles; they have not been discarded.
Near an isotropic Fermi surface,
The last relation defines the local isotropic Fermi-liquid effective mass. The Landau quasiparticle energy functional expresses how changing occupations shifts other quasiparticle energies:
with normalization of the sum absorbed into . Its angular components give Landau parameters describing collective density and spin response.
The Pauli exclusion principle severely restricts the available final states for low-energy collisions. In an ordinary three-dimensional Fermi liquid with regular interactions,
At zero temperature the decay width decreases faster than , making quasiparticles increasingly sharp near the Fermi surface. Their coherent spectral weight is the quasiparticle residue , and the true electron momentum distribution has a discontinuity of height at zero temperature. The rest of the spectral weight is incoherent. The quasiparticle occupation itself still follows a Fermi-Dirac distribution.
The low-temperature electronic heat capacity is linear, , where is the spin-summed quasiparticle density of states. In an isotropic spin-one-half model, the uniform spin susceptibility and compressibility take the forms
Thus enhanced effective mass and residual interactions can separately change observable responses. The familiar electrical resistivity requires a mechanism relaxing electrical current, such as allowed Umklapp processes; collisions in a perfectly clean Galilean-invariant continuum conserve total momentum and do not alone produce a finite dc resistivity.
Zero sound is a collisionless collective density oscillation sustained by the quasiparticle interaction; ordinary first sound occurs when collisions establish local equilibrium during a cycle. Metals and normal liquid helium-3 supply important examples. The description fails in phases with no sufficiently sharp low-energy quasiparticles, such as generic interacting one-dimensional systems, and must be modified after pairing or another ordering transition.
Band edges and optical transitions. A direct band gap has the conduction minimum and valence maximum at the same crystal momentum. A near-edge photon then connects the two almost vertically on an energy-versus-wavevector plot, because its momentum is small on the Brillouin zone scale. An indirect band gap has the extrema at different momenta, so an additional momentum transfer, commonly supplied by phonon absorption or emission, is required for a near-edge radiative transition. This reduces its probability relative to an allowed direct transition.
Near the zone centre, an isotropic direct-gap model has and downward valence dispersions . The heavy-hole and light-hole bands are the upper branches with large and small positive hole masses: the heavy branch is flatter. In common cubic semiconductors these arise from the manifold, while spin-orbit coupling places the split-off valence band lower by . Its thermally generated hole population is suppressed by approximately relative to upper-valence holes, apart from density-of-states factors. Neglect is justified when ; it is not a universal exclusion of this band at arbitrary temperature.
Figure 1.
Original schematic direct and indirect band edges, including heavy, light and split-off hole branches
. Details.
For an indirect material such as Si, the conduction minimum need not be at . In applying the carrier-density formulae, must mean the actual minimum of the included conduction valley. The supplied scalar effective masses and densities of states constitute a simplified band-edge model; they do not specify Si's full valley anisotropy or multiplicity.
Carrier populations. Let be the Fermi-Dirac distribution. In the nondegenerate conduction tail, . Set and integrate the given density of states:
The integral is . Defining the effective density of states gives
The factor two is already included in the supplied density of states; it must not be added a second time.
A hole is an unoccupied valence state. For sufficiently below , . Substitution of gives in exactly the same way
Their product eliminates the chemical potential:
This is the semiconductor mass-action law at thermal equilibrium. In an intrinsic semiconductor, charge neutrality additionally gives
In a doped semiconductor the product still applies in this nondegenerate equilibrium model, but the two populations need not be equal.
Intrinsic chemical potential. Taking logarithms of yields the intrinsic-semiconductor chemical potential
For unequal band degeneracies, use the first form with the full . For the supplied mass ratio , the shift from midgap at is approximately , toward the valence edge.
Three temperature regimes. The original plot runs from high temperature on the left to low temperature on the right. At high temperature, thermal excitation across the gap dominates and rises rapidly with temperature. In the middle, donors and acceptors are effectively ionized, the material is extrinsic, and makes the plateau. At low temperature, dopant freeze-out binds increasing numbers of electrons to donor levels, reducing the mobile-electron population as the sample cools. Its activation scale is a donor binding energy, not the full band gap. Compensation and incomplete ionization affect that regime.
Gap estimate from the original figure. On the dashed intrinsic branch, . Hence band-gap extraction from intrinsic carrier activation uses the slope of versus . With and expressed in eV/K, two points give
The very compressed intrinsic branch is only approximately readable: near and near . These give , rather than a precision determination. Neglecting the thermal prefactor gives about from those same rough points. Thus a gap of order – is justified by this schematic plot; its resolution does not support distinguishing a few hundredths of an electronvolt. The declining freeze-out branch must not be used to infer the full gap.
Room-temperature densities. At , the supplied masses give
In the fully ionized extrinsic plateau, put . Charge neutrality and the mass-action law imply
Since , is insensitive to the imprecise gap. The minority density is much more sensitive:
Using gives and . Using instead gives and . These are explicitly conditional estimates within the stated mass model. The plot alone cannot fix a precise minority density because the uncertain gap occurs inside an exponential.
Noninteracting spectral line. Treat as a wavevector, so physical crystal momentum is . Write ; an energy zero at the chemical potential can be chosen throughout. The temporal factor of a noninteracting state is . With the Fourier transform convention ,
Thus its normalized frequency spectral function is
One should not square the delta distribution as if it were an ordinary finite-time transition amplitude. More precisely, the fermionic retarded Green function is . Its Fourier transform is , whose negative imaginary part divided by is the same delta peak.
Interactions and lifetime. Interactions change the dispersion through the real part of the self-energy and permit decay into other excitations through its imaginary part. In the simplest quasiparticle-pole model, exponential decay of the retarded amplitude gives
Performing the causal integral,
The corresponding Lorentzian quasiparticle spectral peak is
Its total frequency integral is one. The exponential amplitude lifetime is ; the probability lifetime is . Describing as a scattering rate therefore requires the lifetime convention to be specified. The energy half-width is .
A general interacting electron is more accurately described by near the pole. The coherent peak then integrates to the quasiparticle residue and an incoherent continuum carries the remaining weight. In an energy-unit self-energy convention,
Linearizing the real denominator gives and at the pole. The requested unit-weight Lorentzian is the normalized quasiparticle model, not the most general exact interacting spectrum.
Peak position and width. At fixed , put and . Differentiation gives
It is positive for and negative for , so the only maximum is at , with height . Setting to half this height gives . Therefore
In energy units the maximum is and the half-width is . Frequency dependence of the self-energy would require reconsidering both the shape and peak position.
Fermi gas and Fermi liquid. A noninteracting Fermi gas has a delta-function spectral line at its bare dispersion for every momentum. A normal three-dimensional Fermi liquid has a shifted coherent peak plus incoherent spectral weight; close to the Fermi surface its width scales as . At zero temperature the peak becomes arbitrarily sharp at the Fermi surface, while it broadens farther away. Spectral functions include addition and removal weight; occupied-only photoemission contains an extra Fermi factor.
Figure 1.
A free spectral line and broadened quasiparticle peaks, with a separate incoherent background
. Details.
Photoemission measurement. A photon of known energy ejects an electron. Measure its kinetic energy and emission direction. For binding energy below the Fermi level and the appropriate calibrated work function , conservation of energy gives
Conservation of momentum parallel to a periodic surface gives
where is a surface reciprocal-lattice vector. At suitable photon energies the photon momentum can be neglected; at higher energies it must be retained. The surface breaks perpendicular translation symmetry, so the perpendicular crystal momentum generally needs a final-state model. A two-dimensional band can be reconstructed directly from the parallel components.
Within the sudden approximation and electric-dipole approximation, photoemission intensity and spectral function are related by
Here is measured relative to the chemical potential. Instrumental resolution convolves the product and inelastic scattering adds background. Thus ARPES maps the occupied electron removal spectral function, weighted by a matrix element; it does not measure the full unoccupied spectrum or an unweighted spectral function automatically. Bright ridges give dispersions, and intrinsic widths give decay information after instrumental effects are removed.
Reading the supplied image. Its bright ridges form a graphene-like cone. The crossing near is the Dirac point, not in general the Fermi level. In the paper's plotted energy coordinates its energy is roughly –. A velocity follows from a branch slope independently of the energy zero:
For example, the lower branch moves by about between and . This gives and
Reading other points changes the estimate by tens of percent. The original exam image has unequal energy increments on nearly equally spaced labelled ticks, so it does not justify a precise fitted velocity. The value is an estimate of the cone slope extrapolated to the Fermi crossing, rather than a measured derivative exactly at a visibly identified Fermi cutoff.
To determine a Fermi energy, first calibrate the Fermi edge, then subtract the chosen reference energy. Neither a Fermi-level marker nor an explicit energy-zero convention is supplied in this crop. If its upper termination near is intended to represent the Fermi cutoff, the conditional graph answer is on that axis, or relative to the Dirac crossing. Identifying the crossing itself as the Fermi energy would incorrectly assume undoped graphene.
The research article cited by the exam identifies this as linearly polarized graphene ARPES and independently states for the substrate-doped sample. That resolves the relative energy when this external calibration is admitted.[ref] The crop alone does not uniquely fix an absolute Fermi energy; the source-supported relative value and the conditional axis reading should be distinguished.
Work at zero temperature and fixed total carrier density. Write and . Let be the single-spin density of states per unit volume in the unpolarised state. Since and , the linear kinetic-energy shifts cancel and the quadratic cost is
The repulsive interaction instead decreases under spin polarization:
Hence the Stoner criterion for an opposite-spin contact repulsion is
The two forms differ only in whether the stated density of states counts one spin or both. For the three-dimensional free electron model, , so
Equivalently, for fractional polarisation , the kinetic-energy density is ; its quadratic coefficient is , while the interaction contributes . With number densities, has units of energy times volume and the displayed interaction is an energy density; multiplying by the sample volume does not change the threshold.
Ordinary doped semiconductor devices need a regime where their donor impurities or acceptor impurities supply mobile charge carriers while intrinsic thermal excitation remains sufficiently small.
At low temperature, donor impurities or acceptor impurities may no longer be ionized. Dopant freeze-out reduces the carrier density, raises electrical resistance and changes the intended carrier injection and p-n junction response. At high temperature the intrinsic carrier concentration grows as
when band effective masses and band gap are treated as constant. Once it becomes comparable with the net doping density, it overwhelms the intended electron-hole imbalance; reverse leakage and dissipation increase and useful p/n control is lost.
The usual operating window lies between freeze-out and intrinsic conduction, with mobility and thermal management placing further practical limits. This is a statement about conventional doped-device operation, not a fundamental prohibition on engineered cryogenic or wide-bandgap devices.
A current-carrying sample in a transverse magnetic field develops a transverse electric field because the Lorentz force deflects its charge carriers. In steady state this Hall voltage balances the transverse deflection. With Hall coefficient , one isotropic electron band gives and one hole band gives .
When electrons and electron holes both conduct, the weak-field two-carrier result is
The two-carrier Hall coefficient weights each species by mobility squared. Its sign changes when
Temperature changes dopant ionization, intrinsic electron-hole generation, and the two charge carrier mobilities, so it can move the system across this condition. For example, a p-type semiconductor can be hole-dominated in its extrinsic regime but become negative in its intrinsic regime if and . A sign change is possible, not mandatory: in many n-type semiconductors electrons dominate at both temperatures. The density majority alone does not determine the Hall sign in a two-carrier material.
Use a real, even interaction . Let and . In the Bogoliubov theory expansion, keep terms with zero or two noncondensed operators. The zero-momentum term gives . The terms with two noncondensed operators give
Here is the direct condensate interaction and is its exchange contribution. At fixed total particle number, , so
The second term cancels the direct normal term. Adding the kinetic energy yields the number-conserving Bogoliubov quadratic Hamiltonian
The constant is the condensate mean-field theory interaction energy; the normal term describes moving a particle out of the condensate, while the anomalous term creates or annihilates an opposite-momentum pair by exchanging two particles with the condensate. Total particle number remains fixed in the underlying theory; the condensate amplitudes in this notation are treated as classical amplitudes.
For condensate fluctuation power counting, put . Formal interaction-amplitude counting gives , cubic terms and quartic terms ; these are not uniform bounds on individual modes. Thus the latter are down by and relative to the quadratic interaction terms. The correction is relative . Homogeneity eliminates terms linear in nonzero-momentum operators. Small depletion and weak interactions, not low temperature alone, justify the truncation. [ref]
In the Bogoliubov quasiparticle vacuum, , but the inverse Bogoliubov transformation gives . Consequently the quantum depletion of a Bose condensate is
This expression assumes zero temperature. At nonzero temperature, thermal Bogoliubov quasiparticles add to the depletion:
Thus the printed formula is the quantum contribution, or a low-temperature approximation when the thermal contribution is negligible. It is nonnegative for a stable quadratic Hamiltonian and vanishes in the noninteracting zero-temperature limit.
A periodic harmonic chain with masses , spacing and nearest-neighbour spring constant has
Hamilton's equations give . The discrete Fourier transform , , with in one Brillouin zone, diagonalizes the quadratic form:
The harmonic-chain normal-mode reality condition means , not that and are independent real coordinates. For canonical quantization, impose , equivalently . For set
Then the quantum harmonic chain is a set of independent phonon oscillators:
The centre-of-mass zero mode of a harmonic chain is the free coordinate , with ; it cannot be normalized by dividing by . Retaining it gives
with the corresponding momentum expansion and uniform term . The nonzero modes alone give the site commutator ; the uniform term supplies the missing . Each oscillator occupation is a nonnegative integer, and adding one phonon costs and crystal momentum modulo a reciprocal-lattice vector. [ref]
In the continuum limit of a harmonic chain, , and
The wave equation and linear low- dispersion relation follow. The lattice gives an ultraviolet cutoff; the acoustic approximation does not reproduce the zone-edge flattening.
With the centre of mass fixed, the oscillator partition function and Bose-Einstein distribution are
For a long chain, above its finite-size level spacing and below the zone-edge scale, integrating the two acoustic propagation directions gives . At high temperature each retained oscillator instead contributes to the heat capacity.
Figure 1.
Original monatomic-chain mode sketch and its acoustic lattice dispersion
. Source. Details.
The Holstein–Primakoff transformation converts a spin algebra into bosonic operators while retaining the finite physical spin Hilbert space. With , write
The Holstein–Primakoff occupation constraint is . Acting on , the raising matrix element is and the lowering element is . Their squared difference is , reproducing . Operator ordering inside the square roots is essential. For low boson occupation relative to , expand them to obtain and . Keeping quadratic terms gives the linear spin-wave approximation; higher terms describe magnon interactions. Large helps control the expansion, but a low-dimensional infrared divergence can invalidate it even at large finite .
For the ferromagnet choose , , with unit lattice spacing and . Expand around the all-up state. A longitudinal bond contributes , and the transverse exchange contributes . 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 to ; 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 , with infrared behaviour . 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 , still , on an even bipartite chain. The classical reference state alternates up and down. Rotate one sublattice through so that both local quantization axes point along their reference spins. In unrotated physical components, to leading order , , and . Choosing Fourier phases using the actual site positions gives
The reduced Brillouin zone has momenta and length . The pair creation and annihilation terms express the fact that an antiparallel reference product state is not an eigenstate. A bosonic Bogoliubov transformation with , , diagonalizes each two-mode block. With the sign convention and , the result is
Thus 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 or excitations but contains original spin bosons, with . Its staggered magnetization is . This quantum depletion of Néel order diverges logarithmically at small 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]
Figure 1.
Original harmonic chain dispersions; the antiferromagnetic curve is a spin-wave approximation, not the exact spectrum for every spin
. Details.
For one coordinate with Hamiltonian , the propagator is , where . Divide into intervals , use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:
Multiplying the short-time kernels gives the time-sliced configuration-space path integral,
with , . Continuum notation abbreviates this limit as , where . The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.[ref]
When actions vary rapidly relative to , stationary phase selects classical paths. Put , with zero at the endpoints. Vanishing first variation gives . The quadratic term is
The resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,
The branch convention makes for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example gives the exact free-particle kernel . Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.
Quantum statistical mechanics follows by Wick rotation . The imaginary-time evolution operator is , with Euclidean action
Taking a thermal trace identifies and integrates the endpoint positions:
This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.
The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in spatial dimensions becomes a classical-looking field model in dimensions, with anisotropic couplings and extra extent . Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient , illustrating the mathematical correspondence.
The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in dimensions becomes
For indistinguishable particles there is a factor in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power is small. For a harmonic oscillator, the exact quantum result tends to as , recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.
For a uniform time-independent saddle, extremizing the action gives . Its grand-potential density is , with . Since , the minimum at positive density for is
For the uniform mean-field minimum is the vacuum, rather than a positive-density condensate. At fixed nonzero density in the weakly interacting condensed regime, the saddle relation sets . The mean-field one-body density matrix factorizes as , giving macroscopic occupation of the uniform mode. Depletion corrections are omitted at this stage.
The action has the global particle-number symmetry , . Choosing one saddle phase breaks this symmetry; the degenerate mean-field minima form the circle . There is no energetic preference for a uniform value of , but spatially varying phase costs gradient energy. Consequently a neutral condensate has a gapless Goldstone mode and a nonzero phase stiffness. The mode becomes sound once density dynamics is included. With the superfluid velocity is , so a single-valued field allows circulation in units around a closed loop that avoids zeros of the amplitude. This phase coherence and rigidity are the physical content of the broken-symmetry description.
The claim is a mean-field saddle statement, not an assertion that every and temperature in the question supports true long-range order. In a finite number-conserving system ; the symmetry-breaking order parameter is defined with the thermodynamic limit taken before an infinitesimal source is removed. In a homogeneous infinite short-range gas at nonzero temperature, infrared phase fluctuations prohibit true condensation in . In one dimension they also destroy true condensate long-range order at zero temperature. A two-dimensional gas can nevertheless have superfluidity and algebraically decaying correlations below a Berezinskii–Kosterlitz–Thouless transition. A finite trap can have substantial coherence without contradicting these thermodynamic statements.[ref]
At wavelengths with , neglect the density-gradient term in the quadratic density-phase action. Complete the square in its remaining density terms:
The convergent Gaussian integral can be performed by a contour shift or by its Fourier-transform identity. Its determinant is independent of at this order. Combining that normalization with the saddle action into a constant gives
The saddle contribution to is ; the additional measure-dependent constant need not be assigned an invented numerical value. The density Gaussian extends the fluctuation integration beyond its exact positivity domain, which is justified only while fluctuations around are small. The smooth zero-winding assumption is also required as discussed above. A constant global phase zero mode gives a finite compact-group volume, handled separately from nonzero phonon modes.
Fourier transforming the Bose-gas phase-only action gives a kernel . Analytically continuing gives its real-frequency poles,
In the present units the excitation energy equals . Each real Fourier coordinate is a harmonic oscillator with inertia coefficient and spring coefficient , hence frequency equal to the square root of their ratio. This is the direct relation to a quantum harmonic chain. The mode is a gapless sound phonon; its finite speed reflects both the compressibility and the phase stiffness .
Retaining the density-gradient term provides a useful check. Its integration instead gives a temporal coefficient , so
This is the Bogoliubov spectrum. Its small- limit is the derived sound branch; at large within the contact model it crosses towards particle-like behaviour. The crossover in the squared dispersion is , precisely the scale beyond which the discarded density gradient is not small.
The same action checks the dimensional restriction in the saddle analysis. At nonzero temperature its zero Matsubara mode gives phase fluctuations proportional to ; at zero temperature integrating the frequency gives . Their infrared divergences invalidate uniform true long-range order for at nonzero temperature and for at zero temperature. The long-wave collective description can remain useful for correlations even where a global condensate saddle is not a literal symmetry-broken state. The contact coupling here is an effective low-energy parameter; this calculation does not determine ultraviolet renormalization or trap-specific surface modes.
Figure 1.
Original dimensionless Bogoliubov dispersion and its sound and particle limits
. Details.