$^1$S$_0$ pairing gaps, chemical potential and entrainment matrix in superfluid neutron-star cores for the Brussels-Montreal functionals

Imagine the universe's most extreme playground: a neutron star. These are the collapsed, super-dense cores of dead stars, so heavy that a single teaspoon of their material would weigh a billion tons on Earth. Deep inside, the pressure is so immense that atoms are crushed flat, leaving behind a soup of subatomic particles. But this isn't just a static, cold lump; it's a dynamic, churning fluid. In fact, the core is expected to be a "superfluid," a state of matter where particles flow with zero friction, like a liquid that never stops moving.

To understand how these stars spin, wobble, or even glitch (suddenly speed up), scientists need to know how the different particles inside interact. Think of the core as a crowded dance floor with two main groups: neutrons and protons. In a superfluid, these particles pair up and dance in perfect sync. However, they don't dance alone; they influence each other. If the neutrons start to spin, they drag the protons along with them, and vice versa. This "drag" is called entrainment. It's like two skaters holding hands; if one tries to turn, the other has to turn too, even if they didn't want to. The paper you are about to read dives deep into the math that describes exactly how strong this grip is, how temperature changes the dance, and how fast the particles can move before the pairing breaks apart.

---

The Cosmic Dance Floor: A New Map for Neutron Stars

Neutron stars are nature's ultimate laboratories, but they are also incredibly tricky to model. The core of a neutron star is a dense mixture of neutrons and protons (along with some electrons and muons) that behave like a superfluid. In this state, the particles form pairs and flow without resistance. But here's the catch: the neutrons and protons don't just flow independently. Because of the strong nuclear forces between them, they are "entrained." This means if you try to push the neutrons, the protons get dragged along, and if you push the protons, the neutrons follow. It's a cosmic version of a tangle of wet spaghetti where moving one strand pulls the whole mess with it.

For a long time, scientists have tried to map out exactly how this dragging works. They needed to know three main things: how much energy it takes to break the particle pairs (the pairing gap), how much "push" is needed to keep them moving (the chemical potential), and exactly how tightly the neutrons and protons are linked (the entrainment matrix). The problem is that these values change depending on how hot the star is and how fast the particles are flowing. Previous models often had to make guesses or use simplified rules that didn't account for all these changing conditions at once.

This paper, written by Valentin Allard and Nicolas Chamel, takes a giant leap forward by solving the full, complicated math equations (known as the time-dependent Hartree–Fock–Bogoliubov equations) without cutting corners. They used a specific set of rules for how nuclear matter behaves, called the Brussels–Montreal functional BSk24, which is known to be very accurate. They simulated the entire range of temperatures and speeds where superfluidity can exist in the outer core of a neutron star.

What They Found: The Rules of the Dance

The team's main discovery is a complete, self-consistent map of how neutrons and protons behave in this extreme environment. They didn't just guess; they calculated the exact numbers for the pairing gaps, chemical potentials, and the entrainment matrix across the whole range of conditions.

One of their most interesting findings concerns the pairing gaps. This is the energy "glue" holding the particle pairs together. They found that for protons, this glue is surprisingly weak—much weaker than what many other studies have assumed. In fact, the proton pairing gaps they calculated are significantly smaller than the neutron ones. This matters because if the glue is weak, the protons might stop superfluid behavior at lower temperatures or speeds than previously thought. Their results align with recent, very complex calculations that consider how particles interact with each other in the medium, suggesting that the "weak glue" for protons is likely the real deal.

They also discovered that the behavior of these particles is surprisingly universal. When you look at the data using the right "rulers" (normalizing the temperature and speed by their critical limits), the behavior of neutrons and protons looks almost identical. Whether it's a neutron or a proton, the way the pairing gap shrinks as the star gets hotter or the flow gets faster follows the same pattern. This suggests that the underlying physics is the same for both, governed by the same weak-coupling rules.

Another key result involves the entrainment matrix, which tells us how much the neutrons drag the protons. The authors found that as long as the flow speed is below a certain critical limit (Landau's critical velocity), the dragging effect doesn't change much, even if the particles are moving fast. However, once the flow gets too fast or the temperature gets too high, the superfluidity starts to break down. At that point, the neutrons and protons stop dragging each other and become "dynamically uncoupled." The neutrons might stop being superfluid entirely, while the protons keep dancing on their own, no longer tethered to the neutrons.

What They Ruled Out and What's Still Unknown

The paper is very clear about what it doesn't cover. The authors explicitly state that they are only looking at the outer core of the neutron star, where the density is high but not extreme enough to create exotic particles like hyperons or quarks. They also focused only on the 1S0 pairing phase (a specific way neutrons and protons pair up). They ruled out the 3PF2 neutron superfluidity, noting that in the regions where both could exist, the 1S0 phase completely pushes the 3PF2 phase out of the picture unless there are incredibly strong magnetic fields present. So, for the outer core, the 3PF2 phase is effectively non-existent.

They also didn't just assume the results; they tested them against various approximations. They found that while some older, simpler models (like Landau's theory) work well for the entrainment matrix if you know the critical temperatures, they aren't enough to calculate the chemical potentials accurately. To get the true chemical potentials, you really need to solve the full, heavy-duty equations they used.

Why This Matters

This work provides the "ingredients" needed to build better models of neutron stars. Just like a chef needs accurate measurements of flour and sugar to bake a perfect cake, astrophysicists need precise numbers for pairing gaps and entrainment to understand how neutron stars spin, how they cool down, and how they react to gravitational waves. By providing these consistent, microscopic inputs based on the BSk24 functional, Allard and Chamel have given the scientific community a more reliable recipe for simulating the hearts of these mysterious, dense stars. They haven't solved every mystery of the neutron star, but they have cleared up the fog around how the two main dancers in the core's ballroom move together.