Solar neutrino flavor conversion in the Sun is described by the Mikheyev–Smirnov–Wolfenstein (MSW) effect, and for the standard smooth, spherically symmetric solar density profile the adiabatic approximation is known to be excellent. The framework paper by Han et al. (2026) flags as unresolved the quantitative assessment of how three-dimensional, time-dependent density perturbations in the solar interior — such as those associated with solar gravity modes — modify neutrino propagation through nonadiabatic transitions and parametric-resonance-type responses, and what time variation they induce in the predicted solar neutrino flux at the precision level of current and next-generation detectors. The production-side channel (g-mode temperature fluctuations modulating neutrino genesis) has recently been quantified and is undetectable at present, while the propagation-side channel has only been analyzed in pre-modern two-flavor, one-dimensional idealizations; the full three-flavor treatment with realistic 3D time-dependent perturbations and current oscillation parameters remains quantitatively unassessed.

Background

Solar neutrinos are electron neutrinos produced by nuclear fusion reactions in the Sun's core (the pp, pep, 7^7Be, 8^8B, and hep reactions of the proton–proton chain, together with the carbon–nitrogen–oxygen (CNO) cycle), and they reach Earth as flavor mixtures because of neutrino oscillations, enhanced in matter by the Mikheyev–Smirnov–Wolfenstein (MSW) effect: coherent forward scattering off solar electrons adds an effective potential to the neutrino Hamiltonian, modifying the effective mass-squared splitting and mixing angle during propagation. Because the solar electron density decreases smoothly and monotonically outward, high-energy neutrinos (such as 8^8B neutrinos above roughly 5 MeV) undergo adiabatic flavor conversion — the propagation eigenstate amplitudes follow the instantaneous matter eigenstates, and the standard treatment reduces to evaluating the survival probability PeeP_{ee} along a static, spherically symmetric density profile; Han et al. (2026) verified numerically that a phase-incoherent adiabatic treatment misses the full three-flavor nighttime survival probability only at the O(10−5)\mathcal{O}(10^{-5}) level for standard 8^8B neutrinos. Solar gravity modes (g-modes) are oscillation modes of the Sun's radiative interior restored by buoyancy rather than pressure, with periods of hours; they probe the energy-generating core but have never been unambiguously observed — the claimed detection of asymptotic g-modes by Fossat et al. (2017) remains contested. If g-modes of sufficient amplitude exist, they perturb the solar electron density three-dimensionally and periodically in time, and neutrinos propagating through such a medium need not follow the adiabatic picture: a neutrino traversing a spatially periodic, moving density modulation can undergo parametric resonance — a resonant amplification of flavor conversion driven by the periodic variation of the matter potential along the path, conceptually distinct from the static MSW resonance. Neutrino propagation through stochastic or periodic solar density fluctuations was studied around 1996–2003 (Balantekin–Fetter–Loreti; Bamert–Burgess–Michaud; Burgess et al.) in the two-flavor framework with then-current parameters, with the general conclusion that small fluctuations negligibly affect high-energy 8^8B neutrinos while low-energy components (7^7Be, pp) are more susceptible to spectral distortion. Independently, g-modes modulate neutrino production through core temperature fluctuations; Lopes & Turck-Chièze (2014) showed this sensitivity is large (amplification factor of order 170 for 8^8B) and bounded the core temperature variation, and Hatta et al. (2026) recently computed that the time-varying part of this production-side modulation is of relative order 10−910^{-9}, far below detectability. Han et al. (2026) built a high-precision numerical framework (a full 3×33\times3 flavor propagator based on Strang splitting, a standard second-order operator-splitting scheme that evolves the Hamiltonian pieces in small alternating steps, together with three-dimensional Earth electron-density models and site-specific day–night flux predictions for underground laboratories) and point out that the adiabatic approximation, which discards path-dependent quantum phases, cannot capture nonadiabatic transitions or parametric responses induced by three-dimensional, time-dependent density perturbations such as g-modes. As solar neutrino experiments enter the sub-percent precision era — the Super-Kamiokande and Hyper-Kamiokande water-Cherenkov detectors in Japan, the Jiangmen Underground Neutrino Observatory (JUNO) liquid-scintillator detector in China, and the Deep Underground Neutrino Experiment (DUNE) in the United States — the question is whether such perturbations can measurably alter the solar neutrino flux and its time dependence, or whether they are guaranteed to be negligible.

Problem Statement

Quantitatively determine how three-dimensional, time-dependent electron-density perturbations in the solar interior — such as those associated with solar gravity modes — affect solar neutrino flavor conversion and the observable flux, by propagating neutrinos through the perturbed solar medium with the full three-flavor (coherent, phase-preserving 3×33\times3) treatment rather than the phase-incoherent adiabatic approximation. Specifically: (i) map physically motivated g-mode (or comparable helioseismic) perturbation fields — mode spectra, geometries, and amplitudes consistent with observational bounds — onto time-dependent electron-density perturbations δNe(r,t)\delta N_e(\mathbf{r},t) (with NeN_e the electron number density) along neutrino trajectories from production points in the core to the solar surface; (ii) compute the induced corrections to the survival probability and their dependence on neutrino energy, mode degree and radial order, amplitude, and oscillation phase, identifying whether nonadiabatic transitions or parametric-resonance-type enhancement can occur for any physically admissible perturbation or whether such enhancement is absent throughout the observationally allowed parameter range; (iii) propagate these probability-level effects through production-point, trajectory, and time averaging to obtain the resulting time variation (and any time-average shift) of the predicted solar neutrino fluxes for the relevant components, including 8^8B and the low-energy 7^7Be/pp components where earlier work indicated the greatest susceptibility; and (iv) state whether the induced signals are in principle observable at current or next-generation detectors (Super-Kamiokande/Hyper-Kamiokande, successors to the Borexino low-energy liquid-scintillator detector, JUNO, DUNE) or constrained by existing data. A complete solution must address all four parts: even a negative determination in (ii) requires the corresponding bound on the resulting flux time variation in (iii) and the comparison with detector sensitivity in (iv).

References

  1. Keyu Han, Isabella Yin, Kaoru Yagi, Kevin Yifan Jiang, Xiangpan Ji, Shaomin Chen (2026), 'A High-Precision Numerical Framework for Time-Varying Solar Neutrino Flux with Full Earth Matter Oscillation Corrections for Global Underground Laboratories', arXiv:2607.24421, DOI 10.48550/arxiv.2607.24421, https://arxiv.org/abs/2607.24421
  2. Ilídio Lopes, Sylvaine Turck-Chièze (2014), 'Detecting gravity modes in the solar 8^8B neutrino flux', The Astrophysical Journal Letters 792, L35, DOI 10.1088/2041-8205/792/2/L35, arXiv:1408.6671, https://arxiv.org/abs/1408.6671
  3. Yoshiki Hatta, Yuuki Nakano, Sho Sugama, Masanobu Kunitomo, Hiroshi Ito, Takashi Sekii (2026), 'Solar Neutrino Flux Fluctuations Caused by Solar Gravity Modes', arXiv:2604.06535, https://arxiv.org/abs/2604.06535
  4. A. B. Balantekin, J. M. Fetter, F. N. Loreti (1996), 'MSW effect in a fluctuating matter density', Physical Review D 54, 3941, DOI 10.1103/PhysRevD.54.3941, https://doi.org/10.1103/PhysRevD.54.3941
  5. P. Bamert, C. P. Burgess, D. Michaud (1998), 'Neutrino propagation through helioseismic waves', Nuclear Physics B 513, 319, DOI 10.1016/S0550-3213(97)00672-X, https://doi.org/10.1016/S0550-3213(97)00672-X
  6. C. Burgess, N. S. Dzhalilov, M. Maltoni, T. I. Rashba, V. B. Semikoz, M. A. Tórtola, J. W. F. Valle (2003), 'Large Mixing Angle Oscillations as a Probe of the Deep Solar Interior', The Astrophysical Journal 588, L65, DOI 10.1086/375482, arXiv:hep-ph/0209094, https://arxiv.org/abs/hep-ph/0209094
  7. A. B. Balantekin, H. Yüksel (2003), 'Do the KamLAND and solar neutrino data rule out solar density fluctuations?', Physical Review D 68, 013006, DOI 10.1103/PhysRevD.68.013006, arXiv:hep-ph/0303169, https://arxiv.org/abs/hep-ph/0303169
  8. E. Fossat, P. Boumier, T. Corbard, J. Provost, D. Salabert, F. X. Schmider, A. H. Gabriel, G. Grec, C. Renaud, J. M. Robillot, T. Roca-Cortés, S. Turck-Chièze, R. K. Ulrich, M. Lazrek (2017), 'Asymptotic g modes: Evidence for a rapid rotation of the solar core', Astronomy & Astrophysics 604, A40, DOI 10.1051/0004-6361/201730460, https://doi.org/10.1051/0004-6361/201730460
  9. S. J. A. J. Salmon, G. Buldgen, A. Noels, P. Eggenberger, R. Scuflaire, G. Meynet (2021), 'Standard solar models: Perspectives from updated solar neutrino fluxes and gravity-mode period spacing', Astronomy & Astrophysics 651, A106, DOI 10.1051/0004-6361/202140769, arXiv:2105.00911, https://arxiv.org/abs/2105.00911
  10. Liangliang Ren, Ze-Lin Zhang, Bing Xu, Tian-Cheng Huang, Ran Wang, Jia-Xin Dong, Yan-Ping Wang (2026), 'Comparative Periodogram Analysis of 22 Years of Super-Kamiokande Solar 8^8B Neutrino Data: Classical, Phase-Based, and Information Theoretic Methods', arXiv:2607.27979, https://arxiv.org/abs/2607.27979