Physicists prove that neutron star vibrations follow harmonic oscillator rules, even in Einstein's gravity

What is inside a neutron star? We know the outer layers contain neutrons, protons, and electrons crushed to densities several times that of an atomic nucleus. But the core remains, as physicists put it, a giant question mark. Some theories predict a quark-gluon plasma. Others suggest quantum superfluid or superconducting phases. The only place in the universe where such matter exists at the relevant temperatures is inside neutron stars themselves -- and we cannot exactly crack one open to look.

Gravitational waves offer a potential window. When two neutron stars spiral toward each other in a binary system, each star's gravity tidally deforms its partner, exciting internal vibrations that leave imprints on the gravitational waves they emit. If physicists can decode those imprints, they might be able to determine what the stars are made of.

But decoding requires a solid theoretical foundation, and a key piece of that foundation has been missing. Until now.

The completeness problem

In Newtonian gravity, a tidally deformed body responds by vibrating in modes -- oscillation patterns that behave like damped springs. These modes form a complete set, meaning any possible tidal response can be expressed entirely in terms of these modes, with nothing left out. This completeness is what makes Newtonian tidal models reliable: you know you are not missing part of the physics.

Scientists have long hoped the same would hold in Einstein's general relativity. But inspiraling neutron stars are profoundly relativistic objects. They are extraordinarily dense, approach speeds of 40% the speed of light before merging, and strongly distort the spacetime around them. The Einstein field equations governing their behavior are far more complex than Newton's, and three specific problems had blocked any proof of mode completeness.

First, separating the tidal effect of one star from the other in a binary system is technically difficult. Second, a star's own gravity modifies the equations inside and outside its boundary -- a complication absent in Newtonian theory. Third, the system constantly loses energy to gravitational radiation, which means the modes should not, strictly speaking, be complete at all.

Matched asymptotic expansion

A team led by Nicolas Yunes at the University of Illinois Urbana-Champaign, with collaborators at UC Santa Barbara, Montana State University, and the Tata Institute of Fundamental Research, solved these problems by dividing the space around the neutron star into zones. Near the star, gravity is strong. Far away, gravity is weak. The researchers found approximate solutions in each zone separately, then stitched them together using a technique called matched asymptotic expansion.

The far-zone decomposition handled the radiation problem: by restricting their analysis to the near zone and treating radiation as a small correction, they eliminated it from the mode analysis. The near-zone analysis captured the tidal field. And by manipulating the Einstein-Euler equations (which describe how matter generates gravitational fields and evolves in spacetime), they showed that the interior tidal field drives oscillations that produce damped harmonic oscillator modes -- exactly as in Newtonian theory.

"We showed two major things," said lead author Abhishek Hegade, now a postdoctoral scholar at Princeton. "First, we were able to subtract off radiation, finding that a neutron star's modes do indeed form a complete set. Second, we found that if you consistently solve a certain set of equations using a tidal field that's sufficiently 'smooth,' it's a solution to the interior of the star, and you can do all the same things in general relativity as in Newtonian gravity."

The result was published as an Editors' Suggestion in Physical Review Letters on February 18, 2026.

What this makes possible

With the proof in hand, researchers can now build gravitational wave models that fully capture neutron star tidal responses without worrying about missing contributions. The mode frequencies and decay times carry information about the star's internal composition -- including, potentially, whether quark cores or phase transitions exist deep inside.

"If we can understand the mode frequencies of oscillation and their decay times, we might be able to determine the composition of neutron stars in a regime not accessible on Earth," Yunes said.

Waiting for better detectors

There is a practical barrier, though. Current gravitational wave detectors, including LIGO, do not have sufficient sensitivity at the high frequencies where most neutron star oscillation mode information resides. The signal-to-noise ratios from LIGO's most recent neutron star merger detection (in 2017) are too low to reveal the features this framework captures.

Next-generation detectors, expected to come online within the next few years, should improve both sensitivity and frequency range. Until then, the team plans to extend their framework to rotating neutron stars (most spin rapidly), nonlinear tidal forces, and non-gravitational fields like magnetic fields.

"We've figured out the hard part -- gravity," Hegade said. "Now it's just a matter of applying our models to more realistic configurations."