Background: The microscopic composition and properties of infinite hadronic matter at a wide range of densities and temperatures have been subjects of intense investigation for decades. The equation of state (EoS) relating pressure, energy density, and temperature at a given particle number density is essential for modeling compact astrophysical objects such as neutron stars, core-collapse supernovae, and related phenomena, including the creation of chemical elements in the universe. The EoS depends not only on the particles present in the matter, but, more importantly, also on the forces acting among them. Because a realistic and quantitative description of infinite hadronic matter and nuclei from first principles in not available at present, a large variety of phenomenological models has been developed in the past several decades, but the scarcity of experimental and observational data does not allow a unique determination of the adjustable parameters.

The uncertainties in neutron star radii and crust properties due to our limited knowledge of the equation of state are quantitatively analyzed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core equation of state based on models with different properties at nuclear matter saturation, the uncertainties can be as large as $\sim 30$ % for the crust thickness and 4% for the radius. Necessary conditions for causal and thermodynamically consistent matchings between the core and the crust are formulated and their consequences examined. A large set of unified equations of state for purely nucleonic matter is obtained based on twenty-four Skyrme interactions and nine relativistic mean-field nuclear parametrizations. In addition, for relativistic models fifteen equations of state including a transition to hyperonic matter at high density are presented. All these equations of state have in common the property of describing a $2M_{\odot }$ star and of being causal within stable neutron stars. Spans of $\sim 3$ and $\sim 4$ km are obtained for the radius of, respectively, $1.0M_{\odot }$ and $2.0M_{\odot }$ stars. Applying a set of nine further constraints from experiment and ab initio calculations the uncertainty is reduced to $\sim 1$ and 2 km, respectively. These residual uncertainties reflect lack of constraints at large densities and insufficient information on the density dependence of the equation of state near the nuclear matter saturation point. The most important parameter to be constrained is shown to be the symmetry energy slope $L$. Indeed, this parameter exhibits a linear correlation with the stellar radius, which is particularly clear for small mass stars around $1.0M_{\odot }$. The other equation-of-state parameters do not show clear correlations with the radius, within the present uncertainties. Potential constraints on $L$, the neutron star radius, and the equation of state from observations of thermal states of neutron stars are also discussed. The unified equations of state are made available in the Supplemental Materials and via the CompOSE database.

Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and nonrelativistic mean-field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star $\approx 2M_{\odot }$. The tidal deformability shows a strong correlation with specific linear combinations of the isoscalar and isovector nuclear matter parameters associated with the EoS. Such correlations suggest that a precise value of the tidal deformability can put tight bounds on several EoS parameters, in particular on the slope of the incompressibility and the curvature of the symmetry energy. The tidal deformability obtained from the GW170817 and its UV, optical and infrared counterpart sets the radius of a canonical $1.4M_{\odot }$ neutron star to be $11.82\le R_{1.4}\le 13.72\mathrm{km}$.

In the present work we use the large-$N_c$ approximation to investigate quark matter described by the SU(2) Nambu–Jona-Lasinio model subject to a strong magnetic field. The Landau levels are filled in such a way that usual kinks appear in the effective mass and other related quantities. $\beta$ equilibrium is also considered and the macroscopic properties of a magnetar described by this quark matter is obtained. Our study shows that the magnetar masses and radii are larger if the magnetic field increases but only very large fields ($⩾{10}^{18}$ G) affect the equation of state in a non-negligible way.

The QCD phase diagram at zero chemical potential and finite temperature subject to an external magnetic field is studied within the three-flavor Nambu\char21{}Jona-Lasinio (NJL) model and the NJL model with the Polyakov loop. A scalar coupling parameter dependent on the magnetic field intensity is considered. The scalar coupling has been fitted so that the lattice QCD pseudocritical chiral transition temperatures are reproduced and in the limit of large magnetic field decreases with the inverse of the magnetic field intensity. This dependence of the coupling allows us to reproduce the lattice QCD results with respect to the quark condensates and Polyakov loop: due to the magnetic field the quark condensates are enhanced at low and high temperatures and suppressed for temperatures close to the transition temperatures and the Polyakov loop increases with the magnetic field.