RESEARCH

Nuclear physics

Microscopic modeling of compact stars

Since the recent (2017) detection of gravitational waves (GW) during the coalescence of two neutron stars (NS) leading to the emergence of  multi-messenger astronomy, it is clear that nuclear physics has a key role to play in the understanding of these spectacular astrophysical phenomena, by providing microscopic theoretical and experimental data that can be directly connected with the observations. 

The LPC theory team has been involved for several years in the microscopic modeling of dense matter as found in all astrophysical phenomena concerning compact stars, namely neutron stars and proto-neutron stars (pulsars, magnetars, gamma and X-ray sources), core-collapse supernovae and neutron star binaries, with various theoretical methods (variational approaches, extended Thomas-Fermi theory, relativistic mean-field, compressible liquid droplet model, statistical models). More specifically, our expertise concerns the theoretical modeling of the equation of state (EoS), composition and phase transitions of baryonic matter under conditions of density and temperature where it can be described by nucleonic or hyperonic degrees of freedom.  In the context of the current discoveries, our research aims (i) on the one hand to establish a direct and quantitative connection between microphysical data and astrophysical observables using advanced statistical techniques, (ii) and on the other hand to develop original theoretical tools to describe the EoS at finite temperature, a subject that is still treated with strong approximations.

Degenerate regime

We have developed an EoS meta-modeling technique [3] that allows to extrapolate with Bayesian techniques the uncertainties on microscopic data (experimental or theoretical) into controlled predictions of astrophysical observables, such as masses and radii of compact stars, tidal polarizability measurable by the GW signal, moment of inertia of the crust, amplitude of pulsar glitches.

We have extended the formalism initially valid only for homogeneous matter in the liquid phase, to the solid phase of the crust and the associated inhomogeneous matter [7]. This allowed us to have a unified treatment of the EoS which we showed to be essential for a correct evaluation of the NS radius [1] as can be estimated by tidal polarizability measurements, and to show for the first time the key role of the surface energy of very neutron-rich nuclei in the determination of the crustal properties.

Figure tirée de DOI: 10.1103/PhysRevC.100.055803

Estimation bayesienne du moment d’inertie (en ordonnées) et de l’épaisseur (en abscisses) d’une étoile à neutron de masse M=1.4Mo, avec prise en compte des contraintes provenant des données expérimentales en physique nucléaire, des calculs ab-initio, et des contraintes observationnelles. Les symboles donnent les prédictions de quelques modèles populaires d’équation d’état.

Bayesian estimation of the moment of inertia (ordinate) and thickness (abscissa) of a neutron star of mass M=1.4Mo, taking into account constraints from experimental nuclear physics data, ab-initio calculations, and observational constraints. The symbols give the predictions of some popular equation of state models.

Composition de la croute d’une étoile à neutrons à différentes températures, où Tm représente la température de cristallization, obtenue avec l’équation d’état unifiée BSK24.

Composition of the crust of a neutron star at different temperatures, where Tm is the crystallization temperature, obtained with the BSK24 unified equation of state.

The obvious continuation of this project concerns the proposal of waveforms for the analysis of GW signals incorporating the current microscopic knowledge of the EoS. For this purpose, the team is part of the VIRGO/LIGO collaboration since December 2020.

Finite temperature

The additional difficulty of finite temperature modeling is that, if the density is below the nuclear saturation density, the matter is organized into clusters of nucleons and the statistical distribution of microstates must be explicitly considered. We have produced [5] the only equation-of-state table including the complete cluster distributions and based on modern non-relativistic functionals that is publicly available for direct use in core collapse calculations. This table is available on the CompOse archive and partially allows to test the effect of microphysics on the collapse, by comparison with relativistic models computed by other teams and available on the same archive.

In order to be able to test separately the different physical ingredients, we have developed [4] an original perturbative theoretical approach that allows to compute the cluster distribution from any standard EoS model, based on the single nucleus approximation.   Thanks to this approach, we have been able to deconvolute the effect of the different ingredients, and to show that the key variable is given by the electron capture rates on very neutron rich nuclei [2]. Very recently [8], we have extended this approach to the study of crustal crystallization and produced the first self-consistent microscopic calculation of the impurity factor that governs the neutrino scattering rates in proto-NS.

Distribution de noyaux normalisée en fonction de la densité dans la croûte interne, calculée à la température de crystallization.

Distribution of nuclei at different densities and temperatures. The simplified estimate of the single nucleus approximation is given by the vertical lines.

Distribution de noyaux à différentes densités et température. L’estimation simplifiée du noyau unique est donnée par les lignes verticales.
Distribution de noyaux à différentes densités et température. L’estimation simplifiée du noyau unique est donnée par les lignes verticales.

Finally we started to tackle the problem of determining the mass change of light clusters in a dense medium [6] and suggested a novel method of analyzing nuclear multi-fragmentation data to be able to constrain the parameters of the theory.

Read more

Scroll to top