Interaction between nuclear clusters and superfluid phonons in the neutron-star inner crust

This paper employs nuclear density functional theory to derive a microscopic description of the interaction between nuclear clusters and superfluid phonons in the neutron-star inner crust, revealing a significantly smaller coupling constant than previous hydrodynamical estimates due to the suppression of phonon amplitude within the clusters.

The Gist

The Setting: A Cosmic "Cheese"

Imagine the inside of a neutron star (the incredibly dense, dead core of a collapsed star). The "inner crust" of this star isn't just a solid rock; it's a bizarre, quantum kitchen.

Think of it like a giant block of Swiss cheese:

• The Holes: These are actually nuclear clusters (tiny, heavy balls of protons and neutrons) arranged in a perfect crystal grid.
• The Cheese: The space between the holes is filled with a neutron superfluid. This isn't normal water or milk; it's a quantum liquid that flows with zero friction, like a super-conductor but for matter.
• The Electrons: Floating around everywhere are electrons, acting like the air in the room.

The Problem: Two Types of "Waves"

In this cosmic cheese, there are two different ways things can wiggle or vibrate:

1. Lattice Phonons (The Cheese Wiggle): If you poke the grid of nuclear clusters, the whole grid shakes. It's like shaking a crystal lattice.
2. Superfluid Phonons (The Liquid Wiggle): The neutron superfluid can also ripple and wave, like sound moving through water.

The Big Question: How do these two wiggles talk to each other? If the grid shakes, does it make the liquid ripple? If the liquid ripples, does it shake the grid?

Scientists have known for a long time that they do interact, and this interaction is crucial for understanding "glitches" (sudden speed-ups) in spinning pulsars. However, previous estimates of how strong this connection is were based on big, rough approximations (like looking at the cheese from a satellite). They assumed the liquid and the grid were perfectly mixed and that the waves moved through everything equally.

The New Discovery: The "Force Field" Effect

The authors of this paper decided to look much closer, using a microscopic microscope (a super-computer simulation based on quantum physics). They asked: What actually happens to a wave when it tries to pass right next to one of those heavy nuclear clusters?

The Analogy: The Crowd and the VIP
Imagine a wave of people (the superfluid phonon) trying to walk through a hallway.

• Old View (Macroscopic): We assumed the people walk smoothly past a VIP (the nuclear cluster) sitting in a chair. We thought the VIP just sat there, and the people flowed around them easily.
• New View (Microscopic): The authors found that the VIP is actually surrounded by a force field. When the wave of people gets close to the VIP, they get scared or repelled. They don't just flow around; they flatten out and avoid the VIP entirely. The "amplitude" (the height of the wave) drops to almost zero right near the VIP.

The Result: The Connection is Much Weaker

Because the wave avoids the VIP, the "handshake" between the grid and the liquid is much weaker than we thought.

• The Math: The paper calculates a "coupling constant" (a number that tells us how strong the link is).
• The Surprise: Their number is significantly smaller (about 5 to 7 times smaller) than previous estimates.
• Why? Because the superfluid wave is "suppressed" (squashed down) right where it needs to touch the nuclear cluster to transfer energy. It's like trying to shake hands with someone, but they are wearing a giant, invisible bubble suit that pushes your hand away before you can touch them.

Why Does This Matter?

This isn't just about math; it changes how we understand neutron stars.

1. Pulsar Glitches: When a pulsar suddenly speeds up, it's often blamed on the "unpinning" of these quantum waves from the grid. If the connection is weaker than we thought, the physics of how these glitches happen needs to be rewritten.
2. Heat and Sound: The way heat and sound travel through the star's crust depends on this mixing. If the mixing is weak, the star might cool down or conduct heat differently than models predicted.

The Bottom Line

The authors built a detailed, microscopic model of a single nuclear cluster sitting in a sea of superfluid neutrons. They found that the superfluid waves are afraid of the clusters and stay away from them. This "fear" (or physical suppression) means the two parts of the star's crust are not as tightly coupled as we previously believed.

In short: We thought the grid and the liquid were holding hands tightly. It turns out they are actually holding hands through a thick, bouncy mattress, making the connection much weaker than expected.

Technical Summary

1. Problem Statement

The inner crust of neutron stars consists of a lattice of nuclear clusters immersed in a superfluid neutron liquid and a background of relativistic degenerate electrons. A critical aspect of the dynamics in this region is the interaction between lattice phonons (vibrations of the nuclear clusters) and superfluid phonons (collective excitations of the neutron superfluid).

• Context: This coupling influences the heat capacity, thermal conductivity, elastic properties, and potentially the "glitch" phenomenon (sudden changes in pulsar rotation) via vortex pinning.
• The Gap: Previous studies have relied on macroscopic approaches (hydrodynamics and effective field theory) to estimate the coupling strength. These models treat parameters phenomenologically based on bulk thermodynamic properties.
• The Issue: The microscopic origin of this coupling and the precise value of the effective coupling constant have remained unclear. Previous microscopic studies suggested the coupling might be weak due to the non-penetration of phonons into clusters, but a rigorous quantitative derivation linking microscopic structure to the macroscopic coupling constant was missing.

2. Methodology

The authors developed a framework to derive the coupling constant from first principles using Nuclear Density Functional Theory (DFT) and the Quasiparticle Random-Phase Approximation (QRPA).

A. Theoretical Framework

1. Microscopic Derivation: They started with a microscopic description where nuclear clusters are treated as a lattice. The interaction between lattice vibrations and the neutron superfluid is mediated by the self-consistent mean field potential.
2. Effective Hamiltonian: They formulated an effective Hamiltonian in the long-wavelength limit ($k \to 0$) describing the mixing of lattice and superfluid phonons. This Hamiltonian includes a coupling term characterized by a constant $g_0$.
3. Matching Condition: To determine $g_0$, they established a matching condition between the microscopic response (involving the actual cluster potential) and the effective theory response (involving a simplified effective potential).

B. Computational Strategy

Since a full linear-response calculation for a periodic lattice is computationally prohibitive and not yet available, the authors employed a single-cluster approximation:

• Setup: A single nuclear cluster (specifically $Z=28$) is embedded in a large spherical box ($R=50$ fm) containing the neutron superfluid.
• Microscopic Calculation: They used the Time-Dependent Kohn-Sham-Bogoliubov-de Gennes (TDKS-BdG) formalism with the Skyrme energy density functional (SLy4) and a density-dependent delta pairing interaction.
• Linear Response: They calculated the transition density ($\delta\rho$) and the response to the dipole operator ($rY_{1M}$) and the potential gradient ($\partial v_0 / \partial r$) using QRPA.
• Effective Theory Matching: They compared the microscopic matrix elements with those derived from an effective potential $v_{eff}$ (parameterized as a Gaussian) acting on a uniform superfluid. The strength of $v_{eff}$ was tuned to reproduce the microscopic matrix elements.

3. Key Contributions

1. Microscopic Determination of Coupling: The paper provides the first rigorous microscopic derivation of the coupling constant between lattice and superfluid phonons in the neutron star crust, moving beyond phenomenological estimates.
2. Quantification of Suppression: The study quantifies the strong suppression of the superfluid phonon amplitude inside and immediately around the nuclear cluster. This suppression is the primary physical mechanism reducing the coupling strength.
3. Resolution of Discrepancies: The work explains why previous macroscopic estimates (hydrodynamics) overestimated the coupling: they failed to account for the distortion of the phonon wavefunction near the cluster surface.
4. Effective Potential Construction: The authors successfully constructed an effective interaction potential ($\bar{v}_{eff}(0)$) that bridges the gap between the microscopic DFT description and the long-wavelength effective field theory.

4. Key Results

• Coupling Strength: The derived coupling constant $g_0$ (or equivalently the effective interaction volume integral $\bar{v}_{eff}(0)$) is significantly smaller than previous estimates.
  • Microscopic Result: $\bar{v}_{eff}(0) \approx 7.1 \times 10^3 \text{ MeV fm}^3$.
  • Comparison with Naive Estimate: This is roughly 7 times smaller than a naive estimate using the bare Hartree-Fock potential ($\approx 4.7 \times 10^4 \text{ MeV fm}^3$).
  • Comparison with Hydrodynamics: It is roughly 5 times smaller than estimates based on hydrodynamic models (e.g., Kobyakov & Pethick) when entrainment effects are neglected.
• Physical Mechanism: The reduction is attributed to the distortion of the phonon wavefunction.
  • The transition density (phonon amplitude) is almost zero inside the nuclear cluster ($r \lesssim 6$ fm).
  • The amplitude is also significantly suppressed near the cluster surface.
  • Since the coupling matrix element depends on the spatial overlap between the phonon wavefunction and the cluster potential, this suppression drastically reduces the effective interaction.
• Weak Mixing: The resulting coupling strength ($g_0/c \approx 4.3 \times 10^{-3}$) is much smaller than the superfluid phonon velocity ($v_s/c \approx 4.9 \times 10^{-2}$), indicating that the mixing between lattice and superfluid phonons is weak.

5. Significance

• Astrophysical Implications: The reduced coupling strength suggests that the interaction between the lattice and the superfluid is weaker than previously thought in macroscopic models. This has direct consequences for:
  • Glitch Mechanics: The efficiency of angular momentum transfer between the superfluid and the crust.
  • Thermal Evolution: The thermal conductivity and heat capacity of the inner crust, which rely on phonon scattering.
  • Elastic Properties: The rigidity and vibrational modes of the crust.
• Theoretical Advancement: The study validates the use of microscopic DFT/QRPA to determine parameters for effective field theories. It demonstrates that "naive" macroscopic approximations can fail significantly when the microscopic structure (specifically the exclusion of phonons from the cluster interior) is ignored.
• Future Directions: The authors note that the entrainment effect (drag between normal and superfluid components) was not explicitly included in this specific calculation but is known to enhance coupling in hydrodynamic models. Future work will aim to incorporate entrainment within this microscopic framework to provide a complete picture of crust dynamics.

In summary, this paper establishes that the interaction between nuclear clusters and superfluid phonons is governed by a microscopic suppression mechanism, leading to a coupling constant significantly smaller than previously assumed in macroscopic models. This finding necessitates a re-evaluation of neutron star crust dynamics based on more accurate microscopic inputs.