The concept of tidal Love numbers dates back more than a century, when the British mathematician A. E. H. Love introduced dimensionless parameters describing how the Earth yields to the Moon’s pull — how much its surface bulges vertically, laterally, and how its gravitational potential warps in response. T. Shida later added a third parameter for lateral deformation. For decades these numbers belonged to geophysics. That changed in 2008, when Éanna Flanagan and Tanja Hinderer showed that the same mathematics applies to neutron stars in full general relativity, and that the tidal deformation leaves a measurable imprint on the gravitational waves emitted by merging binaries, entering the waveform phase at the fifth post-Newtonian order. Suddenly, Love numbers became a tool for reading the interior structure of stars light-years away.

Black holes, by contrast, are eerily rigid. Subsequent calculations demonstrated that the tidal Love numbers of Schwarzschild black holes vanish exactly, and the result extends to slowly rotating black holes as well. Lacking internal structure or any elastic mechanism to respond, a black hole simply refuses to bulge. This vanishing is more than a curiosity: it cleanly separates black holes from every other kind of compact object, and it means any nonzero tidal signal in a gravitational waveform is evidence that the object involved is not a black hole. That makes exotic alternatives — boson stars, gravastars, axion stars, Proca stars — prime targets for tidal analysis, because their Love numbers differ in both magnitude and sign from those of neutron stars.

Boson stars themselves trace back to the late 1960s, when David Kaup showed that a complex scalar field coupled to Einstein gravity could form a self-gravitating, stable lump — a Klein-Gordon geon — that resists gravitational collapse. Ruffini and Bonazzola obtained equivalent solutions from a quantized real scalar field around the same time. In the decades since, the boson-star family has grown to include self-interacting, charged, rotating, and oscillating varieties. These objects are more than mathematical toys: they can mimic black holes in some observational channels, they are serious candidates for dark matter, and their tidal Love numbers are known to be smaller in magnitude than those of neutron stars, offering a potential observational discriminator.

What the Lanzhou group — Xin-Lei Zhao, Jun-Ru Chen, and Yong-Qiang Wang — did differently was to consider a boson star made of two complex scalar fields simultaneously: one in its ground state, with no radial nodes in its wavefunction, and one in the first excited state, carrying a single node. Such multi-state boson stars were first constructed by Bernal and collaborators in 2010. On their own, excited-state boson stars are generally believed to be unstable, decaying into the ground state or collapsing into black holes under perturbation. But mixing the two states can produce stable configurations — a genuinely multi-field compact object, in contrast to the single-field objects that dominate the literature.

The team solved the coupled Einstein-Klein-Gordon equations numerically under two scenarios. In the synchronized case, both fields share the same frequency while their masses differ; in the nonsynchronized case, the frequencies differ while the masses are equal. In both settings the solutions split into single-branch and double-branch families. In the single-branch synchronized case, each frequency corresponds to exactly one configuration, and as the frequency varies the star smoothly interpolates between a pure excited-state star at the minimum frequency and a pure ground-state star at the maximum. In a narrow window of the mass ratio — between roughly 0.7922 and 0.7976 — the structure bifurcates, and each frequency supports two distinct solutions, a phenomenon reminiscent of multi-state Dirac stars. The authors assessed stability using the binding energy, the difference between the star’s total ADM mass and the sum of the Noether charges weighted by the field masses. Negative binding energy signals stability; positive signals decay. The verdict: double-branch solutions are always unstable, while stable configurations survive only in the single-branch families, and only when the two field masses are sufficiently close or the ground-state frequency is sufficiently high.

The centerpiece of the study is the calculation of the quadrupolar, or l = 2, tidal Love numbers, computed separately for the electric-type response — induced by even-parity perturbations of the metric — and the magnetic-type response, induced by odd-parity perturbations. Following Thorne’s framework for multipole moments of an arbitrary spacetime, the team extracted the induced mass and current multipoles from the asymptotic expansion of the perturbed metric, solving the coupled perturbation equations for the metric functions and both scalar-field perturbations on a 10,000-point finite-element grid with relative errors below one part in a hundred thousand. Only branches containing stable solutions were analyzed; the unstable families were set aside.

The results are striking. For stable single-branch stars, the electric tidal Love numbers start positive, grow with the star’s mass, and then — at a critical mass — abruptly jump to negative values, with their magnitudes then decreasing as the mass continues to rise. The authors call this a peak: near the transition, a tiny change in mass produces a dramatic swing in deformability, complete with a sign reversal from positive to negative feedback. The transition appears whenever the excited-state field mass exceeds roughly 0.891 times the ground-state mass, or the ground-state frequency exceeds roughly 0.777; below those thresholds the electric Love numbers stay positive throughout. The magnetic Love numbers tell a calmer story: they are always negative, and their absolute values remain smaller than those of the electric numbers. In the double-branch nonsynchronized case, the electric Love numbers remain positive but develop spiral-like tails near the branch-turning point, while the magnetic numbers climb monotonically.

Two broader conclusions stand out. First, the Love numbers of multi-state boson stars are larger in absolute value than those of ordinary ground-state boson stars — the presence of the excited-state component makes the star measurably softer, more easily deformed by an external tidal field. Second, the sign flip itself is a signature tied directly to the excited state and to how close its mass sits to the ground state’s. Because black holes have identically vanishing Love numbers and neutron stars show positive, smoothly varying responses, a gravitational-wave signal carrying this kind of abrupt sign transition would point squarely at exotic, multi-field matter. The authors suggest natural extensions: adding rotation and self-interactions, applying the same machinery to multi-state Dirac stars and other multi-field models. As gravitational-wave detectors grow more sensitive, such exotic tidal fingerprints may move from theoretical prediction to observable test — turning the century-old mathematics of Earth’s ocean tides into a probe of what quantum fields, if any, hide inside the universe’s darkest objects.

Subject of Research: Tidal Love numbers of multi-state boson stars composed of ground-state and first-excited-state complex scalar fields

Article Title: Tidal Love numbers of multi-state Boson stars

Article References: Zhao, X.-L., Chen, J.-R., & Wang, Y.-Q. (2026). Tidal Love numbers of multi-state Boson stars. The European Physical Journal C, 86(9), Article 1121. https://doi.org/10.1140/epjc/s10052-026-16369-6

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16369-6

Keywords: boson stars, tidal Love numbers, gravitational waves, compact objects, scalar fields, general relativity, black hole mimickers, dark matter, excited states, binding energy, tidal deformability, theoretical physics