Eur. Phys. J. A (2019) 55: 188
https://doi.org/10.1140/epja/i2019-12884-1

Regular Article - Theoretical Physics

Bayesian analysis of the crust-core transition with a compressible liquid-drop model

¹ Laboratoire de Physique Corpusculaire, CNRS, ENSICAEN, UMR6534, Université de Caen Normandie, F-14050, Caen Cedex, France
² Institut de Physique Nucléaire de Lyon, CNRS/IN2P3, Université de Lyon, Université Claude Bernard Lyon 1, F-69622, Villeurbanne Cedex, France

^* e-mail: carreau@lpccaen.in2p3.fr

Received: 25 February 2019
Accepted: 17 September 2019
Published online: 30 October 2019

Abstract

The crust-core (CC) phase transition of neutron stars is studied within a unified meta-modeling of the nuclear Equation of State (EoS). The variational equations in the crust are solved within a Compressible Liquid-Drop (CLD) approach, with surface parameters consistently optimized for each EoS set on experimental nuclear mass data. When EoS parameters are taken from known Skyrme or RMF functionals, the transition point of those models is nicely reproduced as expected. The probability distribution of EoS parameters and of the transition density and pressure is determined with a Bayesian analysis, where the prior is given by an uncorrelated distribution of parameters considering the present empirical uncertainties, and constraints are applied both from neutron star physics and ab initio predictions. We show that the properties of the CC transition point are largely independent of the high density properties of the EoS, while ab initio EoS calculations of neutron and symmetric matter are far more constraining. The surface tension of extremely neutron rich matter, which remains poorly known, is the most influential parameter for the CC transition point. This explains the large dispersion of existing predictions of the CC transition point. Fixing the isospin dependence of the surface tension to a reasonable but somewhat arbitrary value, strong correlations with isovector EoS parameters (, and ) are recovered. Within the present experimental and theoretical uncertainties on those parameters, we estimate the transition density as fm^-3 and the transition pressure as MeV fm^-3.