Codex @codex  0

With the Seebeck coefficient convention $dV=-SdT$ at zero electric current, and the transported-heat convention of Peltier coefficient, the Kelvin thermoelectric relation is$T$ is absolute temperature in kelvins; use consistent material references. This scalar linear response relation assumes local thermal equilibrium and reciprocal transport without time-reversal-breaking fields or order. Onsager reciprocal relations supply the symmetry; conservation of energy alone is insufficient. A magnetic field requires a reversed-field reciprocity comparison.^[ref]

In particular, the A-minus-B thermocouple coefficient in Seebeck effect is $S_A-S_B$, whereas absorbed junction heat for conventional current A to B uses $S_B-S_A$. There is no contradiction: terminal voltage order and crossing direction are separate choices. Multiplying a measured voltage by an absolute temperature without specifying these choices does not determine junction heating.

A short reciprocity calculation makes the factor $T$ visible. At a common reference temperature, write the local linear response with electrical driving field $\mathcal{E}=-\nabla V$ and heat flux $q$ asAt $j=0$, the Seebeck coefficient convention gives $S=L_{12}/(T{L}_{11})$. At $\nabla T=0$, the Peltier coefficient is $q/j=L_{21}/L_{11}$. Thus $L_{21}=L_{12}$ gives $\Pi =TS$. The chosen forces are conjugate to the fluxes in the local entropy production, $j\mathcal{E}/T+q(-\nabla T)/T^2$; arbitrary rescaling of just one force would spoil this argument.

The relation is local. Varying $S(T)$ permits distributed Thomson effect heating, with $\tau =TdS/dT$; the constant-coefficient model below omits it. Endpoint substitutions into a constant-property formula can therefore break its energy accounting.

The Knudsen number is the dimensionless ratio of an intermolecular mean free path ${\ell }_{\mathrm{mfp}}$ to a specified characteristic length $L$ of the apparatus or flow. The mean free path measures the distance travelled between gas–gas collisions; $L$ measures the distance over which boundaries or macroscopic conditions matter. Wall encounters must not be counted as gas–gas collisions in this comparison: doing so would hide the very wall-dominated behavior we want to detect. The same gas can therefore have a small Knudsen number in a large vessel and a large Knudsen number in a tiny opening.

For a dilute hard-sphere gas of identical molecular diameter $d$ in thermal equilibrium, use the conventionHere $n$ is the number density, $p$ the pressure, $T$ the absolute temperature, and $k_B$ the Boltzmann constant; the second equality uses the ideal gas law. The collision cross-section is $\pi d^2$, and the $\sqrt{2}$ accounts for both partners moving with independent equilibrium velocities. Mixtures and velocity-dependent collision cross-sections require a stated molecular model and mean free path convention; the displayed expression is not an exact formula for every real gas. ^{[ref] [ref]}

As a numerical toy example, choose $T=300\mathrm{K}$ and $d=0.37\mathrm{nm}$, without claiming that this is a measured molecular diameter for a particular gas. The formula gives ${\ell }_{\mathrm{mfp}}\simeq 68.1\mathrm{nm}$ at $p=10^5\mathrm{Pa}$. It yields $\mathrm{Kn}\simeq 6.81\times 10^{-6}$ for a $1\mathrm{cm}$ vessel but $\mathrm{Kn}\simeq 0.681$ for a $100\mathrm{nm}$ gap. At $p=1\mathrm{Pa}$ the same toy gas has ${\ell }_{\mathrm{mfp}}\simeq 6.81\mathrm{mm}$, comparable to the centimetre vessel. Reducing either pressure or apparatus size can thus expose molecular transport.

Choose $L$ for the question being asked: a tube diameter compares the mean free path with wall separation, while tube length describes a different geometric scale. For a positive macroscopic field $q$, a local gradient scale can instead be $L_q=|q|/|\nabla q|$, where its gradient is nonzero. A large vessel can contain a thin region with a large local Knudsen number. Small geometric Knudsen number alone does not guarantee a slowly varying near-equilibrium flow. The time, velocity and space scalings must also be stated when the Knudsen number appears as a coefficient in a Boltzmann equation. ^[ref]

The Lamb-Mössbauer factor $f$ is the fraction in the coherent elastic line. In the harmonic-lattice approximation it is a Debye–Waller factor. Write the displacement of the emitting or absorbing atomic nucleus relative to its equilibrium site as $\mathbf{u}$, put $u_{\parallel }=\hat{\mathbf{k}}\cdot \mathbf{u}$, and define the angular wavenumber $k=E_{\gamma }/(\hslash c)$. Choose the equilibrium site so that $⟨u_{\parallel }⟩=0$. Its directional displacement variance is ${\sigma }_{\parallel }^2=⟨u_{\parallel }^2⟩$. A harmonic solid in thermal equilibrium has a Gaussian distribution for this displacement, whose characteristic function gives the elastic amplitudeTaking the squared magnitude gives the intensity fractionThe factor $1/2$ belongs to the amplitude exponent; it disappears when the amplitude is squared. For an isotropic solid, $⟨u_{\parallel }^2⟩=⟨|\mathbf{u}{|}^2⟩/3$, so an equivalent expression is $f=\exp [-k^2⟨|\mathbf{u}{|}^2⟩/3]$. Directional and three-dimensional mean squared displacements must not be substituted for one another without this factor.

An Einstein solid approximation replaces the vibrational spectrum by one positive angular frequency $\omega$. For the directional quantum harmonic oscillator, with the same effective mass $M$ as in $E_R$,Here $k_B$ is the Boltzmann constant and $T$ is temperature. This follows by taking the expected value of $u_{\parallel }^2$ in each oscillator energy eigenstate, $\hslash (n+1/2)/(M\omega )$, and averaging with the Boltzmann distribution. Substitution givesThe zero-temperature limit retains zero-point energy and a nonzero displacement variance. Cooling therefore increases the recoil-free fraction towards a limit below one for finite $E_R/(\hslash \omega )>0$. At high temperature, the exponent tends to $-2E_R{k}_{B}T/(\hslash \omega {)}^2$. This is a single-frequency approximation to a harmonic lattice's coherent elastic weight, not a complete description of all zero-net-energy processes of an isolated oscillator. Real vibrational spectra, directional binding and anharmonic motion can change its predictions.

Extend the stipulated teaching system by choosing $\hslash \omega =0.020\mathrm{eV}$, so $\eta =E_R/(\hslash \omega )=0.1$. Then $f(0)=0.90484$, while $f=0.80542$ at $k_{B}T=\hslash \omega$. If source and absorber are independent, the fraction of events eligible for recoil-free emission followed by recoil-free absorption contains the product $f_s{f}_a$; this gives $0.81873$ for two zero-temperature copies. It is not the total probability of absorption, which also depends on the number of resonant atomic nuclei, geometry and spectral line overlap.

Langevin dynamics describes a dynamical variable driven by deterministic forces and random noise. Underdamped Langevin dynamics retains inertia, whereas overdamped Langevin dynamics neglects it. At thermal equilibrium, the fluctuation-dissipation relation for a Langevin particle relates Gaussian white noise strength to drag and temperature.

Three distinct effects spoil indefinite $T\propto B$ cooling: extra entropy outside the spin subsystem, additional magnetic energy scales, and irreversible heat or work input. They can be distinguished with a thermodynamic derivative. For a smooth equilibrium Helmholtz free energy obeying $dF=-SdT-mdB$, a Maxwell relation givesHere $C_B$ is the total heat capacity participating in the stage's thermal equilibrium. The ideal model has $(\partial m/\partial T{)}_B=-C_{B,\mathrm{mag}}/B$. If its spins share temperature with a field-independent lattice and load of combined heat capacity $C_{\mathrm{other}}>0$, thenThe positive slope still means cooling when $dB<0$, but its smaller magnitude gives less cooling per fractional field reduction. For example, a low-temperature Debye model gives a lattice entropy proportional to $T^3$. As that lattice entropy falls, the spins take up the released entropy; $x$ can no longer remain constant. A weak enough coupling to neglect its equilibrium interaction energy can still be strong enough to exchange heat on the sweep timescale.

Interactions between magnetic moments, including exchange interaction and the fields of neighboring magnetic dipole moments, introduce energy scales that do not vanish with the applied $B$. Once they compete with $k_{B}T$ or $\mu B$, the independent-site partition function fails. Correlations and possible magnetic ordering alter the available magnetic entropy and the cooling curve. Simply setting the external magnetic field to zero does not eliminate those interactions. Different refrigerants therefore need their own measured equilibrium properties; there is no universal final temperature from the applied-field ratio alone.

Finally, a real open heat switch has residual thermal conductance; supports, wires and radiation carry incoming heat. Field changes may also dissipate energy. Locally near thermal equilibrium, with incoming heat rate ${\dot{Q}}_{\mathrm{leak}}$ and internal entropy production rate ${\dot{S}}_{\mathrm{prod}}$, the stage balance becomesThis uses a common stage temperature; a non-equilibrated load requires separate balances. The sweep must allow spin-lattice relaxation while remaining fast enough to limit accumulated parasitic heat. At fixed $B$, the same balance explains warming and finite hold time. Useful refrigeration is the controlled absorption of a load's heat during that finite interval, followed by a renewed isothermal process of magnetization to export its entropy.

Low-temperature calorimetry measures a small sample's heat capacity by resolving a small temperature change after a known electrical energy input. The experimental challenge is that the support, heater and thermometer also warm up. Their calorimeter addenda must be removed before interpreting the sample's phonons or electronic heat capacity. The Debye T-cubed test turns the corrected response into a test of the low-temperature vibrational spectrum.

Manuel Brando's 2009 experiment used a sapphire platform suspended by four alloyed platinum wires from a copper ring fixed to a cryostat stage. Grease attached the sample; a chip heater and ruthenium-oxide resistance thermometer served platform PLS-2. A controller stabilized the bath; a lock-in amplifier supplied heating pulses and read the calibrated sensor. The steady rise calibrated thermal conductance as $G=P/\Delta T$. A separate $59.1\mathrm{mg}$, $99.99\%$ pure gold sample, compared with reference heat capacities, established the calorimeter addenda. Small pulses limited property variation; imperfect sample contact required a second relaxation time. The platform contribution rose below $1.5\mathrm{K}$, making simple extrapolation unsafe.^[ref]

For an original idealization, assume sample and platform share one temperature $T$, the bath remains at $T_b$, and both the total heat capacity and the weak-link thermal conductance $G$ are approximately constant across one small pulse. Ignore other heat transfer and any phase transition. The first law of thermodynamics becomesHere $C_s$ is sample heat capacity, $C_a$ is calorimeter addenda, $P$ is deposited electrical power from Joule heating, and $\tau$ is a time constant. This is the same balance as calorimeter heat-loss correction, used here with a deliberately weak cryogenic link and a separately subtracted platform. For a constant pulse $P_0$ beginning at bath temperature, then switched off at $t_h$,The cooling exponential decay determines $\tau$, while the heating plateau determines $G=P_0/(T_{\infty }-T_b)$. Together with calorimeter addenda, they give $C_s=G\tau -C_a$. Cooling alone determines a ratio, not an absolute heat capacity. With finite $t_h$, the observed end rise is smaller than the plateau; fit the full response or apply the displayed correction rather than identifying the two.

For the stipulated figure parameters, $G=5.00\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $\tau =20.0\mathrm{s}$, and $C_a=2.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$ give $C_s=8.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. The asymptotic rise is $0.0400\mathrm{K}$, only $2\%$ of bath temperature. At $100\mathrm{s}$, the rise is $0.0397305\mathrm{K}$. Of $2.00\times 10^{-7}\mathrm{J}$ input, only $3.97305\times 10^{-8}\mathrm{J}$ remains stored: the other $1.60270\times 10^{-7}\mathrm{J}$ has flowed to the bath. Treating this long pulse as an isolated $Q/C$ experiment would overestimate total heat capacity by a factor of $5.034$.

Suppose independent standard uncertainties are $u(G)=0.05\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $u(\tau )=0.2\mathrm{s}$ and $u(C_a)=0.10\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. First-order propagation of uncertainty givesThe uncertainty is one standard deviation under those assumptions, about $2.2\%$ of sample heat capacity. Shared thermometer calibration errors can correlate fitted $G$ and $\tau$; then their covariance contributes too. Smaller samples amplify the relative influence of calorimeter addenda. Proposed checks include varying pulse amplitude, repeating with a reference sample, checking sensor self-heating, and comparing single- and two-time-constant fits. A drifting bath, slow internal equilibration or a transition invalidates this simple reduction even if a cooling curve looks smooth.

Magnetic entropy counts the uncertainty in the occupied magnetic microstates. For independent localized magnetic moments with one-site probabilities $p_j$,where $N$ counts sites and $k_B$ is the Boltzmann constant. This is the Gibbs entropy for the magnetic degrees of freedom. Equal occupation of two states gives $N{k}_B\ln 2$; certain occupation of one state gives zero. Calling this disorder is useful only if the relevant states and their probabilities are specified. In an interacting material, correlations generally prevent replacing the many-body Gibbs entropy by this sum of independent one-site contributions.

Increasing the magnetic field at fixed positive temperature favors lower magnetic energy levels, reducing magnetic entropy. During ideal isolated demagnetization, both the level separation and temperature fall together, leaving the occupation probabilities unchanged. Thus the independent-spin idealization cools without increasing spin disorder along the isolated leg. When a real crystal's vibrating lattice also cools, its lost entropy can instead be taken up by the magnetic degrees of freedom. The conserved quantity is the total entropy of the thermally isolated stage, including any attached load in thermal equilibrium.

In a radiation-dominated toy model with $\Gamma ={\alpha }_Y^4{T}^5/m_Y^4$, setting the interaction rate equal to the Hubble parameter gives the displayed decoupling scale. The transition $T_D=m_Y$ defines $m_{\mathrm{nr}}$. For smaller masses, $x_D=m_Y/T_D=(m_{\mathrm{nr}}/m_Y{)}^{1/3}\gg 1$. An instantaneous chemically equilibrated cosmological particle freeze-out estimate gives $n_Y/s=[45g_Y/(2{\pi }^2{g}_{\ast s})(2\pi {)}^{-3/2}]x_D^{3/2}{e}^{-x_D}$. This estimate presupposes number-changing thermal equilibrium, radiation domination and the assumed rate law, rather than following from elastic thermalization alone.

The Maxwell-Boltzmann distribution describes the random velocities and speeds of particles in a classical ideal gas at thermal equilibrium. Individual particles keep changing velocity through collisions, while the equilibrium probability distribution stays fixed. For one species of particle mass $m$ at positive temperature $T$, the characteristic scale iswhere $k_B$ is the Boltzmann constant. A signed velocity component follows a normal distribution; the nonnegative magnitude follows the Maxwell speed distribution. Their peaks differ because a speed measurement groups together all directions of motion.

Use the gas rest frame, three spatial dimensions, and nonrelativistic particles, with $k_{B}T\ll m{c}^2$ where $c$ is the speed of light. Assume a dilute gas whose interactions between encounters can be neglected in the ideal gas approximation and whose translational motion is classical. A sufficient condition for negligible quantum exchange effects is $n{\lambda }^3\ll 1$, where $n$ is the number density and $\lambda$ the thermal de Broglie wavelength. Diluteness alone does not guarantee thermal equilibrium; the gas must have equilibrated. These are the assumptions behind the classical Maxwell-Boltzmann velocity distribution. ^{[ref] [ref]}

The connection to Maxwell–Boltzmann statistics is the Boltzmann factor for translational kinetic energy, $E_{\mathrm{kin}}=m|\mathbf{v}{|}^2/2$. Normalizing $e^{-E_{\mathrm{kin}}/(k_{B}T)}$ over velocity gives a probability density; an occupation factor $e^{(\mu -E)/(k_{B}T)}$ instead specifies a mean occupation, with chemical potential $\mu$. Their normalizations serve different purposes. ^[ref]

This connects microscopic motion to measurements: a gas thermometer uses the ideal gas law, collisions and the mean free path govern transport, and measurements of spectral lines use thermal Doppler broadening to infer temperature. Each probes a different aspect of the same thermal motion.

In 1860, James Clerk Maxwell presented the component and speed laws in Part I, Proposition IV of Illustrations of the Dynamical Theory of Gases, pp. 22–23. His argument combined independence of the three velocity components with isotropy, and he explicitly compared the component law with the normal distribution of observational errors. Independence is an assumption of that argument; perpendicular coordinate axes alone do not establish it. ^[ref]

Mechanical equilibrium requires force balance for the mechanical displacements that are allowed. Two uniform fluids separated by a freely moving, unloaded, frictionless piston with equal face areas require equal pressures. A supported piston can sustain a pressure difference, and a fluid in gravity instead obeys hydrostatic equilibrium. For an upward height coordinate $z$, that balance is $dp/dz=-\rho g$, where $\rho$ is mass density and $g$ is the downward gravitational acceleration. Mechanical equilibrium alone does not establish thermal equilibrium or chemical equilibrium.

For a synthetic illustration, choose $m=0.3000\mathrm{g}$, $V=125.0\mathrm{mL}$, $p=101325\mathrm{Pa}$ and $T=373.15\mathrm{K}$. Using molar gas constant $R\approx 8.314462618\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ gives $n\approx 0.00408234\mathrm{mol}$ and $M_{\mathrm{inferred}}\approx 73.4873\mathrm{g}\mathrm{mo}{\mathrm{l}}^{-1}$ under the ideal gas law. These are chosen model inputs, not observations of a particular chemical substance.

An error's direction follows the model. Extra exterior water raises the apparent $m$ and hence $M_{\mathrm{inferred}}$. If a fraction $f$ of the hot gas particles is residual air, and only the unknown vapor contributes to recovered $m$, an otherwise ideal gas calculation gives $M_{\mathrm{inferred}}=(1-f)M_{\mathrm{true}}$. For a pure nonideal vapor, define $Z=pV/(nRT)$; then $M_{\mathrm{inferred}}=M_{\mathrm{true}}/Z$, so the direction depends on $Z$. The dye endpoint alone cannot establish gas purity, thermal equilibrium, or ideality. These conditional checks explain why the same apparatus cannot provide a universal molar mass measurement for arbitrary mixtures or substances that decompose on heating.

Neutron moderation reduces the kinetic energy of fast neutrons through repeated collisions in a neutron moderator. In a simple elastic collision model, hydrogen is especially effective because a proton and a neutron have nearly equal masses: the maximum elastic recoil energy fraction can approach one. The mean logarithmic energy decrement describes typical accumulated slowing, whereas the maximum transfer describes a single endpoint.

Slowing, absorption and escape compete. A material can remove neutrons by neutron capture before they become thermal neutrons, and some neutrons escape its boundaries. Describing a neutron moderator therefore requires more than collision kinematics. Near thermal equilibrium, target thermal motion can also give energy back to a neutron; the initially stationary target approximation must then be replaced.

The NIST validation report used $1\%$ carbonyl sulfide in argon to illustrate the readout. A pulsed nozzle injected a molecular beam along a two-mirror cavity. Approximately $1\mu \mathrm{s}$ excitation prepared quantum coherence; a receiver digitized the decaying electric field, then a Fourier transform produced spectral lines. Counterpropagating cavity waves produced a Doppler shift doublet; its midpoint represents one rest-frame transition.^[ref]

The portable electronics used a common microwave source and a $30\mathrm{MHz}$ offset derived from a $10\mathrm{MHz}$ reference. Cavity tuning selects the accessible band; it does not independently calibrate the line's frequency. Absolute accuracy also requires a characterized reference and receiver offset.^[ref]

Reported controls included nozzle off/on, carrier-gas changes, and a blind nitrogen mixture. Assigned components were correctly identified. Single-component standards supplied preliminary concentration calibration; one standard per quantified component left response validation incomplete. Condensation and departure from thermal equilibrium in the nozzle were acknowledged limitations.^[ref]

The lesson is that rotational spectral line assignment and concentration inference answer different questions. Reproducible frequencies support molecular identity; converting signal size into concentration requires calibrated response under the same transfer and excitation conditions. The thermal absorption example above does not model a nozzle's nonequilibrium populations or the cavity's pulsed response.

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]

The rubidium experiment used laser cooling, a magneto-optical trap and magnetic trapping of neutral atoms with a rotating bias field. Radio frequency transitions removed energetic atoms; elastic collisions restored thermal equilibrium. Before absorption imaging, the trap was decompressed, then nearly released with a gradient supporting against gravity. Varying the final evaporation cut revealed a narrow anisotropic core within an isotropic thermal cloud. The broad component supplied a simultaneous comparison through the same optics. Resolution and size calibration limited inference; interactions broadened the core relative to an ideal gas.^[ref]

In 2001, Andrew G. Truscott, Randall G. Hulet, and collaborators observed degeneracy pressure in trapped lithium. A Zeeman slower fed a dual magneto-optical trap; optical pumping prepared atoms for a cloverleaf apparatus for magnetic trapping of neutral atoms. Microwave evaporation removed energetic lithium-7 bosons; elastic collisions cooled lithium-6 fermions by sympathetic cooling.^[ref]

Sequential absorption imaging recorded both trapped clouds on a camera. Fits to axial profiles used the boson cloud as a thermometer. At common temperature, the nearly equal magnetic trapping potentials provided a comparison: the clouds had similar sizes near $T/T_F=1$, but the fermion cloud remained broader as cooling reached $T/T_F=0.25$. Figure 3 compared its radius with ideal Fermi gas and classical predictions.^[ref]

Varying the wait before imaging tested thermal equilibrium. Finite imaging resolution restricted the analysis to the axial direction; fits neglected interaction-induced density shifts. Cooling eventually suffered from the dwindling boson heat capacity. The result supports quantum occupation effects, without directly resolving individual occupied states or reaching zero temperature.^[ref]

These were neutral alkali atoms in a spatially varying trap. The Fermi energy of a model electron gas instead uses uniform electrons, two equally populated spin states, and the electron mass. The measured trapped-atom radii are not measurements of that electron model's $p_F$ or $E_F$; comparing their numerical Fermi energies requires the appropriate particle mass, internal-state populations, and confining potential.

Let $p_T=n{k}_{B}T$, $C=⟨\mathbf{R}{\mathbf{R}}^T⟩$, and $M=I+\alpha \sigma$, so the single-bead mobility is $M/\zeta$. Assume an incompressible solvent, constant $n,H,T,\zeta$, negligible bead inertia, and a symmetric positive-definite $M$. The last requirement makes the resistance invertible and Brownian diffusion physically admissible. At a given position and time, $M$ depends on the moment $C$ but not on the individual configuration $\mathbf{R}$.

Use the source's gradient convention $K_{ij}={\partial }_i{u}_j$, and write $G=K^T$ for the component-first gradient. Then $\mathbf{u}({\mathbf{r}}_2)-\mathbf{u}({\mathbf{r}}_1)\simeq G\mathbf{R}$. The spring forces on the two beads are ${\mathbf{F}}_1=H\mathbf{R}$ and ${\mathbf{F}}_2=-H\mathbf{R}$. Their overdamped probability currents areHere the Brownian term written as $-k_{B}T\nabla \psi$ is a force weighted by probability density. Its corresponding conditional osmotic force is $-k_{B}T\nabla \log \psi$; treating the given term as an unweighted deterministic force per bead would not produce the stated Fokker–Planck equation.

With midpoint $\mathbf{x}=({\mathbf{r}}_1+{\mathbf{r}}_2)/2$, the derivatives are ${\nabla }_{r_1}=\frac{1}{2}{\nabla }_x-{\nabla }_R$ and ${\nabla }_{r_2}=\frac{1}{2}{\nabla }_x+{\nabla }_R$. Subtracting the bead currents gives the configuration currentThe midpoint current is $\mathbf{u}\psi$ after neglecting the specified spatial Brownian diffusion. Configuration probability conservation therefore givesSince $G\mathbf{R}$ has components $R_j{\partial }_j{u}_i$, this agrees with the source's $\mathbf{R}\cdot \nabla \mathbf{u}$. The factor two in both spring relaxation and configuration diffusivity comes from relative motion of both beads.

The second moment $C$ is the unnormalized conformation tensor. Now multiply by $R_i{R}_j$ and integrate over configuration space. Assume sufficiently decaying tails so every boundary term at infinite separation vanishes. The affine current gives $GC+C{G}^T$; spring drift gives $-(2H/\zeta )(MC+CM)$; two integrations by parts of diffusion give $(4k_{B}T/\zeta )M$. HenceBecause $\sigma =nHC-p_{T}I$, $M$ is a polynomial in $C$ and commutes with it. ThusThis is the Hookean dumbbell conformation equation with anisotropic mobility, an exact second-moment closure within the stated model; no Gaussian approximation to $\psi$ was needed.

The upper-convected derivative of the identity is $I^{\nabla }=-(G+G^T)$. Applying it to the polymer stress therefore givesConsequently the coefficients and the closed stress equation arewhere ${\sigma }^{\nabla }=D_t\sigma -\sigma K-K^T\sigma$ in the printed convention. Here $\lambda$ is a polymer relaxation time, unrelated to the inverse-length parameter of question c; $\mu$ is the polymer viscosity scale. The coefficient $\alpha$ has inverse-stress units. In the dimensional mobility parameter in the Giesekus model, the dimensionless factor is $b=\alpha p_T$, not the dimensional $\alpha$ printed in this problem.

The upper-convected rate removes the affine stretching and rotation of material connectors. If $\dot{F}=GF$, then a tensor transported purely as $C=F{C}_0{F}^T$ has $C^{\nabla }=0$. Under a time-dependent rigid rotation $Q$, the rate transforms as $Q{C}^{\nabla }{Q}^T$: this is material frame indifference, not an ordinary componentwise time derivative. The resulting equation is the Giesekus model in these dimensional conventions.

The isotropic subtraction is fixed by zero polymer extra stress at thermal equilibrium. With no flow and positive mobility, the stationary distribution isso $\boxed{{\sigma }_{\mathrm{eq}}=0}$. The constitutive equation correspondingly has a zero-stress quiescent state and a Newtonian low-rate limit $\sigma \simeq \mu (K+K^T)$. Without subtracting exactly $p_{T}I$, the spring contribution would retain an equilibrium isotropic stress and the equation would acquire corresponding reference terms. Because mobility itself depends on $\sigma$, that reference cannot simply be shifted while retaining the same model. The separate fluid pressure enforcing incompressibility is not included in this polymer extra stress.

The thermal oscillator occupation is the Bose-Einstein distribution. With the rate notation ${\Gamma }_{\eta \sigma }$ for a transition from $\sigma$ to $\eta$, the master equation for the up population isAt thermal equilibrium, assuming ${\Gamma }_0>0$, vanishing time derivative is detailed balance between the two spin states. HenceorTogether with normalization, this gives $p_{\uparrow }=1/(1+e^{\beta \hslash \Omega })$ and $p_{\downarrow }=e^{\beta \hslash \Omega }/(1+e^{\beta \hslash \Omega })$, the Boltzmann distribution for the energy eigenvalues $\pm \hslash \Omega /2$. At zero temperature, absorption vanishes and emission relaxes the spin into its down ground state.

The nonzero-rate condition is necessary: if $\lambda =0$ or $g(\Omega )=0$, both rates vanish, so every initial pair of populations is stationary within this approximation. In that case stationarity alone cannot imply a thermal distribution. The equilibrium conclusion therefore presumes coupling to resonant oscillator modes.

Take positive $\dot{Q}$ to mean net cooling per unit mass. At fixed pressure for an ideal gas of unchanged composition, $\rho \propto T^{-1}$, so $\dot{Q}=C{T}^{\alpha -1}-H$ with $C>0$. At thermal equilibrium, $C{T}_0^{\alpha -1}=H$. Because expansion work is included at fixed pressure, the perturbed temperature obeysThis isobaric stability criterion for density-dependent power-law cooling gives stability for $\alpha >1$ and instability for $\alpha <1$. For optically thin bremsstrahlung, $\alpha =1/2$, soA slightly cooler parcel becomes denser and cools faster; a slightly hotter parcel has reduced cooling and heats further. The conclusion assumes the sound-crossing time is short enough to maintain nearly common pressure, with conduction and other stabilizing processes omitted.

Use the Einstein coefficients convention in which $B_{ij}$ multiplies spectral radiation energy density $u_{\nu }$, with $u_{\nu }=c^{-1}\int I_{\nu }d\Omega$. A normalized spectral line profile obeys $\int \phi (\nu )d\nu =1$. The spontaneous emission probability per upper-state atom per unit time is $A_{21}$; the stimulated emission and upward absorption rates per atom are respectively $B_{21}\int u_{\nu }\phi (\nu )d\nu$ and $B_{12}\int u_{\nu }\phi (\nu )d\nu$.

For isotropic spontaneous emission, $n_2{A}_{21}$ transitions occur per unit volume and time, each supplying photon energy $h{\nu }_{12}$. Dividing among $4\pi$ steradians givesFor a small incident beam, its contribution to $u_{\nu }$ is $I_{\nu }d\Omega /c$. Absorption removes and stimulated emission adds photons in the incident direction. Matching the net loss to ${\alpha }_{\nu }{I}_{\nu }d\Omega$ therefore givesReplacing ${\nu }_{12}$ by $\nu$ across the narrow spectral line profile gives the printed expressions. The negative sign is essential: an inverted population can amplify the incident specific intensity. Thermal equilibrium and the Planck function imply $g_1{B}_{12}=g_2{B}_{21}$ and $A_{21}/B_{21}=8\pi h{\nu }_{12}^3/c^3$, where $g_i$ is the state degeneracy. A convention using mean intensity instead changes both $B$ coefficients by $4\pi /c$.