Short range correlations (SRC) have been known to be an important aspect of nuclear theory for some time. Recent works have re-ignited interest on this topic, particularly due to the fact that it has recently been demonstrated that SRC may be responsible for breaking pairing gaps in nuclear matter. In this work we revisit the concept of SRC for beta equilibrated matter in neutron stars. We construct two equivalent models, with and without SRC and proceed to investigate the thermal evolution of stars described by such models. We show that SRC play a major role in the thermal evolution of neutron stars. It will be shown that while the SRC largely leaves the macroscopic properties of the star unaltered, it significantly alters the proton fraction, thus leading to an early onset of the direct Urca (DU) process, which in turns leads to stars exhibiting much faster cooling.

In this work we revisit the thermal relaxation process for neutron stars. Such process is associated with the thermal coupling between the core and the crust of neutron stars. The thermal relaxation, which takes place at around 10 -- 100 years, is manifested as a sudden drop of the star's surface temperature. Such drop is smooth for slow cooling objects and very sharp for those with fast cooling. In our study we focus particularly on the cooling of neutron stars whose mass is slightly greater than the value above which the direct Urca (DU) process sets in. Considering different mechanisms for neutrino production in each region of the star, and working with equations of state with different properties, we solve the thermal evolution equation and calculate the thermal relaxation time for ample range of neutron star masses. By performing a comprehensive study of neutron stars just above the onset of the direct Urca process we show that stars under these conditions exhibit a peculiar thermal relaxation behavior. We demonstrate that such stars exhibit an abnormally late relaxation time, characterized by a second drop of its surface temperature taking place a later ages. We qualify such behavior by showing that it is associated with limited spatial distribution of the DU process is such stars. We show that as the star's mass increase, the DU region also grows and the start exhibits the expected behavior of fast cooling stars. Finally we show that one can expect high relaxation times for stars in which the DU process takes place in a radius not larger than 3 km.

In the present work, we investigate the effects of short-range correlations (SRC) on the dimensionless deformability of the binary neutron system related to the GW170817 event. We implemented phenomenological SRC in a relativistic mean-field model in which the bulk parameters, namely, incompressibility ($K_0$), effective nucleon mass ratio ($m^*$), symmetry energy ($J$) and its slope ($L_0$), are independently controlled. Our results point out that the SRC favor the model to pass through the constraints, established by the LIGO/Virgo Collaboration, on the values of $\Lambda_{1.4}$ and on the $\Lambda_1\times\Lambda_2$ region. We also found a clear linear correlation between $\Lambda_{1.4}$ with $K_0$ and $L_0$ (increasing dependence), and with $m^*$ and $J$ (decreasing dependence). Finally, we also obtained compatible numbers for $R_{1.4}$ (model with and without SRC) in comparison with recent data from the Neutron Star Interior Composition Explorer (NICER) mission.

We propose a density-dependent function for the attractive interaction in the original van der Waals model to correctly describe the flow constraint at the high-density regime of the symmetric nuclear matter. After a generalization to asymmetric nuclear matter, it was also possible to study the stellar matter regime from this new model. The mass-radius relation for neutron stars under $\beta$-equilibrium is found to agree with recent X-ray observations. The neutron star masses supported against gravity, obtained from some parametrizations of the model, are in the range of $(1.97-2.07)M_{\odot}$, compatible with observational data from the PSR J0348+0432 pulsar. Furthermore, we verify the reliability of the model in predicting tidal deformabilities of the binary system related to the GW170817 neutron star merger event and find a full agreement with the new bounds obtained by the LIGO/Virgo collaboration.

In the present work we used five different versions of the quark-meson coupling (QMC) model to compute astrophysical quantities related to the GW170817 event and to neutron star cooling process. Two of the models are based on the original bag potential structure and three versions consider a harmonic oscillator potential to confine the quarks. The bag-like models also incorporate the pasta phase used to describe the inner crust of neutron stars. We show that the pasta phase always play a minor or negligible role in all studies. Moreover, while no clear correlation between the models that satisfy the GW170817 constraints and the slope of the symmetry energy is found, a clear correlation is observed between the slope and the fact that the cooling is fast or slow, i.e., fast (slow) cooling is related to higher (lower) values of the slope. We did not find one unique model that can describe, at the same time, GW170817 constraints and give a perfect description of the possible cooling processes.

We have obtained the Love number and corresponding tidal deformabilites ($\Lambda$) associated with the relativistic mean-field parametrizations shown to be consistent (CRMF) with the nuclear matter, pure neutron matter, symmetry energy and its derivatives [Dutra, et al., Phys. Rev. C 90, 055203 (2014)]. Our results show that CRMF models present very good agreement with the recent data from binary neutron star merger event GW170817. They also confirm the strong correlation between $\Lambda_{1.4}$ and the radius of canonical stars ($R_{1.4}$). When a recently GW170817 constraint on $\Lambda_{1.4}$ and the corresponding radius $R_{1.4}$ is used, the majority of the models tested are shown to satisfy it.

The present work starts by providing a clear identification of correlations between critical parameters ($T_c$, $P_c$, $\rho_c$) and bulk quantities at zero temperature of relativistic mean-field models (RMF) presenting third and fourth order self-interactions in the scalar field $\sigma$. Motivated by the nonrelativistic version of this RMF model, we show that effective nucleon mass ($M^*$) and incompressibility ($K_o$), at the saturation density, are correlated with $T_c$, $P_c$, and $\rho_c$, as well as, binding energy and saturation density itself. We verify agreement of results with previous theoretical ones regarding different hadronic models. Concerning recent experimental data of the symmetric nuclear matter critical parameters, our study allows a prediction of $T_c$, $P_c$ and $\rho_c$ compatible with such values, by combining them, through the correlations found, with previous constraints related to $M^*$ and $K_o$. An improved RMF parametrization, that better agrees with experimental values for $T_c$, is also indicated.

We have analyzed stellar properties of the relativistic mean-field (RMF) parametrizations shown to be consistent with the recently studied constraints related to nuclear matter, pure neutron matter, symmetry energy and its derivatives [Dutra et al., Phys. Rev. C 90, 055203 (2014)]. Our results show that only two RMF parametrizations do not allow the emergence of the direct Urca process, important aspect regarding the evolution of a neutron star. Moreover, among all approved RMF models, fourteen of them produce neutron stars with maximum masses inside the range $1.93\leqslant M/M_\odot\leqslant 2.05$, with $M_\odot$ being the solar mass. Only three models yield maximum masses above this range and a discussion on the inclusion of hyperons is presented. Finally, we have verified that the models satisfying the neutron star maximum mass constraint do not observe the squared sound velocity bound, namely, $v_s^2 < 1/3$, corroborating recent findings. However, the recently proposed $\sigma$-cut scheme can make the RMF models consistent with both constraints depending on the isoscalar-vector interaction of each parametrization.

In this work, we study the arising of correlations among some isoscalar ($K_o$, $Q_o$, and $I_o$) and isovector ($J$, $L_o$, $K_{sym}^o$, $Q_{sym}^o$, and $I_{sym}^o$) bulk parameters in nonrelativistic and relativistic hadronic mean-field models. For the former, we investigate correlations in Skyrme and Gogny parametrizations, as well as in the nonrelativistic (NR) limit of relativistic point-coupling models. We provide analytical correlations among bulk parameters for the NR limit, discussing the conditions in which they are linear ones. Based on a recent study [B. M. Santos et al., Phys. Rev. C 90, 035203 (2014)], we also show that some correlations presented in the NR limit are reproduced for relativistic models presenting cubic and quartic self-interactions in the scalar field $\sigma$, mostly studied in this work in the context of the relativistic framework. We also discuss how the crossing points, observed in the density dependence of some bulk parameters, can be seen as a signature of linear correlations between the specific bulk quantity presenting the crossing, and its immediately next order parameter.

Relativistic mean-field (RMF) models have been widely used in the study of many hadronic frameworks because of several important aspects not always present in nonrelativistic models, such as intrinsic Lorentz covariance, automatic inclusion of spin, appropriate saturation mechanism for nuclear matter, causality and, therefore, no problems related to superluminal speed of sound. With the aim of identifying the models which best satisfy well known properties of nuclear matter, we have analyzed $263$ parameterizations of seven different types of RMF models under three different sets of constraints related to symmetric nuclear matter, pure neutron matter, symmetry energy, and its derivatives. One of these (SET1) is formed of the same constraints used in a recent work [M. Dutra et al., Phys. Rev. C 85, 035201 (2012)] in which we analyzed $240$ Skyrme parameterizations. The results pointed to $2$ models consistent with all constraints. By using another set of constraints, namely, SET2a, formed by the updated versions of the previous one, we found $4$ models approved simultaneously. Finally, in the third set, named SET2b, in which the values of the constraints are more restrictive, we found $3$ consistent models. Another interesting feature of our analysis is that the results change dramatically if we do not consider the constraint regarding the volume part of the isospin incompressibility ($K_{\tau,\rm v}$). In this case, we have $35$ approved models in SET2a and $30$ in SET2b.

By using point-coupling versions of finite range nuclear relativistic mean field models containing cubic and quartic self interactions in the scalar field $\sigma$, a nonrelativistic limit is achieved. This approach allows an analytical expression for the symmetry energy ($J$) as a function of its slope ($L$) in a unified form, namely, $\,L\,=\,3J\,+f(m^{*},\rho_{o},B_{o},K_{o})$, where the quantities $m^{*}$, $\rho_{o}$, $B_{o}$ and $K_{o}$ are bulk parameters at the nuclear matter saturation density $\rho_{o}$. This result establishes a linear correlation between $L$ and $J$ which is reinforced by exact relativistic calculations. An analogous analytical correlation is also found for $J$, $L$ and the symmetry energy curvature ($K_{sym}$). Based on these results, we propose graphic constraints in $L\times J$ and $K_{sym}\times L$ planes which finite range models must satisfy.

We show that the magnitude of the order parameters in Polyakov-Nambu-Jona-Lasinio (PNJL) model, given by the quark condensate and the Polyakov loop, can be used as a criterium to clearly identify, without ambiguities, phases and boundaries of the strongly interacting matter, namely, the broken/restored chiral symmetry, and confinement/deconfinement regions. This structure is represented by the projection of the order parameters in the temperature-chemical potential plane, which allows a clear identification of pattern changes in the phase diagram. Such a criterium also enables the emergence of a quarkyonic phase even in the two-flavor system. We still show that this new phase diminishes due to the influence of an additional vector-type interaction in the PNJL phase diagrams, and is quite sensitive to the effect of the change of the $T_0$ parameter in the Polyakov potential. Finally, we show that the phases and boundaries constructed by our method indicate that the order parameters should be more strongly correlated, as in the case of entanglement PNJL (EPNJL) model. This result suggests a novel way to pursue further investigation of new interactions between the order parameters in order to improve the PNJL model.

This paper presents a detailed assessment of the ability of the 240 Skyrme interaction parameter sets in the literature to satisfy a series of criteria derived from macroscopic properties of nuclear matter in the vicinity of nuclear saturation density at zero temperature and their density dependence, derived by the liquid drop model, experiments with giant resonances and heavy-ion collisions. The objective is to identify those parameterizations which best satisfy the current understanding of the physics of nuclear matter over a wide range of applications. Out of the 240 models, only 16 are shown to satisfy all these constraints. Additional, more microscopic, constraints on density dependence of the neutron and proton effective mass beta-equilibrium matter, Landau parameters of symmetric and pure neutron nuclear matter, and observational data on high- and low-mass cold neutron stars further reduce this number to 5, a very small group of recommended Skyrme parameterizations to be used in future applications of the Skyrme interaction of nuclear matter related observables. Full information on partial fulfillment of individual constraints by all Skyrme models considered is given. The results are discussed in terms of the physical interpretation of the Skyrme interaction and the validity of its use in mean-field models. Future work on application of the Skyrme forces, selected on the basis of variables of nuclear matter, in Hartree-Fock calculation of properties of finite nuclei, is outlined.

In this work we study the hadron-quark phase transition matching relativistic hadrodynamical mean-field models (in the hadronic phase) with the more updated versions of the Polyakov-Nambu-Jona-Lasinio models (on the quark side). Systematic comparisons are performed showing that the predicted hadronic phases of the matching named as RMF-PNJL, are larger than the confined phase obtained exclusively by the Polyakov quark models. This important result is due to the effect of the nuclear force that causes more resistance of hadronic matter to isothermal compressions. For sake of comparison, we also obtain the matchings of the hadronic models with the MIT bag model, named as RMF-MIT, showing that it presents always larger hadron regions, while shows smaller mixed phases than that obtained from the RMF-PNJL ones. Thus, studies of the confinement transition in nuclear matter, done only with quark models, still need nuclear degrees of freedom to be more reliable in the whole $T\times\mu$ phase diagram.