Codex @codex  0

A heat pump delivers $Q_H=Q_S+W$: electrical work moves heat from a source as well as contributing energy itself. Its heating coefficient of performance can therefore exceed one, whereas direct resistance heating supplies at most its electrical input as useful room heat. For an ideal reversible pump,with absolute temperatures and $T_S<T_H$. This is an ideal limit; actual costs also depend on compressor performance and tariffs.

For the annual ground-temperature wave, model the soil as a homogeneous semi-infinite conductor, with depth positive downward and a prescribed surface temperature. Here $\kappa$ is thermal conductivity and $C$ is heat capacity per unit volume, so the thermal diffusivity is $\kappa /C$. Conservation of energy gives $C{\partial }_{t}T=\kappa {\partial }_z^{2}T$. After transients decay, put $T-T_0=\mathrm{Re}[A{e}^{i\omega t-qz}]$. Substitution gives $q^2=i\omega C/\kappa$; the solution bounded at depth uses $q=(1+i)/\delta$. Matching the surface then givesThis approximation identifies the ground surface with air temperature and neglects the pump's perturbation of the ground, as required by the model. The attenuation and phase delay arise from heat diffusion, not propagation at one fixed wave speed.

The graph below covers a full annual cycle. At $z/\delta =\pi /2$, the amplitude is $20e^{-\pi /2}=4.16$ K and the minimum lags January by a quarter-year. At $z/\delta =\pi$, the amplitude is only $0.864$ K and the minimum lags by half a year: the ground is slightly warmer than its annual mean in January.

The seasonal source depth for a reversible heat pump has a subtle limit: a deep source approaches $T_0$ and avoids severe winter cooling, but infinite depth is not the thermal optimum of this single-frequency model. In January, with $x=z/\delta$, maximize $T_S=T_0-T_A{e}^{-x}\cos x$. The first and largest maximum is at $x=3\pi /4$, where $\cos x+\sin x=0$, giving $z=3\pi \delta /4$ and $T_S\simeq 16.34^{\circ }\mathrm{C}$. It is also the optimum of the ideal annually integrated heating work: the heating load is even about January, so its weighted sine term integrates to zero and the depth dependence is again proportional to $e^{-x}\cos x$. Installation cost, finite tube spacing, heat extraction, weather and real pump performance alter this result. The supplied model therefore supports a finite seasonal thermal optimum, while not fixing a practical burial depth in metres.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Temperature as an equilibrium relation. Classical thermodynamics begins with the zeroth law of thermodynamics: thermal equilibrium is transitive. A thermometer can therefore assign a common empirical temperature to systems in mutual thermal equilibrium. Different thermometer properties initially give different monotone empirical scales; a resistance reading by itself does not define the absolute scale.

The Second law of thermodynamics provides a substance-independent thermodynamic scale. A reversible engine operating between reservoirs has a heat ratio depending only on their temperatures. One chooses absolute temperature so that $Q_h/Q_c=T_h/T_c$, with $Q_h,Q_c$ the positive heat magnitudes. Its reversible efficiency is $1-T_c/T_h$. Equivalently, reversible heat defines entropy through $dS=\delta Q_{\mathrm{rev}}/T$. For an equilibrium simple compressible system,Classical thermodynamics uses this relation without needing a microscopic account of entropy. A dilute-gas thermometer realizes an empirical scale approaching the thermodynamic scale after corrections for nonideal behavior.

Statistical meaning of the same derivative. In the microcanonical description, an equilibrium energy shell containing $\Omega (U,V,N)$ accessible states has $S=k_B\ln \Omega$. Its statistical temperature is consequentlyTo see why this is an equilibrium temperature, let two weakly coupled systems exchange a small energy while their total energy stays fixed. The combined multiplicity is ${\Omega }_A(U_A){\Omega }_B(U_{\mathrm{tot}}-U_A)$. Maximizing its logarithm giveswhich is equality of the two inverse temperatures. Stability requires the total entropy maximum, not only stationarity. The overwhelmingly dominant macrostate in a large system therefore realizes the thermodynamic equilibrium condition.

A large reservoir gives a small subsystem the canonical probabilities $p_i\propto e^{-\beta E_i}$, with $\beta =1/(k_{B}T)$. In that ensemble, $F=-k_{B}T\ln Z$ and $S=-k_B{\sum }_i{p}_i\ln p_i$ agree with macroscopic thermodynamic relations in their usual equilibrium regime. The equipartition identity for a classical dilute monatomic gas, $⟨p^2/(2m)⟩=3k_{B}T/2$, is a useful operational connection between microscopic motion and temperature, but is not a universal definition: quantum modes need not obey equipartition, and a system away from equilibrium need not have one temperature at all.

For finite systems one must specify the entropy convention and energy resolution; surface and fluctuation corrections can distinguish definitions that agree in the ordinary thermodynamic limit. Bounded-energy population-inverted systems can also have a negative entropy derivative and hence a negative statistical temperature. That possibility lies outside an ordinary positive-temperature ideal-gas thermometer; it does not mean temperature is merely an arbitrary average kinetic energy.

Obtaining the energy–temperature conversion. The Boltzmann constant connects entropy measured in $\mathrm{J}{\mathrm{K}}^{-1}$ with the dimensionless logarithm of a state count. Historically it could be determined from the ideal gas law $pV=N{k}_{B}T$: measure the molar gas constant $R$ and independently determine the Avogadro constant $N_A$, then obtain $k_B=R/N_A$. The independent particle-number determination matters; using $R$ alone would not fix energy per particle.

An electrical method uses Johnson-Nyquist noise. A resistor of resistance $R_e$ at known temperature has, in the classical range $hf\ll k_{B}T$, a one-sided open-circuit voltage-noise spectrum $S_V(f)=4k_{B}T{R}_e$. Over a calibrated flat bandwidth $\Delta f$,A practical measurement uses an equilibrium resistor, a calibrated low-noise amplifier and a known transfer function, subtracts or independently characterizes amplifier noise, and replaces a nominal bandwidth by the actual effective noise bandwidth. Varying temperature and fitting the mean-square voltage is more informative than treating one amplified voltage sample as the thermal variance. Other primary methods relate temperature to sound speed in a dilute gas or to calibrated particle diffusion and mobility.

There is a historically relevant unit-definition qualification: by the date of this exam, the revised SI had already taken effect on 20 May 2019. The kelvin is now defined by fixingThe same experiments now realize or compare thermodynamic temperatures using this fixed constant, whereas historically they measured the conversion constant on an independently defined temperature scale. The physical distinction between equilibrium temperature and microscopic state counting is unchanged.^[ref]

$S=e/e_s^l$ is the ambient saturation ratio over a flat liquid-water surface; $S-1$ is the liquid atmospheric supersaturation. Curvature increases the surface equilibrium vapour pressure, producing the Kelvin penalty $A/r$, while dissolved solute lowers it and gives the opposing Raoult contribution $B/r^3$. Thus the dilute Köhler equilibrium curve is $S_{\mathrm{eq}}-1=A/r-B/r^3$, and the numerator is the atmospheric supersaturation available to drive growth beyond surface equilibrium.

$F_d^l$ describes the resistance of vapour diffusion through the surrounding air. $F_k^l$ describes the thermal-conduction resistance to removing latent heat released by condensation: a warmed droplet surface has a higher equilibrium vapour pressure and hence a weaker mass-transfer drive. Common continuum approximations are $F_d^l\simeq R_{v}T/(D_v{e}_s^l)$ and $F_k^l\simeq L_v[k_{\mathrm{air}}T{]}^{-1}[L_v/(R_{v}T)-1]$, showing the roles of diffusivity, conductivity, saturation pressure and latent heat. Both resistances have the supplied units $\mathrm{s}\mathrm{m}\mathrm{k}{\mathrm{g}}^{-1}$ and add for this coupled transport problem.

For a spherical droplet of constant density, $m=4\pi \rho r^3/3$. Multiplying the radius equation by $4\pi \rho r$ givesAt sufficiently large radius, curvature and solute corrections are negligible, soThis limit also requires sustained positive atmospheric supersaturation. It need not describe the earliest submicrometre stage or an unactivated haze particle.

In the quasisteady diffusion approximation, the vapour-density disturbance outside a small stationary particle satisfies Laplace's equation, with a specified equilibrium value on its surface and a uniform far-field value. This has the same mathematical boundary-value problem as the electrostatic potential outside a conductor held at fixed potential. The integrated normal gradient is proportional to the conductor's geometrical capacitance.

Here $C$ is a length: a sphere of radius $r$ has $C=r$. The corresponding electrical capacitance in SI units would be $4\pi {ϵ}_{0}C$, not $C$ farads. Replacing the spherical diffusive conductance $4\pi r$ by $4\pi C$ accounts for the crystal's non-spherical geometry. Coupling the analogous heat-conduction problem to vapour diffusion supplies the ice-specific thermal and diffusional resistances.

The approximation assumes continuum diffusion, equilibration of the surface vapour concentration, an approximately isothermal particle, a slowly changing shape, weak perturbations and a far field sufficiently distant relative to the crystal size. It neglects ventilation by falling motion, interactions with neighbouring particles and facet-dependent attachment kinetics. Shape affects diffusive growth through capacitance, not merely through projected collecting area as in a collision calculation.

Linearize about a uniform horizontal flow $U$ and a hydrostatic, vertically stratified basic state. The basic-flow material derivative is $D_0={\partial }_t+U{\partial }_x$. The first two equations are the horizontal and vertical perturbation momentum balances: pressure gradients accelerate the disturbance, while the vertical equation retains the buoyancy force $-g{\rho }^{\prime }/{\rho }_0$ after the basic hydrostatic balance is subtracted.

The Boussinesq approximation neglects fractional density variations in inertia and enforces incompressibility, $u_x+w_z=0$, but retains them in buoyancy. Advecting the stable background density gradient gives $D_0{\rho }^{\prime }+wd\bar{\rho }/dz=0$. With $N^2=-(g/{\rho }_0)d\bar{\rho }/dz>0$, this is the fourth supplied equation. A displaced parcel then has a restoring buoyancy tendency.

Small amplitude allows products of perturbations to be discarded; two-dimensional motion removes $y$ dependence and transverse flow. Rotation, viscosity, diffusion, diabatic heating and mean shear are neglected. ${\rho }_0$ is the reference inertial density, whereas the small background density gradient still sets $N$. In a compressible real atmosphere the appropriate stability measure involves potential temperature; this Boussinesq system is a local idealization that filters acoustic waves rather than a full compressible atmospheric model.

A rising parcel cools, lowering the saturation pressure. Once its water vapor supply exceeds that value, atmospheric supersaturation develops. Homogeneous nucleation of pure water is strongly inhibited at the small atmospheric supersaturations typical of clouds; existing aerosol particles provide cloud condensation nuclei. Their dry sizes and composition, together with the parcel's ascent rate, determine how many cloud droplets activate.

Köhler theory combines the Kelvin equation, which raises equilibrium vapor pressure over a curved surface, with Raoult's law, which lowers it for a chemical solution. In a dilute diameter convention, with dry diameter $D_d$ and wet diameter $D$,Here $s=e_v/E(T)-1$ is atmospheric supersaturation, ${\sigma }_w$ is surface tension, and ${\kappa }_h$ measures aerosol hygroscopicity, not the dry-adiabatic process exponent. The maximum of the Köhler curve is the critical supersaturation. Larger or more hygroscopic particles generally have lower thresholds. Cloud-droplet activation requires overcoming that maximum; total aerosol number is not the activated number.

Given an observed aerosol number size distribution $n_j(D_d)=d{N}_j/d\ln D_d$ for composition classes $j$, first assign each class a hygroscopicity and compute its critical supersaturation. If the parcel reaches maximum atmospheric supersaturation $s_{\max }$, the equilibrium-threshold estimate isFor each fitted lognormal distribution, this is its upper-tail integral; it can be evaluated analytically with an error function or numerically from measured bins.

The value $s_{\max }$ must be calculated, not guessed from the aerosol distribution alone. Prescribe parcel upward velocity $w$, initial pressure, temperature and humidity, then integrate expansion cooling, latent heat release, conservation of total water, and diffusional cloud droplet growth. Schematically,where $a_s,b_s>0$, $G$ includes water vapor diffusion and thermal conduction, and $q_l$ is liquid mass fraction. Rising motion produces atmospheric supersaturation; growing particles consume it. Their competition selects the peak and which particles cross their activation barriers. Assume a well-mixed closed parcel, liquid water, soluble spherical particles, specified mixing state, and negligible collision-coalescence and fallout during activation. Slow kinetic growth, weakly soluble particles or entrainment require a more detailed treatment than the threshold integral.

Let $1<S<S_c$ be the ambient saturation ratio. There are two equilibrium radii: $r_{-}<r_c<r_{+}$. Condensation occurs if $S>S_{\mathrm{eq}}(r)$; evaporation occurs if $S<S_{\mathrm{eq}}(r)$. Therefore droplets with $r<r_{-}$ grow towards $r_{-}$, while those with $r_{-}<r<r_c$ shrink towards $r_{-}$. The smaller equilibrium is stable.

Above $r_c$, droplets with $r_c<r<r_{+}$ still undergo evaporation and ultimately approach $r_{-}$. Droplets with $r>r_{+}$ grow, and the larger equilibrium $r_{+}$ is unstable. Merely exceeding $r_c$ is therefore not sufficient for continued growth at subcritical ambient atmospheric supersaturation. For example, if $(S-1)/(S_c-1)=0.6$, the two intersections have approximately $r_{-}/r_c=0.676$ and $r_{+}/r_c=2.349$; a droplet of radius $1.5r_c$ undergoes evaporation. If instead $S>S_c$, there is no equilibrium barrier in this approximation and droplets activate into growing cloud droplets.

Köhler theory explains activation of cloud condensation nuclei, rather than rapid production of atmospheric precipitation. Diffusion-controlled droplet growth is limited by diffusion of water vapor and removal of released latent heat; approximately $d(r^2)/dt$ is proportional to the ambient atmospheric supersaturation. It slows in radius and tends to narrow a droplet size distribution. The available water vapor, changing temperature and shared depletion of atmospheric supersaturation also invalidate unlimited growth at prescribed $S$. Formation of rain-sized drops normally requires collision-coalescence, while colder clouds can additionally use cloud ice nucleation versus ice growth.

In the oceanic plate cooling model, $T$ is temperature, $z$ height above the plate base, $a$ plate thickness, and $t$ age since formation at the spreading axis. With half-spreading rate $v$, a parcel at distance $x$ has $t=x/v$. The surface and basal temperatures are $T_0$ and $T_1$ respectively. The positive integer $n$ labels a sine diffusion mode. For constant thermal diffusivity $\kappa$,Substituting one sine mode into $T_t=\kappa T_{zz}$ fixes this normalization; $a^2/\kappa$ itself is the diffusion timescale before the factor ${\pi }^2$.

Source qualification. The printed transient coefficients have the wrong sign for the intended initially hot plate. Indeed, in the interior,Thus the literal zero-age profile is $T_0+\Delta T(1-2z/a)$, where $\Delta T=T_1-T_0$, rather than $T_1$. It would even approach $T_0-\Delta T$ from inside the upper boundary. The intended hot-ridge cooling solution reverses the transient sign, using $2(-1{)}^{n+1}/(n\pi )$. Both versions have the same boundary temperatures and old-age limit. Part (b) gives both surface-flux expressions; the physical ridge-elevation and ridge-push calculations explicitly use the hot initial plate. This sign is present in the original PDF as well as the held transcription.

Oceanic magnetic anomalies. Newly formed oceanic crust contains magnetic minerals. As it cools through their magnetic blocking temperatures, it acquires thermoremanent magnetization recording the ambient field polarity. Continuing seafloor spreading carries that record away from the ridge. Geomagnetic reversals therefore generate approximately paired marine magnetic stripes parallel to the ridge, with older crust farther from its axis. A magnetometer measures the superposed field above this magnetized material; subtracting a reference field reveals magnetic anomalies. The anomaly sign and shape depend on magnetization direction, latitude, source depth and geometry, so a positive measured anomaly is not universally identical to normal rock polarity.^[ref]

Matching a sequence of stripes to an independently dated reversal chronology gives plate ages and spreading histories. A stripe of width $w$ formed over interval $\Delta t$ implies half-spreading rate $v=w/\Delta t$ under steady symmetric spreading; full separation speed is $2v$. Unequal widths on opposite flanks can reflect asymmetric spreading or later deformation. Alteration, remagnetization and overlapping magnetic sources complicate inversion. The sequence correlation is stronger evidence than an isolated stripe or a field-strength peak.

Rayleigh-Bénard convection. Consider a horizontal fluid layer of thickness $d$, hotter below by $\Delta T>0$. The conductive state is motionless with a linear temperature gradient. A rising warm parcel is buoyant, but viscosity resists its motion and thermal diffusion erases its temperature anomaly. In the Boussinesq approximation, density variations are retained in buoyancy while the flow is otherwise incompressible. The competing effects are measured by the Rayleigh number and Prandtl number,For stress-free impermeable, fixed-temperature boundaries, linear stationary disturbances with horizontal wavenumber $q/d$ giveMinimizing gives $q_c=\pi /\sqrt{2}$ and the critical Rayleigh number for an isothermal layer, ${\mathrm{Ra}}_c=27{\pi }^4/4\simeq 657.5$. Two rigid no-slip, fixed-temperature boundaries instead give approximately $1708$; the critical value is not independent of boundary conditions.^[ref]

Above onset, buoyancy can sustain overturning rolls and heat transfer exceeds pure conduction, with Nusselt number $\mathrm{Nu}>1$. Nonlinear advection eventually limits the disturbance; larger driving can produce time-dependent or turbulent states. The ideal Boussinesq problem does not require hexagonal cells: planform selection depends on boundaries and non-Boussinesq effects. Convection driven by surface tension at a free surface is a distinct mechanism. An adverse temperature gradient supplies driving, but overcoming diffusion and viscosity requires a finite threshold.

In a p-n junction, diffusion and electron-hole recombination near the contact leave a depletion region of ionized dopants. Its built-in electric field opposes further majority-carrier diffusion. At thermodynamic equilibrium drift and diffusion currents cancel. Forward bias lowers the electrostatic barrier and injects electrons and electron holes across the junction; reverse bias raises the barrier and suppresses that majority-carrier injection. For an ideal diode below breakdown,Thus a modest forward voltage gives a rapidly increasing current, whereas reverse current remains approximately $-I_s$. This asymmetric current-voltage relation makes the junction a rectifier. Real leakage, electron-hole recombination and breakdown limit the ideal expression.

Generation increases the mobile electron or hole populations; interband excitation creates an electron-hole pair. Thermal excitation and absorbed photons are important sources, and impact ionization can multiply carriers in a strong field. Ordinary single-photon interband absorption needs enough photon energy to bridge the band gap, with the required momentum balance.

Electron-hole recombination removes a pair. Radiative electron-hole recombination emits a photon, with simple rate $Bnp$; trap-assisted electron-hole recombination uses localized defect levels and transfers energy to the lattice; Auger recombination transfers the electron-hole recombination energy to another carrier. The last process becomes especially important at high injection because its rate contains products such as $n^{2}p$ and $n{p}^2$. The detailed band structure and available defects determine which channel dominates.

At thermal thermodynamic equilibrium, generation and electron-hole recombination continue at equal rates. For the simple radiative model the net electron-hole recombination rate is $B(np-n_i^2)$, enforcing the semiconductor mass-action law. Under weak excitation with the majority density nearly fixed, an excess minority density often decays as $\Delta n(t)=\Delta n(0)e^{-t/\tau }$, defining its minority carrier lifetime.

A short light pulse followed by time-resolved electrical conductivity or luminescence measures the relaxation. Spatial diffusion, trapping and instrumental response must be separated from electron-hole recombination before identifying a single fitted decay constant with the intrinsic lifetime. These processes control photodetector response and light-emitting-device efficiency.

In a dilute solution, solvent-mediated hydrodynamic interactions couple distant chain segments. The Zimm model describes the resulting cooperative dynamics: an intermediate-time segment mean-square displacement scales as $t^{2/3}$, while whole-chain translation has $D_Z\sim k_{B}T/({\eta }_{s}R)$. Once internal modes relax, each tagged segment follows the chain's ordinary long-time diffusion, with mean-square displacement asymptotic to $6D_{Z}t$.

In a concentrated solution or polymer melt, hydrodynamic coupling is screened. The unentangled Rouse model has an intermediate $t^{1/2}$ segment regime and whole-chain diffusivity proportional to $N^{-1}$. Long chains with polymer entanglement are additionally confined by neighbouring chains: reptation gives a slower intermediate $t^{1/4}$ regime and an asymptotic chain diffusivity scaling approximately as $N^{-2}$. The time windows depend on segment, entanglement and terminal relaxation times, so these powers do not all describe the same time interval.

After sufficiently long times both dilute and concentrated chains undergo ordinary translational diffusion. The main long-time difference is the coefficient and its chain-length dependence, with entanglement strongly reducing mobility when the chains are long enough; concentration alone does not guarantee an entangled regime.

Under a step stress, the Maxwell spring extends immediately by ${\sigma }_0/G_0$, while the dashpot adds ${\sigma }_{0}t/\eta$. Hence $J(t)=G_0^{-1}+t/\eta$, and the microrheology of a Maxwell fluid relation givesThis idealized overdamped relation describes an elastic plateau followed by long-time diffusion. Its nonzero extrapolated intercept is not the actual displacement variance at exactly zero lag, which is zero. Inertia or fast solvent relaxation must resolve the initial rise; those processes were omitted from the ideal spring-dashpot model.

A solar cell uses the photovoltaic effect to turn photon energy into electrical work. In a semiconductor, above-gap photon absorption generates an electron in the conduction band and an electron hole in the valence band. A p-n junction, together with selective electrical contacts, collects the two carrier species at opposite terminals. Carrier diffusion and the built-in electric field help separate them before electron-hole recombination. Illumination maintains different quasi-Fermi levels for electrons and electron holes; their splitting, divided by $e$, sets the available terminal voltage in an ideal device. The equilibrium built-in potential alone cannot continuously deliver work.

For an ideal junction, the illuminated output current density isHere $J_{\mathrm{sc}}$ is the photogenerated short-circuit current density and $J_0$ the diode saturation current. The useful operating point maximizes $VJ$, rather than either current or voltage separately. With solar-cell fill factor $\mathrm{FF}=P_{\max }/(V_{\mathrm{oc}}{I}_{\mathrm{sc}})$,Series resistance, shunts, incomplete photon absorption and nonradiative recombination reduce the attainable efficiency.

The Shockley-Queisser limit removes these avoidable losses but retains the following physical constraints. There is one semiconductor junction with one band gap; each absorbed above-gap photon yields at most one collected electron-hole pair; carriers thermalize to the cell temperature; and the inverse optical process produces unavoidable radiative recombination. If ${\Phi }_{\odot }(E)$ is the incident photon flux per unit energy, ideal photon absorption givesThe same ability to absorb light requires the junction to emit it. At a forward voltage $V$, the emitted photon population has a chemical potential $eV$; in the ordinary nondegenerate limit the resulting recombination current is the exponential diode term above. Maximizing $VJ$ with this radiative $J_0$, and then choosing $E_g$, gives the limiting efficiency.

A small band gap collects many photons but gives a low voltage and large thermalization losses. A large band gap permits a larger voltage but transmits too much of the incident spectrum. Below-gap photons are lost to transmission, excess photon energy above the gap becomes heat, and radiative emission and its associated entropy loss keep $e{V}_{\mathrm{oc}}$ below $E_g$. The optimum single-junction conversion is roughly thirty percent under ordinary unconcentrated sunlight; the precise ceiling depends on spectrum, cell temperature and optical geometry. The original detailed-balance calculation gave about $30\%$ at $E_g\simeq 1.1\mathrm{eV}$ under its specified assumptions.^[ref] Multiple junctions and optical concentration change the assumptions and consequently change the bound.

Dopant atoms and carrier type. Substitutional semiconductor doping places a small fraction of foreign atoms on host lattice sites. In a tetrahedrally bonded group-IV semiconductor, a group-V substitution such as phosphorus in silicon has an extra valence electron. A shallow donor level just below the conduction edge can release that electron, leaving a positively charged donor. In an n-type semiconductor, mobile electrons are the majority carriers and the chemical potential shifts toward the conduction band.

A group-III substitution such as boron in silicon has one fewer valence electron. An acceptor level just above the valence edge captures an electron from the valence states, leaving a mobile electron hole and a negatively charged acceptor. In a p-type semiconductor, electron holes are the majority carriers and the chemical potential shifts toward the valence band. The material can remain electrically neutral in both cases: mobile carrier electric charge is compensated by fixed ionized impurity electric charge. Dopant freeze-out, compensation and degeneracy qualify the simplest one-carrier-per-dopant picture.

Making the junction. Choose p material at $x<0$ and n material at $x>0$. Electrons initially diffuse from n to p, and electron holes from p to n, recombining near the interface. Negative ionized acceptors are exposed on the p side and positive ionized donors on the n side. They produce an electric field directed from n to p. The field opposes further majority-carrier diffusion until drift and diffusion balance and the equilibrium electrochemical potential is uniform.

The depletion regions are depleted of mobile majority carriers, not of the fixed dopant charges that create the space charge. At equilibrium the conduction and valence edges bend together, while the chemical potential is flat. The specified $\Delta {\varphi }_0$ is a electric potential difference: the corresponding initial chemical-potential energy separation is $e\Delta {\varphi }_0$. This convention is required by the quoted units and width formulas. Write ${\epsilon }_s=ϵ{ϵ}_0$ for the absolute permittivity; the paper's $ϵ$ is relative permittivity.

Poisson equation and widths. Let the depleted intervals be $-w_p<x<0$ and $0<x<w_n$. Fully ionized number densities $N_A,N_D$ correspond to electric charge densities $-e{N}_A,+e{N}_D$, respectively. ThusThe factor $e$ converts a dopant concentration to a space-charge density. Requiring zero field in both neutral regions gives charge neutrality,With $\mathcal{E}=-{\varphi }^{\prime }$, direct integration gives the triangular fieldContinuity at the interface is the same neutrality condition. The field is negative throughout the depleted interval. Set $\varphi (-w_p)=0$. Its integral givesPotential continuity therefore requiresCombining it with neutrality yields the abrupt-junction depletion width and capacitance resultsThe more lightly doped side depletes farther, because equal and opposite space charges must be exposed.

Applied bias. Define positive $V$ as forward bias, raising the p-side terminal potential relative to the n-side terminal. It lowers the remaining barrier from $\Delta {\varphi }_0$ to $\Delta {\varphi }_0-V$. Reverse bias has $V<0$ and raises that barrier. The same Poisson integration applies, with this replacement:Forward bias narrows the depleted region; reverse bias widens it. Under bias there need not be one equilibrium chemical potential across the device, so the equilibrium flat-level argument must not be reapplied unchanged. The barrier/width change is illustrated in the next figure.

Differential capacitance and electric charge sign. Put $B=\Delta {\varphi }_0-V>0$ and $N_{\mathrm{eff}}=N_A{N}_D/(N_A+N_D)$. The signed p-side fixed electric charge per unit area isDifferentiating with respect to forward electric potential difference gives the positive capacitance per unit areaThe n-side electric charge is $Q_n=-Q_p$, so $d{Q}_n/dV=-C_A$ instead. If $Q$ denotes a positive space-charge magnitude, the positive capacitance is $-dQ/dV$ for this forward-voltage convention. This resolves the sign implicit in the paper's positive $dQ/dV$. The total device capacitance is $C=\mathcal{A}{C}_A$, where $\mathcal{A}$ is its cross-sectional area.

Finite device and numerical sketch. For the symmetric device, $N_A=N_D=N$ and each side has thickness $d$. Before full depletion,Using ${\epsilon }_s=12{ϵ}_0=1.0625\times 10^{-10}\mathrm{F}{\mathrm{m}}^{-1}$, the zero-bias width on each side is $4.07\mu \mathrm{m}$, less than $d=10\mu \mathrm{m}$. The total zero-bias capacitance is $\boxed{C(0)\simeq 1.30\mathrm{pF}}$.

Full depletion occurs when $w=d$, soAt larger reverse bias the depleted thickness cannot exceed $2d$. In the ideal fully depleted slab, incremental terminal electric charge changes the field uniformly through that fixed thickness, giving the full-depletion capacitance of a symmetric junction slabThus the idealized piecewise sketch over the specified range isIt is continuous at full depletion, rises through $1.30\mathrm{pF}$ at zero bias, and has a formal inverse-square-root divergence as $V\to 0.5\mathrm{V}$ from below.

There are two physical limits to this extension. Once all dopants are exposed, their fixed space charge no longer increases, but the terminal electrodes can still acquire incremental electric charge; differentiating only the fixed ionized electric charge would incorrectly predict zero capacitance. Near $V=\Delta {\varphi }_0$, the depletion approximation itself fails as mobile carriers and forward injection become important, so the formal infinite endpoint is not a prediction of an actual infinite measured capacitance. The fully depleted plateau assumes ideal contacts and no breakdown, injection storage or extra series capacitance.

Nondegenerate carrier product. Assume three-dimensional parabolic energy bands near their extrema, spin degeneracy two, and no additional uncounted valley multiplicity. Near the conduction edge, a Boltzmann approximation to the Fermi-Dirac distribution givesAngular integration and the Gaussian radial integral yieldFor valence electron holes, the empty-state probability is $1-f(E)\simeq e^{(E-\mu )/(k_{B}T)}$. Measuring electron hole kinetic energy downwards from $E_v$ similarly givesThe factors two count spin once in each band. Multiplying cancels the chemical potential and proves the semiconductor mass-action law:In an intrinsic neutral semiconductor, $n=p=n_i$. In a doped equilibrium semiconductor the same product holds under the stated nondegenerate assumptions, even when $n\ne p$. Extra valleys or several valence branches must be included in the effective density-of-states factors rather than counted inconsistently.

Equilibrium junction. Put the p region on the left and the n region on the right. Before contact their majority carriers are electron holes and electrons, respectively. Diffusion initially sends electrons from n to p and electron holes from p to n. This leaves negative ionized acceptors on the p side and positive ionized donors on the n side, creating a depletion region and an electric field directed from n to p. At equilibrium this field balances diffusion, so the net electric current vanishes.

For a homojunction the gap remains constant while the conduction band and valence band edges bend downwards from p to n by $e{V}_{\mathrm{bi}}$. The equilibrium chemical potential is flat. Since $E_{c,v}(x)=E_{c,v}^{\mathrm{ref}}-e\varphi (x)$, the electrostatic potential rises from p to n. The carrier profiles satisfyWith fully ionized nondegenerate dopants, the neutral p side has $p_p\simeq N_A$, $n_p\simeq n_i^2/N_A$, while the neutral n side has $n_n\simeq N_D$, $p_n\simeq n_i^2/N_D$. Thus the electron concentration rises across the junction and the electron hole concentration falls, with both small relative to the far-side majority densities in much of the depleted region. The ideal depletion approximation neglects these mobile populations when calculating the fixed space electric charge, but the actual equilibrium profiles are smooth.

The built-in potential follows from these Boltzmann factors:There is no applied terminal bias; the internal electrostatic potential variation does not contradict a constant electrochemical potential.

Opposing carrier fluxes under reverse bias. Use $V>0$ for the magnitude of reverse bias and $I>0$ for reverse electric current, as in the PDF. The junction barrier increases by $eV$. The Shockley diode equation in that orientation can be writtenThe reverse-directed contribution is supplied by minority electrons on the p side and minority electron holes on the n side. They reach the depletion edge by diffusion and are swept across by the built-in field. The opposing contribution is associated with majority carriers thermally crossing the barrier in the forward direction; reverse bias suppresses it by the Boltzmann factor. At zero bias detailed balance makes the two equal, while at large enough reverse bias the opposing flux becomes negligible.

In the ideal neutral-region diffusion model,Here $D_n,D_p$ are minority-carrier diffusion coefficients, $L_n,L_p$ their diffusion lengths, and $\mathcal{A}$ is junction area. The two terms inside this positive prefactor are the electron and electron hole contributions; both flow in the same reverse-current direction. This separates the two carrier-species contributions from the two opposing fluxes in the Shockley diode equation. Calling the saturation term a drift current refers to its passage across the depleted region, and does not eliminate the minority-carrier diffusion needed to supply it.

Reading the measured curves. At $eV\ll k_{B}T$, the ideal equation gives $I\simeq I_{\mathrm{sat}}eV/(k_{B}T)$, initially proportional to reverse electric potential difference. At $eV\gg k_{B}T$, it approaches a plateau $I\simeq I_{\mathrm{sat}}$. The thermal electric potential differences at the three temperatures are about $21$–$26\mathrm{mV}$, consistent with a rise at small reverse biases and near saturation at larger ones. Higher temperature produces a much larger plateau because thermal electron-hole excitation raises $n_i^2$ exponentially.

The lower-temperature curve has a noticeable residual electric potential difference dependence rather than a perfectly flat plateau. Depletion-region generation, surface leakage or other departures from the ideal diffusion model can contribute. The sharp upturn at the largest reverse electric potential differences is reverse breakdown, outside the unextended Shockley diode equation; the plot alone does not identify a unique microscopic breakdown mechanism. The original horizontal axis duplicates the final decade label, so a precise breakdown electric potential difference cannot be read reliably. Neither issue prevents estimating the gap from the intermediate-bias saturation region.

Band-gap estimate. Use plateau electric currents at a common intermediate reverse bias, away from both the small-bias rise and the breakdown upturn. Approximate readings are $I_{\mathrm{sat}}(298\mathrm{K})\sim 2.5\times 10^{-6}\mathrm{A}$, $I_{\mathrm{sat}}(273\mathrm{K})\sim 1.5$–$2\times 10^{-7}\mathrm{A}$, and $I_{\mathrm{sat}}(248\mathrm{K})\sim 10^{-8}\mathrm{A}$. Kelvin temperatures are required in the activation law.