Time-dependent nuclear energy-density functional theory toolkit for neutron star crust: Dynamics of a nucleus in a neutron superfluid

Daniel Pȩcak,Agata Zdanowicz,Nicolas Chamel,Piotr Magierski,Gabriel Wlazłowski

2024-03-26

Abstract:We present a new numerical tool designed to probe the dense layers of neutron star crusts. It is based on the Time-Dependent Hartree-Fock-Bogoliubov theory with generalized Skyrme nuclear energy density functionals, such as the Brussels-Montreal ones. We use it to study the time evolution of a nucleus accelerating through superfluid neutron medium in the inner crust of a neutron star. We extract an effective mass in the low velocity limit. We observe a threshold velocity and specify mechanisms of dissipation: phonon emission, Cooper pairs breaking, and vortex rings creation. The microscopic effects we study have impact on neutron star. Moreover, the mechanisms, we described, are general and apply also to other fermionic superfluid mixtures like liquid helium, or ultracold gases.

Nuclear Theory,High Energy Astrophysical Phenomena,Superconductivity

What problem does this paper attempt to address?

The main focus of this paper is on the dense inner layers of neutron stars, particularly the dynamic properties of the neutron star crust. The research team developed a new numerical tool based on the time-dependent Hartree-Fock-Bogoliubov theory (TDHFB), utilizing the Skyrme nuclear energy density functional, specifically the Brussels-Montreal series. They used this tool to study the time evolution of nuclear clusters (referred to as impurities) accelerating in neutron superfluids, extracted the effective mass at low velocities, and analyzed threshold speeds and dissipation mechanisms such as phonon emission, Cooper pair breaking, and vortex formation. These microphysical effects are crucial for understanding the properties of neutron stars and are also applicable to other fermionic superfluid mixtures.

The background of the study revolves around neutron stars as remnants of stellar explosions, which have extremely high internal densities and may exhibit various quantum phase transitions, including the formation of neutron superfluidity through 1S0 pairing. As experimental conditions cannot replicate the environment inside neutron stars, theoretical models and indirect astronomical observations are needed to explore their properties. Density functional theory (DFT) becomes a powerful tool in this context as it can handle both electron and nuclear systems, although with lower accuracy in the case of nuclear systems. Nonetheless, with improved functional quality, it remains an important method for describing nuclear properties.

The W-BSk toolkit introduced in the paper is a versatile software for three-dimensional simulations of dynamic processes within the neutron star core. It can handle static and dynamic phenomena at zero and finite temperatures, without restrictions on symmetry. Leveraging high-performance computing (HPC) capabilities, it has the ability to simulate regions exceeding (100 fm)^3 in volume, encompassing thousands of neutrons and protons. The study uses a nuclear cluster in motion within a neutron superfluid as an example to investigate its effective mass and dissipation channels, providing new insights into the complex dynamics within neutron stars.