
2
are neutral elementary particles, and it is intriguing that
this particle could be a Majorana particle, where neu-
trino and anti-neutrino are indistinguishable in the con-
ventional sense of the CPT invariant paradigm. The neu-
trino nature under the CPT-violating scenario has been
explored in Ref. [12]. Here we focus on the phenomeno-
logical consequence of the CPT violation in the observ-
able neutrino oscillation.
In context of three-flavor PMNS framework [13,14], for
a given propagation distance Land matter density ρ, the
probabilities (Pνα→νβ, Pνα→νβ) for a neutrino and anti-
neutrino at a specific energy (Eν, Eν) oscillating from one
flavor (να, να) to another flavor (νβ, νβ) are completely
and commonly described with six oscillation parameters
including three leptonic mixing angles (θ12, θ13 , θ23),
one Dirac CP-violating phase δCP , and two mass-squared
differences (∆m2
21,∆m2
31). Under CPT symmetry, the
neutrino and anti-neutrino oscillation probabilities are
well connected as follows:
Pνα→νβ
CPT
−−−→Pνβ→να=Pνα→νβ
=f(θ12, θ13, θ23, δCP,∆m2
21,∆m2
31).
If the CPT is violated in neutrino sector, the under-
lying sets of oscillation parameters in neutrino and anti-
neutrino may differ. Empirically, we assume
Pνα→νβ=f(θ12, θ13, θ23, δCP,∆m2
21,∆m2
31),(3)
for describing the neutrino oscillations, and
Pνβ→να=f(θ12, θ13, θ23, δCP ,∆m2
21,∆m2
31),(4)
for anti-neutrino oscillations.
If there are observable differences in the parameters
of the two sets, it may indicate a CPT violation in the
lepton sector. Since the discovery of neutrino oscilla-
tions [15,16] at the end of the twentieth century, neu-
trino oscillation experiments [8] using both natural and
man-made neutrino sources have transitioned into the
precision measurement phase of three mixing angles and
two mass-squared differences, and being explored three
remained known unknowns including the neutrino mass
ordering, whether CP is violated, and whether the mix-
ing angle θ23 is maximal (θ23 = 45◦) or belong to a
lower (θ23 <45◦) or higher (θ23 >45◦) octant. Each
experiment is typically sensitive to a subset of the os-
cillation parameters but not the entire set. The ex-
periments with solar neutrinos provide the most con-
straints on the (θ12,∆m2
21) parameters while the reactor-
based long-baseline neutrino (R-LBL) experiments can
measure precisely the (θ12,∆m2
21) parameters. The
reactor-based short-baseline (order of 1 km) neutrino (R-
SBL) experiments play a central role in measuring the
(θ13,∆m2
31) parameters. The under-developing reactor-
based medium-baseline neutrino (R-MBL) experiment
JUNO, which will be discussed later, takes advantage of
interference of oscillations at different wavelengths, huge
statistics, and good energy resolution to achieve sub-
percent precision in measuring the (θ12,∆m2
21,∆m2
31)
parameters. Experiments with the atmospheric neu-
trino and accelerator-based neutrino sources can pre-
cisely measure the (θ23, θ23,∆m2
31,∆m2
31) parameters.
Besides, this type of experiment is also sensitive to the
(θ13,θ13) parameters, but the precision of these param-
eters is much lower in comparison to the R-SBL experi-
ment due to the statistical limit and their strong correla-
tion with two known unknowns, CP-violating phase and
neutrino mass ordering. Although there is some hint [17]
of non-zero CP-violating phase δCP , precise measure-
ment on this parameter is not possible until the next gen-
eration of the accelerator-based long-baseline (A-LBL)
experiments. It is provided in Ref. [18] the most recent
update at 3σconfidence level (C. L.) on the bounds of
CPT violation on each individual parameter with global
neutrino data.
|δνν (∆m2
21)|<4.7×10−5eV2,
|δνν (∆m2
31)|<2.5×10−4eV2,
|δνν (sin2θ12)|<0.14,
|δνν (sin2θ13)|<0.029,
|δνν (sin2θ23)|<0.19,(5)
where δνν (X) = X−Xfor the X neutrino oscillation
parameter and the Xanti-neutrino oscillation param-
eter. In this study, we focus on the synergy between
two on-going A-LBL experiments (T2K and NOνA) and
one under-developing R-MBL experiment (JUNO) to ex-
plore the potential sensitivity to the measurement of
δνν (∆m2
31) and δνν (sin2θ23 ) parameters. The A-LBL
experiments utilize the highly intense beam of the al-
most pure muon neutrinos νµand muon anti-neutrinos
νµfor measuring the four transitions categorized into two
channels, appearance channels (νµ→νe, νµ→νe), and
disappearance channels (νµ→νµ, νµ→νµ). While the
appearance channels are sensitive to a wider subset of pa-
rameters and being explored for searching the CP viola-
tion in the lepton sector, measuring (νµ→νe, νµ→νe)
is not sufficient to test CPT directly since the corre-
sponding CPT-mirrored processes are missing. The dis-
appearance channels, on the other hand, are well-suited
for testing CPT since they are two CPT-mirrored pro-
cesses. We characterize the difference in the probabili-
ties of the muon neutrino disappearance and muon anti-
neutrino disappearance,ACPT
µµ =Pνµ→νµ−Pνµ→νµas an
observable measure of the CPT-violating effect.
The observable asymmetry ACPT
µµ may consist of two
parts: intrinsic CPT asymmetry and extrinsic CPT
asymmetry caused by differences in interactions between
neutrinos and anti-neutrinos with the matter of the prop-
agation medium [19–23]. Fig. 1illustrates the CPT
asymmetries ACPT
µµ calculated in vacuum and in the mat-
ter presence at baselines of the T2K experiment (L =
295 km) and of the NOνA experiment (L = 810 km).
Here we take the best-fit values of the mainly involved
(∆m2
31,∆m2
31, θ23, θ23) parameters from the recent T2K
results [24] and of the others from the global data anal-
ysis [25], which are summarized in Table I. It is wor-