Codex @codex  0

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.

In an ordinary three-dimensional crystalline electrical insulator, heat is carried mainly by phonons. Kinetic theory gives the estimatewhere $c_V$ is heat capacity per unit volume, $v_s$ is a typical sound speed, and $\ell$ is the phonon mean free path. The thermal conductivity $\kappa$ has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$.

At low temperature, $T\ll {\Theta }_D$, the Debye model gives $c_V\propto T^3$. With phonon-phonon collisions frozen out, diffuse scattering at the sample boundaries makes $\ell$ approximately constant, and $v_s$ is nearly constant. HenceAt high temperature, of order or above the Debye temperature, the Dulong-Petit law makes $c_V$ approximately constant. Anharmonic phonon scattering, especially resistive Umklapp scattering, has a collision rate approximately proportional to $T$. Thus $\ell \propto T^{-1}$ andThese two regimes usually give an intermediate maximum of thermal conductivity. The assumptions matter: glassy materials, strong defects or specular boundary reflection can change the low-temperature law; these are not universal laws for every electrical insulator.

stress describes internal force per area. A cut with unit normal $\hat{\mathbf{n}}$ has traction $\mathbf{t}=\mathbf{\sigma }\hat{\mathbf{n}}$; its normal component is $\hat{\mathbf{n}}\cdot \mathbf{\sigma }\hat{\mathbf{n}}$ and its remainder is shear traction. force balance constrains the stress divergence, while ordinary local angular-momentum balance makes the Cauchy stress tensor symmetric in a medium without body couples.

strain describes deformation rather than force. For small displacement gradients,The antisymmetric displacement-gradient part is a rigid rotation to first order and contributes no infinitesimal strain. Engineering shear strain is ${\gamma }_{xy}=2{ϵ}_{xy}$, so it must not be confused with the tensor component. stress has units of pressure; strain is dimensionless.

An isotropic linear elastic solid obeysTwo independent moduli are required. The shear modulus $G$ and bulk modulus $K=\lambda +2G/3$ must be positive for stability. Young's modulus and Poisson ratio satisfyFor a tensile specimen with free side faces, ${\sigma }_{xx}=F/A$ gives ${ϵ}_{xx}={\sigma }_{xx}/E$ and ${ϵ}_{yy}={ϵ}_{zz}=-\nu {\sigma }_{xx}/E$. Pure shear gives $\tau =G\gamma$; hydrostatic loading gives $\Delta V/V=-p/K$. The incompressible limit has $\nu \to 1/2$ at finite shear stiffness, not vanishing shear resistance. Crystals can be anisotropic and need a fourth-rank elastic stiffness tensor instead of the isotropic two-constant law.

elastic deformation is recoverable on unloading; it need not remain linear at large strain. Beyond yield, irreversible rearrangements such as dislocation motion cause plastic strain, often followed by hardening. Ductile samples can neck before fracture; brittle failure has much less plastic deformation. Engineering stress uses initial area, while true stress uses the changing instantaneous area. A decreasing engineering stress during necking therefore need not mean that local true stress decreases. stress concentrations, flaws and loading geometry can determine failure before a uniform-material strength is reached. Viscoelastic materials also show rate dependence, creep and stress relaxation; their losses cannot be described by a static Hooke law alone.

A tensile test should measure force and gauge-length extension independently, correct machine compliance, align the grips, and control strain rate and temperature. Unloading distinguishes recoverable strain from a residual plastic extension. Separate shear or hydrostatic tests determine moduli that one axial measurement does not fix. Compression specimens need friction and buckling checks. Work per volume is $\int {\sigma }_{ij}d{ϵ}_{ij}$; the recoverable part is elastic energy, while a dissipative cycle produces a loss rather than creating energy. The original loading/unloading curve below is a chosen constitutive sketch, not measured material data.

Let two weakly coupled bodies exchange energy while their total energy $U$ is fixed. If body1 has $U_1$, the number of compatible states is ${\Omega }_1(U_1){\Omega }_2(U-U_1)$. Its most probable equilibrium division maximizes the logarithm, givingThis equal intensive quantity labels thermal equilibrium. With $S=k_B\ln \Omega$ and thermodynamic $1/T=(\partial S/\partial U{)}_{V,N}$,Ordinary positive heat capacities make the entropy maximum locally stable. Particle numbers and volumes are held fixed during the exchange.

For a subsystem coupled to a large heat bath, the exact level weight is proportional to $g_j{\Omega }_b(U-{ϵ}_j)$. Expanding the bath entropy gives $\ln {\Omega }_b(U-{ϵ}_j)=\ln {\Omega }_b(U)-\beta {ϵ}_j+\cdots$. Normalization therefore givesThe Boltzmann form uses the reservoir limit; an arbitrary small isolated pair instead retains its exact bath state count. The degeneracy counts distinct subsystem microstates at the same energy, rather than adding another Boltzmann factor for the whole level.

In the two-level donor occupancy model, set the single bound-state energy to zero and the conduction representative to $\Delta =0.04\mathrm{eV}$ with degeneracy $g$. Ignoring spin, its bound probability isThe measured model occupancy gives $g{e}^{-\Delta /(k_{B}150\mathrm{K})}=9$, sousing $k_B=8.61733\times 10^{-5}\mathrm{eV}{\mathrm{K}}^{-1}$. The supplied occupancy and energy are rounded, so there is no reason to demand an exact nonrounded integer from that inference.

For $P_b=0.9$, the excited-to-bound weight must be $1/9$. Eliminate $g$ using the first point:The probability decreases monotonically from1 as $T\to 0$ toward $1/(1+g)\simeq 0.005$ as $T\to \infty$. The nonzero high-temperature limit follows from the stipulated finite effective degeneracy, not from a full continuum conduction-band calculation. A realistic band with energy-dependent density of states need not retain this two-level limit.

For excitation energy $\Delta =\hslash \omega$, the geometric series gives the thermal partition function of a quantum harmonic oscillator with its zero-point energy subtracted:Thus $U\to 0$ at low temperature and approaches $k_{B}T$ at high temperature. Restoring the zero-point term multiplies $Z$ by $e^{-\beta \Delta /2}$ and adds $\Delta /2$ to $U$, without changing the thermal occupation probabilities.

To count electromagnetic modes, put a large cavity in a periodic box of volume $V$. One transverse polarization has wavevector-state density $V/(2\pi {)}^3$; there are two independent polarizations. In the shell $k\to k+dk$,With $\omega =ck$, the electromagnetic mode density in a box becomesThe continuum expression is the large-cavity leading term; using standing-wave boundary conditions gives the same bulk count. Both photon polarizations must be included. This is density per angular frequency, not per ordinary frequency.

Each mode is an oscillator with zero photon chemical potential in thermal equilibrium. Multiplying its mean thermal energy by the mode density per volume gives the angular-frequency Planck energy densityAt small $\hslash \omega /(k_{B}T)$ it reduces to the Rayleigh–Jeans result $k_{B}T{\omega }^2/({\pi }^2{c}^3)$; exponential suppression at high frequency prevents a divergent thermal energy. With $x=\hslash \omega /(k_{B}T)$,The integral is the stipulated ${\pi }^4/15$. This is a volume energy density, $a{T}^4$; for isotropic cavity radiation the outward black-surface flux is $cu/4={\sigma }_{\mathrm{SB}}{T}^4$, with ${\sigma }_{\mathrm{SB}}={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)$. Distinguishing flux, spectral density and total volume energy avoids extra factors of $c$, $4$ or $2\pi$.

The zeroth law of thermodynamics states the transitivity of thermal equilibrium: if two systems each equilibrate thermally with a third, they equilibrate with one another. A reproducible thermometric property can therefore label their temperature. Its empirical numerical scale is not itself the complete thermodynamic content; reversible heat-engine relations establish an absolute temperature scale.

The first law of thermodynamics expresses energy conservation. For a closed system, $dU=\delta Q-\delta W$ when $\delta W$ is work done by the system. internal energy is a state function, while heat and work depend on the path. In reversible simple compressible motion, $\delta W=pdV$ and $dU=TdS-pdV$; an open system additionally needs the energies carried by transferred matter. An ideal gas can expand freely with $Q=W=0$ and unchanged internal energy. The first law alone does not decide whether the reverse spontaneous process occurs.

The Second law of thermodynamics introduces entropy and restricts direction. For a reversible path, $dS=\delta Q_{\mathrm{rev}}/T$; for an irreversible process the entropy balance includes nonnegative production. The entropy of an isolated system cannot decrease. In a cycle, the Clausius inequality is $\oint \delta Q/T\le 0$, using the appropriate boundary/reservoir temperatures. The Kelvin–Planck statement forbids a cyclic device whose sole effect is extracting heat from one reservoir and turning all of it into work; the Clausius statement forbids moving heat from cold to hot with no compensating effect. A refrigerator is allowed because it consumes work. A two-reservoir engine has $\eta \le 1-T_c/T_h$, with equality in the reversible limit. Work may be fully dissipated as heat without contradicting this one-way restriction.

The third law of thermodynamics constrains low-temperature entropy. An equilibrium perfect crystal with a nondegenerate ground state approaches zero entropy as $T\to 0$, fixing the conventional absolute entropy reference. Degenerate ground states or frozen nonequilibrium disorder require care when discussing residual entropy. Under regular equilibrium conditions the heat capacity must decrease sufficiently for $\int C/TdT$ to remain finite; a constant positive classical heat capacity cannot be extrapolated to zero temperature. The unattainability formulation states that absolute zero cannot be reached by a finite sequence of admissible cooling operations. These statements concern the limiting thermodynamic regime, rather than predicting that every substance stays in one classical phase during cooling.

The laws work together: conservation sets a device's energy budget, entropy production restricts the attainable budget, and constitutive equations give the properties of its working material. For fixed temperature and volume, Helmholtz free energy decreases toward equilibrium; at fixed temperature and pressure the corresponding potential is Gibbs free energy. Stable equilibrium, phase coexistence and Maxwell relations follow from these potentials with appropriate constraints. Calorimetry, mechanical-work measurements, thermometry and reversible-path comparisons test different parts of this framework. None of the laws by itself supplies a particular gas equation of state or microscopic heat capacity.

The microscopic magnetic response becomes a macroscopic field problem through the magnetization $\mathbf{M}$, the magnetic dipole moment per unit volume. Definein magnetostatics. Bound currents are represented by ${\mathbf{J}}_b=\nabla \times \mathbf{M}$ and surface current ${\mathbf{K}}_b=\mathbf{M}\times \hat{\mathbf{n}}$. Normal $B$ is continuous across an interface, while $\hat{\mathbf{n}}\times ({\mathbf{H}}_2-{\mathbf{H}}_1)={\mathbf{K}}_{\mathrm{free}}$. It is generally incorrect to impose continuity of every component of both fields. A linear isotropic material has $\mathbf{M}=\chi \mathbf{H}$ and ${\mu }_r=1+\chi$, but a ferromagnet usually requires a nonlinear, history-dependent constitutive relation.

Internal field differs from externally applied field. For uniform magnetization along a principal axis of an ellipsoid, the demagnetizing field is $-N_{d}M$ and $H_{\mathrm{int}}=H_{\mathrm{applied}}-N_{d}M$. A linear specimen therefore has apparent susceptibility $\chi /(1+N_d\chi )$ relative to the applied field. Shape is especially important for large susceptibility. A long closed magnetic circuit reduces stray fields; a bar magnet still has demagnetizing fields despite the absence of magnetic monopoles. Equivalent magnetic-pole descriptions are bookkeeping for magnetization divergence, not free magnetic charge.^[ref]

In paramagnetism, permanent localized moments are preferentially aligned by the field, against thermal disorder. For independent spin-one-half moments with projections $\pm \mu$, the two Zeeman energies givein the dilute weak-field limit with $B\simeq {\mu }_{0}H$. This is a Curie law; the response saturates at $n\mu$ when alignment dominates. More generally, spin $J$ gives $\chi ={\mu }_{0}n(g{\mu }_B{)}^{2}J(J+1)/(3k_{B}T)$. Interactions can produce a Curie–Weiss response $C/(T-\Theta )$ over a suitable range. Conduction-electron Pauli paramagnetism need not show the localized-moment $1/T$ law.

In diamagnetism, the field induces a magnetic moment opposing it, so the linear susceptibility is negative. The microscopic cause is the change in orbital electronic motion, and a closed-shell insulating material has no permanent atomic moment to orient. This weak response is often only mildly temperature dependent away from other transitions. It also exists as a background in materials with a stronger positive response. Perfect flux expulsion in a superconductor is a distinct collective phenomenon and should not be identified with ordinary weak atomic diamagnetism.

In ferromagnetism, exchange interactions stabilize spontaneous magnetic order below the Curie temperature. A specimen can have zero net moment while its magnetic domains remain locally magnetized in different directions. Domain-wall motion and moment rotation produce a large response; anisotropy and pinning lead to remanence, coercivity and hysteresis. Increasing field eventually approaches saturation. Heating above the ordering transition removes spontaneous order, although an applied field still induces a response. Soft magnetic cores favor easy reversal and small cyclic loss; hard permanent-magnet materials favor retained magnetization and resistance to reversal. These properties concern collective order and barriers, not merely independent moments all following a weak field.^[ref]

A practical distinction requires a controlled field sweep and a separate magnetic readout. A toroidal core with drive winding $N_1$, path length $l$ and current $I$ has approximately $H=N_{1}I/l$. A pickup winding measures ${\mathcal{E}}_2=-N_{2}AdB/dt$; integrating after calibration gives $B$, and $M=B/{\mu }_0-H$. Reversing the sweep reveals a loop whose positive loss area per volume is $\oint HdB$. An air-core/empty-holder reference, known geometry and turn count, offset correction, temperature control and slower sweeps test background, heating and rate-dependent eddy-current losses. A single force or susceptibility reading cannot establish spontaneous order or hysteresis. The original curves below are illustrative responses, not measured material data.

A fluid is described on a continuum scale by density, pressure and velocity fields. An Eulerian field is measured at fixed position; a material element changes at rate $D/Dt={\partial }_t+\mathbf{v}\cdot \nabla$. mass conservation requiresFor incompressible flow this becomes $\nabla \cdot \mathbf{v}=0$. A streamline is instantaneously tangent to velocity, whereas a pathline follows a particle through time; they coincide in a steady flow but need not in an unsteady one.

For an ideal inviscid fluid, pressure and body forces give the Euler equation. In the constant-density gravitational case,A stationary impermeable boundary requires zero normal velocity, but inviscid theory does not impose no slip. Taking its steady projection along a streamline gives $p/\rho +v^2/2+gz=\mathrm{constant}$. If the flow is also irrotational in the connected region, that constant is shared between streamlines. Inviscid does not mean irrotational: rotating ideal-fluid flows are possible. A velocity potential exists locally for an irrotational field, and in incompressible flow it obeys Laplace's equation. The Bernoulli assumptions must be checked before converting a pressure difference into a speed.^[ref]

Real fluids transmit shear stress. An incompressible Newtonian fluid with constant dynamic viscosity $\eta$ obeysThe constitutive shear stress is proportional to strain rate, and stationary solid boundaries usually impose no slip. viscosity diffuses vorticity and converts organized mechanical energy into heat; its local dissipation $2\eta e_{ij}{e}_{ij}$ is nonnegative. It cannot generally be included merely by changing the constant in a lossless Bernoulli equation.

In a long straight circular tube, steady fully developed laminar flow hasThis Hagen-Poiseuille equation follows by balancing the axial pressure gradient with viscous shear and setting $v_z(a)=0$. The parabolic profile, the fourth-power radius sensitivity and the pumping power $\Delta pQ$ are experimentally distinct consequences. Entrance regions, rough walls, non-Newtonian response or a changing temperature require additional modelling.

The Reynolds number $\mathrm{Re}=\rho UL/\eta$ compares inertia with viscosity. Small values favor creeping-flow balances; high values permit thin boundary layers, separation and turbulence, although transition depends on geometry and perturbations, not one universal critical number. At high Reynolds number an outer flow can be approximately inviscid while near-wall viscosity still determines drag. This resolves why a zero-drag ideal potential-flow model cannot predict the resistance of an ordinary moving body. In turbulent flow, fluctuating Reynolds stresses redistribute momentum and increase the required pressure drop; averaging the equations introduces closure requirements.

A controlled tube experiment measures pressure difference and collected volume per time while changing forcing. Use tracers or a velocity probe to test profile shape, allow sufficient entrance length, and monitor radius and temperature independently. Repeat with reversed flow and at several rates to test offsets and linearity. Dye mixing alone does not measure viscosity or establish a fully developed velocity profile. The original figure below separates a specified laminar velocity profile from an independent pressure/discharge 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.

Keep the problem's one-mole convention throughout. The excluded-volume parameter $b$ represents the loss of available translational volume due to short-range repulsion; it is an effective covolume, not simply the literal molecular material volume. The attraction parameter $a$ reduces the pressure below its purely repulsive value. In a mean-field picture, interaction energy is proportional to pair density, so the attractive free-energy contribution scales as $-a/V$ for one mole and the pressure correction as $-a/V^2$. Equivalently, molecular cohesive forces and the number striking a wall each supply a density factor. These parameters approximate microscopic correlations rather than measuring every pair potential exactly.

At fixed temperature, $dF=-pdV$. Substituting $p=RT/(V-b)-a/V^2$ and integrating gives the isothermal free-energy change of a Van der Waals gasThe arbitrary temperature-only term cancels. The internal-energy derivative is insteadExpansion raises the attractive internal energy toward zero. As a check, the entropy change is $R\ln [(V_2-b)/(V_1-b)]$, so $\Delta F=\Delta U-T\Delta S$. Constant $a,b$ and a homogeneous branch with $V>b$ are assumed. If a subcritical equilibrium path crosses coexistence, a Maxwell construction replaces the unstable homogeneous continuation; the displayed integral is the direct equation-of-state result, not a claim that such a continuation is physically stable.

At the Van der Waals critical point, the critical isotherm has a stationary inflection:Equating the resulting expressions for $R{T}_c$ gives $2V_c=3(V_c-b)$, henceIn reduced variables $x=V/V_c$, the critical curve is $p/p_c=8/(3x-1)-3/x^2$: it decreases through $(1,1)$ with zero first and second derivative, rather than having separate liquid/gas turning points at the critical temperature.

In the Carnot cycle, let $1\to 2$ be reversible isothermal expansion at $T_h$, $2\to 3$ reversible adiabatic expansion, $3\to 4$ isothermal compression at $T_c$, and $4\to 1$ adiabatic compression. The pressure–volume curve is clockwise. The entropy–temperature figure below uses temperature on the horizontal axis and entropy on the vertical axis, literally an $S$–$T$ plot: isotherms are vertical and adiabats horizontal. The hot isotherm increases entropy by $\Delta S$ and absorbs $Q_h=T_h\Delta S$; the cold one rejects $Q_c=T_c\Delta S$. Thus $W=(T_h-T_c)\Delta S$ and ${\eta }_C=1-T_c/T_h$.

For the second engine, choose the working direction with heat addition during $A\to B$ at $p_1$ and rejection during $C\to D$ at $p_2$. A reversible ideal-gas adiabat obeys $T\propto p^{(\gamma -1)/\gamma }$. DefineThe two constant-pressure heat transfers are $Q_{\mathrm{in}}=c_p(T_B-T_A)$ and $Q_{\mathrm{out}}=c_p(T_C-T_D)=Q_{\mathrm{in}}/r$, where rejected heat is quoted positive. Therefore the Brayton cycle efficiency isfor $\gamma =5/3$. All heat-capacity amount conventions cancel in the ratio.

Since $T_B>T_A$ and $T_C=T_B/r$, $T_D=T_A/r$, the maximum is at B and the minimum at D. The Brayton-Carnot efficiency comparison isThe inequality is strict for a nondegenerate cycle with positive heat input. Although each process can be reversible, the isobaric heat transfer spans a range of temperatures; making the whole cycle reversible requires a succession of matched-temperature reservoirs, rather than just the hottest and coldest reservoirs.

For an ideal diatomic gas with active rotation and frozen vibration, $\gamma =7/5$ andHere $0<p_2/p_1<1$. This comparison holds at the same pressure ratio and assumes temperature-independent heat capacities in the stated mode regime. Frozen rotation or excited vibration would require their actual heat capacities.

A microstate specifies the microscopic state: classically the positions and momenta of all particles, or quantum mechanically an allowed many-body state. A macrostate specifies a few coarse variables such as energy, volume and particle number. Many microstates can give the same macrostate. A statistical ensemble is a probability assignment over hypothetical copies prepared with the same macroscopic constraints; it does not require a sample to be literally divided into many physical copies. An observable's ensemble mean is $⟨A⟩={\sum }_i{p}_i{A}_i$, with an integral replacing the sum in classical phase space.

For an isolated system at fixed $U,V,N$, the microcanonical ensemble assigns equal probability to accessible microstates in a narrow energy shell. If their number is $\Omega (U,V,N)$, the Boltzmann entropy is $S=k_B\ln \Omega$. Thermal contact between two initially isolated systems allows energy to be shared while their sum stays fixed. The most probable division maximizes $S_1(U_1)+S_2(U-U_1)$, so equilibrium requires equal $\partial S/\partial U$. This gives the statistical temperature definitionFor positive temperature, energy tends to flow from lower entropy slope to higher entropy slope until those slopes agree. The narrow energy-shell width represents macroscopic resolution; it must not exclude the physically accessible states.

For a small subsystem exchanging energy with a much larger bath, the bath has energy $U_{\mathrm{tot}}-E_i$ when subsystem microstate $i$ has energy $E_i$. Expanding its entropy gives $S_b(U_{\mathrm{tot}}-E_i)=S_b(U_{\mathrm{tot}})-E_i/T+\cdots$. Counting compatible bath states therefore yields the Boltzmann distribution and canonical ensemble:The exponential weights individual microstates. An energy level of degeneracy $g_j$ has total probability $g_j{e}^{-\beta E_j}/Z$; a more degenerate excited level can outweigh one ground-state level even though each of its microstates is less probable. The Gibbs entropy $-k_B{\sum }_i{p}_i\ln p_i$ gives $S=k_B\ln Z+U/T$. The equilibrium energy variance is $⟨(\Delta E{)}^2⟩=k_B{T}^2{c}_V$, linking measurable fluctuations to a macroscopic response.

If particle exchange is also allowed, a reservoir fixes $T$ and chemical potential $\mu$. The grand canonical ensemble weights a state by $e^{-\beta (E_i-\mu N_i)}/\Xi$, with $\Xi ={\sum }_{N,i}{e}^{-\beta (E_i-\mu N)}$. Mean particle number is $k_{B}T\partial \ln \Xi /\partial \mu$, while volume remains a constraint in this version. Thus the ensemble follows the physical exchange conditions: fixed energy, thermal contact, or thermal and particle contact.

In ordinary large short-range systems, relative fluctuations often decrease as $N^{-1/2}$ and the ensembles give the same bulk thermodynamics away from singular limits. Finite systems, critical fluctuations and long-range interactions require care. Interpreting a time average as an ensemble average additionally requires appropriate sampling or mixing dynamics; it is not guaranteed solely by naming an ensemble. thermal equilibrium probability assignments do not describe every driven or relaxing state.

Use peak-amplitude phasors with time factor $e^{-i\omega t}$ and $\mu ={\mu }_0$. Faraday's law of induction and the Ampere-Maxwell law giveIn the good-conductor limit $\sigma \gg \omega ϵ$, take the root with positive attenuation into the slab:The real physical field is $\mathrm{Re}[{\tilde{E}}_y{e}^{-i\omega t}]$. A thick slab excludes a returning wave from its rear face. The skin depth is the electric-field amplitude decay length; the intensity decay length is $\delta /2$.

Faraday's law of induction fixes ${\tilde{H}}_z=k{\tilde{E}}_y/({\mu }_0\omega )$. The surface impedance of a good conductor and the complete transmitted magnetic field are consequentlyThe sign of the imaginary part belongs to the stated time convention. Average the Poynting vector using complex conjugation and the factor one half appropriate to peak phasors:As an independent energy check, $-{\partial }_x⟨N_x⟩=\sigma |E_0{|}^2{e}^{-2x/\delta }/2$, the mean Joule heating per unit volume.

Ohm's law gives ${\tilde{J}}_y=\sigma {\tilde{E}}_y$ and ${\tilde{B}}_z={\mu }_0{\tilde{H}}_z$. Hence the averaged Lorentz force per unit volume isIntegrating through the illuminated half-space gives the positive inward radiation pressureFinally let $Z_0=\sqrt{{\mu }_0/{ϵ}_0}$ be the vacuum wave impedance. Tangential-field boundary conditions at normal incidence give the electric-amplitude reflection coefficient $r=(Z-Z_0)/(Z+Z_0)$. Since $|Z|\ll Z_0$,Here $N_0=(1-R)I_{\mathrm{inc}}$ by surface energy balance. Combining the leading formulas gives $p\simeq 2I_{\mathrm{inc}}/c$, the nearly perfect-reflector limit. The displayed reflection expansion also assumes a good conductor relative to the incident vacuum wave, $\omega {ϵ}_0/\sigma \ll 1$.

Use the Euler-Bernoulli beam equation for a uniform slender elastic beam, small deflection and small slope, neglecting axial tension, shear deformation, rotational inertia of a cross-section and damping. Let $E$ be Young's modulus, $I={\int }_A{z}^{2}dA$ the second moment of area about the neutral axis, and $\rho$ the mass per unit length. If density instead denotes mass per volume, replace this $\rho$ by ${\rho }_{\mathrm{bulk}}A$. The bending stiffness is $EI$; the small-slope curvature is $y_{xx}$. The bending and kinetic energies areTwice integrating the variation of $U$ by parts gives the restoring load $-EI{y}_{xxxx}$ per unit length. Newton's second law therefore givesAt a pivot, displacement is fixed but slope is free; zero end bending moment gives $EI{y}_{xx}=0$. This also follows from the natural boundary term $EI{y}_{xx}\delta y_x$ in the variation.

With separation of variables, $y=XT$ givesChoose the time normalization $T(0)=1$. Initial rest removes the sine time term, so $T(t)=\cos \omega t$. The spatial characteristic roots are $\pm ik,\pm k$, givingThe two conditions at $x=0$ give $B+D=0$ and $-B+D=0$, hence $B=D=0$. At $L$, addition and subtraction of $X(L)=X^{\prime \prime }(L)=0$ give $C\sinh kL=0$ and $A\sin kL=0$. Nontrivial positive-frequency modes therefore have $C=0$ and $k_n=n\pi /L$. Thus the simply supported beam normal modes areThe zero-frequency cubic solution also vanishes under these four endpoint conditions, so $n=0$ is not an additional mode.

For the specified release, small slope requires $|a|L\ll 1$. Expand the initial deflection in the Fourier sine series. Its coefficients areWriting ${\omega }_1={\pi }^2\sqrt{EI/\rho }/L^2$, the resulting motion isOnly odd normal modes occur because the initial profile is symmetric about the midpoint, and their amplitudes decrease as $n^{-3}$. The frequencies grow as $n^2$, so the motion is dispersive and is not a single shape multiplied by one cosine. Nevertheless the ideal undamped solution is exactly periodic with period $2\pi /{\omega }_1$. At $t=\pi /(2{\omega }_1)$ all odd-mode cosines vanish, so the bar is instantaneously straight; at $t=\pi /{\omega }_1$ it has the opposite initial deflection. The first mode dominates displacement, with third-mode amplitude only $1/27$ of the first.

There is initial curvature incompatibility at a moment-free beam end: the prescribed parabola has $y_{xx}(x,0)=2a$, whereas a moment-free pivot requires zero curvature at an endpoint. Therefore the data do not define a classical solution with curvature continuous through the two end-time corners. They do define a finite-energy weak solution after release of the forces or end moments that prepared the initial shape. The sine expansion represents that solution: its displacement reproduces the initial parabola, each mode satisfies the released endpoint conditions, and ${\sum }_n{n}^4|A_n{|}^2<\infty$ gives finite bending energy. One must not differentiate the initial series twice at its endpoints and assign the interior parabola's curvature there. This qualification leaves the modal displacement prediction intact.

Use peak phasors with time dependence $e^{-i\omega t}$, real positive electrical conductivity $\sigma$, and absolute permeability ${\mu }_a=\mu {\mu }_0$. Ampère-Maxwell equation and Faraday equations giveso ${ϵ}^{\prime }=ϵ+i\sigma /(\omega {ϵ}_0)$. The sign of the imaginary term follows from the stated time convention. For a transverse plane wave in a homogeneous conductor,In the good-conductor limit, choose the root that decays into $z>0$:The electric phasor and physical field are consequently ${\tilde{E}}_x=E_0{e}^{(i-1)z/\delta }$ and $E_x=\mathrm{Re}[E_0{e}^{i(z/\delta -\omega t)-z/\delta }]$, as required. Tangential boundary matching sets $E_0$; it need not equal the incident amplitude.

Faraday's law fixes the magnetic field, including its phase:For real $E_0$ and $\psi =z/\delta -\omega t$, this means $H_y=E_0{e}^{-z/\delta }(\cos \psi -\sin \psi )/(\omega {\mu }_a\delta )$, with other components zero. The time-averaged Poynting vector isThe energy entering a thin slab exceeds that leaving it by $-d⟨N_z⟩/dz$ per unit volume, givingDirectly, $\mathbf{J}=\sigma \mathbf{E}$ and $\mathbf{J}\cdot \mathbf{E}=J^2/\sigma$ at each time. Therefore $⟨J^2/\sigma ⟩=\sigma |{\tilde{E}}_x{|}^2/2$ is exactly this loss. The factors one half use peak rather than root mean square phasors; amplitude decays over $\delta$, while intensity decays over $\delta /2$.