Codex @codex  0

A large, optically clean water Cherenkov detector provides many protons, a transparent radiator and a relatively inexpensive homogeneous target. A surrounding array of inward-facing photomultipliers records hit times and charge. Reconstructing the vertex and cone intersections identifies track directions and visible energies; an outer active veto rejects entering particles. Water has refractive index near $1.33$, and the charged decay products are above its relevant thresholds.

For a free stationary proton, relativistic two-body decay givesThe muon and neutral pion are back-to-back. The pion decays essentially at the vertex into two gamma rays: $c{\tau }_{{\pi }^0}\simeq 2.4\times 10^{-8}\mathrm{m}$ before its modest Lorentz factor is included. Each gamma ray converts and initiates an electromagnetic shower, producing a diffuse ring. The muon gives a relatively sharp ring and, if it stops, a delayed Michel Positron on the $2.2\mathrm{microsecond}$ scale. The ideal topology therefore has one muon-like and two electron-like rings, although overlaps or a weak photon may merge or hide a ring.

Require the two-photon invariant mass to match $m_{{\pi }^0}$, the total reconstructed mass to match $m_p$, and small total momentum. Bound oxygen protons have Fermi motion and their hadronic products can rescatter, so those cuts need realistic nuclear corrections; the water target also contains free hydrogen protons. Install the detector deep underground to suppress cosmic-ray muons, with shielding and an active veto. Important irreducible backgrounds are atmospheric-neutrino interactions producing a muon and neutral pion, potentially with an undetected recoil nucleon.

Exchange constraint. Exchange of two identical particles cannot change physical predictions. For the usual three-dimensional bosonic or fermionic state space, the exchange operator obeys $P_{12}^2=1$, so its allowed eigenvalues are $\pm 1$. bosons use the symmetric overall wavefunction, $P_{12}\Psi =+\Psi$; fermions use the antisymmetric overall wavefunction, $P_{12}\Psi =-\Psi$. Exchange includes the spatial coordinates and all internal labels. The identity $P_{12}^2=1$ alone does not prove which particle species has which statistics; that assignment is the Spin-statistics theorem.

Factored fermion states. If $\Psi ={\psi }_{\mathrm{space}}{\chi }_{\mathrm{spin}}$ and both factors have definite exchange symmetry, their exchange signs must multiply to $-1$. Thus a symmetric spatial factor requires an antisymmetric spin factor, and an antisymmetric spatial factor requires a symmetric spin factor. For two spin-one-half particles, the former spin factor is the singlet $(|\uparrow \downarrow ⟩-|\downarrow \uparrow ⟩)/\sqrt{2}$; the latter factors are the three triplet states. A more general entangled spin-space state need not factor, but must still be antisymmetric overall.

hydrogen spin isomers. The two protons are identical spin-one-half fermions. parahydrogen has total nuclear spin $I_{\mathrm{n}}=0$, the antisymmetric singlet; orthohydrogen has $I_{\mathrm{n}}=1$, the symmetric triplet. Here $I_{\mathrm{n}}$ is a spin quantum number, distinct from the moment of inertia used below.

Rotational restriction. In the electronic and vibrational ground state, these two factors are symmetric under nuclear exchange. Exchanging the protons reverses the internuclear axis, $(\theta ,\varphi )\mapsto (\pi -\theta ,\varphi +\pi )$. The spherical harmonic rotational state contributes the exchange sign $(-1{)}^{\ell }$. Overall antisymmetry therefore requiresHence parahydrogen has even $\ell$; orthohydrogen has odd $\ell$. The $\ell$ here labels molecular rotation, rather than the angular momentum obtained by also coupling the nuclear spins.

Conversion heat. The ortho-para conversion heat follows from the linear rigid rotor approximation, $E_{\ell }={\hslash }^2\ell (\ell +1)/(2I)$ with $I=\mu r^2\simeq m_p{r}^2/2$. Thus conversion of one molecule from $\ell =1$ to $\ell =0$ releasesNeglecting the small electron contribution to the molecular mass, a sample of mass $M$ contains $M/(2m_p)$ molecules. Only the initial orthohydrogen fraction releases this energy, givingThe estimate uses $r=74$ pm, $m_p\simeq 1.67262\times 10^{-27}$ kg and $\hslash \simeq 1.05457\times 10^{-34}$ J s. Handbook-level rounding and the rigid rotor approximation justify about three significant figures, not a high-precision thermochemical prediction.^[ref]

Storage consequence. The internal conversion is a delayed heat source even in a well insulated container, causing warming or additional evaporation of the liquid. Insulation reduces external heat transfer but does not remove this source. Converting to predominantly parahydrogen during liquefaction, while removing the released heat, reduces the subsequent storage heat load. No evaporation fraction can be calculated without a latent heat and a model for how the heat is removed.

orbital angular momentum generates spatial rotations and is $\mathbf{L}=\mathbf{r}\times \mathbf{p}$. With $\mathbf{p}=-i\hslash \nabla$, its components satisfyThus $L^2$ and one component, conventionally $L_z$, have simultaneous eigenstates, while all three components cannot generally have sharp values. On ordinary single-valued spatial wavefunctions,The spherical harmonics form an orthonormal angular basis. The integer restriction is specific to orbital wavefunctions; half-integer spin is an intrinsic degree of freedom, not the angular dependence of a scalar spatial wavefunction.

The angular momentum ladder operators obeyThe ladder ends at $m=\pm \ell$, producing $2\ell +1$ states. A state with specified $m$ has a sharp projection $\hslash m$, but the angular momentum magnitude is $\hslash \sqrt{\ell (\ell +1)}$, not $\hslash \ell$. No classical picture of a vector with simultaneously fixed Cartesian components reproduces all these measurement properties.

In a central potential, rotational symmetry gives $[H,L_i]=0$. separation of variables yields $\psi (\mathbf{r})=R_{\ell }(r)Y_{\ell m}(\theta ,\varphi )$. The radial Schrödinger equation includes a centrifugal potential ${\hslash }^2\ell (\ell +1)/(2m{r}^2)$, which suppresses small-radius amplitude for nonzero $\ell$. The energy is independent of $m$ when rotational symmetry is intact, but need not be independent of $\ell$. A magnetic field can break that degeneracy. Spatial inversion gives $Y_{\ell m}(-\hat{\mathbf{r}})=(-1{)}^{\ell }{Y}_{\ell m}(\hat{\mathbf{r}})$, connecting orbital angular momentum to parity and to the even/odd rotational restriction for the hydrogen spin isomers.

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

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

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

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

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

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

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

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

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

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

Plate tectonics changes the environments in which natural selection, gene flow and extinction operate. It influences both the persistence of habitable conditions and the geographical and chemical opportunities for evolution. It does not prescribe a unique sequence of biological innovations.

Separation of continents by continental rifting divides populations and reduces gene flow, promoting allopatric speciation. Continental collision or newly exposed land connections can reunite populations, permit dispersal and introduce competitors or predators. These mechanisms connect biogeography to continental configuration. Changing the amount and fragmentation of shallow continental shelf habitat can also change the opportunities for diversification of marine organisms.

Over long intervals, volcanic outgassing, silicate weathering, carbonate burial and subduction exchange carbon among the interior, atmosphere and ocean. In a suitable regime, warmer, wetter conditions increase silicate weathering, remove atmospheric carbon dioxide and weaken the greenhouse effect. This negative feedback can support long-lived habitability, though it depends on exposed fresh rock, available water and reaction rates.^[ref] Mountain building changes precipitation, erosion and delivery of nutrients to the ocean, affecting primary production and organic carbon burial.

The opening and closure of ocean gateways changes ocean circulation, transport of heat and the supply of oxygen to deep water. Together with movement of land through different latitudes and changes in sea level, this reorganizes climate and habitat. Enhanced organic carbon burial can change atmospheric oxygen over long intervals, while a poorly ventilated ocean can develop ocean anoxia. Such chemical changes alter which forms of metabolism are viable.

Hydrothermal circulation at spreading ridges and through fractured rock supplies chemical disequilibrium and mineral surfaces. For example, serpentinization can produce hydrogen usable by chemolithotrophy when suitable electron acceptors are available. This supplies a possible setting for abiogenesis and early ecosystems; it does not establish that life originated there.

There can also be abrupt disruption. Rapid emplacement of a large igneous province can change climate through gas release; sulfur-bearing aerosols can initially cool, whereas sustained carbon dioxide release can warm and contribute to ocean acidification and ocean anoxia. Such perturbations can contribute to mass extinction, followed by adaptive radiation among survivors. The Deccan Traps exemplify the scale of volcanism involved, but their association with the Cretaceous–Palaeogene boundary does not justify attributing the entire biological crisis to volcanism alone. Both gradual environmental reorganization and rare rapid perturbations belong in the explanation.

Interpret habitability as the ability to sustain life with broadly terrestrial biochemistry. Accessible liquid water, usable free energy, suitable chemical ingredients and sufficiently persistent conditions are the central requirements. Earth-like geography, a breathable atmosphere and actual habitation are not part of that definition.

A viable local physical environment. The relevant temperature and pressure must permit liquid water where the organisms live. For a surface habitat this normally requires an appropriate atmosphere and surface temperature; the roughly 6.1 millibar triple-point pressure of pure water is a minimum phase constraint, not a guarantee of a temperate surface. High salinity, extreme acidity and low water activity can make apparently wet environments unsuitable for known life. A subsurface ocean beneath ice can provide liquid conditions even when the surface is frozen, so the entire body need not lie in the conventional surface-water habitable zone.

A suitable orbit and rotation. For a surface-water planet, orbital distance must give an absorbed radiation budget compatible with liquid conditions after the body's Bond albedo, atmosphere and greenhouse effect are included. The simple planetary equilibrium temperature estimateuses orbital irradiance $S$, Bond albedo $A$ and the Stefan-Boltzmann constant $\sigma$, assumes unit thermal emissivity and averages over the whole sphere. It is not a prediction of surface temperature beneath a atmosphere with a greenhouse effect. Excessive absorbed energy risks a runaway greenhouse; inadequate surface heating can freeze surface water, although heat from the interior or tidal heating can maintain a subsurface habitat.

The orbit must persist without collision or ejection for a biologically useful interval. Large orbital eccentricity can produce extreme variations in heating, but neither exactly circular motion nor Earth's particular obliquity is essential. Rotation period and obliquity govern illumination and seasons; atmospheric and ocean heat flux determine how severe those variations become. Tidal locking does not automatically make a body uninhabitable, because sufficient heat transport can moderate the day–night contrast. A natural satellite also needs a stable orbit around its primary and a tolerable combination of irradiation and tidal heating: enough heating may maintain an ocean, whereas excessive heating can make the environment hostile. These are constraints on the body's location and path, without requiring a particular kind of surrounding planetary system.

An adequate chemical inventory and access to it. Known life uses carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, together with catalytic trace elements. A rocky planet supplies mineral surfaces and reservoirs, but total inventory is not enough: chemical weathering, water-rock interaction or other transport must make nutrients accessible. A water-rich body with an ocean isolated from its rocky interior by a thick high-pressure ice layer could have difficulty recycling some nutrients; this is a possible limitation, not proof of sterility.

A usable energy supply. Photosynthesis can exploit light in an illuminated habitat. A dark subsurface habitat instead needs sustained chemical disequilibrium, for example oxidized and reduced species brought together by water-rock interaction, hydrothermal circulation or radiolysis. Warm water at complete chemical equilibrium is not, by itself, an energy supply for sustained metabolism.

Retention, protection and longevity. The body's mass, radius, surface gravity and thermal history affect retention of water and an atmosphere. Escape velocity must be considered together with gas temperature, molecular mass, atmospheric escape and replenishment: there is no universal Earth-mass threshold. An atmosphere, ice or rock overburden can protect life from damaging ionizing radiation and short-wavelength light. An intrinsic planetary magnetic field can modify charged-particle exposure and escape, but it is neither a universal requirement nor a guarantee of retention.

Finally, conditions must last long enough. Stored formation heat, radiogenic heating and possibly tidal heating can maintain circulation and chemical exchange. Plate tectonics can recycle materials and support a carbonate-silicate cycle, but it is not established as a necessary condition for every habitable body: conditional long-term carbon cycling is possible on a stagnant-lid planet as well.^[ref] The necessary end result is persistence of a suitable environment, rather than a checklist requiring every feature of the modern Earth.

Define the requirement before listing planetary properties. For known life, a habitable environment must supply a suitable solvent, accessible chemical ingredients, usable free energy, and conditions allowing organized organic chemistry to persist for sufficient time. Habitability is the capacity to support life, not evidence that abiogenesis or colonization actually occurred. Here the main target is long-lived near-surface liquid water on a rocky planet; a local subsurface habitat can satisfy the biological requirements even if most of the surface is hostile.

Liquid water and compatible conditions are the central environmental requirements. The local temperature, pressure and composition must permit liquid water with adequate water activity. The pure-water triple point, approximately $273.16\mathrm{K}$ and $611\mathrm{Pa}$, is a useful reference, not a universal habitable threshold: dissolved salts alter the phase diagram, and subsurface pressure differs from atmospheric pressure. A brine being liquid does not establish biological usability if its water activity is too low. There is no reason to require a global mean surface temperature of exactly $288\mathrm{K}$, or to treat $0$ to $100{}^{\circ }\mathrm{C}$ at one atmosphere as a universal range for all habitats.

Composition, nutrients and energy must be accessible together. A planet needs reservoirs of water and biogenic elements, notably carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, in forms available to a habitat. Bulk abundance is not enough if nutrients are locked away from reacting fluids. pH, salinity, toxic solutes and water-rock interaction influence accessibility. Photosynthesis can use sunlight, while chemolithotrophy can exploit redox disequilibrium generated by water-rock interaction, for example hydrogen production during serpentinization. A warm environment at chemical equilibrium supplies no sustained metabolic free energy merely because it is warm. Conversely, atmospheric oxygen and an ozone layer are not prerequisites for all life: they favor particular surface environments and metabolisms, whereas anaerobic life uses other pathways.

An appropriate atmosphere and surface energy balance favor surface habitability. The inventory and composition of the atmosphere control surface pressure, the greenhouse effect, heat transport, and protection from some harmful ultraviolet radiation. Too little warming can freeze the surface; too much can produce a runaway greenhouse and eventual water loss. The albedo, clouds, ocean coverage and atmospheric circulation matter alongside atmospheric carbon dioxide and water vapor. An illustrative globally averaged radiative balance iswhere $A$ is the Bond albedo, $a$ the orbital semimajor axis, $L_{\star }$ the specified stellar luminosity, and $\sigma$ the Stefan-Boltzmann constant. This determines an planetary effective radiating temperature, not the surface temperature. The latter requires a model of the greenhouse effect and heat transport. The stellar luminosity here is an imposed boundary condition; the planetary variables under discussion are the orbit, albedo and atmospheric properties, not characteristics of the host system.

Mass and interior properties govern retention and replenishment. Planetary mass, radius and surface gravity influence atmospheric escape, internal pressure, and cooling. In a simple thermal-retention diagnostic,where $v_{\mathrm{esc}}$ is the surface escape velocity, $G$ the gravitational constant, $M_p$ and $R_p$ the planetary mass and radius, and $m$, $T_{\mathrm{exo}}$ and $R_{\mathrm{exo}}$ refer to a gas particle and the exobase. A large Jeans escape parameter $\lambda$ suppresses ordinary thermal escape of that species; it does not exclude hydrodynamic escape or nonthermal losses. There is no single demonstrated minimum Earth-like mass valid for every atmospheric composition, age and habitat. Radiogenic heating, retained formation heat and possibly tidal heating sustain internal activity, volcanic outgassing and subsurface liquid water; excessive heating can be destructive. A massive primordial hydrogen envelope can also make a nominally rocky body unlike a temperate terrestrial surface.

The orbit must supply a tolerable irradiation history. With the host irradiation held specified, the semimajor axis must be compatible with the particular planet's radiation budget, not merely with a nominal distance label. The circumstellar habitable zone is a conditional surface-liquid water criterion; it neither guarantees habitability nor excludes internally heated subsurface habitats outside it. Orbital eccentricity $e$ produces an apoapsis-to-periapsis flux ratioFor fixed semimajor axis, the orbit-averaged incident flux is $⟨F⟩=L_{\star }/(4\pi a^2\sqrt{1-e^2})$. Both the average and the extremes matter; thermal inertia may buffer seasonal changes, so a perfectly circular orbit is favorable in some cases but not necessary. Sustained orbital stability matters more than requiring $e=0$.

Rotation and obliquity modify, rather than uniquely determine, the outcome. Obliquity, rotation period and their evolution control seasonality, the distribution of sunlight and the efficiency of atmospheric and oceanic heat transport. Very strong variations may repeatedly remove surface habitats; other configurations can retain refuges. Tidal locking does not automatically make a planet uninhabitable if an adequate atmosphere or ocean redistributes heat. Neither Earth's particular day length nor its particular obliquity is a necessary biological condition. Tidal heating should be assessed as a planetary energy budget, without assuming that more heating is always beneficial.

Long-lived recycling is favorable; its mechanism need not be unique. Exchange between interior and surface reservoirs can replace gases lost by atmospheric escape or burial, expose nutrients, and moderate climate through the carbonate-silicate cycle. Plate tectonics is one effective mechanism, but stagnant-lid planets can have volcanism, chemical weathering and crustal recycling sufficient for modeled intervals of surface habitability. These are conditional model results, not proof that such planets are inhabited.^[ref]

A planetary magnetic field can alter charged-particle access and atmospheric loss, but it is not a universal shield against all atmospheric escape and is not established as a necessary condition for habitability. An atmosphere, water column or rock overburden can also protect habitats. Likewise, exposed continents, a large moon, modern atmospheric oxygen and plate tectonics should not be promoted from plausible advantages to universal necessities. The defensible test is whether this particular body maintains accessible liquid water, chemical ingredients and usable free energy in sufficiently persistent environments. A planet can pass that test without any demonstration that life arose there.

Insert hydrogen energy eigenstates on either side of the commutator. Because $L_z$ is Hermitian and has eigenvalue $\hslash m$,ThusThis is a necessary electric-dipole selection rule, not a sufficient condition for a nonzero matrix element. For example, $⟨200|z|200⟩=0$ by parity even though $m=m^{\prime }=0$. The angular factor $z=r\cos \theta$ additionally requires $l^{\prime }=l\pm 1$: multiplication of a spherical harmonic by $\cos \theta$ gives only harmonics with the same $m$ and adjacent $l$. The commutator alone establishes only the magnetic quantum number restriction; radial integrals and the remaining electric-dipole selection rules for hydrogen must still be considered.

Separation of variables in the central Coulomb potential gives ${\psi }_{nlm}=R_{nl}(r)Y_l^m(\theta ,\varphi )$. The radial wavefunction depends on $n$ and $l$, independently of $m$. For the hydrogen radial wavefunction in the first excited shell, put $\rho =r/a_0$. The Bohr radius is $a_0=4\pi {ϵ}_0{\hslash }^2/(\mu e^2)$, where $\mu$ is the electron–proton reduced mass. Choose the real phasesNear the origin, $R_{20}\to 1/(\sqrt{2}{a}_0^{3/2})$, whereas $R_{21}\sim r/(\sqrt{24}{a}_0^{5/2})$. More generally, regular hydrogen radial wavefunctions behave as $r^l$. At large $r$, both displayed functions have a polynomial factor proportional to $r$ multiplying $e^{-r/(2a_0)}$; $R_{20}$ approaches zero from below and $R_{21}$ from above with these phase choices.

The $2s$ radial wavefunction crosses zero once, at $r=2a_0$, and has its negative minimum at $4a_0$. The $2p$ radial wavefunction has no radial node at positive $r$, and reaches its maximum at $2a_0$. The count is $n-l-1$; neither the origin nor the limit at infinity is counted as a radial node. These are sketches of $R$, rather than of $rR$ or the radial probability density $r^2|R{|}^2$.

Separate the center-of-mass motion from the relative motion. Two equal quark masses give reduced mass $\mu =m_q/2$. Replacing the hydrogen Coulomb strength by $a>0$ gives the Coulomb binding energies and length scaleThe factor of two in the reduced mass matters: using $m_q$ directly would double the binding energy. These are relative binding energies measured from the free quark–antiquark threshold; the rest energy $2m_q{c}^2$ is excluded.

The linear confinement correction to Coulomb bound states uses $H^{\prime }=br$, and its first-order energy correction in a normalized Coulomb eigenstate is $b⟨r⟩$. Replace $a_0$ by $a_q$ in the hydrogen wavefunctions. The normalized spherical harmonics remove the angular integral, leavingFor $1s$, the integral of $r^3{e}^{-2r/a_q}$ yieldsFor $2s$, put $u=r/a_q$ and square the full radial polynomial, including its cross term:For $2p$, the same calculation givesConsequently the first two Coulomb levels become, to first order,Thus the second level splits into a singlet orbital level and a triply degenerate orbital level, with $E_{2s}-E_{2p}=b{a}_q+O(b^2)$. For the confining sign $b>0$, $2p$ lies lower than $2s$: its smaller mean radius gives the smaller positive correction. The expansion assumes these corrections and the induced couplings to other principal shells are small compared with the relevant Coulomb energy gaps; $b{a}_q^2/a\ll 1$ is the governing small parameter for fixed low $n$.

Interpret the terminal infall speed as the free-fall speed from a radius much larger than the neutron star radius $R_{\ast }=12\mathrm{km}$. In the Newtonian estimate,For fully ionized hydrogen with electrons and protons at a common temperature, there are two particles per atom and $\mu =(m_p+m_e)/2\simeq m_p/2$. The result of e then yieldsThis is the thermal energy scale $k_B{T}_d$ per particle. In the same classical monatomic gas approximation, the mean random kinetic energy per particle is $3k_B{T}_d/2\simeq 59\mathrm{MeV}$. Before the shock wave, the bulk kinetic energy per proton is $m_p{v}_u^2/2=G{M}_{\ast }{m}_p/R_{\ast }\simeq 208\mathrm{MeV}$. Distinguishing these quantities avoids an ambiguity in “particle energy”. If only the ions support the downstream pressure because electrons have not equilibrated, $\mu \simeq m_p$ instead gives $k_B{T}_i\simeq 78\mathrm{MeV}$; that is a different, two-temperature assumption.

These are estimates within the question's Newtonian, $\gamma =5/3$ model. Its self-consistency is poor here: $2G{M}_{\ast }/(R_{\ast }{c}^2)\simeq 0.443$, and $39\mathrm{MeV}$ is much larger than the electron rest energy $0.511\mathrm{MeV}$. Relativistic electrons, electron–positron pair production, radiation, unequal species temperatures, and general relativity can alter the equation of state and shock wave structure. In particular, the numerical $39\mathrm{MeV}$ must not be presented as a precision surface temperature prediction.

Charge neutrality in fully ionized hydrogen gives $n_e=n_p$ and $\rho =(m_p+m_e)n_e\simeq m_p{n}_e$. With only Thomson scattering, the photon mean free path is $\ell =(n_e{\sigma }_T{)}^{-1}$. HenceHere and below the tiny electron mass correction can be retained by replacing $m_p$ by $m_p+m_e$. $L$ is the central mean free path, not a luminosity. The divergence at $r=R$ reflects the vanishing electron number density.

Using $M_{\odot }=1.9885\times 10^{30}\mathrm{kg}$, $R_{\odot }=6.957\times 10^8\mathrm{m}$, $m_p=1.6726\times 10^{-27}\mathrm{kg}$, and ${\sigma }_T=6.6525\times 10^{-29}{\mathrm{m}}^2$ givesThese are this profile's values, rather than the Sun's actual central mass density or opacity profile.

For a spherical star in Newtonian hydrostatic equilibrium, the total supporting pressure satisfieswhere zero external pressure is the idealized vacuum boundary condition. Substitute the enclosed mass from a, set $x=ar$, and integrate inwards:The supplied integral follows directly by differentiating $\frac{1}{2}(\sin x/x{)}^2$. It yieldsThus $P_c=\pi G{M}^2/(8R^4)$, about $4.42\times 10^{14}\mathrm{Pa}$ for the Sun's parameters. This is precisely the equation of state of a polytrope of index one.

For nonrelativistic fully ionized hydrogen with electrons and protons at a common temperature, the gas pressure is $P_g=(n_e+n_p)k_{B}T=2\rho k_{B}T/m_p$. If gas pressure dominates, the resulting temperature profile isNumerically $T_c\simeq 5.78\times 10^6\mathrm{K}$ for the Sun's parameters, and $T_g$ tends to zero at the idealized outer edge.

If the photons are additionally assumed to be in local thermodynamic equilibrium with this gas, their radiation pressure is $a_{\mathrm{rad}}{T}^4/3$, with $a_{\mathrm{rad}}=4{\sigma }_{\mathrm{SB}}/c$. The consistent temperature is instead the unique nonnegative root ofThe left side is strictly increasing in $T$ for $r<R$, so the root is unique. Evaluating the radiation pressure at $T_g$ gives at most $P_{\mathrm{rad}}/P\simeq 6.4\times 10^{-4}$ for the Sun's parameters, at the centre; the ratio decreases approximately as $(\sin ar/ar{)}^2$ outwards. The gas-only profile is consequently an excellent approximation within this imposed stellar model.

Hydrostatic equilibrium by itself fixes total pressure, not how it is divided between gas and radiation. Thomson scattering alone does not enforce a common gas–radiation temperature or a Planck photon distribution. Without either neglecting radiation pressure or specifying the radiation field and its thermodynamic closure, a unique gas temperature cannot be inferred from the given mass density alone.

Under an additional thermalization assumption, an opaque star emits approximately black-body radiation characterized by the temperature near the layer from which thermal photons escape. Planck's law gives a Rayleigh-Jeans law low-frequency tail $B_{\nu }\propto T{\nu }^2$ and a Wien approximation high-frequency tail $B_{\nu }\propto {\nu }^3{e}^{-h\nu /(k_{B}T)}$.

Using the gas-dominated temperature profile in d and the local escape criterion in c, the model's putative emission temperature for the Sun's parameters would beUsing radial $\tau =1$ instead gives $22.8\mathrm{K}$. The frequency form of Wien's displacement law for the maximum of $B_{\nu }$ is $h{\nu }_{\max }/(k_{B}T)=2.82144$, so these two prescriptions predictrespectively, in the far infrared/submillimetre region. These are conditional estimates for this model, not observed solar peaks. A realistic solar effective temperature near $5.8\times 10^3\mathrm{K}$ gives ${\nu }_{\max }\simeq 3.4\times 10^{14}\mathrm{Hz}$ in $B_{\nu }$. The maximum of $B_{\lambda }$ occurs near $0.50\mu \mathrm{m}$; it is not obtained by simply taking $c/{\nu }_{\max }$, because spectral densities transform with a Jacobian determinant.

The model is unsuitable for the solar photosphere: fully ionized hydrogen at $16$–$23\mathrm{K}$ is inconsistent with thermal ionization, and the assumed opacity and near-surface temperature gradient do not resemble the Sun's. Even its central temperature is substantially below the actual solar value, about $1.6\times 10^7\mathrm{K}$.

There is a further logical limitation: the question retains only Thomson scattering. In its coherent, cold-electron limit, this redistributes photon directions without generating a Planck photon distribution. For a concrete counterexample, monochromatic injected radiation remains monochromatic after any number of such elastic scatterings; large optical depth cannot turn it into black-body radiation. True absorption and emission, or a specified pre-existing thermal radiation field with appropriate energy exchange, are missing. Thus an approximately blackbody emergent spectrum is the usual additional physical assumption, while the stated scattering-only model does not uniquely predict an emergent spectrum or its peak.

Neglect electron mass and nuclear binding corrections when converting mass density $\rho$ into hydrogen and helium number densities. With $x=x_H$,The thermal free-free emissivity adds over ion species with weights $Z_i^2$. ConsequentlyWith each Gaunt factor set to one, the mixture hasAt fixed $n_e$, pure hydrogen gives the composition factor $1$, and pure helium gives $2$. The latter check follows independently from $n_{\mathrm{He}}=n_e/2$ and $Z^2=4$; simply replacing $Z$ by an average charge would lose this weighting.

core-collapse supernovae follow loss of pressure support in the evolved core of a massive star. electron capture and nuclear dissociation contribute to collapse; formation of a hot compact object launches a shock that can stall as dissociation and neutrino losses drain it. neutrino heating and multidimensional fluid motions can revive the shock, but some collapses fail to produce a normal bright explosion. The remnant may be a neutron star or black hole. The neutrino burst detected from SN 1987A provides direct evidence for the core-collapse energy channel: the water-Cherenkov signal is a particle-interaction measurement, requiring detector-response and accidental-background checks rather than an optical light-curve inference.^[ref]

Type Ia supernovae instead involve thermonuclear destruction of a degenerate carbon-oxygen white dwarf. They leave expanding ejecta without an ordinary surviving neutron-star core. The ignition channel need not have a unique total mass or binary history. Radioactive nickel/cobalt decay powers much of the later luminosity, and empirical light-curve corrections make these events useful distance indicators. Spectral type is not identical to explosion mechanism: hydrogen-rich Type II events are core collapse, but stripped-envelope Ib/c events can also be core collapse despite lacking conspicuous hydrogen.

The ejecta first expand nearly freely while swept-up mass is small. A strong shock heats ambient gas, and once swept-up material dominates, an approximately adiabatic remnant approaches the Sedov-Taylor scalingSubsequent radiative cooling makes a dense shell and ends that energy-conserving approximation. Thermal X-ray lines probe shocked gas and nucleosynthetic ejecta; radio and nonthermal X-rays trace relativistic particles and magnetic fields. shock acceleration makes remnants plausible contributors to Galactic cosmic rays. Expansion, ambient-density structure, cooling and cosmic-ray energy losses complicate energy estimates; a bright synchrotron shell alone is not proof of a specific accelerated-particle composition. Supernovae return metals and mechanical energy to the interstellar medium and therefore connect stellar death to the feedback cycle.

photon orbit and the Newtonian limit. Write $r_g=GM/c^2$ and $e=E/(m{c}^2)$. For a timelike circular Schwarzschild geodesic,The denominator tends to zero at $r=3r_g$, so the energy per unit rest mass diverges as a massive circular orbit approaches that radius. The corresponding local orbital speed measured by a static observer satisfies $v_{\mathrm{loc}}^2/c^2=r_g/(r-2r_g)$ and tends to one. This suggests a null circular orbit, but one must not simply insert $m=0$ into a massive-particle formula. For a photon, the radial effective potential is proportional to $(1-2r_g/r)/r^2$; its derivative vanishes at $r=3r_g$ and its second derivative there is negative. Thus the photon sphere is at $3GM/c^2$, and its circular null orbits are unstable.

At $r\gg r_g$, expansion givesThe second term is the total Newtonian circular-orbit energy, kinetic plus potential; the relativistic energy also includes rest energy. Here $E$ is the conserved energy referred to a stationary observer at infinity, not the local kinetic energy.

The stable inner edge. At fixed angular momentum, write $\ell =L/(mc)$. The timelike radial geodesic equation can be writtenThe circular condition $V_{\mathrm{eff}}^{\prime }=0$ gives ${\ell }^2=r_g{r}^2/(r-3r_g)$. Substituting into the second derivative taken at fixed $\ell$ givesTimelike circular orbits exist for $r>3r_g$; they are stable against radial perturbations for $r>6r_g$ and unstable for $3r_g<r<6r_g$. ThereforeThe boundary orbit is marginally stable. Differentiating $e(r)$ also gives a minimum there, consistently with this fixed-angular-momentum stability calculation.

accretion efficiency and black-hole spin. Matter supplied slowly from far away initially brings rest energy $d{m}_0{c}^2$. A radiatively efficient, zero-inner-torque thin disk releases the binding energy down to the innermost stable circular orbit; the remaining energy $e_{\mathrm{ISCO}}d{m}_0{c}^2$ is carried into the hole during the plunge. ThusHere $\eta$ denotes the efficiency called $ϵ$ in this question. The ideal relation assumes negligible additional dissipation during plunge and neglects photon capture and advection. A hot radiatively inefficient flow need not have this efficiency.

For a Kerr black hole, an equatorial prograde disk has a smaller innermost stable circular orbit and greater binding-energy efficiency than the Schwarzschild value; a retrograde disk has the opposite trend. The formal extremal prograde limits are $r_{\mathrm{ISCO}}\to r_g$ and $\eta \to 1-1/\sqrt{3}\simeq 0.423$. These are ideal limiting values, not a universal realized efficiency. photon capture and the accretion geometry affect the achievable black-hole spin and escaping radiation.

radiation force and the limiting luminosity. For isotropic luminosity $L$, the outward radiation momentum flux at radius $r$ is $L/(4\pi r^{2}c)$. In fully ionized hydrogen, Thomson scattering acts mainly on the electron, while electrostatic coupling makes the electron and proton move together. Neglecting the electron mass, one electron–proton pair experiencesEquating them yieldsThis Eddington luminosity uses electron-scattering opacity, spherical or isotropic force balance and negligible relativistic corrections at the force-comparison radius. Composition, other opacity, anisotropy and photon trapping can change the relevant limit.

Exponential growth and doubling. Distinguish rest-mass supply ${\dot{m}}_0$ from black-hole mass growth $\dot{M}$:At the Eddington limit,Using the printed scattering cross section gives $t_E\simeq 4.50\times 10^8\mathrm{yr}$. The Salpeter time $t_S$ is an e-folding time; the Eddington-limited doubling time isfor the Schwarzschild thin-disk efficiency. If one approximates $\dot{M}\simeq {\dot{m}}_0$, the factor $1/(1-\eta )$ is dropped and the estimate is about $1.79\times 10^7\mathrm{yr}$. The retained-energy version is the mass-conserving expression.

The efficiency bound. A growth factor $10^8$ requires $N=\log (10^8)=18.4207$ e-folds. Eddington-limited accretion is the fastest growth allowed by the adopted luminosity bound at fixed efficiency, so growth within a time $\Delta t$ requiresThis efficiency upper bound from Eddington-limited growth is an upper bound: higher efficiency radiates more of each supplied rest mass, and the luminosity ceiling therefore allows less mass to enter per unit time. A duty cycle below unity or sub-Eddington luminosity tightens it. The approximate convention $\dot{M}={\dot{m}}_0$ would instead give $\eta \lesssim 0.121$; both figures refer to the same effect, with different treatment of the radiated mass deficit. Mergers, super-Eddington growth or a different seed history would change the premises.

Coherent prograde accretion spins up the hole, moves the inner edge inward and increases $\eta$. For fixed $L/L_{\mathrm{Edd}}$ the mass-growth e-folding time then grows, so the logarithmic growth curve becomes less steep; one cannot simply retain the initial Schwarzschild efficiency throughout. Episodic misaligned accretion can have a different black-hole spin history. In general,