Transient LMXBs that host neutron stars (NSs) provide excellent laboratories for probing the dense matter physics present in NS crusts. During accretion outbursts in LMXBs, exothermic reactions may heat the NS crust, disrupting the crust-core equilibrium. When the outburst ceases, the crust cools to restore thermal equilibrium with the core. Monitoring this evolution allows us to probe the dense matter physics in the crust. Properties of the deeper crustal layers can be probed at later times after the end of the outburst. We report on the unexpected late-time temperature evolution (>2000 days after the end of their outbursts) of two NSs in LMXBs, XTE J1701-462 and EXO 0748-676. Although both these sources exhibited very different outbursts (in terms of duration and the average accretion rate), they exhibit an unusually steep decay of ~7 eV in the observed effective temperature (occurring in a time span of ~700 days) around ~2000 days after the end of their outbursts. Furthermore, they both showed an even more unexpected rise of ~3 eV in temperature (over a time period of ~500-2000 days) after this steep decay. This rise was significant at the 2.4{\sigma} and 8.5{\sigma} level for XTE J1701-462 and EXO 0748-676, respectively. The physical explanation for such behaviour is unknown and cannot be straightforwardly be explained within the cooling hypothesis. In addition, this observed evolution cannot be well explained by low-level accretion either without invoking many assumptions. We investigate the potential pathways in the theoretical heating and cooling models that could reproduce this unusual behaviour, which so far has been observed in two crust-cooling sources. Such a temperature increase has not been observed in the other NS crust-cooling sources at similarly late times, although it cannot be excluded that this might be a result of the inadequate sampling obtained at such late times.

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object's predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spindown of a pulsar, or thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 year old neutron star. Furthermore, models of cooling neutron stars within the Minimal Cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n-1S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations don't overlap, then pulsar spindown or accretion might be more likely, but the pulsar's period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.

A neo-neutron star is a hot neutron star that has just become transparent to neutrinos. In a core collapse supernova or accretion induced collapse of a white dwarf the neo-neutron star phase directly follows the proto-neutron star phase, about 30 to 60 seconds after the initial collapse. It will also be present in a binary neutron star merger in the case the "born-again" hot massive compact star does not immediately collapse into a black hole. Eddington or even super-Eddington luminosities are present for some time. A neo-neutron star produced in a core collapse supernova is not directly observable but the one produced by a binary merger, likely associated with an off-axis short gamma-ray burst, may be observable for some time as well as when produced in the accretion induced collapse of a white dwarf. We present a first step in the study of this neo-neutron star phase in a spherically symmetric configuration, thus neglecting fast rotation, and also neglecting the effect of strong magnetic fields. We put particular emphasis on determining how long the star can sustain a near-Eddington luminosity and also show the importance of positrons and contraction energy during neo-neutron star phase. We finally discuss the observational prospects for neutron star mergers triggered by LIGO and for accretion-induced collapse transients.

Transient neutron star (NS) LMXBs undergo episodes of accretion, alternated with quiescent periods. During an accretion outburst, the NS heats up due to exothermic accretion-induced processes taking place in the crust. Besides the long-known deep crustal heating of nuclear origin, a likely non-nuclear source of heat, dubbed 'shallow heating', is present at lower densities. Most of the accretion-induced heat slowly diffuses into the core on a timescale of years. Over many outburst cycles, a state of equilibrium is reached when the core temperature is high enough that the heating and cooling (photon and neutrino emission) processes are in balance. We investigate how stellar characteristics and outburst properties affect the long-term temperature evolution of a transiently accreting NS. For the first time the effects of crustal properties are considered, particularly that of shallow heating. Using our code NSCool, we tracked the thermal evolution of a NS undergoing outbursts over a period of $10^5$ yr. The outburst sequence is based on the regular outbursts observed from Aql X-1. For each model, we calculated the timescale over which equilibrium was reached and we present these timescales along with the temperature and luminosity parameters of the equilibrium state. We find that shallow heating significantly contributes to the equilibrium state. Increasing its strength raises the equilibrium core temperature. We find that if deep crustal heating is replaced by shallow heating alone, the core would still heat up, reaching only a 2% lower equilibrium core temperature. Deep crustal heating may therefore not be vital to the heating of the core. Additionally, shallow heating can increase the quiescent luminosity to values higher than previously expected.

The structure and composition of the crust of neutron stars plays an important role in their thermal and magnetic evolution, hence in setting their observational properties. One way to study the crust properties is to measure how it cools after it has been heated during an accretion outburst in a low-mass X-ray binary (LMXB). Such studies have shown that there is a tantalizing source of heat, of currently unknown origin, that is located in the outer layers of the crust and has a strength that varies between different sources and different outbursts. With the aim of understanding the mechanism behind this "shallow heating", we present Chandra and Swift observations of the neutron star LMXB Aql X-1, obtained after its bright 2016 outburst. We find that the neutron star temperature was initially much lower, and started to decrease at much later time, than observed after the 2013 outburst of the source, despite the fact that the properties of the two outbursts were very similar. Comparing our data to thermal evolution simulations, we infer that the depth and magnitude of shallow heating must have been much larger during the 2016 outburst than during the 2013 one. This implies that basic neutron star parameters that do not change between outbursts, do not play a strong role in shallow heating. Furthermore, it suggests that outbursts with a similar accretion morphology can give rise to very different shallow heating. We also discuss alternative explanations for the difference in quiescent evolution after the 2016 outburst.

We present a new Chandra observation (performed in July 2016) of the neutron star X-ray transient IGR J17480-2446, located in the globular cluster Terzan 5. We study the continued cooling of the neutron star crust in this system that was heated during the 2010 outburst of the source. This new observation was performed two years after the last observation of IGR J17480-2446, hence, significantly extending the cooling baseline. We reanalysed all available Chandra observations of the source (but excluding observations during which one of the known transients in Terzan 5 was in outburst) and fitted the obtained cooling curve with our cooling code NSCool, which allows for much improved modelling than what was previously performed for the source. The data and our fit models indicate that the crust was still cooling ~5.5 years after the outburst ended. The neutron star crust has likely not reached crust-core thermal equilibrium yet, and further cooling is predicted (which can be confirmed with additional Chandra observations in >5 years). Intriguingly, we find indications that the thermal conductivity might be relatively low in part of the crust compared to what has been inferred for other crust-cooling sources and tentatively suggest that this layer might be located around the neutron drip. The reason for this difference is unclear, but might be related to the fact that IGR J17480-2446 harbours a relatively slowly rotating neutron star (with a spin of 11 Hz) that has a relatively strong inferred surface magnetic field ($10^{9-10}$ Gauss) compared to what is known or typically assumed for other cooling sources.

The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star.

MASTER OT 075353.88+174907.6 was a blue optical transient reported by the MASTER-Net project on 2017 Oct 31. This source was previously detected by {\it GALEX} in its NUV band but not by the Sloan Digital Sky Survey (in the optical). We carried out multiwavelength follow-up observations of this source during its 2017 outburst using {\it Swift} and RATIR. The source was found to be $\gtrsim$4.4 mag above its quiescent level during the peak of the outburst and the outburst lasted $\gtrsim$19 days. Our observations suggest that it was a superoutburst of a long orbital period U Geminorum type dwarf nova system. The spectral energy distribution during the initial slow decay phase of the outburst was consistent with a disk-dominated spectra (having spectral indices $\Gamma \! \sim$1.5--2.3). After this phase, the UV flux decreased slower than the optical and the spectral energy distribution was very steep with indices $\Gamma \! \sim$3.7$\pm$0.7. This slow decay in the UV may be the emission from a cooling white dwarf heated during the outburst. The spectral shape determined from the assumed pre-outburst quiescent level was also steep ($\Gamma \! \gtrsim$2.5) indicating that the white dwarf is still hot in quiescence (even after the cooling due to the potential accretion-induced heating has halted). No X-ray emission was detected from the source since it is likely located at a large distance $>$2.3 kpc.

A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns y$_H$ > 10$^{7}$ g cm$^{-2}$, while DNB can make helium envelopes more transparent to heat flux for surface temperatures $T_s$ > 2 $\times 10^6$ K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the long-term cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of $\approx$ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.

Monitoring the cooling of neutron-star crusts heated during accretion outbursts allows us to infer the physics of the dense matter present in the crust. We examine the crust cooling evolution of the low-mass X-ray binary MXB 1659-29 up to ~505 days after the end of its 2015 outburst (hereafter outburst II) and compare it with what we observed after its previous 1999 outburst (hereafter outburst I) using data obtained from the Swift, XMM-Newton, and Chandra observatories. The observed effective surface temperature of the neutron star in MXB 1659-29 dropped from ~92 eV to ~56 eV from ~12 days to ~505 days after the end of outburst II. The most recently performed observation after outburst II suggests that the crust is close to returning to thermal equilibrium with the core. We model the crust heating and cooling for both its outbursts collectively to understand the effect of parameters that may change for every outburst (e.g., the average accretion rate, the length of outburst, the envelope composition of the neutron star at the end of the outburst) and those which can be assumed to remain the same during these two outbursts (e.g., the neutron star mass, its radius). Our modelling indicates that all parameters were consistent between the two outbursts with no need for any significant changes. In particular, the strength and the depth of the shallow heating mechanism at work (in the crust) were inferred to be the same during both outbursts, contrary to what has been found when modelling the cooling curves after multiple outburst of another source, MAXI J0556-332. This difference in source behaviour is not understood. We discuss our results in the context of our current understanding of cooling of accretion-heated neutron-star crusts, and in particular with respect to the unexplained shallow heating mechanism.

If the thermal evolution of the hot young neutron star in the supernova remnant HESS J1731-347 is driven by neutrino emission, it provides a stringent constraint on the coupling of light (mass $\ll 10$ keV) axion-like particles to neutrons. Using Markov-Chain Monte Carlo we find that for the values of axion-neutron coupling $g_{ann}^2 > 7.7 \times 10^{-20}$ (90% c.l.) the axion cooling from the bremsstrahlung reaction $n+n\rightarrow n+n +a$ is too rapid to account for the high observed surface temperature. This implies that the Pecci-Quinn scale or axion decay constant $f_a > 6.7 \times 10^7$ GeV for KSVZ axions and $f_a > 1.7 \times 10^9$ GeV for DFSZ axions. The high temperature of this neutron star also allows us to tighten constraints on the size of the nucleon pairing gaps.

A unique signature for the presence of massive black holes in very dense stellar regions is occasional giant-amplitude outbursts of multiwavelength radiation from tidal disruption and subsequent accretion of stars that make a close approach to the black holes. Previous strong tidal disruption event (TDE) candidates were all associated with the centers of largely isolated galaxies. Here we report the discovery of a luminous X-ray outburst from a massive star cluster at a projected distance of 12.5 kpc from the center of a large lenticular galaxy. The luminosity peaked at ~10^{43} erg/s and decayed systematically over 10 years, approximately following a trend that supports the identification of the event as a TDE. The X-ray spectra were all very soft, with emission confined to be <3.0 keV, and could be described with a standard thermal disk. The disk cooled significantly as the luminosity decreased, a key thermal-state signature often observed in accreting stellar-mass black holes. This thermal-state signature, coupled with very high luminosities, ultrasoft X-ray spectra and the characteristic power-law evolution of the light curve, provides strong evidence that the source contains an intermediate-mass black hole (IMBH) with a mass of a few ten thousand solar mass. This event demonstrates that one of the most effective means to detect IMBHs is through X-ray flares from TDEs in star clusters.

With our neutron star crust cooling code {\tt NSCool} we track the thermal evolution of the neutron star in Aql X-1 over the full accretion outburst history from 1996 until 2015. For the first time, we model many outbursts (23 outbursts were detected) collectively and in great detail. This allows us to investigate the influence of previous outbursts on the internal temperature evolution and to test different neutron star crust cooling scenarios. Aql X-1 is an ideal test source for this purpose, because it shows frequent, short outbursts and thermally dominated quiescence spectra. The source goes into outburst roughly once a year for a few months. Assuming that the quiescent {\it Swift}/XRT observations of Aql X-1 can be explained within the crust cooling scenario (Waterhouse et al. 2016), we find three main conclusions. Firstly, the data are well reproduced by our model if the envelope composition and shallow heating parameters are allowed to change between outbursts. This is not the case if both shallow heating parameters (strength and depth) are tied throughout all accretion episodes, supporting earlier results that the properties of the shallow heating mechanism are not constant between outbursts. Second, from our models shallow heating could not be connected to one specific spectral state during outburst. Third, and most importantly, we find that the neutron star in Aql X-1 does not have enough time between outbursts to cool down to crust-core equilibrium and that heating during one outburst influences the cooling curves of the next.

We studied the transient neutron-star low-mass X-ray binary GRS 1747-312, located in the globular cluster Terzan 6, in its quiescent state after its outburst in August 2004, using an archival XMM-Newton observation. A source was detected in this cluster and its X-ray spectrum can be fitted with the combination of a soft, neutron-star atmosphere model and a hard, power-law model. Both contributed roughly equally to the observed 0.5-10 keV luminosity (~4.8X10^33 erg/s). This type of X-ray spectrum is typically observed for quiescent neutron-star X-ray transients that are perhaps accreting in quiescence at very low rates. Therefore, if this X-ray source is the quiescent counterpart of GRS 1747-312, then this source is also accreting at low levels in-between outbursts. Since source confusion a likely problem in globular clusters, it is quite possible that part, if not all, of the emission we observed is not related to GRS 1747-312, and is instead associated with another source or conglomeration of sources in the cluster. Currently, it is not possible to determine exactly which part of the emission truly originates from GRS 1747-312, and a Chandra observation (when no source is in outburst in Terzan 6) is needed to be conclusive. Assuming that the detected emission is due to GRS 1747-312, we discuss the observed results in the context of what is known about other quiescent systems. We also investigated the thermal evolution of the neutron star in GRS 1747-312, and inferred that GRS 1747-312 can be considered a typical quiescent system under our assumptions.

We report on two new quiescent {\it XMM-Newton} observations (in addition to the earlier {\it Swift}/XRT and {\it XMM-Newton} coverage) of the cooling neutron star crust in the low-mass X-ray binary 1RXS J180408.9$-$342058. Its crust was heated during the $\sim$4.5 month accretion outburst of the source. From our quiescent observations, fitting the spectra with a neutron star atmosphere model, we found that the crust had cooled from $\sim$ 100 eV to $\sim$73 eV from $\sim$8 days to $\sim$479 days after the end of its outburst. However, during the most recent observation, taken $\sim$860 days after the end of the outburst, we found that the crust appeared not to have cooled further. This suggested that the crust had returned to thermal equilibrium with the neutron star core. We model the quiescent thermal evolution with the theoretical crustal cooling code NSCool and find that the source requires a shallow heat source, in addition to the standard deep crustal heating processes, contributing $\sim$0.9 MeV per accreted nucleon during outburst to explain its observed temperature decay. Our high quality {\it XMM-Newton} data required an additional hard component to adequately fit the spectra. This slightly complicates our interpretation of the quiescent data of 1RXS J180408.9$-$342058. The origin of this component is not fully understood.

We present a quantitative analysis of superfluidity and superconductivity in dense matter from observations of isolated neutron stars in the context of the minimal cooling model. Our new approach produces the best fit neutron triplet superfluid critical temperature, the best fit proton singlet superconducting critical temperature, and their associated statistical uncertainties. We find that the neutron triplet critical temperature is likely $2.09^{+4.37}_{-1.41} \times 10^{8}$ K and that the proton singlet critical temperature is $7.59^{+2.48}_{-5.81} \times 10^{9}$ K. However, we also show that this result only holds if the Vela neutron star is not included in the data set. If Vela is included, the gaps increase significantly to attempt to reproduce Vela's lower temperature given its young age. Further including neutron stars believed to have carbon atmospheres increases the neutron critical temperature and decreases the proton critical temperature. Our method demonstrates that continued observations of isolated neutron stars can quantitatively constrain the nature of superfluidity in dense matter.

We show that the neutron star in the transient system MXB~1659-29 has a core neutrino luminosity that substantially exceeds that of the modified Urca reactions (i.e., $n+n\to n+p+e^{-}+\bar{\nu}_{e}$ and inverse) and is consistent with the direct Urca reactions ($n\to p+e^{-}+\bar{\nu}_{e}$ and inverse) occurring in a small fraction of the core. Observations of the thermal relaxation of the neutron star crust following 2.5 years of accretion allow us to measure the energy deposited into the core during accretion, which is then reradiated as neutrinos, and infer the core temperature. For a nucleonic core, this requires that the nucleons are unpaired and that the proton fraction exceed a critical value to allow the direct Urca reaction to proceed. The neutron star in MXB~1659-29 is the first with a firmly detected thermal component in its X-ray spectrum that needs a fast neutrino cooling process. Measurements of the temperature variation of the neutron star core during quiescence would place an upper limit on the core specific heat and serve as a check on the fraction of the neutron star core in which nucleons are unpaired.

In accreting neutron star X-ray transients, the neutron star crust can be substantially heated out of thermal equilibrium with the core during an accretion outburst. The observed subsequent cooling in quiescence (when accretion has halted) offers a unique opportunity to study the structure and thermal properties of the crust. Initially crust cooling modelling studies focussed on transient X-ray binaries with prolonged accretion outbursts (> 1 year) such that the crust would be significantly heated for the cooling to be detectable. Here we present the results of applying a theoretical model to the observed cooling curve after a short accretion outburst of only ~10 weeks. In our study we use the 2010 outburst of the transiently accreting 11 Hz X-ray pulsar in the globular cluster Terzan 5. Observationally it was found that the crust in this source was still hot more than 4 years after the end of its short accretion outburst. From our modelling we found that such a long-lived hot crust implies some unusual crustal properties such as a very low thermal conductivity (> 10 times lower than determined for the other crust cooling sources). In addition, we present our preliminary results of the modelling of the ongoing cooling of the neutron star in MXB 1659-298. This transient X-ray source went back into quiescence in March 2017 after an accretion phase of ~1.8 years. We compare our predictions for the cooling curve after this outburst with the cooling curve of the same source obtained after its previous outburst which ended in 2001.

The transient neutron star (NS) low-mass X-ray binary MAXI J0556$-$332 provides a rare opportunity to study NS crust heating and subsequent cooling for multiple outbursts of the same source. We examine {\it MAXI}, {\it Swift}, {\it Chandra}, and {\it XMM-Newton} data of MAXI J0556$-$332 obtained during and after three accretion outbursts of different durations and brightness. We report on new data obtained after outburst III. The source has been tracked up to $\sim$1800 d after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from $\sim$333 eV to $\sim$146 eV was observed during the first $\sim$1200 d. Outburst III reheated the crust up to $\sim$167 eV, after which the crust cooled again to $\sim$131 eV in $\sim$350 d. We model the thermal evolution of the crust and find that this source required a different strength and depth of shallow heating during each of the three outbursts. The shallow heating released during outburst I was $\sim$17 MeV nucleon$^{-1}$ and outburst III required $\sim$0.3 MeV nucleon$^{-1}$. These cooling observations could not be explained without shallow heating. The shallow heating for outburst II was not well constrained and could vary from $\sim$0--2.2 MeV nucleon$^{-1}$, i.e., this outburst could in principle be explained without invoking shallow heating. We discuss the nature of the shallow heating and why it may occur at different strengths and depths during different outbursts.

Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.

We present a brief, observational review about the study of the cooling behaviour of accretion-heated neutron stars and the inferences about the neutron-star crust and core that have been obtained from these studies. Accretion of matter during outbursts can heat the crust out of thermal equilibrium with the core and after the accretion episodes are over, the crust will cool down until crust-core equilibrium is restored. We discuss the observed properties of the crust cooling sources and what has been learned about the physics of neutron-star crusts. We also briefly discuss those systems that have been observed long after their outbursts were over, i.e, during times when the crust and core are expected to be in thermal equilibrium. The surface temperature is then a direct probe for the core temperature. By comparing the expected temperatures based on estimates of the accretion history of the targets with the observed ones, the physics of neutron-star cores can be investigated. Finally, we discuss similar studies performed for strongly magnetized neutron stars in which the magnetic field might play an important role in the heating and cooling of the neutron stars.

We have monitored the transient neutron star low-mass X-ray binary 1RXS J180408.9-342058 in quiescence after its ~4.5 month outburst in 2015. The source has been observed using Swift and XMM-Newton. Its X-ray spectra were dominated by a thermal component. The thermal evolution showed a gradual X-ray luminosity decay from ~18 x 10^32 to ~4 x 10^32 (D/5.8 kpc)^2 erg s^{-1} between ~8 and ~379 d in quiescence, and the inferred neutron star surface temperature (for an observer at infinity; using a neutron star atmosphere model) decreased from ~100 to ~71 eV. This can be interpreted as cooling of an accretion-heated neutron star crust. Modelling the observed temperature curve (using NSCOOL) indicated that the source required ~1.9 MeV per accreted nucleon of shallow heating in addition to the standard deep crustal heating to explain its thermal evolution. Alternatively, the decay could also be modelled without the presence of deep crustal heating, only having a shallow heat source (again ~1.9 MeV per accreted nucleon was required). However, the XMM-Newton data statistically required an additional power-law component. This component contributed ~30 per cent of the total unabsorbed flux in 0.5-10 keV energy range. The physical origin of this component is unknown. One possibility is that it arises from low-level accretion. The presence of this component in the spectrum complicates our cooling crust interpretation because it might indicate that the smooth luminosity and temperature decay curves we observed may not be due to crust cooling but due to some other process.

Crustal cooling of accretion-heated neutron stars provides insight into the stellar interior of neutron stars. The neutron star X-ray transient, KS~1731$-$260, was in outburst for 12.5 years before returning to quiescence in 2001. We have monitored the cooling of this source since then through {\it Chandra} and {\it XMM-Newton} observations. Here, we present a 150 ks {\it Chandra} observation of KS~1731$-$260 taken in August 2015, about 14.5 years into quiescence, and 6 years after the previous observation. We find that the neutron star surface temperature is consistent with the previous observation, suggesting that crustal cooling has likely stopped and the crust has reached thermal equilibrium with the core. Using a theoretical crust thermal evolution code, we fit the observed cooling curves and constrain the core temperature (T$_c = 9.35\pm0.25\times10^7$ K), composition (Q$_{imp} = 4.4^{+2.2}_{-0.5}$) and level of extra shallow heating required (Q$_{sh} = 1.36\pm0.18$ MeV/nucleon). We find that the presence of a low thermal conductivity layer, as expected from nuclear pasta, is not required to fit the cooling curve well, but cannot be excluded either.

The neutron star low-mass X-ray binary and intermittent millisecond X-ray pulsar HETE J1900.1-2455 returned to quiescence in late 2015, after a prolonged accretion outburst of ~10 yr. Using a Chandra observation taken ~180 d into quiescence we detect the source at a luminosity of ~4.5E31 (D/4.7 kpc)^2 erg/s (0.5-10 keV). The X-ray spectrum can be described by a neutron star atmosphere model with a temperature of ~54 eV for an observer at infinity. We perform thermal evolution calculations based on the 2016 quiescent data and a <98 eV temperature upper limit inferred from a Swift observation taken during an unusually brief (<2 weeks) quiescent episode in 2007. We find no evidence in the present data that the thermal properties of the crust, such as the heating rate and thermal conductivity, are different than those of non-pulsating neutron stars. Finding this neutron star so cold after its long outburst imposes interesting constraints on the heat capacity of the stellar core; these become even stronger if further cooling were to occur.

We show that observations of the core temperature of transiently-accreting neutron stars combined with observations of an accretion outburst give a lower limit to the neutron star core heat capacity. For the neutron stars in the low mass X-ray binaries KS 1731-260, MXB 1659-29, and XTE J1701-462, we show that the lower limit is a factor of a few below the core heat capacity expected if neutrons and protons in the core are paired, so that electrons provide the dominant contribution to the heat capacity. This limit rules out a core dominated by a quark color-flavor-locked (CFL) phase, which would have a much lower heat capacity. Future observations of or limits on cooling during quiescence will further constrain the core heat capacity.

Using a theoretical model, we track the thermal evolution of a cooling neutron star crust after an accretion induced heating period with the goal of constraining the crustal parameters. We present for the first time a crust cooling model $-\text{ } NSCool\text{ } -$ that takes into account detailed variability during the full outburst based on the observed light curve. We apply our model to KS 1731-260. The source was in outburst for $\sim$12 years during which it was observed to undergo variations on both long (years) and short (days-weeks) timescales. Our results show that KS 1731-260 does not reach a steady state profile during the outburst due to fluctuations in the derived accretion rate. Additionally, long time-scale outburst variability mildly affects the complete crust cooling phase, while variations in the final months of the outburst strongly influence the first $\sim$40 days of the calculated cooling curve. We discuss the consequences for estimates of the neutron star crust parameters, and argue that detailed modelling of the final phase of the outburst is key to constraining the origin of the shallow heat source.

An accretion outburst in an X-ray transient deposits material onto the neutron star primary; this accumulation of matter induces reactions in the neutron star's crust. During the accretion outburst these reactions heat the crust out of thermal equilibrium with the core. When accretion halts, the crust cools to its long-term equilibrium temperature on observable timescales. Here we examine the accreting neutron star transient MAXI J0556-332, which is the hottest transient, at the start of quiescence, observed to date. Models of the quiescent light curve require a large deposition of heat in the shallow outer crust from an unknown source. The additional heat injected is $\approx 4\textrm{-}10\,\mathrm{MeV}$ per accreted nucleon; when the observed decline in accretion rate at the end of the outburst is accounted for, the required heating increases to $\approx 6\textrm{-}16\,\mathrm{MeV}$. This shallow heating is still required to fit the lightcurve even after taking into account a second accretion episode, uncertainties in distance, and different surface gravities. The amount of shallow heating is larger than that inferred for other neutron star transients and is larger than can be supplied by nuclear reactions or compositionally driven convection; but it is consistent with stored mechanical energy in the accretion disk. The high crust temperature ($T_b \gtrsim 10^{9} \, {\rm K}$) makes its cooling behavior in quiescence largely independent of the crust composition and envelope properties, so that future observations will probe the gravity of the source. Fits to the lightcurve disfavor the presence of Urca cooling pairs in the crust.

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.

Neutron stars provide a fertile environment for exploring superfluidity under extreme conditions. It is not surprising that Cooper pairing occurs in dense matter since nucleon pairing is observed in nuclei as energy differences between even-even and odd-even nuclei. Since superfluids and superconductors in neutron stars profoundly affect neutrino emissivities and specific heats, their presence can be observed in the thermal evolution of neutron stars. An ever-growing number of cooling neutron stars, now amounting to 13 thermal sources, and several additional objects from which upper limits to temperatures can be ascertained, can now be used to discriminate among theoretical scenarios and even to dramatically restrict properties of nucleon pairing at high densities. In addition, observations of pulsars, including their spin-downs and glitch histories, additionally support the conjecture that superfluidity and superconductivity are ubiquitous within, and important to our understanding of, neutron stars.

We explore the spectrum of low-energy collective excitations in the crust of a neutron star, especially in the inner region where neutron-proton clusters are immersed in a sea of superfluid neutrons. The speeds of the different modes are calculated systematically from the nuclear energy density functional theory using a Skyrme functional fitted to essentially all experimental atomic mass data.

Our current description of the large-scale Universe is now known with a precision undreamt of a generation ago. Within the simple standard cosmological model only six basic parameters are required. The usual parameter set includes quantities most directly probed by the cosmic microwave background, but the nature of these quantities is somewhat esoteric. However, many more numbers can be derived that quantify various aspects of our Universe. Using constraints from the Planck satellite, in combination with other data sets, we explore several such quantities, highlighting some specific examples.

It is now possible to model thermal relaxation of neutron stars after bouts of accretion during which the star is heated out of equilibrium by nuclear reactions in its crust. Major uncertainties in these models can be encapsulated in modest variations of a handful of fudge parameters that change the crustal thermal conductivity, specific heat, and heating rates. Observations of thermal relaxation constrain these fudge parameters and allow us to predict longer term variability in terms of the neutron star core temperature. We demonstrate this explicitly by modeling ongoing thermal relaxation in the neutron star XTE J1701-462. Its future cooling, over the next 5 to 30 years, is strongly constrained and depends mostly on its core temperature, uncertainties in crust physics having essentially been pinned down by fitting to the first three years of observations.

A plasma made out of strange-quark matter (SQM) and electrons, has a rather high plasma frequency (>20 MeV). Thus, a compact star made of such material all the way up to its surface, i.e., a bare strange star, would be unable to radiate away its thermal emission. We use the MIT-bag model and assume that SQM is the ground state of nuclear matter at high density. We investigate whether the presence of a magnetic field will allow propagation of radiation at frequencies below the SQM plasma frequencies. Hence, we study the presence of gyrofrequencies in a SQM plasma permeated by a strong magnetic field (B > 10^{12} G). We find that small regions in the frequency spectrum allow radiation propagation due to the presence of the magnetic fields. It is likely that narrow bands of radiation would likely be observable from magnetized bare strange stars .

We study the low-energy collective excitations of the neutron star inner crust, where a neutron superfluid coexists with a Coulomb lattice of nuclei. The dispersion relation of the modes is calculated systematically from a microscopic theory including neutron band structure effects. These effects are shown to lead to a strong mixing between the Bogoliubov-Anderson bosons of the neutron superfluid and the longitudinal crystal lattice phonons. In addition, the speed of the transverse shear mode is greatly reduced as a large fraction of superfluid neutrons are entrained by nuclei. Not only does the much smaller velocity of the transverse mode increase the specific heat of the inner crust, but it also decreases its electron thermal conductivity. These results may impact our interpretation of the thermal relaxation in accreting neutron stars. Due to strong mixing, the mean free path of the superfluid mode is found to be greatly reduced. Our results for the collective mode dispersion relations and their damping may also have implications for neutron star seismology.

It is assumed that accreting neutron stars (NSs) in LMXBs are heated due to the compression of the existing crust by the accreted matter which gives rise to nuclear reactions in the crust. It has been shown that most of the energy is released deep in the crust by pycnonuclear reactions involving low-Z elements. We discuss if NSs in very-faint X-ray transients (VFXTs; those which have peak X-ray luminosities < 1E36 erg/s) can be used to test this model. Unfortunately we cannot conclusively answer this because of the large uncertainties in our estimates of the accretion rate history of those VFXTs, both the short-term (less than a few tens of thousands of years) and the one throughout their lifetime. The latter is important because it can be so low that the NSs might not have accreted enough matter to become massive enough that enhanced cooling processes become active. Therefore, they could be relatively warm compared to other systems for which such enhanced cooling processed have been inferred. However, the amount of matter can also not be too low because then the crust might not have been replaced significantly by accreted matter and thus a hybrid crust of partly accreted and partly original, albeit further compressed matter, might be present. This would inhibit the full range of pycnonuclear reactions to occur and thus very likely decreasing the amount of heat deposited in the crust. Furthermore, better understanding is needed how a hybrid crust affects other properties such as the thermal conductivity. We also show that some individual NS LMXBs might have hybrid crusts as well as the NSs in HMXBs. This has to be taken into account when studying the cooling properties of those systems when they are in quiescence. We show that the VFXTs are likely not the dominate transients that are associated with the brightest low-luminosity X-ray sources in globular clusters as was hypothesized.

We present magnetohydrodynamic numerical simulations of the late post-supernova hypercritical accretion to understand its effect on the magnetic field of the new-born neutron star. We consider as an example the case of a magnetic field loop protruding from the star's surface. The accreting matter is assumed to be non magnetized and, due to the high accretion rate, matter pressure dominates over magnetic pressure. We find that an accretion envelope develops very rapidly and once it becomes convectively stable the magnetic field is easily buried and pushed into the newly forming neutron star crust. However, for low enough accretion rates the accretion envelope remains convective for an extended period of time and only partial submergence of the magnetic field occurs due to a residual field that is maintained at the interface between the forming crust and the convective envelope. In this latter case, the outcome should be a weakly magnetized neutron star with a likely complicated field geometry. In our simulations we find the transition from total to partial submergence to occur around dotM ~ 10 M_sun/yr. Back-diffusion of the submerged magnetic field toward the surface, and the resulting growth of the dipolar component, may result in a delayed switch-on of a pulsar on time-scales of centuries to millenia.

In this review, I present a brief summary of the impact of nucleon pairing at supra-nuclear densities on the cooling of neutron stars. I also describe how the recent observation of the cooling of the neutron star in the supernova remnant Cassiopeia A may provide us with the first direct evidence for the occurrence of such pairing. It also implies a size of the neutron 3P-F2 energy gap of the order of 0.1 MeV.

We review the nuclear and condensed matter physics underlying the thermal and transport properties of the neutron star inner crust. These properties play a key role in interpreting transient phenomena such as thermal relaxation in accreting neutron stars, superbursts, and magnetar flares. We emphasize simplifications that occur at low temperature where the inner crust can be described in terms of electrons and collective excitations. The heat conductivity and heat capacity of the solid and superfluid phase of matter is discussed in detail and we emphasize its role in interpreting observations of neutron stars in soft X-ray transients. We highlight recent theoretical and observational results, and identify future work needed to better understand a host of transient phenomena in neutron stars.

The supernova remnant Cassiopeia A contains the youngest known neutron star which is also the first one for which real time cooling has ever been observed. In order to explain the rapid cooling of this neutron star, we first present the fundamental properties of neutron stars that control their thermal evolution with emphasis on the neutrino emission processes and neutron/proton superfluidity/superconductivity. Equipped with these results, we present a scenario in which the observed cooling of the neutron star in Cassiopeia A is interpreted as being due to the recent onset of neutron superfluidity in the core of the star. The manner in which the earlier occurrence of proton superconductivity determines the observed rapidity of this neutron star's cooling is highlighted. This is the first direct evidence that superfluidity and superconductivity occur at supranuclear densities within neutron stars.

We propose that the observed cooling of the neutron star in Cassiopeia A is due to enhanced neutrino emission from the recent onset of the breaking and formation of neutron Cooper pairs in the 3P2 channel. We find that the critical temperature for this superfluid transition is ~0.5x10^9 K. The observed rapidity of the cooling implies that protons were already in a superconducting state with a larger critical temperature. Our prediction that this cooling will continue for several decades at the present rate can be tested by continuous monitoring of this neutron star.

The properties of a new-born neutron star, produced in a core-collapse supernova, can be strongly affected by the possible late fallback which occurs several hours after the explosion. This accretion occurs in the regime dominated by neutrino cooling, explored initially in this context by Chevalier (1989). Here we revisit this approach in a 1D spherically symmetric model and carry out numerical simulations in 2D in an accretion column onto a neutron star considering detailed microphysics, neutrino cooling and the presence of magnetic fields in ideal MHD. We compare our numerical results to the analytic solutions and explore how the purely hydrodynamical as well as the MHD solutions differ from them, and begin to explore how this may affect the appearance of the remnant as a typical radio pulsar.

Several authors have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere or chromosphere around the black hole due to QED or QCD interactions between the emitted particles. Analyzing these models we identify a number of physical and geometrical effects which invalidate them. In all cases, we find that the observational signatures of a cosmic or Galactic background of black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability.

The minimal cooling paradigm for neutron star cooling assumes that enhanced cooling due to neutrino emission from any direct Urca process, due either to nucleons or to exotica such as hyperons, Bose condensates, or deconfined quarks, does not occur. This scenario was developed to replace and extend the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation processes that occur near the critical temperature for superfluid/superconductor pairing. Recently, it has been found that Cooper-pair neutrino emission from the vector channel is suppressed by a large factor compared to the original estimates that violated vector current conservation. We show that Cooper-pair neutrino emission remains, nevertheless, an efficient cooling mechanism through the axial channel. As a result, the elimination of neutrino emission from Cooper-paired nucleons through the vector channel has only minor effects on the long-term cooling of neutron stars within the minimal cooling paradigm. We further quantify precisely the effect of the size of the neutron 3P2 gap and demonstrate that consistency between observations and the minimal cooling paradigm requires that the critical temperature T_c for this gap covers a range of values between T_c^min < 0.2 x 10^9 K up to T_c^max > 0.5 \times 10^9 K in the core of the star. In addition, it is required that young neutron stars have heterogenous envelope compositions: some must have light-element compositions and others must have heavy-element compositions. Unless these two conditions are fulfilled, about half of the observed young cooling neutron stars are inconsistent with the minimal cooling paradigm and provide evidence for the existence of enhanced cooling.

We discuss a variety of bremsstrahlung processes associated with charged particles emitted by evaporating black holes. We show that such particles produce a negligible number of bremsstrahlung photons from their scattering off each other, though at low frequencies inner bremsstrahlung photons dominate over the direct Hawking emission of photons. This analysis and the further analysis of the accompanying paper invalidate Heckler's claim that sufficiently hot evaporating black holes form QED photospheres.

Several authors, most notably Heckler, have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere around the black hole due to QED or QCD interactions between the emitted particles. In this paper we analyze these claims and identify a number of physical and geometrical effects which invalidate these scenarios. We point out two key problems. First, the interacting particles must be causally connected to interact, and this condition is satisfied by only a small fraction of the emitted particles close to the black hole. Second, a scattered particle requires a distance ~ E/m_e^2 for completing each bremsstrahlung interaction, with the consequence that it is improbable for there to be more than one complete bremsstrahlung interaction per particle near the black hole. These two effects have not been included in previous analyses. We conclude that the emitted particles do not interact sufficiently to form a QED photosphere. Similar arguments apply in the QCD case and prevent a QCD photosphere (chromosphere) from developing when the black hole temperature is much greater than Lambda_QCD, the threshold for QCD particle emission. Additional QCD phenomenological arguments rule out the development of a chromosphere around black hole temperatures of order Lambda_QCD. In all cases, the observational signatures of a cosmic or Galactic halo background of primordial black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability. We also consider the possibility, as proposed by Belyanin et al. and D. Cline et al., that plasma interactions between the emitted particles form a photosphere, and we conclude that this scenario too is not supported.

Observations that we are highly unlikely to be vacuum fluctuations suggest that our universe is decaying at a rate faster than the asymptotic volume growth rate, in order that there not be too many observers produced by vacuum fluctuations to make our observations highly atypical. An asymptotic linear e-folding time of roughly 16 Gyr (deduced from current measurements of cosmic acceleration) would then imply that our universe is more likely than not to decay within a time that is less than 19 Gyr in the future.

Ultra-high-energy cosmic rays (UHECRs) have been tried to be related to the most varied and powerful sources known in the universe. Gamma-ray bursts (GRBs) are natural candidates. Here, we argue that cosmic rays can be accelerated by large amplitude electromagnetic waves (LAEMWs) when the MHD approximation of the field in the wind generated by the GRB's magnetized central engine breaks down. The central engine considered here is a strange star born with differential rotation from the accretion induced conversion of a neutron star into a strange star in a low-mass X-ray binary system. The LAEMWs generated this way accelerate light ions to the highest energies $E = q\eta\Delta\Phi_{max}$ with an efficiency $\eta \sim 10^{-1}$ that accounts for all plausible energy losses. Alternatively, we also consider the possibility that, once formed, the LAEMWs are unstable to creation of a relativistically strong electromagnetic turbulence due to an overturn instability. Under this assumption, a lower limit to the efficiency of acceleration is estimated to be about $\eta \sim 10^{-2.5}$. Due to their age, low mass X-ray binary systems can be located in regions of low interstellar medium density as, e.g., globular clusters or even intergalactic medium in case of high proper motion systems, and cosmic ray energy losses due to proton collisions with photons at the decelerating region are avoided, thus opening the possibility for particles to exploit the full voltage available, well beyond that currently observed.

RX J1856.5-3754 has been proposed as a strange star candidate due to its very small apparent radius measured from its X-ray thermal spectrum. However, its optical emission requires a much larger radius and thus most of the stellar surface must be cold and undetectable in X-rays. In the case the star is a neutron star such a surface temperature distribution can be explained by the presence of a strong toroidal field in the crust (Perez-Azorin et al. 2006, Geppert et al. 2006). We consider a similar scenario for a strange star with a thin baryonic crust to determine if such a magnetic field induced effect is still possible.

Reformulating the instantons in a complex plane for tunneling or transmitting states, we calculate the pair-production rate of charged fermions in a spatially localized electric field, illustrated by the Sauter electric field E_0 sech^2 (z/L), and in a temporally localized electric field such as E_0 sech^2 (t/T). The integration of the quadratic part of WKB instanton actions over the frequency and transverse momentum leads to the pair-production rate obtained by the worldline instanton method, including the prefactor, of Phys. Rev. D72, 105004 (2005) and D73, 065028 (2006). It is further shown that the WKB instanton action plus the next-to-leading order contribution in spinor QED equals the WKB instanton action in scalar QED, thus justifying why the WKB instanton in scalar QED can work for the pair production of fermions. Finally we obtain the pair-production rate in a spatially localized electric field together with a constant magnetic field in the same direction.

We present models of temperature distribution in the crust of a neutron star in the presence of a strong toroidal component superposed to the poloidal component of the magnetic field. The presence of such a toroidal field hinders heat flow toward the surface in a large part of the crust. As a result, the neutron star surface presents two warm regions surrounded by extended cold regions and has a thermal luminosity much lower than in the case the magnetic field is purely poloidal. We apply these models to calculate the thermal evolution of such neutron stars and show that the lowered photon luminosity naturally extends their life-time as detectable thermal X-ray sources.

Unless our universe is decaying at an astronomical rate (i.e., on the present cosmological timescale of Gigayears, rather than on the quantum recurrence timescale of googolplexes), it would apparently produce an infinite number of observers per comoving volume by thermal or vacuum fluctuations (Boltzmann brains). If the number of ordinary observers per comoving volume is finite, this scenario seems to imply zero likelihood for us to be ordinary observers and minuscule likelihoods for our actual observations. Hence, our observations suggest that this scenario is incorrect and that perhaps our universe is decaying at an astronomical rate.

It is shown how pair production itself would almost certainly prevent the astrophysical formation of macroscopic dyadospheres, hypothetical regions, extending many electron Compton wavelengths in all directions, where the electric field exceeds the critical value for microscopically rapid Schwinger pair production. Pair production is a self-regulating process that would discharge a growing electric field, in the example of a hypothetical collapsing charged stellar core, before it reached 6% of the minimum dyadosphere value, keeping the pair production rate more than 26 orders of magnitude below the dyadosphere value, and keeping the efficiency below 0.0002 (M/M_sun)^{1/2}.

We review theoretical developments in studies of dense matter and its phase structure of relevance to compact stars. Observational data on compact stars, which can constrain the properties of dense matter, are presented critically and interpreted.

Pair production itself prevents the development of dyadospheres, hypothetical macroscopic regions where the electric field exceeds the critical Schwinger value. Pair production is a self-regulating process that would discharge a growing electric field, in the example of a hypothetical collapsing charged stellar core, before it reached 6% of the minimum dyadosphere value, keeping the pair production rate more than 26 orders of magnitude below the dyadosphere value.

It is shown how pair production itself would prevent the astrophysical formation of dyadospheres, hypothetical regions where the electric field exceeds the critical value for rapid Schwinger pair production. Pair production is a self-regulating process that would discharge a growing electric field, in the example of a hypothetical collapsing charged stellar core, before it reached 6% of the minimum dyadosphere value, keeping the pair production rate more than 26 orders of magnitude below the dyadosphere value, and keeping the efficiency below 0.0002(M/M_sun)^{1/2}.

We continue the study of the effects of a strong magnetic field on the temperature distribution in the crust of a magnetized neutron star (NS) and its impact on the observable surface temperature. Extending the approach initiated in Geppert et al.(2004), we consider more complex and, hence, more realistic, magnetic field structures but still restrict ourselves to axisymmetric configurations. We put special emphasis on the heat blanketing effect of a toroidal field component. We show that asymmetric temperature distributions can occur and a crustal field consisting of dipolar poloidal and toroidal components will cause one polar spot to be larger than the opposing one. These two warm regions can be separated by an extended cold equatorial belt. We present an internal magnetic field structure which can explain both the X-ray and optical spectra of the isolated NS RXJ 1856-3754. We investigate the effects of the resulting surface temperature profiles on the observable lightcurve which an isolated thermally emitting NS would produce for different field geometries. The lightcurves will be both qualitatively (deviations from sinusoidal shape) and quantitatively (larger pulsed fraction for the same observational geometry) different from those of a NS with an isothermal crust. This opens the possibility to determine the interna} magnetic field strengths and structures in NSs by modeling their X-ray lightcurves and spectra. The striking similarities of our model calculations with the observed spectra and pulse profiles of isolated thermally emitting NSs is an indication for the existence of strong magnetic field components maintained by crustal currents.

Current observations of the fraction of dark energy and a lower limit on its tension, coupled with an assumption of the non-convexity of the dark energy potential, are used to derive a lower limit of 26 billion years for the future age of the universe. Conversely, our ordered observations, coupled with an assumption that observers are smaller than the universe, are used to argue for an upper limit of about e^10^50 years if the universe eventually undergoes power-law expansion, and an upper limit of only about 10^60 years left for our universe if it continues to expand exponentially at the current rate.

We present a calculation of neutrino pair bremsstrahlung due to electron-electron scattering in a relativistic degenerate plasma of electrons. Proper treatment of the in-medium photon propagator, i.e., inclusion of Debye screening of the longitudinal part and Landau damping of the transverse part, leads to a neutrino emissivity which is several orders of magnitude larger than when Debye screening is imposed for the tranverse part. Our results show that this in-medium process can compete with other sources of neutrino radiation and can, in some cases, even be the dominant neutrino emission mechanism. We also discuss the natural extension to quark-quark bremsstrahlung in gapped and ungapped quark matter.

Superbursts are very energetic Type I X-ray bursts discovered in recent years by long term monitoring of X-ray bursters, believed to be due to unstable ignition of carbon in the deep ocean of the neutron star. A number of "intermediate duration" bursts have also been observed, probably associated with ignition of a thick helium layer. We investigate the sensitivity of these long X-ray bursts to the thermal profile of the neutron star crust and core. We first compare cooling models of superburst lightcurves with observations, and derive constraints on the ignition mass and energy release, and then calculate ignition models for superbursts and pure helium bursts, and compare to observations. The superburst lightcurves and ignition models imply that the carbon mass fraction is approximately 20% or greater in the fuel layer, constraining models of carbon production. However, the most important result is that when Cooper pairing neutrino emission is included in the crust, the temperature is too low to support unstable carbon ignition at the observed column depths. Some additional heating mechanism is required in the accumulating fuel layer to explain the observed properties of superbursts. If Cooper pair emission is less efficient than currently thought, the observed ignition depths for superbursts imply that the crust is a poor conductor, and the core neutrino emission is not more efficient than modified URCA. The observed properties of helium bursts support these conclusions, requiring inefficient crust conductivity and core neutrino emission.

Recent models of carbon ignition on accreting neutron stars predict superburst ignition depths that are an order of magnitude larger than observed. We explore a possible solution to this problem, that the compact stars in low mass X-ray binaries that have shown superbursts are in fact strange stars with a crust of normal matter. We calculate the properties of superbursts on strange stars, and the resulting constraints on the properties of strange quark matter. We show that the observed ignition conditions exclude fast neutrino emission in the quark core, for example by the direct Urca process, which implies that strange quark matter at stellar densities should be in a color superconducting state. For slow neutrino emission in the quark matter core, we find that reproducing superburst properties requires a definite relation between three poorly constrained properties of strange quark matter: its thermal conductivity, its slow neutrino emissivity and the energy released by converting a nucleon into strange quark matter.

The cooling of a compact star depends very sensitively on the state of dense matter at supranuclear densities, which essentially controls the neutrino emission, as well as on the structure of the stellar outer layers which control the photon emission. Open issues concern the hyperon population, the presence of meson condensates, superfluidity and superconductivity, and the transition of confined hadronic matter to quark matter. This paper describes these issues and presents cooling calculations based on a broad collection of equations of state for neutron star matter and strange matter. These results are tested against the body of observed cooling data.

We present a detailed, two dimensional numerical study of the microphysical conditions and dynamical evolution of accretion disks around black holes when neutrino emission is the main source of cooling. Such structures are likely to form after the gravitational collapse of massive rotating stellar cores, or the coalescence of two compact objects in a binary (e.g., the Hulse--Taylor system). The physical composition is determined self consistently by considering two regimes: neutrino--opaque and neutrino--transparent, with a detailed equation of state which takes into account neutronization, nuclear statistical equilibrium of a gas of free nucleons and alpha particles, blackbody radiation and a relativistic Fermi gas of arbitrary degeneracy. Various neutrino emission processes are considered, with electron/positron capture onto free nucleons providing the dominant contribution to the cooling rate. We find that important temporal and spatial scales, related to the optically thin/optically thick transition are present in the disk, and manifest themselves clearly in the energy output in neutrinos. This transition produces an inversion of the lepton gradient in the innermost regions of the flow which drives convective motions, and affects the density and disk scale height radial profiles. The electron fraction remains low in the region close to the black hole, and if preserved in an outflow, could give rise to heavy element nucleosynthesis. Our specific initial conditions arise from the binary merger context, and so we explore the implications of our results for the production of gamma ray bursts.

An inexhaustive review of Hawking radiation and black hole thermodynamics is given, focusing especially upon some of the historical aspects as seen from the biased viewpoint of a minor player in the field on and off for the past thirty years.

A new classification of neutron star cooling scenarios, involving either ``minimal'' cooling or ``enhanced'' cooling is proposed. The minimal cooling scenario replaces and extends the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation process. This emission dominates that due to the modified Urca process for temperatures close to the critical temperature for superfluid pairing. Minimal cooling is distinguished from enhanced cooling by the absence of neutrino emission from any direct Urca process, due either to nucleons or to exotica. Within the minimal cooling scenario, theoretical cooling models can be considered to be a four parameter family involving the equation of state of dense matter, superfluid properties of dense matter, the composition of the neutron star envelope, and the mass of the neutron star. Consequences of minimal cooling are explored through extensive variations of these parameters. Results are compared with the inferred properties of thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. All stars for which thermal emissions have been clearly detected are at least marginally consistent with the lack of enhanced cooling. The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced cooling in case their ages and/or temperatures are on the lower side of their estimated values whereas the four stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 may require some source of internal heating in case their age and/or luminosity are on the upper side of their estimated values. The new upper limits on the thermal luminosity of PSR J0205+6449 and RX J0007.0+7302 are indicative of the occurrence of some enhanced neutrino emission beyond the minimal scenario.

Observations to date cannot distinguish neutron stars from self-bound bare quark stars on the basis of their gross physical properties such as their masses and radii alone. However, their surface luminosity and spectral characteristics can be significantly different. Unlike a normal neutron star, a bare quark star can emit photons from its surface at super-Eddington luminosities for an extended period of time. We present a calculation of the photon bremsstrahlung rate from the bare quark star's surface, and indicate improvements that are required for a complete characterization of the spectrum. The observation of this distinctive photon spectrum would constitute an unmistakable signature of a strange quark star and shed light on color superconductivity at stellar densities.

A general overview of the main physical processes driving the cooling of an isolated neutron star is presented. Among the most important ones are the various possible neutrino emission processes and the occurrence of baryon pairing. Special emphasis is also put on the importance of the chemical composition of the upper layers of the star. A detailed analysis of a Minimal Scenario, which explicitly postulates that no ``exotic'' form of matter be present, is summarized and compared with presently available observational data. No striking incompatibility of the data with the predictions of the Minimal Scenario is found. Nevertheless, two, possibly three, conspicuous stars are identified which may, when better data are available, constitute strong astrophysical evidence for the occurrence of a new state of matter at high density.

Cosmological gamma ray bursts (GRBs) are thought to occur from violent hypercritical accretion onto stellar mass black holes, either following core collapse in massive stars or compact binary mergers. This dichotomy may be reflected in the two classes of bursts having different durations. Dynamical calculations of the evolution of these systems are essential if one is to establish characteristic, relevant timescales. We show here for the first time the result of dynamical simulations, lasting approximately one second, of post--merger accretion disks around black holes, using a realistic equation of state and considering neutrino emission processes. We find that the inclusion of neutrino optical depth effects produces important qualitative temporal and spatial transitions in the evolution and structure of the disk, which may directly reflect upon the duration and variability of short GRBs.

We investigate the influence of different magnetic field configurations on the temperature distribution in neutron star crusts. We consider axisymmetric dipolar fields which are either restricted to the stellar crust, ``crustal fields'', or allowed to penetrate the core, ``core fields''. By integrating the two-dimensional heat transport equation in the crust, taking into account the classical (Larmor) anisotropy of the heat conductivity, we obtain the crustal temperature distribution, assuming an isothermal core. Including quantum magnetic field effects in the envelope as a boundary condition, we deduce the corresponding surface temperature distributions. We find that core fields result in practically isothermal crusts unless the surface field strength is well above $10^{15}$ G while for crustal fields with surface strength above a few times $10^{12}$ G significant deviations from isothermality occur at core temperatures inferior or equal to $10^8$ K. At the stellar surface, the cold equatorial region produced by the quantum suppression of heat transport perpendicular to the field in the envelope, present for both core and crustal fields, is significantly extended by the classical suppression at higher densities in the case of crustal fields. This can result, for crustal fields, in two small warm polar regions which will have observational consequences: the neutron star has a small effective thermally emitting area and the X-ray pulse profiles are expected to have a distinctively different shape compared to the case of a neutron star with a core field. These features, when compared with X-ray data on thermal emission of young cooling neutron stars, will open a way to provide observational evidence in favor, or against, the two radically different configurations of crustal or core magnetic fields.

The photon emissivity from the bremsstrahlung process ee-> ee\gamma occuring in the electrosphere at the bare surface of a strange quark star is calculated. For surface temperatures T<10^9K, the photon flux exceeds that of e+e- pairs that are produced via the Schwinger mechanism in the presence of a strong electric field that binds electrons to the surface of the quark star. The average energy of photons emitted from the bremsstrahlung process can be 0.5 MeV or more, which is larger than that in e+e- pair annihilation. The observation of this distinctive photon spectrum would constitute an unmistakable signature of a strange quark star and shed light on color superconductivity at stellar densities.

By combining a renormalization group argument relating the charge e and mass m of the proton by e^2 ln m ~ 0.1 pi (in Planck units) with the Carter-Carr-Rees anthropic argument that gives an independent approximate relation m ~ e^20 between these two constants, both can be crudely estimated. These equations have the factor of 0.1 pi and the exponent of 20 which depend upon known discrete parameters (e.g., the number of generations of quarks and leptons, and the number of spatial dimensions), but they contain NO continuous observed parameters. Their solution gives the charge of the proton correct to within about 8%, though the mass estimate is off by a factor of about 1000 (16% error on a logarithmic scale). When one adds a fudge factor of 10 previously given by Carr and Rees, the agreement for the charge is within about 2%, and the mass is off by a factor of about 3 (2.4% error on a logarithmic scale). If this 10 were replaced by 15, the charge agrees within 1.1% and the mass itself agrees within 0.7%.

Recent Chandra observations have raised expectations that the objects RX J185635-3754 and 3C58 are either bare strange quark stars or stars with extended quark cores. However, these observations can also be interpreted in terms of normal neutron stars. Essential requirements for either explanation are that they simultaneously account for the observed (i) spectral features (i.e., thermal or nonthermal, possible lines), (ii) bounds on the inferred mass (M) and radius (R), and (iii) the cooling curves (effective temperatures T and luminosities L vs. age). Compact stars offer the promise of delineating QCD at finite baryon density (Lattimer and Prakash, 2001) in a fashion complementary to relativistic heavy ion collider experiments, which largely address the finite temperature, but baryon poor regime.

The cooling of a young bare strange star is studied numerically by solving the equations of energy conservation and heat transport for both normal and superconducting strange quark matter inside the star. We show that the thermal luminosity from the strange star surface, due to both photon emission and e+e- pair production, may be orders of magnitude higher than the Eddington limit, for about one day for normal quark matter but possibly for up to a hundred years for superconducting quark matter, while the maximum of the photon spectrum is in hard X-rays with a mean energy of ~ 100 keV or even more. This differs both qualitatively and quantitatively from the photon emission from young neutron stars and provides a definite observational signature for bare strange stars. It is shown that the energy gap of superconducting strange quark matter may be estimated from the light curves if it is in the range from ~ 0.5 MeV to a few MeV.

We explore, using an exact cooling code, the thermal evolution of a neutron star undergoing episodes of intense accretion, alternated by long periods of quiescence (e.g. Soft X-Ray Transients). We find that the soft component of the quiescent luminosity of Aql X-1, 4U 1608-522 and of SAX J1808.4-3658 can be understood as thermal emission from a cooling neutron star with negligible neutrino emission. In the case of Cen X-4 strong neutrino emission from the inner core is necessary to explain the observation: this may indicate that the neutron star of Cen X-4 is heavier than 1.4 Msun. This study opens the possibility of using the quiescent emission of Soft X-Ray Transients as a tool for probing the core superfluidity in relation to the mass of the neutron star.

If one surrounds a black hole with a perfectly reflecting shell and adiabatically squeezes the shell inward, one can increase the black hole area A to exceed four times the total entropy S, which stays fixed during the process. A can be made to exceed 4S by a factor of order unity before the one enters the Planck regime where the semiclassical approximation breaks down. One interpretation is that the black hole entropy resides in its thermal atmosphere, and the shell restricts the atmosphere so that its entropy is less than A/4.

We explore the thermal evolution of a neutron star undergoing episodes of intense accretion, separated by long periods of quiescence. By using an exact cooling code we follow in detail the flow of heat in the star due to the time-dependent accretion-induced heating from pycno-nuclear reactions in the stellar crust, to the surface photon emission, and the neutrino cooling. These models allow us to study the neutron stars of the Soft X-Ray Transients. In agreement with Brown, Bildsten and Rutledge (1998) we conclude that the soft component of the quiescent luminosity of Aql X-1, 4U 1608-522, and of the recently discovered SAX J1808.4, can be understood as thermal emission from a cooling neutron star with negligible neutrino emission. However, we show that in the case of Cen X-4, despite its long recurrence time, strong neutrino emission from the neutron star inner core is necessary to understand the observed low ratio of quiescent to outburst luminosity. This result implies that the neutron star in Cen X-4 is heavier than the one in the other systems and the pairing critical temperature in its center must be low enough (well below 10^9 K) to avoid a strong suppression of the neutrino emission.

We investigate the thermal, magnetic and rotational evolution of isolated neutron stars assuming that the dipolar magnetic field is confined to the crust. Our treatment, for the first time, uses a fully general relativistic formalism not only for the thermal but also for the magnetic part, and includes partial general relativistic effects in the rotational part. Due to the fact that the combined evolution depends crucially upon the compactness of the star, three different equations of state have been employed in the calculations. In the absence of general relativistic effects, while upon increasing compactness a decrease of the crust thickness takes place leading into an accelerating field decay, the inclusion of general relativistic effects intend to ``decelerate this acceleration''. As a consequence we find that within the crustal field hypothesis, a given equation of state is compatible with the observed periods $P$ and period derivative $\dot P$ provided the initial field strength and current location as well as the magnitude of the impurity content are constrained appropriately. Finally, we access the flexibility of the soft, medium and stiff classes of equations of state as candidates in describing the state of the matter in the neutron star interiors. The comparison of our model calculations with observations, together with the consideration of independent information about neutron star evolution, suggests that a not too soft equation of state describes neutron star interiors and its cooling proceeds along the `standard' scenario.

An analysis of the role of general relativistic effects on the decay of neutron star's magnetic field is presented. At first, a generalized induction equation on an arbitrary static background geometry has been derived and, secondly, by a combination of analytical and numerical techniques, a comparison of the time scales for the decay of an initial dipole magnetic field in flat and curved spacetime is discussed. For the case of very simple neutron star models, rotation not accounted for and in the absence of cooling effects, we find that the inclusion of general relativistic effects result, on the average, in an enlargement of the decay time of the field in comparison to the flat spacetime case. Via numerical techniques we show that, the enlargement factor depends upon the dimensionless compactness ratio ${\epsilon}={2GM \over c^{2}R}$, and for ${\epsilon}$ in the range $(0.3~,~0.5)$, corresponding to compactness ratio of realistic neutron star models, this factor is between 1.2 to 1.3. The present analysis shows that general relativistic effects on the magnetic field decay ought to be examined more carefully than hitherto. A brief discussion of our findings on the impact of neutron stars physics is also presented.

Baryon and quark superfluidity in the cooling of neutron stars are investigated. Observations could constrain combinations of the neutron or Lambda-hyperon pairing gaps and the star's mass. However, in a hybrid star with a mixed phase of hadrons and quarks, quark gaps larger than a few tenths of an MeV render quark matter virtually invisible for cooling. If the quark gap is smaller, quark superfluidity could be important, but its effects will be nearly impossible to distinguish from those of other baryonic constituents.

We confront theoretical models for the rotational, magnetic, and thermal evolution of an ultra-magnetized neutron star, or magnetar, with available data on the Anomalous X-ray Pulsars (AXPs). We argue that, if the AXPs are interpreted as magnetars, their clustering of spin periods between 6 and 12 seconds, observed at present in this class of objects, their period derivatives, their thermal X-ray luminosities, and the association of two of them with young supernova remnants, can only be understood globally if the magnetic field in magnetars decays significantly on a time scale of the order of 10^4 years.

Soft Gamma-ray Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are interpreted as young highly magnetized neutron stars (NSs). Their X-ray luminosity in quiescence, exceeding 10^{35} erg s^{-1} cannot be explained as due to cooling of a highly magnetized NS, but requires as an extra heat source the decay of its magnetic field (MF). We study numerically the coupled evolution of the MF, temperature and spin period under the assumption that the currents maintaining the field are confined in the crust of the star. The decay of the field depends on the field strength itself (Hall-drift), on the temperature and injects heat into the star, but is controlled by neutrino emission. Finally we consider the spin down from magnetic dipole braking with this decaying field to track the long term evolution. We find reasonable initial conditions for the MF strength and structure to explain their current observational values both of their rotational period, its time derivative and the X-ray luminosity of AXPs and SGRs.the X-ray luminosity of AXPs and SGRs.

After a brief summary of neutron star cooling theory I present results which emphasize the importance of baryon pairing in the neutron star core. I show how the thermal evolution may be totally controlled by pairing for models which include only nucleons and models containing also hyperons. Finally, I consider the thermal evolution of ultramagnetized neutron stars whose existence has been inferred from the soft gamma repeaters.

I give a general overview of the theory of neutron star cooling, emphasizing the intuitive understanding of the effects of the various physical ingredients, including the recently proposed accreted envelopes and neutrino emission from Cooper pair breaking and formation. I describe how the present data may be compatible with both the `standard model' and the fast cooling scenarios, with or without invoking extensive presence of baryon pairing.

It is generally considered that the neutron star cooling scenarios involving fast neutrino emission, from a kaon or pion condensate, quark matter, or the direct Urca process, require the presence of baryon pairing in the central core of the star to control the strong neutrino emission and produce surface temperatures compatible with observations. I show here that within the kaon condensate scenario pairing is not necessary if: 1) the equation of state is stiff enough for the star to have a thick crust in which sufficient friction can occur to heat the star and 2) a thin layer, of mass larger than 10^{-12} Msol, of light elements (H and He) is present at the stellar surface. The effect of the light elements is to increase the heat flow and thus produce a higher surface temperature. Both the occurrence of heating and the presence of H and/or He at the surface (deposited during the late post-supernova accretion) can possibly be confirmed or infirmed by future observations.

Witten's conjecture about strange quark matter (`Strange Matter') being the ground state of QCD at finite baryon number is presented and stars made of strange matter (`Strange Stars') are compared to neutron stars. The only observable way in which a strange star differs from a neutron star is in its early thermal history and a detailed study of strange star cooling is reported and compared to neutron star cooling. One concludes that future detection of thermal radiation from the compact object produced in the core collapse of SN 1987A could present the first evidence for strange matter.

We complete our study of pulsars' non-uniform surface temperature and of its effects on their soft X-ray thermal emission. Our previous work had shown that, due to gravitational lensing, dipolar fields cannot reproduce the strong pulsations observed in Vela, Geminga, PSR 0656+14, and PSR 1055-52. Assuming a standard neutron star mass of 1.4 Msol, we show here that the inclusion of a quadrupolar component, if it is suitably oriented, is sufficient to increase substantially the pulsed fraction, Pf, up to, or above, the observed values if the stellar radius is 13 km or even 10 km. For models with a radius of 7 km the maximum pulsed fraction obtainable with (isotropic) blackbody emission is of the order of 15% for orthogonal rotators (Vela, Geminga and PSR 1055-52) and only 5% for an inclined rotator as PSR 0656+14. Given the observed values this indicates that the neutron stars in Geminga and PSR 0656+14 have radii significantly larger than 7 km and, given the very specific quadrupole components required to increase Pf, even radii of the order of 10 km may be unlikely in all four cases. We confirm our previous finding that the pulsed fraction always increases with photon energy, below about 1 keV, when blackbody emission is used and show that it is due to the hardenning of the blackbody spectrum with increasing temperature. The observed decrease of pulsed fraction may thus suggest that the emitted spectrum softens with increasing temperature. Finally, we apply our model to reassess the magnetic field effect on the outer boundary condition used in neutron star cooling models and show that, in contradistinction to several previous claims, it is very small and most probably results in a slight reduction of the heat flow through the envelope.

We show that the expected inhomogeneous temperature distribution induced at the surface of a neutron star by the anisotropy of heat transport in the magnetized envelope allows us to understand quite well the observed pulse profiles of the four nearby pulsars for which surface thermal emission has been detected. However, due to gravitational lensing, dipolar magnetic fields are not adequate and the observed high pulsed fractions force us to include a quadrupolar component.

We use models of magnetized hydrogen atmospheres to fit the soft X-ray thermal spectrum of the Vela pulsar obtained by ROSAT. The distance and hydrogen column density required by our fits are in good agreement with other independents estimates and can be considered as an argument for the presence of hydrogen at the surface of this pulsar. The low temperature obtained, T = 7.85 +/- 0.25 10^5 K, is below the predictions of the `standard' model of neutron star cooling and strengthens previous claims for the necessity of the occurrence of fast neutrino cooling in this neutron star.

After Michel (1994) introduced a phenomenological picture of `rapid magnetization' of newly born neutron stars (NSs), Muslimov & Page (1995) suggested that the physical conditions accompanying the formation of a NS may result in the surface magnetic field of the NS having a low value (~ 10^8 - 10^9 G) while the bulk of the flux is submerged under the stellar surface. The subsequent growth of the surface magnetic field occurs owing to ohmic diffusion of a strong internal field (10^12-10^13 G) and proceeds over a time scale of 10^2-10^3 yrs, depending on the early thermal history of the NS, the initial distribution of magnetic flux and the electrical resistivity in the stellar crust, as well as on the details of the stellar structure. Referring to this suggestion, we perform numerical computations to demonstrate its possible relevance to young radio pulsars. In our calculations, we use different NS structures based on a model of dense matter presenting a phase transition to kaon condensation which implies both slow and fast neutrino cooling, depending on the stellar mass. We present sequences of secular evolution of the surface magnetic field, spin-down luminosity, spin period, `braking index', and spin-down age best matching the corresponding measured quantities for PSR 0531+21 (Crab), PSR 0540-69, and PSR 1509-58. We conclude that the effect under consideration reproduces remarkably well both the observed rotational characteristics and derived values of the surface magnetic field strength for these radio pulsars, the only ones having reliable measurements of the braking indices. In addition, our analysis suggests that the NSs in the Crab and PSR 0540-69 have experienced fast neutrino cooling while the NS in PSR 1509-58 has undergone slow cooling and hence has a mass below the critical mass.

We present a model to explain the decrease in the amplitude of the pulse profile with increasing energy observed in Geminga's soft X-ray surface thermal emission. We assume the presence of plates surrounded by a surface with very distinct physical properties: these two regions emit spectra of very distinct shapes which present a crossover, the warm plates emitting a softer spectrum than the colder surrounding surface. The strongly pulsed emission from the plates dominates at low energy while the surroundings emission dominates at high energy, producing naturally a strong decrease in the pulsed fraction. In our illustrative example the plates are assumed to be magnetized while the rest of the surface is field free. This plate structure may be seen as a schematic representation of a continuous but very nonuniform distribution of the surface magnetic field or as a quasi realistic structure induced by past tectonic activity on Geminga.

We consider, in general terms, the early thermal evolution of an isolated neutron star, i.e., during the first $10^5$ years after the supernova explosion when the cooling is driven by neutrino emission from the core. It is shown that, if fast neutrino emission is occuring, the surface temperature is actually determined by the critical temperature $T_c$ at which the core baryons become superfluid, and almost nothing else. In particular, the dependence on the actual neutrino emission process is very weak. In other terms, fast neutrino cooling neutron stars are thermometers for the highest $T_c$ superfluid (superconductor) in the Universe. The observational data are briefly presented.

Surface thermal emission has been detected by ROSAT from four nearby young neutron stars. Assuming black body emission, the significant pulsations of the observed light curves can be interpreted as due to large surface temperature differences produced by the effect of the crustal magnetic field on the flow of heat from the hot interior toward the cooler surface. However, the energy dependence of the modulation observed in Geminga is incompatible with blackbody emission: this effect will give us a strong constraint on models of the neutron star surface.

We model the temperature distribution at the surface of a magnetized neutron star and study the effects on the observable X-ray spectra and light curves, taking fully into account gravitational lensing. Only dipolar fields are considered here. General features are studied and compared with the ROSAT data from the pulsars 0833-45 (Vela), 0656+14, 0630+178 (Geminga), and 1055-52. The composite blackbody spectra obtained are not very different from a black- body spectrum at the star's effective temperature: temperature estimates using single temperature models are practically equivalent to composite model results. Gravitational lensing strongly suppresses the light curve pulsations: this may be a new tool to constrain the size of neutron stars. Dipolar surface fields yield many features somewhat similar to the observed ones but are not sufficient for an accurate description. The pulsed fractions obtained in all our models increase with photon energy: the strong decrease observed in Geminga at energies 0.3 - 0.5 keV is definitely a genuine effect of the atmospheric magnetic field on the locally emitted spectrum in contradistinction to the magnetic effects on the surface temperature. A complete model will have to consider complex surface field structures and include effects of the atmospheric magnetic field on the emitted spectrum.

We show that Geminga's surface temperature can be understood by both types of neutrino cooling scenarios, i.e, slow neutrino cooling by the modified Urca process or fast neutrino cooling by the direct Urca process or by some exotic matter, and thus does not allow us to discriminate between these two competing schemes. However, for both types of scenarios, agreement with the observed temperature can only be obtained if BARYON PAIRING IS PRESENT IN MOST, IF NOT ALL, OF THE CORE OF THE STAR. We also comment on the recent temperature estimates for PSR 0656+14 and PSR 1055-52, which pertain to the same photon cooling era. We argue though that observational evidence for the slow neutrino cooling model (the ``standard'' model) is in fact very dim and that the interpretation of the surface temperature of all neutron stars could be done with a reasonable theoretical a priori within the fast neutrino cooling scenarios only. In this case, Geminga, PSR 0656+14, and PSR 1055-52 all show evidence of baryon pairing down to their very centers. xxx----To be published in Ap.J. (June 10 isue).---xxx