Emergence of the geometric contribution to the superfluid density in the inner crust of neutron stars

This paper clarifies the origin of the geometric contribution to the superfluid density in neutron star inner crusts by deriving its dependence on the pairing gap within a multi-band framework and demonstrating that accounting for corrections to Bogoliubov quasi-particle states, rather than just Hartree-Fock single-particle states, is essential for its emergence.

The Gist

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Picture: A Cosmic Dance Floor

Imagine a neutron star not as a solid rock, but as a giant, cosmic dance floor. Inside the "crust" (the outer shell) of this star, there are two distinct groups of dancers:

1. The Lattice: A rigid, crystalline grid of atomic nuclei (like a fixed pattern of chairs).
2. The Superfluid: A sea of neutrons that can flow without any friction, like a ghostly, slippery fluid sliding between the chairs.

For decades, scientists have been trying to figure out exactly how much of this neutron fluid is actually "superfluid" (able to flow freely) versus how much is "stuck" to the chairs. This is crucial because the amount of free-flowing fluid determines how the star spins and why it sometimes suddenly speeds up (a phenomenon called a "pulsar glitch").

The Problem: The "Stuck" Neutrons

In the past, scientists thought the neutrons were heavily "entrained." Think of this like trying to run through a crowded room where everyone is holding hands. Even if you want to run fast, the crowd drags you down. The lattice (the chairs) was thought to drag the neutrons so much that very little of them could actually flow freely.

However, recent calculations suggested something surprising: There is a hidden "geometric" force that helps the neutrons flow more freely than we thought. This paper explains why that happens.

The Key Discovery: The "Band Theory" Dance

To understand this, we have to look at the neutrons through the lens of Band Theory. Imagine the energy levels of the neutrons as a multi-story building with many floors (bands).

• In normal materials (like copper wire), only the top floor is partially full, and the rest are empty.
• In a neutron star's crust, it's a chaotic skyscraper where many floors are cut in half by the "Fermi energy" (the current energy level of the system). It's like having a building where the water level is halfway up the 10th, 20th, and 30th floors simultaneously.

The author, Giorgio Almirante, shows that when you have this complex, multi-floor situation, a new type of contribution to the superfluid density appears. He calls it the Geometric Contribution.

The Analogy: The "Mixing" Dance Move

Here is the core of the paper, explained with a metaphor:

The Old View (Conventional Contribution):
Imagine the neutrons are dancers on a single floor. If they want to move, they just slide across the floor. The amount they can slide depends on how "sticky" the floor is. This is the "conventional" way we calculated things before.

The New View (Geometric Contribution):
Now, imagine the dancers are on a complex, multi-level stage. Because the neutrons are "superfluid," they have a special ability to mix with dancers on other floors.

• The paper explains that if you only look at the dancers on their own floor (ignoring the connection to other floors), you get the wrong answer.
• The "Geometric Contribution" comes from the fact that the superfluid state allows neutrons to borrow energy and momentum from empty floors and mix with filled floors.
• It's like a dancer who, instead of just walking, can suddenly teleport or slide onto a different level of the stage to help the whole group move faster.

The "Gap" Factor:
The paper finds a specific rule for this: The more "sticky" the pairing is (the stronger the superfluid "glue" holding the neutrons together, called the pairing gap), the more this geometric mixing happens.

• Surprise: In the past, scientists thought the superfluid density wouldn't change much if you changed the "glue."
• New Finding: In the neutron star crust, the superfluid density grows linearly with the strength of the glue. The stronger the pairing, the more the neutrons can use this "geometric mixing" to escape the drag of the lattice.

Why Did We Miss This Before?

The paper points out a mathematical "blind spot" in previous studies.

• The Mistake: Previous calculations treated the neutrons like simple, independent particles. They corrected the "position" of the particles but forgot to correct the "quantum state" (the way the particles are mixed together).
• The Fix: The author shows that to see the geometric contribution, you must account for how the quantum states themselves change when the system is moving. If you ignore this "mixing" of states, you only see the "conventional" drag. If you include it, you see the extra "geometric" flow.

Why Does This Matter?

1. Solving the "Glitch" Mystery: Pulsars (spinning neutron stars) sometimes suddenly speed up, like a figure skater pulling in their arms. This is called a glitch. Scientists have struggled to explain how the crust has enough "free" fluid to cause this. This paper suggests that because of the geometric contribution, there is much more free fluid in the crust than we thought. This might mean the crust alone is enough to cause glitches, without needing to involve the star's core.
2. Connecting the Universe to the Lab: The math used here is similar to what is used in "flat-band" superconductors on Earth (like in some exotic lab materials). This paper suggests that if we find a material on Earth where the superfluid density grows linearly with the pairing strength, we could use it as a miniature model to study neutron stars right here in the lab.

The Takeaway

The inner crust of a neutron star is a complex, multi-layered quantum system. The neutrons aren't just sliding on a floor; they are performing a complex dance where they mix between different energy levels. This "geometric" mixing allows them to flow much more freely than previously believed, potentially solving long-standing mysteries about how these stars spin and glitch.

In short: The neutrons aren't just stuck to the chairs; they are using a secret "quantum elevator" to bypass the traffic, and the stronger their bond, the faster they can take the elevator.

Technical Summary

Here is a detailed technical summary of the paper "Emergence of the geometric contribution to the superfluid density in the inner crust of neutron stars" by Giorgio Almirante.

1. Problem Statement

The inner crust of neutron stars consists of a periodic lattice of nuclear clusters (protons and neutrons) immersed in a gas of unbound neutrons and relativistic electrons. A critical parameter for understanding the hydrodynamical and thermodynamical properties of these stars (including pulsar glitches) is the superfluid density ($\rho_S$) of the neutron component.

Previous theoretical models, particularly those based on linear response theory within the BCS approximation, suggested that the superfluid fraction in the inner crust was weakly dependent on the pairing gap ($\Delta$) and that neutrons were strongly "entrained" (dragged) by the nuclear lattice. However, recent self-consistent calculations revealed a much higher superfluid fraction than previously expected.

The core problem addressed in this work is the neglect of the "geometric contribution" to the superfluid density in earlier derivations. While this contribution has been studied in condensed matter systems (specifically flat-band superconductors), its emergence and behavior in the inner crust of neutron stars—which features a highly non-trivial band structure with many energy bands crossing the Fermi energy—remained unclear. The author seeks to clarify how this geometric term arises, its dependence on the pairing gap, and why previous approximations failed to capture it.

2. Methodology

The author employs linear response theory within a mean-field framework to analyze a 3D periodic system (the neutron star crust) subjected to a stationary flow (Galilean boost).

• Formalism:

  • The system is described using the Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) formalisms.
  • The unperturbed system is defined by the HF eigenvalue problem $h|\alpha_k\rangle = \epsilon_{\alpha_k}|\alpha_k\rangle$, where $\alpha$ represents band indices and $k$ is the Bloch momentum.
  • A perturbation $V' = -\mathbf{v} \cdot \mathbf{p}$ is introduced to simulate a stationary flow.
  • The superfluid density is extracted by comparing the induced average momentum density $\langle \mathbf{j} \rangle$ with the two-fluid model expression: $\langle \mathbf{j} \rangle = (\langle \rho \rangle - \rho_S)m\mathbf{v}$.

• Perturbation Theory Approaches:
  The paper compares three distinct levels of approximation to compute the correction to the density matrix ($\rho'$):

  1. Full HFB Linear Response: Perturbation theory is applied to the full Bogoliubov quasi-particle states, accounting for corrections to both the single-particle states and the Bogoliubov coefficients ($u, v$).
  2. HF Linear Response: Perturbation theory applied to normal HF states (no pairing).
  3. Hybrid Approximation: Perturbation theory applied to HFB equations, but only correcting the single-particle states while keeping the Bogoliubov coefficients unperturbed (effectively neglecting the mixing of bands induced by the gap).

• Band Structure Analysis:
  The study focuses on the regime where many bands cut the Fermi energy (conduction bands), a characteristic of the neutron star inner crust, contrasting it with condensed matter systems where only a few bands near the Fermi surface are relevant.

3. Key Contributions and Results

A. Derivation of the Geometric Contribution

The author derives the total superfluid density $\rho_S$ as the sum of two terms:

1. Conventional Contribution ($\rho_{S}^{conv}$): A sum over single bands, dependent on the derivative of the energy dispersion.
2. Geometric Contribution ($\rho_{S}^{geom}$): A term arising from the mixing of different bands ($\alpha \neq \beta$), involving off-diagonal momentum matrix elements and a geometric tensor $g_{\alpha\beta}^{ij}$ related to the manifold of single-particle states.

$$\rho_{S}^{geom} \propto \sum_{\alpha \neq \beta} \frac{\Delta^2}{E_\alpha E_\beta (E_\alpha + E_\beta)} p_{\alpha\beta}^i p_{\beta\alpha}^j$$

B. Dependence on the Pairing Gap ($\Delta$)

A crucial finding is the behavior of the geometric contribution for small pairing gaps in systems with many bands crossing the Fermi energy:

• Linear Dependence: Unlike the conventional contribution (which is often constant or weakly dependent), the geometric contribution exhibits a linear dependence on the magnitude of the pairing gap ($\rho_{S}^{geom} \propto \Delta$) when many bands are involved.
• Physical Mechanism: This linear dependence arises because the pairing gap mixes empty and filled bands. In the limit where single-particle energies $|\xi| \ll \Delta$ or $|\xi| \approx \Delta$, the contribution scales linearly with $\Delta$.
• Numerical Verification: Figure 2 in the paper confirms this linear relationship numerically for baryon densities typical of the inner crust ($\rho_b = 0.030$ and $0.055 \text{ fm}^{-3}$).

C. Identification of the Source of Previous Errors

The paper clarifies why earlier works (e.g., Carter et al., Chamel) missed the geometric contribution:

• In the Hybrid Approximation (Section IV.B), where corrections are applied only to the single-particle states but the Bogoliubov coefficients ($u, v$) are treated as unperturbed, the geometric term vanishes.
• The geometric contribution strictly requires the perturbative correction to the Bogoliubov coefficients, which accounts for the mixing of bands induced by the pairing interaction.
• Previous derivations effectively neglected the non-diagonal terms in the average momentum density that arise from this band mixing.

D. Interpretation of Conduction vs. Geometric Density

• In the normal phase (HF), the response to a flow splits into "conduction density" (particles free to move, bands crossing the Fermi level) and a rigid contribution (particles moving collectively).
• In the superfluid phase, even if the conduction density is negligible (e.g., in flat-band systems or if bands are fully filled), a non-zero superfluid density can exist solely due to the geometric contribution, provided pairing correlations are strong enough to mix the bands.

4. Significance and Implications

• Resolution of the Pulsar Glitch Problem: The inclusion of the geometric contribution leads to a significantly higher superfluid fraction in the inner crust. This supports the hypothesis that the superfluidity of the crust alone is sufficient to explain the angular momentum transfer observed in pulsar glitches (specifically the Vela pulsar), removing the need to invoke superfluidity in the core or modify existing glitch models.
• Astrophysical Modeling: A reliable estimate of the superfluid density is essential for modeling:
  • Low-lying phonon modes and thermal evolution of neutron stars.
  • Frequencies of star oscillation modes.
  • Hydrodynamical properties of the crust.
• Bridge to Condensed Matter: The work establishes a theoretical link between neutron star physics and condensed matter systems (multi-band superconductors and flat-band systems). It suggests that experimental platforms mimicking the inner crust's band structure could be used to study these geometric effects in the lab.
• Theoretical Clarification: The paper provides a rigorous derivation showing that the geometric contribution is not an artifact of the BCS approximation but a fundamental consequence of band mixing in periodic superfluids, which must be accounted for in any self-consistent calculation involving many bands.

Conclusion

Giorgio Almirante demonstrates that the geometric contribution to the superfluid density is a dominant factor in the inner crust of neutron stars due to the complex, multi-band nature of the system. By showing that this contribution scales linearly with the pairing gap and arises from the perturbative mixing of Bogoliubov states, the paper corrects previous underestimations of the superfluid fraction, offering a more robust explanation for pulsar glitches and a refined framework for neutron star hydrodynamics.