In a dilute thermal helium gas with rapid formation–decay balance of beryllium-8, write its small equilibrium abundance schematically as . A narrow isolated Hoyle state gives a second-step coefficient proportional to . Multiplication yields
with the stabilizing decay strength absorbed into the coefficient. Here is the Boltzmann constant. The two positive offsets count energy relative to their respective entrance thresholds, not relative to the carbon-12 ground state. A resonance enhances capture at the relevant energy; the Boltzmann factor still selects how many encounters reach it. This is a sequential narrow-level approximation with fixed nuclear widths, excluding plasma screening and changes of composition or mass density.
William Alfred Fowler lists the Hoyle state above the carbon-12 ground state, above beryllium-8 plus an alpha particle, and above three alpha particles.[ref] Subtraction gives for the diagram.
Normalize the calculated rate at and . Holding alpha particle number density and the decay strength fixed gives the original sensitivity calculation
At , a deliberately chosen shift multiplies the rate by ; gives . These offsets are hypothetical perturbations, not alternative measured nuclear levels. The reference curve has local temperature exponent of a thermonuclear reaction : raising temperature by gives a factor in this fixed-environment calculation. A real star changes its stellar structure and composition as it burns, so these factors do not directly predict final carbon abundance or a constraint on life.
Figure 1.
Triple-alpha energy thresholds and calculated resonance sensitivity
. Original energy bookkeeping using Fowler's rounded offsets (1983 lecture, p. 179). The normalized narrow-resonance rate calculation holds alpha density and decay strength fixed. The plus/minus 100 keV shifts are deliberately chosen sensitivity tests. No plotted curve is laboratory data. Bound carbon levels lie below the diagram's energy window; the downward arrow indicates stabilization rather than a measured lifetime or branching fraction. Details.
Richardson-Dushman law Created 2026-10-18 Updated 2026-10-20
The Richardson-Dushman law models the magnitude of emitted electric current density as
Here is absolute thermodynamic temperature, is the work function, is the Boltzmann constant, and is the Thermionic Richardson constant. The signed electric current density points opposite to the outgoing electron motion; here is positive. For uniform emission over area , the emitted electric current magnitude is .
The Richardson-Dushman law assumes an emitter close to thermal equilibrium, a surface with a well-defined and approximately constant work function over the fitted temperature range, and predominantly thermal escape over the barrier. Interpreting a measured collector electric current with this law additionally requires an emission-limited thermionic current, known emitting area, and negligible return emission from the collector. Strong electric fields can lower the surface barrier or enable quantum tunnelling; then this zero-field expression is insufficient.
For a deliberately chosen teaching model, take and , expressing the barrier in electronvolts. These are illustrative parameters, not measured properties of a named material. Using , at the Boltzmann factor is about , yet
A 20 percent increase in temperature raises the predicted electric current density by about 69 times. For a chosen emitting area , these correspond to approximately and if collection is complete. The small escaping fraction is replenished by the external electrical circuit.
The sensitivity is especially clear from the logarithmic derivative at fixed :
In the chosen model at , a rise changes by about 2.5 percent to first order. Accurate temperature measurement matters even when electric current measurement is easy.
The rotational partition function sums Boltzmann factors over rotational energy levels, including their degeneracies and the appropriate indistinguishability convention. For a linear rigid rotor at high temperature, , with moment of inertia and rotational symmetry number .
Second virial coefficient Created 2026-09-24 Updated 2026-10-18
The second virial coefficient is the leading interaction correction to the ideal gas law: at fixed temperature, , where and is number density. For identical classical particles in a homogeneous three-dimensional gas, with an isotropic pair potential and the centre-to-centre separation,
The function is the Mayer f-function. The factor counts the volume of a thin spherical shell; the factor counts each unordered pair once. The result has units of volume. Repulsive regions, , contribute positively, while attractive regions, , contribute negatively. [ref]
This is a striking bridge from a microscopic pair potential to a macroscopic pressure reading. At sufficiently small number density, raises pressure relative to an ideal gas at the same temperature and number density, while lowers it. A single value of , however, cannot determine the entire pair potential: many different functions have the same weighted integral. Measurements over several temperatures supply more information.
The integral must exist. For example, a bounded short-range attraction outside a hard core has no difficulty. If instead as , then and the radial tail behaves like : it converges only for . A bare unscreened Coulomb potential therefore falls outside this formula's ordinary short-range setting. This condition on the tail does not by itself establish the stability or convergence of the full virial expansion.
For an exactly integrable toy attraction, choose the following square-well pair potential, with diameter , attraction depth , and range multiplier :
Set and , the second virial coefficient of a hard-sphere gas. Integrating the core and attractive shell separately gives the exact toy result
The core contributes ; the shell contributes . Heating weakens the relative Boltzmann factor enhancement of attractive configurations, so approaches the positive hard-core value as . For the illustrative choice , the zero occurs at
At this temperature, : cancellation removes the linear correction, not every correction. At low , the negative shell contribution becomes large, requiring a correspondingly smaller number density for a useful truncated virial expansion. The square well is a deliberately chosen model, not a fit to the carbon dioxide experiment.
Consider an original counting example: a square two-dimensional XY model of side , a microscopic core radius , and a fixed phase stiffness . This deliberately neglects screening by other vortices. Away from a unit phase vortex, and . Integrating the Goldstone-mode effective free energy over an annulus gives
Boundary shape changes finite terms. There are order core positions, so the positional entropy is , where is the Boltzmann constant. With , subtracting gives the displayed free energy.
For and , the dimensionless free energy cost is about at and at . The same logarithm appears in the cost and the positional reward: increasing the system size strengthens whichever has the larger coefficient. The corresponding one-defect Boltzmann factors, about and , are relative statistical weights, not normalized probabilities. The latter signals failure of a dilute single-defect estimate.
A single net phase vortex requires compatible boundaries: on a periodic sample use neutral pairs, as required by vortex charge neutrality. This argument identifies a possible threshold, but predicting the universal stiffness jump requires the renormalized phase stiffness, after smaller pairs have responded. It does not justify setting the transition temperature equal to a bare lattice coupling divided by a universal constant.
Spin model 2026-09-28
A spin model assigns a discrete or continuous spin variable to each site of a lattice or vertex of a graph and gives interacting configurations a Boltzmann factor. Spin models describe collective order and phase transitions in magnetic and other many-body systems.
A bead in a liquid keeps undergoing Brownian motion even when an optical trap holds its average position fixed. For a stationary harmonic optical trap in thermal equilibrium at absolute temperature , the equipartition theorem assigns mean potential energy to this one quadratic coordinate. With the Boltzmann constant and at equilibrium,
At the same temperature, quadrupling the spring constant divides the position variance by four and the standard deviation by two. The Boltzmann factor gives the complete equilibrium probability density function, a Gaussian distribution:
The prefactor makes ; its second moment is . The probability density function has units of inverse length, and its area gives probability. This density is exact for the ideal quadratic model. In a finite optical trap it approximates the occupied well when thermal excursions remain within the harmonic region and escape is negligible during observation. Driven motion, appreciable nonconservative forces or laser heating need not give this equilibrium density.
Figure 1.
Thermal force calibration model: exact Gaussian position densities for relative trap stiffnesses 1, 2 and 4 at a common temperature
. Details.
Here , and . The plotted density per unit is . Each curve has total area one; their dimensionless standard deviations are , and . The curves are exact model densities, not noisy experimental data.
Brownian motion becomes a measuring tool: a calibrated position measurement supplies , and the known temperature supplies . Once is known, a slowly applied constant external force shifts the mean by , so in the same harmonic range. Fluctuations therefore calibrate the conversion from displacement to force. Svoboda and Block used this equipartition theorem method when calibrating their bead traps. [ref].
For the example stiffness and , and . Using a variance measured in uncalibrated detector units would not determine that physical stiffness. Measurement noise, finite exposure time, drift and insufficient sampling of the correlated Brownian motion can bias the estimate; record positions over many optical trap relaxation times and account for the detector response. The deterministic relaxation curve above describes relaxation without thermal forcing; an individual thermal trajectory continues to fluctuate.
Transfer matrix for the one-dimensional Ising model Created 2026-10-03 Updated 2026-10-05
For nearest-neighbour interactions, a two-by-two transfer matrix assigns one Boltzmann factor to every adjacent spin pair. With periodic boundary conditions, summing all spin configurations gives the partition function as the trace of a power of this matrix.
Equating the site chemical potential with that of an internally structureless classical vapour gives . The coefficient includes the bound vibrational partition factor, the translational thermal de Broglie wavelength, and the positive binding Boltzmann factor. It has inverse-pressure units. Half coverage occurs at ; an unspecified prevents a numerical pressure scale, but not the functional isotherm.
Villain model Created 2026-10-14 Updated 2026-10-25
A Villain model replaces the cosine Boltzmann factor of the XY model by a periodic Gaussian weight. The Poisson summation formula makes its bond weight equal to . The periodic Gaussian retains integer winding defects and permits an explicit separation into spin waves and vortices in two dimensions. Its replacement of the cosine model is an approximation; subsequent transformations within the Villain model can be exact.