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

The new work, published in The European Physical Journal C by Kazem Bitaghsir Fadafan of Shahrood University of Technology, Jesús Cruz Rojas of the National Autonomous University of Mexico, and Jonas Mager of the Vienna University of Technology, applies holography, the correspondence that links strongly coupled quantum field theories to weakly coupled gravitational theories in higher dimensions. The central obstacle in describing quark matter inside neutron stars is that the relevant regime of quantum chromodynamics, or QCD, is strongly coupled. First-principles lattice simulations fail there because of the notorious fermion sign problem at large densities, while traditional phenomenological approaches such as nucleon effective field theories struggle under the extreme conditions found in stellar cores. Holographic models offer a way around this impasse: instead of computing the equation of state of dense quark matter directly, one solves a comparatively tractable problem in a higher-dimensional gravitational dual and reads off the properties of the quantum field theory from the geometry.

The researchers built their model on the well-studied D3/D7-brane configuration of type IIB string theory, a bottom-up construction in which D7-branes embedded in the curved background of D3-branes describe quarks moving in a gauge theory. The key innovation lies in the dilaton, a scalar field that controls the running of the gauge coupling in the dual theory. Rather than keeping the dilaton constant, as in the simplest version of the model, the team adopted a phenomenological profile that interpolates smoothly and monotonically from its ultraviolet value to a finite, regular value in the infrared. Three parameters, denoted A, lambda and kappa, govern the precise shape of this profile, controlling the infrared value of the dilaton, the scale at which conformal symmetry is broken, and the steepness of the transition. This tailored profile captures essential features of QCD such as chiral symmetry breaking and produces a deconfined yet massive quark phase at finite density.

In the holographic setup, the embedding function of the D7-brane encodes the quark mass and condensate, while a gauge field living on the brane describes the quark chemical potential and density. Two distinct phases emerge: a Minkowski, or vacuum, phase with broken chiral symmetry and zero baryon density, and a quark phase in which chiral symmetry is restored and the density is nonzero. By computing the on-shell action of the brane, which the holographic dictionary identifies with the grand canonical potential, the researchers obtained the pressure as a function of chemical potential and hence the equation of state of the quark phase. At asymptotically large densities the model correctly reproduces the perturbative QCD behavior, with the pressure scaling as the fourth power of the chemical potential.

The team also attempted to describe the baryonic phase within the same holographic framework, modeling nucleons as D5-branes wrapped on a five-sphere. The nontrivial dilaton profile stabilizes these wrapped branes against collapse, which is itself a notable achievement. However, in the homogeneous, smeared approximation used here, the resulting nuclear equation of state turned out to be unrealistic, predicting pressures in sharp conflict with phenomenology just above the onset of the baryonic phase. The authors showed analytically that near the transition the chemical potential scales as the one-third power of the density, forcing a fourth-order transition that cannot reproduce the first-order behavior expected of isospin-symmetric QCD. They therefore discarded the holographic baryon phase and instead adopted the phenomenological equations of state of Hebeler and collaborators, derived from chiral effective field theory and constrained by nuclear physics and observation.

With the baryonic phase fixed phenomenologically and the quark phase supplied by holography, the researchers scanned the parameter space of the model, varying the AdS radius and the t’Hooft coupling. They found a region of parameters, with the AdS radius between roughly 0.015 and 0.02 inverse megaelectronvolts and the t’Hooft coupling between 1.9 and 3, that produces rather stiff equations of state with only weakly first-order transitions between baryonic and quark matter. This softness of the transition is crucial: it allows a smooth conversion from nuclear matter to quark matter inside a star rather than a violent discontinuity that would destabilize the configuration. When combined with the stiff Hebeler equation of state, three parameter choices produced equations of state supporting stable quark cores, with phase transitions occurring at chemical potentials around 400 megaelectronvolts.

To connect the microscopic equation of state to observable stellar properties, the team solved the Tolman-Oppenheimer-Volkoff equations of general relativity, which determine the structure of non-rotating compact stars. The resulting mass-radius curves show stars that begin their lives as ordinary nucleonic objects and develop quark cores as the central density increases. Stability requires that the mass increase with central energy density, and the stars with quark cores satisfy this criterion all the way up to the maximum of the mass-radius curve, beyond which a radial mode becomes unstable and the star collapses into a black hole. The maximum masses obtained for the three quark-star-supporting parameter choices are 2.17, 2.13 and 1.90 solar masses, the largest of which is broadly consistent with the heaviest precisely measured neutron star mass of about 2.35 solar masses, albeit somewhat favoring lower values.

The authors also computed the tidal deformability, the quantity that measures how easily a star is deformed by the gravitational field of a companion and that leaves a characteristic imprint on gravitational wave signals from binary neutron star inspirals. Once a star becomes heavy enough to develop a quark core, the tidal deformability drops rapidly, falling by roughly 80 units between the onset of the quark phase and the maximum allowed mass. This rapid decrease is a potentially observable signature of quark matter formation. However, the very stiffness of the baryonic phase required to obtain stable quark cores causes the model to overshoot the tidal deformability constraint from the GW170817 event at 1.4 solar masses, predicting a value near 950 compared with the observational bound of about 190 with large uncertainties. The authors note that a baryonic equation of state intermediate between the stiff and medium phenomenological cases would likely ease this tension while still permitting stable quark cores.

Several other results add nuance to the picture. The speed of sound in the quark phase remains below the speed of light, although the stiff baryonic phase violates a recently derived transport bound. More strikingly, the polytropic index of the holographic quark phase exceeds the value of 1.75 that has been proposed in the literature as a criterion for the onset of quark matter, reaching values as high as 2.5 just after the transition. This suggests that the criterion may be too restrictive, and it contrasts with predictions from the V-QCD holographic model, which places transitions at much larger energy densities with smaller polytropic indices. The authors emphasize that their model cannot definitively establish whether stable quark stars exist in nature, since the answer depends on parameter choices, but it demonstrates that such objects can in principle emerge from holographic models, opening a computational window onto their phenomenology. Future work will explore localized D5-brane configurations, which may improve the low-density baryonic description, along with neutrino transport and rotating stars with quark cores, as multi-messenger observations continue to tighten the constraints on the densest matter in the universe.

Subject of Research: Holographic modeling of stable massive quark stars and the equation of state of dense quark matter in compact star cores

Article Title: Properties of stable massive quark stars in holography

Article References: Bitaghsir Fadafan, K., Cruz Rojas, J., & Mager, J. (2026). Properties of stable massive quark stars in holography. The European Physical Journal C, 86(9), Article 1094. https://doi.org/10.1140/epjc/s10052-026-16338-z

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16338-z

Keywords: quark stars, holographic QCD, neutron stars, D3/D7-branes, equation of state, chiral symmetry breaking, tidal deformability, compact stars, gauge-gravity duality, quark matter, TOV equations, dense nuclear matter