1. 1. Introduction
  2. 2. Neutrinos & Neutrino Oscillations
    1. 2.1. Neutrino States
    2. 2.2. Neutrino Oscillations
    3. 2.3. Neutrino Masses
    4. 2.4. Measuring Neutrino Oscillation Parameters
  3. 3. Neutrino Interactions
  4. 4. Coherent Pion Production
  5. 5. The T2K Experiment
  6. 6. A Search for Coherent Pion Production at T2K
  7. 7. Conclusions

1. Introduction

2. Neutrinos & Neutrino Oscillations

The existence of a light, weakly-interacting neutral particle was first suggested by Wolfgang Pauli in 1930 to explain the energy spectrum observed for electrons/positrons emitted in beta decays. Although not universally popular at first (even Pauli himself was sceptical) the idea eventually became widely accepted, but it was not until 1956 that the first experimental detection of “neutrinos” was made [1].

In the decades which followed discovery, their use in experiments was as probes of other physics. For example, the study of solar neutrinos was mainly motivated by understanding the fusion processes present in the Sun. Likewise, accelerator-based neutrino beams were made to study the workings of the weak force at lower energies than would be possible with electromagnetically or strongly interacting particles.

It wasn’t until observation of the solar and atmospheric neutrino problems, and the eventual discovery of neutrino oscillations, that neutrinos became of interest in their own right.

2.1. Neutrino States

The leptons, like all matter particles in the Standard Model, exhibit the peculiar feature that the states in which they undergo charged-current weak interactions, “flavour states”, are not mass eigenstates 1. Instead the three flavour states are a superposition, or mixture, of the three mass states [2].

Experimentally we have the freedom to define the reference point from which the flavour states are measured. Conventionally in quarks, we define the flavour states in terms of the up-type quarks, such that all the mixing occurs in the down-type quarks. Similarly in neutrinos, the convention is to define the flavour states in terms of the charged leptons with all the mixing confined to the neutrinos. There is good reason for this choice: the fact that neutrinos interact only weakly and that the differences between their masses are tiny, make it impossible to measure which mass state has been created.

One will sometimes find one or other of the flavour/mass states being described as the “physical” states (usually mass in quarks and flavour in neutrinos) but such a description is misleading. It is the measurement being made which determines which states are important, and any perception that one set is more “real” than another is simply an artefact of the view the experimenter has.

It is natural to assume that like the other fermions each neutrino flavour να has an anti-matter equivalent να, i.e. that they are Dirac particles. However because of their neutral charge it is possible that neutrinos could be Majorana particles. If this were to be the case then neutrinos would be their own anti-particles and they could have both Majorana and Dirac masses.

Today it is thought there are three neutrino flavour states which are defined by the charged lepton they interact with at charged-current weak vertices: νe, νμ and ντ. There are also three mass states ν1, ν2, ν3 with masses m1, m2 and m3 respectively:

The mass and flavour states are related by a mixing matrix:

νe νμ ντ = Ue1 Ue2 Ue3 Uμ1 Uμ2 Uμ3 Uτ1 Uτ2 Uτ3 ν1 ν2 ν3

This unitary matrix is known as the “Pontecorvo-Maki-Nakagawa-Sakata” matrix [3] [4]. The expression can also be inverted to express the mass states, i, in terms of the flavour states, α, highlighting the fact that they are both mixtures of each other:

νi = Σα U*αi να

It is convenient for experimental purposes to parametrise the PMNS matrix as the product of four sub-matrices [2]:

UPMNS = 100 0 cosθ23 sinθ23 0 -sinθ23 cosθ23 cosθ13 0 sinθ13e-iδ 010 -sinθ13e-iδ 0 cosθ13 cosθ12 sinθ12 0 -sinθ12 cosθ12 0 001 e½iα1 0 0 0 e½iα2 0 0 0 1

This parametrisation consists of three mixing angles (θ12, θ23 and θ13) two Majorana CP-violating phases (α1 and α2) and one Dirac CP-violating phase (δ).

The value of the three mixing angles effectively determines the degree to which the mass and flavour states are mixed.

The phases α1 and α2 are only physically observable if neutrinos are Majorana particles. Even if they exist it will be very difficult to measure their values experimentally, but in either case they do not contribute at all to neutrino oscillations (so they will be neglected from here on in).

However the final phase δ, which is commonly referred to as δcp, is directly observable and, if its value is non-zero, represents the source of leptonic CP-violation which most strongly motivates our interest in oscillations. As can be seen in Equation 2.3 however, it only appears with terms of sinθ13 so, regardless of the value of δcp, we also require θ13 ≠ 0 for any CP-violating effects to exist.

2.2. Neutrino Oscillations

Neutrino oscillations is a phenomenon which occurs as a result of the mixing of neutrino flavour and mass states discussed above. Experimentally it is the observation of a change in the flavour composition of neutrinos as they travel away from their source.

When a neutrino is created, it is created by the weak force, and therefore exists in a weak (flavour) eigenstate. So the wavefunction describing a pure source of να at creation looks like:

ψ 0,0 = να = Uα1 ν1 + Uα2 ν2 + Uα3 ν3

Later, after travelling a distance L through a vacuum the wavefunction becomes:

ψ L,t = Uα1 ν1 e-iφ1 + Uα2 ν2 e-iφ2 + Uα3 ν3 e-iφ3

where, assuming a plane-wavefunction 2, each mass state has acquired a phase φi = Eit - |pi|L. Since each of these states represents a different mass, the phases are also slightly different. So as the neutrinos travel away from the source, the mass states get out of phase with one another and the resulting interference means that the wavefunction evolves to contain components from all three flavour states. This can be seen from Equation 2.5 by substituting the mass states for their flavour mixtures from Equation 2.2:

ψ L,t = Uα1 U*e1 e-iφ1 + Uα2 U*e2 e-iφ2 + Uα3 U*e3 e-iφ3 νe + Uα1 U*μ1 e-iφ1 + Uα2 U*μ2 e-iφ2 + Uα3 U*μ3 e-iφ3 νμ + Uα1 U*τ1 e-iφ1 + Uα2 U*τ2 e-iφ2 + Uα3 U*τ3 e-iφ3 ντ

After travelling a distance L through a vacuum the neutrino wavefunction, which was initially pure να, is now a mixture of all three neutrino types and the proportion of each flavour varies as the neutrinos travel. This is neutrino oscillations.

Expanding Equation 2.6 out and calculating an example probability gets cumbersome so we make the assumption, valid at the resolution of current experiments, that δcp is zero to simplify things. Then for a beam of initially pure νμ with energy E, the probability of a neutrino being a νe at a distance L from the source is:

P νμ νe = sin2 θ23 sin2 13 sin2 Δm232 Φ + cos2 θ23 sin2 13 sin2 Δm221 Φ + cosθ13 sin12 sin13 sin23 sin Δm232 Φ sin Δm221 Φ cos Δm232 Φ

Where, for the sake of brevity:

Φ = L4E

and we meet for the first time the mass-splittings:

Δm2ij = mi2 - mj2

We can plot this probability as a function of L/E (Figure 2.1), making the origin of the name “oscillations” more apparent.

Figure 2.1: The detection probability of the three neutrino flavours from an initially pure νμ beam, as a function of L/Eν (using values from Table 2.1). Over longer distances, or at lower energies, “beat” patterns emerge - the neutrinos essentially behave as three coupled oscillators. The fact that the νμ component can reach both 1 and 0 is a consequence of θ23 = 45 °.

Oscillations are therefore a function of the three mixing angles, the distance travelled (often called the baseline) and the neutrino energy. They are also a function of the size of the differences between the mass states, Δm2ij, which is unsurprising given that it is the differences in how the states propagate, and hence their phases, that cause the oscillations. There are three of these mass splittings but since:

Δm232 + Δm221 = Δm231

there are only two which are independent. It is also worth emphasizing that it is only the differences between the masses which matter and neutrino oscillations are insensitive to the absolute scale of the three masses.

The Particle Data Group [2] provide a review of current world knowledge of the values of all the fundamental parameters that contribute to neutrino oscillations, which are summarised in Table 2.1.

Parameter Value
θ12 33.9 ± 1.0 °
θ23 39 ° < θ23 < 51 °
θ13 9.1 ± 0.6 °
Δm221 (7.50 ± 0.20) ×10-5  eV2
|Δm232| (2.32 +0.12-0.08) ×10-3  eV2
δcp unknown
Table 2.1: Current world knowledge of the neutrino oscillation parameters [2].

The equations above all apply to oscillations in a vacuum, however the situation is complicated further when a neutrino beam passes through matter. Because matter contains electrons but none of the other charged leptons/anti-leptons, νe have an additional interaction mode (CC forward scattering) available to them which is unavailable to the other neutrino flavours. This gives a larger effective mass to the νe component, which in turn affects the propagation of the associated states and hence the oscillations. The magnitude of the effect increases with the distance travelled in matter and the electron density encountered. Matter effects are not significant in the T2K experiment so won't be discussed in any more depth, but they are important in longer baseline experiments.

Finally, it is possible that there could be more than three neutrino flavour states, though the additional states would have to be “sterile”, i.e. not interact via the weak force, or it would conflict with other experimental evidence (the decay width of the Z0 [6], for example). If this were the case the mixing matrix would have to grow to accommodate it, and there may or may not be more neutrino mass states as well. However at present there is no convincing or consistent evidence to suggest that steriles exist [2].

2.3. Neutrino Masses

In the Standard Model neutrinos have always been treated as massless particles and it is only the discovery of neutrino oscillations which has given us any evidence to suggest otherwise.

The fact that any oscillations occur means that one of the Δm2ij ≠ 0, implying that at least one of the neutrino masses must also be non-zero. In fact, since all of the mass splittings are measured to be non-zero at least two of the masses are required to be non-zero. However, as discussed in Section 2.2, oscillations are insensitive to the absolute values of the neutrino masses and can only observe the differences between them.

At present the masses have proven to be so small that they are inaccessible to current experiments. The end-point of tritium beta-decay in the Troitzk experiment sets the world's most stringent limit: mνe < 2.05 eV at 95% CL [7]. Though it should be noted that the mass measured here, mνe, is that of a flavour state, and therefore represents a combination of the three mass states it contains.

Even with the mass differences we can measure, thus far it has not been possible to determine the sign of Δm231. The two possibilities are referred to as the “normal” (Δm231 > 0) and “inverted” (Δm231 < 0) hierarchies, and amount to a choice between which of the mass states is largest:

Normal Hierarchy: m3 > m2 > m1 Inverted Hierarchy: m2 > m1 > m3

The difficulty in determining the sign is essentially caused by the uncertainties on Δm232 and Δm231 being larger than the size of Δm221, and therefore the sign has a negligible effect on the oscillation probabilities compared to current experimental uncertainties. It is only matter effects in solar neutrino oscillations that have allowed us to determine the sign of Δm221.

2.4. Measuring Neutrino Oscillation Parameters

Equation 2.7 tells us that neutrino oscillations depend on the mixing angles, the mass splittings, the distance travelled and the neutrino energy. The mixing angles and mass splittings are properties of the Universe, we can't change them, we can only measure them. The distance travelled in a conventional accelerator-based oscillation experiment is a fixed quantity which we know.

The final parameter, the neutrino energy, is a property of the neutrino source - but that's an oversimplification. There is no practical way to create a mono-energetic neutrino source 3 and in practice experiments produce neutrinos over a range of energies.

The values of the oscillation parameters then, are determined by measuring the Eν spectrum for a given flavour before and after traversing the baseline. Any deficits/excesses found are the result of oscillations and the energy, size and shape of those changes can be used to measure the oscillation parameters. We can see this by considering the example of an initially pure νμ beam, measured after traversing a baseline of 295 km - approximately equivalent to the T2K experiment (Chapter 5) 4.

Determination of the mixing angles comes from measurement of the amplitude of oscillations (Figure 2.2). For example, the value of θ23 can be determined from the νμ spectrum, where the depth of the minimum resulting from oscillations away from νμ indicates the size of θ23. Likewise, measuring the height of the peak in oscillations to νe allows the determination of θ13.

Figure 2.2: (left) The effect of θ23 on P(νμ → νμ), its value can be determined by measuring the depth of the minimum. (right) The effect of θ13 on P(νμ → νe), in this case the value can be determined from the height of the peak.

The mass splittings on the other hand affect the energy at which these oscillation features occur (Figure 2.3). The value of Δm232 can be measured from the Eν of the minimum in the νμ spectrum, and the choice of normal or inverted mass hierarchy results in a shift in the Eν of the peak in the νe spectrum.

Figure 2.3: (left) The effect of Δm232 on P(νμ → νμ), the value can be determined by measuring the energy of the minimum. (right) The dominant effect on the νe spectrum of changing the mass hierarchy is a shift of the peak's energy.

The long-term goal is to determine whether or not neutrinos violate CP, and perhaps therefore contribute to the matter – anti-matter asymmetry of the Universe. In other words, to measure the value of δcp.

While the other oscillation parameters tend to affect one or other of the amplitude or energy of oscillations, δcp affects both (Figure 2.4). As the value of δcp is rotated, the position of the νe peak varies both in height and energy, traversing a circular path which reaches extremes in amplitude at δcp = 90/270 ° and energy at δcp = 0/180 °. It also affects νe and νe differently, with νe traversing the same circular path but in the opposite direction. Looking at lower Eν (Figure 2.5), equivalent to higher L/Eν, there are additional oscillation peaks which are much more strongly affected by the value of δcp.

Figure 2.4: The effect of δcp on P(νμ → νe). (left) As you rotate the angle of δcp, the position of the peak follows a circular path, changing both its height and energy. (right) The peak for P(νμνe) follows the same circular path, but in the opposite direction.

This leaves three ways through which such an experiment could determine δcp: a precise measurement of the position of the first peak in P(νμ → νe) spectrum, a comparison of that peak in neutrino and anti-neutrinos or, if the experiment is capable of resolving it, measuring the position of the second P(νμ → νe) peak.

Figure 2.5: The effect of δcp on P(νμ → νe). If an experiment can resolve the second oscillation maximum at lower Eν, the effects of δcp can be much larger.

However, determining the value of δcp through any of these methods will be particularly challenging. Because they only interact weakly, experiments cannot observe neutrinos directly. Instead, they only observe them via the secondary particles produced when they interact. In order to determine the proportion of neutrinos of a given flavour, the experiment must measure the number of events in which the corresponding charged lepton is produced, then use knowledge of the probability of that neutrino interacting to infer the number of neutrinos that actually passed through the detector. In order to measure the energy of the neutrino, the experiment must infer it from measurements of the secondary particles produced, and an understanding of how the neutrino's interaction produced these particles.

Not understanding the probability of a given neutrino flavour interacting will lead to incorrectly determining the number of those neutrinos present, and hence the oscillation probability. Not understanding the kinematics of secondary particles produced by a neutrino's interaction, will result in the energy being mis-reconstructed, and hence the neutrino energy spectrum being distorted. So in order to make precise measurements of the oscillation parameters, a precise understanding of neutrino interactions is needed.

This is particularly true for the determination of δcp. The variations in the height and position of the νe peak resulting from changes in δcp could also be attributed to changes in the values of the mixing angles and mass-splittings. And the size of those variations resulting from δcp are much smaller in scale than the current precision in those parameters will allow us to resolve. Furthermore, neutrinos and anti-neutrinos have quite different interaction probabilities, and result in a different range of secondary particles.

So our ability to measure δcp is limited by the need for precise knowledge of the other oscillation parameters, and of the interactions of neutrinos and anti-neutrinos. The precision of our knowledge of the other oscillation parameters, is also limited by our understanding of neutrino interactions, and how they relate to the energy and number of neutrinos in our experiment.

In order to increase the precision to which we understand neutrino oscillations, measure δcp, and perhaps solve the matter – anti-matter asymmetry of the Universe, we must first understand neutrino interactions.

3. Neutrino Interactions

4. Coherent Pion Production

5. The T2K Experiment

6. A Search for Coherent Pion Production at T2K

7. Conclusions


  1. This is unlike the electromagnetic and strong interactions, whose interaction states are mass eigenstates.
  2. Assuming plane-wavefunctions gives the same results as a more complete treatments in all but extreme cases [5].
  3. There are options for creating mono-energetic neutrinos but none of them are practical. Pion decay at rest would provide a mono-energetic νμ source, but their ∼30 MeV energy would be too low for νμ CC interactions. Alternatively, Z0 decay at rest would make an interesting neutrino source. Around 20% would decay to give an even mix of all six neutrinos and anti-neutrinos at 45.6 GeV. However the oscillation baseline required for such a high energy would be difficult, and the problems in making a high-flux source from stationary Z0s are many and large.
  4. Unless specified otherwise, for all the plots in this chapter the oscillation parameters are set to the values in Table 2.1, the initial flux is 100% νμ and L = 295 km.


  1. C.L. Cowan, F. Reines et al.
    Detection of the Free Neutrino: A Confirmation
    Science, Volume 124, p103-104, 1956
  2. The Particle Data Group
    Review of Particle Physics - Neutrino Mass, Mixing and Oscillations
    Physical Review D, Volume 86, p177-198, 2012
  3. B. Pontecorvo
    Neutrino Experiments and the Problem of Conservation of Leptonic Charge
    Journal of Experimental and Theoretical Physics, Volume 26, p984, 1968
  4. Z. Maki, M. Nakagawa, and S. Sakata
    Remarks on the Unified Model of Elementary Particles
    Progress of Theoretical Physics, Volume 28, p870-880, 1962
  5. B. Kayser
    On the quantum mechanics of neutrino oscillation
    Physical Review D, Volume 24, p110-116, 1981
  6. The ALEPH, DELPHI, L3, OPAL and SLD Collaborations, the LEP Electroweak Working Group, and The SLD Electroweak and Heavy Flavour Groups
    Precision electroweak measurements on the Z resonance
    Physics Reviews, Volume 427, p257-454, 2006
  7. V.N. Aseev et al.
    Upper limit on the electron antineutrino mass from the Troitsk experiment
    Physical Review D, Volume 84, p112003, 2011