Codex @codex  0

Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path $\lambda$ small compared with the gradient length, and relaxation times short compared with the observation time.

For a dilute hard-sphere gas,where $d$ is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.

A simple kinetic estimate makes the common structure explicit. For a quantity $a(y)$ carried per particle, a molecule moving at angle $\theta$ to $+y$ brings approximately $a(y-\lambda \cos \theta )$. Expanding to first order and averaging an isotropic velocity distribution givessince $⟨\cos \theta ⟩=0$ and $⟨{\cos }^2\theta ⟩=1/3$. This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.

Momentum transport and viscosity. For a shear flow with mean velocity $U_x(y)$, the transported directed momentum per molecule is $m{U}_x$. The net flux of $x$ momentum in the $+y$ direction isHere $\eta$ is dynamic viscosity in $\mathrm{Pa}\mathrm{s}$. The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is $\nu =\eta /(nm)$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.

Heat transport. Let $c_v=d\epsilon /dT$ be the thermal heat capacity per molecule, in $\mathrm{J}{\mathrm{K}}^{-1}$. Substituting $a=\epsilon (T)$ gives Fourier's law,Thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a monatomic gas, translational equipartition gives $c_v=3k_B/2$. Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.

Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,The diffusivity $D$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.

Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as $T^{1/2}$; at fixed pressure, $n=p/(k_{B}T)$ andKinematic viscosity and thermal diffusivity also grow roughly as $T^{3/2}/p$ at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change $c_v$. Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.

Heat diffusion equation. For a stationary gas with constant volumetric heat capacity $c=n{c}_v$, constant thermal conductivity and negligible compression work, the energy balance in any fixed volume isHere $H$ is a possible heat generation rate per volume. The divergence theorem and Fourier's law yieldThermal diffusivity $\alpha$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is of order $\overline{v}\lambda /3$. In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.

Representative solutions. With no source, a steady one-dimensional slab between temperatures $T_a,T_b$ at $x=0,L$ has $T=T_a+(T_b-T_a)x/L$ and uniform heat flux $-\kappa (T_b-T_a)/L$. Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.

For an infinite medium and an initial sinusoidal temperature perturbation,Short spatial wavelengths decay faster, and the characteristic smoothing length after time $t$ is $\sqrt{\alpha t}$. A finite slab with both ends held at $T_b$ has a sine expansion with factors $e^{-\alpha (n\pi /L{)}^{2}t}$; its slowest relaxation time is $L^2/({\pi }^2\alpha )$.

An instantaneous heat input $Q_A$ per unit transverse area at $x=0$ in an otherwise uniform infinite medium gives the Gaussian solutionIts integral $\int c(T-T_b)dx=Q_A$ conserves the deposited heat. The peak falls as $t^{-1/2}$ while its width grows as $\sqrt{t}$; diffusion redistributes energy instead of transporting a sharp front at a single velocity.

Finally, a semi-infinite body initially at $T_i$ and held at surface temperature $T_s$ from $t=0$ hasIt has $T(0,t)=T_s$ and approaches $T_i$ far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

Origin of diffusion. In a material near local thermal equilibrium, microscopic carriers transfer energy from hotter regions to colder regions. On distances long compared with their mean free path and times long compared with relaxation times, the local constitutive law is Fourier's law,Combining this temperature-gradient-driven flux with local conservation of energy produces a diffusion equation. This description assumes diffusive thermal conduction; it is not generally valid for ballistic transport, and flow would add advective transport.

Energy balance with heating. Choose any fixed volume $V$, let $c$ be its uniform volumetric heat capacity and $H$ its heat generation rate per unit volume. The rate of increase of internal energy equals heat generation minus outward heat flux:Since $V$ is arbitrary, for constant thermal conductivityThe ratio $\kappa /c$ is the thermal diffusivity in ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Here $c$ is in $\mathrm{J}{\mathrm{m}}^{-3}{\mathrm{K}}^{-1}$, $\kappa$ in $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$ and $H$ in $\mathrm{W}{\mathrm{m}}^{-3}$, so both right-hand terms have units $\mathrm{K}{\mathrm{s}}^{-1}$. If thermal conductivity varies with position or temperature, the divergence $\nabla \cdot (\kappa \nabla T)$ must be retained rather than replacing it by $\kappa {\nabla }^{2}T$.

Replacing the fixed surface temperature by Newton's law of cooling changes the boundary condition but not the radial derivative in steady state. At the surface,HenceHere $\gamma$ is a surface heat transfer coefficient in $\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$, unrelated to the surface tension symbol used in C10. An independent energy-balance check uses a length $\ell$: the electrical input $I^2\ell /(\sigma \pi a^2)$ must equal the surface loss $2\pi a\ell \gamma [T(a)-T_{\mathrm{air}}]$, which gives the same answer.

The complete temperature profile isThermal conductivity controls the internal gradient, but cancels from the surface-temperature rise because all generated power must eventually cross the prescribed surface. This assumes positive constant $\gamma$, uniform electrical conductivity and negligible end or additional radiative losses; temperature-dependent material properties could make the heating and boundary problem nonlinear.

Take positive heat flux upwards, along increasing $z$, and denote thermal conductivity by $k$. Fourier's law gives $q_s=-k{T}_z{|}_{z=a}$. Since $\cos (n\pi )=(-1{)}^n$, differentiation of the literal printed expression givesIts negative young-age flux is another sign of the unphysical initial profile. For the explicitly corrected hot initial plate, the surface heat flux of a cooling plate is insteadThe old-age limit is $k\Delta T/a$. For $t\ll {\tau }_0$, replacing the Gaussian sum by its leading integral gives $q_s\sim k\Delta T/\sqrt{\pi \kappa t}$, recovering the half-space cooling model. Neither $k$ nor $\kappa$ is numerically specified, so the age-dependent flux remains symbolic as requested.

Use the supplied surface temperature, $420^{\circ }\mathrm{C}\simeq 693\mathrm{K}$, as the premise. With otherwise comparable material and geothermal gradient, this hotter boundary increases thermally activated material creep at shallow depth and places a chosen mechanical cutoff temperature closer to the surface. That effect alone would favour a thinner elastic layer than in a cool terrestrial oceanic plate.

Water weakening of rocks makes dryness act in the opposite direction. Water can enhance material creep; sufficiently dry material can remain strong at temperatures at which its wet equivalent relaxes. Experiments on dry diabase provide a concrete basis for a strong dry crust and possible mechanical coupling between crust and mantle, rather than requiring a weak crustal layer solely because the surface is hot.^[ref]

For an illustrative steady conductive geotherm with thermal conductivity $k_T$ and heat flux $q$, a fixed strength cutoff givesTherefore the comparison would scale aswhere all temperatures are compared consistently and a dry rheology may have a higher effective cutoff $T_{\ast ,V}$. This is a conditional strength-layer estimate, not a universal definition of $T_e$.

Venus can retain appreciable elastic thickness despite its hot surface; temperature favours thinning, while dryness favours strengthening. The supplied facts alone determine neither a numerical $T_e$ nor a definite ordering relative to Earth's oceanic plates. The comparison requires heat flow, subsurface thermal history, material rheology, interior water content and loading timescale. Little atmospheric water is evidence relevant to the hypothesis of dry rocks, not a direct measurement of deep water activity.

I would test whether the ancient farside crust grew locally with an early thermal/compositional gradient or acquired a substantial later accreted layer. The research programme should estimate a globally consistent crustal mass and chronology before tuning either mechanism to a single attractive map.

First, map crustal geometry with independent density constraints. Jointly invert GRAIL gravity and laser-altimeter surface topography, allowing spatially varying density and porosity rather than converting every anomaly to thickness. A new broad-band seismic network should include farside stations, long-lived polar power systems and a communications relay. Natural impacts and moonquakes, supplemented by calibrated active sources near landing sites, would constrain layer velocities, attenuation, crustal interfaces and the crust-mantle boundary. Crossing ray paths and receiver-side conversions are essential to distinguish a thick crust from an incorrect density or velocity assumption. A coherent buried layer or age-correlated discontinuity across the farside would favor addition or stacking; a gradual hemispheric compositional gradient would better support asymmetric local growth. Neither pattern is guaranteed to survive later impacts.

Second, obtain samples that address the initial contrast. Collect geographically distributed ancient highland rocks from both hemispheres, using basin ejecta, central peaks and drilling to sample beneath the impact-gardened regolith. Return samples with mapped geological context. Mineral chemistry and bulk Ca–Al–Mg–Fe abundances would test a primary refractory gradient; crystallization ages would date crustal construction, while shock-reset ages would date overprinting. Isotopic and trace-element measurements could identify an accreted population, but similar isotopes are possible for material derived from the same early debris system. A young deposit draped over older crust would favor accretion; coeval flotation-related crust with a systematic chemical gradient would favor local asymmetric differentiation. Multiple sites are necessary because a single basin's ejecta can mix distant sources.

Third, measure the thermal and volcanic consequences. Place shielded, sufficiently deep heat flow probes in representative highlands and mare terrains. Combine their geothermal gradients and measured thermal conductivity with orbital gamma-ray constraints on heat-producing elements. Map mare volumes and eruption ages, rather than only their exposed areas, to separate melt production from storage and eruption efficiency. A strong correlation between late volcanism, subsurface heat producers and high inferred mantle temperatures would support radiogenic amplification, while old crustal asymmetry without that correlation would show that later heating is not the entire cause.

Finally, compare forward models against withheld observations. Couple irradiation, vapor transport, condensation, lunar magma ocean mixing and crystallization for the thermal scenario; couple orbital evolution, impact deposition, mixing and later viscous relaxation for accretion. Account for later true polar wander before treating present-day hemispheric orientation as the original formation geometry. Require each model to conserve mass and energy and to reproduce the independent thickness, chemistry, age and volcanic data within their uncertainties. Vary unknown initial conditions and compare predicted observables at sites not used for calibration. The discriminating outcome is a conjunction such as a coherent accreted unit with a consistent emplacement age, or an ancient chemical gradient quantitatively consistent with transport under Earth irradiation. A lower gravity misfit alone cannot establish the historical mechanism, because gravity is a nonunique inverse problem.

Specify what “onset” means. I would distinguish the first local, perhaps short-lived subduction event from persistent, self-sustaining recycling in a connected global network of tectonic plates. Modern plate tectonics combines relatively rigid moving lithosphere, localized boundary deformation, seafloor spreading and sustained convergence. A sinking drip or a single imposed descending sheet is not sufficient evidence of that global regime. Different processes may become established at different times, so a transition interval may be a better hypothesis than an instantaneous worldwide switch.

Construct competing physical models. Solve conservation of mass, momentum and energy for evolving mantle convection and lithosphere, including temperature- and stress-dependent rheology, brittle yielding, plastic deformation, melting, compositional buoyancy and relevant phase changes. Use geodynamic modelling to compare a stagnant lid, episodic overturn, drip-dominated recycling and a sustained mobile lid. In a simplified incompressible creeping-flow framework,Here $\mathbf{u}$ is velocity, $p$ pressure, $\eta$ effective viscosity, $\dot{\mathbf{ϵ}}$ the strain-rate tensor, $\rho$ density, $\mathbf{g}$ gravitational acceleration, $c_p$ specific heat capacity, $k$ thermal conductivity, $H$ radiogenic power per unit mass, and $\Phi$ dissipative heating per unit volume. These equations are a framework: compressibility, latent heat and chemical transport may require additional terms. A simple yield stress law ${\tau }_y=C+\mu (\rho gz-P_f)$ represents cohesion $C$, friction coefficient $\mu$ and pore pressure $P_f$, illustrating why hydration and interface weakening can change the outcome.

Vary the early-planet controls rather than assuming the answer. Explore plausible initial mantle potential temperature, radiogenic heating, crust thickness and composition, water inventory, grain-size evolution, surface temperature, and boundary strength. Hotter mantle lowers viscosity but can produce thicker buoyant basaltic crust; it does not automatically make sustained subduction easier. A cold strong plate can transmit slab pull yet resist bending and initiation. Water weakening of rocks, melt-induced weakening, inherited faults and mechanically weak sediments can change this balance. Test proposed triggers such as mantle plume impingement or gravitational instability without presupposing an established global plate-force system.

Useful diagnostics include the Rayleigh number,with thermal expansion coefficient $\alpha$, layer thickness $D$, thermal diffusivity $\kappa$ and temperature contrast $\Delta T$, together with surface mobility, the localization of deformation, plate longevity, melt production and material recycling. A large Rayleigh number establishes vigorous convection under suitable conditions, not necessarily a mobile lid. A convincing model should sustain one-sided subduction, connected boundary evolution and seafloor spreading after any artificial trigger is removed. Compare initiation thresholds with continuation thresholds: hysteresis can allow an established regime to persist under conditions in which it would not spontaneously start.

As an example of the type of experiment, three-dimensional models have tested mantle plume-induced subduction initiation, finding sensitivity to plate buoyancy, magmatic weakening and hydration of the interface. The scientific use of such results here is to identify a feasible mechanism and its predicted signatures, not to adopt their success as a direct measurement of an ancient onset date.^[ref]

Connect a model regime to a date through independent evidence. Couple the experiments to candidate planetary cooling histories, then predict observable pressure-temperature-time paths, magmatic compositions, deformation patterns and recycling rates for each regime. If independently dated rocks constrain a thermal or compositional state, ask which regimes can reproduce that state and the observed geological associations. Repeat over uncertain initial conditions and constitutive parameters, check numerical resolution, and compare two- and three-dimensional results where geometry matters. Models that prescribe a descending slab, fix plate velocity, or insert a weak boundary may be excellent forward models of a working subduction zone but do not demonstrate spontaneous initiation.

The result is a conditional statement: persistent plate tectonics is physically feasible under a range of conditions, and dated geological observations may favor that range at a particular time. Uncertain cooling histories, several feasible regimes and uncertain rheology prevent numerical feasibility alone from dating onset. The strongest constraint comes when the same parameter range explains multiple independently dated observations and alternative regimes fail those tests.

Let $n=N/\mathcal{V}$ and $c_e=C_e/\mathcal{V}$ be the carrier density and electronic heat capacity per unit volume. In isotropic kinetic theory, the angular average $⟨v_x^2⟩=v_F^2/3$ gives the electronic thermal conductivitywhere ${\ell }_Q=v_F{\tau }_Q$ is the mean free path associated with the heat-current relaxation time. With $v_F^2=2E_F/m^{\ast }$,The Drude model gives the electrical conductivity $\sigma =n{e}^2{\tau }_{\sigma }/m^{\ast }$, where ${\tau }_{\sigma }$ is the momentum relaxation time. ThereforeWhen the same approximately elastic collisions relax heat and charge currents, ${\tau }_Q={\tau }_{\sigma }$. Using $A={\pi }^2/2$ gives the Wiedemann-Franz law$L_0$ is the Lorenz number. Degenerate electrons move at the Fermi velocity, not the classical thermal velocity; using classical equipartition would give the wrong coefficient.

The Wiedemann-Franz law applies to a normal degenerate metal when the same itinerant electrons carry heat and charge, scattering is approximately elastic over the thermal shell, and the transport probability varies slowly across the energy interval $k_{B}T$ around the Fermi level. The comparison uses the electronic thermal conductivity $K_e$, with any lattice contribution removed.

It works particularly well as $T\to 0$ when elastic impurity or defect scattering dominates. It also often holds approximately for ${\Theta }_D\ll T\ll T_F$, where thermally populated phonons act as a nearly elastic classical scattering field. Room temperature need not lie above the Debye temperature for every metal.

At intermediate temperatures, inelastic small-angle electron-phonon scattering can relax heat much more efficiently than charge. Then ${\tau }_Q\ne {\tau }_{\sigma }$ and $K_e/(\sigma T)$ can be smaller than $L_0$. Strong energy dependence of scattering, additional neutral heat carriers, or a transition out of the normal metallic state also invalidates the simple derivation. An observed deviation alone does not identify which assumption failed.

Assume a stationary temperature profile, spatially uniform mass density and radioactive heat generation rate, constant thermal conductivity, isotropic solid thermal conduction, and blackbody emission into effectively cold surroundings. The volumetric source isTotal radiogenic power equals the surface blackbody radiation loss. The Stefan–Boltzmann law givesFor steady radiogenic heating of a conducting sphere, the interior heat equation isRegularity at the centre removes the $1/r^2$ integration term and yields $dT/dr=-q_{v}r/(3\kappa )$. Matching $T(R)=T_S$ givesThe high central temperature is compatible with a cold surface because finite thermal conductivity requires a large temperature gradient to carry the accumulated interior heat. A nonzero background radiation temperature $T_{\mathrm{bg}}$ replaces $T_S^4$ by $T_S^4-T_{\mathrm{bg}}^4$ in the surface balance. The calculation neglects convection, stored transient heat, changing radioactive abundances, phase changes and any substantial external irradiation.

By Fourier's law, $q_x=-k{T}_x$. With nonzero thermal conductivity $k$, zero outward heat flux at both endpoints gives homogeneous Neumann boundary conditions:The outward normal has opposite signs at the two ends, but a zero-flux condition gives zero derivative at each.

At uniform temperature, define the material Peltier coefficient by transported heat current along positive conventional current:Its units are $W/A=J/C=\mathrm{V}$. In one dimension, $q=\Pi j-\kappa dT/dx$, with electric current density $j$ and thermal conductivity $\kappa$. Here heat excludes the electrochemical potential contribution to total energy flux.^[ref]

Now orient positive conventional current $I$ from material A into material B, at a common junction temperature $T$. Define ${\dot{Q}}_{\mathrm{abs}}>0$ as heat taken from the junction's thermal surroundings. The outgoing transported heat minus the incoming transported heat givesThe local released heat is ${\dot{Q}}_{\mathrm{rel}}=-{\dot{Q}}_{\mathrm{abs}}$. Authors defining a junction coefficient by released heat therefore use the opposite sign. Exchanging A and B reverses the coefficient; reversing $I$ at fixed A and B reverses the junction's thermal response. An individual material coefficient and a junction difference are different quantities.

For an original stipulated example, choose $S_A=-100\mu V/K$ and $S_B=+150\mu V/K$ at $T=300\mathrm{K}$. The Kelvin thermoelectric relation gives ${\Pi }_A=-0.030\mathrm{V}$, ${\Pi }_B=+0.045\mathrm{V}$. A positive $2\mathrm{A}$ flowing A to B absorbs $(0.075\mathrm{V})(2\mathrm{A})=0.150\mathrm{W}$. A separate $0.010\Omega$ junction electrical resistance would dissipate $0.040\mathrm{W}$ for either direction. If all that dissipation heats the junction, its net absorbed power is $0.110\mathrm{W}$ forwards and $-0.190\mathrm{W}$ backwards. These numbers illustrate a chosen thermal accounting boundary, not measured material data or a universal allocation of contact heating.

A radiometer thermal link is the intentional path for thermal conduction from receiver to reference block. Near a chosen operating point, write $Q_{\mathrm{link}}\simeq G(T-T_b)$ with local conductance $G$ in $\mathrm{W}{\mathrm{K}}^{-1}$. A weaker link increases the steady temperature rise per added power, $\delta T\simeq \delta P/G$, but also increases the uncontrolled relaxation time $\tau =C/G$ when the link dominates other heat transfer. These are local linear approximations: actual thermal conductivity depends on temperature, and active feedback control introduces its own response. NIST explicitly describes the tradeoff between weaker links and smaller receiver heat capacity.^[ref]

A single-category model represents local ice by one thickness $h$ and one surface temperature $T_s$. With basal freezing temperature $T_m$, atmospheric temperature $T_A$, thermal conductivity $k$ and atmospheric heat transfer coefficient ${\lambda }_A$, quasistatic thermal conduction giveswhile $T_s<T_m$. If $F_O$ is upward ocean heat flux, the Stefan condition gives $\rho L\dot{h}={\lambda }_A(T_m-T_A)/(1+{\lambda }_{A}h/k)-F_O$. During surface melting, $T_s=T_m$ and $\rho L\dot{h}=-{\lambda }_A(T_A-T_m)-F_O$ when $T_A\ge T_m$. Thickness is constrained to be nonnegative, and the open-water energy balance applies after complete melt.

A one-dimensional thermal conduction model with constant thermal conductivity, source $q^{\prime \prime \prime }(x)=q_0{e}^{-\mu x}$ and prescribed face temperatures $T(0)=T(a)=T_0$. Its steady state solution isPositive internal heating makes the profile concave. Its maximum lies toward the more strongly heated face, at $x_m={\mu }^{-1}\ln [\mu a/(1-e^{-\mu a})]$. The uniform-source limit is $T-T_0=q_{0}x(a-x)/(2k)$.

Snow is solid atmospheric precipitation made of ice crystals and their aggregates. Crystals grow in cold clouds through transfer of water vapor to ice; collisions and the collection of supercooled droplets can further change their form. Air temperature and moisture conditions influence the resulting particles.

A settled snowpack contains both ice and pore space, so its bulk mass density and thermal conductivity differ from those of solid glacier ice. Compaction and recrystallization can turn old snow into firn. Its high albedo also changes the absorbed solar energy, while its insulating effect alters heat exchange with the ground or underlying sea ice. These optical and thermal roles can act in different directions in a surface energy budget.

In one-dimensional steady heat conduction without internal heat sources, conservation of energy gives constant heat flux $q=-k{T}^{\prime }$. Through a layer of thickness $L_j$ and positive thermal conductivity $k_j$, the temperature drop is $q{L}_j/k_j$. If adjoining layers are in perfect thermal contact, temperature and heat flux are continuous at each interface, so summing the layer drops proves the displayed formula. Temperature is piecewise linear and its slope may jump when thermal conductivity changes.

For constant volumetric heat generation rate $q_v$ and thermal conductivity $\kappa$, spherical steady thermal conduction obeys $\kappa r^{-2}(r^2{T}^{\prime }{)}^{\prime }+q_v=0$. Regularity at the centre and a specified surface temperature give the displayed parabolic profile. For a black surface radiating into negligible background radiation, global energy conservation gives ${\sigma }_{\mathrm{SB}}{T}_S^4=q_{v}R/3$. The central excess scales as $R^2$, while the surface temperature scales only as $R^{1/4}$.

For a long homogeneous circular wire of radius $a$, uniform volumetric heating $H$, constant thermal conductivity $\kappa$ and negligible end losses, the steady heat equation is $\kappa {\rho }^{-1}(\rho T^{\prime }{)}^{\prime }+H=0$. Regularity at the axis excludes a logarithm term, giving $T=T_s+H(a^2-{\rho }^2)/(4\kappa )$. If the surface obeys Newton's law of cooling, $-\kappa T^{\prime }(a)=\gamma (T_s-T_{\mathrm{air}})$, then $T_s=T_{\mathrm{air}}+Ha/(2\gamma )$. Here $\gamma$ has units $\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$ and $H$ has units $\mathrm{W}{\mathrm{m}}^{-3}$. For uniform electric current $I$ and constant electrical conductivity $\sigma$, Joule heating gives $H=I^2/(\sigma {\pi }^2{a}^4)$.

The 1935 primary paper The Spin-Isomerization of Hydrogen at Paramagnetic and Diamagnetic Surfaces, pp. 1251–1252, describes evacuated catalyst samples cooled before admitting molecular hydrogen. Gas samples were analyzed in a glass cell in liquid air: a heated platinum wire and Wheatstone bridge converted changes in thermal conductivity into electrical resistance readings, calibrated against normal hydrogen.^[ref]

Its comparison of oxide surfaces supplied a control: specimens exhibiting strong paramagnetism converted the gas much faster than the comparison exhibiting diamagnetism despite all adsorbing molecular hydrogen. Contact time, pressure, and temperature were recorded. This supports a role for surface magnetic properties beyond adsorption alone; it does not identify every microscopic surface site or directly measure nuclear-spin orientations.^[ref]

The inference chain is preparation, contact with a specified surface, then a calibrated bulk transport readout of composition. It is distinct from a direct calorimetry measurement of heat capacity: transport and heat capacity can both depend on the spin populations, but they are different measured quantities. Changing the catalyst's own surface structure, introducing magnetic impurities, or changing the detector's gas pressure can confound an uncalibrated comparison. The free-rotor partition functions predict equilibrium populations and stored energy, not the conversion rate on a particular surface.