 Solution

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.