Neutrino oscillations are one of the most remarkable phenomena in modern particle physics. Neutrinos come in three flavors, electron, muon and tau, but they do not travel as pure flavors. Instead, each flavor is a quantum superposition of three mass states, mixed together by the Pontecorvo-Maki-Nakagawa-Sakata matrix. As a neutrino races through space, these mass states accumulate slightly different quantum phases, and the relative phase difference determines the probability that the particle will be detected as a different flavor than the one it started as. The oscillation phase depends on the neutrino’s energy, the squared mass differences between the eigenstates, and the distance traveled. In flat spacetime the formula is elegant and well tested by experiments such as Super-Kamiokande, IceCube and Borexino. Curved spacetime, however, changes the game, because the phase must now be computed along the curved geodesic that the particle actually follows.

The team modeled a static, spherically symmetric Schwarzschild black hole immersed in what they call DMPF, dark matter within a phantom field background. This geometry, built from an Einstein-Hilbert action augmented with a scalar field, adds a logarithmic correction to the standard black hole metric, governed by a parameter a that encodes the density of the entire dark sector. The energy density of this dark component falls off as one over the cube of the distance from the hole, a profile consistent with the nearly flat rotation curves observed in spiral galaxies. By varying the dimensionless ratio a/M, where M is the black hole mass, the researchers could dial the strength of the dark matter influence up and down and watch how the neutrino oscillation patterns responded.

The first major result is a striking cancellation. For neutrinos traveling purely radially, straight outward from the black hole, the gravitational contributions to the oscillation phase exactly cancel. The phase difference between mass eigenstates ends up depending only on the radial distance, the mass-squared splitting and the neutrino energy, precisely as it would in flat spacetime. In other words, a neutrino beaming directly away from a black hole carries no memory of the gravitational maelstrom it just escaped, at least not in its oscillation phase. This result follows from expanding the covariant phase integral, built from the canonical momentum along the geodesic, in the ultra-relativistic limit where the neutrino mass is tiny compared with its energy, and it holds even in the presence of the dark matter logarithmic term.

The story changes dramatically for neutrinos whose paths are bent. When a trajectory is deflected by the gravitational field, characterized by a nonzero angular momentum and an impact parameter b, extra phase terms survive. The researchers derived analytical expressions showing that the dark matter parameter enters both as a small correction to the effective mass-splitting coefficient and as an additive logarithmic contribution to the total accumulated phase. When gravitational lensing splits the journey into two segments, from source to the point of closest approach and from there to the detector, the phase acquires further corrections involving the lens mass and the logarithm of the source and detector distances. Gravitational lensing thus does more than bend the neutrino’s path; it imprints measurable modifications on the interference pattern that determines which flavor arrives at the detector.

To make these abstract calculations concrete, the team computed two-flavor transition probabilities, electron to muon, for realistic astrophysical configurations. They considered three lensing objects: the Sun, Sagittarius A*, the four-million-solar-mass black hole at the center of our galaxy, and M87*, the 6.5-billion-solar-mass giant famously imaged by the Event Horizon Telescope. For Sagittarius A*, Event Horizon Telescope shadow observations constrain the dark sector parameter to values up to roughly 0.15 at two sigma, while for M87* the bound reaches about 0.33. The corresponding neutrino energies span an enormous range: solar neutrinos carry around 10 megaelectronvolts, astrophysical neutrinos detected by IceCube from sources such as the blazar TXS 0506+056 and the Seyfert galaxy NGC 1068 reach 10^8 to 10^10 megaelectronvolts, and hypothetical cosmogenic neutrinos, produced when ultra-high-energy cosmic rays collide with cosmic microwave background photons, could climb to the exaelectronvolt scale near M87*.

The numerical results reveal a rich structure. Plotting the transition probability as a function of the azimuthal angle of the source, the curves shift noticeably as the dark sector parameter is increased from zero, the pure Schwarzschild case, to a/M values of 0.05 and 0.1. The probability is highly sensitive to the dark sector density parameter, and it also depends on whether the neutrino mass hierarchy is normal or inverted. Intriguingly, the team uncovered a degeneracy: different combinations of lens mass and dark matter parameter can produce the same oscillation probability. For a fixed probability, achieving the same flavor transition under an inverted hierarchy requires a slightly larger a/M than under a normal hierarchy, and this difference grows with mass. This degeneracy is both a challenge and an opportunity, since breaking it with independent mass measurements could turn neutrino detectors into dark matter probes.

The study also tackled a subtler question: what happens when neutrinos are treated realistically as Gaussian wave packets rather than infinite plane waves? Wave packets of different mass states gradually separate in space, causing the oscillation pattern to wash out, a process called decoherence. Because gravity dilates proper time, neutrinos in curved spacetime must travel farther before their wave packets lose overlap. The researchers derived the effective damping factor for lensed trajectories and found that the dark matter background exerts only a minimal influence on decoherence. Instead, the absolute neutrino mass is the dominant factor determining the coherence length and the duration over which oscillations persist. Once decoherence sets in, the dark sector parameter has a negligible effect on the transition probability envelope, meaning that any dark matter signature must be extracted from neutrinos that are still coherent when detected.

The implications are tantalizing. In the Solar System, the dark matter density bound at one astronomical unit is so stringent that the parameter a/M is capped at roughly 7.5 times 10 to the minus 12, making dark matter effects on solar neutrino oscillations utterly unobservable with current techniques. But around supermassive black holes like Sagittarius A* and M87*, the allowed dark sector parameters are many orders of magnitude larger, and the extreme gravitational lensing amplifies the phase corrections. Future neutrino astronomy, whether through next-generation IceCube extensions, kilometer-cubic detectors in the Mediterranean, or hypothetical cosmogenic neutrino observatories, could in principle compare observed flavor ratios with the predictions of standard general relativity and detect the fingerprint of the dark halo surrounding these monsters.

Of course, this remains a theoretical edifice, resting on a two-flavor approximation and idealized geometries, and no one has yet observed a neutrino oscillation modified by dark matter. But the work demonstrates that the intersection of quantum interference and strong-field gravity is not merely a curiosity; it is a potential measurement channel. Any departure from the predictions of pure general relativity in neutrino flavor ratios arriving from the vicinity of compact objects would constitute compelling evidence for modified gravity or for the dark environments in which black holes are embedded. As detectors grow more sensitive and theoretical models more refined, the humble neutrino, the particle that barely interacts with anything, may become one of humanity’s sharpest tools for mapping the invisible universe that surrounds every black hole.

Subject of Research: Effects of dark matter in a phantom field background on neutrino flavor oscillations near black holes

Article Title: Probing dark matter and phantom field effects on neutrino oscillations around black holes

Article References: Ergashov, I., Narzilloev, B., Hussain, I., & Ahmedov, B. (2026). Probing dark matter and phantom field effects on neutrino oscillations around black holes. The European Physical Journal C, 86(9), Article 1032. https://doi.org/10.1140/epjc/s10052-026-16278-8

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16278-8

Keywords: neutrino oscillations, dark matter, black holes, phantom field, gravitational lensing, general relativity, neutrino astronomy, decoherence, Sagittarius A*, M87*, particle astrophysics, spacetime curvature

The study, published in The European Physical Journal C by Adel Awad, Dmitry Chirkov, Alexey Golovnev and Alexey Toporensky, examines the Bianchi IX cosmology, the mathematically richest of the homogeneous but anisotropic universe models, whose spatial sections are three-dimensional spheres. The authors set out to explain, from first principles, a puzzling phenomenon that had previously been observed only in numerical simulations: the volume of such a universe does not monotonically decrease toward the singularity, but instead oscillates, repeatedly flipping between contraction and expansion. Their analysis reveals that this behavior is not a numerical artifact but a direct and inevitable consequence of the mathematics of phantom-dominated dynamics.

To appreciate the result, it helps to recall the standard BKL framework. In the vacuum Bianchi IX universe, the approach to the singularity can be described as a sequence of Bianchi I Kasner epochs, in which each scale factor grows or shrinks as a power of time, a ∝ t^p1, b ∝ t^p2, c ∝ t^p3, with the exponents satisfying p1 + p2 + p3 = 1 and p1² + p2² + p3² = 1. These constraints force at least one exponent to be negative but never allow any exponent to fall below −1/3. When one direction becomes dominant, curvature terms act like an exponential potential wall, and the universe bounces into a new Kasner epoch with new exponents. This discrete transformation, called the Mixmaster map, has been shown to be chaotic, and the BKL conjecture holds that this oscillatory regime is the generic behavior of spacetime near a singularity.

Matter, remarkably, barely matters in this regime. As the original BKL authors already understood, ordinary matter becomes negligible near the singularity, with one crucial exception: a stiff fluid with equation-of-state parameter w = 1, whose most natural representation is a massless scalar field. A standard massless scalar field actually kills the chaos entirely, replacing the oscillations with monotonic Kasner-like behavior. But the new work considers a massless scalar field with the wrong sign of its kinetic term, the so-called phantom field, whose energy density is negative. This single sign flip rewrites the rules of the game in profound ways.

The technical heart of the analysis lies in the constraint equation governing the Kasner exponents. With a phantom field, the constraint becomes p1² + p2² + p3² = 1 + q², where q is the conserved scalar momentum, a quantity protected by the shift symmetry of the massless field. Because the right-hand side is now larger than one, the exponents are no longer bounded. Two consequences follow immediately. First, two of the three Kasner indices can be negative simultaneously, meaning that two spatial dimensions grow while approaching the singularity in the past, a situation impossible in vacuum. Second, individual exponents can take arbitrarily large values, positive or negative, as q grows.

The possibility of two negative indices has a striking consequence for the Mixmaster map itself. In the standard picture, only one scale factor can trigger a bounce, so the exponents alone determine the next epoch. With two negative indices, either of the two corresponding dimensions can drive the bounce, and which one actually does depends on the actual sizes of the scale factors, information that is not encoded in the exponents at all. The authors show that when the two dimensions have very different sizes, the usual bounce formulas apply, but with the roles of the indices assigned according to which dimension dominates. When the two scale factors are nearly equal, the exponential potential walls disappear altogether, and the system undergoes a smooth, continuous evolution described by Bessel functions, an analog of the small oscillations near the Taub solution known from the vacuum case, but now accompanied by a reversal of the volume expansion rate.

The most dramatic finding concerns the volume of the universe. After each bounce, the rate of change of the volume in logarithmic time is proportional to 1 + 2p1, where p1 is the index of the bouncing dimension. In vacuum, p1 is always greater than −1/3, so this factor stays positive and the volume keeps shrinking toward the singularity. But with a phantom field, each bounce with a mildly negative index increases q², and once q² exceeds 3/8, the smallest index can drop below −1/2. At that point the denominator of the Mixmaster map changes sign, and the volume, followed backward in time, switches from contraction to expansion. The authors’ numerical integrations confirm this: after a finite sequence of bounces, the volume derivative flips sign, and the universe that was heading toward a singularity instead begins to re-expand.

This mechanism explains the volume oscillations observed in earlier numerical studies of phantom-dominated Bianchi IX cosmologies. The expansion phase cannot last forever, however. The Bianchi IX universe has the topology of a sphere, and its positive isotropic curvature falls off only as the inverse square of the mean scale factor, much more slowly than the anisotropies and the phantom energy, which both scale as the inverse sixth power. At large volumes, the curvature inevitably dominates and forces recollapse, a feature already known for vacuum Bianchi IX from the work of Lin and Wald. The result is a universe whose volume oscillates between finite bounds, bouncing off a curvature-driven maximum at large sizes and a phantom-driven minimum at small ones.

The oscillating volume has a deep implication for the BKL conjecture itself. Because volume maxima are unavoidable and the Mixmaster approximation breaks down near them, the map cannot be iterated infinitely. In the phantom case, the isotropic curvature term, though negligible near the singularity, always retains a finite contribution as the volume never approaches zero. The authors therefore conclude that the chaotic Mixmaster regime, rather than being a global attractor of the dynamics, becomes a transient phenomenon realized only with finite accuracy over a finite number of bounces. This challenges the expectation, still lacking a rigorous proof even in the vacuum case, that BKL oscillations describe the generic approach to a cosmological singularity.

The work also uncovers a genuinely new exact solution of the Bianchi I equations, possible only with phantom matter: an exponentially expanding or contracting universe of constant volume rate, in which the negative phantom energy exactly balances the anisotropic kinetic energy. Formally, such solutions correspond to infinite Kasner exponents, illuminating why the Mixmaster map can produce arbitrarily large indices. While a strict mathematical proof that phantom-filled Bianchi IX universes never hit a singularity remains to be constructed, the connection between the condition p1 < −1/2 and volume bounces makes that conjecture eminently reasonable. If phantom fields, which appear naturally in dark energy models and modified gravity theories such as Horndeski gravity, played a role in the early universe, the violent chaotic singularity foreseen by classical general relativity may be replaced by an eternal cosmic heartbeat, a universe that forever oscillates between expansion and collapse without ever reaching an end.

Subject of Research: Bianchi IX cosmological dynamics with a phantom scalar field and its effect on BKL oscillations near the singularity

Article Title: Bianchi IX dynamics with a phantom field

Article References: Awad, A., Chirkov, D., Golovnev, A., & Toporensky, A. (2026). Bianchi IX dynamics with a phantom field. The European Physical Journal C, 86(9), Article 1106. https://doi.org/10.1140/epjc/s10052-026-16386-5

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16386-5

Keywords: Bianchi IX, BKL oscillations, phantom field, Mixmaster map, Kasner exponents, cosmological singularity, anisotropic cosmology, general relativity, scalar field, volume oscillations, dark energy, theoretical physics