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}$.

We present a Bayesian analysis of the Landau mass within the extended $\sigma$-$\omega$ model for neutron star matter. To this purpose, we consider the mass measurement of the object PSR 0740+6620, the tidal deformability estimation from the GW170817 and the mass-radius estimate of PSR J0030+0451 by NICER. Using Landau mass as free parameter of the theory, we rely on the prediction power of the Bayesian method to find the best value for this nuclear quantity.

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.

We present a new method to interpolate between two matter phases that allows for a description of mixed phases and can be used, e.g., for mimicking transitions between pasta structures occuring in the crust as well as in the inner core of compact stars. This interpolation method is based on assuming switch functions that are used to define a mixture of subphases while fulfilling constraints of thermodynamic stability. The width of the transition depends on a free parameter, the pressure increment relative to the critical pressure of a Maxwell construction. As an example we present a trigonometric function ansatz for the switch function together with a pressure increment during the transition. We note that the resulting mixed phase equation of state bears similarities with the appearance of substitutional compounds in neutron star crusts and with the sequence of transitions between different pasta phases in the hadron-to-quark matter transition. We apply this method to the case of a hadron-to-quark matter transition and test the robustness of the compact star mass twin phenomenon against the appearance of pasta phases modeled in this way.

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 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 present a systematic investigation of the cosmological constant effects in compact stars interiors in the framework of Einstein's gravity. Consideration of a cosmological constant $\Lambda$ in compact stars is motivated by the mechanism of acceleration of the observable universe, where $\Lambda$ is usually related to the dark energy. In particular, we consider compact star mass twins, hybrid neutron stars that populate both the second and third branch of the mass-radius diagram. For those models, the need of consideration of excluded volume effects in the equation of state, resulting from the finite size volume of nucleons, leads to a stiffening of matter causing compact stars to acquire higher mass and radius values. We demonstrate that certain values of the cosmological constant can also modify the compact star properties but in an opposite way. In addition, we find that the inclusion of $\Lambda$ can have a similar effect to the existence of pasta phases at the hadron-quark interface.

In the past few years, cosmic-rays beyond the GZK cut-off ($E > 5 \times 10^{19}$ eV) have been detected by leading collaborations such as Pierre Auger Observatory. Such observations raise many questions as to how such energies can be reached and what source can possibly produce them. Although at lower energies, mechanisms such as Fermi acceleration in supernovae front shocks seem to be favored, top-down scenarios have been proposed to explain the existence of ultra-high energy cosmic-rays: the decay of super-massive long-lived particles produced in the early Universe may yield to a flux of ultra-high energy photons. Such photons might be presently generating so called super-preshowers, an extended cosmic-ray shower with a spatial distribution that can be as wide as the Earth diameter. The Cosmic Ray Extremely Distributed Observatory (CREDO) mission is to find such events by means of a network of detectors spread around the globe. CREDO's strategy is to connect existing detectors and create a worldwide network of cosmic-ray observatories. Moreover, citizen-science constitutes an important pillar of our approach. By helping our algorithms to recognize detection patterns and by using smartphones as individual cosmic-ray detectors, non-scientists can participate in scientific discoveries and help unravel some of the deepest mysteries in physics.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may overspread a significant fraction of the Earth, is a scientific terra incognita. Instead of developing and commissioning a completely new global detector infrastructure, CREDO proposes approaching the global cosmic-ray analysis objectives with all types of available detectors, from professional to pocket size, merged into a worldwide network. With such a network it is possible to search for evidences of correlated cosmic-ray ensembles. One of the observables that can be investigated in CREDO is a number of spatially isolated events collected in a small time window which could shed light on fundamental physics issues. The CREDO mission and strategy requires active engagement of a large number of participants, also non-experts, who will contribute to the project by using common electronic devices (e.g. smartphones). In this note the status and perspectives of the project is presented.

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.

Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, that give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star--neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.

We present a multi-polytrope approach to describe high-mass twins fulfilling chiral effective field theory estimations of the neutron star equation state and test it against the appearance of mixed phases at the hadron-quark interface. In addition, we discuss astrophysical applications of this method and expected future measurements that shall further constrain neutron star matter and the understanding of the QCD phase diagram.

Recently it has been found that a certain class of hybrid star equations of state with a large latent heat (strong first order phase transition obtained by a Maxwell construction) between stiff hadronic hadronic and stiff quark matter phases allows for the appearance of a third family of compact stars (including "twins") at high mass of $\sim 2 M_\odot$. We investigate how robust this high-mass twin phenomenon is against a smoothing of the transition which would occur, e.g., due to pasta structures in the mixed phase. To this end we propose a simple construction of a pasta-like equation of state with a parameter that quantifies the degree of smoothing of the transition and could eventually be related to the surface tension of the pasta structures. It is interesting to note that the range of energy densities for the transition as well as the pressure at the onset of the transition of this class of hybrid star matter at zero temperature corresponds well to values of the same quantities found in finite temperature lattice QCD simulations for the 1 $\sigma$ region at the pseudocritical temperature $T_c=154 \pm 9$ MeV. The pattern of the speed of sound as a function of energy density is very different.

We show that in the 3-polytropes model of Hebeler et al. \cite{Hebeler:2013nza} for the neutron star equation of state at supersaturation densities a third family of compact stars can be obtained which confirms the possibility of high-mass twin stars that have coincident masses $M_1=M_2\approx 2~M_\odot$ and significantly different radii $|R_1-R_2|>2-3 $ km. We show that the causality constraint puts severe limitations on the maximum mass of the third family sequence which can be relaxed when this scheme is extended to four polytropes thus mimicking a realistic high-density matter EoS.

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 point out that the very existence of a "horizontal branch" in the mass-radius characteristics for neutron stars indicates a strong first-order phase transition and thus supports the existence of a critical endpoint (CEP) of first order phase transitions in the QCD phase diagram. This branch would sample a sequence of hybrid stars with quark matter core, leading to the endpoint of stable compact star configurations with the highest possible baryon densities. Since we know of the existence of compact stars with $2M_\odot$, this hypothetical branch has to lie in the vicinity of this mass value, if it exists. We report here a correlation between the maximal radius of the horizontal branch and the pressure at the onset of hadron-to-quark matter phase transition which is likely to be a universal quantity of utmost relevance to the upcoming experiments with heavy-ion collisions at NICA and FAIR.

Both the symmetry energy part and excluded volume corrections to the equation of state play an important role for the neutron star interior structure and composition, namely for the profile of the baryon density and the proton fraction. While the symmetry energy uniquely determines the proton fraction, excluded volume effects control the maximum density values inside neutron stars. Observations of cooling neutron stars indicate that the fast direct Urca cooling is not operative for the typical, low mass stars, pointing at proton fractions that lie below the threshold for the onset of direct Urca cooling process. This in turn, restricts the density range admissible in neutron star interiors and may require an excluded volume correction. In this contribution we discuss the interplay between fast cooling, symmetry energy and excluded volume corrections to the equation of state that would be required to fulfil the direct Urca cooling constraint.

We discuss the observation that under neutron star conditions of charge neutrality and $\beta-$equilibrium the contribution from the symmetry energy to the equation of state (EoS) follows a universal behaviour. We call this behaviour the conjecture of a Universal Symmetry Energy Contribution (USEC). We find that an USEC holds provided the density dependence of the symmetry energy $E_s(n)$ follows a behaviour that limits the proton fraction $x(n)$ to values below the threshold for the direct Urca (DU) cooling process. The absence of DU cooling in typical mass neutron stars appears to be supported by the phenomenology of neutron star cooling data and allows to constrain the behaviour of $E_s(n)$ at high densities. Two classes of symmetry energy functions are investigated more in detail to elucidate the USEC. We derive an analytic formula for the USEC to the neutron star EoS based on the result for the symmetry energy extracted from isobaric analog states of nuclei.

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.

In this work we study the case of deconfinement in compact star interiors in the presence of a strong first order phase transition associated to a critical end point in the QCD phase diagram. Neutron stars fulfilling these conditions show a third branch in the mass-radius diagram with the first and second branches being the white dwarfs and neutron stars configurations. The transition to the third branch can be reached by a pure hadronic neutron star through an induced collapse releasing energy that corresponds to a mass-energy difference between the second and third branch configurations. Physical outcomes of this phenomenon that can potentially explain the already detected astrophysical signals are discussed. In particular we present energy estimations for the case of a fast radio burst, seen as a double-peak structure in the object's light curve.

In this contribution, we outline a new 2-phase description of the quark-nuclear matter hybrid equation of state that takes into account effects of phase space occupation (excluded volume) in both, the hadronic and the quark matter phases. For the nuclear matter phase, the reduction of the available volume at increasing density leads to a stiffening, while for the quark matter phase a reduction of the effective string tension in the confining density functional is obtained. The deconfinement phase transition in the resulting hybrid equation of state is sensitive to both excluded volume effects. As an application, we consider matter under compact star constraints of electric neutrality and $\beta$-equilibrium. We obtain mass-radius relations for hybrid stars that fulfill the $2M_\odot$ constraint and exhibit the high-mass twin phenomenon. Both features depend sensitively on the excluded volume.

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.

We estimate the energy reservoir available in the deconfinement phase transition induced collapse of a neutron star to its hybrid star mass twin on the "third family" branch, using a recent equation of state of dense matter. The available energy corresponding to the mass-energy difference between configurations is comparable with energies of the most violent astrophysical burst processes. An observational outcome of such a dynamical transition might be fast radio bursts, specifically a recent example of a FRB with a double-peak structure in its light curve.

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.

Aims: We present a new microscopic hadron-quark hybrid equation of state model for astrophysical applications, from which compact hybrid star configurations are constructed. These are composed of a quark core and a hadronic shell with a first-order phase transition at their interface. The resulting mass-radius relations are in accordance with the latest astrophysical constraints. Methods: The quark matter description is based on a quantum chromodynamics (QCD) motivated chiral approach with higher-order quark interactions in the Dirac scalar and vector coupling channels. For hadronic matter we select a relativistic mean-field equation of state with density-dependent couplings. Since the nucleons are treated in the quasi-particle framework, an excluded volume correction has been included for the nuclear equation of state at suprasaturation density which takes into account the finite size of the nucleons. Results: These novel aspects, excluded volume in the hadronic phase and the higher-order repulsive interactions in the quark phase, lead to a strong first-order phase transition with large latent heat, i.e. the energy-density jump at the phase transition, which fulfils a criterion for a disconnected third-family branch of compact stars in the mass-radius relationship. These twin stars appear at high masses ($\sim$ 2 M$_\odot$) that are relevant for current observations of high-mass pulsars. Conclusions: This analysis offers a unique possibility by radius observations of compact stars to probe the QCD phase diagram at zero temperature and large chemical potential and even to support the existence of a critical point in the QCD phase diagram.

In order to prove the existence of a critical end point (CEP) in the QCD phase diagram it is sufficient to demonstrate that at zero temperature $T=0$ a first order phase transition exists as a function of the baryochemical potential $\mu$, since it is established knowledge from ab-initio lattice QCD simulations that at $\mu=0$ the transition on the temperature axis is a crossover. We present the argument that the observation of a gap in the mass-radius relationship for compact stars which proves the existence of a so-called third family (aka "mass twins") will imply that the $T=0$ equation of state of compact star matter exhibits a strong first order transition with a latent heat that satisfies $\Delta\epsilon/\epsilon_c \gt 0.6$. Since such a strong first order transition under compact star conditions will remain first order when going to symmetric matter, the observation of a disconnected branch (third family) of compact stars in the mass-radius diagram proves the existence of a CEP in QCD. For the equation of state of the twins the quark matter description is based on a QCD-motivated chiral approach with higher-order quark interactions in the Dirac scalar and vector coupling channels. For hadronic matter we select a relativistic mean-field equation of state with density-dependent couplings. Since the nucleons are treated in the quasi-particle framework, an excluded volume has been included for the nuclear equation of state at super-saturation density that takes into account the finite size of the nucleons. Furthermore we show results of a Bayesian analysis (BA) using disjunct M-R constraints for extracting probability measures for cold, dense matter equations of state. This study reveals that measuring the radii of neutron star twins has the potential to support the existence of a first order phase transition for compact star matter.

We aim at contributing to the resolution of three of the fundamental puzzles related to the still unsolved problem of the structure of the dense core of compact stars (CS): (i) the hyperon puzzle: how to reconcile pulsar masses of $2\,$M$_\odot$ with the hyperon softening of the equation of state (EoS); (ii) the masquerade problem: modern EoS for cold, high density hadronic and quark matter are almost identical; and (iii) the reconfinement puzzle: what to do when after a deconfinement transition the hadronic EoS becomes favorable again? We show that taking into account the compositeness of baryons (by excluded volume and/or quark Pauli blocking) on the hadronic side and confining and stiffening effects on the quark matter side results in an early phase transition to quark matter with sufficient stiffening at high densities which removes all three present-day puzzles of CS interiors. Moreover, in this new class of EoS for hybrid CS falls the interesting case of a strong first order phase transition which results in the observable high mass twin star phenomenon, an astrophysical observation of a critical endpoint in the QCD phase diagram.

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 study the possibility that the transition from hadron matter to quark matter at vanishing temperatures proceeds via crossover, similar to the crossover behavior found with lattice QCD studies at high temperatures. The purpose is to examine astrophysical consequences of this postulate by constructing hybrid star sequences fulfilling current experimental data.

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.

We present two types of models for hybrid compact stars composed of a quark core and a hadronic mantle with an abrupt first order phase transition at the interface which are in accordance with the latest astrophysical measurements of two 2 M_sun pulsars. While the first is a schematic one, the second one is based on a QCD motivated nonlocal PNJL model with density-dependent vector coupling strength. Both models support the possibility of so called twin compact stars which have the same mass but different radius and internal structure at high mass (~2 M_sun), provided they exhibit a large jump \Delta \epsilon in the energy density of the first order phase transition fulfilling \Delta \epsilon/\epsilon_crit > 0.6. We conclude that the measurement of high-mass twin stars would support the existence of a first order phase transition in symmetric matter at zero temperature entailing the existence of a critical end point in the QCD phase diagram.

We present a model for hybrid compact stars composed of a quark core and a hadronic mantle with an abrupt first order phase transition at the interface and in accordance with the latest astrophysical measurements of two 2 M_sun pulsars. We demonstrate the possibility of a disconnected mass-radius sequence (third family) of high-mass pulsars as a distinct feature due to a large jump $\Delta \epsilon$ in the energy density of the first order phase transition setting in at $\epsilon_{\rm crit}\approx 500$ MeV/fm$^3$ and fulfilling $\Delta \epsilon/\epsilon_{\rm crit} > 0.6$. We conclude that the measurement of so called \textit{twin} compact stars at high mass (~ 2 M_sun) would support the existence of a first order phase transition in symmetric matter at zero temperature entailing the existence of a critical end point (CEP) in the QCD phase diagram.

The functional form of the nuclear symmetry energy has only been determined in a very narrow range of densities. Uncertainties concern both the low as well as the high density behaviour of this function. In this work different shapes of the symmetry energy, consistent with the experimental data, were introduced and their consequences for the crustal properties of neutron stars are presented. The resulting models are in agreement with astrophysical observations.

The functional form of the nuclear symmetry energy in the whole range of densities relevant for the neutron stars is still unknown. Discrepancies concern both the low as well as the high density behaviour of this function. By use of Bezier curves three different families of the symmetry energy shapes, relevant for different density range were introduced. Their consequences for the crustal properties of neutron stars are presented.

Models for the internal composition of Dense Compact Stars are reviewed as well as macroscopic properties derived by observations of relativistic processes. Modeling of pure neutron matter Neutron Stars is presented and crust properties are studied by means of a two fluid model.

The form of the nuclear symmetry energy $E_s$ around saturation point density leads to a different crust-core transition point in the neutron star and affect the crust properties. We show that the knowledge about $E_s$ close to the saturation point is not sufficient, because the very low density behaviour is relevant. We also claim that crust properties are strongly influenced by the very high density behavior of $E_s$, so in order to conclude about the form of low density part of the symmetry energy one must isolate properly the high density part.