A new statistical method for the structure of the inner crust of neutron stars

Alessandro Pastore*, M. Shelley, Simone Baroni, C. A. Diget

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review


We investigated the structure of the low density regions of the inner crust of neutron stars using the Hartree-Fock-Bogoliubov (HFB) model to predict the proton content Z of the nuclear clusters and, together with the lattice spacing, the proton content of the crust as a function of the total baryonic density . The exploration of the energy surface in the configuration space and the search for the local minima require thousands of calculations. Each of them implies a HFB calculation in a box with a large number of particles, thus making the whole process very demanding. In this work, we apply a statistical model called Gaussian process emulator that makes the exploration of the energy surface ten times faster. We also present a novel treatment of the HFB equations that leads to an uncertainty on the total energy of per particle. Such a high precision is necessary to distinguish neighbour configurations around the energy minima.

Original languageEnglish
Article number094003
JournalJournal of physics g-Nuclear and particle physics
Issue number9
Publication statusPublished - 11 Aug 2017

Bibliographical note

Invited article for special issue "emerging leaders".
This is an author-produced version of the published paper. Uploaded in accordance with the publisher’s self-archiving policy. Further copying may not be permitted; contact the publisher for details


  • Gaussian process emulator
  • Hartree-Fock-Bogoliubov
  • inner crust
  • neutron stars

Cite this