For $h\nu \ll k_{B}T$, expand the exponential in the thermal free-free emissivity:The energy emitted per unit frequency is therefore approximately flat, with the weak Gaunt factor dependence superposed. At these low photon energies, most of the thermal electrons can participate in emission; there is no exponentially small population threshold.

For $h\nu >k_{B}T$, the spectrum turns down approximately as $e^{-h\nu /(k_{B}T)}$, apart from the Gaunt factor. Conservation of energy requires an emitting electron to supply at least $h\nu$. The high-energy tail of its Maxwell-Boltzmann distribution is exponentially rare, which produces this cutoff.

“Flat” refers to $j_{\nu }$, an energy spectrum per unit frequency. The emitted photon number spectrum instead contains $1/(h\nu )$, while energy per logarithmic frequency interval is $\nu j_{\nu }$. At sufficiently low frequency the plasma can also become optically thick to free-free absorption; the given optically thin emissivity is then not the escaping spectrum.

Thermal bremsstrahlung is radiation from the acceleration of electrons in the Coulomb force of ions. At fixed temperature, the encounter rate per electron is proportional to ion number density $n_i$, and the number of radiating electrons per unit volume is $n_e$. Therefore $j_{\nu }\propto n_e{n}_i$ in a dilute, optically thin plasma. This is a binary-encounter rate, not a one-particle emission process.

For a fixed encounter geometry, the ionic electric field and hence the electron acceleration are proportional to ionic charge $Ze$. The Larmor formula is quadratic in acceleration, giving $j_{\nu }\propto Z^2$ at fixed species number densities. The dimensionless Gaunt factor $g(\nu ,T)$ corrects the simple classical result for the detailed quantum scattering and radiation calculation, after averaging over the thermal electron distribution. It depends on frequency and temperature and, in a more complete notation, ionic charge.

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.

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.

Put $\Delta =R-r\ll R$. The near-surface mass density and mean free path from a and b areA simple local escape estimate sets the remaining distance equal to the mean free path, $\ell \simeq \Delta$. It givesFor the Sun's parameters this is $\boxed{\Delta R_{\mathrm{esc},\mathrm{local}}\simeq 1.94\mathrm{km},\Delta R_{\mathrm{esc},\mathrm{local}}/R\simeq 2.79\times 10^{-6}}$.

A more precise operational definition uses the outward radial optical depth. In this varying mass density,The chance of an outward radial photon escaping without another scattering is $e^{-\tau }$. Thus the conventional $\tau =1$ last-scattering depth givesThere is no sharp escape surface: $\tau \ll 1$ gives nearly free propagation, and oblique rays have different optical depths. The local $\ell =\Delta$ estimate actually has $\tau \simeq 1/2$ in this model, not $1$; the factor $\sqrt{2}$ is a convention-dependent order-one correction, not a discrepancy in the Thomson cross-section. Neither scattering criterion establishes a thermalization surface.

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.

Write $x=ar=\pi r/R$. The mass density has a finite central value because $\sin x/x\to 1$ as $x\to 0$, is monotonically decreasing on $0<x<\pi$, and first vanishes at $r=R$. Near the centre, $\rho /{\rho }_c=1-{\pi }^2{r}^2/(6R^2)+O(r^4/R^4)$; near the surface, $\rho /{\rho }_c\simeq (R-r)/R$. It is zero outside. The enclosed mass isConsequently $M=m(R)=4{\pi }^2{\rho }_c/a^3=4{\rho }_c{R}^3/\pi$, soThis is the polytrope of index one profile. Its derivative is negative because $\sin x-x\cos x>0$ for $0<x<\pi$: the latter function vanishes at zero and has derivative $x\sin x>0$.

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.

For the stipulated Strong-shock Rankine-Hugoniot conditions, ${\rho }_d=4{\rho }_u$ and $v_d=v_u/4$. The ideal gas law is $P_d={\rho }_d{k}_B{T}_d/\mu$, with $\mu$ a dimensional mean particle mass. The specific enthalpy of a nonrelativistic monatomic gas is $5k_{B}T/(2\mu )$. Dividing the energy flux equation in d by the common mass flux givesThereforeAs an independent check, conservation of momentum gives $P_d={\rho }_u{v}_u^2-4{\rho }_u(v_u/4{)}^2=3{\rho }_u{v}_u^2/4$, and division by ${\rho }_d=4{\rho }_u$ gives exactly the same temperature. The numerical coefficient relies on $\gamma =5/3$ and the absence of dynamically important radiative or magnetic stresses.

Work in the locally planar shock frame, taking $v_u,v_d>0$ towards the surface. For a stationary normal shock wave with negligible radiative loss inside the discontinuity, the Rankine-Hugoniot conditions for a perfect gas express conservation of mass, conservation of momentum, and conservation of energy per unit area asHere $h$ is specific enthalpy, including the pressure work; using only internal energy would omit that work. The gravitational potential is continuous across an infinitesimally thin front. For cold upstream gas, $P_u$ and $h_u$ are negligible, while a monatomic gas has $\gamma =5/3$ and $h=5P/(2\rho )$.

The hydrodynamic equations apply to flow parallel to the guiding magnetic field when its stress contributions cancel across the front. An oblique magnetohydrodynamic shock requires magnetic stresses and Poynting flux; a radiation-mediated front or appreciable escaping radiation requires the corresponding radiative momentum and energy flux terms. These are additional physics, not terms hidden in $h$ here.

In a sufficiently tight binary star, the main sequence companion can fill its Roche lobe. Roche-lobe overflow then transfers gas through the inner Lagrange point towards the neutron star. The gas carries angular momentum, so it generally forms an accretion disk rather than falling radially onto the surface. Internal stresses carry angular momentum outwards and permit accretion inwards. Capture of a stellar wind is another possible supply mechanism if the companion does not fill its Roche lobe.

Near the neutron star, the rapidly increasing magnetic dipole field, with $B\propto r^{-3}$, becomes dynamically important when its magnetic pressure $B^2/(2{\mu }_0)$ is comparable to the inflow stresses. The accretion disk is truncated around a magnetospheric radius. Ionized material attaches to magnetic field lines and streams along them in magnetically channelled accretion, forming an accretion column above one or both magnetic poles. The Lorentz force impedes motion across the magnetic field, but does not forbid motion parallel to it. Gas can therefore reach the polar surface, with its gravitational potential energy released in a hot shock wave and radiation.

This description presumes the neutron star rotates slowly enough for material to reach it: sufficiently rapid rotation can centrifugally expel the inflow instead. A strong magnetic field by itself does not guarantee surface accretion.

In the ultrarelativistic regime of a, retaining a constant mean fractional energy gain and ignoring losses gives, after $N$ individual crossings,After $n$ complete return cycles, $N=2n$ and $E_n\simeq E_0(1+u/c{)}^n$, to the same first-order accuracy. These products approximate the mean amplification; individual electron crossing angles give a distribution around it.

The downstream electron distribution is nearly isotropic. Its crossing particle flux per unit area is $n_{\mathrm{a}}c/4$, where $n_{\mathrm{a}}$ is the accelerated-particle number density. The downstream flow permanently advects particles away at particle flux $n_{\mathrm{a}}{u}_2$. Their ratio estimates the escape probability after a downstream visit:Upstream particles are swept back towards the shock wave in this idealized infinite-upstream model. Consequently this return probability applies once per complete cycle, not once per individual crossing or local scattering event.

Let $A\simeq 1+u/c$ be the amplification per complete cycle and $P=P_{\mathrm{ret}}$. A steady injection of ${\mathcal{N}}_0$ particles leaves ${\mathcal{N}}_0{P}^n$ surviving to reach $E=E_0{A}^n$. Eliminating $n$ gives the cumulative power law and differential energy distributionThe escape distribution has the same differential exponent when $P$ is independent of energy. Thus the standard Strong-shock Rankine-Hugoniot conditions giveEquivalently the phase-space one-particle distribution function is $f(p)\propto p^{-4}$ for ultrarelativistic particles. Equal logarithmic intervals then contain equal particle energy. A finite source age, upstream escape, and radiative losses eventually introduce cutoffs or breaks.

A particle exponent is not automatically the radiation spectral index. For optically thin uncooled synchrotron radiation from $d\mathcal{N}/dE\propto E^{-s}$, the characteristic frequency scales as $\nu \propto E^2$ and the emitted power as $E^2$. Hence $L_{\nu }\propto E^{-s}{E}^2/(d\nu /dE)\propto {\nu }^{-(s-1)/2}$. With the convention $L_{\nu }\propto {\nu }^{-\alpha }$, this gives $\boxed{\alpha \simeq 1/2}$ and $\nu L_{\nu }\propto {\nu }^{1/2}$. Strong synchrotron cooling steepens a continuously injected particle distribution by one exponent, giving $\alpha \simeq 1$ above its cooling break. A photon number distribution $d{\dot{N}}_{\gamma }/d\epsilon \propto {\epsilon }^{-\Gamma }$ has $\Gamma =\alpha +1$, or $3/2$ in the uncooled case.

Convention qualification. If the printed $P=1-u/c$ is instead imposed independently at every individual crossing, its combination with $A_{\mathrm{crossing}}=1+u/(2c)$ gives $s\simeq 3$ and uncooled $\alpha \simeq 1$. For example $u/c=0.01$ gives $s=1-\ln (0.99)/\ln (1.005)\simeq 3.015$, whereas the consistent cycle model gives $s\simeq 2.010$. The standard $s=2$ estimate requires matching the probability and amplification to the same complete cycle.

Use the shock frame, with upstream and downstream fluid velocities $u_1=u$ and $u_2=u/4$ pointing in the same direction. These are the Strong-shock Rankine-Hugoniot conditions for a nonrelativistic monatomic gas with adiabatic index $5/3$. The two scattering media approach one another at relative speed $U=u_1-u_2=3u/4$, to first order in $u/c$. An energetic electron repeatedly crosses the shock wave, scatters from magnetic turbulence carried by the local plasma, and returns. Each scattering approximately conserves its energy in the local fluid rest frame while changing its direction; comparing successive fluid frames produces a systematic energy and momentum gain. The energy comes from the converging bulk flow, through the moving magnetic field irregularities and their associated electric fields. This is diffusive shock acceleration, a first-order Fermi acceleration process.

Let $\eta$ be the positive cosine of the crossing angle with the shock wave's normal. The Lorentz transformation between fluid frames gives, for either crossing direction with the corresponding sign of the normal reversed,For an almost isotropic population, the number crossing a unit area is weighted by its normal velocity. Thus the normalized crossing distribution is $2\eta d\eta$, rather than a uniform distribution in $\eta$, andFor an ultrarelativistic electron, $E\simeq pc$ and its speed $v\simeq c$, soMore generally $dE/dp=v$ gives $⟨\Delta p⟩/p=2U/(3v)$, while a nonrelativistic kinetic energy has $⟨\Delta K⟩/K=4U/(3v)$. A complete upstream–downstream–upstream cycle contains two crossings and has ultrarelativistic fractional energy gain $u/c$.

The initially cold electrons need an injection stage: shock wave heating and small-scale plasma processes must first produce particles with speed much larger than $u$ that can return across the front. The $v\simeq c$ result is a description of that accelerated population; a cold electron cannot simply be assigned $E\simeq pc$ at its first encounter.