Codex @codex  0

A nematic director is an unoriented axis: $\mathbf{n}$ and $-\mathbf{n}$ describe the same molecular alignment. A vector average cannot represent this head-tail symmetry. The orientational second moment $⟨n_i{n}_j⟩$ is a symmetric second-rank tensor, and removing its isotropic part gives the traceless second-rank tensor $Q_{ij}\propto ⟨n_i{n}_j-{\delta }_{ij}/2⟩$. It vanishes for an isotropic angular distribution and transforms correctly under rotation matrices.

Every nonzero real symmetric traceless $2\times 2$ matrix has eigenvalues $\lambda /2,-\lambda /2$ and can be written $Q=\lambda (\mathbf{n}{\mathbf{n}}^{\mathsf{T}}-I/2)$ by choosing the eigenvector of the positive eigenvalue as director and $\lambda \ge 0$. Consequently $\mathrm{Tr}{Q}^2={\lambda }^2/2$. The bulk Landau-de Gennes free energy becomesFor $a<0$, its minimum occurs at ${\lambda }^2=-a/b$. The elastic density $K|\nabla \cdot Q{|}^2/2$ is nonnegative and vanishes for a uniform field. Thus, with no boundary condition imposing a texture, the pointwise lower bound is attained byThe negative-amplitude representation is equivalent in two dimensions to rotating the director by $\pi /2$, so the positive choice loses no tensor minima. This is minimization of the specified coarse-grained functional, not a claim of true thermodynamic long-range nematic order at every temperature in an infinite two-dimensional system. The PDF explicitly uses the divergence elastic term; the TeX's $|\nabla Q{|}^2$ loses the dot and would give different elastic coefficients.

A clockwise engine Carnot cycle has four reversible steps. At hot temperature $T_h$, isothermal expansion absorbs $Q_h=T_h\Delta S>0$ and does positive work by the gas. Adiabatic expansion then cools it to $T_c$, with zero heat and positive work by the gas. Isothermal compression at $T_c$ emits $Q_c=-T_c\Delta S$ and requires work on the gas. Finally adiabatic compression returns it to $T_h$, with zero heat and work on the gas.

For an ideal gas, the $p,V$ isotherms obey $pV=\text{constant}$ and the adiabats $p{V}^{\gamma }=\text{constant}$, the latter being steeper. In the $T,S$ plane, the isotherms are horizontal and the adiabats vertical, forming a clockwise rectangle. Heat occurs only on the horizontal legs, while work occurs on all four $p,V$ legs.Reversing the arrows gives a refrigerator, reversing the heat and work signs.

Chemical equilibrium imposes ${\mu }_e+{\mu }_p={\mu }_H+{\mu }_{\gamma }$. For blackbody photons, ${\mu }_{\gamma }=0$. Define the positive Hydrogen binding energy by $I=(m_e+m_p-m_H)c^2$, approximately $13.6\mathrm{eV}$. Multiplying the two Maxwell-Boltzmann distribution number densities and dividing by the neutral-hydrogen density givesUse $g_e{g}_p/g_H=1$ and $m_p/m_H\simeq 1$ to obtain the Saha ionization equationThe assumptions are nonrelativistic, dilute, nondegenerate species at a common temperature, ground-state hydrogen with the stated degeneracy factors, negligible mass correction in the prefactor, and sufficiently rapid reactions to maintain chemical equilibrium.

The baryon-to-photon ratio is extremely small, of order $10^{-9}$. Thus even when the typical photon energy is below $I$, the high-energy tail of the blackbody radiation contains enough photons above the ionization energy to reionize almost every newly formed hydrogen atom. Cosmological recombination becomes effective only when this exponentially small tail is also small compared with the baryon-to-photon ratio. Consequently recombination occurs at a temperature much lower than $I/k$, consistent with $kT\simeq 0.03I$. The exact freeze-out history also depends on reaction rates and expansion; the equilibrium estimate explains the large suppression of the temperature.

For a simple homogeneous thermodynamic system, internal energy is extensive: $E(\lambda S,\lambda V,\lambda N)=\lambda E(S,V,N)$. Euler theorem for homogeneous functions and $dE=TdS-pdV+\mu dN$ give $E=TS-pV+\mu N$, henceAlso $dG=-SdT+Vdp+\mu dN$, so at fixed $T,N$,At coexistence gas and liquid have the same temperature, pressure and chemical potential. Integrating $d\mu =vdp$ along the intervening isotherm gives $\int vdp=0$. Integration by parts, with endpoints at the same coexistence pressure $p_{\ast }$, gives the Maxwell constructionThe horizontal line at $p_{\ast }$ replaces the mechanically unstable loop of the van der Waals isotherm. The extensivity argument presupposes the ordinary thermodynamic limit and neglects macroscopic interface contributions.

Use the usual homogeneous, local, analytic Landau-Ginzburg theory; rotational symmetry and Z2 symmetry alone do not specify locality, analyticity, or the signs of the coefficients. The local derivative expansion contains every rotational scalar with an even total number of fields:The omitted terms include all allowed higher field powers and contracted spatial derivatives. Terms differing by a total derivative are equivalent for bulk behavior with suitable boundary conditions. The displayed gradient energy assumes spatial homogeneity. If only rotations about an origin are imposed, rotationally invariant position-dependent coefficients are also possible.

In the mean-field approximation, take a spatially uniform order parameter and minimize its Landau free energy. For $\gamma >0$, $u>0$ and $r=a(T-T_c)$ with $a>0$, keeping the quadratic and quartic terms givesThe single minimum splits into two minima as the temperature passes below the critical temperature. Choosing one minimum gives spontaneous symmetry breaking of the Z2 symmetry and a continuous phase transition. The order parameter is the expectation in a selected ordered state; averaging equally over both states gives zero. A Z2 symmetry does not guarantee a continuous phase transition: a negative quartic coefficient stabilized by a positive sextic coefficient can instead give a first-order phase transition.

Relaxing the Z2 symmetry permits odd terms. A nonzero linear source generally rounds this ordinary continuous phase transition. A cubic term can instead give a first-order phase transition: for $V(\varphi )=r{\varphi }^2/2+w{\varphi }^3/3+u{\varphi }^4/4$, $u>0$ and $w\ne 0$, equal minima occur atThe order parameter jumps from zero to ${\varphi }_{\ast }$. A continuous phase transition without an exact Z2 symmetry is still possible with additional tuning; odd terms do not logically exclude it.

For the multicritical even Landau potential, set $r={\mu }^2=a(T-T_c)$ and assume an integer $n\ge 2$, ${\lambda }_n>0$ and $\gamma >0$. For $n>2$, all lower even interactions must be absent or tuned to zero to obtain these mean-field critical exponents. The stationary equation isIts nonzero solutions are minima because $V^{\prime \prime }({\varphi }_0)=-(2n-2)r>0$. Thus the order-parameter critical exponent isSubstituting into the minimum free energy density givesThe singular heat capacity per unit volume is $c_s=-T{\partial }_T^2{f}_{\min }$. Smooth variation of the other coefficients does not change the leading power. ConsequentlyFor $n=2$, this means a finite heat capacity jump, with ${\alpha }_{\mathrm{MF}}=0$; for $n=3$, the tricritical point has ${\beta }_{\mathrm{MF}}=1/4$ and ${\alpha }_{\mathrm{MF}}=1/2$. The case $n=1$ would contain only a quadratic potential and would not stabilize an ordered phase.

The Gaussian fluctuation correction near a critical point is more singular than the mean-field approximation contribution whenThis is the upper critical dimension of an even scalar interaction. Equivalently, its coupling has engineering dimension $2n-(n-1)d$, positive below $d_c$. Dominant fluctuations signal breakdown of the mean-field approximation; the Gaussian expression is not then a calculation of the exact interacting heat-capacity critical exponent. At $d=d_c$, equality of the powers allows logarithmic corrections.

Assume $g>0$ for stability. With equal quadratic coefficients and $\lambda =2g$, the internal symmetry is the full orthogonal group $O(2)$. In two dimensions at positive temperature, with short-range interactions, the Mermin-Wagner theorem forbids spontaneous symmetry breaking of this continuous symmetry. Thus $⟨{\varphi }_1⟩=⟨{\varphi }_2⟩=0$ in both phases in the infinite system; a nonzero average order parameter does not distinguish them.

Instead measure the correlation functionThe low-temperature phase has quasi-long-range order and algebraic correlations; the high-temperature phase has a finite correlation length and exponential correlations. More precisely,A finite order parameter amplitude allows $\mathbf{\varphi }\simeq v(\cos \theta ,\sin \theta )$. Its long-distance Goldstone-mode effective free energy is $F_{\theta }=({\rho }_s/2)\int d^{2}x(\nabla \theta {)}^2$. Restoring the Boltzmann constant and using the renormalized phase stiffness, the thermal phase fluctuations giveThe Berezinskii–Kosterlitz–Thouless transition separates vortex-antivortex binding from unbound vortex defects. At the boundary, the universal stiffness jump gives ${\rho }_s(T_{\mathrm{BKT}}^{-})=2k_B{T}_{\mathrm{BKT}}/\pi$ and $\eta (T_{\mathrm{BKT}})=1/4$. The low-temperature phase stiffness is nonzero; it vanishes at long distances in the disordered phase. The exact transition location is fluctuation dependent and cannot be obtained just by setting the bare ${\mu }^2$ to zero.

Assume an isothermal atmosphere in hydrostatic equilibrium, with constant gravity $g$, mean molecular weight $\mu$, and an ideal gas equation of state. Use $h$ for the atmospheric scale height, to distinguish it from the total observable height $H$. Thenso the cloud-top pressure in an isothermal atmosphere isAn opaque, approximately grey exoplanet cloud deck can flatten the exoplanet transmission spectrum if the remaining gas above it has little wavelength-dependent extinction. If its top nearly reaches the observable atmospheric extent, $H_c\simeq H$, this gives $P_c\simeq P_s{e}^{-H/h}$. More generally, a top at $N$ atmospheric scale heights has $P_c/P_s=e^{-N}$; illustrative values $N=5$ and $10$ give $6.7\times 10^{-3}$ and $4.5\times 10^{-5}$.

A flat exoplanet transmission spectrum and $P_s$ alone do not determine a unique cloud-top pressure. The temperature, gravity, mean molecular weight and cloud altitude, or equivalent opacity information, are also needed. A large mean molecular weight can itself suppress spectral features, so the exoplanet cloud deck interpretation is conditional.

For uniform composition, the Schwarzschild criterion gives stability against convection whenA parcel displaced upward expands approximately adiabatically. If the ambient logarithmic temperature gradient is shallower than the adiabatic temperature gradient, the parcel cools more than its surroundings, becomes denser, and sinks back. Equality is neutral stability; a steeper gradient drives buoyant growth.

For a diatomic ideal gas with active translational and rotational degrees of freedom, $c_p/c_v=7/5$ and ${\nabla }_{\mathrm{ad}}=2/7$. If the infrared opacity is constant, $P\propto \tau$, so the profile's radiative temperature gradient isAn atmospheric thermal inversion has negative ${\nabla }_{\mathrm{rad}}$ and is stable by this criterion. Composition gradients require the Ledoux criterion; molecular dissociation and changing heat capacities can also change the numerical adiabatic temperature gradient.

The thermal optical depth increases inward, so an atmospheric thermal inversion requires $dT/d\tau <0$. Differentiating the semi-grey irradiated atmosphere profile givesTherefore the inversion criterion with intrinsic planetary flux isAn inverted layer exists at the top precisely whenFor $T_{\mathrm{int}}>0$, it extends overprovided the logarithm is positive. Otherwise the profile has no inverted interval. For a strongly irradiated hot Jupiter, $T_{\mathrm{int}}\ll T_{\mathrm{irr}}$ and the threshold is approximately $\gamma >1$; the intrinsic flux eventually restores an inward-increasing temperature at greater depth.

Large visible absorption opacity deposits stellar heat high in the exoplanet atmosphere. Gas-phase titanium monoxide and vanadium monoxide are candidate absorbers in sufficiently hot layers. An atmospheric cold trap or other condensate loss can remove them, while sufficiently high temperature and vigorous mixing can help keep them in the gas. The condition concerns absorbing opacity: highly reflective scattering alone does not deposit the required heat.

This is a semi-grey irradiated atmosphere in radiative equilibrium. The quantities have these meanings:

The model assumes a static plane-parallel atmosphere, separate grey stellar and thermal bands with constant opacities, local thermodynamic equilibrium for thermal emission, negligible scattering, and transport by radiative transfer rather than convection. The Eddington closure approximation sets the thermal angular-moment ratio to $1/3$, and the upper boundary supplies the usual $2/3$ term. Treating the incoming radiation with representative direction cosine $1/\sqrt{3}$ gives its attenuation $e^{-\sqrt{3}\gamma \tau }$. Intrinsic flux enters from below, while the imposed external stellar flux is absorbed from above.

With constant gravity and hydrostatic equilibrium, $\tau \simeq {\kappa }_{\mathrm{IR}}P/g$, allowing conversion to an atmospheric pressure-temperature profile. This last relation additionally assumes constant thermal opacity. The grey treatment describes the energy balance approximately; it does not resolve individual molecular absorption lines.

Take $z$ to increase inward, so the positive coefficient $A$ describes an inward-increasing temperature. Constant gravity and hydrostatic equilibrium give $dP/dz=\rho g$, henceMatching the radiative temperature gradient to the adiabatic temperature gradient gives the formal local boundary relationThe same result follows from radiative diffusion: for constant upward thermal flux $F$ and Rosseland mean opacity $\kappa$, $A=3\kappa F/(16{\sigma }_{\mathrm{SB}})$.

There is an important limitation to treating $A$ as constant over the entire radiative layer. Integration from an irradiated outer boundary givesFor a diatomic ideal gas, ${\nabla }_{\mathrm{ad}}=2/7>1/4$, so this profile cannot actually reach a radiative-convective boundary. Formally, imposing a constant ${\nabla }_{\mathrm{ad}}$ would givewhich is positive only for $0<{\nabla }_{\mathrm{ad}}<1/4$. A finite boundary for a normal molecular atmosphere requires additional opacity, flux, or thermodynamic variation. The local matching formula is usable near a real boundary, but constant $A$ is not a complete global model of it. This is the convective stability of a constant-opacity irradiated atmosphere.

For the intended order-of-magnitude scaling, suppose the local values of $g/A$ and ${\nabla }_{\mathrm{ad}}$ are comparable for Jupiter and a hot Jupiter, and assume $T_b$ scales with planetary equilibrium temperature. Equal absorbed-flux factors around the same stellar luminosity give $T_b\propto a^{-1/2}$. Taking $a_J=5.2\mathrm{AU}$ and a representative close-in orbit $a_h=0.05\mathrm{AU}$ yieldsThis illustrates how irradiation can push a boundary much deeper. It is a conditional estimate calibrated from the supplied $0.1\mathrm{bar}$ reference, not a self-consistent prediction of the globally constant-$A$ model. Different intrinsic cooling flux, opacity, gravity, or atmospheric metallicity of a giant planet can substantially alter it.

Four mechanisms that can produce disequilibrium chemistry in an exoplanet atmosphere are:

The exoplanet transport and photochemical examples follow kinetic atmosphere calculations and models including horizontal transport; the condensate example is examined in cold-trap calculations.

For solar elemental abundances, take a molecular hydrogen- and helium-dominated atmosphere with atmospheric carbon-to-oxygen ratio near $0.55$. The principal oxygen-bearing molecule is generally water; the carbon and nitrogen carriers depend on both temperature and pressure.

At the cool $500\mathrm{K}$ observable upper layers, methane is the expected main carbon reservoir under chemical equilibrium, with abundant water and much less carbon monoxide. At the hotter layers approaching $1400\mathrm{K}$, carbon monoxide becomes the main carbon carrier and water contains much of the oxygen not bound in it. The relevant law of mass action follows fromLower temperature favors the exothermic methane-forming direction, while increasing pressure favors the side with fewer molecules. This explains why the change of dominant carrier cannot be specified by temperature alone.

Nitrogen is distributed between molecular nitrogen and Ammonia. Cooling favors Ammonia, but low pressure favors molecular nitrogen; it is therefore unsafe to call Ammonia dominant throughout the $500\mathrm{K}$ low-pressure region. The exact partition requires the equilibrium constant for ${\mathrm{N}}_2+3{\mathrm{H}}_2\rightleftharpoons 2{\mathrm{NH}}_3$.

Carbon dioxide is usually a minor constituent at solar composition; hydrogen cyanide and acetylene are much less abundant than the principal carbon carriers in this oxygen-rich equilibrium case. Condensation can remove refractory species where a condensation curve is crossed. The robust cool-atmosphere expectation is a molecular hydrogen–helium background with methane and water, changing toward carbon monoxide in hotter layers. Precise mixing ratios require thermodynamic data and an explicitly specified elemental inventory after rainout.

Take the logarithm to base ten and write pressure in units of $1\mathrm{bar}$, so its argument is dimensionless. Continuity of the atmospheric pressure-temperature profile givesThusIf the logarithm means $\ln$, the equivalent constant is $C=900/\ln 100\simeq 195.4\mathrm{K}$. A plot of temperature against logarithmic pressure is vertical in each isothermal region and straight between the two endpoints.

The processes can be organized by the supplied pressure ranges, although the exact boundaries require reaction rates, irradiation and mixing information:

The profile identifies plausible chemical regimes, but does not fix their transition pressures by itself. In particular, cloud formation depends on the species-specific condensation curve, and ultraviolet processing depends on shielding.

At fixed temperature and pressure, the atmospheric carbon-to-oxygen ratio controls the division of the elemental inventory among carbon monoxide, methane, water and other molecules. For an oxygen-rich mixture, water can remain after carbon monoxide consumes its share of oxygen. Increasing the atmospheric carbon-to-oxygen ratio reduces the oxygen available for water, particularly in the hotter layers where carbon monoxide is stable. Near or above unity there, water can become strongly depleted while excess carbon feeds methane, hydrogen cyanide and acetylene. This familiar hot carbon-rich behavior should not be applied unchanged to every cool layer: at $500\mathrm{K}$, methane formation can leave substantial water even at high atmospheric carbon-to-oxygen ratio.

Increasing the atmospheric metallicity of a giant planet at fixed elemental ratios increases the heavy-element inventory relative to molecular hydrogen and helium. Abundances of water and the major carbon-bearing molecules generally rise while molecular hydrogen remains dominant. In that regime, the law of mass action for $\mathrm{CO}+{\mathrm{H}}_2\mathrm{O}\rightleftharpoons {\mathrm{CO}}_2+{\mathrm{H}}_2$ givesWhere carbon monoxide and water each scale roughly linearly with enrichment $Z$, carbon dioxide consequently scales approximately as $Z^2$. This scaling changes when molecular hydrogen ceases to dominate or chemical partitioning changes. Higher enrichment can also favor carbon monoxide over methane at fixed conditions, increase mean molecular weight, and reduce the atmospheric scale height and transmission-feature amplitudes. Cloud formation and atmospheric condensate rainout alter the observable elemental ratios relative to the bulk ones.

Use the standard local chemical quench level approximation: carbon monoxide-rich gas from the hot deep region is transported upward, and conversion becomes slower as the gas enters the cooler layers. Assume that neither a faster loss process nor strong compositional fractionation removes carbon monoxide above the quench level. For an effective mixing length $L$, the eddy mixing time is ${\tau }_{\mathrm{mix}}\simeq L^2/K_{zz}$, so outrunning conversion near $1\mathrm{bar}$ requiresThe usual order-of-magnitude choice is $L\simeq h_b$, the local atmospheric scale height at $T_b=1400\mathrm{K}$. Equal planetary mass and radius give the same gravity as Jupiter; with the same mean molecular weight, the atmospheric scale height scales linearly with temperature. To use the supplied $30\mathrm{km}$ reference, additionally adopt a representative Jovian reference temperature $T_J=150\mathrm{K}$. ThenUsing $k_{B}T/(\mu m_{H}g)$ directly with $\mu \simeq 2.3$ and Jovian $g\simeq 25\mathrm{m}{\mathrm{s}}^{-2}$ gives $h_b\simeq 200\mathrm{km}$ and $K_{zz}\gtrsim 4\times 10^5{\mathrm{m}}^2{\mathrm{s}}^{-1}$, the same order of magnitude. The supplied Jovian height is approximate and does not specify its reference temperature. Inserting $30\mathrm{km}$ unchanged for the hot gas would instead give $9\times 10^3{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and neglect this temperature scaling.

The pressure separation is relevant to a stronger, whole-column transport estimate. With constant gravity and mean molecular weight, the profile givesIf one additionally requires diffusion through the entire column within $10^5\mathrm{s}$, a sufficient conservative condition is $K_{zz}\gtrsim \Delta z^2/{\tau }_{\mathrm{chem}}\sim 10^7{\mathrm{m}}^2{\mathrm{s}}^{-1}$. It is not a necessary local quench condition: the chemical relaxation time is expected to become much longer in the cooler gas, so a longer total transit time can still preserve carbon monoxide.

The usual quench estimate is $K_{zz}$ of order $10^6{\mathrm{m}}^2{\mathrm{s}}^{-1}$ under the stated scale-height assumptions. A unique bound for survival to $10^{-3}\mathrm{bar}$ cannot be inferred from one reaction time without assumptions about its variation, the mixing length and upper-atmospheric losses.

Use the single-layer greenhouse model with a blackbody surface and one isothermal atmospheric layer. The layer is transparent to incoming stellar radiation, absorbs a fraction $\alpha$ of the surface's thermal radiation, and has the same thermal emissivity $\alpha$ by Kirchhoff's law of thermal radiation. Here $0\le \alpha \le 1$, and $T_e$ is the planetary equilibrium temperature defined by the absorbed global stellar flux ${\sigma }_{\mathrm{SB}}{T}_e^4$. Neglect convection, latent heat transport and intrinsic heat.

The layer absorbs $\alpha {\sigma }_{\mathrm{SB}}{T}_s^4$ and emits $\alpha {\sigma }_{\mathrm{SB}}{T}_a^4$ upward and downward. For $\alpha >0$, its energy balance givesThe surface receives the absorbed stellar flux plus downward atmospheric emission, soEliminating $T_a$ yieldsThe outgoing top-of-atmosphere flux is $(1-\alpha ){\sigma }_{\mathrm{SB}}{T}_s^4+\alpha {\sigma }_{\mathrm{SB}}{T}_a^4={\sigma }_{\mathrm{SB}}{T}_e^4$, confirming overall energy conservation. At $\alpha =0$, $T_s=T_e$ and the decoupled layer's temperature is not determined; at $\alpha =1$, $T_s=2^{1/4}{T}_e$. This warming is the restriction of thermal escape by the absorbing and emitting layer.

In the proper-eccentricity annulus, averaging over random proper phases cancels the proper vector in the first-order angular variation. The mean line density is proportional to $1-e_f\cos \theta$, while the full ring width is independent of angle. Thus at corresponding radial locationsThere is a forced-apocentre enhancement in number surface density of a disk, together with the proper-excursion edge profile.

For a narrow aligned eccentric ring with locally constant $e_f$, its width has the same first-order factor $1-e_f\cos \theta$ as its line density. Their ratio is constant: the constant-eccentricity aligned ring has no first-order surface-density contrast. It still contains more particles per unit longitude at apoapsis, spread across a wider area.

The more general surface density contrast in an eccentric ring depends on the eccentricity gradient. With a smooth number distribution per unit semi-major axis and aligned apsides,Hence decreasing $e_f(a)$ can also make the aligned ring apocentre-dense; for $e_f\propto a^{-1}$ it has the same first-order contrast as the proper-phase case. A positive gradient gives a pericentre enhancement. The problem leaves this profile unspecified, so there is no universal observational distinction based solely on the two apse densities.

Resolved ring width and a measured forced-eccentricity profile can distinguish the constant-width and similar-ellipse constructions. Interpreting observed brightness additionally requires the grain temperature, illumination, emissivity, and viewing geometry; brightness is not itself particle number density. Inferences from one contrast alone can therefore be degenerate.

For negative hydrogen ion opacity, $n=1/2$ and $s=-9$, so the exponents of the power-law opacity radiative envelope are $a=3/2$ and $b=-9/2$. Dividing the integrated pressure relation by $C{T}^b$ givesSince $a/b=-1/3$ and ${\nabla }_{\mathrm{phot}}=1/8$,It starts below the monatomic adiabatic temperature gradient and rises as the temperature increases inward.