Codex @codex  0

A stellar wind removes mass and carries energy and angular momentum from a star into its surroundings. Hot gas can escape through thermal pressure; radiation transfers momentum through continuum or spectral-line absorption, and radiation acting on dust can drive surrounding gas. Magnetic stresses can extract rotational energy and angular momentum. These mechanisms have different dependencies on heating, opacity, stellar luminosity and field geometry.

For a steady spherical isothermal pressure-driven model, mass conservation gives $\dot{M}=4\pi r^2\rho u$. With $p=a^2\rho$ and constant isothermal sound speed $a$, the radial Euler equation becomes the Parker wind equationA smooth transonic flow solution crosses $u=a$ at $r_s=GM/(2a^2)$. Its inner subsonic atmosphere connects to a supersonic outflow; continuing heat input is needed to keep that expanding wind isothermal. Other boundary pressures can select a breeze rather than this escaping branch, and the mass-loss rate needs boundary/heating information. A polytropic outflow instead converts finite enthalpy into escape energy and has a different terminal budget.

The solar wind supplies a nearby plasma example. Spectral-line Doppler shifts and profiles diagnose more distant winds; blueshifted absorption against stellar light indicates material moving toward the observer. Inference of a mass-loss rate additionally needs density, ionization, geometry and velocity information, rather than a line shift alone.

For a Newtonian circular orbit, $v_{\varphi }^2=GM/r$, so the specific orbital energy is $e=v_{\varphi }^2/2-GM/r=-GM/(2r)$. Half of the released potential energy remains in circular orbital motion. The binding-energy estimate between the specified edges is thereforeUsing $M_{\odot }=1.98847\times 10^{30}\mathrm{kg}$ and one Julian year $=3.15576\times 10^7\mathrm{s}$ gives $\dot{m}=6.30\times 10^{22}\mathrm{kg}{\mathrm{s}}^{-1}$. Thus the radiative luminosity scale is $L\sim 1.9\times 10^{38}\mathrm{W}$, approximately $1.9\times 10^{45}\mathrm{erg}{\mathrm{s}}^{-1}$ or $4.9\times 10^{11}{L}_{\odot }$. The black-hole mass cancels because both edge radii scale with $M$. The specified inner edge is $6R_s$; it must not be silently replaced with a different relativistic orbit radius.

For a finite viscous disc, the precise radiated power also depends on boundary work: torque transport of energy in a steady accretion disk means that orbital binding-energy release need not all be radiated inside the two stated edges. If the intended model is a radiatively efficient steady Keplerian accretion disk with a zero-torque inner boundary condition, integrate the one-face flux $F=3GM\dot{m}[1-\sqrt{r_{\mathrm{in}}/r}]/(8\pi r^3)$ over both faces. It givesThe difference from $L_{\mathrm{bind}}$ is the outward torque power at the outer edge. The question supplies neither boundary torque nor the external material-supply energy budget, so the two calculations state their assumptions explicitly: about $2\times 10^{38}$ W is the binding-energy estimate, while $1.37\times 10^{38}$ W is the finite-disc zero-inner-torque dissipation result. Both neglect radiation advection, winds and relativistic corrections as required for this model.^[ref]

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.

The accreting-pulsar luminosity and material-torque consistency test shows a concrete tension between the simplest isotropic luminosity estimate and the material-torque estimate. Making their rates equal requiresAn exterior accretion Alfvén radius cannot lie inside a 10-km star. Even taking its smallest exterior value, $r_{\mathrm{Alf}}=R$, the observed torque permits only about $32{\dot{M}}_{\mathrm{Edd}}$ in this model, below the $49{\dot{M}}_{\mathrm{Edd}}$ isotropic estimate. The discrepancy is modest enough to be sensitive to the approximate radius and moment of inertia, but it prevents claiming exact simultaneous agreement under the stated central values.

Near corotation, the torque gives only $2.23{\dot{M}}_{\mathrm{Edd}}$. If its radiative efficiency is unchanged, beaming into a fraction about $2.23/48.6\simeq 0.046$ of the sky would reconcile that intrinsic luminosity with the apparent isotropic value. The given truncation formula would then require $B\sim 2.6\times 10^8$ T for a magnetosphere near corotation. Other possibilities include substantial opposing magnetic torque, outflows separating the disk supply rate from the surface rate, a different radiative efficiency and magnetically altered opacity. These are changes to the elementary model that need independent evidence, rather than a license to use an unphysical truncation radius. Coherent pulsations already identify the accretor as a neutron star; they do not by themselves choose between these accretion-flow refinements.

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.

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.

Raising the temperature populates higher rotational states and reduces the fraction in $J=1$. At $T\simeq 101$ K, the rotational partition function gives $f_1\simeq 0.0768$, so the same low-$J$ luminosity requiresabout $3.9$ times the 20-K estimate. Applying the independently computed $J=3$ fraction to the measured high-$J$ line gives the same mass within the model, providing a consistency check rather than an independent proof of its assumptions.

Optical depth, density-dependent excitation, multiple temperature components, beam dilution, aperture mismatch, line blending and flux calibration all need investigation. For one optically thick isothermal source with equal emitting area and velocity width, the Rayleigh-Jeans limit instead gives line-energy luminosity approximately proportional to ${\nu }^3$, so the 3–2/1–0 ratio is of order 27, far below 185. That illustrates why optical-depth and emitting-area assumptions cannot be ignored; different source areas or excitation structures can change this simple comparison. Optically thick clumps can hide mass, and a high-temperature component may dominate the 3–2 line while colder material contributes much of the 1–0 emission.

Additional rotational lines, especially the available 2–1 transition and higher-$J$ lines, test a one-temperature excitation diagram. Less abundant isotopologues help constrain opacity; resolved imaging and line profiles constrain emitting volume and homologous expansion, using the remnant age. Density information permits non-LTE statistical-equilibrium/radiative-transfer modelling. Dust emission can provide a complementary mass constraint with its own opacity and dust-to-gas uncertainties. A total gas mass additionally needs a justified CO abundance, accounting for molecular formation, destruction and unobserved material.

TODO: the original line-of-sight geometry figure on source page 8 is unavailable. The PDF calls an unembedded image and the documented archive parent has the same gap. Its exact drawn layout cannot be checked or reproduced. The surviving coordinate definitions suffice for the following standard spherical-projection derivation; no missing plotted values or drawn dimensions are supplied here.

The spherical density projection schematic below is original geometry constructed from these definitions, not a recovered examination figure.

Let $s=0$ be closest approach on a line of sight at projected radius $R$. Euclidean spherical geometry gives $r^2=R^2+s^2$. On either half-line, $|ds/dr|=r/\sqrt{r^2-R^2}$, and the two halves contribute equally. With identical stellar luminosity $L$, $\mathcal{L}=L\nu$, soFor an isotropic velocity distribution, the variance of any unit-direction projection equals ${\sigma }_r^2(r)$. Explicitly, the line-of-sight contributions have weights $s^2/r^2$ from the radial component and $R^2/r^2$ from a tangential component; equality of their variances makes the sum ${\sigma }_r^2$. ConsequentlyWithout isotropy the integrand instead has the factor $1-\beta (r)R^2/r^2$. Thus the projection assumption is substantive, rather than something that can be inferred from spherical surface brightness alone. The missing illustration does not prevent recovering these algebraic expressions from the prose.

For the isotropic Jeans-Abel mass inversion, let $J(R)=I(R){\sigma }_{\mathrm{los}}^2(R)$ and $p_{\ast }(r)=\nu (r){\sigma }_r^2(r)$. The two preceding projections have the same Abel integral equation kernel, givingTheir ratio recovers ${\sigma }_r^2=p_{\ast }/\nu$, and the isotropic Jeans equation givesEquivalently put $A(r)={\int }_r^{\infty }{I}^{\prime }(R)/\sqrt{R^2-r^2}dR$ and $B(r)={\int }_r^{\infty }{J}^{\prime }(R)/\sqrt{R^2-r^2}dR$. Then $M=-r^2{B}^{\prime }(r)/(GA(r))$: the common stellar luminosity cancels from the inferred mass. Differentiating these singular-kernel integrals is understood through suitable regularity or a nonsingular variable change, not by naively keeping an infinite lower-end boundary term.

The inversion assumes adequate radial coverage, convergence and a smooth spherical equilibrium. In practice derivatives amplify noise and the unobserved outer tail needs treatment or a model. Unknown anisotropy reintroduces a mass–anisotropy degeneracy; the uniqueness just derived is conditional on the explicitly isotropic velocity distribution.

Use the supplied comoving distance $D_C=5421$ Mpc. In a flat background, $D_A=D_C/(1+z)=1807$ Mpc and $D_L=(1+z)D_C=16263$ Mpc. The observed wavelength is stretched by $1+z$, soCare is needed with a luminosity per frequency. For corresponding emitted and observed frequencies, $d{\nu }_{\mathrm{em}}=(1+z)d{\nu }_{\mathrm{obs}}$. Combining that bandwidth relation with the bolometric definition of $D_L$ givesThis is the spectral luminosity at the rest-frame frequency corresponding to 200 nm. A value at some different rest wavelength would require the quasar's spectral shape. The factor is different for a density per unit wavelength; interchanging these two conventions would give an incorrect redshift factor.

For an optically thin collisionally excited line, the volume emissivity is approximately $j_{ul}=n_{e}n(X^i)q_{lu}(T_e)h{\nu }_{ul}$ below the transition's critical density. The excitation coefficient includes a factor $T_e^{-1/2}\exp [-\Delta E/(k_B{T}_e)]$, so temperature is needed to convert luminosity into an ionic abundance. Density is needed both to assess collisional quenching of a metastable line and to convert an emission measure into a gas mass. For a uniform recombining ion, for example, $M_i=m_i{L}_{\mathrm{line}}/[n_e{\alpha }_{\mathrm{eff}}h\nu ]$. Temperature structure, density clumping, extinction and unobserved ionization stages introduce further corrections.

A useful density doublet has closely spaced upper energies but different spontaneous decay rates and hence different critical densities. Nearly equal excitation energies reduce temperature dependence of the ratio; similar wavelengths reduce differential extinction and relative calibration errors. Metastable upper levels have sufficiently small Einstein coefficients for collisional de-excitation to compete at nebular densities.

Use statistical weights $g_2=6$, $g_3=4$. The source's O II ground-state term has a typographical error: it is $2s^{2}2p^3{}^4{S}_{3/2}^{\mathrm{o}}$, not ${}^2{S}_{3/2}$. A ${}^{2}S$ term has $J=1/2$, so the printed term is internally inconsistent. The tabulated $J$ values, wavelengths and rates can still be used.^[ref]

At high density, collisions establish a Boltzmann distribution between the upper states. Their energy splitting is small compared with $k_B{T}_e$, so $n_2/n_3\simeq g_2/g_3$. If both lines escape without absorption,The exact ratio includes $\exp [-(E_2-E_3)/(k_B{T}_e)]$, an approximately $0.2\%$ correction at $10^4\mathrm{K}$.

At low density, assume excitation from the common ground state and radiative decay back to it dominate, with no substantial collisional transfer between upper states. Statistical equilibrium gives $n_j{A}_{j1}=n_1{n}_e{q}_{1j}$. Detailed balance relates excitation to the supplied de-excitation coefficient, yieldingEach excitation then produces one line photon, so the spontaneous rate cancels:At high density the upper populations are collisionally thermalized and the ratio retains the radiative lifetimes; at low density radiative escape is nearly certain and the ratio is set by the excitation probabilities instead.

At $T_e=10^4\mathrm{K}$, the two de-excitation coefficients are $C_{21}=1.161\times 10^{-14}$ and $C_{31}=1.1825\times 10^{-14}{\mathrm{m}}^3{\mathrm{s}}^{-1}$. ThusThe useful density range lies between and around these scales. In the independent-upper-level approximation,The density sensitivity of a metastable doublet depends on the error measure. The fractional response $|d\ln R/d\ln n_e|$ is maximal atIf absolute ratio errors are constant, maximizing $|dR/d\ln n_e|$ instead gives $n_e=n_{c,2}$. An unspecified experimental noise model therefore does not define a unique sharper optimum than this transition-density range.

For Ne IV, the same calculation gives critical densities $4.23\times 10^{10}$ and $3.93\times 10^{11}{\mathrm{m}}^{-3}$, with fractional sensitivity peaking near $1.29\times 10^{11}{\mathrm{m}}^{-3}$. This is a useful higher-density diagnostic when O II is saturated. However, Ne${}^{3+}$ and O${}^{+}$ need not occupy the same ionization zone, so their inferred densities need not describe the same gas. The ultraviolet Ne IV pair also needs suitable observations and extinction correction. The proposal is useful with these qualifications, not a guaranteed replacement for the density of the O II emitting zone.

The main evidence combines a large luminosity with a compact source and a massive central gravitating object. Rapid continuum variability bounds the emitting size by a light-crossing scale, while broad emission lines and reverberation delays constrain the gas orbital speeds and characteristic radii. Their virial masses, together with stellar and gas dynamics in galactic nuclei, imply compact central masses far larger than ordinary stars. Thermal optical/ultraviolet disk emission, a high-energy coronal component, jets and relativistically broadened inner-disk features are naturally associated with an accretion flow. No one signature alone is a unique proof, but the combined compactness, mass and energy requirements strongly favor accretion onto supermassive black holes. Gravitational accretion releases a much larger fraction of rest energy than ordinary stellar burning, without requiring the luminosity to come from the black-hole interior.

For a fully ionized pure-hydrogen plasma, each electron receiving Thomson radiation force is coupled electrostatically to a proton carrying almost all the mass. With spherical isotropic luminosity $L$, the outward force per electron is ${\sigma }_{T}L/(4\pi r^{2}c)$, while gravity per electron–proton pair is approximately $GM{m}_p/r^2$. Equality defines the Eddington luminosity:Here $G$ is the gravitational constant, $m_p$ the proton mass, $c$ the speed of light, $M$ the central mass and ${\sigma }_T$ the Thomson cross-section. The estimate assumes electron-scattering opacity, tight electron–ion coupling and no other dominant force. Geometry, beaming, other opacities and photon trapping can change the relevant luminosity limit.

If ${\dot{M}}_0$ is the supplied rest-mass rate and a fraction $ϵ$ of its rest energy is radiated, $L=ϵ{\dot{M}}_0{c}^2$. The associated Eddington accretion rate is consequentlyThis is a supply rate. The black-hole growth rate is ${\dot{M}}_{\mathrm{BH}}=(1-ϵ){\dot{M}}_0$ if the remaining supplied energy increases the hole's mass. Confusing these rates loses a factor relevant to the growth-time calculation.

The order-of-magnitude efficiency $ϵ\sim 0.1$ is motivated by orbital binding energy released near the innermost stable orbit. A Newtonian orbit at $6GM/c^2$ has binding energy per unit mass $c^2/12\simeq 0.083c^2$; the exact nonspinning thin-disk value is $1-\sqrt{8/9}\simeq 0.057$. Black-hole spin changes the innermost stable orbit and can raise or lower the efficiency. Thus $0.1$ is a useful radiatively efficient accretion estimate, not a universal constant; an advective or jet-dominated flow can have a different radiative efficiency.

Using the stipulated luminosity and ${\sigma }_T=6.65\times 10^{-29}{\mathrm{m}}^2$,The source's redshift and age are treated as its supplied data; its record wording is historical, not a present-day claim about the most distant known quasar.

At continuous Eddington luminosity and fixed efficiency, growth is exponential:This is the Salpeter time, including retained-mass efficiency. Even granting all $0.77\mathrm{Gyr}$ since the big bang permits only about 15.4 e-folds, requiring a seed of at least roughly $400M_{\odot }$. A $100M_{\odot }$ seed would take about $0.84\mathrm{Gyr}$ under these assumptions, before allowing any formation delay; a $10^5{M}_{\odot }$ seed takes about $0.49\mathrm{Gyr}$. These numbers motivate early massive seeds, sustained high duty-cycle accretion, lower radiative efficiencies, super-Eddington phases or contributions from mergers. They do not prove a unique seed channel, since the luminosity-to-mass conversion and accretion history are model-dependent.

An additional atomic input is needed: the total Case B recombination coefficient does not specify the number of H-alpha photons per net radiative recombination. Let $y_{\alpha }={\alpha }_{\mathrm{H}\alpha }^{\mathrm{eff}}/{\alpha }_r$, where ${\alpha }_{\mathrm{H}\alpha }^{\mathrm{eff}}$ is the effective recombination coefficient. At the approximately $10^4\mathrm{K}$ temperature implicit in the supplied ${\alpha }_r$, take $y_{\alpha }\simeq 0.45$. The coefficient 2.2 used later in the paper corresponds to $y_{\alpha }=1/2.2\simeq 0.455$, which we use consistently here.

For an unextinguished, fully filled Strömgren sphere with negligible central cavity, the H-alpha luminosity and observed radiative flux satisfySubstitution eliminates $r_S$ and givesConvert the angular diameter to radians: $\theta =56\times 10^{-3}/206265=2.7150\times 10^{-7}$. Also $hc/{\lambda }_{\mathrm{H}\alpha }=3.0268\times 10^{-19}\mathrm{J}$. ThusThe implied hydrogen-ionizing photon production rate is $Q\simeq 1.00\times 10^{39}{\mathrm{s}}^{-1}$. Using $y_{\alpha }=0.45$ instead gives $D\simeq 505\mathrm{pc}$. These values use the printed milliarcseconds and integrated line observed radiative flux. Assuming one H-alpha photon per recombination would give a different, unjustified distance. An appreciable cavity or a nonunit filling factor must also enter the emitting volume.

The Brackett series consists of transitions ending at $n=4$, at infrared wavelengths. They suffer much less interstellar extinction than H-alpha, so their corrected observed radiative fluxes can give a more reliable version of the same inference. At fixed measured angular size and number density, the formula has $D\propto F$: uncorrected interstellar extinction therefore makes this particular distance estimate too small. Infrared lines are weaker and need their own effective recombination coefficients; their advantage is reduced sensitivity to dust.

Let $q_{\mathrm{in}}$ be the photoionization rate per neutral hydrogen atom just outside the empty cavity. With an isotropic central source and no attenuation inside that cavity,The cross section is zero below threshold; this permits the lower integration limit zero. Pure hydrogen with hydrogen ionization fraction $\zeta$ has $n_{\mathrm{HI}}=n_H(1-\zeta )$ and $n_e=n_p=n_H\zeta$. Local photoionization equilibrium is thereforeDefine the dimensionless quantity ${\Theta }_{\mathrm{phys}}={\alpha }_r{n}_H/q_{\mathrm{in}}$. Solving the resulting quadratic equation, ${\Theta }_{\mathrm{phys}}{\zeta }_{\mathrm{in}}^2+{\zeta }_{\mathrm{in}}-1=0$, gives the physical rootThe rationalized form avoids numerical cancellation for small ${\Theta }_{\mathrm{phys}}$. Strong illumination gives ${\zeta }_{\mathrm{in}}\simeq 1-{\Theta }_{\mathrm{phys}}$; weak illumination gives ${\zeta }_{\mathrm{in}}\simeq {\Theta }_{\mathrm{phys}}^{-1/2}$.

Eliminating $Q$ with the H-alpha luminosity, $Q=2.2(4\pi D^2{F}_{\mathrm{H}\alpha }){\lambda }_{\mathrm{H}\alpha }/(hc)$, yieldsThe printed definition of $\Theta$ is inconsistent with the defined radii. The source calls $r_{\ast }$ a physical stellar radius, but its factor $(2\pi r_{\ast }D{)}^2$ gives $\Theta$ units of length to the fourth power, rather than a dimensionless value. It also omits the physically necessary cavity radius $r_{\mathrm{in}}$. For a concrete counterexample, double $r_{\mathrm{in}}$ while holding the source, $D$, and $n_H$ fixed: the incident photoionization rate falls by four, so ${\Theta }_{\mathrm{phys}}$ rises by four; the printed expression does not change. The root has the requested form with the corrected ${\Theta }_{\mathrm{phys}}$ above. There is no supplied relation allowing $r_{\mathrm{in}}$ to be replaced by $r_{\ast }$.

TODO: assign the physical differences to the actual two panels once the missing evolved-star sequences are recovered.

The main intrinsic parameters for old globular clusters are age, metallicity, helium abundance and the distribution of envelope masses after red giant branch mass loss. Distance and foreground interstellar extinction alter the appearance and must first be removed before assigning intrinsic differences.

Age controls the main-sequence turnoff: a more massive star has a shorter main-sequence lifetime, approximately fuel divided by luminosity. An older coeval population therefore has a smaller turnoff mass and a fainter, generally redder main-sequence turnoff. The separation in astronomical magnitude between the main-sequence turnoff and horizontal branch is useful because a common distance modulus cancels.

At comparable age, larger metallicity increases stellar-envelope opacity and atmospheric line blanketing. This generally makes the red giant branch and main sequence redder; metal-poor populations tend to have hotter, bluer red giant branches. Metallicity also affects hydrogen burning through the CNO cycle, so a quantitative age inference needs composition-dependent stellar models.

The horizontal branch is especially sensitive to envelope mass: at similar helium-core mass, a thinner hydrogen envelope gives a higher effective temperature and hence a bluer location. Red giant branch mass loss, age and helium abundance all influence that envelope and produce the familiar “second parameter” complexity. Enhanced helium also changes mean molecular weight, luminosity and evolutionary lifetime, and can shift main sequences and horizontal branches.

Interstellar extinction dims the diagram, while wavelength-dependent interstellar extinction reddens $B-V$; distance primarily supplies a vertical distance modulus shift. A systematic bend or different horizontal branch colour distribution cannot be explained by distance alone. Because the source hides the evolved branches, it is not possible to assign “older”, “more metal rich”, or “bluer horizontal branch” to a particular panel reliably. The preceding mechanisms give the warranted physical interpretation without inventing the lost observational differences.

On the main sequence, central hydrogen burning converts four protons into helium, positrons, neutrinos and released energy. The net reaction is $4p\to {}^4\mathrm{He}+2e^{+}+2{\nu }_e$, with positron annihilation contributing to the total roughly $26.7\mathrm{MeV}$ release before subtracting escaping neutrino energy. In the low-mass main-sequence stars of an old globular cluster, the proton–proton chain generally dominates. The CNO cycle makes the same net conversion and becomes important in hotter burning regions.

After central hydrogen exhaustion, a subgiant contracts its helium core and establishes a hydrogen-burning shell. On the red giant branch, that shell supplies most of the luminosity, commonly through the temperature-sensitive CNO cycle, while the inert helium core grows and becomes supported mainly by electron degeneracy pressure. The envelope expands and its surface cools. The helium core is not yet powered by helium nuclear fusion.

Helium ignition in a sufficiently low-mass degenerate core produces the helium flash, removing the strong electron degeneracy pressure dominance. A horizontal-branch star subsequently burns helium in its core through the Triple-alpha process, $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}+\gamma$, and through ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$. A surrounding hydrogen-burning shell can continue simultaneously.

After central helium exhaustion, the Asymptotic giant branch has an inert carbon-oxygen core with a helium-burning shell and an outer hydrogen-burning shell. Alternation between quiescent burning and helium-shell thermal pulses drives mixing and substantial mass loss. Subsequent post-asymptotic giant branch evolution exposes a hot core that may ionize its expelled envelope as a planetary nebula. Finally a white dwarf cools mainly by loss of stored thermal energy, rather than sustained central nuclear fusion.

A dense cloud core first collapses when self-gravity overwhelms its effective support. While radiation escapes readily, compression can remain approximately isothermal and the Jeans mass decreases as density rises, encouraging gravitational fragmentation. Angular momentum usually leads to a disk rather than unrestricted spherical infall.

As the central optical depth increases, cooling becomes slower than compression and pressure produces a first hydrostatic core. molecular hydrogen dissociation near $2000\mathrm{K}$ absorbs compressional energy and lowers the effective adiabatic response, permitting a second collapse of a protostar. After dissociation the gas stiffens and a much more compact stellar core forms. An accretion shock and surrounding disk then convert infall energy into radiation. A small central core forms before most of the envelope arrives, as demonstrated in radiation-hydrodynamic collapse calculations.^[ref]

The isothermal inside-out collapse of a singular isothermal sphere model provides one idealized supply law, $\dot{M}\sim a_T^3/G$, but rotation, magnetic fields and an unsteady reservoir modify real accretion. Its luminosity scale is $L_{\mathrm{acc}}\sim GM\dot{M}/r_{\ast }$. Dust can absorb the emerging short-wavelength radiation and re-emit it in the infrared, so an embedded source need not have a directly visible stellar photosphere.

As envelope accretion declines, gravitational contraction becomes important on the Kelvin-Helmholtz cooling time $t_{\mathrm{KH}}\sim G{M}^2/(r_{\ast }L)$. deuterium-burning energy-generation rate can temporarily influence the radius; sustained core hydrogen burning eventually establishes the main sequence if the final mass is sufficient. A low-mass convective object approaches along a Hayashi-like track, while more massive objects need not share that track. accretion, contraction and nuclear burning are distinct energy sources, and a protostar is not generally a fixed-mass object.

mass and energy signs. The enclosed mass is $M(r)={\int }_0^{r}4\pi s^2\rho (s)ds$, so $dM/dr=+4\pi r^2\rho$. The negative sign printed in that equation cannot describe positive enclosed mass. Use the physical sign and define gravitational potential energy by

Boundary-pressure virial theorem. Multiply hydrostatic equilibrium by $4\pi r^3$ and integrate. The left side is $[4\pi r^{3}p{]}_0^{r_0}-3{\int }_0^{r_0}4\pi r^{2}pdr$, while the right side is $\Omega$. The central boundary contribution vanishes when $r^{3}p\to 0$ and the outer pressure is $p(r_0)=p_0$. ThusA finite external pressure is compatible with a gravitationally self-contained cloud, but is not the zero-surface-force notion of isolation. External ram stress or mass flow would require additional terms.

radiation budget of monatomic cloud collapse. For an ideal gas with fixed $\gamma$, the specific internal energy is $u=p/[\rho (\gamma -1)]$. Therefore the pressure integral is $3(\gamma -1)U$. Cold monatomic hydrogen has $\gamma =5/3$ when only translational degrees of freedom matter. With negligible final surface pressure, $2U+\Omega =0$, so $E_f=U+\Omega =\Omega /2$. If the initial cloud is sufficiently extended and cold that $E_i\simeq 0$, the emitted energy isBeing initially at rest alone does not set its binding energy to zero. If the initial total energy is appreciable, retain $E_{\mathrm{rad}}=E_i-\Omega /2$; a different equation of state also changes the half-energy result.

Inside-out accretion. The combination of $a_T$ and $G$ with dimensions of mass per time iswhere dimensional analysis leaves the dimensionless coefficient $C$ undetermined. An isothermal gas keeps $a_T$ constant, giving constant supply in the self-similar idealization. A front at radius $a_{T}t$ and an isothermal mass scale $M(r)\sim a_T^{2}r/G$ give the same scaling. The front's swept mass is not automatically the accreted central mass; the coefficient requires a dynamical similarity solution.^[ref]

constant-radius accretion energy budget. If all available accretion power is radiated at the fixed radius, $L_{\mathrm{acc}}=GM\dot{M}/r_0$. Neglecting the seed mass,This equals $-\Omega /2$ precisely whenThe positive sign in the final printed energy condition is inconsistent with the earlier use of negative gravitational potential energy. Equivalently its positive right side must denote the binding-energy magnitude $|\Omega |$. A general structural coefficient $\Omega =-qG{M}^2/r_0$ leaves a mismatch unless $q=1$ under these other idealizations. Radiation trapped in the object, outflows or release over a different radius also change the budget.

Which temperature? With $M=M_s+\dot{M}t$ and fixed $r_0$, $L_{\mathrm{acc}}\propto M$. If this power emerges thermally from area $4\pi r_0^2$, the accretion-powered effective-temperature scaling iswhen the seed mass is negligible. This is the radiating surface's effective temperature, inferred from the preceding accretion-luminosity calculation. A homologous ideal-gas interior instead has virial temperature $T_c\sim \mu m_{H}GM/(k_B{r}_0)\propto t$. Thus the surface-effective and interior-virial scalings are $t^{1/4}$ and $t$, respectively, under their stated assumptions. Inferring a detailed interior temperature profile would additionally require a structural and transport model.

virial pressure-radius curve with uniform-sphere binding energy. Put $c_s^2=k_{B}T/\mu$. Isothermality makes $\int (p/\rho )dM=c_s^2{M}_0$, so the virial theorem givesThe curve crosses zero at $r_z=B/A=G{M}_0/(5c_s^2)$, rises to its maximum atthen decreases toward zero at large radius. Since $p_0^{\prime }=(4B-3A{r}_0)/r_0^5$, a cloud displaced outward at fixed external pressure has a restoring pressure deficit only on the descending branch. Thus the homologous, isothermal radial mode is stable for $r_0>r_c$, unstable for $r_z<r_0<r_c$, and marginal at $r_c$. The negative-pressure branch is not a physical equilibrium with a positive confining pressure. This conclusion uses the stipulated gravitational-energy approximation, not the exact isothermal density profile.

critical mass from a virial pressure turning point. At the turning point, $B=3A{r}_c/4$; henceUsing the cloud's mean density $\bar{\rho }=3M_0/(4\pi r_0^3)$, the condition $r_0>4G{M}_0/(15c_s^2)$ becomesUsing the boundary density instead changes the numerical coefficient, not this scaling. The coefficient should not be identified with the exact Bonnor-Ebert sphere threshold because the assumed potential energy is different.

mass and luminosity weighting of a stellar initial mass function. Impose finite birth-mass cutoffs $m_l,m_u$. The stellar birth mass and main-sequence luminosity scale asFor a broad mass interval, the low-mass end controls the first integral and the high-mass end the second. Per logarithmic mass interval the weights are respectively $m^2\xi \propto m^{-0.35}$ and $mL\xi \propto m^{2.65}$. Finite lifetimes remove massive stars from an aging population, so the luminosity conclusion applies while those stars remain on the main sequence.

exponential closed-box gas consumption. The gross formation rate is $\Psi =ϵM_g$, but the permanent locking rate is $\alpha \Psi$. The instantaneous recycling approximation and mass conservation give ${\dot{M}}_g=-\alpha ϵM_g$, ${\dot{M}}_{\ast }=\alpha ϵM_g$. ThereforeHere $M_{\ast }$ means surviving locked stellar mass in the galactic gas fraction, not the cumulative gross mass that has passed through star formation.

locked-mass and returned-mass yield conventions. Pre-existing metals are locked at rate $\alpha Z\Psi$; the stipulated newly produced metal source is $b(1-\alpha )\Psi$. ThusFor $\alpha >0$ and initially pristine gas,The effective yield $y$ is per unit permanently locked mass; it is not the source parameter $b$. For $\alpha =0$, $M_g=M_0$ and the direct time equation gives $Z=bϵt$. For $\alpha =1$, the explicitly stipulated metal source vanishes, so this particular enrichment model gives $Z=0$. Constant yield is a trace-enrichment approximation; the formal logarithmic divergence as gas is exhausted must not be interpreted as a physical metal fraction exceeding one.

backward extrapolation of exponential closed-box star formation. Denote the present age of the universe by $t_{\mathrm{now}}$. With the last paragraph's $\alpha =1$,Adopting a present age of the universe $t_{\mathrm{now}}=13.8\mathrm{Gyr}$ for the estimate gives about $365M_{\odot }\mathrm{y}{\mathrm{r}}^{-1}$, requiring $7.3\times 10^{11}{M}_{\odot }$ of gas at that epoch in the same model. The boxed expression makes the dependence on the adopted present age explicit. This is an extreme history for a present Milky-Way-like system, although intense high-redshift star formation by itself is not forbidden.

The extrapolation forces a fixed gas reservoir to decline exponentially over many gas-depletion times. Gas accretion, outflows, evolving gas density and star-formation efficiency, delayed recycling, and galaxy assembly break that assumption. The empirical disk-averaged relation is not a universal constant-efficiency law for all gas: the original measurements combined Balmer series H-alpha, atomic gas, carbon monoxide and infrared data, with tracer calibration and dust corrections entering the inference.^[ref] Cosmological calculations explicitly support continuing cold-gas supply to rapidly forming high-redshift galaxies, providing a physical alternative to this closed-box backward extrapolation.^[ref]