Tidal deformations of general-relativistic multifluid compact stars

This paper presents a fully general-relativistic multifluid framework for modeling adiabatic tidal deformations in compact stars, demonstrating that nondissipative mutual entrainment between fluid species leaves tidal deformabilities unchanged and thus produces no measurable effect on gravitational-wave signals during the inspiral phase.

The Gist

Imagine two cosmic dancers, a pair of neutron stars, spiraling toward each other in a slow, graceful waltz. As they get closer, they don't just pull on each other with gravity; they actually stretch and squish one another, like two blobs of dough being pulled apart by invisible hands. This stretching is called a tidal deformation.

For the last decade, scientists have been listening to the "music" of these dances using gravitational wave detectors. By measuring exactly how much the stars squish, we can figure out what they are made of inside. It's like trying to guess the ingredients of a cake just by listening to how it bounces when you poke it.

However, there's a problem. The standard recipe we use to describe these stars assumes they are made of a single, perfect, smooth fluid (like water). But in reality, the inside of a neutron star is more like a complex cocktail or a layered smoothie. It might have:

• Superfluid neutrons (neutrons that flow without friction).
• Superconducting protons.
• Maybe even dark matter mixed in.
• Different layers of "crust" and "core."

When you have multiple fluids interacting, they don't just flow past each other; they get "entangled." In physics, this is called entrainment. Imagine two people running on a treadmill. If they are holding hands (entrainment), when one speeds up, the other is dragged along, even if they aren't trying to. This interaction changes how the fluids move.

The Big Question

Scientists wondered: Does this "holding hands" (entrainment) change how the star squishes?
If it does, then our current models for reading gravitational waves are missing a huge piece of the puzzle. Some previous studies suggested that superfluidity (the frictionless flow) could change the squishiness by up to 20%. That would be a massive difference!

The Discovery: The "Ghost" Effect

In this paper, the authors (Ethan Carlier and Nicolas Chamel) built a super-advanced mathematical model to simulate these multi-fluid stars. They used a sophisticated framework developed by physicist Brandon Carter to handle all the complex interactions between different types of fluids.

They asked: "If we include all these messy, entangled interactions, does the star's shape change differently when it gets stretched?"

The answer is a surprising "No."

Here is the analogy:
Imagine you are stretching a piece of Jell-O that has different flavors mixed inside (strawberry, lime, grape).

• The Old View: If the flavors are "entangled" (stuck together), stretching the Jell-O should feel different than stretching plain Jell-O.
• The New Discovery: The authors proved that for the specific type of stretching that happens slowly during the early part of the dance (the "adiabatic" phase), the entangled flavors don't matter at all.

The "squishiness" of the star depends only on the total energy and pressure of the mixture, not on how the different fluids are holding hands with each other. The entrainment effect is like a ghost: it exists and is real, but it leaves no fingerprint on the gravitational wave signal we detect during the early inspiral.

Why This Matters

1. Simplifying the Search: Scientists can stop worrying about calculating the complex "entanglement" between superfluids or dark matter when trying to interpret gravitational waves from the early stages of a merger. They can use simpler models and still get the right answer.
2. Resolving Conflicts: Some previous studies said superfluidity changes the result; others said it doesn't. This paper settles the debate: if the star is in a stable, slow-moving state, superfluidity does not change the tidal deformability. The differences seen in earlier studies were likely due to calculation errors or different assumptions.
3. Dark Matter: Even if neutron stars are hiding dark matter in their cores, and even if that dark matter is "entangled" with normal matter, it won't change the gravitational wave signal in a way we can detect right now.

The Bottom Line

The universe is full of complex, multi-layered fluids. But when it comes to the slow, gentle stretching of two stars before they crash, nature simplifies things. The "messy" interactions between the different fluids cancel out, leaving us with a signal that depends only on the star's overall weight and pressure.

So, while the inside of a neutron star is a chaotic, superfluid, multi-fluid party, the gravitational waves it sends out during the early dance tell us a much simpler story. We don't need to know who is holding hands to know how the star will squish.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the need for increasingly accurate theoretical models of compact stars (neutron stars, white dwarfs, and exotic objects) in the context of gravitational-wave (GW) astronomy.

• Context: Third-generation GW detectors (e.g., Einstein Telescope, Cosmic Explorer) will measure tidal deformabilities with unprecedented precision, potentially revealing subleading tidal contributions.
• The Gap: Most existing models assume a single, barotropic, perfect fluid. However, realistic compact stars often contain multiple interacting fluid components (e.g., superfluid neutrons, superconducting protons, dark matter admixtures, or hyperons).
• The Specific Issue: Previous studies on superfluid neutron stars yielded contradictory results regarding the impact of mutual entrainment (a non-dissipative coupling where the momentum of one fluid species depends on the velocity of another) on tidal deformabilities. Some studies suggested a ~20% change, while others argued for no effect.
• Goal: To develop a fully general-relativistic framework for adiabatic tidal deformations of compact stars composed of an arbitrary number ($N$) of interacting fluids using Carter's variational formalism, and to definitively determine how entrainment affects tidal Love numbers.

2. Methodology

The authors employ Carter's multifluid variational formalism, a powerful and general approach to relativistic hydrodynamics.

• Formalism:

  • The system is described by a master function $\Lambda(n^\mu_x, g_{\mu\nu})$, which acts as the Lagrangian density. It depends on the metric $g_{\mu\nu}$ and the density currents $n^\mu_x$ for $N$ fluid species.
  • The master function encapsulates the microphysics, including entrainment effects via off-diagonal terms in the current couplings.
  • The authors use a constrained action principle with Lagrangian displacements to derive the equations of motion, ensuring conservation of particle fluxes ($\nabla_\mu n^\mu_x = 0$).

• Equilibrium Configuration:

  • They derive the hydrostatic equilibrium equations for a static, spherically symmetric star.
  • These equations generalize the Tolman-Oppenheimer-Volkoff (TOV) equations to the multifluid case, resulting in a system of $N+2$ coupled differential equations for the metric potentials and density profiles.

• Perturbation Theory:

  • The authors introduce stationary (zero-frequency) tidal perturbations to the equilibrium configuration.
  • They decompose perturbations into Polar (Gravitoelectric) and Axial (Gravitomagnetic) sectors using spherical harmonics.
  • They derive the linearized Einstein equations and the perturbed Euler equations for an arbitrary number of fluids ($N$).
  • Crucially, they derive a master differential equation for the metric perturbation function $H_\ell$ (polar) and $h_\ell$ (axial) that includes a specific term $J_N$, which encapsulates the microphysical response of the fluid mixture.

• Analytical Equation of State (EOS):

  • To analyze the entrainment effect, they introduce an analytical representation of the master function $\Lambda$ as a power series expansion around the static limit ($n^2_{xy} \approx n_x n_y$).
  • This allows them to separate static thermodynamic terms from dynamic entrainment terms (coefficients $X^k_{xy}$ for $k>0$).

3. Key Contributions

1. Generalized Formalism: The paper provides the first fully general-relativistic description of adiabatic tidal deformations for compact stars with an arbitrary number ($N$) of interacting fluids. Explicit expressions for the perturbation equations are provided for $N=1, 2,$ and $3$.
2. Derivation of $J_N$: The authors derive a general expression for the quantity $J_N$, which appears in the polar perturbation equation. This term determines the effective stiffness of the star against tidal deformation.
3. Resolution of Entrainment Controversy: By evaluating the perturbation equations in the static background limit, the authors prove that mutual entrainment has no effect on adiabatic tidal deformabilities.
   • They demonstrate that the term $J_N$ reduces to the standard barotropic form $J_N = (E+P)/(dP/dE)$, where $E$ is the static energy density and $P$ is the pressure.
   • The entrainment coefficients (which depend on relative velocities) vanish in the static limit relevant to adiabatic tides.
4. Unified Treatment: The framework unifies the treatment of superfluid neutron stars, dark matter admixed stars, and hot compact stars under a single mathematical structure.

4. Key Results

• Independence from Entrainment: The adiabatic tidal Love numbers ($k_\ell$) and deformabilities ($\Lambda_\ell$) are uniquely determined by the static equation of state $E(n_x)$. They are completely independent of the non-dissipative mutual entrainment between fluid species.
• Implication for Superfluidity: For cold superfluid neutron stars in beta equilibrium, the presence of superfluidity (and the associated entrainment) leaves no imprint on the gravitational wave signal during the early inspiral. The tidal response is identical to that of a single-fluid barotropic model with the same static EOS. This resolves the discrepancy with previous numerical studies (e.g., Refs. [51, 52]), suggesting those results were likely due to numerical artifacts or misinterpretations of the beta-equilibrium condition.
• Implication for Dark Matter: For compact stars admixed with dark matter, any entrainment between baryonic and dark matter components does not affect the tidal Love numbers. Current models that ignore this coupling are therefore not introducing an approximation error regarding adiabatic tides.
• Gravitomagnetic Sector: The result holds for both polar (electric) and axial (magnetic) tidal responses. The perturbation equations for the gravitomagnetic sector also reduce to the single-fluid form when evaluated on the static background.

5. Significance and Conclusion

• Theoretical Clarity: The paper establishes a rigorous theoretical baseline: for adiabatic tides, the complex microphysics of entrainment is irrelevant. Only the static thermodynamic properties (the relationship between energy density and pressure) matter.
• Observational Impact: This simplifies the interpretation of future GW data. When third-generation detectors measure tidal deformabilities, they will constrain the static Equation of State of dense matter directly, without needing to disentangle complex entrainment effects.
• Limitations and Future Work: The authors note that this result applies strictly to static backgrounds and adiabatic (slowly varying) tidal fields. If the tidal field becomes dynamical (exciting internal oscillation modes) or if the star is rotating, entrainment effects may become non-negligible.
• Conclusion: The development of more complex multifluid models is essential for understanding the internal structure of compact stars, but for the specific purpose of modeling adiabatic tidal deformations in binary inspirals, the standard single-fluid perfect-fluid approximation (with the correct static EOS) is sufficient. The "richness" of the multifluid description does not alter the tidal Love numbers.