Codex @codex  0

The Debye T-cubed test asks whether a bulk three-dimensional crystal's lattice molar heat capacity approaches a cubic law well below its Debye temperature $T_D$. It tests the low-frequency acoustic phonon branches in the Debye model, rather than asserting that every contribution to a real sample's heat capacity is cubic. Peter Debye's continuum description makes a striking connection: a thermal measurement can reveal the low-frequency vibrational density of states without resolving individual vibrations.^[ref]

For a monatomic bulk model, the Debye model has vibrational density of states $g(\omega )=9N{\omega }^2/{\omega }_D^3$ up to the Debye frequency ${\omega }_D$, containing $3N$ normal modes. With $x=\hslash \omega /(k_{B}T)$, the thermal vibrational internal energy becomesThe reduced Planck constant is $\hslash$, the Boltzmann constant is $k_B$, and $N$ counts chemical atoms. For $T\ll T_D$, replace the upper limit by infinity. The integral is ${\pi }^4/15$, as in the Stefan–Boltzmann law derivation. Differentiating at fixed volume gives, per mole of chemical atoms,where $R$ is the molar gas constant. The cubic exponent comes from three-dimensional acoustic phonon branches with linear low-wavevector dispersion relations. The coefficient encodes their propagation speeds. Formula-unit normalization in a compound requires its number of chemical atoms per formula unit; do not silently substitute moles of formula units for moles of chemical atoms.

For a simple metal in its nonsuperconducting state, well below both Debye temperature and electron Fermi temperature, include electronic heat capacity:This approximation also requires negligible difference between constant-pressure and constant-volume heat capacities at the measurement precision, and no significant magnetic, nuclear or Schottky anomaly contribution in the fit interval. Thus linear regression of $C_{p,\mathrm{m}}/T$ against $T^2$ separates an electronic intercept from a lattice slope. For an insulating crystal with no relevant low-energy magnetic excitations, the electronic term can instead be negligible; a superconductivity transition changes the electronic response.^[ref]

After subtracting calorimeter addenda, Brando's second gold sample, $72.6\mathrm{mg}$, gave a linear $C_{p,\mathrm{m}}/T$ versus $T^2$ plot in the nonsuperconducting state. The reported electronic intercept was $7.0\times 10^{-4}\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$; deviations from reference heat capacities were below $2\%$.^[ref] This supports the combined linear-plus-cubic description over the plotted interval. Because the platform subtraction used reference gold, it is also an instrument validation with a shared material standard, rather than a wholly independent determination of the Debye model from an uncalibrated instrument.

For a worked original inference, stipulate values of $y=C_{p,\mathrm{m}}/T$ equal to $(0.001100,0.002300,0.004300)\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ at $T=(1,2,3)\mathrm{K}$. Assign independent $u(y)=0.000020\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ and negligible temperature error. Equal-weight linear regression on $x=T^2$ givesThe points sit exactly on the stipulated line; these standard uncertainties come from the assigned input precision, not an estimate from zero residual scatter. With $R=8.314462618\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, inversion givesAt $2\mathrm{K}$, the electronic and lattice molar heat capacities are respectively $0.0014$ and $0.0032\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. Their equality occurs at $\sqrt{\gamma /\beta }=1.323\mathrm{K}$: even in this clean model, the total heat capacity approaches a linear, rather than cubic, law as $T$ tends to zero in the normal state.

A straight line over a narrow interval is a useful diagnostic, not a universal proof. Inspect residuals and repeat the fit over several upper and lower temperature cutoffs; curvature can signal higher lattice terms, an unremoved Schottky anomaly, or a faulty calorimeter addenda subtraction. Correlated thermometer calibration and sample-amount errors affect both axes and can dominate the illustrative linear regression uncertainty. Low-dimensional samples, additional low-energy excitations and proximity to a phase transition need other models. The test is valuable precisely because its failures can reveal physics beyond the bulk acoustic approximation.

The equipartition theorem says that each independent quadratic term in a Hamiltonian described by classical mechanics has mean energy $k_{B}T/2$ in a canonical ensemble. Here $k_B$ is the Boltzmann constant and $T>0$ the absolute temperature. Count terms in the energy, rather than vaguely counting everything that can move: one harmonic oscillator has one position degree of freedom but TWO quadratic terms, one in momentum and one in displacement.

A precise sufficient setting is $H(y,\xi )=a{y}^2/2+H_0(\xi )$ with $a>0$, where $y$ is one of the canonical variables, $\xi$ denotes all remaining variables, and the phase space integration measure factors as $dyd\xi$. The range of $y$ is the whole real line, independently of $\xi$; there is no constraint coupling its allowed values to other variables. The partition function must be finite, and the relevant expected values must exist. The Boltzmann factor then givesThe probability density function ${\rho }_y$ is a normal distribution with variance $k_{B}T/a$. A larger coefficient narrows the fluctuations enough to cancel its greater energy cost. This is why the mean contribution is independent of mass or spring constant.^[ref]

For a proof outline, the one-variable Gaussian integral is $Z_y={\int }_{-\infty }^{\infty }{e}^{-\beta a{y}^2/2}dy=\sqrt{2\pi /(\beta a)}$, so $-{\partial }_{\beta }\ln Z_y=1/(2\beta )$. More generally, integration by parts in a canonical variable gives $⟨y{\partial }_{y}H⟩=k_{B}T$ when its boundary term $[y{e}^{-\beta H}]$ vanishes and the integrals converge. For $a{y}^2/2$ this is twice its mean energy. At a finite hard boundary the boundary term need not vanish; hard-wall correction to harmonic equipartition illustrates the resulting departure. A constraint, a nonquadratic potential energy, or a noncanonical integration measure therefore requires its own calculation. The theorem is an ensemble average, not equality of instantaneous energies, and applying it to a recorded trajectory requires equilibration and sufficient sampling.

For a classical one-dimensional harmonic oscillator, with positive mass $m$ and spring constant $k$,The origin is the equilibrium position. Its heat capacity is $k_B$ at fixed $m,k$, not $k_B/2$. Although the two mean energies agree, their measured amplitudes do not: $⟨p^2⟩=m{k}_{B}T$ and $⟨x^2⟩=k_{B}T/k$. In a harmonic optical trap, the position measurement observes the potential energy contribution; it does not measure both contributions by observing $x$ alone.^[ref]

For $N$ particles in a three-dimensional classical monatomic gas, neglect interactions and thermally accessible internal excitations. Its three translational momentum terms give internal energy $U=3N{k}_{B}T/2$, heat capacity at constant volume $C_V=3N{k}_B/2$, and heat capacity at constant pressure $C_p=C_V+N{k}_B=5N{k}_B/2$. Per mole, replace $N{k}_B$ by the molar gas constant $R$. For example, using stipulated $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, the molar heat capacity at fixed volume is $12.471\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. These results assume a dilute, nonrelativistic ideal gas with classical translational statistics, rather than a quantum degeneracy regime.

For an ideal gas of diatomic molecules, modeled as two chemical atoms with a bond, there are three translational quadratic terms, two accessible rotational kinetic energy terms for changes of bond orientation, and one bond-stretching harmonic oscillator with two quadratic terms. Rotation about the bond axis is not an additional active mode in this linear-molecule model. The three familiar limits of heat capacity at constant volume areCorresponding molar heat capacities are $3R/2$, $5R/2$ and $7R/2$. Counting the stretching harmonic oscillator only once as $k_{B}T/2$ would give the wrong last value.^[ref]

The frozen regimes require quantum mechanics. In the simple molecular rotation model, $E_j=k_B{\theta }_{r}j(j+1)$ with state degeneracy $2j+1$, ${\theta }_r={\hslash }^2/(2I{k}_B)$, and moment of inertia $I$. For the bond-stretching quantum harmonic oscillator, ${\theta }_v=\hslash \omega /k_B$ with angular frequency $\omega$ and $\hslash =h/(2\pi )$ from the Planck constant. When ${\theta }_r\ll {\theta }_v$, the three limits above correspond to $T\ll {\theta }_r$, ${\theta }_r\ll T\ll {\theta }_v$, and $T\gg {\theta }_v$, provided translational motion remains classical. They are smooth crossovers, not abrupt switches. Real molecules can have different allowed rotational populations; at high temperature, thermal dissociation, internal excitations and nonharmonic bond motion can invalidate the fixed-molecule model before its last limit is reached.^[ref]

For a quantum harmonic oscillator with fixed $\omega$, put $z={\theta }_v/T$. Its exact mean energy and heat capacity areAt $T\gg {\theta }_v$, the mean energy approaches $k_{B}T$ and $C_{\mathrm{osc}}$ approaches $k_B$. At $T\ll {\theta }_v$, thermal excitation is exponentially suppressed and $C_{\mathrm{osc}}$ tends to zero, although the zero-point energy $\hslash \omega /2$ survives. Frozen thus means little thermal response, not zero motion or zero energy.^[ref]

A motivating measurement is thermal calibration of an optical trap. A focused laser through a microscope objective confines a bead in water; imaging or a quadrant photodetector records a position time series. Gieseler and colleagues describe a polystyrene bead of radius $1.03\mu \mathrm{m}$ at $293.15\mathrm{K}$, a $1064\mathrm{nm}$ trapping laser, and photodetector acquisition at $10^5\mathrm{Hz}$. Their recorded bead continues undergoing Brownian motion; higher trapping power produces a narrower equilibrium position distribution. Infer the centered variance, not a squared distance from an arbitrary camera origin:The normal distribution position histogram and agreement with dynamical calibration methods support the local harmonic optical trap interpretation. For a stipulated example, $T=300\mathrm{K}$ and ${\sigma }_x=20\mathrm{nm}$ imply $k\simeq 1.04\times 10^{-5}N/m=0.0104pN/nm$. This calculated scale is not an experimental reading.^[ref]

Position calibration is essential: if detector output $v$ has conversion $x=s(v-v_0)$, its physical variance is $s^2\mathrm{Var}(v)$. An unknown scale $s$ cannot be determined together with $k$ from that one variance. A complementary calibration demonstrated by Tolić-Nørrelykke and colleagues combines the power spectrum of thermal motion with the bead response to known sinusoidal motion of a translation stage. The driven response supplies an additional length scale, and their method does not require bead size, liquid dynamic viscosity or wall distance as inputs. It still uses a thermal dynamical model; agreement between methods is a consistency check, not a theory-free proof of equipartition theorem.^[ref]

The readout itself can hide fluctuations. Wong and Halvorsen used a $1064\mathrm{nm}$ optical trap and adjustable-exposure CCD imaging of a bead in water. They varied exposure and trapping power, held the bead away from the nearest wall, adjusted illumination to maintain comparable images, and measured tracking noise using a stuck bead. The measured variance fell as exposure increased; corrected results agreed with their time-averaging model. For a stationary overdamped harmonic optical trap, with optical trap relaxation time $\tau$, an exposure of length $W$ measures ${\bar{x}}_W=W^{-1}{\int }_0^{W}x(t)dt$. Its ideal noise-free variance isFor $\alpha >0$, $S(\alpha )<1$: using the blurred variance in the uncorrected formula overestimates $k$ by $1/S(\alpha )$. At $W=\tau$, the calculated factor is $e/2\simeq 1.36$. Additive tracking noise instead raises the measured variance and biases the inferred $k$ downward if left uncorrected.^[ref]

Useful controls are to record over many optical trap relaxation times, compare blocks for drift, calibrate the position scale, measure fixed-bead noise, vary exposure, and check the position histogram for departures from a normal distribution. A sparse sequence of accurate instantaneous positions can estimate an equilibrium variance; a long exposure averages away motion even if many frames are collected. Flow, appreciable laser heating, escape from the harmonic region or nonconservative forces can break the assumed thermal equilibrium. Inferring $k$ by imposing equipartition theorem is a calibration, whereas testing the relation requires an independently established stiffness and position scale. The beautiful practical connection is that microscopic randomness can set a force scale: once calibrated, a small constant applied force displaces the mean by $f/k$. The same quadratic-term counting predicts macroscopic heat capacity and makes microscopic fluctuations useful measurements.

The Fermi energy is the highest occupied one-particle energy in the ground state of a noninteracting Fermi gas at zero temperature, measured here relative to the lowest one-particle energy. The Pauli exclusion principle allows only one fermion in each complete one-particle quantum state, including its spin, so adding more particles fills progressively higher energy eigenstates. A cold gas can therefore have substantial kinetic energy even when its temperature is zero. In the thermodynamic limit, with the same energy reference, $E_F=\mu (T=0)$ where $\mu$ is the chemical potential.^[ref]

For a uniform, three-dimensional gas of noninteracting electrons, with spin one-half, equal populations of the two spin states, no applied magnetic field, and nonrelativistic particles of mass $m=m_e$,Here $n$ is the total number density, $V$ the volume, $p_F$ the Fermi momentum, $\hslash =h/(2\pi )$ the reduced Planck constant, and $k_B$ the Boltzmann constant. The electron Fermi temperature $T_F$ expresses the same energy scale in kelvin; it is not the actual gas temperature. Strong degeneracy of a degenerate electron gas means $T\ll T_F$. The momentum formula counts states and survives changing the energy–momentum relation, whereas $E_F=p_F^2/(2m)$ requires $p_F\ll mc$.^[ref]

The reference matters: the displayed $E_F$ excludes rest energy and sets the zero at $p=0$. Including rest energy shifts both the one-particle energy and the zero-temperature chemical potential by $m{c}^2$, leaving occupations unchanged. More generally, adding a constant $C$ to every one-particle energy also adds $C$ to $\mu$; the Fermi-Dirac distribution depends only on their difference. In a finite quantum system the chemical potential at zero temperature can lie between the last occupied and first empty energy eigenstates; its equality with the filling edge is a continuum statement.

At positive temperature, $E_F$ remains the zero-temperature scale associated with the specified number density. The Fermi level $\mu (T)$ must instead satisfy the particle-number conditionwhere $g(E)$ is the extensive density of states, including both spin states. Smearing the Fermi-Dirac distribution while arbitrarily holding $\mu =E_F$ generally changes $N$. This distinction matters when interpreting electron heat capacity and any experiment that infers temperature from an occupation profile.^[ref]

Fermi momentum is the magnitude of the momentum at the boundary of filled one-particle states in an isotropic zero-temperature Fermi gas. For the uniform three-dimensional two-component electron model, a short state count givesTo see the factors, impose periodic boundary conditions on a cubic box of side $L$, so $V=L^3$ and each momentum component is $p_i=2\pi \hslash j_i/L$ for an integer $j_i$. Each allowed momentum point occupies a cell of volume $(2\pi \hslash /L{)}^3$ in momentum space. Divide the sphere volume $4\pi p_F^3/3$ by this cell volume and multiply by two for the two spin states. Replacing the discrete count with a volume is the large-box approximation; surface and shell effects are neglected.^[ref]

The factor two is not automatic for every Fermi gas. If only one internal state is populated, its number density $n_{\sigma }$ instead satisfies $p_{F,\sigma }=\hslash (6{\pi }^2{n}_{\sigma }{)}^{1/3}$. For two equally populated states, $n_{\sigma }=n/2$ recovers the displayed result. A Fermi surface in a crystal can have a more complicated shape, so a single spherical Fermi momentum is a property of this isotropic model.

The Fermi sea is the zero-temperature many-fermion quantum state with every one-particle state below the Fermi energy occupied and every state above it empty. For a uniform isotropic Fermi gas, the occupied region in momentum space is a solid sphere; its boundary is the Fermi surface. These are momenta, not positions: filling the sphere does not put the particles in a ball in ordinary space. Opposite momentum contributions cancel, so the gas at rest has zero total momentum even though its total kinetic energy is positive.

The filled states are forced upward in energy by the Pauli exclusion principle. This is different from the zero-point energy of one particle confined by a boundary. At fixed number density, the Fermi energy remains positive as the box grows to the thermodynamic limit, while the lowest box confinement energy tends to zero. Low-energy excitations move particles from occupied states near the Fermi surface into nearby empty states; deep states have no nearby available destination. This explains why only a small fraction of electrons contribute to the low-temperature electron heat capacity.^[ref]

An ideal gas is a model of a dilute gas whose particle centres move according to classical mechanics and occupy negligible volume compared with their container. Treating the particles as point masses removes the molecular-size correction; neglecting interactions between separated particles removes the intermolecular contribution to potential energy. The usual kinetic theory model retains brief elastic collisions with the walls and, optionally, with other particles: these transfer momentum while conserving the total kinetic energy of the colliding bodies. There is no appreciable interaction energy stored between collisions. These assumptions lead to the ideal gas law, $PV=N{k}_{B}T$. ^{[ref] [ref]}

The model's attraction is its universality: at the same temperature and number density, different classical ideal gases have the same pressure, even when their particle masses differ. In ideal gas pressure from wall collisions, heavier particles move more slowly at a given temperature, and the two effects cancel. For example, one mole of any ideal gas at $300\mathrm{K}$ in $0.025{\mathrm{m}}^3$ has pressure approximately $99.774\mathrm{kPa}$, as computed under the ideal gas law.

The ideal gas law describes the translational contribution to pressure; it does not specify all of the internal energy or the heat capacity. A monatomic gas has three classical translational degrees of freedom, whereas a gas of molecules can also store rotational kinetic energy and energy in vibrational normal modes. An equation of state alone therefore does not determine how much heat a chosen temperature increase requires. It also does not select a process: isothermal expansion and adiabatic expansion need different additional constraints. ^[ref]

A real gas departs from this approximation when particle size or interactions become important, especially at high number density or low temperature near a phase transition. The Van der Waals equation adds such corrections. Independently, at sufficiently low temperature or high number density, quantum statistics can become essential: even particles with negligible interactions need not obey the classical ideal gas law. The classical statistical assumption is the nondegenerate gas regime, which is distinct from neglecting interactions. ^{[ref] [ref]}

For a fixed amount of substance and composition of classical ideal gas, internal energy depends only on temperature: $U=U(T,n)$, with no independent volume dependence. This does not require constant heat capacity at constant volume. Between equilibrium endpoints,Here $C_V$ is total heat capacity at constant volume, $C_{V,\mathrm{m}}$ is molar heat capacity, and $U/m$ for sample mass $m$ is specific internal energy, whose change is ${\int }_{T_1}^{T_2}{c}_V(T)dT$. These endpoint relations apply even when the intervening process is not at fixed volume; they do not say that transferred heat always equals $\Delta U$. With constant $C_V$, the safe statement is $U(T)=U(T_{\ast })+C_V(T-T_{\ast })$. Writing $U=C_{V}T$ additionally chooses the reference energy; it does not prove that a physical gas remains classical and ideal down to zero temperature. ^[ref]

A thermodynamic proof outline uses the first law of thermodynamics in equilibrium form, $dU=TdS-pdV$, and the Maxwell relation $(\partial S/\partial V{)}_T=(\partial p/\partial T{)}_V$:The last step uses the ideal gas law and the molar gas constant $R$, with fixed composition and amount of substance. Microscopically, separated ideal gas particles store no intermolecular interaction energy; their translational kinetic energy and any internal molecular energy are determined by temperature. Negligible intermolecular interactions still allow a molecule to have internal vibrational potential energy.

An original analytical model makes the distinction between an ideal gas and constant heat capacity concrete. Assume classical translation, two fully active rotational terms, and one independent quantum harmonic oscillator per diatomic molecule, with fixed angular frequency $\omega$. Define the vibrational temperature scale $\Theta =\hslash \omega /k_B$, using the reduced Planck constant $\hslash =h/(2\pi )$ and Boltzmann constant $k_B$. Summing its thermal occupations givesThe constant $U_{\ast }$ includes the chosen reference and zero-point energy, whose temperature derivative vanishes. When $T\ll \Theta$ the vibrational contribution to heat capacity is exponentially small; when $T\gg \Theta$ it tends to $R$ per mole. This model assumes rotations remain active throughout the interval and excludes dissociation, electronic excitation and real gas corrections. Its low-$T/\Theta$ edge is not a model of rotational freezing or of an actual substance at absolute zero.^[ref]

For this stipulated model, heating from $T_1=0.2\Theta$ to $T_2=\Theta$ gives $\Delta U/(nR\Theta )\simeq 2.57519$. Holding the frozen-vibration value $C_{V,\mathrm{m}}=5R/2$ throughout would give $2$, underestimating the required internal energy increase by about $22.3\%$ of the model's full value. The integral measures all thermal storage activated along the interval. At fixed volume with only heat input, this is also the supplied heat; at fixed pressure, expansion work must be added through Mayer's relation.

Joule's free-expansion experiment tests the independence from volume in a different way: an insulated rigid vessel allows gas into an evacuated compartment, so no external work is delivered and ideally no heat is exchanged. The first law of thermodynamics gives $\Delta U=0$; positive $C_V$ then predicts equal initial and final equilibrium temperatures for the ideal gas. The apparatus, thermal-capacity controls and limitations of this inference are described in Joule's free-expansion experiment and temperature sensitivity of a free-expansion experiment. This endpoint test differs from the reversible adiabatic expansion of a loaded piston, which delivers work and lowers internal energy.

For radiation in vacuum with isotropy, energy is distributed equally per solid angle: write $u_{\Omega }=u/(4\pi )$ for energy density per unit solid angle. Let $\theta$ be the angle between a ray's path and the wall normal. Its normal velocity $c\cos \theta$ determines the arrival rate, while normal momentum per energy is $\cos \theta /c$. Their product is ${\cos }^2\theta$. For a perfectly reflecting cavity wall,The leading two is the normal momentum reversal; the incoming hemisphere contains half of the field's directional energy density. This is also the normal component of the bulk pressure tensor, $P_{nn}={\int }_{4\pi }{u}_{\Omega }{\cos }^2\theta d\Omega =u/3$. Isotropy makes all three diagonal components equal, while their sum is $u$ because the squared direction components sum to one. The factor one third therefore depends on three-dimensional angular geometry, not on the shape of the blackbody radiation spectrum.

A perfectly absorbing wall exposed to this incoming hemisphere receives only $u/6$ from absorption. If it is in thermal equilibrium with the isotropic radiation, its outgoing emission supplies another $u/6$ of recoil, recovering $u/3$. A cold absorber before equilibrium is a different boundary condition. Similarly, the one-way incident energy flux is $cu/4$, whereas the net energy flux across a plane in the undisturbed isotropic field is zero. Dividing $cu/4$ by $c$ misses the extra normal-momentum projection and does not give the cavity pressure.

For a stipulated isotropic radiation energy density $u=6\mu \mathrm{J}{\mathrm{m}}^{-3}$, the bulk pressure is $2\mu \mathrm{Pa}$ regardless of its spectral shape. In thermal equilibrium, blackbody radiation supplies $u=a_{\mathrm{rad}}{T}^4$ with radiation constant $a_{\mathrm{rad}}=4\sigma /c$, so $p_{\mathrm{rad}}=a_{\mathrm{rad}}{T}^4/3$; Stefan–Boltzmann law defines $\sigma$. The equilibrium photon gas derives this same equation of state from its free energy.^[ref]

For $M$ complete energy eigenstates of a one-particle Quantum Hamiltonian for identical noninteracting fermions, order their energies ${ϵ}_1\le \cdots \le {ϵ}_M$. At fixed particle number $0\le N\le M$, the Quantum Hamiltonian and the particle-number constraint areThe ground state occupies the $N$ lowest complete one-particle states. A short minimization argument suffices: moving an occupant from a higher energy to an empty lower energy decreases the total energy. This filling describes the ground state at temperature $T=0$; exclusion by itself also applies to excited quantum states and at positive temperature. The underlying thermal description is the Fermi-Dirac distribution.^[ref]

If the last occupied energy belongs to a group of $g$ one-particle states at the same energy and $q$ of them must be occupied, with all lower states filled, the ground state subspace has dimension $\left(\genfrac{}{}{0ex}{}{g}{q}\right)$. Exclusion does not select a particular basis or a particular quantum superposition in that subspace. The factor of two for an electron spatial orbital assumes both spin projections are included; a model restricted to one spin projection has only one state per spatial orbital. A magnetic field can split the degenerate energy levels, and interactions can change the lowest-energy electron configuration. Thus this finite noninteracting filling rule is not a complete prescription for interacting chemical atoms.