The detailed knowledge of the inner crust properties of neutron stars might be important to explain different phenomena such as pulsar glitches or the possibility of an {\it r-process} site in neutron star mergers. It has been shown in the literature that quantal effects like shell correction or pairing may play a relevant role to determine the composition of the inner crust of the neutron star. In this paper we construct the equation of state of the inner crust using the finite-range Gogny interactions, where the mean field and the pairing field are calculated with same interaction. We have used the semiclassical Variational Wigner-Kirkwood method along with shell and pairing corrections calculated with the Strutinsky integral method and the BCS approximation, respectively. Our results are compared with those of some popular models from the literature. We report a unified equation of state of the inner crust and core computed with the D1M* Gogny force, which was specifically fabricated for astrophysical calculations.

We perform a systematic study of the dependence of the r-mode phenomenology in normal fluid pulsar neutron stars on the symmetry energy slope parameter $L$. An essential ingredient in this study is the bulk viscosity, which is evaluated consistently for several equations of state corresponding to different values of the slope parameter $L$. Direct Urca processes, which are allowed from a critical $L$-value onwards, enhance the bulk viscosity and have large influence on the $r$-mode features, such as the instability boundary and spin-down properties of newborn neutron stars. The magnitude of the changes in the $r$-mode properties induced by the direct Urca processes are driven by the $L$-value of the equation of state and the mass of the pulsar. The study has been done by using a family of equations of state of $\beta$-equilibrated neutron star matter obtained with the finite range simple effective interaction, which provides realistic results for nuclear matter and finite nuclei properties. These equations of state predict the same properties in symmetric nuclear matter and have the same value of the symmetry energy parameter, $E_s(\rho_0)$, but differ in the slope parameter $L$. The range chosen for the variation of $L$ is decided from the tidal deformability data extracted from the GW170817 event and the maximum mass constraint.

We analyze the effects of including $\Delta(1232)$ isobars in an equation of state (EoS) for cold, $\beta$-stable neutron star matter, employing relativistic nuclear mean field theory. The selected EoS reproduces the properties of nuclear matter and finite nuclei and, in the astrophysical context, allows for the presence of hyperons in neutron stars having masses larger than 2$M_{\odot}$. We find that the composition and structure of neutron stars is critically influenced by the addition of the $\Delta$ isobars, which allows us to constrain their interaction with the meson fields taking into account astrophysical information. Imposing that the EoS is stable and ensures the existence of 2$M_{\odot}$ neutron stars, as well as requiring agreement with data of $\Delta$ excitation in nuclei, we find that, in the absence of other mechanisms stiffening the EoS at high densities, the interaction of the $\Delta$ isobars with the sigma and omega meson fields must be at least 10\% stronger than that of the nucleons. Moreover, the neutron star moment of inertia turns out to be sensitive to the presence of $\Delta$ isobars, whereas the inclusion of $\Delta$ isobars in the EoS allows for smaller stellar radii and for a lower value of the tidal deformability consistent with the analysis of the GW170817 merger event.

We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2$M_{\odot}$ neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above $1.8 M_{\odot}$.

We perform cooling simulations for isolated neutron stars using recently developed equations of state for their core. The equations of state are obtained from new parametrizations of the FSU2 relativistic mean-field functional that reproduce the properties of nuclear matter and finite nuclei, while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions, measurements of massive 2$M_{\odot}$ neutron stars, and neutron star radii below 13 km. We find that two of the models studied, FSU2R (with nucleons) and in particular FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including extensive nucleon pairing. This suggests that the cooling observations are more compatible with an equation of state that produces a soft nuclear symmetry energy and, hence, generates small neutron star radii. However, both models favor large stellar masses, above $1.8 M_{\odot}$, to explain the colder isolated neutron stars that have been observed, even if nucleon pairing is present.

The characteristic physical properties of rotating neutron stars under the r-mode oscillation are evaluated using the finite-range simple effective interaction. Emphasis is given on examining the influence of the stiffness of both the symmetric and asymmetric parts of the nuclear equation of state on these properties. The amplitude of the r-mode at saturation is calculated using the data of particular neutron stars from the considerations of "spin equilibrium" and "thermal equilibrium". The upper limit of the r-mode saturation amplitude is found to lie in the range 10^{-8}-10^{-6}, in agreement with the predictions of earlier work.

We study the symmetry energy and the core-crust transition in neutron stars using the finite-range Gogny nuclear interaction and examine the deduced crustal thickness and crustal moment of inertia. We start by analyzing the second-, fourth- and sixth-order coefficients of the Taylor expansion of the energy per particle in powers of the isospin asymmetry for Gogny forces. These coefficients provide information about the departure of the symmetry energy from the widely used parabolic law. The neutron star core-crust transition is evaluated by looking at the onset of thermodynamical instability of the liquid core. The calculation is performed with the exact (i.e., without Taylor expansion) Gogny EoS for the core, and also with its Taylor expansion in order to assess the influence of isospin expansions on locating the inner edge of neutron star crusts. It is found that the properties of the core-crust transition derived from the exact EoS differ from the predictions of the Taylor expansion even when the expansion is carried through sixth order in the isospin asymmetry. Gogny forces, using the exact EoS, predict the ranges $0.094 \text{ fm}^{-3} \lesssim \rho_t \lesssim 0.118\text{ fm}^{-3}$ for the transition density and $0.339 \text{ MeV fm}^{-3} \lesssim P_t \lesssim 0.665 \text{ MeV fm}^{-3}$ for the transition pressure. The transition densities show an anticorrelation with the slope parameter $L$ of the symmetry energy. The transition pressures are not found to correlate with $L$. Neutron stars obtained with Gogny forces have maximum masses below $1.74M_\odot$ and relatively small moments of inertia. The crustal mass and moment of inertia are evaluated and comparisons are made with the constraints from observed glitches in pulsars.

The D1 family of parametrizations of the Gogny interaction commonly suffers from a rather soft neutron matter equation of state that leads to maximal masses of neutron stars well below the observational value of two solar masses. We propose a reparametrization scheme that preserves the good properties of the Gogny force but allows one to tune the density dependence of the symmetry energy, which, in turn, modifies the predictions for the maximum stellar mass. The scheme works well for D1M, and leads to a new parameter set, dubbed D1M*. In the neutron-star domain, D1M* predicts a maximal mass of two solar masses and global properties of the star in harmony with those obtained with the SLy4 Skyrme interaction. By means of a set of selected calculations in finite nuclei, we check that D1M* performs comparably well to D1M in several aspects of nuclear structure in nuclei.

We reexamine the equation of state for the nucleonic and hyperonic inner core of neutron stars that satisfies the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km, while fulfilling the saturation properties of nuclear matter and finite nuclei together with the constraints on the high-density nuclear pressure coming from heavy-ion collisions. The recent nucleonic FSU2R and hyperonic FSU2H models are fine tuned improving the density dependence of pure neutron matter at subsaturation densities. The corresponding nuclear matter properties at saturation, the symmetry energy and its slope turn out to be compatible with recent experimental and theoretical determinations. We obtain the mass, radius and composition of neutron stars for the two updated models and study the impact on these properties of the uncertainties in the hyperon-nucleon couplings estimated from hypernuclear data. We find that the onset of appearance of each hyperon strongly depends on the hyperon-nuclear uncertainties, whereas the maximum masses for neutron stars differ by at most 0.1 $M_{\odot}$, although a larger deviation should be expected tied to the lack of knowledge of the hyperon potentials at the high densities present in the center of $2 M_\odot$ stars. For easier use, we provide tables with the results from the FSU2R and FSU2H models for the equation of state and the neutron star mass-radius relation.

We obtain a new equation of state for the nucleonic and hyperonic inner core of neutron stars that fulfills the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km. The nucleonic equation of state is obtained from a new parametrization of the FSU2 relativistic mean-field functional that satisfies these latest astrophysical constraints and, at the same time, reproduces the properties of nuclear matter and finite nuclei while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions. On the one hand, the equation of state of neutron star matter is softened around saturation density, which increases the compactness of canonical neutron stars leading to stellar radii below 13 km. On the other hand, the equation of state is stiff enough at higher densities to fulfill the 2$M_{\odot}$ limit. By a slight modification of the parametrization, we also find that the constraints of 2$M_{\odot}$ neutron stars with radii around 13 km are satisfied when hyperons are considered. The inclusion of the high magnetic fields present in magnetars further stiffens the equation of state. Hyperonic magnetars with magnetic fields in the surface of $ \sim 10^{15}$ G and with values of $\sim 10^{18}$ G in the interior can reach maximum masses of 2$M_{\odot}$ with radii in the 12-13 km range.

Although ab-initio calculations of relativistic Brueckner theory lead to large scalar isovector fields in nuclear matter, at present, successful versions of covariant density functional theory neglect the interactions in this channel. A new high precision density functional DD-ME$\delta$ is presented which includes four mesons $\sigma$, $\omega$, $\delta$, and $\rho$ with density dependent meson-nucleon couplings. It is based to a large extent on microscopic ab-initio calculations in nuclear matter. Only four of its parameters are determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, in particular the density dependence of the meson-nucleon vertices, are adjusted to non-relativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass $m_{p}^{\ast}-m_{n}^{\ast}$ derived from relativistic Brueckner theory is used to determine the coupling strength of the $\delta$-meson and its density dependence.

We derive a new equation of state (EoS) for neutron stars (NS) from the outer crust to the core based on modern microscopic Brueckner-Hartree-Fock (BHF) calculations using the Argonne $v_{18}$ potential plus three-body forces computed with the Urbana model. To deal with the inhomogeneous structures of matter in the NS crust, we use the recent Barcelona-Catania-Paris-Madrid (BCPM) nuclear energy density functional that is directly based on the same microscopic BHF calculations, and which is able to reproduce the ground-state properties of nuclei along the periodic table. The EoS of the outer crust requires the masses of neutron-rich nuclei, which are obtained through Hartree-Fock-Bogoliubov calculations with the BCPM functional when they are unknown experimentally. To compute the inner crust, Thomas-Fermi calculations in Wigner-Seitz cells are performed with the same functional. Existence of nuclear pasta is predicted in a range of average baryon densities between $\simeq$0.067 fm$^{-3}$ and $\simeq$0.0825 fm$^{-3}$, where the transition to the core takes place. The NS core is computed from the nuclear EoS of the BHF calculation assuming non-exotic constituents (core of $npe\mu$ matter). In each region of the star, we discuss the comparison of the new EoS with previous EoSes for the complete NS structure, in particular, with the Lattimer-Swesty EoS and with the Shen et al. EoS widely used in astrophysical calculations. The new microscopically derived EoS fulfills at the same time a NS maximum mass of 2~$M_\odot$ with a radius of 10 km, and a 1.5~$M_\odot$ NS with a radius of 11.7 km.

Within a microscopic approach the structure of Neutron Stars is usually studied by modelling the homogeneous nuclear matter of the core by a suitable Equation of State, based on a many-body theory, and the crust by a functional based on a more phenomenological approach. We present the first calculation of Neutron Star overall structure by adopting for the core an Equation of State derived from the Brueckner-Hartree-Fock theory and for the crust, including the pasta phase, an Energy Density Functional based on the same Equation of State, and which is able to describe accurately the binding energy of nuclei throughout the mass table. Comparison with other approaches is discussed. The relevance of the crust Equation of state for the Neutron Star radius is particularly emphasised.

The accurate characterization of the nuclear symmetry energy and its density dependence is one of the outstanding open problems in nuclear physics. A promising nuclear observable in order to constrain the density dependence of the symmetry energy at saturation is the neutron skin thickness of medium and heavy nuclei. Recently, a low-energy peak in the isovector dipole response of neutron-rich nuclei has been discovered that may be correlated with the neutron skin thickness. The existence of this correlation is currently under debate due to our limited experimental knowledge on the microscopic structure of such a peak. We present a detailed analysis of Skyrme Hartree-Fock (HF) plus random phase approximation (RPA) predictions for the dipole response in several neutron-rich nuclei and try to elucidate whether models of common use in nuclear physics confirm or dismiss its possible connection with the neutron skin thickness. Finally, we briefly present theoretical results for parity violating electron scattering on ${}^{208}$Pb at the conditions of the PREx experiment and discuss the implications for the neutron skin thickness of ${}^{208}$Pb and the slope of the symmetry energy.

We review and extend with nonrelativistic nuclear mean field calculations a previous study of the impact of the nuclear symmetry energy on the structure and composition of the outer crust of nonaccreting neutron stars. First, we develop a simple "toy model" to understand the most relevant quantities determining the structure and composition of the outer crust: the nuclear symmetry energy and the pressure of the electron gas. While the latter is a well determined quantity, the former ---specially its density dependence--- still lacks an accurate characterization. We thus focus on the influence of the nuclear symmetry energy on the crustal composition. For that, we employ different nuclear models that are accurate in the description of terrestrial nuclei. We show that those models with stiffer symmetry energies ---namely, those that generate thicker neutron skins in heavy nuclei and have smaller symmetry energies at subnormal nuclear densities--- generate more exotic isotopes in the stellar crust than their softer counterparts.

The symmetry energy coefficients, incompressibility, and single-particle and isovector potentials of clusterized dilute nuclear matter are calculated at different temperatures employing the $S$-matrix approach to the evaluation of the equation of state. Calculations have been extended to understand the aforesaid properties of homogeneous and clusterized supernova matter in the subnuclear density region. Comparison of the results in the $S$-matrix and mean-field approach reveals some subtle differences in the density and temperature region we explore.