Search

We present the discovery of Aquarius III, an ultra-faint Milky Way satellite galaxy identified in the second data release of the DECam Local Volume Exploration (DELVE) survey. Based on deeper follow-up imaging with DECam, we find that Aquarius III is a low-luminosity ($M_V = -2.5^{+0.3}_{-0.5}$; $L_V = 850^{+380}_{-260} \ L_{\odot}$), extended ($r_{1/2} = 41^{+9}_{-8}$ pc) stellar system located in the outer halo ($D_{\odot} = 85 \pm 4$ kpc). From medium-resolution Keck/DEIMOS spectroscopy, we identify 11 member stars and measure a mean heliocentric radial velocity of $v_{\rm sys} = -13.1^{+1.0}_{-0.9} \ \rm km \ s^{-1}$ for the system and place an upper limit of $\sigma_v < 3.5 \rm \ km \ s^{-1}$ ($\sigma_v < 1.6 \rm \ km \ s^{-1}$) on its velocity dispersion at the 95% (68%) credible level. Based on Calcium-Triplet-based metallicities of the six brightest red giant members, we find that Aquarius III is very metal-poor ([Fe/H]$ = -2.61 \pm 0.21$) with a statistically-significant metallicity spread ($\sigma_{\rm [Fe/H]} = 0.46^{+0.26}_{-0.14}$ dex). We interpret this metallicity spread as strong evidence that the system is a dwarf galaxy as opposed to a star cluster. Combining our velocity measurement with $Gaia$ proper motions, we find that Aquarius III is currently situated near its orbital pericenter in the outer halo ($r_{\rm peri} = 78 \pm 7$ kpc) and that it is plausibly on first infall onto the Milky Way. This orbital history likely precludes significant tidal disruption from the Galactic disk, notably unlike other satellites with comparably low velocity dispersion limits in the literature. Thus, if further velocity measurements confirm that its velocity dispersion is truly below $\sigma_v \lesssim 2 \rm \ km \ s^{-1}$, Aquarius III may serve as a useful laboratory for probing galaxy formation physics in low-mass halos., Comment: 25 pages, 7 figures. Submitted to AAS Journals

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems., Comment: In press for Astrophysical Journal

The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of M33 (M_HI(tot) =1.4 x 10^9 Msun) is found beyond the star forming disk as traced in the far-ultraviolet (FUV). The most distinct features are extended warps, an arc from the northern warp to the disk, diffuse gas surrounding the galaxy, and a southern cloud with a filament back to the galaxy. The features extend out to 22 kpc from the galaxy center (18 kpc from the edge of the FUV disk) and the gas is directly connected to M33's disk. Only five discrete clouds (i.e., gas not directly connected to M33 in position-velocity space) are catalogued in the vicinity of M33, and these clouds show similar properties to Galactic and M31 halo clouds. M33's gaseous features most likely originate from the tidal disruption of M33 by M31 1-3 Gyr ago as shown from an orbit analysis which results in a tidal radius < 15 kpc in the majority of M33's possible orbits. M33 is now beyond the disruptive gravitational influence of M31 and the gas appears to be returning to M33's disk and redistributing its star formation fuel. M33's high mean velocity dispersion in the disk (~18.5 km/s) may also be consistent with the previous interaction and high rate of star formation. M33 will either exhaust its star formation fuel in the next few Gyrs or eventually become star formation fuel for M31. The latter represents the accretion of a large gaseous satellite by a spiral galaxy, similar to the Magellanic Clouds' relationship to the Galaxy., Comment: published in ApJ (13 figures, several with color), see published version for full resolution figures; data available at http://sites.google.com/site/galfahi/data

We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and superfluid properties of neutron star cores by comparing the cooling theory with observations of thermal radiation from isolated neutron stars., Comment: 9 pages, 4 figures, 3 tables, to appear in the proceedings of "40 Years of Pulsars" held in Montreal, Canada, August 12-17, 2007, eds. C. Bassa, Z. Wang, A. Cumming, V. Kaspi, AIP, in press (v.2 - minor bibliography corrections)

We present five new satellites of the Milky Way discovered in Sloan Digital Sky Survey (SDSS) imaging data, four of which were followed-up with either the Subaru or the Isaac Newton Telescopes. They include four probable new dwarf galaxies -- one each in the constellations of Coma Berenices, Canes Venatici, Leo and Hercules -- together with one unusually extended globular cluster, Segue 1. We provide distances, absolute magnitudes, half-light radii and color-magnitude diagrams for all five satellites. The morphological features of the color-magnitude diagrams are generally well described by the ridge line of the old, metal-poor globular cluster M92. In the last two years, a total of ten new Milky Way satellites with effective surface brightness mu_v >~ 28 mag/sq. arcsec have been discovered in SDSS data. They are less luminous, more irregular and appear to be more metal-poor than the previously-known nine Milky Way dwarf spheroidals. The relationship between these objects and other populations is discussed. We note that there is a paucity of objects with half-light radii between ~40 pc and ~ 100 pc. We conjecture that this may represent the division between star clusters and dwarf galaxies., Comment: 10 pages, 8 figures, submitted to the Astrophysical Journal

In this Letter, we study a localized stellar overdensity in the constellation of Ursa Major, first identified in Sloan Digital Sky Survey (SDSS) data and subsequently followed up with Subaru imaging. Its color-magnitude diagram (CMD) shows a well-defined sub-giant branch, main sequence and turn-off, from which we estimate a distance of ~30 kpc and a projected size of ~250 x 125 pc. The CMD suggests a composite population with some range in metallicity and/or age. Based on its extent and stellar population, we argue that this is a previously unknown satellite galaxy of the Milky Way, hereby named Ursa Major II (UMa II) after its constellation. Using SDSS data, we find an absolute magnitude of M_V \~ -3.8, which would make it the faintest known satellite galaxy. UMa II's isophotes are irregular and distorted with evidence for multiple concentrations; this suggests that the satellite is in the process of disruption., Comment: 4 pages, 4 figures, submitted to ApJ Letters

We study thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in a spherical internal layer. We explore the location of this layer as well as the heating rate that could explain high observable thermal luminosities of magnetars and would be consistent with the energy budget of neutron stars. We conclude that the heat source should be located in an outer magnetar's crust, at densities rho < 5e11 g/cm^3, and should have the heat intensity of the order of 1e20 erg/s/cm^3. Otherwise the heat energy is mainly emitted by neutrinos and cannot warm up the surface., Comment: 8 pages, 5 figures, submitted to MNRAS

The host galaxies of the five local, z<0.25, long-duration gamma-ray bursts, each of which had a well-documented associated supernova, are all faint and metal-poor compared to the population of local star-forming galaxies. We quantify this statement by using a previous analysis of star-forming galaxies (z<0.2) from the Sloan Digital Sky Survey to estimate the fraction of local star formation as a function of host galaxy oxygen abundance. We find that only a small fraction (<25%) of current star formation occurs in galaxies with oxygen abundance 12+log(O/H)<8.6, i.e., about half that of the Milky Way. However, all five low-z GRB hosts have oxygen abundance below this limit, in three cases very significantly so. If GRBs traced local star formation independent of metallicity, the probability of obtaining such low abundances for all five hosts would be P~0.1%. We conclude that GRBs trace only low-metallicity star formation, and that the Milky Way has been too metal rich to host long GRBs for at least the last several billion years. This result has implications for the potential role of GRBs in mass extinctions, for searches for recent burst remnants in the Milky Way and other large galaxies, for non-detections of late radio emission from local core-collapse supernovae, and for the production of cosmic rays in the local Universe. We also find that the isotropic energy release of these five GRBs, E_iso, steeply decreases with increasing host oxygen abundance, suggesting an upper metallicity limit for ``cosmological'' GRBs at ~0.15 Z_solar. (Abridged), Comment: To be published in Acta Astronomica. 16 pages, 3 figures; several small changes

We study the "minimal" cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and the enhanced cooling is produced by the neutrino emission due to Cooper pairing of neutrons. Extending our previous consideration (Gusakov et al. 2004a), we include the effects of accreted envelopes of light elements. We employ phenomenological density-dependent critical temperatures T_{cp}(\rho) and T_{cnt}(\rho) of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature T_{cns}(\rho) of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of Gusakov et al. (2004a) and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations., Comment: 19 pages, 4 figures; accepted for publication in MNRAS

We study the cooling of superfluid neutron stars whose cores consist of nucleon matter with the Akmal-Pandharipande-Ravenhall equation of state. This equation of state opens the powerful direct Urca process of neutrino emission in the interior of most massive neutron stars. Extending our previous studies (Gusakov et al. 2004a, Kaminker et al. 2005), we employ phenomenological density-dependent critical temperatures T_{cp}(\rho) of strong singlet-state proton pairing (with the maximum T_{cp}^{max} \sim 7e9 K in the outer stellar core) and T_{cnt}(\rho) of moderate triplet-state neutron pairing (with the maximum T_{cnt}^{max} \sim 6e8 K in the inner core). Choosing properly the position of T_{cnt}^{max} we can obtain a representative class of massive neutron stars whose cooling is intermediate between the cooling enhanced by the neutrino emission due to Cooper pairing of neutrons in the absence of the direct Urca process and the very fast cooling provided by the direct Urca process non-suppressed by superfluidity., Comment: 9 pages, 6 figures; accepted for publication in MNRAS

We have derived a set of equations to describe the thermal evolution of a neutron star which undergoes small-amplitude radial pulsations. We have taken into account, in the frame of the General Theory of Relativity, the pulsation damping due to the bulk and shear viscosity and the accompanying heating of the star. The neutrino emission of a pulsating non-superfluid star and its heating due to the bulk viscosity are calculated assuming that both processes are determined by the non-equilibrium modified Urca process. Analytical and numerical solutions to the set of equations of the stellar evolution are obtained for linear and strongly non-linear deviations from beta-equilibrium. It is shown that a pulsating star may be heated to very high temperatures, while the pulsations damp very slowly with time (a power law damping for 100-1000 years), as long as the damping is determined by the bulk viscosity. The contribution of the shear viscosity to the damping becomes important in a rather cool star with a low pulsation energy., Comment: 10 pages, 3 figures, an important reference to the paper by Finzi & Wolf (1968) is added; analytical consideration of the problem (Section 5) is essentially extended

We discuss pycnonuclear burning of highly exotic atomic nuclei in deep crusts of neutron stars, at densities up to 1e13 g/cc. As an application, we consider pycnonuclear burning of matter accreted on a neutron star in a soft X-ray transient (SXT, a compact binary containing a neutron star and a low-mass companion). The energy released in this burning, while the matter sinks into the stellar crust under the weight of newly accreted material, is sufficient to warm up the star and initiate neutrino emission in its core. The surface thermal radiation of the star in quiescent states becomes dependent of poorly known equation of state (EOS) of supranuclear matter in the stellar core, which gives a method to explore this EOS. Four qualitatively different model EOSs are tested against observations of SXTs. They imply different levels of the enhancement of neutrino emission in massive neutron stars by (1) the direct Urca process in nucleon/hyperon matter; (2) pion condensates; (3) kaon condensates; (4) Cooper pairing of neutrons in nucleon matter with the forbidden direct Urca process. A low level of the thermal quiescent emission of two SXTs, SAX J1808.4-3658 and Cen X-4, contradicts model (4). Observations of SXTs test the same physics of dense matter as observations of thermal radiation from cooling isolated neutron stars, but the data on SXTs are currently more conclusive., Comment: 4 pages, 3 figures, style file included. Proc. of ENAM-2004, Callaway Gardens, Sept. 2004; Europ. Phys. J. A, 2005, in press

We investigate the combined effect of neutron and proton superfluidities on the cooling of neutron stars whose cores consist of nucleons and electrons. We consider singlet-state pairing of protons and triplet-state pairing of neutrons in the cores of neutron stars. The critical superfluid temperatures T_c are assumed to depend on the density of matter. We study two types of neutron pairing with different components of the total angular momentum of Cooper pairs along the quantization axis (|m_J| =0 or 2). Our calculations are compared with observations of thermal emission from isolated neutron stars. We show that the observations can be interpreted by using two classes of superfluidity models: (1) strong proton superfluidity with a maximum critical temperature in the stellar core T_c^{max} > 4 \times 10^9 K and weak neutron superfluidity of any type (T_c^{max} < 2 \times 10^8 K); (2) strong neutron superfluidity (pairing with |m_J|=0) and weak proton superfluidity. The two types of models reflect an approximate symmetry with respect to an interchange of the critical temperatures of neutron and proton pairing., Comment: 20 pages, 8 figures

The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. Physical properties of neutron star matter important for cooling are reviewed such as composition, the equation of state, superfluidity of various baryon species, neutrino emission mechanisms. Theoretical results are compared with observations of thermal radiation from neutron stars. Current constraints on theoretical models of dense matter, derived from such a comparison, are formulated., Comment: 10 pages, 8 figures, to appear in Proceedings of the International Nuclear Physics Conference (INPS2004) (Goeteborg, Sweden, June 27 - July 2, 2004)

We show that the direct Urca process of neutrino emission is allowed in two possible phases of nonspherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3 orders of magnitude higher than the emissivity of the modified Urca process in the stellar core. The process may noticeably accelerate the cooling of low-mass neutron stars., Comment: 7 pages, 3 figures, submitted to A&A

We simulate cooling of superfluid neutron stars with nucleon cores where direct Urca process is forbidden. We adopt density dependent critical temperatures $T_{cp}(\rho)$ and $T_{cn}(\rho)$ of singlet-state proton and triplet-state neutron pairing in a stellar core and consider a strong proton pairing (with maximum $T_{cp}^{max} \ga 5 \times 10^9$ K) and a moderate neutron pairing ($T_{cn}^{max} \sim 6 \times 10^8$ K). When the internal stellar temperature $T$ falls below $T_{cn}^{max}$, the neutrino luminosity $L_{CP}$ due to Cooper pairing of neutrons behaves $\propto T^8$, just as that produced by modified Urca process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing neutrino emission acts like an enhanced cooling agent. By tuning the density dependence $T_{cn}(\rho)$ we can explain observations of cooling isolated neutron stars in the scenario in which direct Urca process or similar process in kaon/pion condensed or quark matter are absent., Comment: 9 pages, 3 figures, submitted to A&A

During the APPLES parallel campaign, the HST Advanced Camera for Surveys has resolved a distant stellar system, which appears to be an isolated dwarf galaxy. It is characterized by a circularly symmetric distribution of stars with an integrated magnitude m(F775W) = 20.13 +- 0.02, a central surface brightness of ~ 21.33 +- 0.18 mag/arcsec^2 and a half-light radius of ~ 1.8 arcsec. The ACS and VLT spectra show no evidence of ionized gas and appear dominated by a 3 Gyr old stellar population. The OB spectral type derived for two resolved stars in the grism data and the systemic radial velocity of ~ 670 km/s measured from the VLT data give a fiducial distance of ~ 9 +- 2 Mpc. These findings, with the support of the spatial morphology, would classify the system among the dwarf spheroidal (dSph) galaxies. Following the IAU rules, we have named this newly discovered galaxy APPLES 1. An intriguing peculiarity of APPLES 1 is that the properties (age and metallicity) of the stellar content so far detected are similar to those of dSph galaxies in the Local Group, where star formation is thought to be driven by galaxy interactions and mergers. Yet, APPLES 1 seems not to be associated with a major group or cluster of galaxies. Therefore, APPLES 1 could be the first example of a field dSph galaxy with self-sustained and regulated star formation and, therefore, would make an interesting test case for studies of the formation and evolution of unperturbed dSph galaxies., Comment: 46 pages, 13 figures, AJ accepted

We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer (of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic weight A > 100) on the top (Haensel & Zdunik 1990, 2003, astro-ph/0305220). The relation between the internal and surface stellar temperatures is obtained and fitted. The quiescent luminosity of the hottest (low-mass) and coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. The observations of SAX J1808.4-3658 in a quiescent state on March 24, 2001 (Campana et al. 2002, astro-ph/0206376) can be explained only if this SXT contains a massive neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal emission of SAX J1808.4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers and can serve for testing such models., Comment: 12 pages, 5 figures, 4 tables, LaTeX2e with aa.cls v.5.3 (included). Accepted by A&A

The cooling theory of isolated neutron stars is reviewed. The main cooling regulators are discussed, first of all, operation of direct Urca process (or similar processes in exotic phases of dense matter) and superfluidity in stellar interiors. The prospects to constrain gross parameters of supranuclear matter in neutron-star interiors by confronting cooling theory with observations of isolated neutron stars are outlined. A related problem of thermal states of transiently accreting neutron stars with deep crustal heating of accreted matter is discussed in application to soft X-ray transients., Comment: 10 pages, 3 figures, Proceedings of the 34th COSPAR Scientific Assembly (Adv. Sp. Res., accepted)

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The relation between the internal and local surface temperatures is calculated and fitted by an analytic function of the internal temperature, magnetic field strength, angle between the field lines and the normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron superfluidity in the stellar crust. Our results are important for interpretation of observations of isolated neutron stars hottest for their age, such as RX J0822-43 and PSR B1055-52., Comment: 15 pages, 12 figures, 2 tables. Corrected title only (v2)

We extend our modeling of cooling of superfluid neutron stars by including the production of muons in the core, in addition to neutrons, protons, and electrons. The results are confronted with observations of middle-aged isolated NSs. Muons have little effect on the hydrostatic structure of NSs, on the slow cooling of low-mass NSs (RX J0822--43 and PSR 1055--52 in our model) and on the rapid cooling of massive NSs. They affect, however, the moderately fast cooling of medium-mass NSs (1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga) and shift appreciably the mass range of these NSs to lower masses, which is important for correct interpretation of the observations. Moreover, the effects of muons can accurately be reproduced by a simple renormalization of NS models with no muons in the NS cores., Comment: 10 pages, 3 figures included, accepted for publication in Acta Physica Polonica

We consider cooling of neutron stars (NSs) with superfluid cores composed of neutrons, protons, and electrons (assuming singlet-state pairing of protons, and triplet-state pairing of neutrons). We mainly focus on (nonstandard) triplet-state pairing of neutrons with the $|m_J| = 2$ projection of the total angular momentum of Cooper pairs onto quantization axis. The specific feature of this pairing is that it leads to a power-law (nonexponential) reduction of the emissivity of the main neutrino processes by neutron superfluidity. For a wide range of neutron critical temperatures $T_{cn}$, the cooling of NSs with the $|m_J| = 2$ superfluidity is either the same as the cooling with the $m_J = 0$ superfluidity, considered in the majority of papers, or much faster. The cooling of NSs with density dependent critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ can be imitated by the cooling of the NSs with some effective critical temperatures $T_{cn}$ and $T_{cp}$ constant over NS cores. The hypothesis of strong neutron superfluidity with $|m_J| = 2$ is inconsistent with current observations of thermal emission from NSs, but the hypothesis of weak neutron superfluidity of any type does not contradict to observations., Comment: 10 pages, 6 figures

The upper limit of the effective surface temperature of the neutron star (NS) PSR J0205+6449 in the supernova remnant 3C 58 obtained recently by Slane et al. (2002) is analyzed using a modern theory of NS cooling (Kaminker et al. 2002). The observational limit can be explained by cooling of a superfluid NS with the core composed of neutrons, protons and electrons, where direct Urca process is forbidden. However, combined with the data on the surface temperatures of other isolated NSs, it gives evidence (emphasized by Slane et al.) that direct Urca process is open in the inner cores of massive NSs. This evidence turns out to be less stringent than the evidence provided by the well known observations of Vela and Geminga., Comment: 3 pages, 2 figures, AA Letters (submitted)

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is analyzed. The main cooling regulators are discussed: opening of direct Urca process in a NS central kernel; superfluidity of nucleons in NS interiors; surface layers of light (accreted) elements; strong surface magnetic fields. An emphasis is paid on the cooling scenario with strong $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, as well as strong $^1$S$_0$ pairing of neutrons in the NS crust. The theory predicts three types of isolated cooling middle-aged NSs with distinctly different properties: low-mass, slowly cooling NSs; medium-mass, moderately cooling NSs; massive, rapidly cooling NSs. The theory is compared with observations of isolated NSs -- pulsars and radio quiet NSs in supernova remnants. The constraints on physical properties of NSs which can be inferred from such a comparison are outlined., Comment: 14 pages, 13 figures, he_symp; to be published in Proceedings of the 270. Heraeus Seminar on Neutron Stars, Pulsars and Supernova Remnants, Bad Honnef, Jan. 21-25, 2002

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is simulated assuming $^1$S$_0$ pairing of neutrons in the NS crust, and also $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, and using realistic density profiles of the superfluid critical temperatures $T_{\rm c}(\rho)$. The theoretical cooling models of isolated middle-aged NSs can be divided into three main types. (I) {\it Low-mass}, {\it slowly cooling} NSs where the direct Urca process of neutrino emission is either forbidden or almost fully suppressed by the proton superfluidity. (II) {\it Medium-mass} NSs which show {\it moderate} cooling via the direct Urca process suppressed by the proton superfluidity. (III) {\it Massive} NSs which show {\it fast} cooling via the direct Urca process weakly suppressed by superfluidity. Confronting the theory with observations we treat RX J0822--43, PSR 1055--52 and RX J1856--3754 as slowly cooling NSs. To explain these sufficiently warm sources we need a density profile $T_{\rm c}(\rho)$ in the crust with a rather high and flat maximum and sharp wings. We treat 1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga as moderately cooling NSs. We can determine their masses for a given model of proton superfluidity, $T_{\rm cp}(\rho)$, and the equation of state in the NS core. No rapidly cooling NS has been observed so far., Comment: 12 pages, 10 figures, Astron. Astrophys., submitted

As shown recently by Kaminker et al. (2001), current observations of thermal emission of isolated middle-aged neutron stars (NSs) can be explained by cooling of NSs of different masses with the cores composed of neutrons, protons and electrons, assuming rather strong superfluidity (SF) of protons, weak triplet-state SF of neutrons and neglecting singlet-state SF of neutrons. We show that this explanation remains correct in the presence of singlet-state SF of neutrons in the NS crust and outermost core but under stringent constraints on the density profile of the SF critical temperature $T_{\rm cns}(\rho)$. In order to explain observations of (young and hot) RX J0822--43 and (old and warm) PSR 1055--52 and RX J1856--3754 as cooling not too massive NSs, the maximum $T_{\rm cns}^{\rm max}$ should be rather high ($\ga 5 \times 10^9$ K) and the decrease of $T_{\rm cns}(\rho)$ outside the maximum should be sharp. These results place important constraints on the models of nucleon SF in dense matter., Comment: 4 pages, 3 figures, AA Letters, accepted

Matter of subnuclear density in the inner crusts of neutron stars consists of neutron-rich atomic nuclei immersed in strongly degenerate relativistic gas of electrons and strongly nonideal liquid of neutrons. We give a brief account of thermodynamic and kinetic properties of this matter which are greatly affected by Coulomb and nuclear interactions (a companion paper on the outer envelopes is astro-ph/0012316) and show how they can be studied from observations of thermal radiation of young (age < 100 yr) neutron stars (for a discussion of the latter possibility, see astro-ph/0012306)., Comment: 4 pages, 1 figure, LaTeX2e using cpp2e.cls (included). Proc. PNP-10 Workshop, Greifswald, Germany, 4-9 Sept. 2000. Accepted for publication in Contrib. Plasma Phys. 41 (2001) no. 2-3

We review the main neutrino emission mechanisms in neutron star crusts and cores. Among them are the well-known reactions such as the electron-positron annihilation, plasmon decay, neutrino bremsstrahlung of electrons colliding with atomic nuclei in the crust, as well as the Urca processes and neutrino bremsstrahlung in nucleon-nucleon collisions in the core. We emphasize recent theoretical achievements, for instance, band structure effects in neutrino emission due to scattering of electrons in Coulomb crystals of atomic nuclei. We consider the standard composition of matter (neutrons, protons, electrons, muons, hyperons) in the core, and also the case of exotic constituents such as the pion or kaon condensates and quark matter. We discuss the reduction of the neutrino emissivities by nucleon superfluidity, as well as the specific neutrino emission produced by Cooper pairing of the superfluid particles. We also analyze the effects of strong magnetic fields on some reactions, such as the direct Urca process and the neutrino synchrotron emission of electrons. The results are presented in the form convenient for practical use. We illustrate the effects of various neutrino reactions on the cooling of neutron stars. In particular, the neutrino emission in the crust is critical in setting the initial thermal relaxation between the core and the crust. Finally, we discuss the prospects of exploring the properties of supernuclear matter by confronting cooling simulations with observations of the thermal radiation from isolated neutron stars., Comment: review, 165 pages, Physics Reports, 2001 in press