The equation of state of an ideal Fermi gas is expressed in terms of Fermi-Dirac integrals. We give formulae for evaluation the Fermi-Dirac integrals of orders 1/2, 3/2, and 5/2 and their derivatives in various limits of non- and extreme degeneracy and relativity. We provide tables and a Fortran subroutine for numerical evaluation of the integrals and derivatives when a limit does not apply. The functions can be evaluated to better than 1% accuracy for any temperature and density using these methods.

We have performed Monte Carlo calculations of gamma-ray transport in models of the clumpy debris cloud of the LMC supernova, SN1987A, to study the influence of composition mixing and heterogeneity on its emergent gamma-ray and X-ray fluxes. In particular, we have focused on the problematic Ginga band (16 -- 28 keV) flux at day 600, whose measured value was an order of magnitude higher than predicted by previous theory. We find that the hydrogen of the envelope could not have been intimately mixed with the heavy elements of the core and that the hydrogen/helium volume filling factor interior to 4000 km/s must have been large ($\ge 40\%$). A physical separation of the scattering region and the regions occupied by the high-Z elements is required. The 600-day models that best fit both the line data at 847 keV and 1238 keV and the measured Ginga band fluxes suggest that as much as 50\% of the explosively produced $^{56}$Ni stayed interior to 1000 km s$^{-1}$ and 2 M$_\odot$. The $^{56}$Ni may have been more centrally-concentrated than in the benchmark models. $^{56}$Ni filling factors greater than 60\% are not preferred, since such models are too good at absorbing photons below 100 keV. Furthermore, a total envelope mass between 10 $\mdot$ and 15 $\mdot$ is favored.

Differential rotation between the neutron star crust and a more rapidly rotating interior superfluid leads to frictional heating that affects the star's long-term thermal evolution and resulting surface emission. Here we present the results of thermal evolution simulations based on two models of the vortex pinning forces that bracket a range of plausible pinning strengths. We include the effects of superfluidity, magnetic fields, and temperature gradients. As representative standard and accelerated neutrino emission processes taking place in the core, we consider the modified Urca process in normal baryonic matter, and the much faster quark Urca process. Comparison of our results with neutron star surface temperature data, including the recent temperature measurement of the Geminga pulsar, shows that stars with soft equations of state and modest frictional heating are in closest agreement with the data; stars with stronger frictional heating have temperatures inconsistent with the upper limit of PSR 1929+10. Stiffer stars undergoing standard cooling generally have temperatures lying above the Vela detection, a situation worsened by the inclusion of frictional heating. Stars undergoing accelerated cooling without frictional heating have temperatures that fall far below most temperature measurements; the Vela and Geminga detections being the most compelling examples. Only in stiff stars, which have thick crusts, can the inclusion of strong frictional heating raise the temperature at late stages in the evolution to a level consistent with the data. However, such a large amount of heating leads to a temperature at $\sim 1000$ yr in excess of the Crab upper limit. Suppression of accelerated neutrino emission processes, perhaps by superfluid pairing in the core, may yield acceptable cooling models.