Search

The implications of including the scalar isovector $\delta$-meson in a relativistic mean-field description of nuclear matter are discussed. A Bayesian inference approach is used to determine the parameters that define the isovector properties of the model. The properties of nuclear matter and neutron stars are discussed. The inclusion of the $\delta$-meson has only a small effect on the maximum mass of the neutron star (NS) and on the speed of sound in its interior, but it has a strong effect on the radius and the tidal deformability of low and medium mass stars. This is mainly due to the effect of the $\delta$-meson on the symmetry energy and its slope and curvature at saturation, increasing the range of possible values of these three properties, and in particular allowing positive values of the symmetry energy curvature. Due to the effect of the $\delta$-meson on the symmetry energy, the proton content of the star is also strongly affected. The inclusion of the $\delta$-meson in the relativistic mean-field description of nuclear matter extends the phase space spanned by the model, allowing for a more flexible density dependence of the symmetry energy compatible with experimental, observational, and ab initio constraints.

The CompactObject package is an open-source software framework developed to constrain the neutron star equation of state (EOS) through Bayesian statistical inference. It integrates astrophysical observational constraints from X-ray timing, gravitational wave events, and radio measurements, as well as nuclear experimental constraints derived from perturbative Quantum Chromodynamics (pQCD) and Chiral Effective Field Theory ($\chi$EFT). The package supports a diverse range of EOS models, including meta-model like and several physics-motivated EOS models. It comprises three independent components: an EOS generator module that currently provides seven EOS choices, a Tolman-Oppenheimer-Volkoff (TOV) equation solver, that allows the determination of the Mass Radius and Tidal deformability as observables, and a comprehensive Bayesian inference workflow module, including a complete pipeline for implementing EOS Bayesian inference. Each component can be used independently in different scientific research contexts, such as nuclear physics and astrophysics. In addition, CompactObject is designed to work in synergy with existing software such as CompOSE, allowing the use of the CompOSE EOS database to extend the EOS options available., Comment: submitted to JOSS, Github page: https://github.com/ChunHuangPhy/CompactObject , Documentation: https://chunhuangphy.github.io/CompactObject/

A covariant relativistic approach based on the Vlasov equation is applied to the study of infinite asymmetric nuclear matter. We use several Walecka-type hadronic models and obtain the dispersion relations for the longitudinal modes. The isovector and isoscalar collective modes are determined for a wide range of densities as a function of isospin asymmetry and momentum transfer within a set of eleven relativistic mean field models with different nuclear matter properties. Special attention is given to beta-equilibrium matter. It is shown that the possible propagation of isoscalar and isovector-like modes depends directly on the density dependence of the symmetric nuclear matter equation of state and of the symmetry energy, with a stiff equation of state favouring the propagation of isoscalar like collective modes at high densities, and a stiff symmetry energy defining the behavior of the isovector like modes which propagate for densities below two times saturation density. The coupling of the nuclear modes to the electron plasmon is also discussed., Comment: 18 pages, 11 figures, 3 tables, submitted to Phys Rev C

We investigate whether measurements of the neutron star mass and radius or the tidal deformability can provide information about the presence of hyperons inside a neutron star. This is achieved by considering two inference models, with and without hyperons, based on a field-theoretical approach. While current observations do not distinguish between the two scenarios, we have shown that data simulating expected observations from future large area X-ray timing telescopes could provide some information through Bayes factors. Inference using simulated data generated from an EOS containing hyperons decisively favours the hyperonic model over the nucleonic model. However, a 2\% uncertainty in the mass and radius determination may not be sufficient to constrain the parameters of the model when only six neutron star mass-radius measurements are considered., Comment: Submitted to MNRAS

We present a neural network classification model for detecting the presence of hyperonic degrees of freedom in neutron stars. The models take radii and/or tidal deformabilities as input and give the probability for the presence of hyperons in the neutron star composition. Different numbers of observations and different levels of uncertainty in the neutron star properties are tested. The models have been trained on a dataset of well-calibrated microscopic equations of state of neutron star matter based on a relativistic mean-field formalism. Real data and data generated from a different description of hyperonic matter are used to test the performance of the models., Comment: 13 pages, 8 figures

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon-oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB observational properties and taxonomy. We perform three-dimensional, smoothed-particle-hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state (EOS), with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, $t_{\rm col}$. We include a non-zero angular momentum of the NSs and find that $t_{\rm col}$ ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about five minutes orbital period), promptly form a BH and release $\gtrsim 10^{52}$ erg. They form NS-BH binaries with tens of kyr merger timescale by gravitational-wave emission. BdHNe II and III do not form BHs, release $\sim 10^{50}$-$10^{52}$ erg and $\lesssim 10^{50}$ erg. They form NS-NS binaries with a range of merger timescales larger than for NS-BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a non-rotating NS. Magnetic braking by a $10^{13}$ G field can delay BH formation, leading to BH-BH or NS-BH of tens of kyr merger timescale., Comment: Version accepted for publication in The Astrophysical Journal

The equation of state (EOS) for neutron stars is modeled using the Relativistic Mean Field (RMF) approach with a mesonic nonlinear (NL) interaction, a modified sigma cut potential (NL-$\sigma$ cut), and the influences of dark matter in the NL (NL DM). Using a Bayesian analysis framework, we evaluate the plausibility and impact of each scenario. Experimental constraints on the general properties of finite nuclei and heavy ion collisions, along with astrophysical observational data on neutron star radii and tidal deformation, have been taken into account. It was shown that all models, including the PREX-II data, were less favored, indicating that this experimental data seemed to be in tension with the other constraints included in the inference procedure, and were incompatible with chiral effective field theoretical calculations of pure neutron matter. Considering the models with no PREX-II constraints, we find the model NL-$\sigma$ cut with the largest Bayes evidence, indicating that the constraints considered favor the stiffening of the EOS at large densities. Conversely, the neutron star with a dark matter component is the least favorable case in light of recent observational constraints, among different scenarios considered here. The $f$ and $p$ modes were calculated within the Cowling approximation, and it can be seen that $f$ modes are sensitive to the EOS. An analysis of the slopes of the mass-radius curves and $f$-mode mass curves has indicated that these quantities may help distinguish the different scenarios.We also analyzed the impact of new PSR J0437-4715 measurements on neutron star mass-radius estimates, noting a $\sim$ 0.2 km reduction in the 90\% CI upper boundary across all models and a significant Bayes evidence decrease, indicating potential conflicts with previous data or the necessity for more adaptable models., Comment: 18 pages, 10 figures and 6 tables

A microscopic nuclear matter formalism with explicit chiral symmetry based on the Nambu Jona-Lasinio model is considered to describe nuclear matter. To reproduce nuclear matter properties adequately at the saturation density, four-point and eight-point interactions are introduced. Within a Bayesian inference approach, the parameters of the model are determined by imposing nuclear matter, both experimental and from {\it ab-initio} calculations, and neutron star observational constraints. Nuclear matter properties are well reproduced with an effective mass of 0.75 to 0.8 nucleon mass at the saturation density. At 90% confidence level, the radius of a $1.4 ~\rm M_\odot$ star varies between 11.48 km and 13.20 km, masses as large as $\sim 2.2 ~\rm M_\odot$ are predicted and the radius of a 2 M$_\odot$ star is above 10.5 km. High-density perturbative QCD (pQCD) results exclude equations of state that predict larger maximum masses and radii. The speed of sound increases monotonically with density and reaches values as large as $\sqrt{0.7}c$-$\sqrt{0.8}c$ in the center of massive stars. Several properties such as the polytropic index or the renormalized trace anomaly, that have been proposed to identify the deconfined phase transition, are analyzed. Interestingly, the radius of the obtained posterior that also meets pQCD constraints aligns closely with the mass-radius measurement of the recent PSR J0437-4715, which contrasts with other relativistic mean field model results., Comment: 15 pages, 7 figures

The density dependence of the symmetry energy in relativistic mean-field models with density dependent couplings is discussed in terms of the possible opening of nucleonic direct Urca processes inside neutron stars, which induce a very rapid cooling of the star. The modification of the parametrization of the isospin channel of two models, DD2 and DDMEX, keeping the same isoscalar properties is considered and the implications are discussed. Within the models discussed it is not possible the onset of nucleonic direct Urca processes in stars with a mass below $\sim1.6\,M_\odot$ if chiral effective field theory constraints for neutron matter are imposed. A Bayesian inference calculation confirms the low probability that nucleonic direct Urca processes occur inside stars with masses below 1.8$M_\odot$, considering the isoscalar channel of the equation of state described by DD2 or DDMEX and the same symmetry energy at saturation. The lowest masses allowing direct Urca processes are associated with a slope of the symmetry energy above $60$ MeV and most likely a positive symmetry energy incompressibility. It is shown that the parametrization of the isospin channel proposed destroys the correlation between symmetry energy slope and incompressibility previously identified in several works., Comment: 18 pages, 7 figures, 6 tables, accepted for publication in Phys. Rev. C

In this work, we study the effect of strong magnetic fields on the equation of state (EoS) of warm, homogeneous, Neutron Star (NS) matter in beta equilibrium. NS matter is described within a relativistic mean field (RMF) approximation, including both models with non-linear meson terms or with density dependent nucleon-meson couplings. We first study the effect of magnetic fields and finite temperature on the EoS separately, finding that the effect of the latter to be significantly stronger than the one of the former. We then study the combined effect of magnetic fields and temperature on the internal composition. We show how both factors cause an increase in the proton fraction at low density and that, as long as the temperatures considered are not higher than 10 MeV, the effect of the magnetic field on the proton fraction is not small enough to be neglected., Comment: 10 pages, 11 figures 1 table. Submitted to Phys. Rev. C

We propose a Bayesian inference estimation of in-medium modification of the cluster self-energies from light nuclei multiplicities measured in selected samples of central $^{136,124}$Xe$+^{124,112}$Sn collisions with the INDRA apparatus. The data are interpreted with a relativistic quasi-particle cluster approach in the mean-field approximation without any prior assumption on the thermal parameters of the model. An excellent reproduction is obtained for H and He isotope multiplicities, and compatible posterior distributions are found for the unknown thermal parameters. We conclude that the cluster-$\sigma$-meson coupling is temperature dependent, becoming weaker when the temperature increases, in agreement with microscopic quantum statistical calculations. This implies a faster decrease of the light cluster abundances with temperature than previously estimated.

The composition of the core of neutron stars is still under debate. One possibility is that because of the high densities reached in their cores, matter could be deconfined into quark matter. The possible existence of hybrid stars is investigated using microscopic models to describe the different phases of matter. Within these microscopic models we aim at calculating the properties of neutron stars and matter. We want to probe the influence of pQCD calculations and analyze the properties that identify a transition to deconfined matter. A Bayesian approach using a Markov Chain Monte Carlo sampling process is applied to generate 8 sets of equations of state. A Maxwell construction describes the deconfinement transition. For the hadronic phase, we consider a stiff and a soft EOS obtained from the Relativistic Mean Field model with non-linear meson terms. For the quark phase, we use two different models: the NJL model with multiquark interactions and the Mean Field Theory of QCD, similar to the MIT bag model with a vector term. The model parameters were determined by Bayesian inference imposing observations from NICER and the phase transition density range. We have also applied restrictions from the pQCD calculations. Hybrid stars are compatible with current observational data. The restrictions of pQCD reduce the value of the maximum mass. However, even when applying this restriction, the models were able to reach 2.1 to 2.3 solar masses. The conformal limit was not attained at the center of the most massive stars. The vector interactions are essential to describe hybrid stars with a mass above two solar masses. The multiquark interactions affect the limits of some quantities considered as indicators of the presence of a deconfined phase. It is possible to find a set of EOS, that predict that inside NS the renormalized matter trace anomaly is always positive., Comment: There are 17 pages, 11 figures

We analyse the implications of information about local derivatives from the mass-radius diagram in neutron star matter. It is expected that the next generation of gravitational wave and electromagnetic detectors will allow the determination of the neutron star radius and mass with a small uncertainty. Observations of neutron stars clustered around a given neutron star mass allow the estimation of local derivatives in the $M(R)$ diagram, which can be used to constrain neutron star properties. From a model-independent description of the neutron star equation of state, it is shown that a $M(R)$ curve with a negative slope at 1.4$M_\odot$ predicts a $2M_\odot$ neutron star radius below 12 km. Furthermore, a maximum mass below 2.3$M_\odot$ is obtained if the $M(R)$ slope is negative in the whole range of masses above $1M_\odot$, and a maximum mass above 2.4$M_\odot$ requires the $M(R)$ slope to be positive in some range of masses. Constraints on the mass-radius curve of neutron stars will place strong constraints on microscopic models., Comment: 10 pages, 11 figures, title changed, published version

We investigate how vector-isoscalar and vector-isovector interactions can be determined within the density regime of neutron stars (NSs), while fulfilling nuclear and astrophysics constrains. We make use of the Chiral Mean Field (CMF) model, a SU(3) nonlinear realization of the sigma model within the mean-field approximation, for the first time within a Bayesian analysis framework. We show that neutron-matter $\chi$EFT constraints at low density are only satisfied if the vector-isovector mixed interaction term is included, e.g., a $\omega^2\rho^2$ term. We also show the behavior of the model with respect to the conformal limit. We demonstrate that the CMF model is able to predict a value for the parameter $d_c$ related to the trace anomaly and its derivative takes values below 0.2 above four times saturation density within a hadronic version of the model that does not include hyperons or a phase transition to deconfined matter. We compare these effects with results from other (non-chiral) Relativistic Mean Field models to assess how different approaches to incorporating the same physical constraints affect predictions of NS properties and dense matter equations of state. We also include data from the gravitation wave event GW230529 detected by the LIGO-Virgo-Kagra collaboration and the most recent radius measurement of PSR J0437-4715 from the NASA NICER mission. Our analysis reveals that this new NICER measurement leads to an average reduction of approximately $\sim 0.1$ km radius in the posterior of the NS mass-radius relationship., Comment: 8 pages, 5 figures

We explore the anomaly detection framework based on Normalizing Flows (NF) models introduced in \cite{PhysRevC.106.065802} to detect the presence of a large (destabilising) dense matter phase transition in neutron star (NS) observations of masses and radii, and relate the feasibility of detection with parameters of the underlying mass-radius sequence, which is a functional of the dense matter equation of state. Once trained on simulated data featuring continuous $M(R)$ solutions (i.e., no phase transitions), NF is used to determine the likelihood of a first-order phase transition in a given set of $M(R)$ observations featuring a discontinuity, i.e., perform the anomaly detection. Different mock test sets, featuring two branch solutions in the $M(R)$ diagram, were parameterized by the NS mass at which the phase transition occurs, $M_c$, and the radius difference between the heaviest hadronic star and lightest hybrid star, $\Delta R$. We analyze the impact of these parameters on the NF performance in detecting the presence of a first-order phase transition. Among the results, we report that given a set of 15 stars with radius uncertainty of $0.2$ km, a detection of a two-branch solution is possible with 95\% accuracy if $\Delta R > 0.4$ km., Comment: 9 pages, 8 figures

Recently, in Malik23, a dataset of several EoS for purely nucleonic stellar matter based on a non-linear RMF model prescription, and constrained to properties of nuclear matter, to state-of-the-art chiral effective field theory calculations for low-density neutron matter, and to astrophysical data, were proposed. In this work, twenty one unified neutron star EoS were chosen from that dataset, in such a way that a large range of values of the slope of the symmetry energy at saturation is covered. Several quantities are calculated and discussed, such as the the proton fraction and the direct Urca behavior, the density dependence of the speed of sound and the trace anomaly, the crust-core transition properties, the compatibility with astrophysical observations, and the neutron matter properties from $\chi$EFT calculations and pQCD constraints. We construct unified EoS, where the outer crust is given by the BSk22 functional, and the inner crust is calculated from a CLD approximation. The core is purely nucleonic, made of protons, neutrons, electrons and muons, under charge neutrality and in $\beta-$equilibrium conditions. The correlation of the slope of the symmetry energy at saturation with the crust-core transition density and proton fraction is analysed, and equations that translate these relations are proposed. Moreover, the spectral representation for all the EOS according to the format proposed in Lindblom10 is given, which is a convenient representation to study quasi-periodic oscillations with realistic EOS. It is shown that several of these EoS have in the center of the most massive NS a speed of sound squared of the order of $\lesssim 0.5$. Most of the EoS predict a maximum central density of the order of about 6 times the nuclear saturation density. Three of the EoS satisfy all of the constraints imposed. All these EoS will be made available in the CompOSE platform., Comment: 14 pages, 12 figures, 7 tables

This study is devoted to the inference problem of extracting the nuclear matter properties directly from a set of mass-radius observations. We employ Bayesian neural networks (BNNs), which is a probabilistic model capable of estimating the uncertainties associated with its predictions. To simulate different noise levels on the $M(R)$ observations, we create three different sets of mock data. Our results show BNNs as an accurate and reliable tool for predicting the nuclear matter properties whenever the true values are not completely outside the training dataset statistics, i.e., if the model is not heavily dependent on its extrapolating capacities. Using real mass-radius pulsar data, the model predicted, for instance, $L_{\text{sym}}=39.80\pm17.52 $ MeV and $K_{\text{sym}}=-101.67\pm62.86 $ MeV ($2\sigma$ interval). Our study provides a valuable inference framework when new NS data becomes available., Comment: 15 pages, 12 figures

The equations of state (EoSs) governing neutron star (NS) matter obtained for both non-relativistic and relativistic mean-field models are systematically confronted with a diverse set of terrestrial data and astrophysical observations within the Bayesian framework. The terrestrial data, spans from bulk properties of finite nuclei to the heavy-ion collisions, constrain the symmetric nuclear matter EoS and the symmetry energy up to twice the saturation density ($\rho_0$= 0.16 fm$^{-3}$). The astrophysical observations encompass the NS radius, the tidal deformability, and the lower bound on maximum mass. Three distinct posterior distributions of EoSs are generated by gradually updating the priors with different constraints: (i) only the maximum NS mass, (ii) incorporating additional terrestrial data, (iii) combining both the terrestrial data and astrophysical observations. These EoS distributions are then compared using the Kullback-Liebler divergence which highlights the significant constraints imposed on the EoSs by the currently available lower bound of NS maximum mass and terrestrial data. The remaining astrophysical observations marginally refine the EoS within the density range $\sim$ 2-3$\rho_0$. It is observed that the relativistic mean field model yields stiffer EoS around the saturation density, but predict smaller values of the speed of sound and proton fraction in the interior of massive stars.

We investigate the impact of asymmetric fermionic dark matter (DM) on the thermal evolution of neutron stars (NSs), considering a scenario where DM interacts with baryonic matter (BM) through gravity. Employing the two-fluid formalism, our analysis reveals that DM accrued within the NS core exerts an inward gravitational pull on the outer layers composed of BM. This gravitational interaction results in a noticeable increase in baryonic density within the core of the NS. Consequently, it strongly affects the star's thermal evolution by triggering an early onset of the direct Urca (DU) processes, causing an enhanced neutrino emission and rapid star cooling. Moreover, the photon emission from the star's surface is modified due to a reduction of radius. We demonstrate the effect of DM gravitational pull on nucleonic and hyperonic DU processes that become kinematically allowed even for NSs of low mass. We then discuss the significance of observing NSs at various distances from the Galactic center. Given that the DM distribution peaks toward the Galactic center, NSs within this central region are expected to harbor higher fractions of DM, potentially leading to distinct cooling behaviors., Comment: 21 pages, 6 figures

We present a combination of two quick guides aimed at summarizing relevant information about the CompOSE nuclear equation of state repository. The first is aimed at nuclear physicists and describes how to provide standard equation of state tables. The second quick guide is meant for users and describes the basic procedures to obtain customized tables with equation of state data. Several examples are included to help providers and users to understand and benefit from the CompOSE database., Comment: Slightly updated from published version: Section 2.2.3. Stellar Information table has one extra line and optional additional columns

We propose the 3-momentum sphere intersection regularization applied to the one and two fermion line integrals at finite temperature and chemical potential. The quark-antiquark polarization function in this new regularization approach is equivalent to the usual 3-momentum regularization, when the absolute value of the external 3-momentum of the polarization is zero. Additionally, it respects the particle-antiparticle symmetry of meson states in the Nambu$-$Jona-Lasinio (NJL) model for all values of temperature and chemical potential. Without this symmetry, in-medium cross sections calculated in the 3-momentum regularized NJL model are not consistent. In order to demonstrate the difference between the usual 3-momentum regularization with the one proposed in this work, we study the quark-quark and quark-antiquark cross sections in both regularization schemes. To this end we use the standard $SU(3)$ NJL model, with four and six quark interactions. We observe major quantitative and qualitative differences when comparing quark-quark cross sections in both schemes. The quark-antiquark cross sections, on the other hand, are very similar in both regularizations, owning to the equivalence between the regularizations when the absolute value of the external 3-momentum is zero., Comment: 65 pages, 25 figures

We study the effect of asymmetric fermionic dark matter (DM) on the thermal evolution of neutron stars (NSs). No interaction between DM and baryonic matter is assumed, except the gravitational one. Using the two-fluid formalism, we show that DM accumulated in the core of a star pulls inwards the outer baryonic layers of the star, increasing the baryonic density in the NS core. As a result, it significantly affects the star's thermal evolution by triggering an early onset of the direct Urca process and modifying the photon emission from the surface caused by the decrease of the radius. Thus, due to the gravitational pull of DM, the direct Urca process becomes kinematically allowed for stars with lower masses. Based on these results, we discuss the importance of NS observations at different distances from the Galactic center. Since the DM distribution peaks towards the Galactic center, NSs in this region are expected to contain higher DM fractions that could lead to a different cooling behavior., Comment: 10 pages, 7 figures

The current observational properties of neutron stars have not definitively ruled out the possibility of dark matter. In this study, we primarily focus on exploring correlations between the dark matter model parameters and different neutron star properties using a rich set of EOSs. We adopt a two-fluid approach to calculate the properties of neutron stars. For the nuclear matter EOS, we employ several realistic EOS derived from the relativistic mean field model (RMF), each exhibiting varying stiffness and composition. In parallel, we look into the dark matter EOS, considering fermionic matter with repulsive interaction described by a relativistic mean field Lagrangian. A reasonable range of parameters is sampled meticulously. Interestingly, our results reveal a promising correlation between the dark matter model parameters and stellar properties, particularly when we ignore the uncertainties in the nuclear matter EOS. However, when introducing uncertainties in the nuclear sector, the correlation weakens, suggesting that the task of conclusively constraining any particular dark matter model might be challenging using global properties alone, such as mass, radius, and tidal deformability. Notably, we find that dark-matter admixed stars tend to have higher central baryonic density, potentially allowing for non-nucleonic degrees of freedom or direct Urca processes in stars with lower masses. There is also a tantalizing hint regarding the detection of stars with the same mass but different surface temperatures, which may indicate the presence of dark matter. With our robust and extensive dataset, we delve deeper and demonstrate that even in the presence of dark matter, the semi-universal C-Love relation remains intact., Comment: 15 Pages, 12 figures, 4 tables

Recently, a resonant state of four neutrons (tetraneutron) with an energy of $E_{4n}=2.37\pm 0.38 \rm{(stat)} \pm 0.44 \rm{(sys)}$ MeV and a width of $\Gamma=1.75\pm 0.22 \rm{(stat)} \pm 0.30 \rm{(sys)}$ MeV was reported. In this work, we analyse the effect of including such an exotic state on the yields of other light clusters, that not only form in astrophysical sites, such as core-collapse supernovae and neutron star mergers, but also in heavy-ion collisions. To this aim, we use a relativistic mean-field formalism, where we consider in-medium effects in a two-fold way, via the couplings of the clusters to the mesons, and via a binding energy shift, to compute the low-density equation of state for nuclear matter at finite temperature and fixed proton fraction. We consider five light clusters, deuterons, tritons, heliums, $\alpha$-particles, and $^6$He, immersed in a gas of protons and neutrons, and we calculate their abundances and chemical equilibrium constants with and without the tetraneutron. We also analyse how the associated energy of the tetraneutron would influence such results. We find that the low-temperature, neutron-rich systems, are the ones most affected by the presence of the tetraneutron, making neutron stars excellent environments for their formation. Moreover, its presence in strongly asymmetric matter may increase considerably the proton and the $\alpha$-particle fractions. This may have an influence on the dissolution of the accretion disk of the merger of two neutron stars., Comment: 10 pages, 6 figures

The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NS), is studied within a relativistic mean field description of nuclear matter. Different formulations, both with density dependent couplings and with non-linear mesonic terms, are considered and their predictions compared and discussed. A special attention is drawn to the effect on the neutron star properties of the inclusion of exotic degrees of freedom as hyperons. Properties such as the speed of sound, the trace anomaly, the proton fraction and the onset of direct Urca processes inside neutron stars are discussed. The knowledge of the general behavior of the hadronic equation of state and the implication it has on the neutron star properties will allow to identify signatures of a deconfinement phase transition discussed in other studies., Comment: 16 pages, 11 figures

Once further confirmed in future analyses, the radius and mass measurement of HESS J1731-347 with $M=0.77^{+0.20}_{-0.17}~M_{\odot}$ and $R=10.4^{+0.86}_{-0.78}~\rm km$ will be among the lightest and smallest compact objects ever detected. This raises many questions about its nature and opens up the window for different theories to explain such a measurement. In this article, we use the information from Doroshenko et al. (2022) on the mass, radius, and surface temperature together with the multimessenger observations of neutron stars to investigate the possibility that HESS J1731-347 is one of the lightest observed neutron star, a strange quark star, a hybrid star with an early deconfinement phase transition, or a dark matter-admixed neutron star. The nucleonic and quark matter are modeled within realistic equation of states (EOSs) with a self-consistent calculation of the pairing gaps in quark matter. By performing the joint analysis of the thermal evolution and mass-radius constraint, we find evidence that within a 1$\sigma$ confidence level, HESS J1731-347 is consistent with the neutron star scenario with the soft EOS as well as with a strange and hybrid star with the early deconfinement phase transition with a strong quark pairing and neutron star admixed with dark matter., Comment: 9 pages, 3 figures

The non-radial oscillations of the neutron stars (NSs) have been suggested as a useful tool to probe the composition of neutron star matter (NSM). With this scope in mind, we consider a large number of equations of states (EOSs) that are consistent with nuclear matter properties and pure neutron matter EOS based on a chiral effective field theory (chEFT) calculation for the low densities and perturbative QCD EOS at very high densities. This ensemble of EOSs is also consistent with astronomical observations, gravitational waves in GW170817, mass and radius measurements from Neutron star Interior Composition ExploreR (NICER). We analyze the robustness of known universal relations (URs) among the quadrupolar $f$ mode frequencies, masses and radii with such a large number of EOSs and we find a new UR that results from a strong correlation between the $f$ mode frequencies and the radii of NSs. Such a correlation is very useful in accurately determining the radius from a measurement of $f$ mode frequencies in the near future. We also show that the quadrupolar $f$ mode frequencies of NS of masses 2.0 M$_\odot$ and above lie in the range $\sim$ 2-3 kHz in this ensemble of physically realistic EOSs. A NS of mass 2M$_{\odot}$ with a low $f$ mode frequency may indicate the existence of non-nucleonic degrees of freedom., Comment: 7 page 3 figure

We exploit the great potential offered by Bayesian Neural Networks (BNNs) to directly decipher the internal composition of neutron stars (NSs) based on their macroscopic properties. By analyzing a set of simulated observations, namely NS radius and tidal deformability, we leverage BNNs as effective tools for inferring the proton fraction and sound speed within NS interiors. To achieve this, several BNNs models were developed upon a dataset of $\sim$ 25K nuclear EoS within a relativistic mean-field framework, obtained through Bayesian inference that adheres to minimal low-density constraints. Unlike conventional neural networks, BNNs possess an exceptional quality: they provide a prediction uncertainty measure. To simulate the inherent imperfections present in real-world observations, we have generated four distinct training and testing datasets that replicate specific observational uncertainties. Our initial results demonstrate that BNNs successfully recover the composition with reasonable levels of uncertainty. Furthermore, using mock data prepared with the DD2, a different class of relativistic mean-field model utilized during training, the BNN model effectively retrieves the proton fraction and speed of sound for neutron star matter., Comment: 16 pages, 15 figures, published version

Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models which deliver no information on composition. In this paper we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current measurements from NICER and gravitational wave observations constrain primarily the symmetric nuclear matter EOS. We then explore what could be delivered by measurements of mass and radius at the level anticipated for future large-area X-ray timing telescopes. These should be able to place very strong limits on the symmetric nuclear matter EOS, in addition to constraining the nuclear symmetry energy that determines the proton fraction inside the neutron star., Comment: Accepted for publication in MNRAS

This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided.

The high-density behavior of nuclear matter is analyzed within a relativistic mean-field description with non-linear meson interactions. To assess the model parameters and their output, a Bayesian inference technique is used. The Bayesian setup is limited only by a few nuclear saturation properties, the neutron star maximum mass larger than 2 M$_\odot$, and the low-density pure neutron matter equation of state (EOS) produced by an accurate N$^3$LO calculation in chiral effective field theory. Depending on the strength of the non-linear scalar vector field contribution, we have found three distinct classes of EOSs, each one correlated to different star properties distributions. If the non-linear vector field contribution is absent, the gravitational maximum mass and the sound velocity at high densities are the greatest. However, it also gives the smallest speed of sound at densities below three times saturation density. On the other hand, models with the strongest non-linear vector field contribution predict the largest radii and tidal deformabilities for 1.4 M$_\odot$ stars, together with the smallest mass for the onset of the nucleonic direct Urca processes and the smallest central baryonic densities for the maximum mass configuration. These models have the largest speed of sound below three times saturation density, but the smallest at high densities, in particular, above four times saturation density the speed of sound decreases approaching approximately $\sqrt{0.4}c$ at the center of the maximum mass star. On the contrary, a weak non-linear vector contribution gives a monotonically increasing speed of sound. A 2.75 M$_\odot$ NS maximum mass was obtained in the tail of the posterior with a weak non-linear vector field interaction. We found that pQCD favors models with a large non-linear vector field contribution or hyperons., Comment: There are 17 pages, 16 figures (The typographical error in the supplemental update)

The recently developed resummation technique known as {\it renormalization group optimized perturbation theory} (RGOPT) is employed in the evaluation of the EoS describing non-strange cold quark matter at NLO. Inspired by recent investigations, which suggest that stable quark matter can be made only of up and down quarks, the mass-radius relation for two flavor pure quark stars is evaluated and compared with the predictions from perturbative QCD (pQCD) at NNLO. This comparison explicitly shows that by being imbued with renormalization group properties, and a variational optimization procedure, the method allows for an efficient resummation of the perturbative series. Remarkably, when the renormalization scale is chosen so as to reproduce maximum mass stars with M=$2-2.3M_\odot$, one obtains a mass-radius curve compatible with the masses and radii of the pulsars PSR J0740+6620, PSR J0030+0451, and the compact object HESS J1731-347. Moreover, the scale dependence of the EoS (and mass-radius relation) obtained with the RGOPT is greatly improved when compared to that of pQCD. This seminal application to the description of quark stars shows that the RGOPT represents a robust alternative to pQCD when describing compressed quark matter., Comment: 13 pages, 5 figures, 2 tables

The possible existence of hybrid stars is studied using several multi-quark interaction channels. The hadronic phase consists of an equation of state (EoS) with presently accepted nuclear matter properties and the quark model is constrained by the vacuum properties of several light mesons. The dependence of several NS properties on the different quark interactions is analyzed. We show that the present constraints from neutron star observations allow for the existence of hybrid stars with large strangeness content and large quark cores., Comment: 6 pages, 5 figures; Proceedings of the XVth Quark Confinement and the Hadron Spectrum conference, August 1st - 6th, 2022, University of Stavanger, Norway

The use of the Audio Spectrogram Transformer (AST) model for gravitational-wave data analysis is investigated. The AST machine-learning model is a convolution-free classifier that captures long-range global dependencies through a purely attention-based mechanism. In this paper a model is applied to a simulated dataset of inspiral gravitational wave signals from binary neutron star coalescences, built from five distinct, cold equations of state (EOS) of nuclear matter. From the analysis of the mass dependence of the tidal deformability parameter for each EOS class it is shown that the AST model achieves a promising performance in correctly classifying the EOS purely from the gravitational wave signals, especially when the component masses of the binary system are in the range $[1,1.5]M_{\odot}$. Furthermore, the generalization ability of the model is investigated by using gravitational-wave signals from a new EOS not used during the training of the model, achieving fairly satisfactory results. Overall, the results, obtained using the simplified setup of noise-free waveforms, show that the AST model, once trained, might allow for the instantaneous inference of the cold nuclear matter EOS directly from the inspiral gravitational-wave signals produced in binary neutron star coalescences., Comment: 11 pages, 11 figures

We propose a model of asymmetric bosonic dark matter (DM) with self-repulsion mediated by the vector field coupled to the complex scalar particles. By adopting the two-fluid formalism, we study different DM distribution regimes, either, fully condensed inside the core of a star or, otherwise, distributed in a dilute halo around a neutron star (NS). We show that DM condensed in a core leads to a decrease of the total gravitational mass, radius and tidal deformability compared to a pure baryonic star with the same central density, which we will perceive as an effective softening of the equation of state (EoS). On the other hand, the presence of a DM halo increases the tidal deformability and total gravitational mass. As a result, an accumulated DM inside compact stars could mimic an apparent stiffening of strongly interacting matter equation of state and constraints we impose on it at high densities. From the performed analysis of the effect of DM particles in a MeV-GeV mass-scale, interaction strength, and relative DM fractions inside NSs we obtained a rigorous constraint on model parameters. Finally, we discuss several smoking guns of the presence of DM that are free from the above mentioned apparent modification of the strongly interacting matter equation of state. With this we could be probed with the future astrophysical and gravitational wave (GW) surveys., Comment: 14 pages, 5 figures

The extraction of the nuclear matter properties from neutron star (NS) observations is nowadays an important issue, in particular, the properties that characterize the symmetry energy which are essential to describe correctly asymmetric nuclear matter. We use deep neural networks (DNNs) to map the relation between cold $\beta$-equilibrium NS matter and the nuclear matter properties. Assuming a quadratic dependence on the isospin asymmetry for the energy per particle of homogeneous nuclear matter and using a Taylor expansion up to fourth order in the iso-scalar and iso-vector contributions, we generate a dataset of different realizations of $\beta$-equilibrium NS matter and the corresponding nuclear matter properties. The DNN model was successfully trained, attaining great accuracy in the test set. Finally, a real case scenario was used to test the DNN model, where a set of 33 nuclear models, obtained within a relativistic mean field approach or a Skyrme force description, were fed into the DNN model and the corresponding nuclear matter parameters recovered with considerable accuracy, in particular, the standard deviations $\sigma(L_{\text{sym}})= 12.85$ MeV and $\sigma(K_{\text{sat}})= 41.02$ MeV were obtained, respectively, for the slope of the symmetry energy and the nuclear matter incompressibility at saturation., Comment: 10 pages, 5 figures; published version

Equations of state for a cold neutron star's interior are presented in three-column tables that relate the baryonic density, the energy density, and the pressure. A few analytical expressions for those tables have been established these past two decades, as a convenient way to present a large number of nuclear models for neutron star matter. Some of those analytical representations are based on nonunified equations of state, in the sense that the high and the low density part of the star are not computed with the same nuclear model. Fits of equations of state based on a piecewise polytropic representation are revised by using unified tables of equations of state, that is to say models which have been calculated consistently for the core and the crust. A set of 52 unified equations of state is chosen. Each one is divided in seven polytropes via an adaptive segmentation, and two parameters per polytrope are fitted to the tabulated equation of state. The total mass, radius, tidal deformability and moment of inertia of neutron stars are modelled from the fits and compared with the quantities calculated from the original tables to ensure the accuracy of the fits on macroscopic parameters. We provide the polytropes parameters for 15 nucleonic relativistic mean field models, seven hyperonic relativistic mean field models, five hybrid relativistic mean-field models, 24 nucleonic Skyrme models, and one ab initio model. The fit error on the macroscopic parameters of neutron stars is small and well within the estimated measurement accuracy from current and next generation telescopes., Comment: 28 pages, 11 figures, 7 tables and 3 ancillary materials (tables in ASCII format and routine in Python). Accepted for publication with PRC

We demonstrate the application of the YOLOv5 model, a general purpose convolution-based single-shot object detection model, in the task of detecting binary neutron star (BNS) coalescence events from gravitational-wave data of current generation interferometer detectors. We also present a thorough explanation of the synthetic data generation and preparation tasks based on approximant waveform models used for the model training, validation and testing steps. Using this approach, we achieve mean average precision ($\text{mAP}_{[0.50]}$) values of 0.945 for a single class validation dataset and as high as 0.978 for test datasets. Moreover, the trained model is successful in identifying the GW170817 event in the LIGO H1 detector data. The identification of this event is also possible for the LIGO L1 detector data with an additional pre-processing step, without the need of removing the large glitch in the final stages of the inspiral. The detection of the GW190425 event is less successful, which attests to performance degradation with the signal-to-noise ratio. Our study indicates that the YOLOv5 model is an interesting approach for first-stage detection alarm pipelines and, when integrated in more complex pipelines, for real-time inference of physical source parameters., Comment: 11 pages, 9 figures

An ambitious goal of the astrophysical community is not only to constrain the equation of state (EOS) of neutron star (NS) matter by confronting it with astrophysics observations, but ultimately also to infer the NS composition. Nevertheless, the composition of the NS core is likely to remain uncertain unless we have an accurate determination of the nuclear symmetry energy at supra saturation density ($\rho>\rho_0$). We investigate how the nucleonic direct Urca (dUrca) processes can be used as an effective probe to constraint the high density nuclear symmetry energy. A large number of minimally constrained EOSs has been constructed by applying a Bayesian approach to study the correlations of the symmetry energy at different densities with a few selected properties of a NS. The nuclear symmetry energy above the baryon density 0.5 fm$^{-3}$ ($\sim 3 \rho_0$) is found to be strongly correlated with NS mass at which the onset of nucleonic dUrca neutrino cooling takes place in the core. This allows us to constrain the high density behavior of nuclear symmetry energy within narrow bounds. {The pure neutron matter pressure constraint from chiral effective field theory rules out the onset of nucleonic dUrca in stars with a mass $\lesssim$ 1.4 $M_\odot$.} The onset of dUrca inside 1.6 M$_\odot$ to 1.8 M$_\odot$ NS implies a slope of the symmetry energy $L$ at $\sim 2.5~\rho_0$, respectively, between 54 and 48 MeV., Comment: 6 pages, 3 figures and 1 supplemental material

By applying a relativistic mean-field description of neutron star matter with density dependent couplings, we analyse the properties of two different matter compositions: nucleonic matter with delta baryons and nucleonic matter with hyperons and delta baryons. The delta-meson couplings are allowed to vary within a wide range of values obtained by experimental data, while the hyperon-meson couplings are fitted to hypernuclear properties. Neutron star properties with no deconfinement phase transition are studied. It is verified that many models are excluded because the effective nucleon mass becomes zero before the maximum mass configuration is attained. Hyperon-free with delta-dominated composition compact stars are possible, the deltic stars. It is found that with a convenient choice of parameters the existence of deltic stars with 80% of delta baryons at the center of the star is possible. However, the presence of hyperons lowers the delta baryon fraction to values below 20% at the center and below 30% at 2-3 saturation densities. It is discussed that in the presence of delta baryons, the hyperon softening is not so drastic because deltas couple more strongly to the $\omega$-meson, and the stiffness of the equation of state is determined by the $\omega$-dominance at high densities. The speed of sound reflects very well this behavior. The compactness of the pulsar RX J0720.4-3125 imposes $x_{\sigma\Delta}>x_{\omega\Delta}>1$ and favors $x_{\rho\Delta}>1$., Comment: 14 pages, 12 figures

The possible signatures of the presence of hyperons inside neutron stars are discussed within a Bayesian inference framework applied to a set of models based on a density-dependent relativistic mean-field description of hadronic matter. Nuclear matter properties, hypernuclei properties, and observational information are used to constrain the models. General properties of neutron stars such as the maximum mass, radius, tidal deformability, proton fraction, hyperon fraction, and speed of sound are discussed. It is shown that the two solar mass constraint imposes that neutron stars described by equations of state that include hyperons have on average a larger radius, $\gtrsim 0.5$km, and a larger tidal deformability, $\gtrsim 150$, than the stars determined from a nucleonic equation of state, while the speed of sound at the center of the star is more than 25\% smaller. If a 1.4 M$_\odot$ star with a radius $\lesssim 12.5$ km is measured it is quite improbable that a massive star described by the same model contains hyperons. A similar conclusion is drawn if a two solar mass star with a radius $\lesssim 11.5$ km or a neutron star with a mass above 2.2 M$_\odot$ is observed: the possible hyperon content of these stars is ruled out or very reduced. The hyperon presence inside neutron stars is compatible with the present NICER mass-radius observations of the pulsars PSR J0030+0451 and PSR J0740+6620 and the gravitational wave detection GW170817. {It is shown that if the polytropic index $\gamma=\partial\ln p/\partial\ln\epsilon$ takes values of the order of 1.75 at not too large densities, it may indicate the onset of some kind of exotic matter, but not necessarily of deconfined quark matter., Comment: Two new fig. were added. Accepted in Phys. Rev. D

The production of light nuclei and hypernuclei together with heavy baryons, both hyperons and $\Delta$-baryons, in low density matter as found in stellar environments such as supernova or binary mergers is studied within relativistic mean-field models. Five light nuclei were considered together with three light hypernuclei. The presence of both hyperons and $\Delta$-baryons shift the dissolution of clusters to larger densities and increase the abundance of clusters. This effect is larger the smaller the charge fraction and the higher the temperature. The couplings of the $\Delta$-baryons were chosen imposing that the nucleon effective mass remains finite inside neutron stars., Comment: 11 pages, 8 figures, submitted to Phys. Rev. C

In the present work, we consider nuclear matter in the innermost crust of neutron stars under the presence of a strong magnetic field within the framework of a relativistic mean-field description. Two models with a different slope of the symmetry energy are considered in order to discuss the density-dependence of the equation of state on the crust structure. The non-homogeneous matter in $\beta$-equilibrium is described within the coexisting phases method, and the effect of including the anomalous magnetic moment is discussed. Five different geometries for the pasta structures are considered. It is shown that strong magnetic fields cause an extension of the inner crust of the neutron stars, with the occurrence of a series of disconnected non-homogeneous matter regions above the one existing for a null magnetic field. Moreover, we observed that in these disconnected regions, for some values of the magnetic field, all five different cluster geometrical shapes occur, and the gas density is close to the cluster density. Also, the pressure at the neutron star crust-core transition much larger than the pressure obtained for a zero magnetic field. Another noticeable effect of the presence of strong magnetic fields is the increase of the proton fraction, favoring the appearance of protons in the gas background., Comment: 13 pages, 8 figures. Accepted for publication in PRD

Nuclear matter properties based on a relativistic approach suitable for the description of multi-component systems are calculated. We use a set of nuclear relativistic mean-field models that satisfy acceptable nuclear matter properties and neutron star observations. The effects of the density dependence of the symmetry energy and of the Landau quantization due to the presence of a strong external magnetic field are discussed. Properties such as the proton fraction, the Landau mass, Landau parameters and entrainment matrix, the adiabatic index and speed of sound are calculated for cold $\beta$-equilibrium matter. A large dispersion on the calculated properties is obtained at two to three times saturation density $\rho_0 $. The proton Landau mass can be as low as one third of the vacuum nucleon mass at 2-3$~\rho_0 $. Similar effects are obtained for the Landau parameters, in particular, the ones involving protons, where the relative dispersion of $F^0_{pp}$ and $F^1_{pp}$ is as high as 30\% to 50\% at 2-3$~\rho_0 $. These parameters are particularly sensitive to the symmetry energy. The effect of the magnetic field on the nuclear properties is small for fields as high as 10$^{18}$G except for a small range of densities just above the crust-core transition. Tables with the EoS, and the parameters, are provided in the Supplementary Material section., Comment: 26 pages, 10 figures

The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NS), is studied within a Bayesian approach applied to a set of models based on a density dependent relativistic mean field description of nuclear matter. The EOS is subjected to a minimal number of constraints based on nuclear saturation properties and the low density pure neutron matter EOS obtained from a precise next-to-next-to-next-to-leading order (N$^{3}$LO) calculation in chiral effective field theory ($\chi$EFT). The posterior distributions of the model parameters obtained under these minimal constraints are employed to construct the distributions of various nuclear matter properties and NS properties such as radii, tidal deformabilites, central energy densities and speeds of sound etc. We found that 90% confidence interval (CI) for allowed NS mass - radius relationship and tidal deformabilites are compatible with GW170817 and recent NICER observations, without invoking the exotic degrees of freedom. A central speed-of-sound of the order of $\sqrt{2/3}$ $c$ is obtained. The maximum neutron star mass allowed by the model is 2.5 $M_\odot$., Comment: 20 pages, 11 figures (Accepted in ApJ)

We analyze the uncertainties introduced in the determination of the neutron star matter proton fraction, in a range of densities close to the saturation density, if the cold $\beta$-equilibrium neutron star matter equation of state (EoS) is known. In particular, we discuss the effect of neglecting the muon contribution and of considering that the energy density of nuclear matter is well described by taking only terms until second order in the proton-neutron asymmetry. It is shown that two types of uncertainties may be associated with the extraction of the symmetry energy from the $\beta$-equilibrium equation of state: an overestimation if terms above the parabolic approximation on the asymmetry parameter are neglected, or an underestimation if the muon contribution is not considered. The effect of the uncertainty on the symmetric nuclear matter EoS on the determination of the proton fraction is discussed. It could be shown that the neutron star mass-radius curve is sensitive to the parabolic approximation on the asymmetry parameter., Comment: 13 pages, 8 figures, published version

We determine, within a meta-model, the properties of the nuclear matter equation of state (EoS) that allow for a phase transition to deconfinement matter. It is shown that the properties that define the isoscalar channel are the ones that are affected, in particular, a phase transition implies much larger values of the skewness and kurtosis. The effect of multi-quark interaction channels in the description of the quark phase in hybrid stars is also studied. NS properties, such as the mass and radius of the quark core, show an interplay dependence between the 8-quark vector and the 4-quark isovector-vector interactions. We show that low mass NS, $M\sim 1.4 M_\odot$, may already contain a quark core, and satisfy all existing NS observational constraints. We discuss the strangeness content of the quark core and its influence on the speed of sound., Comment: 8 pages, 5 figures, Contribution to the Virtual Tribute to Quark Confinement and the Hadron Spectrum 2021 conference proceedings (vConf21)

The low density nuclear matter equation of state is strongly constrained by nuclear properties, however, for constraining the high density equation of state it is necessary to resort to indirect information obtained from the observation of neutron stars, compact objects that may have a central density several times nuclear matter saturation density, $n_0$. Taking a meta-modelling approach to generate a huge set of equation of state that satisfy nuclear matter properties close to $n_0$ and that do not contain a first order phase transition, the possibility of constraining the high density equation of state was investigated. The entire information obtained from the GW170817 event for the probability distribution of $\tilde{\Lambda}$ was used to make a probabilistic inference of the EOS, which goes beyond the constraints imposed by nuclear matter properties. Nuclear matter properties close to saturation, below $2n_0$, do not allow us to distinguish between equations of state that predict different neutron star (NS) maximum masses. This is, however, not true if the equation of state is constrained at low densities by the tidal deformability of the NS merger associated to GW170817. Above $3n_0$, differences may be large, for both approaches, and, in particular, the pressure and speed of sound of the sets studied do not overlap, showing that the knowledge of the NS maximum mass may give important information on the high density EOS. Narrowing the maximum mass uncertainty interval will have a sizeable effect on constraining the high density EOS., Comment: 11 pages, 6 figures

The abundance of light nuclei and hyperons, that are produced in stellar environments such as supernova or binary mergers, is calculated within a relativistic mean-field model with density dependent couplings in low-density matter. Five light nuclei are considered, together with three light hyper-nuclei. We show that the presence of hyperons shifts the dissolution of clusters to larger densities, and increases the amount of clusters. This effect is larger the smaller the charge fraction, and the higher the temperature. The abundance of hyperons is also affected by the cluster formation: neutral and positively charged hyperons suffer a reduction, and the negatively charged ones an increase. We also observe that the dissolution of the less-abundant clusters occurs at larger densities due to smaller Pauli-blocking effects. Overall, hyper-nuclei set in at temperatures above 25 MeV, and depending on the temperature and chemical composition, they may be more abundant than $\alpha$-particles, or even more abundant than other heavier clusters., Comment: 12 pages, 8 figures