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