In this work we study the parameters of the extended $\sigma$-$\omega$ model for neutron star matter by a Bayesian analysis on state-of-the-art multi-messenger astronomy observations, namely mass, radius and tidal deformabilities. We have considered three parameters of the model, the Landau mass $m_L$, the nuclear compressibility $K_0$, and the value of the symmetry energy $S_0$, all at saturation density $n_0$. As a result, we are able to estimate the values of the Landau mass of f $m_L = 739\pm17$ MeV, whereas the values of $K_0$ and $S_0$ fall within already known empirical values. Furthermore, for neutron stars we find the most probable value of 13 km $<R_{1.4}<$ 13.5 km and the upper mass limit of $M_{max} \approx 2.2$ M$_{\odot}$.

The relativistic mean field (RMF) model of the nuclear matter equation of state has been modified by including the effect of Pauli-blocking owing to quark exchange between the baryons. Different schemes of a chiral enhancement of the quark Pauli blocking have been suggested according to the adopted density dependence of the dynamical quark mass. The resulting equations of state for the pressure are compared to the RMF model DD2 with excluded volume correction. On the basis of this comparison a density-dependent nucleon volume is extracted which parametrises the quark Pauli blocking effect in the respective scheme of chiral enhancement. The dependence on the isospin asymmetry is investigated and the corresponding density dependent nuclear symmetry energy is obtained in fair accordance with phenomenological constraints. The deconfinement phase transition is obtained by a Maxwell construction with a quark matter phase described within a higher order NJL model. Solutions for rotating and nonrotating (hybrid) compact star sequences are obtained which show the effect of high-mass twin compact star solutions for the rotating case.

We investigate the possibility that GW170817 has not been the merger of two conventional neutron stars (NS) but involved at least one if not two hybrid stars with a quark matter core which might even belong to a third family of compact stars. To this end, we develop a Bayesian analysis method for selecting the most probable equation of state (EoS) under a set of constraints from compact star physics, which now also include the tidal deformability from GW170817 and the first result for the mass and radius determination for PSR J0030+0451 by NICER. We apply this method for the first time to a two-parameter family of hybrid EoS based on the DD2 model with nucleonic excluded volume for hadronic matter and the color superconducting generalized nlNJL model for quark matter. The model has a variable onset of deconfinement and can mimic the effects of pasta phases with the possibility of a third family of hybrid stars in the mass-radius ($M-R$) diagram. The main findings of this study are that: 1) the presence of multiple configurations for a given mass (twins) corresponds to a set of disconnected lines in the diagram of tidal deformabilities for binary mergers, so that merger events from the same mass range may result in a probability landscape with different peak positions; 2) the Bayesian analysis with the above observational constraints favors an early onset of the deconfinement transition, at masses of $M_{\rm onset}\le 0.8~M_\odot$ with a $M-R$ relationship that in the range of observed neutron star masses is almost indistinguishable from that of a soft hadronic APR EoS; 3) a few yet fictitious measurements of the NICER experiment with a $1\sigma$ range that is half of the present value and different mass and radius would change the posterior likelihood so that hybrid EoS with a phase transition onset in the range $M_{\rm onset} = 1.1 - 1.6~M_\odot $ would be favored.

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.

The effect of pasta phases on the quark-hadron phase transition is investigated for a set of relativistic mean-field equations of state for both hadron and quark matter. The results of the full numerical solution with pasta phases are compared with those of an interpolating construction used in previous works, for which we demonstrate an adequate description of the numerical results. A one-to-one mapping of the free parameter of the construction to the physical surface tension of the quark-hadron interface is obtained for which a fit formula is given. For each pair of quark and hadron matter models the critical value of the surface tension is determined, above which the phase transition becomes close to the Maxwell construction. This result agrees well with earlier theoretical estimates. The study is extended to neutron star matter in beta equilibrium with electrons and muons and is applied to investigate the effect of pasta phases on the structure of hybrid compact stars and the robustness of a possible third family solution.

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.

We develop a Bayesian analysis method for selecting the most probable equation of state under a set of constraints from compact star physics, which now include the tidal deformability from GW170817. We apply this method for the first time to a two-parameter family of hybrid equations of state that is based on realistic models for the hadronic phase (KVORcut02) and the quark matter phase (SFM$\alpha$) which produce a third family of hybrid stars in the mass-radius diagram. One parameter ($\alpha$) characterizes the screening of the string tension in the string-flip model of quark matter while the other ($\Delta_P$) belongs to the mixed phase construction that mimics the thermodynamics of pasta phases and includes the Maxwell construction as a limiting case for $\Delta_P=0$. We present the corresponding results for compact star properties like mass, radius and tidal deformabilities and use empirical data for them in the newly developed Bayesian analysis method to obtain the probabilities for the model parameters within their considered range.

We demonstrate that the existing neutron-star cooling data can be appropriately described within "the nuclear medium cooling scenario" including hyperons under the assumption that different sources have different masses. We use a stiff equation of state of the relativistic mean-field model MKVORH$\phi$ with hadron effective couplings and masses dependent on the scalar field. It fulfills a large number of experimental constraints on the equation of state of the nuclear matter including the $2\,M_{\odot}$ lower bound for the maximum predicted neutron-star mass and the constraint for the pressure from the heavy-ion particle flow. We select appropriate $^1S_0$ proton and $\Lambda$ hyperon pairing gap profiles from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. The $^3P_2$ neutron pairing gap is assumed to be negligibly small in our scenario. The possibility of the pion, kaon and charged $\rho$-meson condensations is for simplicity suppressed. The resulting cooling curves prove to be sensitive to the value and the density dependence of the pp pairing gap and rather insensitive to the values of the $^1S_0$ neutron pairing gaps.

First-order phase transitions, like the liquid-gas transition, proceed via formation of structures such as bubbles and droplets. In strongly interacting compact star matter, at the crust-core transition, but also at the hadron-quark transition in the core, these structures form different shapes dubbed "pasta phases". We describe two methods to obtain one-parameter families of hybrid equations of state (EoS) which mimic the thermodynamic behavior of pasta phases in between a low-density hadron and a high-density quark matter phase, thus generalizing the Maxwell construction. The first method replaces the behavior of pressure vs. chemical potential in a finite region around the critical %chemical potential pressure of the Maxwell construction by a polynomial interpolation. The second method uses extrapolations of the hadronic and quark matter EoS beyond the Maxwell point to define a mixing of both with weight functions bounded by finite limits around the Maxwell point. We apply both methods to the case of a hybrid EoS with a strong first order transition that entails the formation of a third family of compact stars and the corresponding mass twin phenomenon. We investigate for both models the robustness of this phenomenon against variation of the single parameter, the pressure increment at the critical chemical potential which quantifies the deviation from the Maxwell construction. We also show sets of results for other compact star observables than mass and radius, namely the moment of inertia and the baryon mass.

We investigate the robustness of third family solutions for hybrid compact stars with a quark matter core that correspond to the occurrence of high-mass twin stars against a softening of the phase transition by means of a construction that mimics the effects of pasta structures in the mixed phase. We consider a class of hybrid equations of state that exploits a relativistic mean-field model for the hadronic as well as for the quark matter phase. We present parametrizations that correspond to branches of high-mass twin star pairs with maximum masses between $2.05~M_\odot$ and $1.48~M_\odot$ having radius differences between 3.2 km and 1.5 km, respectively. When compared to a Maxwell construction with a fixed value of critical pressure $P_c$ the effect of the mixed phase construction consists in the occurrence of a region of pressures around $P_c$ belonging to the coexistence of hadronic and quark matter phases between the onset pressure at $P_H$ and the end of the transition at $P_Q$. The maximum broadening which would still allow mass twin compact stars is found to be $(P_Q-P_H)_{\rm max} \approx P_c$ for all parametrizations within the present class of models. At least the heavier of the neutron stars of the binary merger GW170817 could have been a member of the third family of hybrid stars.

The existence of high mass ($\sim 2M_{\odot}$) pulsars PSR J1614-2230 and PSR J0348-0432 requires the compact star matter to be described by a stiff equation of state (EoS). Presence of hyperons in neutron stars leads to a softening of the EoS that results in a decrease of the maximum neutron-star mass below the measured values of masses for PSR J1614-2230 and PSR J0348-0432 pulsars if one exploits ordinary relativistic mean-field (RMF) models (hyperon puzzle). However, within an RMF EoS with a sigma-scaled hadron effective masses and coupling constants, the maximum neutron-star mass remains above $2M_{\odot}$ even when hyperons are included. Also, other important constraints on the equation of state, e.g. the flow constraint from heavy-ion collisions are to be fulfilled. We demonstrate how a satisfactory explanation of all existing observational data for the temperature-age relation is reached within the "nuclear medium cooling" scenario with a relativistic-mean-field EoS with a sigma-scaled hadron effective masses and coupling constants including hyperons.

A simple mixed phase model mimicking so-called "pasta" phases in the quark-hadron phase transition is developed and applied to static neutron stars for the case of DD2 type hadonic and NJL type quark matter models. The influence of the mixed phase on the mass-radius relation of the compact stars is investigated. Model parameters are chosen such that the results are in agreement with the observational constraints for masses and radii of pulsars.

Recent evidence for high masses ($\sim 2~ M_\odot$) pulsars PSR J1614-2230 and PSR J0348-0432 requires neutron star matter to have a stiff equation of state (EoS). The thermal evolution of compact stars (CS) with stiff hadronic EoS necessitates the application of the "nuclear medium cooling" scenario with a selection of appropriate proton gap profiles together with in-medium effects (like pion softening) on cooling mechanisms in order to achieve a satisfactory explanation of all existing observational data for the temperature-age relation of CS. Here we focus on two examples from \cite{Grigorian:2016leu} for a stiff hadronic EoS without (DD2 EoS) and with (DD2vex) excluded volume correction.

We present a Bayesian analysis for a new class of realistic models of two-phase equations of state (EoS) for hybrid stars and demonstrate that the observation of a pair of high-mass twin stars would have a sufficient discriminating power to favor hybrid EoS with a strong first order phase transition over alternative EoS. Such a measurement would provide evidence for the existence of a critical endpoint in the QCD phase diagram.

We explore systematically a new class of two-phase equations of state (EoS) for hybrid stars that is characterized by three main features : (1) stiffening of the nuclear EoS at supersaturation densities due to quark exchange effects (Pauli blocking) between hadrons, modelled by an excluded volume correction, (2) stiffening of the quark matter EoS at high densities due to multiquark interactions and (3) possibility for a strong first order phase transition with an early onset and large density jump. The third feature results from a Maxwell construction for the possible transition from the nuclear to a quark matter phase and its properties depend on the two parameters used for (1) and (2), respectively. Varying these two parameters one obtains a class of hybrid EoS that yields solutions of the Tolman-Oppenheimer-Volkoff (TOV) equations for sequences of hadronic and hybrid stars in the mass-radius diagram which cover the full range of patterns according to the Alford-Han-Prakash classification following which a hybrid star branch can be either absent, connected or disconnected with the hadronic one. The latter case often includes a tiny connected branch. The disconnected hybrid star branch, also called "third family", corresponds to high-mass twin stars characterized by the same gravitational mass but different radii. We perform a Bayesian analysis and demonstrate that the observation of such a pair of high-mass twin stars would have a sufficient discriminating power to favor hybrid EoS with a strong first order phase transition over alternative EoS.

Measurements demonstrate the existence of compact stars with masses in a broad range from 1.2 to 2 $M_\odot$. The most massive of these objects might be hybrid stars. To fulfill the constraint $M_{\rm max}>2~M_{\odot}$ with a reserve we exploit the stiff DD2 hadronic equation of state (EoS) without and with excluded volume (DD2vex) correction, which produce maximum neutron star masses of $M_{\rm max} = 2.43~M_{\odot}$ and $2.70~M_{\odot}$, respectively. We show that the stiffness of the EoS does not preclude an explanation of the whole set of cooling data within "nuclear medium cooling" scenario for compact stars by a variation of the star masses. We select appropriate proton gap profiles from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. However, we suppress the possibility of pion condensation. In general, the stiffer the EoS the steeper a decrease with density of the effective pion gap is required. Results are compared with previously obtained ones for the HDD EoS for which $M_{\rm max} =2.06~M_{\odot}$. The cooling of the compact star in the supernova remnant Cassiopeia A (Cas A) is explained mainly by an efficient medium modified Urca process. With the DD2vex EoS and using an effective pion gap steeper decreasing with the density and/or a proton gap shifted to smaller densities we are also able to reproduce both a strong decline compatible with ACIS-S data and HRC-S instrument data. The mass of Cas A is estimated as $1.6 - 1.9~M_{\odot}$, above the value $1.5~M_{\odot}$, which we have evaluated with the softer HDD equation of state. Different mass choices for the hottest object XMMU J173203.3-344518 are discussed. We make general remarks also on hybrid star cooling and on its dependence on the stiffness of the hadronic EoS.

We suggest a new Bayesian analysis (BA) using disjunct M-R constraints for extracting probability measures for cold, dense matter equations of state (EoS). One of the key issues of such an analysis is the question of a deconfinement transition in compact stars and whether it proceeds as a crossover or rather as a first order transition. We show by postulating results of not yet existing radius measurements for the known pulsars with a mass of $2M_\odot$ that a radius gap of about 3 km would clearly select an EoS with a strong first order phase transition as the most probably one. This would support the existence of a critical endpoint in the QCD phase diagram under scrutiny in present and upcoming heavy-ion collision experiments.

The recent measurements of the masses of the pulsar J00737-3039B and of the companion J1756-2251 and pulsars PSR J1614-2230, PSR J0348-0432 demonstrate the existence of compact stars with masses in a broad range from 1.2 to 2 $M_\odot$. To fulfill the constraint $M_{\rm max}>2M_{\odot}$ and to demonstrate the possibility of cooling scenarios for purely hadronic and further for hybrid stars we exploit the stiff DD2 hadronic equation of state producing a maximum neutron star mass $M\simeq 2.43 M_{\odot}$. We show that the "nuclear medium cooling" scenario for neutron stars comfortably explains the whole set of cooling curves just by a variation of the star masses without the necessity for the occurrence of the direct Urca reaction. To describe the cooling data with the very stiff DD2 equation of state we select a proton gap profile from those exploited in the literature and allow for a variation of the effective pion gap controlling the efficiency of the medium modified Urca process. Fast cooling of young neutron stars like it is seen in the data for Cas A is explained with the DD2 equation of state when the following conditions are provided: the presence of an efficient medium modified Urca process, and a large proton gap at densities $n\le 2n_0$ vanishing for $n\ge (2.5 - 3) n_0$, where $n_0$ is the saturation nuclear density.

We perform a Bayesian analysis in the basis of a recently developed two-parameter class of hybrid equations of state that allow for high-mass compact star twins. While recently a wide range of radii, from 9 - 15 km, has been inferred for different neutron stars using different techniques, we perform our analysis under the supposition that the radii are towards the large end ($13-15$ km). We use this radius constraint together with the undebated statistically independent constraint for high masses ($\sim 2~M_\odot$) as priors in selecting the most probable hybrid equations of state from a family with two free parameters: the baryon excluded volume in the hadronic phase and the 8-quark vector channel interaction in the quark matter phase.

Within the "nuclear medium cooling" scenario of neutron stars all reliably known temperature - age data, including those of the central compact objects in the supernova remnants of Cassiopeia A and XMMU-J1732, can be comfortably explained by a set of cooling curves obtained by variation of the star mass within the range of typical observed masses. The recent measurements of the high masses of the pulsars PSR J1614-2230 and PSR J0348-0432 on the one hand, and of the low masses for PSR J0737-3039B and the companion of PSR J1756-2251 on the other, provide independent proof for the existence of neutron stars with masses in a broad range from $\sim 1.2$ to 2 $M_\odot$. The values $M>2 M_{\odot}$ call for sufficiently stiff equations of state for neutron star matter. We investigate the response of the set of neutron star cooling curves to a stiffening of the nuclear equation of state so that maximum masses of about $2.4 M_\odot$ would be accessible and to a deconfinement phase transition from such stiff nuclear matter in the outer core to color superconducting quark matter in the inner core. Without readjustment of cooling inputs the mass range required to cover all cooling data for the stiff DD2 equation of state should include masses of $2.426 M_\odot$ for describing the fast cooling of CasA while the existence of a quark matter core accelerates the cooling so that CasA cooling data are described with a hybrid star of mass $1.674 M_\odot$.

We suggest a new Bayesian analysis using disjunct mass and radius constraints for extracting probability measures for cold, dense nuclear matter equations of state. One of the key issues of such an analysis is the question of a deconfinement transition in compact stars and whether it proceeds as a crossover or rather as a first order transition. The latter question is relevant for the possible existence of a critical endpoint in the QCD phase diagram under scrutiny in present and upcoming heavy-ion collision experiments.

We perform a Bayesian analysis of probability measures for compact star equations of state using new, disjunct constraints for mass and radius. The analysis uses a simple parametrization for hybrid equations of state to investigate the possibility of a first order deconfinement transition in compact stars. The latter question is relevant for the possible existence of a critical endpoint in the QCD phase diagram under scrutiny in heavy-ion collisions.

We suggest a new Bayesian analysis using disjunct M-R constraints for extracting probability measures for cold, dense matter equations of state. One of the key issues of such an analysis is the question of a deconfinement transition in compact stars and whether it proceeds as a crossover or rather as a first order transition. The latter question is relevant for the possible existence of a critical endpoint in the QCD phase diagram under scrutiny in present and upcoming heavy-ion collision experiments.

Recently, Elshamounty et al. performed a reanalysis of the surface temperature of the neutron star in the supernova remnant Cassiopeia A on the basis of Chandra data measured during last decade, and added a new data point. We show that all reliably known temperature data of neutron stars including those belonging to Cassiopea A can be comfortably explained in our "nuclear medium cooling" scenario of neutron stars. The cooling rates account for medium-modified one-pion exchange in dense matter, polarization effects in the pair-breaking-formation processes operating on superfluid neutrons and protons paired in the 1S_0 state, and other relevant processes. The emissivity of the pair-breaking-formation process in the 3P_2 state is a tiny quantity within our scenario. Crucial for a successful description of the Cassiopeia A cooling proves to be the thermal conductivity from both, the electrons and nucleons, being reduced by medium effects. Moreover, we exploit an EoS which stiffens at high densities due to an excluded volume effect and is capable of describing a maximum mass of 2.1 M_sun, thus including the recent measurements of PSR J1614-2230 and PSR J0348+0432.

The measurement of the mass 1.97 +/- 0.04 M_sun for PSR J1614-2230 provides a new constraint on the equation of state and composition of matter at high densities. In this contribution we investigate the possibility that the dense cores of neutron stars could contain strange quarks either in a confined state (hyperonic matter) or in a deconfined one (strange quark matter) while fulfilling a set of constraints including the new maximum mass constraint. We account for the possible appearance of hyperons within an extended version of the density-dependent relativistic mean-field model, including the phi-meson interaction channel. Deconfined quark matter is described by the color superconducting three-flavor NJL model.

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.

We discuss an idea for how accreting millisecond pulsars could contribute to the understanding of the QCD phase transition in the high-density nuclear matter equation of state (EoS). It is based on two ingredients, the first one being a ``phase diagram'' of rapidly rotating compact star configurations in the plane of spin frequency and mass, determined with state-of-the-art hybrid equations of state, allowing for a transition to color superconducting quark matter. The second is the study of spin-up and accretion evolution in this phase diagram. We show that the quark matter phase transition leads to a characteristic line in the Omega-M plane, the phase border between neutron stars and hybrid stars with a quark matter core. Along this line a change in the pulsar's moment of inertia entails a waiting point phenomenon in the accreting millisecond X-ray pulsar (AMXP) evolution: most of these objects should therefore be found along the phase border in the Omega-M plane, which may be viewed as the AMXP analog of the main sequence in the Hertzsprung-Russell diagram for normal stars. In order to prove the existence of a high-density phase transition in the cores of compact stars we need population statistics for AMXP's with sufficiently accurate determination of their masses and spin frequencies.

We study the structure and evolution of neutron stars (NS) the interiors of which are modeled using microscopic approaches and constrained by the condition that the equation of state (EoS) of matter extrapolated to high densities should not contradict known observational data from compact stars and experimental data from heavy-ion collisions (HIC). We use modern cooling simulations to extract distributions of NS masses required to reproduce those of the yet sparse data in the Temperature-Age (TA) plane. By comparing the results with a mass distribution for young, nearby NSs used in population synthesis we can sharpen the NS cooling constraints.

A new scheme for testing the nuclear matter (NM) equation of state (EoS) at high densities using constraints from compact star (CS) phenomenology is applied to neutron stars with a core of deconfined quark matter (QM). An acceptable EoS shall not to be in conflict with the mass measurement of 2.1 +/- 0.2 solar mass (1 sigma level) for PSR J0751+1807 and the mass radius relation deduced from the thermal emission of RX J1856-3754. Further constraints for the state of matter in CS interiors come from temperature-age data for young, nearby objects. The CS cooling theory shall agree not only with these data, but also with the mass distribution inferred via population synthesis models as well as with LogN-LogS data. The scheme is applied to a set of hybrid EsoS with a phase transition to stiff, color superconducting QM which fulfills all above constraints and is constrained otherwise from NM saturation properties and flow data of heavy-ion collisions. We extrapolate our description to low temperatures and draw conclusions for the QCD phase diagram to be explored in heavy-ion collision experiments.

We introduce a new tool for 'unmasking' the composition of neutron star (NS) interiors which is based on the fact that the state of matter at high densities determines the statistics of both NS observables, the temperature-age (TA) data as well as the mass distribution. We use modern cooling simulations to extract distributions of NS masses required to reproduce those of the yet sparse data in the TA plane. By comparing the results with a mass distribution for young, nearby NSs from population synthesis we can sharpen two NS cooling problems. The direct Urca (DU) problem consists in a narrowing of the NS population at the mass value for which the DU process as the most effective cooling mechanism in the hadronic layer of the star can occur. The Vela mass problem is a broadening of the population beyond the range of the typical mass window of 1.1 - 1.5 M_sun. Applying this tool to modern EoS we discuss examples for pure hadronic stars which are in conflict with these constraints while hybrid stars with a color superconducting quark matter core can predict a satisfactory mass distribution, provided the smallest diquark pairing gap has a properly defined density dependence.

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

A new scheme for testing the nuclear matter (NM) equation of state (EoS) at high densities using constraints from compact star (CS) phenomenology is applied to neutron stars with a core of deconfined quark matter (QM). An acceptable EoS shall not to be in conflict with the mass measurement of 2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807 and the mass-radius relation deduced from the thermal emission of RX J1856-3754. Further constraints for the state of matter in CS interiors come from temperature-age data for young, nearby objects. The CS cooling theory shall agree not only with these data, but also with the mass distribution inferred via population synthesis models as well as with LogN-LogS data. The scheme is applied to a set of hybrid EsoS with a phase transition to stiff, color superconducting QM which fulfills all above constraints and is constrained otherwise from NM saturation properties and flow data of heavy-ion collisions. We extrapolate our description to low temperatures and draw conclusions for the QCD phase diagram to be explored in heavy-ion collision experiments.

We discuss a set of tests which confront observations of cooling compact objects and theories of their thermal evolution. As an example we apply the recently developed $\mathrm{Log N}$-$\mathrm{Log S}$ test of compact star cooling theories to hybrid stars with a color superconducting quark matter core, we also apply and discuss other existing tests. While there is not yet a microscopically founded superconducting quark matter phase which would fulfill constraints from cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude and density-dependence of the X-gap can be chosen to satisfy a set of tests: the temperature -- age ($\mathrm{T}$-$\mathrm{t}$) test, the $\mathrm{Log N}$-$\mathrm{Log S}$ test, the brightness constraint, and the mass spectrum constraint. Some recent modifications of the population synthesis model used to obtain the $\mathrm{Log N}$-$\mathrm{Log S}$ distribution are briefly discussed. In addition, we propose to use the age-distance diagram as a new tool to study the local population of young isolated neutron stars.

We study the systematics of neutron star cooling curves with three representative masses from the most populated interval of the estimated mass distribution for compact objects. The cooling simulations are made in the framework of the nuclear medium cooling (NMC) scenario using different combinations of possible nucleon-nucleon pairing gaps. Possible heating or enhanced cooling mechanisms in the crust are not considered. We define a constraint on the highest possible temperatures for a given age of young neutron stars and show that this limits the freedom of modeling pairing gaps and crust properties.

We discuss three topics related to the neutron star (NS) mass spectrum. At first we discuss the possibility to form low-mass ($ M \stackrel{<}{\sim} 1 M_{\odot}$) and suggest this is possible only due to fragmentation of rapidly rotating proto-NSs. Such low-mass NSs should have very high spatial velocities which could allow identification. A critical assessment of this scenario is given. Secondly, we discuss mass growth due to accretion for NSs in close binary systems. With the help of numerical population synthesis calculations we derive the mass spectrum of massive ($M > 1.8 M_{\odot}$) NSs. Finally, we discuss the role of the mass spectrum in population studies of young cooling NSs. We formulate a kind of {\it mass constraint} which can be helpful, in our opinion, in discussing different competive models of the thermal evolution of NSs.

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.

We consider the early cooling evolution of strongly magnetized strange stars in a CFL phase with high gap Delta > 100 MeV. We demonstrate how this model may explain main features of the gamma-ray burst phenomena and also yield a strong star kick. The mechanism is based on beaming of neutrino emission along the magnetic vortex lines. We show that for sufficiently high initial temperatures T_0 ~ 30 - 60 MeV and surface magnetic fields B_s ~ 10^{15} - 10^{17} G, the energy release within the narrow beam is up to 10^{52} erg with a magnetic field dependent time scale between 10^{-2} s (for a smaller magnetic field) to 10 s. The above mechanism together with the parity violation of the neutrino-producing weak interaction processes in a magnetic field allow for the strange star kick. The higher the magnetic field the larger is the stars kick velocity. These velocities may cover the same range as observed pulsar kick velocities.

We apply the recently developed Log N - Log S test of compact star cooling theories for the first time to hybrid stars with a color superconducting quark matter core. While there is not yet a microscopically founded superconducting quark matter phase which would fulfill constraints from cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude and density-dependence of the X-gap can be chosen to satisfy a set of tests: temperature - age (T - t), the brightness constraint, Log N - Log S, and the mass spectrum constraint. The latter test appears as a new conjecture from the present investigation.

The study of thermal emission from isolated, cooling neutron stars plays a key role in probing the physical conditions of both the star crust and the core. The comparison of theoretical models for the star thermal evolution with the surface temperature derived from X-ray observations of sources of different age is one of the main tools to investigate the properties of the interior and constrain the equation of state. Here we propose to use population synthesis studies as an independent approach to test the physics governing the star cooling. Theoretical Log N - Log S distributions depend on the assumed neutron star thermal evolution. We have computed distributions for several different cooling scenarios and found that comparison with the observed Log N - Log S of isolated neutron stars is effective in discriminating among cooling models. Among the eleven cooling models considered in this paper, all of which may reproduce the observed temperature vs. age diagram, only at most three can explain the Log N - Log S distribution of close-by cooling neutron stars. The Log N - Log S test, being a ``global'' one and despite some limitations, appears capable of complementing the standard temperature vs. age test used up to now.

We study the dependence of the cooling of isolated neutron stars on the magnitude of the $3P_2$ neutron gap. It is demonstrated that our ``nuclear medium cooling'' scenario is in favor of a suppressed value of the $3P_2$ neutron gap.

The physics of compact objects is one of the very actively developing branches of theoretical investigations, since the careful analysis of different models for the internal structure of these objects, their evolutionary characteristics and the comparison with modern observational data give access to discriminate among various speculations about the state of matter under extreme conditions. The lecture provides an overview to the problem of neutron star cooling evolution. We discuss the scheme and necessary inputs for neutron star cooling simulations from the background of a microscopic modeling in order to present the surface temperature - age characteristics of hybrid neutron stars including the determination of main regulators of the cooling process, the question of neutrino production and diffusion.

We study the dependence of the cooling of isolated neutron stars on the magnitude of the $3P_2$ neutron gap. It is demonstrated that our ``Nuclear medium cooling scenario'' is in favor of a suppressed value of the $3P_2$ neutron gap.

The cooling of compact isolated objects for different values of the gravitational mass has been simulated for two alternative assumptions. One is that the interior of the star is purely hadronic and second that the star can have a rather large quark core. It has been shown that within a nonlocal chiral quark model the critical density for a phase transition to color superconducting quark matter under neutron star conditions can be low enough for these phases to occur in compact star configurations with masses below 1.3 M_sun. For a realistic choice of parameters the equation of state (EoS) allows for 2SC quark matter with a large quark gap ~ 100 MeV for u and d quarks of two colors that coexists with normal quark matter within a mixed phase in the hybrid star interior. We argue that, if in the hadronic phase the neutron pairing gap in 3P_2 channel is larger than few keV and the phases with unpaired quarks are allowed, the corresponding hybrid stars would cool too fast. Even in the case of the essentially suppressed 3P_2 neutron gap if free quarks occur for M < 1.3 M_sun, as it follows from our EoS, one could not appropriately describe the neutron star cooling data existing by today. It is suggested to discuss a "2SC+X" phase, as a possibility to have all quarks paired in two-flavor quark matter under neutron star constraints, where the X-gap is of the order of 10 keV - 1 MeV. Density independent gaps do not allow to fit the cooling data. Only the presence of an X-gap that decreases with increase of the density could allow to appropriately fit the data in a similar compact star mass interval to that following from a purely hadronic model.

We show that within a recently developed nonlocal chiral quark model the critical density for a phase transition to color superconducting quark matter under neutron star conditions can be low enough for these phases to occur in compact star configurations with masses below 1.3 M_solar. We study the cooling of these objects in isolation for different values of the gravitational mass. Our equation of state (EoS) allows for 2SC quark matter with a large quark gap \~100 MeV for u and d quarks of two colors that coexists with normal quark matter within a mixed phase in the hybrid star interior. We argue that, if the phases with unpaired quarks were allowed, the corresponding hybrid stars would cool too fast. If they occured for M < 1.3 M_solar, as it follows from our EoS, one could not appropriately describe the neutron star cooling data existing by today. We discuss a "2SC+X" phase, as a possibility to have all quarks paired in two-flavor quark matter under neutron star constraints, where the X-gap is of the order of 10 keV - 1 MeV. Density independent gaps do not allow to fit the cooling data. Only the presence of an X-gap that decreases with increase of the density could allow to appropriately fit the data in a similar compact star mass interval to that following from a purely hadronic model. This scenario is suggested as an alternative explanation of the cooling data in the framework of a hybrid star model.

We show that within a recently developed nonlocal, chiral quark model the critical densities for a phase transition to color superconducting quark matter under neutron star conditions can be low enough that these phases occur in compact star configurations with masses below 1.4 M_sun. We study the cooling of these objects in isolation for different values of the gravitational mass and thus different composition and structure of the interior. Our equation of state allows for a 2SC phase with a large quark gap Delta ~ 100 MeV for u and d quarks of two colors, a normal quark matter phase and their coexistence in a mixed phase within the hybrid star interior. We argue that, if the phases with unpaired quarks were allowed, the corresponding hybrid stars would cool too fast to describe the neutron star cooling data existing by today. We incorporate other attractive channels permitting a weak pairing of the residual quarks which remained unpaired in the 2SC phase and demonstrate that the model does not contradict the cooling data if the weak pairing gaps are of the order of 0.1 MeV.

The influence of neutrino trapping (NT) on the early cooling evolution of hot proto quark stars (PQS) with initial temperatures in the range T ~ 40 MeV is studied. Within a simplified model for the neutrino transport it is shown that the time for reaching neutrino opacity temperature of T_{opac} ~ 1 MeV is about 10 sec. This is an order of magnitude larger than without NT and of the same order as the duration long gamma ray bursts.

We study the cooling of isolated neutron stars. The main cooling regulators are introduced: EoS, thermal transport, heat capacity, neutrino and photon emissivity, superfluid nucleon gaps. Neutrino emissivity includes main processes. A strong impact of medium effects on the cooling rates is demonstrated. With taking into account of medium effects in reaction rates and in nucleon superfluid gaps modern experimental data can be well explained.

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

The equation of state for quark matter is derived for a nonlocal, chiral quark model within the mean field approximation.We investigate the effects of a variation of the formfactors of the interaction on the phase diagram of quark matter. Special emphasis is on the occurrence of a diquark condensate which signals a phase transition to color superconductivity and its effects on the equation of state under the condition of beta- equilibrium and charge neutrality. We calculate the quark star configurations by solving the Tolman- Oppenheimer- Volkoff equations and obtain for the transition from a hot, normal quark matter core of a protoneutron star to a cool diquark condensed one a release of binding energy of the order of Delta M c^2 ~ 10^{53} erg. We find that this energy could not serve as an engine for explosive phenomena since the phase transition is not first order. Contrary to naive expectations the mass defect increases when for a given temperature we neglect the possibility of diquark condensation.

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

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

For the classification of rotating compact stars with two high density phases a phase diagram in the angular velocity (Omega) - baryon number (N) plane is investigated. The dividing line N_crit(Omega) between configurations with one and two phases is correlated to a local maximum of the moment of inertia and can thus be subject to experimental verification by observation of the rotational behavior of accreting compact stars. Another characteristic line, which also can be measured is the transition line to black holes that of the maximum mass configurations. The positions and the shape of these lines are sensitive to changes in the equation of state (EoS) of stellar matter. A comparison of the regional structure of phase diagrams is performed for polytropic and relativistic mean field type EoS and correlations between the topology of the transition lines and the properties of two-phase EoS are obtained. Different scenarios of compact star evolution are discussed as trajectories in the phase diagram. It is shown that a population gap in the Omega - N plane for accreting compact stars would signal a high density phase transition and could be observed in the distribution of so called Z sources of quasi periodic oscillations in low-mass X-ray binaries.

Neutron star configurations are considered as thermodynamical systems for which a phase diagram in the angular velocity (Omega) - baryon number (N) plane is obtained with a dividing line N_{crit}(Omega) for quark core configurations. Trajectories of neutron star evolution in this diagram are studied for different scenarios defined by the external torque acting on the star due to radiation and/or mass accretion. They show a characteristic change in the rotational kinematics when the star enters the quark core regime. For isolated pulsars the braking index signal for deconfinement has been studied in its dependence on the mass of the star. Model calculations of the spin evolution of accreting low-mass X-ray binaries in the phase diagram have been performed for different values of the initial magnetic field, its decay time as well as initial mass and mass accretion rate. Population clustering of these objects at the line N_{crit}(Omega) in the phase diagram is suggested as an observable signal for the deconfinement phase transition if it exists for spinnning and accreting neutron stars.

We present the thermodynamics of a nonlocal chiral quark model with separable 4-fermion interaction for the case of U(3) flavor symmetry within a functional integral approach. The four free parameters of the model are fixed by the chiral condensate, and by the pseudoscalar meson properties (pion mass, kaon mass, pion decay constant). We discuss the T=0 equation of state (EoS) which describes quark confinement (zero quark matter pressure) below the critical chemical potential mu_c=333 MeV. The new result of the present approach is that the strange quark deconfinement is separated from the light quark one and occurs only at a higher chemical potential of mu_{c,s}=492 MeV. We compare the resulting EoS to bag model ones for two and three quark flavors, which have the phase transition to the vacuum with zero pressure also at mu_c. We study quark matter stars in general relativity theory assuming beta-equilibrium with electrons and show that for configurations with masses close to the maximum of stability at M=1.62 ... 1.64 M_{solar} strange quark matter can occur.

It has been suggested that the QPO phenomenon in LMXB's could be explained when the central compact object is a strange star. In this work we investigate within a standard model for disk accretion whether the observed clustering of spin frequencies in a narrow band is in accordance with this hypothesis. We show that frequency clustering occurs for accreting strange stars when typical values of the parameters of magnetic field initial strength and decay time, accretion rate are chosen. In contrast to hybrid star accretion no mass clustering effect is found.

We study the evolution of the rotation frequency for accreting compact stars. The discontinuous change of the moment of inertia of a rapidly rotating star due to the possible quark core appearance entails a characteristic change in the spin evolution. Numerical solutions have been performed using a model equation of state describing the deconfinement phase transition. Trajectories of spin evolution are discussed in the angular velocity - baryon number plane (phase diagram) for different accretion scenarios defined by the initial values of mass and magnetic field of the star, as well as mass accretion rate and magnetic field decay time. We observe a characteristic increase in the waiting time when a configuration enters the quark core regime. Overclustering of the population of Z sources of LMXBs in the phase diagram is suggested as a direct measurement of the critical line for the deconfinement phase transition since it is related to the behaviour of the moment of inertia of the compact star.

We study the consequences of superconducting quark cores (with color-flavor-locked phase as representative example) for evolution of temperature profiles and the cooling curves in quark-hadron hybrid stars and in hypothetical self-bounded objects having no a hadron shell (quark core neutron stars). The quark gaps are varied from 0 to $\Delta_q =50$ MeV. For hybrid stars we find time scales of $1\div5$, $5\div10$ and $50\div100$ years for the formation of a quasistationary temperature distribution in the cases $\Delta_q =0$, 0.1 MeV and $\gsim$ 1 MeV, respectively. These time scales are governed by the heat transport within quark cores for large diquark gaps ($\Delta \gsim$ 1 MeV) and within the hadron shell for small diquark gaps ($\Delta \lsim 0.1$ MeV). For quark core neutron stars we find a time scale $\simeq 300$ years for the formation of a quasistationary temperature distribution in the case $\Delta \gsim$ 10 MeV and a very short one for $\Delta \lsim$ 1 MeV. If hot young compact objects will be observed they can be interpreted as manifestation of large gap color superconductivity. Depending on the size of the pairing gaps, the compact star takes different paths in the ${lg}T_s $ vs. ${lg} t$ diagram where $T_s$ is the surface temperature. Compared to the corresponding hadronic model which well fits existing data the model for the hybrid neutron star (with a large diquark gap) shows too fast cooling. The same conclusion can be drawn for the corresponding self-bound objects.

Using the formalism of general relativity for axially symmetric gravitational fields and their sources - rotating compact stars - a perturbation theory with respect to angular velocity is developed and physical quantities such as mass, shape, momentum of inertia and total energy of the star are defined. The change of the internal structure of the star due to rotation has been investigated and the different contributions to the moment of inertia have been evaluated separately. Numerical solutions have been performed using a two-flavor model equation of state describing the deconfinement phase transition as constrained by the conservation of total baryon number and electric charge. During the spin down evolution of the rotating neutron star, below critical values of angular velocity a quark matter core can appear which might be detected as a characteristic signal in the pulsar timing. Within the spin-down scenario due to magnetic dipole radiation it is shown that the deviation of the braking index from $n=3$ could signal not only the occurrence but also the size of a quark core in the pulsar. A new scenario is proposed where, due to mass accretion onto the rapidly rotating compact star, a spin-down to spin-up transition might signal a deconfinement transition in its interior.

The occurence of a quark matter core in rotating compact stars has been investigated within general relativity as a function of both the rotational frequency and the total baryon number. We demonstrate that the deviation of the braking index from $n=3$ signals not only the occurence but also the size of a quark matter core in a pulsar. We suggest that in systems with mass accretion onto a rapidly rotating compact star a spin-down to spin-up transition might signal a deconfinement transition in its interior.

We explore the consequences of an equation of state (EOS) obtained in a confining Dyson-Schwinger equation model of QCD for the structure and stability of nonstrange quark stars at finite-T, and compare the results with those obtained using a bag-model EOS. Both models support a temperature profile that varies over the star's volume and the consequences of this are model independent. However, in our model the analogue of the bag pressure is (T,mu)-dependent, which is not the case in the bag model. This is a significant qualitative difference and comparing the results effects a primary goal of elucidating the sensitivity of quark star properties to the form of the EOS.