The modification of the ground state properties of light atomic nuclei in the nuclear and stellar medium is addressed, using chemical equilibrium constants evaluated from a new analysis of the intermediate energy heavy-ion (Xe$+$Sn) collision data measured by the INDRA collaboration. Three different reactions are considered, mainly differing by the isotopic content of the emission source. The thermodynamic conditions of the data samples are extracted from the measured multiplicities allowing for a parametrization of the in-medium modification, determined with the single hypothesis that the different nuclear species in a given sample correspond to a unique common value for the density of the expanding source. We show that this correction, which was not considered in previous analyses of chemical constants from heavy ion collisions, is necessary, since the observables of the analyzed systems show strong deviations from the expected results for an ideal gas of free clusters. This data set is further compared to a relativistic mean-field model, and seen to be reasonably compatible with a universal correction of the attractive $\sigma$-meson coupling.

The equation of state with light clusters for nuclear and stellar matter is determined using chemical equilibrium constants evaluated from the analysis of the recently published (Xe$+$Sn) heavy ion data, corresponding to three reactions with different isotopic contents of the emission source. The measured multiplicities are used to extract the thermodynamic properties, and an in-medium correction to the ideal gas internal partition function of the clusters is included in the analysis. This in-medium correction and its respective uncertainty are calculated via a Bayesian analysis, with the unique hypothesis that the different nuclear species in a given sample must correspond to a unique common value for the density of the expanding source. Different parameter sets for the correction are tested, and the effect of the radius of the clusters on the thermodynamics and on the chemical equilibrium constants is also addressed. It is shown that the equilibrium constants obtained are almost independent of the isospin content of the analysed systems. Finally, a comparison with a relativistic mean field model proves that data are consistent with a universal in-medium correction of the scalar $\sigma$-meson coupling for nucleons bound in clusters. The obtained value, $g_s/g_s^0 = 0.92 \pm 0.02$, is larger than that obtained in a previous study not including in-medium effects in the data analysis. This result implies a smaller effect on the binding energy of the clusters and, as a consequence, larger melting densities, and an increased cluster contribution in supernova matter.

The possible presence of amorphous and heterogeneous phases in the inner crust of a neutron star is expected to reduce the electrical conductivity of the crust, with potentially important consequences on the magneto-thermal evolution of the star. In cooling simulations, the disorder is quantified by an impurity parameter which is often taken as a free parameter. We aim to give a quantitative prediction of the impurity parameter as a function of the density in the crust,performing microscopic calculations including up-to-date microphysics of the crust. A multi-component approach is developed at finite temperature using a compressible liquid drop description of the ions with an improved energy functional based on recent microscopic nuclear models and optimized on extended Thomas-Fermi calculations. Thermodynamic consistency is ensured by adding a rearrangement term and deviations from the linear mixing rule are included in the liquid phase. The impurity parameter is consistently calculated at the crystallization temperature as determined in the one-component plasma approximation for the different functionals. Our calculations show that at the crystallization temperature the composition of the inner crust is dominated by nuclei with charge number around $Z \approx 40$, while the range of the $Z$ distribution varies from about 20 near the neutron drip to about 40 closer to the crust-core transition. This reflects on the behavior of the impurity parameter that monotonically increases with density up to around 40 in the deeper regions of the inner crust. Our study shows that the contribution of impurities is non-negligible, thus potentially having an impact on the transport properties in the neutron-star crust. The obtained values of the impurity parameter represent a lower limit; larger values are expected in the presence of non-spherical geometries and/or fast cooling dynamics.

The interior of a neutron star (NS) is usually assumed to be made of cold catalyzed matter. However, the outer layers are unlikely to remain in full equilibrium during the formation of the star and its cooling, especially after crystallization. We study the cooling and equilibrium composition of the outer layers of a NS down to crystallization. Here the impurity parameter, usually a free parameter in cooling simulations, is calculated self-consistently using a microscopic model for which a unified equation of state has recently been determined. We follow the evolution of the nuclear distributions of the multi-component Coulomb liquid plasma (MCP) fully self-consistently, adapting a general formalism originally developed for the description of supernova cores. We calculate the impurity parameter at the crystallization as determined in the one-component plasma (OCP) approximation. Our analysis shows that the sharp changes in composition obtained in the OCP approximation are smoothed out when a full nuclear distribution is allowed. The Coulomb coupling parameter at melting is found to be reasonably close to the canonical value of 175, except for specific pressures for which supercooling occurs in the OCP approximation. Our MCP treatment leads to non-monotonic variations of the impurity parameter with pressure. Its values can change by several orders of magnitude reaching about 50, suggesting that the crust may be composed of an alternation of pure (highly conductive) and impure (highly resistive) layers. The results presented here complement the recent unified equation of state obtained with the same model. Our self-consistent approach to hot MCP shows that the presence of impurities in the outer crust of a NS is non-negligible and may have a sizeable impact on transport properties. In turn, this may have important implications for the cooling of NS and their magneto-rotational evolution.

Context. In the cooling process of a non-accreting neutron star, the composition and properties of the crust are thought to be fixed at the finite temperature where nuclear reactions fall out of equilibrium. A lower estimation for this temperature is given by the crystallization temperature, which can be as high as $\approx 7\times 10^9$ K in the inner crust, potentially leading to sizeable differences with respect to the simplifying cold-catalyzed matter hypothesis. Aims. We extend the recent work by Fantina et al. (2019) on the outer crust, to the study of the crystallization of the inner crust and the associated composition in the one-component plasma approximation. Methods. The finite temperature variational equations for non-uniform matter in both the liquid and the solid phases are solved using a compressible liquid-drop approach with parameters optimized on four different microscopic models which cover the present uncertainties in nuclear modeling. Results. We consider separately the effect of the different nuclear ingredients with their associated uncertainties, namely the nuclear equation of state, the surface properties in the presence of a uniform gas of dripped neutrons, and the proton shell effects arising from the ion single-particle structure. Our results suggest that the highest source of model dependence comes from the smooth part of the nuclear functional. Conclusions. We show that shell effects play an important role at the lowest densities close to the outer crust, but the most important physical ingredient to be settled for a quantitative prediction of the inner crust properties is the surface tension at extreme isospin values.

The impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates.

Light and heavy clusters are calculated for asymmetric warm nuclear matter in a relativistic mean-field approach. In-medium effects, introduced via a universal cluster-meson coupling, and a binding energy shift contribution, calculated in a Thomas-Fermi approximation, were taken into account. This work considers, besides the standard lightest bound clusters $^4$He, $^3$He, $^3$H, and $^2$H, also stable and unstable clusters with higher number of nucleons, in the range $5\leq A \leq 12$, as it is natural that heavier clusters also form in core-collapse supernova matter, before the pasta phases set in. We show that these extra degrees of freedom contribute with non-negligible mass fractions to the composition of nuclear matter, and may prevail over deuterons and $\alpha$ particles at high density in strongly asymmetric matter, and not too high temperatures. The presence of the light clusters reduces the contribution of heavy clusters to a much smaller density range, and to a smaller mass fraction.

The neutron star crustal EoS and transition point properties are computed within a unified meta-modeling approach. A Bayesian approach is employed including two types of filters: bulk nuclear properties are controlled from low density effective field theory (EFT) predictions as well as the present knowledge from nuclear experiments, while the surface energy is adjusted on experimental nuclear masses. Considering these constraints, a quantitative prediction of crustal properties can be reached with controlled confidence intervals and increased precision with respect to previous calculations: {$\approx 11\%$} dispersion on the crustal width and {$\approx 27\%$} dispersion on the fractional moment of inertia. The crust moment of inertia is also evaluated as a function of the neutron star mass, and predictions for mass and radii are given for different pulsars. The possible crustal origin of Vela pulsar glitches is discussed within the present estimations of crustal entrainment, disfavoring a large entrainment phenomenon if the Vela mass is above $1.4M_\odot$. Further refinement of the present predictions requires a better estimation of the high order isovector empirical parameters, e.g. $K_{sym}$ and $Q_{sym}$, and a better control of the surface properties of extremely neutron rich nuclei.

The influence of the uncertainties of the equation of state empirical parameters on the neutron stars crust-core phase transition is explored within a meta-modeling approach, in which the energy per particle is expanded as a Taylor series in density and asymmetry around the saturation point. The phase transition point is estimated from the intersection of the spinodal instability region for dynamical fluctuations with the chemical equilibrium curve. Special attention is paid to the inclusion of high-order parameters of the Taylor series and their influence on the transition point. An uncorrelated prior distribution is considered for the empirical parameters, with bulk properties constrained through effective field theory predictions, while the surface parameters are controlled from a fit of nuclear masses using the extended Thomas Fermi approximation. The results show that the isovector compressibility $K_{sym}$ and skewness $Q_{sym}$ have the most significant correlations with the transition point, along with the previously observed influence of the $L_{sym}$ parameter. The estimated density and pressure of the crust-core transition are $n_t = (0.071 \pm 0.011) fm^{-3}$ and $P_t = (0.294 \pm 0.102) MeV fm^{-3}$.

The presence of strong magnetic fields in neutron stars, such as in magnetars, may significantly affect their crust-core transition properties and the crust size. This knowledge is crucial in the correct interpretation of astrophysical phenomena involving magnetars, such as glitches in observed rotation frequencies, cooling, bursts and possibly tidal polarizabilities. A recently developed meta-modelling technique allows exploring the model dependence of density functional theory equation of state calculations. In this work, we extend this meta-model to investigate the effect of strong magnetic fields on spinodal instabilities of neutron star matter and the associated crust-core properties. Both Tolman-Oppenheimer-Volkov and a full self-consistent numerical calculations are performed for the neutron star structure, the results being quantitatively different for strong magnetic fields.

A quantitative knowledge of the nuclear Equation of State (EoS) requires an accurate estimation of the uncertainties on the EoS parameters and their mutual correlations. Such correlations are empirically observed in a large set of EoS models by different authors, but they are not always fully understood. We show that some of these correlations can be interpreted from basic physical constraints imposed on a simple Taylor expansion of the binding energy around saturation density. In particular, we investigate the correlations among the following empirical parameters: the symmetry energy $E_{sym}$, the slope and curvature of the symmetry energy $L_{sym}$ and $K_{sym}$, and the curvature and skewness of the binding energy in symmetric matter $K_{sat}$ and $Q_{sat}$. The uncertainties on these correlations is estimated through analytical modelling as well as a meta-modelling analysis of the EoS subject to physical constraints. We show that a huge dispersion of the correlations among low order empirical parameters is induced by the unknown higher order empirical parameters, such as $K_{sym}$ (second order) and $Q_{sat}$ and $Q_{sym}$ (third order). We also propose an explanation of the reason why $Q_{sat}$ is weakly constrained by present experimental data. We conclude that selective observables on high order parameters, such as $K_{sym}$ and $Q_{sat}$, should be better determined before the present uncertainties on the EoS can be further reduced.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

Background : The emergence of hyperon degrees of freedom in neutron star matter has been associated to first order phase transitions in some phenomenological models, but conclusions on the possible physical existence of an instability in the strangeness sector are strongly model dependent. Purpose : The purpose of the present study is to assess whether strangeness instabilities are related to specific values of the largely unconstrained hyperon interactions, and to study the effect of the strange meson couplings on phenomenological properties of neutron stars and supernova matter, once these latter are fixed to fulfill the constraints imposed by hypernuclear data. Method : We consider a phenomenological RMF model sufficiently simple to allow a complete exploration of the parameter space. Results : We show that no instability at supersaturation density exists for the RMF model, as long as the parameter space is constrained by basic physical requirements. This is at variance with a non-relativistic functional, with a functional behavior fitted through ab-initio calculations. Once the study is extended to include the full octet, we show that the parameter space allows reasonable radii for canonical neutron stars as well as massive stars above two-solar mass, together with an important strangeness content of the order of 30\%, slightly decreasing with increasing entropy, even in the absence of a strangeness driven phase transition. Conclusions : We conclude that the hyperon content of neutron stars and supernova matter cannot be established with present constraints, and is essentially governed by the unconstrained coupling to the strange isoscalar meson.

The uncertainties in neutron star (NS) radii and crust properties due to our limited knowledge of the equation of state (EOS) are quantitatively analysed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core EOS based on models with different properties at nuclear matter saturation, the uncertainties can be as large as $\sim 30\%$ for the crust thickness and $4\%$ for the radius. Necessary conditions for causal and thermodynamically consistent matchings between the core and the crust are formulated and their consequences examined. A large set of unified EOS for purely nucleonic matter is obtained based on 24 Skyrme interactions and 9 relativistic mean-field nuclear parametrizations. In addition, for relativistic models 17 EOS including a transition to hyperonic matter at high density are presented. All these EOS have in common the property of describing a $2\;M_\odot$ star and of being causal within stable NS. A span of $\sim 3$ km and $\sim 4$ km is obtained for the radius of, respectively, $1.0\;M_\odot$ and $2.0\;M_\odot$ star. Applying a set of nine further constraints from experiment and ab-initio calculations the uncertainty is reduced to $\sim 1$ km and $2$ km, respectively. These residual uncertainties reflect lack of constraints at large densities and insufficient information on the density dependence of the EOS near the nuclear matter saturation point. The most important parameter to be constrained is shown to be the symmetry energy slope $L$ which exhibits a linear correlation with the stellar radius, particularly for masses $\sim 1.0\;M_\odot$. Potential constraints on $L$, the NS radius and the EOS from observations of thermal states of NS are also discussed. [Abriged]

Though generally agreed that the symmetry energy plays a dramatic role in determining the structure of neutron stars and the evolution of core-collapsing supernovae, little is known in what concerns its value away from normal nuclear matter density and, even more important, the correct definition of this quantity in the case of unhomogeneous matter. Indeed, nuclear matter traditionally addressed by mean-field models is uniform while clusters are known to exist in the dilute baryonic matter which constitutes the main component of compact objects outer shells. In the present work we investigate the meaning of symmetry energy in the case of clusterized systems and the sensitivity of the proto-neutron star composition and equation of state to the effective interaction. To this aim an improved Nuclear Statistical Equilibrium (NSE) model is developed, where the same effective interaction is consistently used to determine the clusters and unbound particles energy functionals in the self-consistent mean-field approximation. In the same framework, in-medium modifications to the cluster energies due to the presence of the nuclear gas are evaluated. We show that the excluded volume effect does not exhaust the in-medium effects and an extra isospin and density dependent energy shift has to be considered to consistently determine the composition of subsaturation stellar matter. The symmetry energy of diluted matter is seen to depend on the isovector properties of the effective interaction, but its behavior with density and its quantitative value are strongly modified by clusterization.

A simple, exactly solvable statistical model is presented for the description of baryonic matter in the thermodynamic conditions associated to the evolution of core-collapsing supernova. It is shown that the model presents a first order phase transition in the grandcanonical ensemble which is not observed in the canonical ensemble. Similar to other model systems studied in condensed matter physics, this ensemble in-equivalence is accompanied by negative susceptibility and discontinuities in the intensive observables conjugated to the order parameter. This peculiar behavior originates from the fact that baryonic matter is subject to attractive short range strong forces as well as repulsive long range electromagnetic interactions, partially screened by a background of electrons. As such, it is expected in any theoretical treatment of nuclear matter in the stellar environment. Consequences for the phenomenology of supernova dynamics are drawn.

Bulk fermionic matter, as it can be notably found in supernova matter and neutrons stars, is subject to correlations of infinite range due to the antisymmetrisation of the N-body wave function, which cannot be explicitly accounted for in a practical simulation. This problem is usually addressed in condensed matter physics by means of the so-called Twist Averaged Boundary Condition method. A different ansatz based on the localized Wannier representation has been proposed in the context of antisymmetrized molecular dynamics. In this paper we work out the formal relation between the two approaches. We show that, while the two coincide when working with exact eigenstates of the N-body Hamiltonian, differences appear in the case of variational approaches, which are currently used for the description of stellar matter. Some model applications with Fermionic Molecular Dynamics are shown.

We develop a phenomenological statistical model for dilute star matter at finite temperature, in which free nucleons are treated within a mean-field approximation and nuclei are considered to form a loosely interacting cluster gas. Its domain of applicability, that is baryonic densities ranging from about $\rho>10^8$ g $\cdot$ cm$^{-3}$ to normal nuclear density, temperatures between 1 and 20 MeV and proton fractions between 0.5 and 0, make it suitable for the description of baryonic matter produced in supernovae explosions and proto-neutron stars. The first finding is that, contrary to the common belief, the crust-core transition is not first order, and for all subsaturation densities matter can be viewed as a continuous fluid mixture between free nucleons and massive nuclei. As a consequence, the equations of state and the associated observables do not present any discontinuity over the whole thermodynamic range. We further investigate the nuclear matter composition over a wide range of densities and temperatures. At high density and temperature our model accounts for a much larger mass fraction bound in medium nuclei with respect to traditional approaches as Lattimer-Swesty, with sizeable consequences on the thermodynamic quantities. The equations of state agree well with the presently used EOS only at low temperatures and in the homogeneous matter phase, while important differences are present in the crust-core transition region. The correlation among the composition of baryonic matter and neutrino opacity is finally discussed, and we show that the two problems can be effectively decoupled.

In this letter we show that the presence of the long-range Coulomb force in dense stellar matter implies that the total charge cannot be associated with a chemical potential, even if it is a conserved quantity. As a further consequence, the analytical properties of the partition sum are modified, changing the order of the phase transitions and affecting the possible occurrence of critical behaviours. The peculiar thermodynamic properties of the system can be understood introducing a model hamiltonian in which each charge is independently neutralized by a uniform background of opposite charge.

Because of the presence of a liquid-gas phase transition in nuclear matter, compact-star matter can present a region of instability against the formation of clusters. We investigate this phase separation in a matter composed of neutrons, protons and electrons, within a Skyrme-Lyon mean-field approach. Matter instability and phase properties are characterized through the study of the free-energy curvature. The effect of beta-equilibrium is also analyzed in detail, and we show that the opacity to neutrinos has an influence on the presence of clusterized matter in finite-temperature proto-neutron stars.