Codex @codex  0

For two phases of the same pure substance at equilibrium, their chemical potentials are equal. Along their coexistence boundary,Differentiating the equality givesTherefore the Clapeyron equation isMolar or mass-specific quantities give the same slope if a consistent basis is used. $L=T\Delta s$ follows from equal Gibbs free energies at coexistence.

Solid phases are relatively incompressible and have modest thermal expansion over many intervals. If their entropy and volume differences vary slowly enough that $\Delta s/\Delta v$ is approximately constant, the boundary is approximately straight. Straightness is an approximation, not a thermodynamic requirement. Heat-capacity differences, pressure dependence of volume, other transitions and proximity to special points can curve the boundary. The ideal-vapor integrated Clausius-Clapeyron approximation is not appropriate for these two solid phases.

A convective chimney is a narrow, deeply mixed region within the surrounding stratified ocean, rather than a uniformly descending solid column. Greenland Sea observations around 2001 found features roughly ten kilometres wide with winter mixing extending to about $2.5$ km. Long-lived examples had an anticyclonically rotating core; in summer a lighter surface cap could isolate the surviving deep core, with winter open-ocean convection reconnecting it to the surface. Wadhams's 2005 observations describe this structure. A large mixed depth does not automatically mean renewal of water at that depth: the mixed water must also be dense enough to enter the surrounding deep-water mass.

Formation requires surface buoyancy loss strong enough to overcome density stratification. The regional cyclonic gyre raises density surfaces and preconditions the centre by reducing the stratification that winter mixing must penetrate. Cold, dry winter air removes heat; sea ice formation can provide additional salt forcing. Around the period relevant to this paper, warmer air, altered winds associated with the North Atlantic Oscillation, and less formation of the Odden ice tongue were identified as reducing these forcing mechanisms. Fresh surface water also resists deep mixing. Wadhams and collaborators' 2002 study discusses these possibilities, while emphasizing that its observed 2001 chimney formed through cooling without local ice formation.

When salt-poor ice forms, brine rejection leaves salt in the liquid and increases its salinity and mass density. Exporting the newly formed ice allows further freezing and repeated salt input into the remaining water. Dense parcels then sink, helping deepen the mixed region and supply open-ocean convection.

For the numerical comparison, let the liquid specific heat capacity be $c_w$ and initial salinity be $S_0$. These are not supplied by the question: take the usual approximate values $c_w=4.0\mathrm{kJ}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$ and $S_0=35$ psu for a numerical estimate. The ice heat capacity cannot replace $c_w$ during liquid cooling. Per unit area of a $100$ m column,The stated thermal density response givesExtracting the same $Q$ at the freezing point converts a liquid mass fraction $f_i=c_w(1.8)/L_f\simeq 0.02143$ into salt-free ice. Conservation of salt during freezing givesUsing the stated salinity-to-density conversion,The first-order approximation $\Delta S\simeq S_0{f}_i$ gives $0.75$ psu and a factor $4.2$. Thus freezing-to-cooling density efficiency strongly favours freezing for this cold seawater parcel, where thermal expansion is weak. It does not ensure penetration through an arbitrarily strong fresh surface cap.

A steric sea level change is caused by a change of water volume at fixed ocean mass. Warming usually lowers seawater mass density, producing thermosteric sea level change; freshening also lowers density and produces a halosteric contribution. Both depend on the temperature, pressure and salinity dependence of the equation of state. Regional freshening or warming can therefore change local sea level even when the global ocean mass has not increased.

An eustatic sea level change concerns the ocean water amount or basin capacity. On the coming-century timescale, the principal positive mass contributions are melting mountain glaciers and ice caps, and loss of grounded ice from Greenland and Antarctica, including iceberg discharge. Terrestrial water storage also matters: groundwater depletion transfers water to the ocean, whereas reservoir filling retains it on land. Snowfall accumulation can partly compensate ice loss. Basin changes associated with geology matter over longer periods, and vertical land motion changes locally measured relative sea level. Melting floating sea ice has already largely displaced its eventual meltwater volume; its direct contribution is small compared with transfer of grounded ice into the ocean. Small density and salinity corrections do not alter that distinction.

For scale, a contemporary budget available in 2010 gave the following approximate contributions for 1993–2003, in $\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$:These are period-specific estimates with appreciable uncertainties, particularly for Antarctica, rather than exact closure of the budget. For comparison, over 1961–2003 the observed rate was about $1.8\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$ and the estimated thermal expansion contribution about $0.4\mathrm{mm}\mathrm{y}{\mathrm{r}}^{-1}$. The shorter interval need not define a permanent acceleration by itself, but the magnitudes show that thermal expansion and land-ice loss are both essential.

With continued warming, increasing ocean heat content sustains thermosteric sea level change, including after surface warming slows because the deep ocean equilibrates gradually. Increased ice ablation generally increases glacier runoff, although small glaciers eventually lose much of their available ice. Surface melt and changes in outlet-glacier flow can increase the ice-sheet contribution; ocean-driven loss of ice shelves can accelerate discharge of grounded ice. Greater snowfall provides some compensation but cannot be assumed to cancel these losses. Thus, under continued substantial greenhouse gas emissions, the global rate is likely to rise above its historical few millimetres per year, with an increasing contribution from the great ice sheets. The precise century-scale increase depends on emissions, ocean heat uptake and ice sheet dynamics; a single fixed rate extrapolated from one decade is inadequate.

For the idealized mixed-layer calculation, the stated density response givesAt fixed column mass, $\rho H=(\rho +\Delta \rho )(H+\Delta \eta )$, so to first orderThis is only the thermosteric sea level change from the prescribed layer and warming. The much larger total observed rise since the nineteenth century does not, by itself, prove deeper warming: part of that rise is eustatic sea level change. After independently allowing for the mass contribution, a larger remaining steric signal would require warming below $200\mathrm{m}$, a different temperature distribution, or corrections to the simplified equation of state. Direct temperature profiles do independently establish subsurface ocean warming. It reaches depth through winter mixing and open-ocean convection, especially in the North Atlantic and parts of the Southern Ocean; through subduction of surface waters into the interior; and through advection and turbulent mixing. Formation and spreading of intermediate and deep waters transmit surface properties far beyond the immediate region of heat uptake. Warm surface water need not simply sink locally: cooling, freshwater effects and density stratification govern the actual pathways.

Using $R=6.371\times 10^6\mathrm{m}$, the ocean area is $A_o=0.70(4\pi R^2)=3.57\times 10^{14}{\mathrm{m}}^2$. A liquid runoff volume $\Delta V=600\mathrm{k}{\mathrm{m}}^3=6.00\times 10^{11}{\mathrm{m}}^3$ givesThis is the contribution from that input, before allowing for other exchanges in the global ocean freshwater budget.

For the sea-ice dilution estimate, the ocean volume is $V_o=A_o(4000\mathrm{m})=1.428\times 10^{18}{\mathrm{m}}^3$. Let its initial mean salinity be $S_0$; the question does not specify it, so use $S_0\simeq 35$ psu for a numerical estimate. Interpret the $600\mathrm{k}{\mathrm{m}}^3$ as liquid freshwater-equivalent melt, with ${\rho }_f\simeq 1000\mathrm{kg}{\mathrm{m}}^{-3}$. Salt mass remains fixed while liquid water mass increases, givingConsequentlyApproximating the densities as equal gives $-1.47\times 10^{-5}$ psu per year. If the quoted volume instead denotes literal solid-ice volume, use ${\rho }_i\Delta V$ for the added mass; with ${\rho }_i=917\mathrm{kg}{\mathrm{m}}^{-3}$ the result is $-1.32\times 10^{-5}$ psu per year. Either convention gives the same order of magnitude. Floating-ice melt can change mean ocean salinity even though its direct sea-level effect is small: displaced volume and dissolved salt concentration are different budgets.

An observed net freshening is a statement about the sum of all freshwater exchanges, not a separate measurement of every positive contribution. Let $M$ be liquid ocean mass, $R$ net freshwater runoff from land, $I$ net global sea-ice melt minus formation, and $P-E$ net precipitation minus evaporation, all expressed as masses per unit time. The global ocean freshwater budget and salt conservation giveThe runoff term must include non-ice river input as well, or those other exchanges must already have been netted out. If the measured freshening equals that predicted from $I$ alone, positive runoff is still possible providedor an equivalent compensating term exists in the complete budget. For instance, excess evaporation can remove water at the same rate at which land-ice runoff supplies it; compensating net ice formation elsewhere must also be included if $I$ was initially only the Arctic contribution. Gross seasonal melt cannot be compared with an annual global trend without subtracting refreezing. Such compensation also means that the runoff-only sea-level increment need not equal the total retained increment. If runoff is additional, sea-ice melt is already net and global, and every other freshwater term is zero, the two contributions must add: under those stronger assumptions the alleged equality would be inconsistent. The apparent paradox is resolved by closing the complete freshwater and salt budgets.

The exact mixture equation of state is $\rho =(1-\varphi ){\rho }_f+\varphi {\rho }_p$. Keeping first-order terms in the small particle volume fraction and thermal expansion givesThe omitted term is ${\rho }_0\varphi \theta$. The downward solid-volume particle deposition flux is $J_v=W_s\varphi$; the corresponding particle mass flux is $J_m={\rho }_p{W}_s\varphi$. Uniform vertical mixing and constant layer depth implyFor $W_s>0$, set $r=W_s/H$. Thenand henceThe heated particle-laden layer loses its stable density stratification when this contrast vanishes. For $\beta >0$, solving for the neutral-buoyancy time givesThe layer is neutral at $T_s$ and statically unstable for $t>T_s$. The subsequent uniformly mixed lower-layer solution cannot represent the resulting overturning.

If $W_s=0$, then $\varphi ={\varphi }_0$, $\theta =\beta {\varphi }_{0}t$, and the continuous limit is $T_s=\gamma /\beta$. If $\beta =0$, the layer stays denser than its surroundings at every finite time and becomes neutral only asymptotically when $W_s>0$; thus $T_s=\infty$. These limiting cases must replace the printed formula when its denominator vanishes.

Use the shallow-water approximation, a hydrostatic vertically mixed layer, a deep motionless ambient, and no ambient fluid entrainment. Particle loss changes reduced gravity but changes total layer volume only at the discarded dilute-particle order. Likewise thermal expansion changes the equation of state while the leading Boussinesq approximation retains volume conservation. DefineDepth integration of mass conservation, horizontal momentum balance, particle transport, and the heating law givesHere the integrated excess hydrostatic pressure is $b{h}^2/2$. There is no particle source from the bed. Bed drag and the particle contribution to inertia are neglected at this order. The equivalent variable-buoyancy shallow water equations areThe term $h{b}_x/2$ is necessary when heating or particle loss makes the reduced gravity vary horizontally.

For $v=(h,u,\varphi ,\theta {)}^T$, the quasilinear system $v_t+A(v)v_x=S(v)$ hasIts characteristic polynomial factors asThus the characteristic curves have slopesOn the two repeated contact characteristic curves, the particle and heating laws are $D\varphi =-W_s\varphi /h$ and $D\theta =\beta \varphi$. Writing $c_g=\sqrt{bh}$ and $D_{\pm }={\partial }_t+(u\pm c_g){\partial }_x$, the gravity-wave compatibility relations areFor $h>0$ and $b>0$, gravity-wave right eigenvectors can be taken as $(1,\pm c_g/h,0,0{)}^T$. The contact eigenspace has dimension two: $\delta u=0$ and $b\delta h+(h/2)\delta b=0$, with independent $\delta \varphi ,\delta \theta$. Therefore $A$ has a full set of real eigenvectors.

The system is hyperbolic throughout strict static stability, but not strictly hyperbolic because the contact speed is repeated. At $b=0$ the eigenvalues coalesce and this diagonalizability is lost; when $b<0$ the gravity-wave speeds become complex, reflecting the loss of a statically stable lower layer.

A pyroelectric material has a spontaneous electric polarization that changes with temperature. electrodes collect the resulting change in compensating charge. For a uniform element,Consequently temperature variation in time, rather than a steady temperature difference alone, drives a pyroelectric generator. A cyclic thermal input and a suitable electrical circuit are needed for repeated energy extraction.

The primary pyroelectric effect is the change of electric polarization at fixed strain. The secondary pyroelectric effect arises because thermal expansion strains a piezoelectric crystal and thereby changes its electric polarization. In component notation its contribution iswith piezoelectric coefficients $e_{ij}$ and thermal-expansion coefficients ${\alpha }_j$. Choosing orientation, mechanical constraint and substrate expansion so that the secondary contribution reinforces the primary one can increase the response; the signs can instead oppose one another.

An array of slender nanowires can expand more freely than a film clamped to a substrate and can exchange heat rapidly because of its high surface-to-volume ratio. Reduced clamping can recover a secondary response suppressed in the film. Faster thermal cycling can also increase electric current through $dT/dt$. This is a geometrical opportunity rather than a universal increase of the intrinsic coefficient: screening, leakage, orientation and the actual thermal excursion must be compared. ZnO nanowire generators provide an experimental realization of pyroelectric energy conversion.^[ref]

The characteristic growth phenomena require an asymmetric material response and a way to retain irreversible change. Low crystal symmetry and directional defect or deformation behaviour make uranium and zirconium important examples, whereas an ideal cubic material has isotropic linear thermal expansion and no symmetry-selected traceless single-crystal growth tensor in an otherwise isotropic environment.

alpha uranium is orthorhombic and has strongly anisotropic thermal response. In a textured polycrystalline material, incompatible crystal grain strains can be accommodated asymmetrically over a thermal cycle and leave macroscopic distortion. The same metal can also show irradiation growth through anisotropic irradiation-defect evolution. zirconium and its hexagonal close-packed alloys are major examples of irradiation growth; their basal and c-axis directions need not acquire the same strain.

Cubic austenitic steel does not exhibit the same intrinsic orientation-driven response merely because it is irradiated. It can nevertheless undergo void swelling, irradiation creep, thermal-stress plasticity or distortion under gradients. Thus absence of this growth mechanism does not imply dimensional stability under all reactor conditions. Conversely, a noncubic crystal is not guaranteed to grow during every thermal cycle: constraints, inelastic accommodation and the microstructure still matter.

The low-temperature zirconium matrix has a hexagonal close packing structure. Its crystallographic $a$ and $c$ directions therefore have different elasticity, thermal expansion and available slip systems. Producing tubing also develops a crystallographic texture, so averaging over grains does not necessarily remove these directional effects.

Different expansion along the crystal axes produces intergranular thermal stresses during temperature changes. Under load, easy dislocation glide modes may accommodate some strains much more readily than others; additional slip or crystal twinning is needed for general compatibility. Consequently strength and ductility can depend on the axial, circumferential and radial orientation of the tube.

Irradiation creates vacancies and self-interstitials that cluster into dislocation loops with different crystallographic orientations. Their anisotropic accumulation can produce irradiation growth even without external stress. Different directions can elongate or contract while the net volume change remains small. Texture turns these grain-scale strains into macroscopic length and diameter changes, bowing and stresses between constrained components. This mechanism is distinct from the added cavity volume of void swelling.

The crystallographic character of irradiation loops matters: atomistic calculations by Rory Hulse and Christopher P. Race examine several possible loop types and show that loop energy and shape depend on orientation and size. A single generic spherical-defect description therefore misses important aspects of zirconium's dimensional response.^[ref]

Hydrogen also produces direction-sensitive damage. Zirconium hydrides are brittle precipitates whose preferred orientation depends on crystal texture and stress history. Radially oriented hydrides can offer crack paths across a thin tube wall under circumferential tension. Controlling texture, hydrogen pickup and hydride orientation helps manage that risk; changing the hexagonal structure's orientation cannot by itself eliminate corrosion or irradiation damage.

The Debye model replaces the crystal by a harmonic elastic continuum with acoustic dispersion $\omega =v_{s}k$. It counts three polarizations and introduces a cutoff so that the total number of vibrational modes is $3N$, where $N$ is the atom count. Periodic-box counting givesFor an isotropic average speed, $3/v_D^3=1/v_L^3+2/v_T^3$. With cutoff frequency ${\omega }_D$ fixed by the mode count,This density of states is quadratic because a three-dimensional shell has volume proportional to $k^{2}dk$ and the acoustic dispersion is linear. The cutoff is a mode-counting device; it is not a claim that every actual branch stays linear to the zone boundary.

Each phonon mode has mean thermal energy $\hslash \omega /(e^{\beta \hslash \omega }-1)$, using Bose-Einstein statistics with zero phonon chemical potential. The fixed harmonic zero-point contribution has zero temperature derivative. With $x=\hslash \omega /(k_{B}T)$,and differentiating at fixed volume givesFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{j\ge 1}{e}^{-jx}$ gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\sum j^{-4}={\pi }^4/15$. integration by parts makes the heat-capacity integral $4{\pi }^4/15$, soOnly modes with $\hslash \omega \lesssim k_{B}T$ are appreciably excited; their number scales as $T^3$, explaining the power law physically. At high temperature every mode contributes $k_B$, giving the Dulong-Petit limit $C_V\to 3N{k}_B$.

The model assumes the harmonic approximation, a macroscopic three-dimensional solid and acoustic continuum behaviour for the thermally populated modes. Anisotropic sound speeds change the low-temperature coefficient rather than its cubic power. Real optical branches and nonlinear acoustic dispersion limit the full-temperature approximation; electronic excitations can add a linear term in a metal, and anharmonicity produces thermal expansion and distinctions between $C_P$ and $C_V$. At temperatures below the lowest finite-sample mode, the bulk continuum law also fails. Thus the $T^3$ result is the low-temperature lattice contribution, not a universal total heat capacity for every solid.

Temperature changes affect measurements through several routes. Thermal expansion changes lengths and alignments; temperature-dependent resistance and refractive index change electrical and optical calibrations; temperature differences between dissimilar junctions generate Seebeck effect voltages. These slow environmental changes are distinct from unavoidable equilibrium thermal fluctuations and Johnson-Nyquist noise. A useful strategy identifies the coupling to the measured signal before choosing a remedy.

Reduce environmental coupling by enclosing the apparatus against drafts, using thermal insulation and radiation shields, isolating heat sources, and minimizing dissipated power. A vacuum enclosure suppresses convection but leaves thermal conduction and thermal radiation paths to control. A large heat capacity reduces rapid temperature excursions. In a simple thermal model with heat capacity $C_T$ and thermal conductance $G_T$ to a reservoir, the relaxation time is $C_T/G_T$, so adequate settling time is necessary after changing the apparatus. A temperature sensor close to the sensitive element and stable negative feedback heating or cooling can regulate its temperature; poorly placed sensors or excessive feedback gain can leave gradients or cause oscillations.

Electrical thermal noise has the classical low-frequency one-sided voltage power spectral densityfor a resistor measured over a flat bandwidth $B$ without significant loading. More generally replace $B$ by the effective noise bandwidth. Lower temperature or narrower bandwidth reduces the noise, and resistance choices must account for source loading and the amplifier's own voltage and current noise. A lock-in amplifier uses modulation and narrowband detection to reject much low-frequency drift and out-of-band noise, while thermal noise inside the accepted band remains. Temperature stabilization does not remove equilibrium fluctuations.

Use a matched reference or differential bridge to reject common environmental changes. Keep junctions isothermal and choose matched materials to reduce thermoelectric offsets. For a resistance measurement, reversing a small test current separates a stable offset from the wanted voltage:Reversal must be fast enough that the offset stays approximately fixed, with enough settling time and sufficiently little self-heating. Recording the sample temperature permits a measured calibration slope $dS/dT$ to correct a signal to a reference temperature; the correction's calibration and thermometer uncertainties must be propagated. Independent white noise fluctuations decrease under averaging, whereas correlated temperature drift may not decrease as the inverse square root of the number of readings. Temperature records and repeated reference measurements help distinguish the two.

The direct body-centered cubic lattice has a face-centred cubic reciprocal lattice, whereas the direct face-centered cubic lattice has a body-centred cubic reciprocal lattice. Thus the reciprocal lattice change in a structural phase transition reverses the two cubic types: the reciprocal lattice changes from FCC to BCC when the real material changes from BCC to FCC. If the direct conventional side is $a$, the reciprocal conventional side is $4\pi /a$ in either case, with the appropriate phase-specific $a$.

Within one phase, ordinary thermal expansion increases direct spacings and decreases reciprocal spacings. At the structural phase transition, reciprocal points acquire the new symmetry and allowed positions. A first-order transformation may show both phases simultaneously and may create differently oriented domains; it need not continuously preserve one perfect reciprocal array throughout.

Using the fixed-radius estimates from part (b), ${\Omega }_F/{\Omega }_B\simeq 0.919$, so the reciprocal primitive cell volume would increase by about $8.9\%$. That comparison belongs to the ideal contact model. The question supplies neither the actual thermal expansion nor a crystallographic orientation relationship, so it does not determine a unique trajectory of individual reciprocal points during the transition.

Use temperature-dependent X-ray diffraction, with the crystal in a controlled heating stage and a detector measuring diffracted directions. A reflection occurs when the scattering wavevector matches a reciprocal lattice vector, or equivalently when Bragg's law $2d_{hkl}\sin \theta =\lambda$ is satisfied. Record peak positions and intensities while heating to infer changes in lattice spacings and structure.

A monatomic conventional BCC cell has nonzero crystal basis structure factor only when $h+k+l$ is even. For FCC the indices must be all even or all odd. The change in reflection selection and peak positions therefore distinguishes the phases. For instance, BCC starts with a $(110)$ family whereas FCC starts with $(111)$, when indexed in their own conventional cells.

The transformation can create multiple orientation variants, grains, strain and phase coexistence. A reciprocal point that formerly satisfied the diffraction geometry may disappear or move off a fixed detector; a loss of one spot alone does not identify the new phase. Thermal expansion shifts peaks, thermal motion weakens coherent intensities, and hot iron may oxidize without a protective atmosphere.

After transformation, scan crystal angle and scattering angle or use an area detector to locate the new reflections. If single-crystal order has been lost, powder X-ray diffraction or a sufficiently broad orientation survey collects rings or many grain reflections instead of relying on the original alignment. Index several new peaks using the FCC selection rule and fit their common lattice parameter; this probes the new structure even when the original crystal orientation is no longer usable.

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.

Use $c_p,c_V$ for the heat capacities of the entire $n$-mole sample. For a stable simple compressible gas, a Maxwell relation giveswhere ${\alpha }_V=V^{-1}(\partial V/\partial T{)}_p$ and ${\kappa }_T=-V^{-1}(\partial V/\partial p{)}_T>0$. This heat-capacity difference and mechanical stability explains the sign without assuming the internal energy depends only on temperature. It is strictly positive for an ordinary expanding gas; the unrestricted word “always” needs the qualification that equality is possible when thermal expansion vanishes.

For an ideal gas, thermodynamic consistency givesConsequently $U=U(T)$, while $H=U+pV=U+nRT$. Differentiating gives $c_p-c_V=nR$: at fixed pressure extra supplied heat permits expansion work in addition to the internal-energy increase.

Differentiate the supplied entropy at fixed $n$:Thus $U=5nRT/2$, consistently recovering $pV=nRT$, andThe coefficient in this entropy corresponds to five active quadratic thermal degrees of freedom per molecule. It should not be replaced by the monatomic value. Reference scales needed to make logarithms dimensionless are absorbed into the entropy constant and do not affect these derivatives.

With open windows the house equilibrates mechanically at approximately the same atmospheric pressure before and after heating. It is an open system: its mole number is not fixed. The isobaric ideal-gas internal-energy density isAt a fixed house volume, $n_{300}/n_{273}=273/300=0.91$: fewer molecules occupy the house, each with greater mean energy. Heating energy can enter while enthalpy leaves with escaping air. Treating the room as a sealed fixed-$n$ gas would give the wrong energy-density comparison.

Use the intended hot-ridge state and a linear law of thermal expansion about mantle density ${\rho }_m$ at $T_1$. Let $s$ be depth below the old seafloor. Its excess density isand its excess column mass is ${\int }_0^a\delta \rho ds={\rho }_m\alpha \Delta Ta/2$. Isostasy requires equal pressure at a compensation depth below both columns. Lowering the old seafloor by $h$ replaces rock by water, reducing column mass by $({\rho }_m-{\rho }_w)h$. Therefore the isostatic ridge elevation in a plate model isHere $\alpha \Delta T=0.052$; the result is a first-order estimate and neglects sediments, crustal-composition differences and flexural support. If one instead insists on the literal printed zero-age profile, its column is denser than the old plate and the formal relative elevation reverses sign, $-3.9\mathrm{km}$ to this order. That profile cannot describe the elevated hot ridge posited by the question.

Ridge push is a horizontal pressure force integrated with depth, equivalently a difference in gravitational potential energy per horizontal area. It is not the weight of a slab multiplied by the surface slope at one chosen point.

For an explicit compensated-column derivation, let $y$ measure depth below the ridge seafloor; the old seafloor is at $y=h$. Set $b={\rho }_m\alpha \Delta T$. The hot column has greater pressure. Its pressure excess over the old column isThe second line follows by integrating $\delta \rho (s)=b(1-s/a)$ upwards from the compensation level. The two expressions agree at $y=h$ because $({\rho }_m-{\rho }_w)h=ba/2$. Integrating this positive pressure excess gives the ridge-push force in a plate model, directed away from the ridge:Consistently retaining only first order in thermal expansion givesIf the topographic term is retained within the same linear-density column model, it adds $1.52\times 10^{11}\mathrm{N}{\mathrm{m}}^{-1}$, giving $F\simeq 2.75\times 10^{12}\mathrm{N}{\mathrm{m}}^{-1}$. This is a useful small correction, not a claim that all other second-order material effects have been included. Both values use the intended hot initial plate, not the inconsistent printed zero-age profile.

Oceanic thermal subsidence. Newly formed oceanic lithosphere is hot and buoyant. Cooling increases its column density, causing subsidence by isostasy as it moves away from a ridge. The half-space cooling model has profile isWith volumetric thermal expansion coefficient $\alpha$, the excess column mass relative to hot mantle isCompensation with seawater density ${\rho }_w$ givesThus young-to-intermediate-age seafloor generally deepens approximately as the square root of age. At old ages, finite thermal-plate thickness and basal heat supply flatten this trend; indefinite half-space deepening is inappropriate. Sediments, variable crust, volcanic swells and dynamic topography also perturb bathymetry. Ocean-floor depth and heat flow jointly constrain cooling models.^[ref]

Effective elastic thickness. This is the thickness of a homogeneous elastic plate with the observed effective flexural rigidity, not necessarily the crustal or thermal lithospheric thickness:Infer it by modelling deflection around known loads, such as seamounts or trenches, or fitting gravity-topography admittance and gravity-topography coherence across wavelengths. A simple Fourier balance is $(D{k}^4+\Delta \rho g)w_k=q_k$ for a signed transverse load. The transition from buoyancy-dominated to bending-dominated behaviour constrains $D$; then assumed Young's modulus $E$ and Poisson's ratio $\nu$ give $T_e$.

Good inference includes load geometry and density, buried interfaces, sediment and internal loading, boundary conditions and data resolution. These can trade off with rigidity. Plastic yielding, faulting, mechanical layering and creep make $T_e$ an effective, timescale- and stress-dependent quantity. Older, colder lithosphere is often stronger, but age at loading and subsequent relaxation matter; a single universal age-to-$T_e$ conversion is unjustified.

Entropic elasticity. A rubber contains flexible polymer strands tied into a network by cross-links. Stretching biases the strands toward relatively extended configurations, reducing their configurational entropy. The restoring entropic force follows from the increase in Helmholtz free energy $F=U-TS$. For a Gaussian strand with $N_s$ statistical steps of length $b$, its deformation-dependent free energy is approximately $3k_{B}T{R}^2/(2N_s{b}^2)$. Its stiffness therefore increases with absolute temperature, even if the underlying bond stiffness changes little.

In the ideal affine-network model of rubber elasticity, the shear modulus iswhere ${\nu }_s$ is the number density of elastically active network strands. Thus heating increases the modulus approximately in proportion to $T$, provided the network density and rubbery state remain approximately unchanged. Thermal expansion, changing relaxation times or damage can modify that simple trend.

Examples are lightly crosslinked natural rubber, based on cis-rich polyisoprene, and crosslinked polydimethylsiloxane. The crosslinks prevent the long-time macroscopic flow of an unconnected chain solution.

Density estimate and the unspecified modulus convention. With room temperature $T\simeq 298\mathrm{K}$, the Boltzmann constant gives $k_{B}T\simeq 4.11\times 10^{-21}\mathrm{J}$. If the stated modulus is $G=2\times 10^5\mathrm{Pa}$, thenThe molar value counts moles of active strands. If “modulus” instead means Young's modulus $E$, a nearly incompressible isotropic rubber has $E\simeq 3G$, and the inferred value isThus the robust order-of-magnitude estimate is $10^{25}{\mathrm{m}}^{-3}$. The cross-link density inferred from rubber elasticity is usually an effective strand density; the actual junction density is ${\nu }_j=2{\nu }_s/f$ for junction functionality $f$ in an ideal network without inactive loops or dangling ends. For four-functional junctions this would be half the strand density. The paper supplies neither modulus type nor junction functionality, so a unique chemical junction count cannot be extracted.

A Michelson interferometer turns displacement into a count of interference fringes: a moving mirror changes both the outward and return paths. This lets a wavelength act as a ruler whose divisions need no engraving. Assume normal incidence, a fixed reference arm, stable refractive index $n$ along the moving arm, and fixed vacuum wavelength ${\lambda }_0$. Write $\lambda ={\lambda }_0/n$ for the wavelength in that medium and take $\Delta L$ positive away from the beam splitter. The change in optical path difference, the accumulated phase difference, and the signed number of complete interference fringe cycles areThe factor two comes from propagation geometry: moving the mirror by $\Delta L$ adds $\Delta L$ before reflection and $\Delta L$ afterward. One bright-to-bright cycle means mirror travel $\lambda /2$; bright-to-dark means $\lambda /4$, half a cycle. All lengths must use compatible units, and $\Delta N$ is dimensionless. A constant reflection phase difference changes the initial fringe position but cancels from this change formula.

For an original hypothetical example, suppose $\lambda =600\mathrm{nm}$ in the actual surrounding medium and a continuous scan gives $40.25$ cycles, including the final fractional cycle. Then $\Delta L=(600\mathrm{nm})(40.25)/2=12075\mathrm{nm}=12.075\mu \mathrm{m}$. With the wavelength treated as exact, a fringe-reading measurement uncertainty of $0.02$ cycle contributes $6\mathrm{nm}$ to the displacement measurement uncertainty. An uncertainty of $0.01\mathrm{rad}$ in a local phase difference instead corresponds to $600(0.01)/(4\pi )\mathrm{nm}\approx 0.477\mathrm{nm}$. This illustrates measurement resolution, not a complete claim of accuracy: wavelength calibration, air changes, mirror tilt and missed cycles can dominate.

To separate the oscillating readout from its optical coherence envelope, use the local modelHere $I$ is detected light intensity, $\bar{I}$ is the mean level, ${\varphi }_b$ is a fixed bias wave phase, and $\mathcal{V}$ is interferometric visibility. A fine scan of $D$ measures the carrier phase difference; a coarse scan can map the temporal optical coherence envelope through $\mathcal{V}(D)$. The envelope maximum locates matched delays only under the source and apparatus assumptions; it is not a bright fringe identifying an arbitrary mirror position. Even with known $\mathcal{V}$, one intensity reading leaves the cosine sign and whole-cycle ambiguities. Continuous motion of known direction, or an independently obtained pair of signals in phase quadrature, is needed to recover signed travel through reversals.

The figure's original control example stipulates equal arm light intensities, matching linear polarizations, and a Lorentzian optical spectrum giving $\mathcal{V}(D)=\exp (-|D|/L_c)$ with $L_c=1\mathrm{mm}$. In vacuum, arm mismatches $|L_2-L_1|=0.60\mathrm{mm}$ and $0.10\mathrm{mm}$ give $\mathcal{V}=e^{-1.2}\approx 0.301$ and $e^{-0.2}\approx 0.819$. Matching the arms improves the phase difference readout without changing its ${\lambda }_0/2$ travel per cycle. A fine mirror scan of one cycle changes $D$ by only $600\mathrm{nm}$ for the stipulated $600\mathrm{nm}$ source, much less than $L_c$. These are calculated values, not historical data.

A proposed measurement sequence follows from this model. Overlap the returning beams, match their linear polarizations, find a high-contrast delay and use a small calibrated mirror scan to determine the sign and fringe period. Read a fixed spatial sample of the interference fringes with a photodiode; record each arm separately and the dark signal, avoid detector saturation, and retain enough time resolution to count every cycle. Approach mechanical settings consistently and repeat the return scan to expose lost motion or drift. For very small changes, bias near phase quadrature, where $|dI/d\varphi |=\bar{I}\mathcal{V}$ is largest; an intensity extremum has zero first-order sensitivity. A known drive voltage alone does not provide measurement calibration for mirror travel: compare the actual travel with the fringe count using an independently known wavelength.

Control temperature, pressure, vibration and air currents because the instrument measures optical path difference, not mirror motion alone. For a fixed reference arm and fixed ${\lambda }_0$, the first-order change of the moving arm's round-trip path is $\Delta D\approx 2n\Delta L+2L\Delta n$. Thus an unrecognized $\Delta n$ mimics $\Delta L_{\mathrm{app}}=L\Delta n/n$. A hypothetical $L=0.10\mathrm{m}$, $n\approx 1$, and $\Delta n=10^{-7}$ give $10\mathrm{nm}$ apparent travel. Holding the mirror still while varying or monitoring the air conditions tests this confounder. Separate checks of alignment, wavelength stability and thermal expansion are necessary before sub-fringe measurement resolution becomes a defensible measurement uncertainty.

A $3\times 1.5\times 1.5\mathrm{m}{\mathrm{m}}^3$ single crystal of gadolinium antimonide carried a $2\mathrm{mm}$ fibre Bragg grating cyanoacrylate-bonded along its $a$ axis. The pulsed magnetic field paralleled that axis. Broadband light, a circulator and a spectrometer resolved reflected-wavelength shifts at $5.4\mathrm{K}$. Pre-pulse readings measured noise; low fibre thermal expansion reduced thermal background. A residual post-pulse shift indicated possible heating. Spectral distortion checked bond uniformity, while the fibre's smaller cross-section limited loading; bonding was not assumed force-free. Longitudinal strain showed a kink near $32.6\mathrm{T}$.^[ref]

The distinction between readout and inference is essential. The fibre Bragg grating senses a change of reflected wavelength; interpreting this as the specimen's strain requires good strain transfer and a controlled temperature. A useful independent control would compare a bonded sensor with a nearby unloaded sensor through the same pulse. Another would vary bond stiffness or specimen cross-section to test loading. These are proposed checks, not additional procedures claimed for the experiment above. The pulsed magnet's winding details were not specified in the inspected account.

This experiment measures one longitudinal component of strain. Establishing the first-order volume change would require three orthogonal normal components under the same conditions. GdSb also has its own antiferromagnetism and field-dependent magnetic structure: this measured curve does not establish the simple saturated isotropic law. Its role here is to make the measurement chain and its limitations concrete.