The observation of thermal emission from isolated neutron stars and the modeling of the corresponding cooling curves has been very useful to get information on the properties of matter at very high densities. More recently, the detection of quiescent thermal emission from neutron stars in low mass X-ray binary systems after active periods opened a new window to the physics of matter at lower densities. Here we analyze a few sources that have been recently monitored and we show how the models can be used to establish constraints on the crust composition and their transport properties, depending on the astrophysical scenarios assumed.

We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in x-rays (outbursts). The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659-29, KS 1731-260, EXO 0748-676, XTE J1701-462 and IGR J17480-2446. We find that the evolution of MXB 1659-29, KS 1731-260 and EXO 0748-676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701-462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480-2446 with an additional heat deposition in the outer crust from shallow sources.

In the magneto-centrifugal mechanism for jet formation, accreting neutron stars are assumed to produce relativistic jets only if their surface magnetic field is weak enough ($B \sim 10^8$ G). However, the most common manifestation of neutron stars are pulsars, whose magnetic field distribution peaks at $B \sim 10^{12}$ G. If the neutron star magnetic field has at least this strength at birth, it must decay considerably before jets can be launched in binary systems. We study the magnetic field evolution of a neutron star that accretes matter from the wind of a high-mass stellar companion so that we can constrain the accretion rate and the impurities in the crust, which are necessary conditions for jet formation. We solved the induction equation for the diffusion and convection of the neutron star magnetic field confined to the crust, assuming spherical accretion in a simpliflied one-dimensional treatment. We incorporated state-of-the-art microphysics, including consistent thermal evolution profiles, and assumed two different neutron star cooling scenarios based on the superfluidity conditions at the core. We find that in this scenario, magnetic field decay at long timescales is governed mainly by the accretion rate, while the impurity content and thermal evolution of the neutron star play a secondary role. For accretion rates $\dot{M}\gtrsim10^{-10}$ M$_{\odot}$ yr$^{-1}$, surface magnetic fields can decay up to four orders of magnitude in $\sim$10$^7$ yr, which is the timescale imposed by the evolution of the high-mass stellar companion in these systems. Based on these results, we discuss the possibility of transient jet-launching in strong wind-accreting high-mass binary systems like supergiant fast X-ray transients.

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully-coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

We report on a new mechanism for heat conduction in the neutron star crust. We find that collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat conduction in magnetized neutron stars. They can dominate the heat conduction transverse to magnetic field when the magnetic field $B \gsim 10^{13}$ G. At density $\rho \simeq 10^{12}-10^{14} $ g/cm$^3$ the conductivity due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron conductivity when temperature $\simeq 10^8$ K. This new mode of heat conduction can limit the surface anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with and without superfluid heat conduction could show observationally discernible differences.

Context: Many thermally emitting isolated neutron stars have magnetic fields larger than 10^13 G. A realistic cooling model that includes the presence of high magnetic fields should be reconsidered. Aims: We investigate the effects of anisotropic temperature distribution and Joule heating on the cooling of magnetized neutron stars. Methods: The 2D heat transfer equation with anisotropic thermal conductivity tensor and including all relevant neutrino emission processes is solved for realistic models of the neutron star interior and crust. Results: The presence of the magnetic field affects significantly the thermal surface distribution and the cooling history during both, the early neutrino cooling era and the late photon cooling era. Conclusions: There is a large effect of the Joule heating on the thermal evolution of strongly magnetized neutron stars. Both magnetic fields and Joule heating play a key role in keeping magnetars warm for a long time. Moreover, this effect is important for intermediate field neutron stars and should be considered in radio-quiet isolated neutron stars or high magnetic field radio-pulsars.

We present two-dimensional simulations for the cooling of neutron stars with strong magnetic fields (B > 1e13 Gauss). We study how the cooling curves are influenced by magnetic field decay. We show that the Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field and predict the surface temperature distribution for the early and late stages of the evolution of isolated neutron stars, comparing our results with available observational data of isolated neutron stars.

The impact of strong magnetic fields B>10e13 G on the thermal evolution of neutron stars is investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D cooling simulations with anisotropic thermal conductivity considering all relevant neutrino emission processes for realistic neutron stars. The standard cooling models of neutron stars are called into question by showing that the magnetic field has relevant (and in many cases dominant) effects on the thermal evolution. The presence of the magnetic field significantly affects the thermal surface distribution and the cooling history of these objects during both, the early neutrino cooling era and the late photon cooling era. The minimal cooling scenario is thus more complex than generally assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain the observations.

We present 2D simulations of the cooling of neutron stars with strong magnetic fields (B \geq 10^{13} G). We solve the diffusion equation in axial symmetry including the state of the art microphysics that controls the cooling such as slow/fast neutrino processes, superfluidity, as well as possible heating mechanisms. We study how the cooling curves depend on the the magnetic field strength and geometry. Special attention is given to discuss the influence of magnetic field decay. We show that Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field for the early and late stages of the evolution of isolated neutron stars.

We consider the color spin locking (CSL) phase of two-flavor quark matter at zero temperature for nonlocal instantaneous, separable interactions. We employ a Lorentzian-type form factor allowing a parametric interpolation between the sharp (Nambu-Jona-Lasinio (NJL) model) and very smooth (e.g. Gaussian) cut-off models for systematic studies of the influence on the CSL condensate the deviation from the NJL model entails. This smoothing of the NJL model form factor shows advantageous features for the phenomenology of compact stars: (i) a lowering of the critical chemical potential for the onset of the chiral phase transition as a prerequisite for stability of hybrid stars with extended quark matter cores and (ii) a reduction of the smallest pairing gap to the order of 100 keV, being in the range of values interesting for phenomenological studies of hybrid star cooling evolution.

We present results of a microscopic calculation using NJL-type model of possible spin-one pairings in two flavor quark matter for applications in compact star phenomenology. We focus on the color-spin locking phase (CSL) in which all quarks pair in a symmetric way, in which color and spin states are locked. The CSL condensate is particularly interesting for compact star applications since it is flavor symmetric and could easily satisfy charge neutrality. Moreover, the fact that in this phase all quarks are gapped might help to suppress the direct Urca process, consistent with cooling models. The order of magnitude of these small gaps (~1 MeV) will not influence the EoS, but their also small critical temperatures (T_c ~800 keV) could be relevant in the late stages neutron star evolution, when the temperature falls below this value and a CSL quark core could form.

We present results for the spin-1 color-spin locking phase (CSL) using a NJL-type model in two flavor quark matter for compact stars applications. The CSL condensate is flavor symmetric and therefore charge and color neutrality can easily be satisfied. We find small energy gaps ~1 MeV, which make the CSL matter composition and the EoS not very different from the normal quark matter phase. We keep finite quark masses in our calculations and obtain no gapless modes that could have strong consequences in the late cooling of neutron stars. Finally, we show that the region of the phase diagram relevant for neutron star cores, when asymmetric flavor pairing is suppressed, could be covered by the CSL phase.

We study the influence of nonlocality in the interaction on two spin one pairing patterns of two-flavor quark matter: the anisotropic blue color paring besides the usual two color superconducting matter (2SCb), in which red and green colors are paired, and the color spin locking phase (CSL). The effect of nonlocality on the gaps is rather large and the pairings exhibit a strong dependence on the form factor of the interaction, especially in the low density region. The application of these small spin-one condensates for compact stars is analyzed: the early onset of quark matter in the nonlocal models may help to stabilize hybrid star configurations. While the anisotropic blue quark pairing does not survive a big asymmetry in flavor space as imposed by the charge neutrality condition, the CSL phase as a flavor independent pairing can be realized as neutral matter in compact star cores. However, smooth form factors and the missmatch between the flavor chemical potential in neutral matter make the effective gaps of the order of magnitude $\simeq 10$ keV, and a more systematic analysis is needed to decide whether such small gaps could be consistent with the cooling phenomenology.

We study a spin-1 single flavor color superconducting phase which results from a color-spin locking (CSL) interaction in two-flavor quark matter. This phase is particularly interesting for compact star cooling applications since the CSL phase may survive under charge neutrality constraints implying a mismatch between up- and down-quark chemical potentials which can destroy the scalar diquark condensate. CSL gaps are evaluated within an NJL model and they are found to be consistent with cooling phenomenology if a density dependent coupling constant is used.

The equation of state for quark matter is derived for a nonlocal, chiral quark model within the mean field approximation. Special emphasis is on the occurrence of a diquark condensate which signals a phase transition to color superconductivity and its effects on the equation of state. We present a fit formula for the Bag pressure, which is density dependent in case when the quark matter is color superconducting. We calculate the quark star configurations by solving the Tolman- Oppenheimer- Volkoff equations and demonstrate the effects of diquark condensation on the stability of hybrid stars for different formfactors of the quark interaction.

An equation of state for 2-flavor quark matter (QM) with diquark condensation under the conditions for compact stars -$\beta$-equilibrium, charge and color neutrality- is presented. Trapped antineutrinos prevent the formation of the diquark condensate at moderate densities above a critical value of the antineutrino chemical potential $\mu_{\bar \nu_e}^c$. The following consequences are presented: 1) The star develops a 2-phase structure ($\mu_{\bar \nu_e}\geq \mu_{\bar \nu_e}^c$): a color superconducting QM core and a normal QM shell. 2)During the cooling, when the temperature is small enough ($T<1$ MeV) the antineutrino mean free path becomes larger than the thickness of the normal QM shell and the antineutrinos get untrapped in a sudden burst. The energy release is estimated as $\simeq 10^{52}$ erg and an antineutrino pulse is expected to be observed.

We study the consequences of antineutrino trapping in hot quark matter for quark star configurations with possible diquark condensation. Due to the conditions of charge neutrality and beta-equilibrium the flavor asymmetry increases with the number density of trapped antineutrinos such that above a critical value of the antineutrino chemical potential of 30 MeV diquark condensation is inhibited. When the quark star cools a two-phase structure occurs: a superconducting quark matter core surrounded by a shell of normal quark matter. Below the critical temperature of about 1 MeV, the antineutrino mean free path becomes larger than the thickness of the normal quark matter shell so that they get untrapped within a sudden process. By comparing the masses of configurations with the same baryon number we estimate that the release of energy due to the antineutrino untrapping transition can be in the range of 10^{51} to 10^{52} erg.