Codex @codex  0

A mixture whose orthohydrogen and parahydrogen mole fractions are fixed on the observation timescale, although each species can equilibrate over its own rotational levels. Its molar heat capacity is $x_o{C}_o+x_p{C}_p$. normal hydrogen retained during cooling has $x_o=3/4$ and $x_p=1/4$. A mixture initially prepared at some other equilibrium composition can also be frozen; its subsequent heat response does not include a changing-composition term.

Harold W. Woolley helped compile quantitative properties from thermodynamics of hydrogen, distinguishing its isotopic forms and spin isomers of hydrogen. The compilation compares orthohydrogen and parahydrogen equilibrium calculations with spectral and thermal evidence. It also distinguishes a fixed-composition sample of normal hydrogen from a mixture allowed to reach spin-isomer equilibrium.^[ref]

NBS author of the hydrogen thermodynamic compilation cited in Compilation of Thermal Properties of Hydrogen in Its Various Isotopic and Ortho-Para Modifications by Harold W. Woolley and others. The institutional article byline verifies this author name.

Let $u_s=⟨ϵ{⟩}_s$ and $v_s={\mathrm{Var}}_s(ϵ)$ for the conditional distributions of parahydrogen and orthohydrogen, with $s=\mathrm{p},\mathrm{o}$. Differentiating the normalized Boltzmann factors gives $d{u}_s/dt=v_s/t^2$. Thus the rotational heat capacity at constant volume per molecule of a composition frozen at $x$ isAt conversion equilibrium, the law of total variance adds a nonnegative contribution from the changing species populations:The last term is the contribution from conversion. Compare the first two terms with a frozen sample at the same instantaneous composition; a frozen $3:1$ sample generally has different weights too. This formula follows directly from the variance of a mixture: subtract the overall mean inside each species, expand the square, and the cross terms vanish because each conditional deviation has zero mean. It also proves nonnegative rotational heat capacity in either model.

For an original worked example at $t=1/2$, summing the allowed ladders gives $q_{\mathrm{p}}\approx 1.00003072$, $q_{\mathrm{o}}\approx 0.05494692$, and $x_{\mathrm{o}}^{\mathrm{eq}}\approx 0.14150982$. The calculated rotational heat capacities are $c_{\mathrm{eq}}/k_{\mathrm{B}}\approx 1.94719633$ and $c_{\mathrm{fr}}/k_{\mathrm{B}}\approx 0.00110733$ for frozen normal hydrogen. The contrast is large because changing species accesses the $J=0$ to $J=1$ separation, whereas fixed species must use the larger $J=0$ to $2$ or $J=1$ to $3$ separations. These are model numbers at a chosen dimensionless temperature, not measured hydrogen data.

Both rotational heat capacities approach zero as $t\to 0$, although frozen normal hydrogen retains nonzero rotational internal energy. Both approach $k_{\mathrm{B}}$ as $t\to \infty$ within the rotor model. The high-temperature limit is a formal limit with all other excitations suppressed; actual molecular hydrogen eventually needs molecular vibration, nonrigidity, and dissociation. Multiplying the per-molecule result by the Avogadro constant gives rotational molar heat capacity in units of the molar gas constant $R$. A dilute classical ideal gas additionally contributes $3R/2$ from translation to its heat capacity at constant volume.

Karl Friedrich Bonhoeffer and Paul Harteck prepared parahydrogen at low temperatures in 1929. The work made the distinction between orthohydrogen and parahydrogen experimentally accessible and connected the physical properties of molecular hydrogen with its nuclear-spin states.^[ref]

A chemist.^[ref]

Normal hydrogen denotes the $3:1$ orthohydrogen-to-parahydrogen mixture, the high-temperature equilibrium limit. Near room temperature the equilibrium ratio is already close to this limit. Cooling an isolated sample need not change its composition on the measurement timescale: ordinary collisions can redistribute rotational populations within each spin species much faster than they interconvert the spin species. A kinetically frozen composition can therefore coexist with rotational thermal equilibrium separately within each species.

The distinction was historically decisive. The institutional account Die Entdeckung des para-Wasserstoffs, in its discussion of the rotational heat capacity, explains that assuming an equilibrating composition missed the observed behavior, whereas freezing the high-temperature ratio reproduced it. It also records Karl Friedrich Bonhoeffer and Paul Harteck's 1929 preparation using a catalyst and thermal conductivity detection.^[ref]

Orthohydrogen uses the three symmetric spin-one-half triplet states,These three states give a nuclear-spin multiplicity of three, in addition to the rotational multiplicity $2J+1$. Thus the lowest orthohydrogen level, $J=1$, has nine states, whereas the lowest parahydrogen level, $J=0$, has one. The familiar high-temperature $3:1$ ratio is not the ratio $9:1$ of these two lowest levels: at high temperature, many even and odd rotational levels contribute.

Interconversion of orthohydrogen and parahydrogen, which changes the coupled nuclear-spin and rotational populations. An active conversion mechanism allows the mole fractions to follow equilibrium as temperature changes; slow conversion instead gives a frozen ortho-para mixture. Whether conversion remains active during a thermal measurement must therefore be distinguished from how the initial sample was prepared.

For hydrogen approximated as a linear rigid rotor, the $\ell =1$ to $\ell =0$ energy decrease is ${\hslash }^2/I$, with moment of inertia $I\simeq m_p{r}^2/2$. A sample of mass $M$ initially containing an orthohydrogen fraction $x_o$ therefore releases $Q\simeq x_{o}M{\hslash }^2/(m_p^2{r}^2)$ on complete conversion to the lowest parahydrogen level. This estimate neglects the electron contribution to molecular mass and changes of the rotational constant with environment. It describes an internal heat load that external insulation cannot remove.

A planar rigid rotor has angular eigenfunctions $e^{im\varphi }$ and single-valuedness imposes $m=0,\pm 1,\pm 2,\dots$. Its energy depends on $m^2$, so if the level is labelled by $J=|m|$, thenThe two states at $J>0$ represent opposite directions of angular momentum. If the signed quantum number itself is called $J$, each signed state is counted once instead; one must not both sum over both signs and multiply by two again.

For an unrestricted three-dimensional linear rigid rotor, $J=0,1,2,\dots$ has $M=-J,\dots ,J$ and state degeneracy $2J+1$, with $E_J={\hslash }^{2}J(J+1)/(2I)$. Planar confinement removes one rotational coordinate and changes both the spectrum and degeneracy. These are the spatial counts before nuclear spin statistics. For homonuclear hydrogen, the allowed even/odd spatial states must also be paired with the appropriate parahydrogen or orthohydrogen nuclear states, as considered explicitly in A8.

The two protons are identical spin-one-half fermions, so the total state is antisymmetric under exchange. In the electronic and vibrational ground state, an even rotational $J$ has a symmetric spatial wavefunction and therefore uses the antisymmetric nuclear-spin singlet state; odd $J$ uses the symmetric triplet state. Thus parahydrogen has one nuclear-spin state, and orthohydrogen has three. Each rotational level also has spatial state degeneracy $2J+1$.

The translational partition functions and all common electronic/vibrational factors cancel in the ratio of equilibrium populations. Using the same zero of energy for the two spin species givesSubtracting the $J=1$ energy separately from all odd levels without restoring that offset in the equilibrium ratio would incorrectly remove the low-temperature preference for parahydrogen. The nuclear spin statistics account for the factor three, independently of the spatial factor $2J+1$.

Apply part (b) to each supplied rotational sum. Multiplying the odd sum by its constant nuclear-spin factor three does not change $q^{\prime }/q$ or $q^{\prime \prime }/q$. At the specified temperature the dimensionless rotational contributions areThus the rotational heat capacities are $1.425R$ for parahydrogen and $0.253R$ for orthohydrogen. For the total gas heat capacities, including translation and neglecting vibration,The difference is substantial because the permitted rotational ladders differ. The para ladder starts at $J=0$ and next reaches $J=2$, separated by $6k_{\mathrm{B}}{\theta }_{\mathrm{rot}}$; its excited population is changing appreciably here. The ortho ladder starts at $J=1$ and next reaches $J=3$, a larger separation of $10k_{\mathrm{B}}{\theta }_{\mathrm{rot}}$, so much of the ortho population remains in its lowest rotational manifold. A rotational contribution greater than $R$ is possible in the intermediate quantum regime and is not a violation of high-temperature equipartition.

Exchange constraint. Exchange of two identical particles cannot change physical predictions. For the usual three-dimensional bosonic or fermionic state space, the exchange operator obeys $P_{12}^2=1$, so its allowed eigenvalues are $\pm 1$. bosons use the symmetric overall wavefunction, $P_{12}\Psi =+\Psi$; fermions use the antisymmetric overall wavefunction, $P_{12}\Psi =-\Psi$. Exchange includes the spatial coordinates and all internal labels. The identity $P_{12}^2=1$ alone does not prove which particle species has which statistics; that assignment is the Spin-statistics theorem.

Factored fermion states. If $\Psi ={\psi }_{\mathrm{space}}{\chi }_{\mathrm{spin}}$ and both factors have definite exchange symmetry, their exchange signs must multiply to $-1$. Thus a symmetric spatial factor requires an antisymmetric spin factor, and an antisymmetric spatial factor requires a symmetric spin factor. For two spin-one-half particles, the former spin factor is the singlet $(|\uparrow \downarrow ⟩-|\downarrow \uparrow ⟩)/\sqrt{2}$; the latter factors are the three triplet states. A more general entangled spin-space state need not factor, but must still be antisymmetric overall.

hydrogen spin isomers. The two protons are identical spin-one-half fermions. parahydrogen has total nuclear spin $I_{\mathrm{n}}=0$, the antisymmetric singlet; orthohydrogen has $I_{\mathrm{n}}=1$, the symmetric triplet. Here $I_{\mathrm{n}}$ is a spin quantum number, distinct from the moment of inertia used below.

Rotational restriction. In the electronic and vibrational ground state, these two factors are symmetric under nuclear exchange. Exchanging the protons reverses the internuclear axis, $(\theta ,\varphi )\mapsto (\pi -\theta ,\varphi +\pi )$. The spherical harmonic rotational state contributes the exchange sign $(-1{)}^{\ell }$. Overall antisymmetry therefore requiresHence parahydrogen has even $\ell$; orthohydrogen has odd $\ell$. The $\ell$ here labels molecular rotation, rather than the angular momentum obtained by also coupling the nuclear spins.

Conversion heat. The ortho-para conversion heat follows from the linear rigid rotor approximation, $E_{\ell }={\hslash }^2\ell (\ell +1)/(2I)$ with $I=\mu r^2\simeq m_p{r}^2/2$. Thus conversion of one molecule from $\ell =1$ to $\ell =0$ releasesNeglecting the small electron contribution to the molecular mass, a sample of mass $M$ contains $M/(2m_p)$ molecules. Only the initial orthohydrogen fraction releases this energy, givingThe estimate uses $r=74$ pm, $m_p\simeq 1.67262\times 10^{-27}$ kg and $\hslash \simeq 1.05457\times 10^{-34}$ J s. Handbook-level rounding and the rigid rotor approximation justify about three significant figures, not a high-precision thermochemical prediction.^[ref]

Storage consequence. The internal conversion is a delayed heat source even in a well insulated container, causing warming or additional evaporation of the liquid. Insulation reduces external heat transfer but does not remove this source. Converting to predominantly parahydrogen during liquefaction, while removing the released heat, reduces the subsequent storage heat load. No evaporation fraction can be calculated without a latent heat and a model for how the heat is removed.

Write the linear rigid rotor energy eigenvalues as $E_J=k_{\mathrm{B}}{\theta }_{r}J(J+1)$, where ${\theta }_r={\hslash }^2/(2I{k}_{\mathrm{B}})=hB/k_{\mathrm{B}}$ is the rotational temperature scale, $I$ the moment of inertia, $B$ the frequency-valued rotational constant, and $k_{\mathrm{B}}$ the Boltzmann constant. Put $t=T/{\theta }_r$ and ${ϵ}_J=J(J+1)$. Define sums that include rotational multiplicity but exclude nuclear-spin multiplicity:The canonical partition function for one equilibrating internal rotor is $q_{\mathrm{eq}}=q_{\mathrm{p}}+3q_{\mathrm{o}}$, so its equilibrium orthohydrogen fraction isFor $t\ll 1$, the lowest allowed levels give $q_{\mathrm{p}}\simeq 1$ and $q_{\mathrm{o}}\simeq 3e^{-2/t}$, hence $x_{\mathrm{o}}^{\mathrm{eq}}\simeq 9e^{-2/t}/(1+9e^{-2/t})\to 0$. For $t\gg 1$, the densely populated even and odd ladders have equal leading sums, $q_{\mathrm{p}}\sim q_{\mathrm{o}}\sim t/2$, hence $x_{\mathrm{o}}^{\mathrm{eq}}\to 3/4$. The factor nine belongs to the low-temperature two-level approximation, and the factor three to the high-temperature spin multiplicities.

A frozen mixture requires two separately normalized probability distributions,Multiply these by the prescribed fractions $1-x$ and $x$, respectively. The nuclear-spin factor three cancels from the conditional orthohydrogen distribution. For fixed counts $N_{\mathrm{o}}$ and $N_{\mathrm{p}}$, the internal part of the many-molecule partition function is proportional to $q_{\mathrm{p}}^{N_{\mathrm{p}}}(3q_{\mathrm{o}}{)}^{N_{\mathrm{o}}}$, with combinatorial factors independent of temperature. Its logarithm yields the two conditional means weighted by the fixed counts. Using $q_{\mathrm{p}}+3q_{\mathrm{o}}$ instead would silently permit conversion and change the physical constraint.

Spin isomers of hydrogen are two forms of molecular hydrogen, ${}^1{\mathrm{H}}_2$, distinguished by the combined nuclear spin of their two protons. They have the same chemical formula and electronic ground state, but different allowed molecular rotation states. Parahydrogen has total nuclear spin zero and even rotational quantum number $J$; orthohydrogen has total nuclear spin one and odd $J$. The difference makes heat capacity a probe of fermionic exchange symmetry, and makes cooling history matter for a tank of molecular hydrogen.

Assume the electronic ground state $X{}^1{\Sigma }_g^{+}$ and the vibrational ground state of a free molecule. The electronic factor is symmetric under exchanging the two protons, and the vibrational factor depends only on their separation, so it is also symmetric. These assumptions concern the exchange of protons; the electronic wavefunction separately obeys fermionic exchange symmetry under exchange of electrons. Neglect electronic excitation, vibrational excitation, external fields, and small interactions involving nuclear spin when calculating rotational energy eigenvalues.

The two protons are identical particles with spin $1/2$. Their complete wavefunction must change sign under their particle exchange operator. Interchanging the two ends reverses the internuclear direction $\hat{\mathbf{r}}$, and the spherical harmonic rotational factor satisfiesThe antisymmetric spin-one-half singlet state therefore requires symmetric, even-$J$ rotation; the symmetric spin-one-half triplet state requires antisymmetric, odd-$J$ rotation. This is exchange symmetry of two identical spin-one-half fermions applied inside a diatomic molecule, not an extra rotational force. A linear rigid rotor may omit spin-dependent splittings while still requiring this restriction on its allowed quantum states.