ν1 ν2 ν3 νe νμ ντ SOURCE ILLUSTRATIVE PROPAGATION DISTANCE →

Research

Neutrino phenomenology

Neutrino phenomenology · Particle physics beyond the Standard Model · Astroparticle physics · Cosmology

Research statement

Overview

I have been a researcher in theoretical physics from 2005 to 2015 and associate professor until 2023. I am currently Full Professor of Theoretical Physics at the University of Bari "Aldo Moro".

My research lies at the interface between neutrino physics, particle phenomenology and cosmology. With the Bari group I have contributed to a long line of global analyses of neutrino oscillation data, combining solar, atmospheric, reactor and accelerator experiments to extract the fundamental parameters of leptonic mixing with state-of-the-art precision.

As measurements enter the subpercent regime, these results sharpen the picture of the neutrino mass ordering, of leptonic CP violation and of the absolute neutrino mass scale, while probing possible signatures of physics beyond the three-neutrino paradigm.

At a glance

Where this work sits

Full Professor of Theoretical Physics, Università degli Studi di Bari “Aldo Moro” and INFN Sezione di Bari. Researcher from 2005, associate professor to 2023.

Research theme 01

Global analysis of neutrino oscillations

Three-flavour fits to the world neutrino data, with a self-consistent treatment of solar, atmospheric, reactor and long-baseline accelerator experiments.

Question
What are the values of the three mixing angles, the two mass-squared splittings and the CP phase, and how precisely do the world data fix them?
Method
One self-consistent three-flavour fit of solar, atmospheric, reactor and long-baseline accelerator data, with the full matter-effect treatment and a careful statistical combination of datasets with very different systematics.
Our part
The Bari group's long line of global analyses, released and updated as new data arrive.
Main result
In the March 2025 release |Δm²| is known to 0.8 %, the first three-neutrino parameter below 1 %. The November 2025 update, dominated by the first JUNO results, takes δm² from 2.3 % to 1.3 %.
Read
Go deeper
Analysis results · full tables on global-nu.org · the two figures below

The observables are the two mass-squared splittings $\Delta m^2_{21}$ and $|\Delta m^2_{3\ell}|$ and the three mixing angles $\theta_{12}$, $\theta_{13}$, $\theta_{23}$ of the PMNS matrix, together with the CP-violating phase $\delta_{\rm CP}$. In vacuum the survival probability of a flavour eigenstate takes the familiar form

$$ P(\nu_\alpha \to \nu_\alpha) = 1 - 4\sum_{i<j} |U_{\alpha i}|^2 |U_{\alpha j}|^2 \sin^2\!\left(\frac{\Delta m^2_{ij} L}{4E}\right), $$

but the precision now attainable requires the full matter-effect treatment and a careful statistical combination of datasets with very different systematics.

Solar Atmospheric Reactor Accelerator 3ν global fit δm² Δm² θ12 θ13 θ23 δCP DATA PARAMETERS ONE SELF-CONSISTENT ANALYSIS
Schematic. Four classes of experiment constrain overlapping combinations of the same six parameters, which is why they are fitted together rather than one at a time.

Current best fits

The six parameters, as we measure them

What the global analysis returns for the three mixing angles, the two mass-squared splittings and the CP phase, from the Bari Release 2026.a after the second JUNO data release — preliminary, the paper is in preparation.

Solar mixing angle
0.3047 3σ: 0.2888 – 0.3217
sin²θ₁₂ · 1.8 % precision (1σ)
Preliminary Bari Release 2026.a
Reactor mixing angle
0.0222 3σ: 0.0205 – 0.0238
sin²θ₁₃ · 2.5 % precision (1σ)
Preliminary Bari Release 2026.a
Atmospheric mixing angle
0.491 3σ: 0.444 – 0.573
sin²θ₂₃ · 4.4 % precision · both octants allowed up to 3σ
Preliminary Bari Release 2026.a
Solar splitting
7.465 ×10⁻⁵ eV² · 3σ: 7.239 – 7.681
δm² · 1.0 % precision (1σ)
Preliminary Bari Release 2026.a
Atmospheric splitting
2.472 ×10⁻³ eV² · 3σ: 2.432 – 2.517
|Δm²| · 0.5 % precision (1σ)
Preliminary Bari Release 2026.a
Leptonic CP phase
1.19 3σ: 0.79 – 1.95
δ/π · δ = 0 disfavoured at 3.2σ, δ = π at 1.4σ
Preliminary Bari Release 2026.a

Preliminary. Bari Release 2026.a, after the second JUNO data release: the paper is in preparation and these values may still change. Final combination of long-baseline accelerator, solar, reactor and atmospheric data; normal ordering, favoured over inverted at 3.9σ (Δχ²IO−NO = 15.34). δm² = m₂² − m₁², Δm² = m₃² − (m₁² + m₂²)/2. The 1σ precision is the 1σ error divided by the best fit; for sin²θ₂₃, whose 3σ range spans both octants, one sixth of the 3σ range divided by the best fit.

Research theme 02

Mass ordering and leptonic CP violation

Question
Is the neutrino mass ordering normal or inverted, and is CP violated in the lepton sector?
Method
Statistical extraction of δCP and of the sign of Δm² from T2K, NOvA, reactor and atmospheric data together, and forecasts for JUNO, DUNE and Hyper-Kamiokande.
Our part
Combining datasets that pull in partly different directions, with their systematics treated jointly, since the significance depends delicately on it.
Main result
In the March 2025 release normal ordering is favoured at 2.2σ and, in normal ordering, CP violation at 1.3σ: both questions are still open.
Read

Statistical extraction of the leptonic CP phase $\delta_{\rm CP}$ and of the sign of $\Delta m^2_{3\ell}$ from the interplay of T2K, NOvA, reactor and atmospheric data. The two datasets pull in partly different directions, and the resulting significance depends delicately on how their systematics are combined.

A related strand is forecasting the discovery reach of the next generation of experiments — JUNO, DUNE and Hyper-Kamiokande — and identifying which combinations of measurements would resolve the ordering and the $\theta_{23}$ octant with the least model dependence.

Diagram comparing the normal and inverted neutrino mass orderings, with the flavour composition of each mass eigenstate shown as coloured fractions.
The neutrino mass spectrum, with the flavour composition $|U_{\alpha i}|^2$ of each eigenstate coded by colour. Whether $\nu_3$ sits at the top (normal) or at the bottom (inverted) is one of the open questions our global fits aim to resolve.

Research theme 03

Absolute neutrino masses and cosmology

Question
How heavy are neutrinos? Oscillations measure only differences of squared masses.
Method
A joint statistical treatment of the three probes of the absolute scale — the β-decay endpoint (mβ), neutrinoless double-β decay (mββ) and cosmology (Σ) — together with the oscillation results.
Our part
Allowed regions in the (mβ, mββ, Σ) space, and how the tension between cosmological bounds and laboratory limits evolves as data improve.
Main result
Global constraints on the absolute masses and on their ordering from all these data at once, and Majorana-mass constraints that carry the spread of nuclear matrix elements through.
Read

Oscillation experiments are sensitive only to mass-squared differences. The absolute scale is constrained instead by three complementary probes: the endpoint of the β-decay spectrum, which measures the effective mass $m_\beta$; searches for neutrinoless double-β decay, sensitive to the Majorana mass $\langle m_{\beta\beta} \rangle$; and cosmological observations of the CMB and of large-scale structure, which bound the sum $\Sigma = \sum_i m_i$.

Combining these with the oscillation results is a genuinely statistical problem, since the three observables depend on different combinations of the same underlying parameters and carry very different systematics. Our analyses map the allowed regions in the $(m_\beta,\ \langle m_{\beta\beta}\rangle,\ \Sigma)$ space, and quantify how the tension between cosmological bounds and laboratory limits evolves as data improve.

m1, m2, m3 β-decay endpoint mβ 0νββ decay ⟨mββ⟩ Cosmology Σ THE SAME THREE MASSES · THREE DIFFERENT COMBINATIONS
Schematic. Each probe measures a different combination of the same three masses, so the three bounds have to be combined statistically rather than compared directly.

Current bounds

Where the absolute scale stands

The three probes, as they constrain the scale today — one laboratory limit each from β decay and from neutrinoless double-β decay, and the cosmological sum.

Direct kinematic mass
< 0.45 eV
mβ · KATRIN, 90 % C.L.
Majorana mass
≲ 28–122 meV
⟨mββ⟩ · KamLAND-Zen, ¹³⁶Xe
Cosmological sum
Σ mν
CMB, large-scale structure and DESI

Research theme 04

Supernova neutrinos

Flavour evolution in the one place where neutrinos are dense enough to refract off each other.

Question
How do neutrinos change flavour near a collapsing core, where the flux is dense enough for neutrinos to refract off each other — and what will the next galactic supernova's events tell us?
Method
Non-linear, collective flavour evolution: multi-angle geometry instead of a single-angle approximation, three flavours together, and the dispersion relation of the flavour field.
Our part
With the Bari group: trajectory averaging, three-flavour collective evolution, spectral splits, and the low-energy features in inverted ordering.
Main result
Spectral-split patterns mapped against the relative luminosities of the species, and a classification of fast flavour conversions through dispersion relations.
Read

A core-collapse supernova releases almost all of its gravitational binding energy as neutrinos, and it does so before anything is visible at the surface. They leave the collapsing core while the shock is still buried, which makes them the only messengers that carry direct information about the proto-neutron star and about the mechanism that revives the explosion. SN 1987A remains the sole detection: a couple of dozen events, and still the observational basis of most of what can be said. The next galactic supernova will be recorded by detectors that did not exist then, and the question is what those events can be made to say.

What makes the problem hard — and what our work has been about — is that the neutrino flux near the core is dense enough that neutrinos refract off each other. The flavour evolution stops being a one-particle problem and becomes a non-linear, collective one: the state of the ensemble enters the Hamiltonian that evolves it. With the Bari group we studied how these collective flavour transitions behave once the geometry is treated honestly rather than in a single-angle approximation, and how the multi-angle structure of the emission changes the outcome; how the transitions look when all three flavours evolve together instead of two; and how the resulting spectral splits depend on the relative luminosities of the species. We also worked out the consequences for the low-energy end of the spectrum in inverted ordering, and, more recently, on classifying the fast flavour conversions — instabilities that develop on scales far shorter than the ones set by the mass splittings — through the dispersion relation of the flavour field.

Selected work

Papers on collective flavour transitions

Research theme 05

Beyond the three-neutrino paradigm

Tests of the standard 3ν framework against possible extensions.

Question
Is the three-neutrino picture complete, or do the data leave room for something more?
Method
Fits that add light sterile states (new splittings and mixing angles), non-standard interactions (a modified matter potential), and checks of the consistency among datasets.
Main result
The three-neutrino picture describes the data remarkably well; tensions among datasets are followed as statistics grow, since a persistent one may be the first sign of new physics rather than of underestimated systematics.
Read

The three-neutrino picture describes the data remarkably well, but it is not guaranteed to be complete. We test it against light sterile states, which would appear as additional mass-squared splittings and mixing angles; against non-standard neutrino interactions, which modify the matter potential and therefore the effective mixing in the Sun and in the Earth; and against tensions among datasets, which — if they persist as statistics grow — may be the first indication of new physics rather than of underestimated systematics.

THREE NEUTRINOS 3 + 1 ν1 ν2 ν3 ν1 ν2 ν3 ν4 mostly sterile MASS SCHEMATIC · NOT TO SCALE
Schematic, not to scale. A light sterile state would add a fourth mass eigenstate far above the three active ones, and with it new splittings and mixing angles for the data to constrain.

Research theme 06

Modular invariance and the flavour puzzle

Question
Is there a symmetry behind the masses and mixings of quarks and leptons?
Method
Modular-invariant flavour models, in which Yukawa couplings are modular forms of a complex modulus τ, confronted with data in a global fit.
Our part
The first joint fit of quark and lepton observables — 22 quantities with 14 real parameters, in a model based on the binary octahedral group 2O — and modular invariance as an answer to the strong CP problem.
Main result
Strong correlations, invisible in separate analyses, between quark mass ratios and the leptonic CP phase, which become sharp predictions for KATRIN, neutrinoless double-β decay, JUNO and DUNE.
Read

Following the proposal by F. Feruglio that the Yukawa couplings of leptons could be modular forms of a complex modulus $\tau$ living in the upper half-plane, modular invariance has emerged as a powerful candidate symmetry for the flavour structure of quarks and leptons.

With G.-J. Ding (USTC), E. Lisi (INFN Bari) and S. T. Petcov (SISSA/IPMU) we have carried out the first joint fit of quark and lepton observables — 22 quantities in total — within a modular flavour model based on the binary octahedral group $2O$, with just 14 real parameters. The fit reveals strong correlations among observables, notably between quark mass ratios and the leptonic CP phase $\delta_{\rm CP}$, that are invisible in separate analyses and that translate into sharp predictions for KATRIN, neutrinoless double-β decay, JUNO and DUNE.

In a complementary direction, with F. Feruglio (Padova), A. Strumia and A. Titov (Pisa), we have shown that modular invariance can also address the strong CP problem in string-inspired settings, where quarks have positive modular weights and gauge kinetic functions are non-trivial.

The upper half-plane with the shaded fundamental domain of the modular group SL(2,Z), bounded by the unit circle and the lines Re(tau) = ±1/2, with the fixed points i, omega and −omega-bar marked.
The complex modulus $\tau$ lives in the shaded fundamental domain $\mathcal{F}$ of $SL(2,\mathbb{Z})$. Special points $i$, $\omega$, $-\bar\omega$ and $i\infty$ are fixed by residual symmetries and play a distinguished role in flavour model building.

Interactive

Interactive figures

Eight figures that compute what they show, each on its own page with its assumptions, sources, a PNG for slides and a licence to reuse. The full collection, with filters by topic and level, is on the Interactive page.

Oscillations

How do three neutrino flavours oscillate?

Exact three-flavour oscillation probabilities in vacuum as a function of L/E, with the parameters of the Bari global analysis.

Advanced undergraduate: quantum mechanics and the idea of flavour mixing.

Open →
Detectors

What does a neutrino look like in a water-Cherenkov detector?

The neutrino itself is invisible: the ring of light comes from the charged lepton it produces.

Undergraduate: special relativity and the refractive index.

Open →
Cosmology

What happened to the neutrinos of the Big Bang?

A scrubbable history of the relic neutrinos, from the primordial plasma to the sea that fills every cubic centimetre of space today.

Graduate: the thermal history of the early Universe.

Open →
Supernovae

What would a galactic supernova look like in neutrinos?

The neutrino burst of a core collapse in time and energy, and what three detector technologies would record at a distance you choose.

Graduate: neutrino astrophysics.

Open →
Mass models

What must we add to the Standard Model to give neutrinos mass?

A map of the theoretical choices, from the minimal Standard Model to SO(10), with the fields, the mass matrix and the observable consequences of each.

Graduate: the Standard Model and basic quantum field theory.

Open →
Seesaw

How do three different seesaws become one operator?

The type I, II and III seesaw diagrams collapse onto the same dimension-five operator as the energy scale runs below the heavy mass.

Graduate: quantum field theory and effective field theory.

Open →
Absolute masses

What do β decay, 0νββ and cosmology each measure?

mβ, mββ and Σ as functions of the lightest mass, the ordering and the Majorana phases, set against the current bounds.

Advanced undergraduate to graduate.

Open →
Mass ordering

How can a reactor 52.5 km away tell the neutrino mass ordering?

The reactor antineutrino spectrum at JUNO in seven scenes: the fast ripple whose phase tells normal from inverted ordering, and what resolution and statistics leave of it.

Advanced undergraduate to graduate.

Open →

Selected recent publications

The complete list is maintained on INSPIRE-HEP.

Conferences

NOW — Neutrino Oscillation Workshop

Otranto, biennial. One of the longest-running European workshops in neutrino physics, gathering theorists and experimentalists for in-depth discussions on the frontiers of the field. I am among the organizers.

Opportunities

For prospective students and postdocs

I welcome motivated students and young researchers interested in neutrino phenomenology, statistical analysis of particle physics data, and the interplay between particle physics and cosmology.

A historical note

Pauli’s letter, 1930

The story of neutrinos began with one of the most famous letters in twentieth-century physics. On 4 December 1930 Wolfgang Pauli postulated a new neutral particle of tiny mass to rescue energy conservation in β-decay — a hypothesis he himself called “a desperate remedy”.

Dear Radioactive Ladies and Gentlemen, as the bearer of these lines will explain to you in more detail, because of the “wrong” statistics of the N and Li⁶ nuclei and the continuous beta spectrum, I have hit upon a desperate remedy to save the “exchange theorem” of statistics and the law of conservation of energy. Namely, the possibility that there could exist in the nuclei electrically neutral particles, that I wish to call neutrons, which have spin ½ and obey the exclusion principle and which further differ from light quanta in that they do not travel with the velocity of light. […]

But only the one who dare can win, and the difficult situation, due to the continuous structure of the beta spectrum, is lighted by a remark of my honoured predecessor, Mr Debye, who told me recently in Bruxelles: “Oh, it is best not to think about this at all, like new taxes”. […]

Unfortunately, I cannot appear in Tübingen personally since I am indispensable here in Zurich because of a ball on the night of 6/7 December. With my best regards to you, and also to Mr Back.

Your humble servant, W. Pauli.

Facsimile of the typewritten German original of Pauli's 1930 open letter to the Tübingen conference.
Facsimile of the German original, addressed to the group of radioactives at the Tübingen meeting.
Historical photograph of Wolfgang Pauli, Werner Heisenberg and Enrico Fermi together.
Pauli, Heisenberg and Fermi.