The accurate determination of the maximum mass of the neutron stars is one of the most important tasks in astrophysics. It is directly related to the identification of the black holes in the universe, the production of neutron stars from the supernovae explosion, and the equation of state (EoS) of dense matter. However, not only the EoS is directly connected with neutron star masses, but also the speed of sound in dense matter is a crucial quantity which characterizes the stiffness of the EoS. The upper bound of the speed of sound imposes strong constraints on the maximum mass of neutron stars. However, this upper bound remains still an open issue. Recent observations, of binary neutron star systems, offer the possibility of measuring with high accuracy both the mass and the tidal polarizability of the stars. We study possible effects of the upper bound of the speed of sound on the upper bound of the mass and the tidal polarizability. We conclude that these kinds of measurements, combined with recent observations of neutron stars with masses close to $2 M_{\odot}$, will provide robust constraints on the equation of state of hadronic matter at high densities.

Nuclear matter and compact neutron stars are studied in the framework of the non-linear derivative (NLD) model which accounts for the momentum dependence of relativistic mean-fields. The generalized form of the energy-momentum tensor is derived which allows to consider different forms of the regulator functions in the NLD Lagrangian. The thermodynamic consistency of the NLD model is demonstrated for arbitrary choice of the regulator functions. The NLD approach describes the bulk properties of the nuclear matter and compares well with microscopic calculations and Dirac phenomenology. We further study the high density domain of the nuclear equation of state (EoS) relevant for the matter in $\beta$-equilibrium inside neutron stars. It is shown that the low density constraints imposed on the nuclear EoS and by the momentum dependence of the Schr\"odinger-equivalent optical potential lead to a maximum mass of the neutron stars around $M \simeq 2 M_{\odot}$ which accommodates the observed mass of the J1614-2230 millisecond radio pulsar.

We study the transition from hadronic matter to a mixed phase of quarks and hadrons at high baryon and isospin densities reached in heavy ion collisions. We focus our attention on the role played by the nucleon symmetry energy at high density.In this respect the inclusion of a scalar isovector meson, the \delta-coupling, in the Hadron Lagrangian appears rather important. We study in detail the formation of a drop of quark matter in the mixed phase, and we discuss the effects on the quark drop nucleation probability of the finite size and finite time duration of the high density region. We find that, if the parameters of quark models are fixed so that the existence of quark stars is allowed, then the density at which a mixed phase starts forming drops dramatically in the range Z/A \sim 0.3--0.4. This opens the possibility to verify the Witten-Bodmer hypothesis on absolute stability of quark matter using ground-based experiments in which neutron-rich nuclei are employed. These experiments can also provide rather stringent constraints on the Equation of State (EoS) to be used for describing the pre-Supernova gravitational collapse. Consistent simulations of neutron rich heavy ion collisions are performed in order to show that even at relatively low energies, in the few AGeV range, the system can enter such unstable mixed phase. Some precursor observables are suggested, in particular a ``neutron trapping'' effect.

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.