A numerical rotating neutron star solver is used to study the temporal evolution of accreting neutron stars using a multi-polytrope model for the nuclear equation of state named ACB5. The solver is based on a quadrupole expansion of the metric, but confirms the results of previous works, revealing the possibility of an abrupt transition of a neutron star from a purely hadronic branch to a third-family branch of stable hybrid stars, passing through an unstable intermediate branch. The accretion is described through a sequence of stationary rotating {stellar} configurations which lose angular momentum through magnetic dipole emission while, at the same time, gaining angular momentum through mass accretion. The model has several free parameters which are inferred from observations. The mass accretion scenario is studied in dependence on the effectiveness of angular momentum transfer which determines at which spin frequency the neutron star will become unstable against gravitational collapse to the corresponding hybrid star on the stable third-family branch. It is conceivable that the neutrino burst which accompanies the deconfinement transition may trigger a pulsar kick which results in the eccentric orbit. A consequence of the present model is the prediction of a correlation between the spin frequency of the millisecond pulsar in the eccentric orbit and its mass at birth.

In the first part of this paper, we investigate the possible existence of a structured hadron-quark mixed phase in the cores of neutron stars. This phase, referred to as the hadron-quark pasta phase, consists of spherical blob, rod, and slab rare phase geometries. Particular emphasis is given to modeling the size othis phase in rotating neutron stars. We use the relativistic mean-field theory to model hadronic matter and the non-local three-flavor Nambu-Jona-Lasinio model to describe quark matter. Based on these models, the hadron-quark pasta phase exists only in very massive neutron stars, whose rotational frequencies are less than around 300 Hz. All other stars are not dense enough to trigger quark deconfinement in their cores. Part two of the paper deals with the quark-hadron composition of hot (proto) neutron star matter. To this end we use a local three-flavor Polyakov-Nambu-Jona-Lasinio model which includes the 't Hooft (quark flavor mixing) term. It is found that this term leads to non-negligible changes in the particle composition of (proto) neutron stars made of hadron-quark matter.

In part one of this paper, we use a non-local extension of the 3-flavor Polyakov-Nambu-Jona-Lasinio model, which takes into account flavor-mixing, momentum dependent quark masses, and vector interactions among quarks, to investigate the possible existence of a spinodal region (determined by the vanishing of the speed of sound) in the QCD phase diagram and determine the temperature and chemical potential of the critical end point. In part two of the paper, we investigate the quark-hadron composition of baryonic matter at zero as well as non-zero temperature. This is of great topical interest for the analysis and interpretation of neutron star merger events such as GW170817. With this in mind, we determine the composition of proto-neutron star matter for entropies and lepton fractions that are typical of such matter. These compositions are used to delineate the evolution of proto-neutron stars to neutron stars in the baryon-mass versus gravitational-mass diagram. The hot stellar models turn out to contain significant fractions of hyperons and $\Delta$-isobars but no deconfined quarks. The latter, are found to exist only in cold neutron stars.

The recent direct observation of gravitational wave event $GW170817$ and its $GRB170817A$ signal has opened up a new window to study neutron stars and heralds a new era of Astronomy referred to as the Multimessenger Astronomy. Both gravitational and electromagnetic waves from a single astrophysical source have been detected for the first time. This combined detection offers an unprecedented opportunity to place constraints on the neutron star matter equation of state. The existence of a possible hadron-quark phase transition in the central regions of neutron stars is associated with the appearance of g-modes, which are extremely important as they could signal the presence of a pure quark matter core in the centers of neutron stars. Observations of g-modes with frequencies between 1 kHz and 1.5 kHz could be interpreted as evidence of a sharp hadron-quark phase transition in the cores of neutron stars. In this article, we shall review the description of the dense matter composing neutron stars, the determination of the equation of state of such matter, and the constraints imposed by astrophysical observations of these fascinating compact objects.

In this chapter we will introduce an effective equation of state (EoS) model based on polytropes that serves to study the so called "mass twins" scenario, where two compact stars have approximately the same mass but (significant for observation) quite different radii. Stellar mass twin configurations are obtained if a strong first-order phase transition occurs in the interior of a compact star. In the mass-radius diagram of compact stars, this will lead to a third branch of gravitationally stable stars with features that are very distinctive from those of white dwarfs and neutron stars. We discuss rotating hybrid star sequences in the slow rotation approximation and in full general relativity and draw conclusions for an upper limit on the maximum mass of nonrotating compact stars that has recently be deduced from the observation of the merger event GW170817.

The thermal evolution of hypernuclear compact stars is studied for stellar models constructed on the basis of covariant density functional theory in Hartree and Hartree-Fock approximation. Parametrizations of both types are consistent with the astrophysical mass constraints on compact stars and available hypernuclear data. We discuss the differences of these density functionals and highlight the effects they have on the composition and on the cooling of hypernuclear stars. It is shown that hypernuclear stars computed with density functional models that have a low symmetry energy slope, $L$, are fairly consistent with the cooling data of observed compact stars. The class of stellar models based on larger $L$ values gives rise to the direct Urca process at low densities, which leads to significantly faster cooling. We conjecture high-density pairing for protons and $\Lambda$'s in the $P$-wave channel and provide simple scaling arguments to obtain these gaps. As a consequence the most massive stellar models with masses $1.8 \le M/M_{\odot} \le2$ experience slower cooling by hyperonic dUrca processes which involve $\Lambda$'s and protons.

The $\Delta$-isobar degrees of freedom are included in the covariant density functional (CDF) theory to study the equation of state (EoS) and composition of dense matter in compact stars. In addition to $\Delta$'s we include the full octet of baryons, which allows us to study the interplay between the onset of delta isobars and hyperonic degrees of freedom. Using both the Hartree and Hartree-Fock approximation we find that $\Delta$'s appear already at densities slightly above the saturation density of nuclear matter for a wide range of the meson-$\Delta$ coupling constants. This delays the appearance of hyperons and significantly affects the gross properties of compact stars. Specifically, $\Delta$'s soften the EoS at low densities but stiffen it at high densities. This softening reduces the radius of a canonical $1.4 M_\odot$ star by up to 2~km for a reasonably attractive $\Delta$ potential in matter, while the stiffening results in larger maximum masses of compact stars. We conclude that the hypernuclear CDF parametrizations that satisfy the 2$M_\odot$ maximum mass constraint remain valid when $\Delta$ isobars are included, with the important consequence that the resulting stellar radii are shifted toward lower values, which is in agreement with the analysis of neutron star radii.

In this work we investigate the effect a crystalline quark-hadron mixed phase can have on the neutrino emissivity from the cores of neutron stars. To this end we use relativistic mean-field equations of state to model hadronic matter and a nonlocal extension of the three-flavor Nambu-Jona-Lasinio model for quark matter. Next we determine the extent of the quark-hadron mixed phase and its crystalline structure using the Glendenning construction, allowing for the formation of spherical blob, rod, and slab rare phase geometries. Finally we calculate the neutrino emissivity due to electron-lattice interactions utilizing the formalism developed for the analogous process in neutron star crusts. We find that the contribution to the neutrino emissivity due to the presence of a crystalline quark-hadron mixed phase is substantial compared to other mechanisms at fairly low temperatures ($\lesssim 10^9$ K) and quark fractions ($\lesssim 30\%$), and that contributions due to lattice vibrations are insignificant compared to static-lattice contributions. There are a number of open issues that need to be addressed in a future study on the neutrino emission rates caused by electron-quark blob bremsstrahlung. Chiefly among them are the role of collective oscillations of matter, electron band structures, and of gaps at the boundaries of the Brillouin zones on bremsstrahlung, as discussed in the summary section of this paper. We hope this paper will stimulate studies addressing these issues.

We study the thermal evolution of hypernuclear compact stars constructed from covariant density functional theory of hypernuclear matter and parameterizations which produce sequences of stars containing two-solar-mass objects. For the input in the simulations, we solve the Bardeen-Cooper-Schrieffer gap equations in the hyperonic sector and obtain the gaps in the spectra of $\Lambda$, $\Xi^0$ and $\Xi^-$ hyperons. For the models with masses $M/M_{\odot} \ge 1.5$ the neutrino cooling is dominated by hyperonic direct Urca processes in general. In the low-mass stars the $(\Lambda p)$ plus leptons channel is the dominant direct Urca process, whereas for more massive stars the purely hyperonic channels $(\Sigma^-\Lambda)$ and $(\Xi^-\Lambda)$ are dominant. Hyperonic pairing strongly suppresses the processes on $\Xi^-$s and to a lesser degree on $\Lambda$s. We find that intermediate-mass $1.5 \le M/M_{\odot} \le 1.8$ models have surface temperatures which lie within the range inferred from thermally emitting neutron stars, if the hyperonic pairing is taken into account. Most massive models with $M/M_{\odot} \simeq 2$ may cool very fast via the direct Urca process through the $(\Lambda p)$ channel because they develop inner cores where the $S$-wave pairing of $\Lambda$s and proton is absent.

In this work, we study the properties of magnetized white dwarfs taking into account possible instabilities due to electron capture and pycnonuclear fusion reactions in the cores of such objects. The structure of white dwarfs is obtained by solving the Einstein-Maxwell equations with a poloidal magnetic field in a fully general relativistic approach. The stellar interior is composed of a regular crystal lattice made of carbon ions immersed in a degenerate relativistic electron gas. The onsets of electron capture reactions and pycnonuclear reactions are determined with and without magnetic fields. We find that magnetized white dwarfs violate the standard Chandrasekhar mass limit significantly, even when electron capture and pycnonuclear fusion reactions are present in the stellar interior. We obtain a maximum white dwarf mass of around $2.14\,M_{\odot}$ for a central magnetic field of $\sim 3.85\times 10^{14}$~G, which indicates that magnetized white dwarfs may play a role for the interpretation of superluminous type Ia supernovae. Furthermore, we show that the critical density for pycnonuclear fusion reactions limits the central white dwarf density to $9.35\times 10^9$ g/cm$^3$. As a result, equatorial radii of white dwarfs cannot be smaller than $\sim 1100$~km. Another interesting feature concerns the relationship between the central stellar density and the strength of the magnetic field at the core of a magnetized white dwarf. For high magnetic fields, we find that the central density increases (stellar radius decrease) with magnetic field strength, which makes ultramagnetized white dwarfs more compact. The opposite is the case, however, if the central magnetic field is less than $\sim 10^{13}$~G. In the latter case, the central density decreases (stellar radius increases) with central magnetic field strengths.

We study the QCD phase diagram using a non-local SU(3) NJL model with vector interactions among quarks. We analyze several thermodynamic quantities such as entropy and specific heat, and study the influence of vector interactions on the thermodynamic properties of quark matter. Upon imposing electric charge neutrality and baryon number conservation on the field equations, we compute models for the equation of state of the inner cores of proto-neutron stars providing a non-local treatment of quark matter for astrophysics.

In this paper, we discuss the impact of rotation on the particle composition of rotating neutron stars (pulsars). Particular emphasis is put on the formation of quark matter during stellar spin-down, driven by continuous gravitational compression. Our study is based on modern models for the nuclear equation of state whose parameters are tightly constrained by nuclear data, neutron star masses, and the latest estimates of neutron star radii.

In this paper, we use a three flavor non-local Nambu--Jona-Lasinio (NJL) model, an~improved effective model of Quantum Chromodynamics (QCD) at low energies, to investigate the existence of deconfined quarks in the cores of neutron stars. Particular emphasis is put on the possible existence of quark matter in the cores of rotating neutron stars (pulsars). In contrast to non-rotating neutron stars, whose particle compositions do not change with time (are frozen in), the type and structure of the matter in the cores of rotating neutron stars depends on the spin frequencies of these stars, which opens up a possible new window on the nature of matter deep in the cores of neutron stars. Our study shows that, depending on mass and rotational frequency, up to around 8% of the mass of a massive neutron star may be in the mixed quark-hadron phase, if the phase transition is treated as a Gibbs transition. We also find that the gravitational mass at which quark deconfinement occurs in rotating neutron stars varies quadratically with spin frequency, which can be fitted by a simple formula.

We compute models for the equation of state (EoS) of the matter in the cores of hybrid stars. Hadronic matter is treated in the non-linear relativistic mean-field approximation, and quark matter is modeled by three-flavor local and non-local Nambu$-$Jona-Lasinio (NJL) models with repulsive vector interactions. The transition from hadronic to quark matter is constructed by considering either a soft phase transition (Gibbs construction) or a sharp phase transition (Maxwell construction). We find that high-mass neutron stars with masses up to $2.1-2.4 M_\odot$ may contain a mixed phase with hadrons and quarks in their cores, if global charge conservation is imposed via the Gibbs conditions. However, if the Maxwell conditions is considered, the appearance of a pure quark matter core either destabilizes the star immediately (commonly for non-local NJL models) or leads to a very short hybrid star branch in the mass-radius relation (generally for local NJL models).

The discovery of pulsars as heavy as 2 solar masses has led astrophysicists to rethink the core compositions of neutron stars, ruling out many models for the nuclear equations of state (EoS). We explore the hybrid stars that occur when hadronic matter is treated in a relativistic mean-field approximation and quark matter is modeled by three-flavor local and non-local Nambu Jona-Lasinio (NJL) models with repulsive vector interactions. The NJL models typically yield equations of state that feature a first order transition to quark matter. Assuming that the quark-hadron surface tension is high enough to disfavour mixed phases, and restricting to EoSes that allow stars to reach 2 solar masses, we find that the appearance of the quark matter core either destabilizes the star immediately (this is typical for non-local NJL models) or leads to a very short hybrid star branch in the mass-radius relation (this is typical for local NJL models). Using the Constant-Sound-Speed parametrization we can see that the reason for the near-absence of hybrid stars is that the transition pressure is fairly high and the transition is strongly first order.

Numerous theoretical studies using various equation of state models have shown that quark matter may exist at the extreme densities in the cores of high-mass neutron stars. It has also been shown that a phase transition from hadronic matter to quark matter would result in an extended mixed phase region that would segregate phases by net charge to minimize the total energy of the phase, leading to the formation of a crystalline lattice. The existence of quark matter in the core of a neutron star may have significant consequences for its thermal evolution, which for thousands of years is facilitated primarily by neutrino emission. In this work we investigate the effect a crystalline quark-hadron mixed phase can have on the neutrino emissivity from the core. To this end we calculate the equation of state using the relativistic mean-field approximation to model hadronic matter and a nonlocal extension of the three-flavor Nambu-Jona-Lasinio model for quark matter. Next we determine the extent of the quark-hadron mixed phase and its crystalline structure using the Glendenning construction, allowing for the formation of spherical blob, rod, and slab rare phase geometries. Finally we calculate the neutrino emissivity due to electron-lattice interactions utilizing the formalism developed for the analogous process in neutron star crusts. We find that the contribution to the neutrino emissivity due to the presence of a crystalline quark-hadron mixed phase is substantial compared to other mechanisms at fairly low temperatures ($\lesssim 10^9$ K) and quark fractions ($\lesssim 30\%$), and that contributions due to lattice vibrations are insignificant compared to static-lattice contributions.

We study the thermal evolution of neutron stars containing deconfined quark matter in their core. Such objects are generally referred to as quark-hybrid stars. The confined hadronic matter in their core is described in the framework of non-linear relativistic nuclear field theory. For the quark phase we use a non-local extension of the SU(3) Nambu Jona-Lasinio model with vector interactions. The Gibbs condition is used to model phase equilibrium between confined hadronic matter and deconfined quark matter. Our study indicates that high-mass neutron stars may contain between 35 and 40 % deconfined quark-hybrid matter in their cores. Neutron stars with canonical masses of around $1.4\, M_\odot$ would not contain deconfined quark matter. The central proton fractions of the stars are found to be high, enabling them to cool rapidly. Very good agreement with the temperature evolution established for the neutron star in Cassiopeia A (Cas A) is obtained for one of our models (based on the popular NL3 nuclear parametrization), if the protons in the core of our stellar models are strongly paired, the repulsion among the quarks is mildly repulsive, and the mass of Cas A has a canonical value of $1.4\, M_\odot$.

Non-rotating neutron stars are generally treated in theoretical studies as perfect spheres. Such a treatment, however, may not be correct if strong magnetic fields are present (such as for magnetars) and/or the pressure of the matter in the cores of neutron stars is non-isotropic (e.g., color superconducting). In this paper, we investigate the structure of non-spherical neutron stars in the framework of general relativity. Using a parameterized metric to model non-spherical mass distributions, we first derive a stellar structure equation for deformed neutron stars. Numerical investigations of this model equation show that the gravitational masses of deformed neutron stars depend rather strongly on the degree and type (oblate or prolate) of stellar deformation. In particular, we find that the mass of a neutron star increases with increasing oblateness but decreases with increasing prolateness. If this feature carries over to a full two-dimensional treatment of deformed neutron stars, this opens up the possibility that, depending on the type of stellar deformation, there may exist multiple maximum-mass neutron stars for one and for the same model for the nuclear equation of state.

This short review aims at giving a brief overview of the various states of matter that have been suggested to exist in the ultra-dense centers of neutron stars. Particular emphasis is put on the role of quark deconfinement in neutron stars and on the possible existence of compact stars made of absolutely stable strange quark matter (strange stars). Astrophysical phenomena, which distinguish neutron stars from quark stars, are discussed and the question of whether or not quark deconfinement may occur in neutron stars is investigated. Combined with observed astrophysical data, such studies are invaluable to delineate the complex structure of compressed baryonic matter and to put firm constraints on the largely unknown equation of state of such matter.

The recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

In this paper, we explore whether or not quark deconfinement may occur in high-mass neutron stars such as J1614-2230 (1.97 \pm 0.04 M_Sun) and J0348+0432 (2.01 \pm 0.04 M_Sun). Our study is based on a non-local extension of the SU(3) Nambu Jona-Lasinio (n3NJL) model with repulsive vector interactions among the quarks. This model goes beyond the frequently used local version of the Nambu Jona-Lasinio (NJL) model by accounting for several key features of QCD which are not part of the local model. Confined hadronic matter is treated in the framework of non-linear relativistic mean field theory. We find that both the local as well as the non-local NJL model predict the existence of extended regions of mixed quark-hadron (quark-hybrid) matter in high-mass neutron stars with masses of 2.1 to 2.4 M_Sun. Pure quark matter in the cores of neutron stars is obtained for certain parametrizations of the hadronic lagrangian and choices of the vector repulsion among quarks. The radii of high-mass neutron stars with quark-hybrid matter and/or pure quark matter cores in their centers are found to lie in the canonical range of 12 to 13 km.

In this work we consider the thermal evolution of rigidly rotating neutron stars. In order to perform such study we first calculate the structure of rotating objects, which is considerably more complicated than that of spherical objects. The structure of rotating neutron stars is obtained by solving Einstein's equation for a rotationally deformed fluid distributions. The numerical method used is based on the the KEH. The equation of state used for computing the neutron star structure and composition is a simple relativistic mean field model, with parameter set G300. With the structure of rotating neutron stars computed we calculate the thermal evolution of these objects. In order to do so, we re-derive the thermal evolution equations to account for the metric of a rotating object. The cooling of neutron stars with different frequencies is then calculated. We show that the cooling of the star strongly depends on the frequency of the object, with higher frequencies stars showing a substantial temperature difference between the equator and poles.

Depending on mass and rotational frequency, gravity compresses the matter in the core regions of neutron stars to densities that are several times higher than the density of ordinary atomic nuclei. At such huge densities atoms themselves collapse, and atomic nuclei are squeezed so tightly together that new particle states may appear and novel states of matter, foremost quark matter, may be created. This feature makes neutron stars superb astrophysical laboratories for a wide range of physical studies. And with observational data accumulating rapidly from both orbiting and ground based observatories spanning the spectrum from X-rays to radio wavelengths, there has never been a more exiting time than today to study neutron stars. The Hubble Space Telescope and X-ray satellites such as Chandra and XMM-Newton in particular have proven especially valuable. New astrophysical instruments such as the Five hundred meter Aperture Spherical Telescope (FAST), the square kilometer Array (skA), Fermi Gamma-ray Space Telescope (formerly GLAST), and possibly the International X-ray Observatory (now Advanced Telescope for High ENergy Astrophysics, ATHENA) promise the discovery of tens of thousands of new non-rotating and rotating neutron stars. The latter are referred to as pulsars. This paper provides a short overview of the impact of rotation on the structure and composition of neutron stars. Observational properties, which may help to shed light on the core composition of neutron stars--and, hence, the properties of ultra-dense matter--are discussed.

The unexpected temperature evolution of the neutron star in the Cassiopeia A supernova remnant (Cas A, for short) has renewed tremendous interest in the cooling mechanisms of neutron stars. In particular, the formation of superconducting protons and superfluid neutrons deep inside the cores of neutron stars have become focal points of the discussion. The purpose of this letter is to add a new aspect to this discussion, which focuses on the connection between proton superconductivity and the masses of neutron stars. Assuming (as is currently the case) that the temperature evolution of Cas A is largely controlled by superconducting protons, we study a series of phenomenological proton-pairing models to determine how deep into the stellar core superconducting protons actually penetrate. This allows us to establish a heretofore unknown relationship between the mass of the neutron star in Cas A and the penetration depth of the superconducting proton phase. This relationship can be used to either predict the depth of the superconducting proton phase, or, conversely, determine the mass of Cas A from a reliable calculation of the size of the proton superconducting phase in superdense neutron star matter. We emphasize that the strategy outlined in this paper can be applied to any other neutron star of similar age, whose temperature might be reliably monitored over a several years period. High-mass neutron stars, such as the recently discovered neutron stars J1614-2230 ($1.97 \pm 0.04\, \msun$) and J0348+0432 ($2.01 \pm 0.04 \, \msun$), appear particularly appealing as a significant fraction of the protons in their cores may be superconducting.

Driven by the loss of energy, isolated rotating neutron stars (pulsars) are gradually slowing down to lower frequencies, which increases the tremendous compression of the matter inside of them. This increase in compression changes both the global properties of rotating neutron stars as well as their hadronic core compositions. Both effects may register themselves observationally in the thermal evolution of such stars, as demonstrated in this Letter. The rotation-driven particle process which we consider here is the direct Urca (DU) process, which is known to become operative in neutron stars if the number of protons in the stellar core exceeds a critical limit of around 11% to 15%. We find that neutron stars spinning down from moderately high rotation rates of a few hundred Hertz may be creating just the right conditions where the DU process becomes operative, leading to an observable effect (enhanced cooling) in the temperature evolution of such neutron stars. As it turns out, the rotation-driven DU process could explain the unusual temperature evolution observed for the neutron star in Cas A, provided the mass of this neutron star lies in the range of 1.5 to 1.9 \msun and its rotational frequency at birth was between 40 (400 Hz) and 70% (800 Hz) of the Kepler (mass shedding) frequency, respectively.

Using a nonlocal extension of the SU(3) Nambu-Jona Lasinio model, which reproduces several of the key features of Quantum Chromodynamics, we show that mixed phases of deconfined quarks and confined hadrons (quark-hybrid matter) may exist in the cores of neutron stars as massive as around 2.1 M_Sun. The radii of these objects are found to be in the canonical range of $\sim 12-13$ km. According to our study, the transition to pure quark matter does not occur in stable neutron stars, but is shifted to neutron stars which are unstable against radial oscillations. The implications of our study for the recently discovered, massive neutron star PSR J1614-2230, whose gravitational mass is $1.97 \pm 0.04 M_Sun$, are that this neutron star may contain an extended region of quark-hybrid matter at it center, but no pure quark matter.

Already more that 40 years ago, it has been suggested that because of the enormous mass densities in the cores of neutron stars, the hadrons in the centers of neutron stars may undergo a phase transition to deconfined quark matter. In this picture, neutron stars could contain cores made of pure (up, down, strange) quark matter which are surrounded by a mixed phase of quarks and hadrons. More than that, because of the competition between the Coulomb and the surface energies associated with the positively charged regions of nuclear matter and negatively charged regions of quark matter, the mixed phase may develop geometrical structures, similarly to what is expected of the sub-nuclear liquid-gas phase transition. In this paper we restrict ourselves to considering the formation of rare phase blobs in the mixed quark-hadron phase. The influence of rare phase blobs on the thermal and transport properties of neutron star matter is investigated. The total specific heat, $c_V$, thermal conductivity, $\kappa$, and electron-blob Bremsstrahlung neutrino emissivities, $\epsilon_{\nu,\text{BR}}$, of quark-hybrid matter are computed and the results are compared with the associated thermal and transport properties of standard neutron star matter. Our results show that the contribution of rare phase blobs to the specific heat is negligibly small. This is different for the neutrino emissivity from electron-blob Bremsstrahlung scattering, which turns out to be of the same order of magnitude as the total contributions from other Bremsstrahlung processes for temperatures below about $10^8$ K.

This paper gives an brief overview of the structure of hypothetical strange quarks stars (quark stars, for short), which are made of absolutely stable 3-flavor strange quark matter. Such objects can be either bare or enveloped in thin nuclear crusts, which consist of heavy ions immersed in an electron gas. In contrast to neutron stars, the structure of quark stars is determined by two (rather than one) parameters, the central star density and the density at the base of the crust. If bare, quark stars possess ultra-high electric fields on the order of 10^{18} to 10^{19} V/cm. These features render the properties of quark stars more multifaceted than those of neutron stars and may allow one to observationally distinguish quark stars from neutron stars.

There are many factors that contribute to the breaking of the spherical symmetry of a neutron star. Most notably is rotation, magnetic fields, and/or accretion of matter from companion stars. All these phenomena influence the macroscopic structures of neutron stars, but also impact their microscopic compositions. The purpose of this paper is to investigate the cooling of rotationally deformed, two-dimensional (2D) neutron stars in the framework of general relativity theory, with the ultimate goal of better understand the impact of 2D effects on the thermal evolution of such objects. The equations that govern the thermal evolution of rotating neutron stars are presented in this paper. The cooling of neutron stars with different frequencies is computed self-consistently by combining a fully general relativistic 2D rotation code with a general relativistic 2D cooling code. We show that rotation can significantly influence the thermal evolution of rotating neutron stars. Among the major new aspects are the appearances of hot spots on the poles, and an increase of the thermal coupling times between the core and the crust of rotating neutron stars. We show that this increase is independent of the microscopic properties of the stellar core, but depends only on the frequency of the star.

If quark stars exist, they may be enveloped in thin electron layers (electron seas), which uniformly surround the entire star. These layers will be affected by the magnetic fields of quark stars in such a way that the electron seas would transmit hydromagnetic cyclotron waves, as studied in this paper. Particular attention is devoted to vortex hydrodynamical oscillations of the electron sea. The frequency spectrum of these oscillations is derived in analytic form. If the thermal X-ray spectra of quark stars are modulated by vortex hydrodynamical vibrations, the thermal spectra of compact stars, foremost central compact objects (CCOs) and X-ray dim isolated neutron stars (XDINSs), could be used to verify the existence of these vibrational modes observationally. The central compact object 1E 1207.4-5209 appears particularly interesting in this context, since its absorption features at 0.7 keV and 1.4 keV can be comfortably explained in the framework of the hydro-cyclotron oscillation model.

We discuss constraints for the equation of state of hybrid star matter which can be obtained from heavy-ion collisions at FAIR and NICA. Particular emphasis is on the planned NICA facility at JINR Dubna which shall provide fixed-target and collider experiments just in the relevant energy ranges.

We demonstrate that the high-quality cooling data observed for the young neutron star in the supernova remnant Cassiopeia A over the past 10 years--as well as all other reliably known temperature data of neutron stars--can be comfortably explained within the "nuclear medium cooling" scenario. The cooling rates of this scenario account for medium-modified one-pion exchange in dense matter and polarization effects in the pair-breaking formations of superfluid neutrons and protons. Crucial for the successful description of the observed data is a substantial reduction of the thermal conductivity, resulting from a suppression of both the electron and nucleon contributions to it by medium effects. We also find that possibly in as little as about ten years of continued observation, the data may tell whether or not fast cooling processes are active in this neutron star.

The structure of a magnetized quark-hybrid stars (QHS) is modeled using a standard relativistic mean-field equation of state (EoS) for the description of hadronic matter. For quark matter we consider a bag model EoS which is modified perturbatively to account for the presence of a uniform magnetic field. The mass-radius (M-R) relationship, gravitational redshift and rotational Kepler periods of such stars are compared with those of standard neutron stars (NS).

This paper provides an overview of the possible role of Quantum Chromo Dynamics (QDC) for neutron stars and strange stars. The fundamental degrees of freedom of QCD are quarks, which may exist as unconfined (color superconducting) particles in the cores of neutron stars. There is also the theoretical possibility that a significantly large number of up, down, and strange quarks may settle down in a new state of matter known as strange quark matter, which, by hypothesis, could be more stable than atomic nuclei. In the latter case new classes of self-bound, color superconducting objects, ranging from strange quark nuggets to strange quark stars, should exist. The properties of such objects will be reviewed along with the possible existence of deconfined quarks in neutron stars. Implications for observational astrophysics are pointed out.

In this paper we investigate the electrodynamic surface properties of bare strange quark stars. The surfaces of such objects are characterized by the formation of ultra-high electric surface fields which might be as high as $\sim 10^{19}$ V/cm. These fields result from the formation of electric dipole layers at the stellar surfaces. We calculate the increase in gravitational mass associated with the energy stored in the electric dipole field, which turns out to be only significant if the star possesses a sufficiently strong {\em net} electric charge distribution. In part two of the paper, we explore the intriguing possibility of what happens when the electron layer (sphere) rotates with respect to the stellar strange matter body. We find that in this event magnetic fields can be generated which, for moderate effective rotational frequencies between the electron layer and the stellar body, agree with the magnetic fields inferred for several Central Compact Objects (CCOs). These objects could thus be comfortably interpreted as strange stars whose electron atmospheres rotate at frequencies that are moderately different ($\sim 10$ Hz) from the rotational frequencies of the strange star itself.

Compact stars made of quark matter rather than confined hadronic matter, are expected to form a color superconductor. This superconductor ought to be threaded with rotational vortex lines within which the star's interior magnetic field is confined. The vortices (and thus magnetic flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at the surface of the star and the prolific production of thermal energy. In this Letter, we show that this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star enters the color superconducting phase, which can be used to estimate the density at which quarks deconfine. We find this density to be five times that of nuclear saturation.

Compact stars made of quark matter rather than confined hadronic matter, are expected to form a color superconductor. This superconductor ought to be threaded with rotational vortex lines, within which the star's interior magnetic field is at least partially confined. The vortices (and thus magnetic flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at the surface of the star and the prolific production of thermal energy. In this paper, we show that this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star enters the color superconducting phase, which can be used to estimate the density at which quarks deconfine. From observations, we find this density to be of the order of five times that of nuclear saturation.

New astrophysical instruments such as skA (square kilometer Array) and IXO (formerly Constellation X) promise the discovery of tens of thousands of new isolated rotating neutron stars (pulsars), neutron stars in low-mass X-ray binaries (LMXBs), anomalous X-ray pulsars (AXPs), and soft gamma repeaters (SGRs). Many of these neutron stars will experience dramatic density changes over their active lifetimes, driven by either stellar spin-up or spin-down, which may trigger phase transitions in their dense baryonic cores. More than that, accretion of matter onto neutron stars in LMXBs is believed to cause pycno-nuclear fusion reactions in the inner crusts of neutron stars. The associated reaction rates may be drastically altered if strange quark matter would be absolutely stable. This paper outlines the investigative steps that need to be performed in order to explore the thermal response of neutron stars to rotationally-driven phase transitions in their cores as well as to nuclear burning scenarios in their crusts. Such research complements the exploration of the phase diagram of dense baryonic matter through particle collider experiments, as performed at RHIC in the USA and as planned at the future Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany.

The possible existence of compact stars made of absolutely stable strange quark matter--referred to as strange stars--was pointed out by E. Witten almost a quarter of a century ago. One of the most amazing features of such objects concerns the possible existence of ultra-strong electric fields on their surfaces, which, for ordinary strange matter, is around $10^{18}$ V/cm. If strange matter forms a color superconductor, as expected for such matter, the strength of the electric field may increase to values that exceed $10^{19}$ V/cm. The energy density associated with such huge electric fields is on the same order of magnitude as the energy density of strange matter itself, which, as shown in this paper, alters the masses and radii of strange quark stars at the 15% and 5% level, respectively. Such mass increases facilitate the interpretation of massive compact stars, with masses of around $2 M_\odot$, as strange quark stars.

This paper presents an investigation into the pycnonuclear reaction rates in dense crustal matter of neutron stars contaminated with strange quark matter nuggets. The presence of such nuggets in the crustal matter of neutron stars would be a natural consequence if Witten's strange quark matter hypothesis is correct. The methodology presented in this paper is a recreation of a recent representation of nuclear force interactions embedded within pycnonuclear reaction processes. The study then extends the methodology to incorporate distinctive theoretical characteristics of strange quark matter nuggets, like their low charge-per-baryon ratio, and then assesses their effects on the pycnonuclear reaction rates. Particular emphasis is put on the impact of color superconductivity on the reaction rates. Depending on whether or not quark nuggets are in this novel state of matter, their electric charge properties vary drastically which turns out to have a dramatic effect on the pycnonuclear reaction rates. Future nuclear fusion network calculations may thus have the potential to shed light on the existence of strange quark matter nuggets and on whether or not they are in a color superconducting state, as suggested by QCD.

The thermal evolution of neutron stars is coupled to their spin down and the resulting changes in structure and chemical composition. This coupling correlates stellar surface temperatures with rotational state as well as time. We report an extensive investigation of the coupling between spin down and cooling for hybrid stars which undergo a phase transition to deconfined quark matter at the high densities present in stars at low rotation frequencies. The thermal balance of neutron stars is re-analyzed to incorporate phase transitions and the related latent heat self-consistently, and numerical calculations are undertaken to simultaneously evolve the stellar structure and temperature distribution. We find that the changes in stellar structure and chemical composition with the introduction of a pure quark matter phase in the core delay the cooling and produce a period of increasing surface temperature for strongly superfluid stars of strong and intermediate magnetic field strength. The latent heat of deconfinement is found to reinforce this signature if quark matter is superfluid and it can dominate the thermal balance during the formation of a pure quark matter core. At other times it is less important and does not significantly change the thermal evolution.

Depending on the density reached in the cores of neutron stars, such objects may contain stable phases of novel matter found nowhere else in the Universe. This article gives a brief overview of these phases of matter and discusses astrophysical constraints on the high-density equation of state associated with ultra-dense nuclear matter.

The three-flavor crystalline color-superconducting (CCS) phase of quantum chromodynamics (QCD) is a candidate phase for the ground state of cold matter at moderate densities above the density of the deconfinement phase transition. Apart from being a superfluid, the CCS phase has properties of a solid, such as a lattice structure and a shear modulus, and hence the ability to sustain multipolar deformations in gravitational equilibrium. We construct equilibrium configurations of hybrid stars composed of nuclear matter at low, and CCS quark matter at high, densities. Phase equilibrium between these phases is possible only for rather stiff equations of state of nuclear matter and large couplings in the effective Nambu--Jona-Lasinio Lagrangian describing the CCS state. We identify a new branch of stable CCS hybrid stars within a broad range of central densities which, depending on the details of the equations of state, either bifurcate from the nuclear sequence of stars when the central density exceeds that of the deconfinement phase transition or form a new family of configurations separated from the purely nuclear sequence by an instability region. The maximum masses of our non-rotating hybrid configurations are consistent with the presently available astronomical bounds. The sequences of hybrid configurations that rotate near the mass-shedding limit are found to be more compact and thus support substantially larger spins than their same mass nuclear counterparts.

Neutron stars contain matter in one of the densest forms found in the Universe. This feature, together with the unprecedented progress in observational astrophysics, makes such stars superb astrophysical laboratories for a broad range of exciting physical studies. This paper gives an overview of the phases of dense matter predicted to make their appearance in the cores of neutron stars. Particular emphasis is put on the role of strangeness. Net strangeness is carried by hyperons, K-mesons, H-dibaryons, and strange quark matter, and may leave its mark in the masses, radii, moment of inertia, dragging of local inertial frames, cooling behavior, surface composition, and the spin evolution of neutron stars. These observables play a key role for the exploration of the phase diagram of dense nuclear matter at high baryon number density but low temperature, which is not accessible to relativistic heavy ion collision experiments.

Because of the tremendous densities that exist in the cores of neutron stars, a significant fraction of the matter in the cores of such stars is likely to exist in the form of hyperons. Depending on spin frequency, the hyperon content changes dramatically in rotating neutron stars, as discussed in this paper.

There has been much recent progress in our understanding of quark matter, culminating in the discovery that if such matter exists in the cores of neutron stars it ought to be in a color superconducting state. This paper explores the impact of superconducting quark matter on the properties (e.g., masses, radii, surface gravity, photon emission) of compact stars.

A forefront area of research concerns the exploration of the properties of hadronic matter under extreme conditions of temperature and density, and the determination of the equation of state--the relation between pressure, temperature and density--of such matter. Experimentally, relativistic heavy-ion collision experiments enable physicists to cast a brief glance at hot and ultra-dense matter for times as little as about $10^{-22}$ seconds. Complementary to this, the matter that exists in the cores of neutron stars, observed as radio pulsars, X-ray pulsars, and magnetars, is at low temperatures but compressed permanently to ultra-high densities that may be more than an order of magnitude higher than the density of atomic nuclei. This makes pulsars superb astrophysical laboratories for medium and high-energy nuclear physics, as discussed in this paper.

With central densities way above the density of atomic nuclei, neutron stars contain matter in one of the densest forms found in the universe. Depending of the density reached in the cores of neutron stars, they may contain stable phases of exotic matter found nowhere else in space. This article gives a brief overview of the phases of ultra-dense matter predicted to exist deep inside neutron stars and discusses the equation of state associated with such matter.

It is generally agreed on that the tremendous densities reached in the centers of neutron stars provide a high-pressure environment in which several intriguing particles processes may compete with each other. These range from the generation of hyperons to quark deconfinement to the formation of kaon condensates and H-matter. There are theoretical suggestions of even more exotic processes inside neutron stars, such as the formation of absolutely stable strange quark matter. In the latter event, neutron stars would be largely composed of strange quark matter possibly enveloped in a thin nuclear crust. This paper gives a brief overview of these striking physical possibilities with an emphasis on the role played by strangeness in neutron star matter, which constitutes compressed baryonic matter at ultra-high baryon number density but low temperature which is no accessible to relativistic heavy ion collision experiments.

We summarize masses and radii for a number of white dwarfs as deduced from a combination of proper motion studies, Hipparcos parallax distances, effective temperatures, and binary or spectroscopic masses. A puzzling feature of these data is that some stars appear to have radii which are significantly smaller than that expected for a standard electron-degenerate white-dwarf equations of state. We construct a projection of white-dwarf radii for fixed effective mass and conclude that there is at least marginal evidence for bimodality in the radius distribution forwhite dwarfs. We argue that if such compact white dwarfs exist it is unlikely that they contain an iron core. We propose an alternative of strange-quark matter within the white-dwarf core. We also discuss the impact of the so-called color-flavor locked (CFL) state in strange-matter core associated with color superconductivity. We show that the data exhibit several features consistent with the expected mass-radius relation of strange dwarfs. We identify eight nearby white dwarfs which are possible candidates for strange matter cores and suggest observational tests of this hypothesis.

A new scheme for testing nuclear matter equations of state (EsoS) at high densities using constraints from neutron star phenomenology and a flow data analysis of heavy-ion collisions is suggested. An acceptable EoS shall not allow the direct Urca process to occur in neutron stars with masses below $1.5~M_{\odot}$, and also shall not contradict flow and kaon production data of heavy-ion collisions. Compact star constraints include the mass measurements of 2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807, of 2.0 +/- 0.1 M_sun from the innermost stable circular orbit for 4U 1636-536, the baryon mass - gravitational mass relationships from Pulsar B in J0737-3039 and the mass-radius relationships from quasiperiodic brightness oscillations in 4U 0614+09 and from the thermal emission of RX J1856-3754. This scheme is applied to a set of relativistic EsoS constrained otherwise from nuclear matter saturation properties with the result that no EoS can satisfy all constraints simultaneously, but those with density-dependent masses and coupling constants appear most promising.

Astrophysicists distinguish between three different types of compact stars. These are white dwarfs, neutron stars, and black holes. The former contain matter in one of the densest forms found in the Universe. This feature, together with the unprecedented progress in observational astronomy, makes such stars superb astrophysical laboratories for a broad range of exciting physical studies. This article studies the role of strangeness for compact star phenomenology. Strangeness is carried by hyperons, mesons, H-dibaryons, and strange quark matter, and may leave its mark in the masses, radii, cooling behavior, surface composition and the spin evolution of compact stars.

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.

Astrophysicists distinguish between three different types of compact stars. These are white dwarfs, neutron stars, and black holes. The former contain matter in one of the densest forms found in the Universe which, together with the unprecedented progress in observational astronomy, make such stars superb astrophysical laboratories for a broad range of most striking physical phenomena. These range from nuclear processes on the stellar surface to processes in electron degenerate matter at subnuclear densities to boson condensates and the existence of new states of baryonic matter--like color superconducting quark matter--at supernuclear densities. More than that, according to the strange matter hypothesis strange quark matter could be more stable than nuclear matter, in which case neutron stars should be largely composed of pure quark matter possibly enveloped in thin nuclear crusts. Another remarkable implication of the hypothesis is the possible existence of a new class of white dwarfs. This article aims at giving an overview of all these striking physical possibilities, with an emphasis on the astrophysical phenomenology of strange quark matter. Possible observational signatures associated with the theoretically proposed states of matter inside compact stars are discussed as well. They will provide most valuable information about the phase diagram of superdense nuclear matter at high baryon number density but low temperature, which is not accessible to relativistic heavy ion collision experiments.

We construct self-consistent equilibrium sequences of general relativistic, rotating neutron star models. Special emphasis in put on the determination of the maximum rotation frequency of such objects. Recently proposed models for the equation of state of neutron star matter are employed, which are derived by describing the hadronic phase within the many-body Brueckner--Bethe--Goldstone formalism, and the quark matter phase within the MIT bag model using a density dependent bag constant. We find that the rotational frequencies of neutron stars with deconfined quark phases in their cores rival those of absolutely stable, self-bound strange quark matter stars. This finding is of central importance for the interpretation of extremely rapidly rotating pulsars, which are the targets of present pulsar surveys.

There has never been a more exciting time in the overlapping areas of nuclear physics, particle physics and relativistic astrophysics than today. Orbiting observatories such as the Hubble Space Telescope, Rossi X-ray Timing Explorer (RXTE), Chandra X-ray satellite, and the X-ray Multi Mirror Mission (XMM) have extended our vision tremendously, allowing us to see vistas with an unprecedented clarity and angular resolution that previously were only imagined, enabling astrophysicists for the first time ever to perform detailed studies of large samples of galactic and extragalactic objects. On the Earth, radio telescopes (e.g., Arecibo, Green Bank, Parkes, VLA) and instruments using adaptive optics and other revolutionary techniques have exceeded previous expectations of what can be accomplished from the ground. The gravitational wave detectors LIGO, LISA VIRGO, and Geo-600 are opening up a window for the detection of gravitational waves emitted from compact stellar objects such as neutron stars and black holes. Together with new experimental forefront facilities like ISAC, ORLaND and RIA, these detectors provide direct, quantitative physical insight into nucleosynthesis, supernova dynamics, accreting compact objects, cosmic-ray acceleration, and pair-production in high energy sources which reinforce the urgent need for a strong and continuous feedback from nuclear and particle theory and theoretical astrophysics. In my lectures, I shall concentrate on three selected topics, which range from the behavior of superdense stellar matter, to general relativistic stellar models, to strange quark stars and possible signals of quark matter in neutron stars.

We explore the impact of in-medium modification of the properties of vector mesons on the nuclear equation of state and neutron star properties. It is found that in-medium modifications stiffen the nuclear equation of state considerably. If this feature has its correspondence in the full treatment of dense hadronic matter, then very little room is left for the existence of exotic phases like quark matter or boson condensates in the centers of neutron stars of canonical mass.

This paper discusses several most intruigung astrophysical implications connected with the possible absolute stability of strange quark matter.This is followed by a discussion of two astrophysical signals that may point at the existence of quark matter in both isolated neutron stars as well as in neutron stars in low-mass x-ray binaries (LMXBs).

We comment on the paper by Tanusri Ghosh and Somenath Chakrabarty, Phys. Rev. D63, 043006 (2001). In that paper it was argued that a first-order transition to quark matter in the core of a magnetar is absolutely forbidden if its magnetic field strength exceeds $10^{15}$ G. However, we show in this comment that if the quark anomalous magnetic moment and the population of higher Landau levels is taken into account, it may still be possible for a first-order phase transition from nuclear matter to quark matter to occur in the core of a magnetar. These effects may also obscure the question of whether beta equilibrium favors strange quark matter in the core of a magnetar.

A neutron star in binary orbit with a low-mass non-degenerate companion becomes a source of x-rays with millisecond variability when mass accretion spins it up. Centrifugally driven changes in density profile may initiate a phase transition in a growing region of the core parallel to what may take place in an isolated millisecond pulsar, but in reverse. Such a star will spend a longer time in the spin frequency range over which the transition occurs than elsewhere because the change of phase, paced by the spinup rate, is accompanied by a growth in the moment of inertia. The population of accreters will exhibit a clustering in the critical frequency range. A phase change triggered by changing spin and the accompanying adjustment of moment of inertia has its analogue in rotating nuclei.

Theoretically, the phase transition between the confined and deconfined phases of quarks can have a remarkable effect on the spin properties of millisecond pulsars and on the spin distribution of the population of x-ray neutron stars in low-mass binaries. In the latter class of stars, the effect has already been observed---a strong clustering in the population in a narrow band of spins. The observed clustering cannot presently be uniquely assigned to the phase transition as cause. However, there is another possible signal---not so far observed---in millisecond pulsars that we also discuss. It would have the same origin, and its discovery would tend to confirm the interpretation in terms of a phase transition in the stellar core and that the high density phase exists in the higer mass canonical neutron stars.

We model the evolution of canonical pulsars from the death line to millisecond pulsars through the X-ray neutron star stage of accretion from a low-mass companion. We trace this evolution in magnetic field strength starting at $B=10^{12}$ to $10^8$ G and in period of about 1 s to milliseconds. Important factors are accretion rate and the decay rate of the magnetic field. A broad swathe is traced in the $B-P$ plane according to the value of these factors, which represent different conditions of the binary pair. An important ingredient is the dependence of the stellar moment of inertia on rotation frequency (time).

We study the spin evolution of X-ray neutron stars in binary systems, which are being spun up by mass transfer from accretion disks. Our investigation reveals that a quark phase transition resulting from the changing central density induced by the changing spin, can lead to a pronounced peak in the frequency distribution of X-ray neutron stars. This finding provides one of several possible explanations available in the literature, or at least a contributor to part of the observed anomalous frequency distribution of neutron stars in low-mass X-ray binaries (LMXBs), which lie in a narrow band centered at about 300 Hz, as found by the Rossi Explorer (RXTE).

It is generally agreed on that the tremendous densities reached in the centers of neutron stars provide a high-pressure environment in which numerous novel particles processes are likely to compete with each other. These processes range from the generation of hyperons to quark deconfinement to the formation of kaon condensates and H-matter. There are theoretical suggestions of even more exotic processes inside neutron stars, such as the formation of absolutely stable strange quark matter, a configuration of matter even more stable than the most stable atomic nucleus, iron. In the latter event, neutron stars would be largely composed of pure quark matter, eventually enveloped in a thin nuclear crust. No matter which physical processes are actually realized inside neutron stars, each one leads to fingerprints, some more pronounced than others though, in the observable stellar quantities. This feature combined with the unprecedented progress in observational astronomy, which allows us to see vistas with remarkable clarity that previously were only imagined, renders neutron stars to nearly ideal probes for a wide range of physical studies, including the role of strangeness in dense matter.

Gravity compresses the matter in the cores of neutron stars to densities which are significantly higher than the density of ordinary atomic nuclei, thus providing a high-pressure environment in which numerous particle processes - from the generation of new baryonic particles to quark deconfinement to the formation of Boson condensates and H-matter - may compete with each other. There are theoretical suggestions of even more `exotic' processes inside pulsars, such as the formation of absolutely stable strange quark matter, a configuration of matter even more stable than the most stable atomic nucleus, iron. In the latter event, neutron stars would be largely composed of pure quark matter, eventually enveloped in nuclear crust matter. No matter which physical processes are actually realized inside neutron stars, each one leads to fingerprints, some more pronounced than others though, in the observable stellar quantities. This feature combined with the tremendous recent progress in observational radio and X-ray astronomy, renders neutron stars to nearly ideal probes for a wide range of dense matter studies, complementing the quest of the behavior of superdense matter in terrestrial collider experiments.

The impact of various competing heating processes on the thermal evolution of neutron stars is investigated. We show that internal heating leads to significantly enhanced surface temperatures for pulsars of middle and old age. The heating due to thermal creep of pinned vortices and due to outward motion of proton vortices in the interior of the star leads to a better agreement with the observed data in the case of enhanced cooling. The strong pinning models are ruled out by a comparison with the cooling data on the old pulsars. For millisecond pulsars, the heating due to thermal creep of pinned vortices and chemical heating of the core have the largest impact on the surface temperatures. The angular dependence of the heating rates require two dimensional cooling simulations in general. Such a simulation is performed for a selected case in order to check the applicability of one-dimensional codes used in the past.

It is shown that the modern equations of state for neutron star matter based on microscopic calculations of symmetric and asymmetric nuclear matter are compatible with the lower bound on the maximum neutron-star mass for a certain range of hyperon coupling constants, which are constrained by the binding energies of hyperons in symmetric nuclear matter. The hyperons are included by means of the relativistic Hartree-- or Hartree--Fock approximation. The obtained couplings are also in satisfactory agreement with hypernuclei data in the relativistic Hartree scheme. Within the relativistic Hartree--Fock approximation hypernuclei have not been investigated so far.

The properties of warm symmetric and asymmetric nuclear matter are investigated in the frame of the Thomas-Fermi approximation using a recent modern parametrization of the effective nucleon-nucleon interaction of Myers and Swiatecki. Special attention is paid to the liquid-gas phase transition, which is of special interest in modern nuclear physics. We have determined the critical temperature, critical density and the so-called flash temperature. Furthermore the equation of state for cold neutron star matter was calculated.

The possibility of a neutron m=2-superfluid in the interior of neutron stars is investigated. This pairing state is energetically favoured in strong magnetic fields ($H\sim 10^{16}-10^{17}$ G). Because of the node in the angular-dependent energy gap along the field direction the neutrino emissivity is only suppressed polynomially as function in $T/T_{c}$ instead of exponentially, as it is obtained for a nodeless pairing state. The effect of this pairing state on the thermal evolution of neutron stars is studied, and its outcome is compared with the evolution of ``normal'', i.e. nodeless, superfluid and non-superfluid neutron stars, and also with observations. We find that particularly the predicted surface temperatures of the enhanced cooling scenario considerably change and come into agreement with temperatures deduced from observational data within the hydrogen atmosphere model. Furthermore the surface temperature depends on the magnetic field strength as an additional parameter aside from the neutron star mass. The latter is however only operative in the case of the intermediate cooling scenario.

The properties of hot matter are studied in the frame of the relativistic Brueckner-Hartree-Fock theory. The equations are solved self-consistently in the full Dirac space. For the interaction we used the potentials given by Brockmann and Machleidt. The obtained critical temperatures are smaller than in most of the nonrelativistic investigations. We also calculated the thermodynamic properties of hot matter in the relativistic Hartree--Fock approximation, where the force parameters were adjusted to the outcome of the relativistic Brueckner--Hartree--Fock calculations at zero temperature. Here, one obtains higher critical temperatures, which are comparable with other relativistic calculations in the Hartree scheme.

We study the properties of static and rotating neutron stars adopting non-relativistic equations of state (EOS) for asymmetric nuclear matter based on the Brueckner-Hartree-Fock (BHF) scheme. The BHF calculation, with the continuous choice for the single particle potential, appears to be very close to the full EOS, which includes the three-hole line contribution calculated by solving the Bethe-Fadeev equations within the gap choice for the single particle potential. Three-body forces are included in order to reproduce the correct saturation point for nuclear matter. A comparison with fully relativistic many-body calculations of nuclear matter EOS is made.

We study the properties of neutron stars adopting relativistic equations of state of neutron star matter, calculated in the framework of the relativistic Brueckner-Hartree-Fock approximation for electrically charge neutral neutron star matter in beta-equilibrium. For higher densities more baryons (hyperons etc.) are included by means of the relativistic Hartree- or Hartree-Fock approximation. The special features of the different approximations and compositions are discussed in detail. Besides standard neutron star properties special emphasis is put on the limiting periods of neutron stars, for which the Kepler criterion and gravitation-reaction instabilities are considered. Furthermore the cooling behaviour of neutron stars is investigated, too. For comparison we give also the outcome for some nonrelativistic equations of state.

The general statement that strange stars cool more rapidly than neutron stars is investigated in greater detail. It is found that the direct Urca process could be forbidden not only in neutron stars but also in strange stars. If so, strange stars would be slowly cooling and their surface temperatures would be more or less indistinguishable from those of slowly cooling neutron stars. The case of enhanced cooling is reinvestigated as well. It is found that strange stars cool significantly more rapidly than neutron stars within the first $\sim 30$ years after birth. This feature could become particularly interesting if continued observation of SN 1987A would reveal the temperature of the possibly existing pulsar at its centre.

The conversion of nuclear matter to quark matter in the core of a rotating neutron star alters its moment of inertia. Hence the epoch over which conversion takes place will be signaled in the spin-down "signal_prl.tex" 581 lines, 22203 characters characteristics of pulsars. We find that an observable called the braking index should be easily measurable during the transition epoch and can have a value far removed (by orders of magnitude) from the canonical value of three expected for magnetic dipole radiation, and may have either sign. The duration of the transition epoch is governed by the slow loss of angular momentum to radiation and is further prolonged by the reduction in the moment of inertia caused by the phase change which can even introduce an era of spin-up. We estimate that about one in a hundred pulsars may be passing through this phase. The phenomenon is analogous to ``bachbending'' observed in the moment of inertia of rotating nuclei observed in the 1970's, which also signaled a change in internal structure with changing spin.

The modern nucleon-nucleon interaction of Myers and Swiatecki, adjusted to the properties of finite nuclei, the parameters of the mass formula, and the behavior of the optical potential is used to calculate the properties of $\beta$--equilibrated neutron star matter, and to study the impact of this equation of state on the properties of (rapidly rotating) neutron stars and their cooling behavior. The results are in excellent agreement with the outcome of calculations performed for a broad collection of sophisticated nonrelativistic as well as relativistic models for the equation of state.

For the past 20 years it had been thought that the coexistence phase of the confined hadronic and quark matter phases, assumed to be a first order transition, was strictly excluded from neutron stars. This, however, was due to a seemingly innocuous idealization which has approximated away important physics. The reason is that neutron stars constitute multi-component bodies rather than single-component ones formerly (and incorrectly) used to describe the deconfinement phase transition in neutron stars. So, contrary to earlier claims, `neutron' stars may very well contain quark matter in their cores surrounded by a mixed-phase region of quark and hadronic matter. Such objects are called hybrid stars. The structure of such stars as well as an observable signature that could signal the existence of quark matter in their cores are discussed in this paper.

The general statement that hypothetical strange (quark matter) stars cool more rapidly than neutron stars is investigated in greater detail. It is found that the direct Urca process could be forbidden not only in neutron stars but also in strange stars. In this case, strange stars are slowly cooling, and their surface temperatures are more or less indistinguishable from those of slowly cooling neutron stars. Furthermore the case of enhanced cooling is reinvestigated. It shows that strange stars cool significantly more rapidly than neutron stars within the first $\sim 30$ years after birth. This feature could become particularly interesting if continued observation of SN 1987A would reveal the temperature of the possibly existing pulsar at its center.

Neutrino emission from the dense hadronic component in neutron stars is subject to strong modifications due to collective effects in the nuclear medium. We implement new estimates of the neutrino emissivities of two processes operating in the nuclear medium into numerical cooling simulations of neutron stars. The first process is the modified Urca process, for which the softening of the pion exchange mode and other polarization effects as well as the neutrino emission arising from the intermediate reaction states are taken into account. The second process concerns neutrino emission through superfluid pair breaking and formation processes. It is found that the medium effects on the emissivity of the modified Urca process result in a strong density dependence, which gives a smooth crossover from the standard to the nonstandard cooling scenario for increasing star masses. For superfluid stars, the superfluid pair breaking and formation processes accelerate mildly both the standard and the nonstandard cooling scenario. This leads to a good agreement between the theoretical cooling tracks and the rather low temperatures observed for objects like PSRs 0833-45 (Vela), 0656+14, and 0630+18 (Geminga). The robustness of our findings against variations in both the underlying equation of state of baryonic matter and the used fast cooling processes is demonstrated. Hence we conclude that the two recalculated neutrino emissivities studied here enable one to reproduce theoretically most of the observed pulsar temperatures by varying the masses of neutron star models.

After giving an overview of the history and idea of neutron stars, I shall introduce, in part one of my lectures, a representative collection of models for the equation of state of superdense neutron star matter, which account for various physical possibilities concerning the unknown behavior of such matter. Part two of my lectures deal with the construction of models of static as well as rapidly rotating neutron stars, and an investigation of the cooling behavior of such objects. Both is performed in the framework of Einstein's theory of general relativity. For this purpose the broad collection of models for the equation of state, derived in part one, will be used as an input. Various computed neutron star properties (e.g., masses, radii, moments of inertia, red and blueshifts, limiting rotational periods, cooling behavior) will be compared with the body of observed pulsar data. Finally, part three deals with an investigation of the properties of stars made up of absolutely stable strange quark matter. Their properties are compared with those of their nonstrange counterparts, i.e., neutron stars.

This paper gives an overview of the properties of all possible equilibrium sequences of compact strange-matter stars with nuclear crusts, which range from strange stars to strange dwarfs. In contrast to their non-strange counterparts, --neutron stars and white dwarfs--, their properties are determined by two (rather than one) parameters, the central star density and the density at the base of the nuclear crust. This leads to stellar strange-matter configurations whose properties are much more complex than those of the conventional sequence. As an example, two generically different categories of stable strange dwarfs are found, which could be the observed white dwarfs. Furthermore we find very-low-mass strange stellar objects, with masses as small as those of Jupiter or even lighter planets. Such objects, if abundant enough in our Galaxy, should be seen by the presently performed gravitational microlensing searches.

This paper gives an overview of the properties of all possible equilibrium sequences of compact strange-matter stars with nuclear crusts, which range from strange stars to strange dwarfs. In contrast to their non-strange counterparts, --neutron stars and white dwarfs--, their properties are determined by two (rather than one) parameters, the central star density and the density at the base of the nuclear crust. This leads to stellar strange-matter configurations whose properties are much more complex than those of the conventional sequence. As an example, two generically different categories of stable strange dwarfs are found, which could be the observed white dwarfs. Furthermore we find very low-mass strange stellar objects, with masses as small as those of Jupiter or even lighter planets. Such objects, if abundant enough in our Galaxy, should be seen by the presently performed gravitational microlensing searches. Further aspects studied in this paper concern the limiting rotational periods and the cooling behavior of neutron stars and their strange counterparts.

A collection of modern, field-theoretical equations of state is applied to the investigation of cooling properties of compact stars. These comprise neutron stars as well as hypothetical strange matter stars, made up of absolutely stable 3-flavor strange quark matter. Various uncertainties in the behavior of matter at supernuclear densities, e.g., hyperonic degrees of freedom, behavior of coupling strengths in matter, pion and meson condensation, superfluidity, transition to quark matter, absolute stability of strange quark matter, and last but not least the many-body technique itself are tested against the body of observed cooling data.