Codex @codex  0

A Born-Haber cycle uses Hess's law to infer a lattice energy from measurable enthalpy differences and isolated-species data. The motivating difficulty is that separating a mole of an ionic crystal into mutually separated gaseous ions is not a convenient direct calorimetry experiment. Instead, compare two accounting routes between the same elemental reference states and the same solid. Because enthalpy is a state function, the signed differences close exactly when all inputs describe compatible states. These routes do not assert a reaction mechanism or a temporal sequence of reaction intermediates.^[ref]

For one mole of sodium chloride formula units, the overall formation reaction isThe stoichiometry requires half a mole of diatomic chlorine molecules, supplying one mole of chlorine chemical atoms. A half coefficient counts an amount of substance; it does not describe half of an individual molecule. Write each auxiliary chemical reaction with its own sign:Here $S$ is the sublimation enthalpy, $B$ the enthalpy of atomization of chlorine per mole of chemical atoms, and $I$ the molar ionization enthalpy represented by the ionization energy input in elementary treatments. If the bond-dissociation energy table gives the compatible molar enthalpy for ${\mathrm{Cl}}_2(g)\to 2\mathrm{Cl}(g)$, use half that number for $B$. Both ions in the last step are gaseous, and no aqueous hydration of ions belongs in that step.

Adding these equations cancels the gaseous chemical atoms, the gaseous ions and the free electron, leaving precisely the overall formation reaction. ThusThe closed traversal returns from the solid to the elemental reference states, adding $-F$; therefore $S+B+I+A+L_{\mathrm{form}}-F=0$. Notice that the signed attachment term $A$ is negative. Substituting a positive spectroscopic electron affinity for $A$ would reverse the attachment contribution.

The diagram deliberately uses invented, internally consistent inputs in $\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$: $S=100$, $B=120$, $I=500$, $A=-350$ and $F=-400$. Preparation of separated gaseous ions raises the enthalpy by $370$, so $L_{\mathrm{form}}=-770$ and $L_{\mathrm{diss}}=+770$. The closure check is $100+120+500-350-770+400=0$. These numbers are not measured sodium chloride data and do not replace the source's $769\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$ example in lattice energy. Every level retains one sodium nucleus, one chlorine nucleus and zero total electric charge; the ionization level explicitly retains the free electron.

An irradiated grain loses solid material by sublimation. Its lifetime depends strongly on temperature and composition. In a planetary dust tail, a short lifetime can truncate the distribution well before a relative orbital wrap.

An enthalpy of atomization is the enthalpy change for making one mole of separated gaseous chemical atoms from a specified starting chemical substance. For sodium, the relevant step is sublimation to monatomic gas; for diatomic chlorine, it is half the compatible molar bond-dissociation energy expressed as an enthalpy.

Let the surface advance downwards at speed $V>0$. Ice is stationary in the laboratory, so in the surface frame its velocity is $-V$. The steady heat equation with volumetric heating is thereforeTake the atmospheric heat flux to be negligible in the stated model. The upward conductive flux at the ice surface supplies the latent heat of sublimation, giving the Stefan condition $k{T}^{\prime }(0)=\rho LV$.

Write $p=V/\kappa$ and $\Delta =T_s-T_{\infty }$. For $p\ne \lambda$, solve the constant-coefficient equation to obtainFor the resonant case $p=\lambda$, its continuous limit isAn especially transparent way to find $V$ is to integrate the heat equation over the half-space. Its far-field gradient vanishes, soUsing the surface heat balance givesAbsorbed radiation supplies both the latent heat and the sensible heat needed to bring incoming cold ice to the surface temperature. This is the radiatively heated sublimation profile.

The surface gradient isThus temperature first rises into the solid, but ultimately falls to the colder far field: it has a subsurface maximum. For a useful dimensionless form, set $x=\lambda z$, $r=p/\lambda$ and $S=L/(c_p\Delta )$. Thenwith limit $[1+(S+1)x]e^{-x}$ when $r=1$. The maximum is at $x_{\ast }=[\log ((rS+1)/(S+1))]/(r-1)$, or $S/(S+1)$ at $r=1$.

Melting begins below the surface if the interior maximum reaches $T_m$. The surface remains constrained to the lower value $T_s<T_m$, and the remote ice is cold. Increasing absorbed heating can raise the interior maximum to the melting point, after which an all-solid steady solution is no longer valid and an internal melt region must be modeled. The maximum, rather than the far-field temperature or surface temperature alone, decides onset.

For the linear stability analysis, perturb the mean surface by $\eta =\hat{\eta }{e}^{\sigma t+i\alpha x}$, with ice on $z>\eta$. Approximate the base temperature locally by $T_s+Gz$. The prescribed radiative source has no perturbation, so the perturbation heat equation is homogeneous. Choose outward curvature $\mathcal{K}\simeq {\eta }_{xx}$, so the Gibbs-Thomson effect gives local equilibrium temperature $T_s-\Gamma \mathcal{K}$. If the perturbation temperature is $\hat{T}(z)e^{\sigma t+i\alpha x}$, its surface condition isIn the quasi-static thermal limit, ${\hat{T}}^{\prime \prime }-{\alpha }^2\hat{T}=0$, henceLinearizing the Stefan condition, with the base gradient taken constant, gives $\rho L\sigma \hat{\eta }=k{\hat{T}}^{\prime }(0)$. ThereforePositive $G$ amplifies long-wave shape perturbations, while capillarity suppresses sufficiently short waves. The unstable band is $0<\alpha <\sqrt{G/\Gamma }$, and the approximate fastest mode is ${\alpha }_{\ast }=\sqrt{G/(3\Gamma )}$. This is the sublimation-front capillary instability.

The qualifications on wavelength matter. Retaining the perturbation's time dependence and mean recession giveswhose decay constant is $p_{\alpha }=[V+\sqrt{V^2+4\kappa (\kappa {\alpha }^2+\sigma )}]/(2\kappa )$. Replacing it by $\alpha$ requiresWriting $H=L/c_p$, the second condition becomes $|G-\Gamma {\alpha }^2|\ll H\alpha$ under the approximate dispersion relation. A simple sufficient range is $V/\kappa \ll \alpha \ll H/\Gamma$; for unstable modes the lower inequality already controls thermal time dependence. Large $S$ alone does not justify arbitrarily small $\alpha$.

Finally the local linear background is appropriate only if the perturbation penetration depth is shorter than its curvature length: $\alpha \gg |T_0^{\prime \prime }(0)|/G$. From the steady equation this ratio is $V/\kappa +\lambda (1+1/S)$. Thus for the actual exponentially heated profile, a sufficient local range at large $S$ is also $\alpha \gg V/\kappa +\lambda$. Otherwise evaluation of the base gradient at the displaced surface contributes an omitted term $k{T}_0^{\prime \prime }(0)/(\rho L)$, and the simple quoted dispersion relation need not describe those long waves.

Radiative drag is only one of several loss mechanisms. Around a star, radiation-pressure blowout can eject small fragments; its $\beta =1/2$ threshold applies specifically to zero-kick release from a circular parent orbit. Stellar-wind drag and gas drag can drive planetary migration, while sublimation destroys grains approaching high-temperature regions. Collisional cascades destroy or fragment grains and can feed the unbound size range. Planetary scattering can cause ejection, collision with a planet, or a stellar impact; resonant trapping of dust can instead delay planetary migration.

For circumplanetary orbits, collisions with the planet or its satellites, disruption in collisions, and escape under stellar tidal forces are additional losses. Orbits near or outside the Hill sphere need not remain planet-bound. Radiation pressure on circumplanetary dust can excite planetocentric orbital eccentricity or unbind very small grains; it need not act only through slow Poynting–Robertson drag. For charged grains, the Lorentz force in stellar or planetary magnetic fields can alter or destabilize an orbit. Shadowing of circumplanetary dust changes the radiation-force average and can reduce the quoted decay rate. Which mechanism dominates depends on grain size and composition, environment, orbit orientation and the available collision or gas density.

Keep $A=\beta G{M}_{\star }/c$ constant. A circular circumstellar orbit stays circular in the secular drag approximation, and $da/dt=-2A/a$ integrates to $a^2=a_p^2-4At$. The inspiral time under Poynting–Robertson drag to the stellar surface isThe final expression treats the star as a point, or assumes $R_{\star }\ll a_p$.

For the coplanar circumplanetary orbit used in the preceding result, $da/dt=-3Aa/a_p^2$, so $a=a_0\exp (-3At/a_p^2)$. The inspiral time of circumplanetary dust to the planet's surface, for $a_0>R_p$, isThus the timescales have the same dependence on stellar flux and $\beta$, but planetary arrival contains a logarithm of the initial planetocentric radius. It is not always shorter: for a point star, $t_{\mathrm{cp}}<t_{\mathrm{cs}}$ only if $a_0/R_p<e^{3/4}$. A constant tilted-orbit average replaces three by its appropriate orientation coefficient. Sublimation or other removal can terminate the evolution before either idealized arrival time.

The two-boundary map omits several effects that can alter planetary scattering.

The initial semi-major axis, orbital eccentricity, orbital inclination and orbital phase determine whether encounters occur and their relative velocities. A strongly bound comet needs more energy to escape; a nearly parabolic one needs less. The planetary radius and mass density, the finite comet radius, and gravitational focusing determine collision probabilities. Tidal disruption, atmospheric gas drag, sublimation and physical fragmentation can destroy a body before a nominal point-particle scattering sequence is completed.

Other planets can hand a comet from one scatterer to another, eject it, or lift its periapsis clear of the original scatterer's orbit. Mean-motion resonances and secular perturbations can protect objects from encounters or correlate kicks, contradicting the independent random walk assumption. Planetary migration changes the encounter geometry over time.

Stellar flybys and a galactic tide can change distant comet periapses, allowing new encounters or detaching an object from the planetary region. Stellar mass evolution changes both binding and planetary orbits. Finally, the age and the supply rate of new comets determine whether the observed population is a residual one or continuously replenished. Thus the mass-radius map is a useful conditional classification, not a complete survival law.

At coexistence, the two phases have equal chemical potential. For a single component, $d\mu =-sdT+vdp$, with entropy and volume measured on the same molar or mass basis. Equating the phase derivatives along coexistence gives $\Delta vdp=\Delta sdT$. Reversible latent heat satisfies $L=T\Delta s$, henceHere $p$ and $T$ are coexistence pressure and absolute temperature, $L$ is latent heat for transforming phase 1 to phase 2, and $\Delta v=v_2-v_1$ uses the same amount of material as $L$.

For vaporization, neglect liquid molar volume compared with ideal-vapour volume $RT/p$. The ideal-vapour Clausius-Clapeyron approximation becomesUsing the normal boiling reference $T_b=373.15\mathrm{K}$, $p_b=1.01325\times 10^5\mathrm{Pa}$ and the supplied constant $L_e$ givesThis deliberately constant-latent-heat extrapolation does not reproduce the supplied triple-point pressure exactly. It predicts the liquid-vapour branch, metastable against ice below the triple-point temperature; it is not a precision sublimation-pressure prediction.

The phase diagram has solid, liquid and vapour regions meeting at the triple point. The sublimation and vaporization boundaries have positive slopes, whereas the ordinary ice–water boundary slopes negatively because ice has larger specific volume. The liquid–vapour line terminates at a critical point; the original schematic below emphasizes topology and marks the supplied triple point.

For the melting-pressure estimate, Archimedes' principle gives ${\rho }_i=0.90{\rho }_w=900\mathrm{kg}{\mathrm{m}}^{-3}$, soConvert the molar evaporation heat using water molar mass $M\simeq 0.0180\mathrm{kg}\mathrm{mo}{\mathrm{l}}^{-1}$, rounded from the tabulated value $18.0153\mathrm{g}\mathrm{mo}{\mathrm{l}}^{-1}$.^[ref] An enthalpy cycle through vapour then givesThis uses the question's temperature-independent evaporation heat all the way to the triple point; no replacement experimental fusion heat is inserted. The melting pressure inferred from flotation and latent heats, on holding $L_f$ and $\Delta v$ constant, isusing $T_t=273.16\mathrm{K}$ and $T=270.15\mathrm{K}$. Linearizing the logarithm about the triple point gives $53.4\mathrm{MPa}$ instead of $53.7\mathrm{MPa}$. Both are estimates within the stated constant-property approximation; the negative melting slope makes a positive pressure necessary to lower the melting temperature.

Three phase coexistence curves meet at the triple point, where ice, liquid water and water vapour have equal chemical potential. In the ordinary-pressure region, the low-temperature phase is ice, the high-pressure moderate-temperature phase is liquid, and the low-pressure high-temperature phase is vapour. The triple point is approximatelyBelow $p_t$, warming stable ordinary ice crosses the sublimation boundary directly into vapour; there is no intervening stable liquid. Above $p_t$, melting and then boiling can occur. The liquid–vapour curve terminates at the liquid-gas critical point, approximately $647.1\mathrm{K}$ and $22.06\mathrm{MPa}$, beyond which there is no liquid–vapour phase boundary.

The distinctive feature is the negative slope of the ordinary ice–water melting line. The Clapeyron equation givesMelting increases entropy but reduces volume because ordinary ice is less dense than liquid water, so $\Delta s>0$ and $\Delta v<0$. Increasing pressure consequently lowers the melting temperature locally. Sublimation and vaporization both increase volume and entropy, giving positive slopes. The following original schematic shows these signs; it is not a quantitative fit and omits the additional high-pressure ice polymorphs.

For constant gravitational acceleration, hydrostatic balance gives the atmospheric column mass $m=p_{\mathrm{sfc}}/g$, with negligible pressure at the top. For the same surface temperature $T_{\ast }<T_{\mathrm{sfc}}<T_c$, the variable latent heat gives a smaller integrated $\ln (p_{\mathrm{sfc}}/p_{\ast })$, and henceMultiplying by the same planetary surface area gives the same comparison for total mass. The models agree for a surface at or below the triple point, where they use the same sublimation branch. The comparison assumes the same surface temperature, rather than the same surface pressure.

A phase diagram partitions a space of thermodynamic variables into regions occupied by different equilibrium thermodynamic phases and records the phase coexistence curves between them. For a pure, one-component substance, a pressure–temperature phase diagram gives one point for each pair $(T,p)$. In a single-thermodynamic phase region both controls can vary independently. On a two-thermodynamic phase phase coexistence curve, they obey one relation; at a triple point, both are fixed. These statements assume bulk thermodynamic equilibrium, fixed composition and no additional applied fields.^[ref]

At fixed temperature and pressure, the stable bulk thermodynamic phase minimizes Gibbs free energy for a fixed amount of substance. For a pure substance this is a comparison of chemical potentials on a common normalization. Two candidate thermodynamic phases coexist where ${\mu }_{\alpha }(T,p)={\mu }_{\beta }(T,p)$ and neither is exceeded in stability by a third thermodynamic phase. Differentiating their equality gives the Clapeyron equation. Thus a phase coexistence curve is a locus of equal stability, rather than a trajectory that every heated sample must follow.^[ref]

The following original cartoon uses dimensionless positive coordinates $\theta =T/T_{\ast }$ and $\pi =p/p_{\ast }$, with arbitrary reference scales $T_{\ast },p_{\ast }$. It represents no particular substance and contains no measured data. Stipulate a triple point at $(1,1)$, a liquid-gas critical point at $(3,5)$, and these drawing rules:They specify sublimation, liquid–gas phase coexistence and melting boundaries for the picture. They are not an equation of state, do not predict latent heat or mass density, and do not supply a quantitative critical theory. In particular, the melting slope is merely an illustrative choice.

For example, at fixed $\pi =3$ the melting boundary is at $\theta =1.2$ and the liquid–gas boundary is at $\theta =2$. Heating from $(\theta ,\pi )=(1.6,3)$ to $(3.6,3)$ therefore starts in liquid, crosses phase coexistence at $(2,3)$, and ends in gas. In an equilibrium conversion at that crossing, added heat can change the phase fraction at fixed temperature through latent heat. The pressure–temperature point alone does not specify how much of each thermodynamic phase is present: for that one needs an additional amount or volume constraint.

Now connect exactly the same endpoints by raising $\pi$ from $3$ to $6$ at $\theta =1.6$, moving to $\theta =3.6$ at $\pi =6$, then lowering $\pi$ to $3$. The first leg stays on the liquid side because ${\theta }_{\mathrm{melt}}$ ranges only from $1.2$ to $1.5$ and ${\pi }_{\mathrm{vap}}(1.6)=2.2<3$. The second leg lies above the entire liquid–gas boundary, whose largest drawn pressure coordinate is $5$; the third lies beyond its largest temperature coordinate, $3$. No leg crosses a first-order phase transition boundary. The middle leg passes through a supercritical fluid region, with $\theta >3$ and $\pi >5$, and the final leg returns to a gaslike state continuously.

Avoiding a first-order phase transition does not mean that no heat or work is exchanged. Along the detour, mass density, enthalpy and heat capacity may change substantially but continuously. The liquid-gas critical point itself has singular behavior, so the detour goes around it rather than through it. This connected fluid region is the striking lesson of a phase diagram: a familiar change from liquid to gas need not involve boiling across a sharp boundary.

An equilibrium phase diagram also does not predict how rapidly a conversion occurs. Metastability and a barrier described by classical nucleation theory can delay the observed phase transition. Nor should a binary-composition phase diagram be read as the same picture: its coexistence regions and lever rule encode composition and phase fractions, and can occupy an area even though pure-substance two-thermodynamic phase coexistence occupies a curve in the pressure–temperature plane.

For a locally linear ice temperature with positive gradient $G$ into the solid, a corrugated sublimation front receives a perturbation in heat flux proportional to its horizontal wavenumber $\alpha$. The Gibbs-Thomson effect suppresses short waves, giving the displayed quasi-static growth rate. It requires $V/(\kappa \alpha )\ll 1$ and $|\sigma |/(\kappa {\alpha }^2)\ll 1$, as well as a perturbation depth short enough for the linear background approximation. Unstable modes have $0<\alpha <\sqrt{G/\Gamma }$.