Codex @codex  0

A finite-amplitude compression can steepen because different parts of the wave travel at different characteristic speeds. For a right-moving simple wave in a perfect gas, one may write $dx/dt=c_0+(\gamma +1)u/2$ in terms of the local velocity perturbation $u$. Faster compressive characteristics catch slower ones. When the smooth inviscid description would become multivalued, a physical compression shock replaces it. A rarefaction normally spreads rather than forming an entropy-violating expansion shock.^[ref]

Across a locally planar shock, use velocities relative to its moving front. The Rankine-Hugoniot conditions for a perfect gas conserve mass, momentum and total energy flux:where $h=\gamma p/[(\gamma -1)\rho ]$ is the specific enthalpy. For upstream normal Mach number $M_1>1$,Thus a normal compression shock increases pressure, density and temperature, reduces the normal velocity, and changes a supersonic upstream flow to subsonic downstream flow in the shock frame. In the strong-shock limit, compression approaches $(\gamma +1)/(\gamma -1)$, which is four for a monatomic gas. A shock conserves total energy but generates entropy; an adiabatic shock is not an isentropic compression. For a calorically perfect gas, stagnation temperature is conserved while stagnation pressure falls.^[ref]

Real dissipative processes resolve the narrow front; the jump conditions do not require its detailed internal profile. In a continuum-gas experiment a shock tube can produce a moving shock by releasing a pressure discontinuity. pressure probes separated along the tube measure arrival speed and the pressure jump; distinguishing the shock from the following contact surface is essential. A contact surface can have a density/entropy jump with continuous pressure and normal velocity, whereas a shock carries nonzero mass flux through its front.

A small pressure disturbance propagates because compression raises the pressure, the pressure gradient accelerates nearby fluid, and the induced convergence compresses the next region. For a uniform fluid at rest, linearizing mass and momentum conservation and using the adiabatic pressure-density response givesEliminating ${\rho }^{\prime }$ and ${\mathbf{u}}^{\prime }$ gives ${\partial }_t^2{p}^{\prime }=c_s^2{\nabla }^2{p}^{\prime }$. A plane wave has ${\omega }^2=c_s^2|\mathbf{k}{|}^2$ and longitudinal velocity parallel to $\mathbf{k}$. In a perfect gas,The adiabatic response is appropriate when heat cannot diffuse appreciably out of each compressed region during a period; an isothermal response would give a different derivative. In a moving uniform medium, the acoustic characteristics travel at $U\pm c_s$ along the flow, rather than at $\pm c_s$ in the laboratory.

For a progressive wave, $p^{\prime }={\rho }_0{c}_s{u}_{\parallel }^{\prime }$; pressure, density and longitudinal velocity are in phase. The acoustic energy density is the sum of kinetic and compression energies,Sound transports energy while fluid elements execute small oscillations about their mean positions. The ideal uniform-fluid acoustic mode is nondispersive. Spherical spreading reduces pressure amplitude as $1/r$ when total power is conserved; viscosity and heat conduction attenuate waves, while sufficiently large compressions can steepen into shocks. In a stratified astrophysical gas the uniform result is local: wavelengths short compared with background scale heights justify it, whereas long waves can couple to buoyancy and acoustic-cutoff behavior.

Treat the gas as a smooth quasi-one-dimensional flow through a rigid passage, neglecting axial body forces and heat input. It is the confining medium that is incompressible; the flowing gas must remain compressible. With area $\mathcal{A}(x)$ and axial speed $u$, continuity and inviscid momentum balance giveCombining their logarithmic differentials yields the nozzle area-velocity relationA smooth crossing of $\mathcal{M}=1$ requires an area extremum. The given radius is positive and has derivative $r^{\prime }(x)=16r_{0}x\exp [-16x^2/L^2]/L^2$, which vanishes only at $x=0$ within the slab and changes from negative to positive there. Thus this station is the unique throat, with $r_s=r_0/2$ and ${\mathcal{A}}_s=\pi r_0^2/4$. The specified subsonic inlet and supersonic outlet select the smooth accelerating branch, soA throat by itself does not guarantee transonic flow for arbitrary reservoir and back pressures; that branch is specified here. Before the throat a converging section accelerates subsonic flow; beyond it a diverging section accelerates supersonic flow.^[ref]

For the temperature task, evaluate the mass flux and sound speed at the throat using $p=K{\rho }^{\gamma }$:The ideal-gas relation then givesThis expresses the sonic temperature in terms of the entropy constant, geometry, mass flux and composition; no numerical value of $K$ or $\dot{M}$ is provided. Equivalently $c_s^{\gamma +1}=\gamma K(\dot{M}/{\mathcal{A}}_s{)}^{\gamma -1}$ and $T_s=\mu c_s^2/(\gamma R^{\ast })$. The slab length drops out because the sonic thermodynamic state depends on the minimum area, not on the axial distance over which a lossless smooth nozzle contracts. As a check, conserved stagnation enthalpy gives $T_s=2T_{\mathrm{stag}}/(\gamma +1)$.^[ref]

The result assumes cross-sectional thermodynamic uniformity and negligible transverse kinetic energy. Exact multidimensional flow through a sharply varying passage can have a curved sonic surface; the radius profile alone does not solve that boundary-value problem. The original numerical plot below computes the standard quasi-one-dimensional branch, rather than reproducing the supplied passage sketch.

An accretion disk forms when infalling material has enough angular momentum to orbit before it can reach the central object. In a thin disc dominated by a central mass $M$, the orbital motion is approximately Keplerian rotation:Vertical hydrostatic equilibrium gives the disk scale height $H\sim c_s/{\Omega }_K$; a geometrically thin disc has $H/R\ll 1$. The surface density is $\Sigma =\int \rho dz$, and a steady inward mass flow rate is $\dot{M}=-2\pi R\Sigma v_R>0$.

Moving inward requires losing specific angular momentum. Internal stresses transport angular momentum outward; much of the mass can then accrete while an outer fraction spreads to larger radius. Molecular viscosity alone is generally inadequate for rapid astrophysical accretion. Turbulent and magnetic stresses, including those generated by the magnetorotational instability, provide effective transport. A common parametrization is $\nu =\alpha c_{s}H$ for an effective kinematic viscosity, with $\alpha$ a dimensionless stress parameter.

For a thin Keplerian accretion disk, the evolution equation isThis viscous evolution of an accretion disk occurs on the viscous timescale $t_{\nu }\sim R^2/\nu$, much longer than an orbital period in a thin disc with moderate effective stress.

The orbital specific energy is $-GM/(2R)$, so inward motion releases binding energy. In a steady, optically thick, radiatively efficient disc with zero torque at its inner edge, the emergent flux from each face isThe Stefan–Boltzmann law identifies the local effective temperature. Far from the inner boundary, $T_{\mathrm{eff}}\propto R^{-3/4}$, so the emitted spectrum combines hotter inner and cooler outer annuli. Integrating both faces gives $L_{\mathrm{disc}}=GM\dot{M}/(2R_{\mathrm{in}})$ under these assumptions. This is not a universal spectrum or efficiency: optical depth, inner torque, energy advection and outflows can change the radiation. The essential connection is outward angular-momentum transport, inward accretion, and release of binding energy.

Interpret $v$ as the normal upstream speed in the shock rest frame and $\rho$ as its upstream mass density. If the quoted jet speed is in another frame, the shock speed must first be supplied to obtain this relative velocity. Use a calorically perfect gas of mean particle mass $\bar{m}=\mu m_u$.

In the strong-shock limit, ${\mathcal{M}}_1\gg 1$. The Strong-shock Rankine-Hugoniot conditions giveThe pressure expression also follows from the supplied ratio by inserting ${\mathcal{M}}_1^2=\rho v^2/(\gamma p_1)$. The ideal gas law then gives the immediate postshock temperatureThe upstream thermal contribution is negligible at this order. This is a temperature in kelvin, rather than a temperature change per unit time.

To turn the postshock thermal energy into a rate, the processed mass per unit time is $\dot{m}=\rho vS$. Its downstream specific internal energy isThus the thermal internal energy carried into the shocked gas per unit time isPer unit area, divide this expression by $S$. Its dimensions are watts; the mean molecular weight cancels because the heat capacity per unit mass is inversely proportional to it, while $T_2$ is proportional to it.

There is a second common convention for shock heating power: the loss of directed bulk kinetic energy flux. The energy fluxes of a strong adiabatic shock giveThis equals the downstream enthalpy flux, and includes pressure work as well as the advected internal energy. For $\gamma =5/3$, the compression ratio is $4$, $T_2=3\mu m_u{v}^2/(16k_B)$, ${\dot{E}}_{\mathrm{int}}=9\rho S{v}^3/32$, and ${\dot{E}}_{\mathrm{kin},\mathrm{loss}}=15\rho S{v}^3/32$. Specifying the energy convention avoids conflating these rates. A rate $dT/dt$ for a fixed slab cannot be inferred without its participating mass and thermal evolution; the supplied data fix the postshock temperature and the processed energy fluxes.

Assume a steady spherical outflow of fixed composition. Let $\bar{m}$ be its mean particle mass, and define the constant isothermal sound speed by $a^2=k_B{T}_0/\bar{m}$. Subsonic here refers to this isothermal closure, so $0<u<a$. Mass conservation and radial momentum balance giveLogarithmically differentiating the first equation yields $d\ln \rho /dr=-2/r-u^{\prime }/u$. Substitution eliminates the density gradient and gives the isothermal Parker wind equationAt a smooth speed maximum, $u^{\prime }=0$, while the left prefactor is finite and negative. The right side must vanish, henceFor $r<r_M$ both the right side and the prefactor are negative, so $u^{\prime }>0$; for $r>r_M$, $u^{\prime }<0$. This is a true maximum on a subsonic isothermal breeze, not the accelerating transonic branch.

There is an underdetermination in the requested density estimate. Introduce the peak Mach number ${\mathcal{M}}_M=u_M/a$, where $0<{\mathcal{M}}_M<1$. The density at the maximum speed of an isothermal breeze is exactlyThe hypotheses do not fix ${\mathcal{M}}_M$. A boundary density or pressure, or an equivalent wind integration constant, is needed. The definite statement following from subsonicity is the lower boundIf a maximum speed of order $a$ is intended as an extra estimate, then ${\rho }_M\sim {\rho }_{\mathrm{char}}\propto \dot{M}{T}_0^{3/2}{M}^{-2}$. A very slow breeze can have a much larger density, so setting $u_M=a$ as an exact equality would contradict the stated always-subsonic assumption. The mean particle mass is another implicit fixed-composition input needed to express a temperature as a squared sound speed.

One can see the free parameter directly by integrating the wind equation:Different $C$ give different subsonic maxima at the same radius. The original plot illustrates this nonuniqueness. These breeze branches require nonzero ambient pressure at large radius; they are not the usual transonic escape into a low-pressure exterior.

Maintaining an isothermal expanding wind generally requires heat exchange; isothermal motion is not adiabatic motion.^[ref]

Write the initial mass density as $\rho =mn$. The ideal gas law gives pressure $p=n{k}_{B}T$, and the sphere has $M=4\pi \rho R^3/3$. Here $k_B$ is the Boltzmann constant and $G$ is the gravitational constant. First obtain the gravitational potential energy by adding concentric shells. The mass inside radius $r$ is $M(r)=M{r}^3/R^3$, soIn a small contraction $dR<0$ at fixed total mass, the gravitational energy released is $-dW=3G{M}^2|dR|/(5R^2)$. The positive work required against the gas pressure is $-pdV=4\pi R^{2}p|dR|$. Equating these givesUse $R=(3M/(4\pi mn){)}^{1/3}$ to eliminate the radius. The uniform-sphere pressure-work collapse criterion isThe units are kilograms: $(k_{B}T/(Gm){)}^{3/2}$ has units $\mathrm{k}{\mathrm{g}}^{3/2}{\mathrm{m}}^{-3/2}$, and ${\rho }^{-1/2}$ supplies $\mathrm{k}{\mathrm{g}}^{-1/2}{\mathrm{m}}^{3/2}$. At fixed initial $n$ and $T$, $M>M_{\mathrm{crit}}$ makes gravity stronger than the pressure-work requirement in this model, allowing contraction.

This is a Jeans mass estimate for homologous contraction of an initially uniform sphere. It is not an exact stability limit for a hydrostatic gas sphere: uniform mass density and uniform temperature would give uniform pressure, rather than the gradient required by hydrostatic equilibrium. The numerical coefficient therefore differs from conventions based on a Jeans instability wavelength or a pressure-confined isothermal sphere. Subsequent contraction also depends on heat exchange; isothermal pressure support during gravitational collapse weakens relative to gravity as the radius decreases.

Following a fluid element. Let $\mathbf{X}(t)$ be a pathline, satisfying $\dot{\mathbf{X}}=\mathbf{u}(\mathbf{X},t)$. The material derivative, also called the Lagrangian time derivative, is the rate observed while following that element. The multivariable chain rule givesThe partial derivative holds position fixed; the advection term accounts for motion through a spatially varying field. For a vector field, the same formula applies component by component in fixed Cartesian axes.

Momentum balance. For an inviscid flow, pressure is the only surface stress. A small material volume $dV$ has net pressure force $-\nabla pdV$, obtained by the divergence theorem. Its mass is $\rho dV$ and its material acceleration is $D\mathbf{u}/Dt$. Applying Newton's second law gives the Euler equations for an inviscid fluid,This is the local momentum conservation equation when viscosity and physical body forces are absent. A prescribed body force per unit mass would add $\rho \mathbf{f}$ on the right.

Changing to rotating coordinates. In a rotating frame with constant angular velocity $\mathbf{\Omega }$, the lab acceleration is the rotating-frame material acceleration plus $2\mathbf{\Omega }\times \mathbf{u}+\mathbf{\Omega }\times (\mathbf{\Omega }\times \mathbf{r})$. ThereforeFor uniform rotation,With constant ${\rho }_0$, this gives the Euler equation in a uniformly rotating frame,The term subtracted from $p$ is the equilibrium centrifugal contribution. Dividing the bracket by ${\rho }_0$ defines a reduced pressure $\Pi =p/{\rho }_0-|\mathbf{\Omega }\times \mathbf{r}{|}^2/2$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-2}$. The Coriolis acceleration remains a velocity-dependent term.

Linear disturbances. Linearize about rest in the rotating coordinates, with the equilibrium part of $\Pi$ removed. The advective acceleration is quadratic in disturbance amplitude, so the linearized Euler equations and incompressible flow condition areTaking the curl removes the pressure gradient. The identitygives the inertial-wave vorticity equation,

Plane-wave frequencies. Let $\mathbf{k}\ne 0$ and use the given plane wave convention. Incompressible flow requires $\mathbf{k}\cdot {\mathbf{u}}_0=0$, and the vorticity equation givesFor a nonzero frequency, apply $\mathbf{k}\times$ again, using $\mathbf{k}\times (\mathbf{k}\times {\mathbf{u}}_0)=-k^2{\mathbf{u}}_0$. It follows thatThis is the inertial-wave dispersion relation. If $\theta$ is the angle of the wavevector to the rotation axis, the positive angular frequency is $2|\Omega \cos \theta |$. It is largest, $2|\Omega |$, for a parallel wavevector, decreases to zero for a perpendicular one, and is independent of the magnitude of the wavenumber at fixed direction. The two signs describe the two oscillatory polarization branches. When $k_z=0$, the zero-frequency result is obtained directly from the linear equations rather than by dividing by $\omega$; it describes the nonoscillating balanced limit.

A close binary star can produce two compact stellar objects if both initial stellar masses are sufficiently large. The initially more massive star evolves first, transfers material or loses its envelope, and leaves a neutron star or black hole after core collapse. The system must remain bound despite mass loss and any natal kick. Subsequent evolution of the companion can produce a common envelope: orbital energy transferred to the envelope can expel it and greatly reduce the separation. A second collapse must again leave a bound pair. A sufficiently close surviving binary star then loses orbital energy and angular momentum through gravitational waves and may merge. Stable transfer, winds and kicks make this a possible channel rather than an inevitable outcome; dynamical encounters in dense stellar systems provide another channel.

For two black holes of equal mass, $M=20M_{\odot }$ and $\mu =5M_{\odot }$. The prescribed instantaneous orbital decay timescale givesThe integrated ten-billion-year criterion gives $51.9$ h. These large initial periods compared with a neutron-star pair reflect the stronger radiation reaction of the more massive system.

Use constant component masses and a circular Keplerian orbit. Then $J=\mu \sqrt{GMA}$ gives $\dot{J}/J=\dot{A}/(2A)$, while $P=2\pi A^{3/2}/\sqrt{GM}$ gives $\dot{P}/P=3\dot{A}/(2A)=3\dot{J}/J$. Dividing the supplied loss law by $J$ therefore yieldsThe positive instantaneous orbital decay timescale is ${\tau }_P=P/|\dot{P}|$, since $\dot{P}<0$. For $t_G=10^{10}$ yr, the requested prescription givesIt is important to distinguish this prescription from the integrated gravitational-wave inspiral time. Integrating $\dot{A}=-C/A^3$ givesFinite contact radii supply a small lower-limit correction when the initial separation is much larger than the compact objects. Thus the periods below satisfy the printed instantaneous-timescale criterion and have a formal merger time $3.75$ Gyr. If instead the requirement is literally $t_{\mathrm{merge}}=10$ Gyr, every period limit is larger by $(8/3{)}^{3/8}=1.445$. Eccentric systems require a different decay calculation.

Infer the component masses and spins from the gravitational wave phasing, using the chirp mass together with binary mass ratio information. The paper's representative pairs have total masses $2.8$, $11.4$ and $20M_{\odot }$, and the more massive mergers terminate their inspirals at lower frequencies. A black-hole ringdown and its quasi-normal modes identify the final compact remnant, although both a neutron-star–black-hole merger and a sufficiently massive neutron-star pair can also leave a black hole.

Tidal deformability and any disruption cutoff provide a more direct signature of a neutron star than total mass alone. A disrupted neutron star may supply a gamma-ray burst, kilonova and afterglow; these help distinguish matter-containing mergers from a clean black-hole pair. Their absence does not establish a black-hole pair, because jets are beamed and ejecta may be faint or swallowed. Measured gravitational wave masses are redshifted detector-frame masses; association with a host can supply the redshift needed for source-frame masses. Classification should combine the waveform and any electromagnetic evidence, allowing for parameter degeneracies.

An ultraluminous X-ray source is a compact, apparently off-nuclear source in another galaxy with an inferred isotropic X-ray luminosity typically above $10^{32}$ W, exceeding the classical Eddington luminosity of an ordinary stellar-mass accretor. Imaging X-ray surveys revealed unusually bright sources away from galactic nuclei; distance association with the host converted the observed fluxes into surprisingly high luminosities. Better angular resolution helps exclude blends and background active nuclei. The 2014 discovery of coherent pulsations established a neutron star in M82 X-2, so the class cannot be identified exclusively with black holes.^[ref]

An intermediate-mass black hole can radiate approximately isotropically below its own Eddington luminosity: this avoids extreme super-Eddington emission, but the inferred mass requires dynamical or suitable spectral/timing support and is not proved by luminosity alone. A stellar-mass accretor with a supercritical accretion disk and a wind funnel can combine genuinely super-Eddington output with geometrical beaming; this is compatible with a neutron-star accretor but requires a consistent radiation/flow/torque model. Relativistic beaming from a jet is another possibility for some objects, but orientation, spectral and variability evidence must support it. These mechanisms need not apply identically to every member of the observed class.

For spherical radiation with flux $F=L/(4\pi r^2)$, the outward acceleration is $\kappa F/c$. Equating it to $GM/r^2$ givesThe limit assumes the stated opacity, approximate spherical symmetry and coupling between the radiatively accelerated plasma and the accreting matter. It is not a universal ceiling on isotropic-equivalent luminosity for an anisotropic magnetized flow.

Write the available accretion power as $L=\eta \dot{M}{c}^2$. ThenAccretion as the energy source alone does not fix $\eta$. For radiatively efficient Newtonian accretion onto a slowly rotating stellar surface, $\eta \simeq GM/(R{c}^2)$ and the corresponding surface accretion Eddington rate isThis is the appropriate convention for the numerical neutron-star estimates below, rather than imposing a separate reference efficiency of $0.1$.

Using $M=1.4M_{\odot }$, $R=10^4$ m and the stated opacity givesIf the observed X-rays represent isotropic, radiatively efficient surface accretion power, thenA bolometric correction, beaming, radiation advected into the flow or energy carried by a wind would change the mapping from the measured band luminosity to the supplied mass rate. This estimate keeps the assumptions needed for comparison with the torque calculation.

Write $R=C{E}^a{\rho }_0^b{t}^d$ for a uniform ambient mass density. Equating mass, length and time dimensions gives $a+b=0$, $2a-3b=1$ and $-2a+d=0$. Thus the Sedov-Taylor blast wave hasThe dimensionless coefficient depends on the gas's adiabatic index and on whether $R$ denotes the shock radius. Solving the similarity profile for a monatomic ideal gas gives $\xi \simeq 1.15$, so taking it as unity is an order-of-magnitude approximation rather than an extra dimensional result.

Timed photographs of the 1945 atomic explosion provide a direct scale test: the measured fireball radius follows $R\propto t^{2/5}$ over the relevant blast-wave interval, and inserting the independently established explosion energy and atmospheric density gives a coefficient of order unity. The published photographic analysis tested the radius-time law and used its normalization to infer an explosion energy of the correct nuclear-test scale. This illustrates both the success and the limits of dimensional reasoning: the power law is fixed dimensionally, whereas the numerical coefficient requires a profile calculation or measured normalization.^[ref]

Neglecting the exterior temperature requires a strong shock, ${\mathcal{M}}_s=V_s/c_{s,u}\gg 1$, or ${\rho }_0{V}_s^2\gg p_u$. The swept-up mass must dominate the initial ejecta for this energy-conserving similarity stage, while radiative losses, gravity, ambient gradients and continuing energy input must be unimportant. A very hot exterior or a late weak disturbance introduces a sound-speed scale and invalidates the cold-upstream similarity assumptions.

Use the local shock rest frame. Let $u$ and $d$ label upstream and downstream gas, with positive normal speeds $v_u,v_d$ through the shock. The unshocked ISM is almost stationary in the stellar frame, so $v_u\simeq V_s$, whereas the shocked gas moves outward at $V_s-v_d$ in that frame. For an ideal monatomic gas, $\gamma =5/3$; neglect upstream pressure and enthalpy in the strong-shock limit. Conservation of mass, momentum and energy per unit mass givesPut $q={\rho }_d/{\rho }_u=v_u/v_d$. The momentum equation gives $p_d={\rho }_u{v}_u^2(1-1/q)$. Substitution into the energy equation givesDiscarding the unshocked root $q=1$ leavesThe downstream speed in the stellar frame is therefore $3V_s/4$, not $V_s/4$. These conditions describe immediate shock heating; subsequent cooling and evolution behind the front require the downstream flow profile. Cosmic-ray pressure, radiative energy loss inside the shock layer or a different equation of state would alter the compression ratio.

Let the mean particle mass be $\bar{m}=\mu m_H$, with $n_u$ counting total particles so that ${\rho }_u=n_u\bar{m}$. The ideal-gas law and the strong-shock pressure give the immediate post-shock temperatureThe robust prediction is $T_d\propto t^{-6/5}$. For a transparent numerical convention take $\mu =1$, $n_u=10^6{\mathrm{m}}^{-3}$ and $\xi =1$. ThenFor an ionized mixture with $\mu =0.6$ under the same total-particle-density convention, these normalizations become $5.13$ pc and $5.49\times 10^7$ K. The paper does not specify composition or whether its quoted density instead counts hydrogen nuclei; that convention must be stated, not hidden in a precise age.

The ambient-temperature limit of an adiabatic blast wave gives a formal matching estimate, $T_d\sim 10^4$ K, givesWith $\mu =0.6$ the age is about $1.3\times 10^6$ yr and the radius stays about 90 pc. This matching criterion already lies outside a very strong shock: $T_d=T_u$ in the strong-shock formula corresponds to ${\mathcal{M}}_s=\sqrt{16/(3\gamma )}\simeq 1.79$. A separate sonic estimate $V_s\sim c_{s,u}$ gives about $4.4\times 10^6$ yr and 133 pc for $\mu =1$; neither extrapolation is a precise transition solution. Using $\xi =1.15$ increases these inferred matching ages and radii by factors ${\xi }^{5/3}$.

In a real supernova remnant, substantial radiation losses commonly end the energy-conserving phase before these formal pressure-merging limits. The relevant comparison is $t_{\mathrm{cool}}\sim (3/2)n{k}_{B}T/(n_e{n}_i\Lambda (T))$ against the expansion time. A cooling function, composition and ionization state are needed to calculate that earlier age and radius; they are not supplied. The numerical values above answer the ambient-temperature breakdown estimate under explicit assumptions and must not be presented as a uniquely determined radiative shell-formation age.

After cooling, an interior-pressure-driven shell can approximately follow $R\propto t^{2/7}$ while pressure remains important. In the limiting momentum-conserving stage, $M_s\dot{R}\simeq \mathrm{const}$ and $M_s\propto {\rho }_u{R}^3$ give $\boxed{R\propto t^{1/4}}$. Ultimately ambient pressure and turbulent motions make the disturbance merge with the ISM. These slower stages respect the reduced energy budget rather than continuing an adiabatic strong-shock law indefinitely.

Stellar radiation and nebular ionization. The planetary-nebula spectral energy distribution is constrained by a shared radiation and thermal balance. The central star's effective temperature, luminosity and evolutionary state set its continuum and ionizing-photon supply. A hotter star shifts the ionizing spectrum and allows higher ionization stages; the total recombination output also depends on how many ionizing photons the gas actually absorbs. In an ionization-bounded nebula, approximate photoionization balance iswhereas a density-bounded nebula can leak photons. Dust absorption of ionizing photons changes this relation further. Stellar-atmosphere and wind effects make a blackbody an approximation, especially in the extreme ultraviolet. A changing core temperature and expanding envelope explain why nebular spectra evolve.

Nebular continua and recombination lines. Their strength depends on the emission measure, roughly $\int n_e{n}_{i}dV$, electron temperature, ionic abundances and radiative branching probabilities. For a recombination line,Case-B line ratios assume appropriate trapping of Lyman-series photons and escape of the lower-series radiation. Density, optical-depth effects and extinction can alter observed ratios. The free-free emissivity has the schematic dependence $j_{\nu }^{\mathrm{ff}}\propto n_e{n}_i{T}_e^{-1/2}\exp [-h\nu /(k_B{T}_e)]g_{\mathrm{ff}}$; the nebular free-bound continuum also carries ionization thresholds and recombination cross-sections. Two-photon emission depends on the population of the hydrogen $2s$ state and can be reduced by collisional mixing into $2p$. Density clumping matters because these rates are weighted approximately by density squared, not simply by total gas mass.

Collisionally excited lines. Below a line's critical density of an atomic transition,Elemental abundance, ionic fraction and electron temperature therefore all enter. The ionizing spectrum and local ionization parameter help determine which ion is present. Above the critical density of an atomic transition, collisional de-excitation competes with radiative decay and can suppress a forbidden line. Low-energy infrared fine-structure lines are generally less exponentially temperature sensitive than high-excitation optical/ultraviolet lines. Metal-line cooling also helps set the temperature, so changing an abundance changes both the number of emitters and the thermal balance; line luminosity need not be proportional to abundance at fixed everything else.

Dust and molecules. Nebular dust spectral features identify the carriers as well as contributing to the integrated luminosity. In an optically thin dust component,The grain mass sets its normalization, while the radiation field, stellar distance, grain size and absorption efficiency determine its temperature distribution. For approximately grey heating, radiative equilibrium suggests $T_d\propto L_{\ast }^{1/4}{r}^{-1/2}$; non-grey grains change this scaling. Long-wavelength opacity is often approximated as ${\kappa }_{\nu }\propto {\nu }^{{\beta }_d}$. Small aromatic grains undergo transient heating, so one equilibrium temperature is not adequate for every infrared feature. Carbon-to-oxygen chemistry, grain processing and shielding affect aromatic versus silicate bands. Molecular line emission further depends on molecule abundance, collisional density, kinetic temperature, ultraviolet pumping and shocks, while shielding and dissociation set whether molecules survive at all.

Finally, interstellar extinction preferentially suppresses ultraviolet and optical light, and dust reradiates some absorbed energy in the infrared. Aperture, angular resolution and distance must be treated consistently before comparing stellar, nebular and dust components. These dependencies make the combined spectrum an energy and abundance diagnostic, but also prevent a unique reconstruction from a single broad-band colour.

Cold molecular hydrogen is difficult to observe directly: its homonuclear symmetry forbids electric-dipole rotational transitions, and its lowest useful quadrupole rotational excitation is high compared with a cold cloud's temperature. Carbon monoxide is relatively abundant, has a permanent dipole moment and has low-lying rotational transitions that can be collisionally excited in cold gas. The low-$J$ lines lie in accessible millimetre bands and suffer much less dust extinction than visible light. They therefore reveal molecular-gas velocities and spatial structure efficiently.

The advantages do not turn CO into a direct mass counter without a model. Common-isotope lines can be optically thick; excitation can be subthermal or have multiple temperatures. Abundance changes with metallicity, chemistry, shielding and photodissociation, and CO can freeze onto grains. Molecular hydrogen can also survive where CO is weak, giving CO-dark gas. A calibrated carbon-monoxide tracer of molecular hydrogen conversion may be useful for appropriate cloud populations but should not be transplanted uncritically to newly formed supernova ejecta. The line calculations here infer CO mass, not total molecular mass unless a CO-to-hydrogen abundance is separately known.

Assume a single-temperature optically thin local thermodynamic equilibrium population, negligible line absorption/stimulated corrections, and that the measured line luminosity is from the remnant's CO. Each spontaneous $J=1\to 0$ decay emits $h{\nu }_{10}$, givingFor a rigid rotor, $E_J/k_B=T_{r}J(J+1)$ with $T_r=h{\nu }_{10}/(2k_B)=2.767$ K. The degeneracy is $g_J=2J+1$, and the rotational partition function isAt 20 K the direct sum gives $Z\simeq 7.57$, $f_1=0.3005$; the approximation $Z\simeq T/T_r+1/3$ is already adequate. Thus, with molecular mass $m_{\mathrm{CO}}\simeq 28m_u$,The remnant age does not enter this optically thin population count. It becomes relevant when estimating emitting volume, density, expansion and optical depth to test the assumptions.