We study an impact of asymmetric dark matter on properties of the neutron stars and their ability to reach the two solar masses limit, which allows us to present a new range of masses of dark matter particles and their fractions inside the star. Our analysis is based on the observational fact of the existence of three pulsars reaching this limit and on the theoretically predicted reduction of the neutron star maximal mass caused by the accumulation of dark matter in its interior. We also demonstrate that light dark matter particles with masses below 0.2 GeV can create an extended halo around the neutron star leading not to decrease, but to increase of its visible gravitational mass. By using recent results on the spatial distribution of dark matter in the Milky Way, we present an estimate of its fraction inside the neutron stars located in the Galaxy center. We show how the detection of a 2Msun neutron star in the most central region of the Galaxy will impose an upper constraint on the mass of dark matter particles of ~ 60 GeV. Future high precision measurements of the neutron stars maximal mass near the Galactic center, will put a more stringent constraint on the mass of the dark matter particle. This last result is particularly important to prepare ongoing, and future radio and X-ray surveys.

In this work we investigate the possible condensation of tetraneutron resonant states in the lower density neutron rich gas regions inside Neutron Stars (NSs). Using a relativistic density functional approach we characterize the system containing different hadronic species including, besides tetraneutrons, nucleons and a set of light clusters ($^3$He, $\alpha$ particles, deuterium and tritium). $\sigma,\omega$ and $\rho$ mesonic fields provide the interaction in the nuclear system. We study how the tetraneutron presence could significantly impact the nucleon pairing fractions and the distribution of baryonic charge among species. For this we assume that they can be thermodynamically produced in an equilibrated medium and scan a range of coupling strengths to the mesonic fields from prescriptions based on isospin symmetry arguments. We find that tetraneutrons may appear over a range of densities belonging to the outer NS crust carrying a sizable amount of baryonic charge thus depleting the nucleon pairing fractions.

We revisit the Polyakov Loop coupled Nambu-Jona-Lasinio model that maintains the Polyakov loop dynamics in the limit of zero temperature. This is of interest for astrophysical applications in the interior of neutron stars. For this purpose we re-examine the form of the potential for the deconfinement order parameter at finite baryonic densities. Since the modification of this potential at any temperature is formally equivalent to assigning a baryonic charge to gluons, we develop a more general formulation of the present model that cures this spurious effect and is normalized to match the asymptotic behaviour of the QCD equation of state given by $\mathcal{O}(\alpha_s^2)$ and partial $\mathcal{O}(\alpha_s^3\ln^2\alpha_s)$ perturbative results.