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.

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.

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.

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.

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.

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

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.