Codex @codex  0

An equation of state relates thermodynamic variables within a material model; a phase diagram records which equilibrium phase or phases are stable under specified controls. For a pure substance at fixed temperature and pressure, stability means minimizing Gibbs free energy. Two phases coexist when their molar chemical potentials agree. The equation of state and the free-energy stability condition must therefore be used together: an algebraic relation between $p,V,T$ alone need not identify the stable branch.

A pressure-temperature phase diagram contains single-phase regions separated by phase coexistence curves. For a one-component system the Gibbs phase rule gives $F=3-P$: a single phase has two independent controls, two-phase coexistence lies on a curve, and three phases coexist at a triple point. The liquid-gas coexistence curve ends at a liquid-gas critical point; beyond it a continuous path connects liquidlike and gaslike states without crossing a first-order coexistence boundary. A fixed $p,T$ point on a coexistence curve does not specify the proportions of the two phases: an additional extensive constraint, such as total volume, fixes the phase fraction.

At a first-order phase transition, first derivatives of Gibbs free energy, such as molar entropy or volume, are discontinuous. The conversion has latent heat $L=T\Delta s$. At a continuous transition the first derivatives match, while response functions associated with higher derivatives can be singular or discontinuous. At an ordinary liquid-gas critical endpoint, the two coexisting densities approach one another and the latent heat tends to zero; the material does not retain a sharp liquid-gas boundary above it.

The coexistence slope follows directly from equality of the two molar chemical potentials. Since $d\mu =-sdT+vdp$, differentiating ${\mu }_{\alpha }={\mu }_{\beta }$ along the boundary givesThis is the exact Clapeyron equation, also commonly called the Clausius-Clapeyron relation. For vaporization, both $L$ and $\Delta v$ are positive, so the saturation pressure rises with temperature. For ordinary ice melting into liquid water, $L>0$ but the molar volume decreases, giving the familiar negative melting-curve slope. The sign follows from the volume change, not from a different law.

For a dilute ideal vapour with $v_g\gg v_l$, $\Delta v\simeq RT/p$. The ideal-vapor Clausius-Clapeyron approximation becomesIf molar latent heat is approximately constant over the interval,These approximations deteriorate close to the critical point, where liquid and gas volumes become comparable, or over a range in which latent heat changes appreciably.^[ref]

The ideal gas law $pV=nRT$ is an equation of state for a dilute gas, but it has no liquid branch and cannot predict condensation. A simple model including repulsion and attraction is the Van der Waals molar equation,where $b>0$ represents excluded molar volume and $a>0$ represents attraction. At its critical point, $\partial p/\partial v={\partial }^{2}p/\partial v^2=0$, givingBelow $T_c$, the homogeneous isotherm has an interval with $\partial p/\partial v>0$. That would have negative isothermal compressibility and is mechanically unstable. Equilibrium liquid-gas coexistence replaces the loop by a constant-pressure segment. Equality of chemical potentials gives the Maxwell construction:The equal-area condition enforces the common tangent to the appropriate free energy; it is a stability construction rather than permission to follow the unstable loop physically. The average molar volume between $v_l$ and $v_g$ varies by changing the phase fraction.

The model captures a liquid-gas critical point and qualitative coexistence, but its simple constants are not a complete quantitative description of real fluids and it has no solid phase. A full phase diagram needs additional solid free energies and can contain several crystal phases. Conversely, an experimentally determined coexistence curve and latent heat constrain equations of state through the Clapeyron equation. This is the useful connection: local thermodynamic response, global free-energy stability and the geometry of phase boundaries describe the same material at different levels.

Entropy at absolute zero. The third law of thermodynamics concerns the low-temperature limit of equilibrium states. In the zero-entropy formulation, the entropy of a perfect crystal with a nondegenerate ground state tends to zero as $T\to 0$. The statistical interpretation is $S=k_B\ln \Omega$: a unique accessible ground state has $\Omega =1$. The entropy-difference formulation of the third law of thermodynamics states that entropy differences associated with regular changes of mechanical parameters tend to zero as absolute zero is approached. These statements require equilibrium and appropriate ground-state assumptions; a frozen disordered material or a degenerate ground state can carry residual entropy.

For a phase of fixed composition and pressure,A finite limiting entropy constrains this integral. Ordinary equilibrium heat capacities therefore vanish as temperature approaches zero, with familiar examples $C\propto T^3$ for three-dimensional acoustic phonons and $C\propto T$ for a normal degenerate Fermi liquid. A constant positive heat capacity extrapolated to zero would give a logarithmically divergent entropy difference and is impossible in this regime. Discontinuous transitions and anomalous limits need separate treatment rather than extending a high-temperature approximation indefinitely.

Gibbs free energy. For fixed particle numbers,Thus, when $S\to 0$, the Gibbs free energy approaches its zero-temperature value with a horizontal tangent:The equality $G(0,p)=H(0,p)$ uses the enthalpy $H=U+pV$; neither limiting value needs to be zero. The low-temperature Gibbs free-energy expansion makes the approach quantitative. For $C_p=a{T}^r$ with $r>0$ and $S(0)=0$,The second derivative is ${\partial }_T^{2}G=-C_p/T\le 0$, although it need not tend to zero. For a reaction between such equilibrium phases, ${\partial }_T\Delta G=-\Delta S\to 0$, and $\Delta G$ and $\Delta H$ approach one another. If a genuine residual entropy $S_0$ is retained instead, $G=G(0,p)-S_{0}T+\cdots$ has limiting slope $-S_0$. A horizontal tangent cannot be asserted without the zero-entropy condition.

Why helium remains a liquid. Helium atoms attract each other weakly and are exceptionally light. Their zero-point energy is large compared with the binding that would localize them in a solid. Both helium-4 and helium-3 can therefore remain liquid down to the lowest temperatures at low pressure; substantial applied pressure is needed to solidify them. This is not a violation of the third law. A quantum liquid can have a unique low-entropy ground state, and zero-point motion does not itself imply thermal entropy. The third law does not require every substance to crystallize at atmospheric pressure.

A helium-4 atom has integer spin zero and is a boson. Below the lambda point near $2.17\mathrm{K}$ at saturated vapor pressure, its liquid is superfluid. Its strong interactions mean that it is not quantitatively an ideal Bose gas, and the condensate fraction need not equal the superfluid fraction. The two-fluid model of superfluidity separates an entropy-carrying normal component of a superfluid from the superfluid component, whose entropy is zero. At sufficiently low temperature the dominant thermal excitations are acoustic phonons: their linear dispersion in three dimensions gives $U-U_0\propto V{T}^4$, $C_V\propto V{T}^3$ and $S\propto V{T}^3$. Higher-energy rotons are exponentially suppressed. Thus the normal component of a superfluid and entropy disappear as $T\to 0$, while superfluid flow and quantum motion can persist.

A helium-3 atom has spin $1/2$ and is a fermion. Well below its Fermi temperature, but above its superfluid transition, liquid helium-3 is a degenerate interacting Fermi liquid. Only quasiparticles close to the Fermi surface are thermally excited, giving $C_V\propto T$ and $S\propto T$. At temperatures of order a millikelvin, helium-3 also becomes superfluid through formation of Cooper pairs of fermions. This is not Bose-Einstein condensation of individual helium-3 atoms. Pairing reorganizes the low-energy excitation spectrum; in the fully gapped bulk B phase, the fermionic quasiparticle contribution to heat capacity is exponentially suppressed at temperatures well below the pairing gap. Gapless collective modes can still give thermal terms following a power law. The entropy nevertheless tends to zero. The normal Fermi-liquid law cannot be extrapolated through its pairing transition unchanged.^[ref]

Helium-3 also illustrates the importance of nuclear spin entropy. Over an intermediate low-temperature range, localized spins in the solid can retain approximately $k_B\ln 2$ of entropy per atom, while the degenerate liquid has much less. Reversible conversion of liquid to the higher-entropy solid then absorbs heat; adiabatic compression along the melting curve can cool the sample. This is Pomeranchuk cooling. At still lower temperature, interactions between the nuclear spins remove this nearly free-spin entropy, so a constant $k_B\ln 2$ cannot be extrapolated to absolute zero.^[ref]

Liquefaction and cooling methods. A gas must be below its liquid-gas critical point to form a distinct liquid phase by compression. Compression above the critical temperature alone cannot achieve ordinary liquid-vapour coexistence. Practical liquefaction removes the compression heat, cools the compressed gas in a heat exchanger and uses expansion to lower its temperature further.

A throttling valve operates approximately at constant specific enthalpy when kinetic-energy changes and heat leakage are negligible. The Joule-Thomson coefficient isusing consistent molar quantities, or consistent total quantities for a fixed amount. A pressure drop cools the gas if ${\mu }_{\mathrm{JT}}>0$. An ideal gas has zero Joule-Thomson coefficient; real intermolecular interactions are essential. In particular, room-temperature helium warms on throttling, so it must first be brought into a region of positive Joule-Thomson coefficient. Expansion through an engine or turbine can instead remove work and produce approximately isentropic cooling.

A regenerative liquefier sends the cold returning gas through a counterflow heat exchanger to precool the incoming compressed gas before its next expansion. Expansion engines and external precooling can supply the initial cooling; repeated operation then produces liquid. Modern helium liquefaction commonly combines work-producing expansion with a final throttle. ^[ref]

Once liquid exists, pumping away vapour lowers its saturation pressure and boiling temperature. The escaping vapour carries latent heat away, cooling the remaining liquid. At still lower temperatures, a dilution refrigerator exploits the different entropy carried by helium-3 in concentrated and dilute helium mixtures; adiabatic demagnetization and Pomeranchuk cooling provide other routes.

Approaching zero. The unattainability formulation of the third law states that regular cooling operations with finite resources cannot reach absolute zero in finitely many stages. As the entropy available to remove becomes small and thermal couplings weaken, each method loses its ability to take the final step. An isentropic cooling trajectory starting at nonzero entropy cannot reach a common zero-entropy endpoint through regular finite changes of parameters. A fully general finite-time claim also requires assumptions about available controls and heat-transfer rates. Quantum liquids and refrigeration methods can approach $0\mathrm{K}$ while remaining consistent with vanishing equilibrium entropy; they do not provide a finite-operation route to absolute zero.

The relation describes the slope of the phase coexistence curve of a pure substance in the pressure-temperature plane:Here $p(T)$ is the common equilibrium pressure of the two phases, $T$ is their common absolute temperature in kelvins, and the derivative is taken along their coexistence curve. Choose a direction of transformation from phase 1 to phase 2. Then $L$ is the reversible latent heat for that transformation and $\Delta V=V_2-V_1$ is the corresponding volume change, using the same amount of substance for both quantities.

In the molar convention used below, $L=h_{m,2}-h_{m,1}$ has units $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$ and $\Delta V=v_{m,2}-v_{m,1}$ has units ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$, so the right-hand side has units $\mathrm{Pa}{\mathrm{K}}^{-1}$. For vaporization, both $L$ and $\Delta V$ are positive. The displayed general relation is also called the Clapeyron equation; the ideal-vapor approximation used in (c) gives its usual Clausius form.

At coexistence, the two phases of a pure substance have equal chemical potential, or equivalently equal molar Gibbs free energy:For either phase, $d\mu =-s_{m}dT+v_{m}dp$. Move an infinitesimal distance along the coexistence curve, where the equality must remain true. Subtracting the two differentials givesA reversible transformation at this common temperature absorbs latent heat $L=T\Delta s_m$. ThusThe derivation assumes thermodynamic equilibrium between two bulk phases of a single pure substance, with the same temperature and pressure, and only ordinary pressure-volume work. Surface-curvature shifts, additional external-field work and composition changes are omitted. It applies to a first-order phase boundary where the latent heat and volume discontinuity are defined. No ideal-gas assumption and no constant-$L$ assumption are needed for this differential relation. At the critical endpoint both discontinuities vanish, so their ratio requires a limiting treatment rather than direct substitution of two zeros.

For liquid-vapor equilibrium, use molar quantities. At low vapor pressure the gas is dilute, approximately ideal, and has a much larger molar volume than the liquid. HenceSubstitution into the Clapeyron equation givesIf the molar latent heat varies sufficiently slowly to be treated as constant over the temperature interval, integration givesThe positive constant $p_0$ has pressure units and is an integration constant. It is not necessarily a real infinite-temperature pressure, since the approximations fail before that limit.

For mercury between $300$ and $500\mathrm{K}$, the original plot places this range above the triple point and below the normal boiling and critical temperatures. Its vapor pressure is low relative to the critical pressure, making a dilute-gas description reasonable; the vapor's molar volume greatly exceeds that of the dense liquid. The interval is also far from the critical region where $L$ tends to zero, so a roughly constant latent heat is a useful first approximation. The approximations are not exact: Kirchhoff's law of thermochemistry gives, at fixed pressure, the temperature derivative of a transformation enthalpy as the difference between the two molar heat capacities, so a precise fit over a wide interval would allow $L$ to vary.

Taking a derivative of the logarithmic form gives a direct way to use the plot's base-ten vertical axis:Alternatively, use two nearby readings and eliminate the unknown $p_0$:Reading the original PDF graph approximately gives ${\log }_{10}(p/\mathrm{Pa})\simeq -0.6$ at $300\mathrm{K}$ and $\simeq 3.8$ at $500\mathrm{K}$. These bracket the requested temperature and give an average latent heat near $400\mathrm{K}$:Thus a suitable graph-based estimate isThe graph does not warrant several significant figures; a local tangent or different reasonable nearby readings give a comparable estimate. The factor $\ln 10$ is essential when using a base-ten logarithmic axis.

Below cloud base, an initially unsaturated closed parcel expands as ambient pressure falls. In the absence of heat exchange its temperature falls approximately at the dry adiabatic lapse rate $g/c_{p,d}$. Its vapour mixing ratio is conserved, but its saturation mixing ratio decreases strongly as it cools. Hence relative humidity rises until the lifting condensation level is reached. The density normally decreases through expansion; the environmental profile determines the parcel's buoyancy and thus whether ascent can continue.

Near saturation, slight atmospheric supersaturation permits suitable cloud condensation nuclei to activate. Köhler theory describes the competition between curvature, which raises equilibrium vapour pressure, and dissolved solute, which lowers it. Cooling generates supersaturation, while condensation onto the new droplets consumes it. Their balance normally keeps in-cloud relative humidity close to, but not identically at, saturation. The activated droplet number depends on the aerosol population and the maximum supersaturation produced by the ascent.

Above cloud base, condensation lowers the vapour content and raises the liquid-water content. The latent heat released partly offsets expansion cooling, so the moist adiabat has a smaller lapse-rate magnitude than the dry one. In a dilute approximation,A closed adiabatic parcel conserves total water $q_t=q_v+q_l$ and approximately its moist entropy; retaining condensate makes this a reversible-cloud idealization, rather than a precipitating pseudoadiabat. In a shallow cloud, the saturation-mixing-ratio decrease and slow density variation yield nearly constant condensation per height, so liquid-water content is approximately linear above base. Once activation is complete, assuming nearly fixed droplet number per volume gives increasing droplet size as condensate accumulates. Entrainment, continued activation, collision-coalescence and atmospheric precipitation can all violate that simplified description; they are excluded from the stated retrieval model.

Rename the Antoine equation coefficients $A_A,B_A,C_A$ to distinguish them from the heat-capacity constants. The logarithmic slope is independent of the pressure unit:Using an ideal vapour and neglecting liquid specific volume, Clausius–Clapeyron gives $d\ln p_{\mathrm{sat}}/dT=L_v/(R_v{T}^2)$. Therefore the positive magnitude of the latent heat released on condensation isFor $C_A=-198.043\mathrm{K}$, $d\ln L_v/dT=2C_A/[T(T+C_A)]<0$ over the quoted range. With $R_v\simeq 461.5\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$, this fit implies approximately $2.984\mathrm{MJ}\mathrm{k}{\mathrm{g}}^{-1}$ at $380\mathrm{K}$ and $1.598\mathrm{MJ}\mathrm{k}{\mathrm{g}}^{-1}$ at $573\mathrm{K}$. The required ratio does not depend on $R_v$:The decreasing trend is qualitatively consistent with the liquid and vapour phases becoming less distinct towards the liquid-gas critical point. However, the full Clapeyron relation is $L_v=T(v_v-v_l)d{p}_{\mathrm{sat}}/dT$: at criticality the two specific volumes merge and the latent heat vanishes. Neglecting $v_l$ and treating the vapour as ideal then fail. The Antoine coefficients are not licensed for extrapolation beyond $573\mathrm{K}$, and derivatives of an approximate pressure fit can be less accurate than the fitted pressures themselves.

A concrete check exposes the extrapolation failure: at the supplied $647\mathrm{K}$ critical temperature the formula would still give $L_v\simeq 1.42\mathrm{MJ}\mathrm{k}{\mathrm{g}}^{-1}$, rather than zero, and $p_{\mathrm{sat}}\simeq 134\mathrm{bar}$ rather than $221\mathrm{bar}$. Above the true liquid-gas critical point there is no liquid–vapour coexistence curve or condensation latent heat in the same sense. The last panel of the original figure in solution makes this limitation visible.

$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.

Let $R_A,R_B$ be the specific specific gas constants and $c_{pA},c_{pB}$ the specific heat capacities of the noncondensible gas and vapour. Write $e_s(T)$ for the saturation vapor pressure, $r$ for vapour mass per unit dry-gas mass, and $ϵ=R_A/R_B$. For a dilute ideal gas mixture,The second equation is the Clausius-Clapeyron equation, neglecting condensed-phase volume and taking constant latent heat $L$. The first law of thermodynamics for saturated ascent, with condensate removed and negligible condensate heat capacity, is to leading orderSince $dr=r[LdT/(R_B{T}^2)-d\ln p]$, the resulting moist adiabat has logarithmic slopeConsequently a sufficient dry adiabat criterion isThe first two conditions justify neglecting the vapour's effects on the gas law and sensible heat capacity; the last two suppress latent heating. With $x=R_A/c_{pA}$ and $\beta =L/(R_{B}T)$, the latent conditions are $\beta e_s/p\ll 1$ and $x{\beta }^2{e}_s/p\ll 1$. Usually $x\beta \gg 1$, so the squared-latent term gives the stronger condition. At sufficiently low surface temperature, the exponential smallness of $e_s(T_s)$ makes these conditions hold andThese are local conditions along the column. They can be checked at the surface when $x{\beta }_s>1+2x$: on this dry profile the strongest latent parameter is proportional to $\exp [(1+2x)\zeta -{\beta }_s(e^{x\zeta }-1)]$ and decreases upwards. Without that restriction a surface-only test need not bound every altitude. If suspended condensate is retained, its own heat capacity and loading must also be negligible; saturation alone does not specify those quantities.

For ordinary condensibles, $L/(R_{B}T)$ is typically of order ten or more near low atmospheric temperatures. At the approximate triple-point temperatures, the supplied substances give $\beta$ about $20$ for water and $10$–$12$ for carbon dioxide, methane and nitrogen. This is consistent with the supplied large ratios of latent heat to sensible heating: for molecular gases $c_p/R$ is only a few. The relevant $x=R_A/c_{pA}$ is typically a few tenths.

The integral is concentrated where ${\beta }_{s}x\zeta =O(1)$. There $x\zeta =O({\beta }_s^{-1})\ll 1$, so $e^{x\zeta }-1\simeq x\zeta$ givesThis is an Laplace's method, rather than an assumption that vapour is uniformly mixed. The effective contributing depth in log-pressure is $\Delta \zeta \sim 1/(x{\beta }_s)$, explaining why the near-surface atmosphere dominates.

An exact representation within the same dilute model makes the error explicit. Set $y={\beta }_s{e}^{x\zeta }$ and define the exponential integral $E_1(b)={\int }_b^{\infty }{e}^{-y}dy/y$. ThenThe saturated vapor column on a dry adiabat estimate overstates the exact dilute-model column by roughly $1/{\beta }_s$, about five to ten per cent for those illustrative values. It is not accurate for a warm, nondilute steam atmosphere.

Notation. The source uses $p_s$ for both the hypothetical dry surface pressure and the actual surface pressure after evaporation. Distinguish them as $P_d$ and $P$, respectively, and interpret the specified pressure thickness $\Delta p$ as fixed. This conserves the original hydrogen inventory and permits $P$ to vary. Reading both occurrences as one fixed number would be inconsistent with adding water vapor to a closed atmospheric column mass budget.

Let $E(T)$ denote the saturation pressure of water, let $x$ be the mole fraction of water vapor in the mixed layer, and let $M_d,M_v$ be the molar masses of hydrogen and water. Surface phase equilibrium and the ideal gas law giveBecause the layer has constant temperature, its relative humidity isThus every level strictly above the surface is subsaturated: hydrostatic balance makes $p<P$, whereas the saturation pressure stays constant. The result does not require a dilute water vapor approximation.

The Clausius-Clapeyron relation, neglecting the specific volume of liquid water and treating the water vapor as an ideal gas, givesHere $T_t,E_t$ are the triple point values, $L_v$ is the specific latent heat, and $R$ is the molar gas constant. All temperatures in these expressions are absolute. As throughout the question, this is a constant-latent heat, ideal gas model, not an extrapolation valid through the real critical point of water.

For a well-mixed ideal gas mixture, its mean molar mass and water vapor mass fraction areBy hydrostatic balance, the mixed layer has total mass per unit area $\Delta p/g$, of which $w_v\Delta p/g$ is water vapor. The remaining atmosphere contains no water vapor. Consequently conservation of hydrogen mass givesEquivalently the requested quadratic equation isFor each $T_s$, calculate $E$, take the positive root of this quadratic equation, and retain it only if $0\le x<1$ and $P>\Delta p$. The last condition leaves a nonempty dry layer. In particular the dry constituent's surface partial pressure is not generally $P_d$: $P_d$ measures its column weight, and the hydrogen and water vapor have different molar masses.

A convenient parameter is the water vapor mole fraction itself. The mass conservation calculation yieldsBoth functions increase for $0\le x\le 1$. At the endpoint,Therefore the physically admissible hydrogen–water vapor mixture requiresBeyond this bound, the formal positive quadratic equation root requires $x>1$, equivalently a negative hydrogen partial pressure. It is the physical constraints, rather than loss of real roots, that end the solution branch. For a concrete example, put $P_d=\Delta p=1$ in consistent pressure units, $ϵ=9$, and $E=3$. The positive root of $P^2+23P-51=0$ is about $2.037$, giving $x\simeq 1.473$, which is impossible.

The limiting temperature is obtained from $E(T_{\ast })=P_d+\Delta p$:This finite limit exists above $T_t$ if $E_t<P_d+\Delta p<E_t\exp [L_v/(R_v{T}_t)]$. The upper qualification matters mathematically: the stipulated constant-latent heat formula has a finite limiting saturation pressure as $T\to \infty$, so for an arbitrarily large unspecified $P_d+\Delta p$ the printed assertion of a finite breakdown temperature is not unconditional. For a concrete formal counterexample, choose $P_d+\Delta p=10^{15}\mathrm{Pa}$: with the usual triple-point constants the model has $E(\infty )\simeq 2.5\times 10^{11}\mathrm{Pa}$, so it never reaches the proposed inventory limit. This is a counterexample inside the stipulated equations, not a physically credible high-pressure ideal gas ocean. Actual liquid water also ceases to exist above its critical point, before such an extrapolation would make sense.

Physically, the imposed finite-pressure-thickness layer can hold no more than $\Delta p/g$ of water vapor per unit area. As $x\to 1$, it expels all its hydrogen into the dry upper layer. A still larger surface saturation pressure cannot coexist with this geometry. The layer must deepen, the temperature profile must change, or the assumed phase equilibrium must fail. If $P_d<\Delta p$, the additional restriction $P>\Delta p$ can also exclude the low-temperature end of the branch.

At the interface the pressure is continuous. Immediately below it the temperature is $T_s$ and mean molar mass is $M_d(1+ax)$; immediately above it the dry hydrogen has temperature $T_1$. The ideal gas law therefore givesAn infinitesimal interchange is unstable to convection when the lower fluid is less dense than the upper fluid. HenceEquality gives neutral interfacial buoyancy; reversing the inequality gives a stable density jump. This is an interfacial criterion, not a claim that the whole dry upper atmosphere is stable. Its vertical temperature profile has not been specified. Also, $\mathrm{RH}(p_1)=p_1/P<1$, so infinitesimal displacement at the interface initially involves no condensation or release of latent heat.

There are two competing effects. Increasing $T_s$ directly reduces the lower density, favouring convection. Increasing the mole fraction of heavy water vapor increases its mean molar mass, opposing convection. With $x\ll 1$, the criterion is approximately $T_s>T_1$. With $x\to 1$, it tends to $T_s>ϵT_1\simeq 9T_1$: substantial surface heating need not overcome the compositional stabilization in hydrogen.

For a specified $P_d,\Delta p,T_1$, solve $T_s/T_1-1-ax(T_s)=0$ on the admissible branch to locate the transitions. Its derivative is $1/T_1-a{x}^{\prime }(T_s)$, so there is no universal monotone trend: rapid growth of the saturation pressure can suppress an instability that was present at lower $T_s$. Without numerical parameter values, neither the number of crossings nor a unique interval of unstable temperatures can be asserted.

The cloud absorbs $e\sigma T_g^4$ from below and emits $e\sigma T_c^4$ from each face. It absorbs no stellar radiation. Its net radiative heating per unit area is thereforeFor $e>0$, a cloud with $q^4<1/2$ warms in the absence of other energy transfers; one with $q^4>1/2$ cools. For the thin pressure thickness $\delta p$, its initial temperature tendency is approximately $g{\mathcal{R}}_c/(c_p\delta p)$ if its sensible heat capacity dominates. Latent heat and changing condensate modify that conversion from heating to temperature change, but not the sign of the radiative input.

With $T_g$ fixed, pure radiative equilibrium requiresThis equals the ambient dry-adiabatic process temperature precisely whenFor ordinary dry air, $\kappa \simeq 2/7$, so the ratio is about $0.545$. If $e=0$, the cloud exchanges no thermal radiation and its radiation budget does not select any equilibrium temperature.

At fixed ground temperature, warming the cloud increases its downward radiant flux $e\sigma T_c^4$, giving a positive contribution to the ground's radiative heating. That is the immediate radiative effect. A self-consistent final answer must also remove the convective heat transfer used to pin the original cloud to the dry adiabatic process.

For pure radiative equilibrium of both bodies, with no net nonradiative energy exchange, their two balances areFor $e>0$ this givesThe $e\to 0^{+}$ limit gives the selected limiting cloud temperature; at exactly $e=0$, only $T_g=T_0$ is fixed radiatively. The maximum greenhouse effect at fixed $Q$ occurs at $e=1$:The pinned dry-adiabatic process case instead has maximum factor $1/q=(p_s/p_c{)}^{\kappa }$. It exceeds $2^{1/4}$ for a high, cold cloud with $q<2^{-1/4}$, equals it at the level found in part b, and is smaller for a low, warm cloud.

In particular, for $q^4<1/2$ the final purely radiative ground is cooler than the pinned solution, even though a warmer cloud increases downward radiation at fixed $T_g$. To see why these statements are compatible, let $H$ be the net upward nonradiative energy transfer from the ground to the cloud in the pinned state. Its conservation of energy equation requiresFor the cold cloud, $H<0$: maintaining it on the assumed dry adiabatic process requires removing its net radiative heating and returning energy elsewhere. Setting $H=0$ changes that energy redistribution. It is therefore incorrect to infer the final ground temperature by changing only the downward radiant flux. Indeed the final cloud is warmer than the pinned one for $q^4<1/2$ and $0<e<1$, while the final ground is cooler. At $e=1$ the cloud temperature is $T_0$ in both whole-column equilibria.

Fixed cloud properties are an artificial constraint. Warming a saturated cloud increases the saturation pressure, tending to evaporate its cloud droplets; the released water vapor, changing latent heat budget and altered buoyancy can change its altitude and emissivity. Cooling favours condensation and potentially atmospheric precipitation. A warm perturbation that survives evaporation can rise by convection, but neither a unique vertical displacement nor survival of the cloud follows without moisture supply and an environmental temperature profile. Finally, pure radiative equilibrium of the two bodies does not require the intervening transparent gas to remain on the original dry adiabatic process; imposing that extra constraint is generally inconsistent except at the matching level in part b.

The dew-point depression is $6\mathrm{K}$, soKeeping $T_{D,s}=276.15\mathrm{K}$, $H\simeq 158.5\mathrm{km}$ and $c_p=1004\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$ gives approximately $747\mathrm{m}$, consistent with the rounded estimate.

The real cloud base can differ because entrainment changes the parcel's atmospheric mixing ratio and temperature, surface humidity varies, and the environmental pressure profile need not match the assumed parcel profile. Finite cloud-droplet activation requires some atmospheric supersaturation, so the first visible cloud droplets need not occur exactly at the thermodynamic lifting condensation level. A parcel must also actually be lifted that far; environmental static stability can prevent this. The constant-latent heat and first-order Taylor expansion approximations produce smaller numerical differences.

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.

For the variable-latent-heat saturation curve of a pure saturated water atmosphere, the total pressure equals the saturation pressure. The exact Clausius-Clapeyron relation along the phase coexistence curve iswhere $v_v,v_c$ are the specific volumes of the vapor and condensed thermodynamic phase. Neglecting $v_c$ and using the ideal gas law, $v_v=R_{\mathrm{wv}}T/p$, givesThe triple point provides the integration constant. Below it, the prescribed latent heat is constant, soAbove it, write $K=L_{\ast }/[R_{\mathrm{wv}}(T_c-T_{\ast })]$. Integrating $K(T_c/T^2-1/T)$ givesThese are the requested implicit temperature profiles $T(p)$. Both branches meet at $(T_{\ast },p_{\ast })$, and each is monotone before the limiting liquid-gas critical point. The lower branch describes equilibrium with ice, rather than a stable liquid surface. A stable liquid ocean surface belongs to the upper branch; the colder upper atmosphere may use the lower branch. Hydrostatic balance is needed to convert this relation into a height profile, but not to obtain $T(p)$.

Let $Y=-\ln (p/p_{\ast })$. For the temperature-dependent latent heat,Keeping the latent heat of vaporization fixed changes only the upper branch toBoth curves fall from $+\infty$ as $T\to 0^{+}$, pass through $Y=0$ at the triple point, and are convex. At the triple point there is a corner because the latent heat changes from $L_s$ to $L_{\ast }$. Their right-hand slopes there agree. Above the triple point,The variable-latent heat curve lies above the fixed-latent heat curve: its downward slope has smaller magnitude at every $T>T_{\ast }$. It reaches a horizontal tangent at $T_c$, whereas the fixed-latent heat curve is still descending. Both limiting ordinates are finite; this model gives $Y_{\mathrm{variable}}(T_c)\simeq -7.30$ and $Y_{\mathrm{fixed}}(T_c)\simeq -11.40$. These are extrapolations of the stipulated model, not accurate water liquid-gas critical point data. Extending the variable branch beyond $T_c$ would give negative latent heat and is physically inadmissible.