We study the thermal evolution of neutron stars described within the equation of state with induced surface tension (IST) that reproduces properties of normal nuclear matter, fulfills the proton flow constraint, provides a high-quality description of hadron multiplicities created during the nuclear-nuclear collision experiments, and is equally compatible with the constraints from astrophysical observations and the GW170817 event. The model features strong direct Urca processes for the stars above 1.91\, M_{\odot}. The IST equation of state shows a very good agreement with the available cooling data, even without introducing nuclear pairing. We also analysed an effect of the singlet proton/neutron and triplet neutron pairing on the cooling of neutron stars of different mass. We demonstrate a full agreement of the predicted cooling curves with the experimental data. Moreover, the IST EoS provides a description of Cas A with both paired and unpaired matter.

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.

We study radial oscillations of non-rotating neutron stars (NSs) in four-dimensional General Relativity. The interior of the NS was modelled within a recently proposed multicomponent realistic equation of state (EoS) with the induced surface tension (IST). In particular, we considered the IST EoS with two sets of model parameters, that both reproduce all the known properties of normal nuclear matter, give a high quality description of the proton flow constraint, hadron multiplicities created in nuclear-nuclear collisions, consistent with astrophysical observations and the observational data from the NS-NS merger. We computed the 12 lowest radial oscillation modes, their frequencies and corresponding eigenfunctions, as well as the large frequency separation for six selected fiducial NSs (with different radii and masses of 1.2, 1.5 and 1.9 solar masses) of the two distinct model sets. The calculated frequencies show their continuous growth with an increase of the NS central baryon density. Moreover, we found correlations between the behaviour of first eigenfunction calculated for the fundamental mode, the adiabatic index and the speed of sound profile, which could be used to probe the internal structure of NSs with the asteroseismology data.

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.

We apply a novel equation of state (EoS) that includes the surface tension contribution induced by interparticle interaction and asymmetry between neutrons and protons, to the study of neutron star (NS) properties. This elaborated EoS is obtained from the virial expansion applied to multicomponent particle mixtures with hard core repulsion. The considered model is in full concordance with all the known properties of normal nuclear matter, provides a high-quality description of the proton flow constraints, hadron multiplicities created during the nuclear-nuclear collision experiments, and equally is consistent with astrophysical data coming from NS observations. The analysis suggests that the best model parameterization gives the incompressibility factor $K_{0}$, symmetry energy $J$, and symmetry energy slope $L$ at normal nuclear density equal to $200$ MeV, $30$ MeV, and $113.28-114.91$ MeV, respectively. The mass-radius relations found for NSs computed with this EoS are consistent with astrophysical observations.

In this work we discuss a novel approach to model the hadronic and nuclear matter equations of state using the induced surface tension concept. Since the obtained equations of state, classical and quantum, are among the most successful ones in describing the properties of low density phases of strongly interacting matter, they set strong restrictions on the possible value of the hard-core radius of nucleons. Therefore, we perform a detailed analysis of its value which follows from hadronic and nuclear matter properties and find the most trustworthy range of its values: the hard-core radius of nucleons is 0.30--0.36 fm. A comparison with the phenomenology of neutron stars implies that the hard-core radius of nucleons has to be temperature and density dependent.

A novel equation of state used for analysis of the heavy ion collision experimental data is generalized to also describe the matter inside neutron stars. This approach differs from others by including an induced surface tension caused by the short range repulsion between particles accounted by presence of their hard cores. This new equation of state respects thermodynamic consistency and is free from the usual causality problems. This makes it applicable over a wide range of temperatures and baryon densities, including the ones inside neutron stars. Despite small number of parameters the present equation of state reproduces properties of the saturated nuclear matter and is consistent with the proton flow constraint as well as with astrophysical data obtained from neutron star observations. Accordingly, we found the sets of parameters that agree with the same ones obtained from the heavy ion collision data analysis.

A novel equation of state with the surface tension induced by particles' interactions was generalized to describe the properties of the neutron stars (NSs). In this equation the interaction between particles occurs via the hard core repulsion by taking into account the proper volumes of particles. Recently, this model was successfully applied to the description of the properties of nuclear and hadron matter created in collisions of nucleons. The new approach is free of causality problems and is fully thermodynamically consistent, which enables us to use it for the investigation of the strongly interacting matter phase-diagram properties in a wide range of temperatures and baryon densities, including NSs. Here, we have calculated the mass-radius relations for a compact star using the Tolmann-Oppenheimer-Volkov (TOV) equation for two sets of parameters that satisfy the existing constraints. Accordingly, we found parameter values that are in good agreement with the same ones obtained from the nuclear-nuclear collision data analysis. The astrophysical constraints suggest that the hard core radius of baryons can vary between 0.425 fm and 0.476 fm.

Short range particle repulsion is rather important property of the hadronic and nuclear matter equations of state. We present a novel equation of state which is based on the virial expansion for the multicomponent mixtures with hard-core repulsion. In addition to the hard-core repulsion taken into account by the proper volumes of particles, this equation of state explicitly contains the surface tension which is induced by another part of the hard-core repulsion between particles. At high densities the induced surface tension vanishes and the excluded volume treatment of hard-core repulsion is switched to its proper volume treatment. Possible applications of this equation of state to a description of hadronic multiplicities measured in A+A collisions, to an investigation of the nuclear matter phase diagram properties and to the neutron star interior modeling are discussed.