Heavy neutrino mixing in the T2HK, the T2HKK and an extension of the T2HK with a detector at Oki Islands

TL;DR

This paper investigates the potential to detect heavy neutrino mixing in extended T2K experiments, showing that adding a detector at Oki Islands improves sensitivity and allows for better measurement of neutrino properties.

Contribution

It introduces a new extension of T2HK with an Oki Islands detector and analyzes its improved sensitivity to heavy neutrino mixing and neutrino oscillation parameters.

Findings

Adding a detector at Oki Islands mitigates sensitivity deterioration.

The extended setup improves measurement accuracy of neutrino mass hierarchy.

Data with varying L/E ratios enhance parameter determination.

Abstract

We study the discovery potential for the mixing of heavy isospin-singlet neutrinos in extensions of the Tokai-to-Kamioka (T2K) experiment: The Tokai-to-Hyper-Kamiokande (T2HK), the Tokai-to-Hyper-Kamiokande-to-Korea (T2HKK), and a plan of adding a new detector at Oki Islands to the T2HK. We parametrize the mixing of heavy neutrinos in terms of a non-unitary mixing matrix acting on active flavors. It is shown that in the T2HK and T2HKK, the sensitivity to heavy neutrino mixing deteriorates for some values of $CP$-violating phases in the standard and the non-unitary mixing matrices, but the deterioration is drastically mitigated if a detector at Oki Islands is added to the T2HK. We also consider the feasibility of measuring the mass hierarchy and the standard $CP$-violating phase $\delta_{CP}$ in the presence of heavy neutrino mixing, by fitting experimental data with the standard oscillation parameters only, i.e., under the assumption of unitary mixing. It is revealed that such measurement can be performed to some accuracy in the T2HKK and the extension of the T2HK with a detector at Oki Islands, owing to the fact that data with largely varying $L/E$ are collected in these experiments, while it cannot be in the T2HK.

Tables

Table 1. Table 1: Parameters used in our simulation.

physical parameter	value in our simulation
${\sin }^2{\theta }_{12}$	0.304
${\sin }^2{\theta }_{13}$	0.0219
${\sin }^2{\theta }_{23}$ (normal hierarchy)	0.51
${\sin }^2{\theta }_{23}$ (inverted hierarchy)	0.50
$\Delta m_{21}^2$	$7.53\times {10}^{-5}$ eV2
$\left|\Delta m_{32}^2\right|$ (normal hierarchy)	$2.44\times {10}^{-3}$
$\left|\Delta m_{32}^2\right|$ (inverted hierarchy)	$2.51\times {10}^{-3}$
${\delta }_{CP}$	free
${\alpha }_{11}$	0.992573
$\left|{\alpha }_{21}\right|$	0.003494
${\alpha }_{22}$	0.999589
${\alpha }_{31}$	0
${\alpha }_{32}$	0
${\alpha }_{33}$	1
${\phi }_{21}=\arg \left({\alpha }_{21}\right)$	free

Table 2. Table 2: Parameters assumed for the detectors used in the experiments.

experiment	site	baseline length $L$	fiducial volume	matter density $\rho$	off-axis angle
T2HK	Kamioka	295 km	374 kton	2.60 g/cm3	2.5∘
T2HKK	Kamioka	295 km	187 kton	2.60 g/cm3	2.5∘
	Korea	1000 km	187 kton	2.90 g/cm3	1.0∘
T2HK+Oki	Kamioka	295 km	374 kton	2.60 g/cm3	2.5∘
	Oki	653 km	10 kton	2.75 g/cm3	1.0∘
T2HK+Kami	Kamoika	295 km	384 kton	2.60 g/cm3	2.5∘

Table 3. Table 3: Flux of neutrino-focusing beam and antineutrino-focusing beam from J-PARC assumed.

	Neutrino-focusing	Antineutrino-focusing
Flux	$27\times {10}^{21}$ POT	$27\times {10}^{21}$ POT

Table 4. Table 4: The total number of neutrino events at each detector, with the detector size of Table  2 and the beam flux of Table  3 . The numbers for the neutrino-focusing beam and for the antineutrino-focusing beam are shown separately.

experiment	detector	number for	number for
		neutrino-focusing beam	antineutrino-focusing beam
T2HK	Kamioka	21$\times {10}^4$	7.5$\times {10}^4$
T2HKK	Kamioka	11$\times {10}^4$	3.7$\times {10}^4$
	Korea	3.1$\times {10}^4$	1.7$\times {10}^4$
T2HK+Oki 10 kton	Kamioka	21$\times {10}^4$	7.5$\times {10}^4$
	Oki	0.39$\times {10}^4$	0.21$\times {10}^4$
T2HK+Kami 10 kton	Kamioka	22$\times {10}^4$	7.7$\times {10}^4$

Equations

$$m_D$$

$$U_{PMNS}$$

$$i⟨{\nu }_e|e^{-i{\int }_0^L\mathrm{d}xH(x)}|{\nu }_{\mu }⟩$$

$$+i{\alpha }_{21}^{\ast }$$

$$-i{V}_{cc}L\left\{\frac{\Delta m_{31}^{2}L}{4E}{e}^{-i{\delta }_{CP}}\sin {\theta }_{13}\sin {\theta }_{23}+\frac{1}{2}\frac{\Delta m_{21}^{2}L}{4E}\sin (2{\theta }_{12})\cos {\theta }_{23}\right\},$$

$$i⟨{\bar{\nu }}_e|e^{-i{\int }_0^L\mathrm{d}xH(x)}|{\bar{\nu }}_{\mu }⟩$$

$$+i{\alpha }_{21}$$

$$+i{V}_{cc}L\left\{\frac{\Delta m_{31}^{2}L}{4E}{e}^{i{\delta }_{CP}}\sin {\theta }_{13}\sin {\theta }_{23}+\frac{1}{2}\frac{\Delta m_{21}^{2}L}{4E}\sin (2{\theta }_{12})\cos {\theta }_{23}\right\},$$

$$|⟨{\nu }_e|e^{-i{\int }_0^L\mathrm{d}xH(x)}|{\nu }_{\mu }⟩{|}^2$$

$$\supset -4\sin \left(\frac{\Delta m_{31}^{2}L}{4E}\right)\sin {\theta }_{13}\sin {\theta }_{23}\mathrm{Im}\left({\alpha }_{21}{e}^{-i\Delta m_{31}^{2}L/(4E)-i{\delta }_{CP}}\right)+2\frac{\Delta m_{21}^{2}L}{4E}\sin (2{\theta }_{12})\cos {\theta }_{23}\mathrm{Im}({\alpha }_{21}),$$

$$|⟨{\bar{\nu }}_e|e^{-i{\int }_0^L\mathrm{d}xH(x)}|{\bar{\nu }}_{\mu }⟩{|}^2$$

$$\supset 4\sin \left(\frac{\Delta m_{31}^{2}L}{4E}\right)\sin {\theta }_{13}\sin {\theta }_{23}\mathrm{Im}\left({\alpha }_{21}{e}^{i\Delta m_{31}^{2}L/(4E)-i{\delta }_{CP}}\right)-2\frac{\Delta m_{21}^{2}L}{4E}\sin (2{\theta }_{12})\cos {\theta }_{23}\mathrm{Im}({\alpha }_{21}),$$

$$|{\nu }_{\alpha }⟩$$

$$P({\nu }_{\alpha }\to {\nu }_{\beta })$$

$$H(x)$$

$$+{\left[U^{†}\mathrm{diag}\left(V_{nc}(x)+V_{cc}(x),V_{nc}(x),V_{nc}(x),0,...,0\right)U\right]}_{3\times 3},$$

$$V_{nc}(x)L$$

$$V_{cc}(x)L$$

$$U$$

$$({\Omega }_{n,m}{)}_{kk}$$

$$({\Omega }_{n,m}{)}_{nm}$$

$$U$$

$${\left[U^{†}\mathrm{diag}\left(V_{nc}(x)+V_{cc}(x),V_{nc}(x),V_{nc}(x),0,...,0\right)U\right]}_{3\times 3}$$

$$=(U_{3\times 3}{)}^{†}\mathrm{diag}\left(V_{nc}(x)+V_{cc}(x),V_{nc}(x),V_{nc}(x)\right)U_{3\times 3}.$$

$$P({\nu }_{\alpha }\to {\nu }_{\beta })$$

$$H(x)$$

$$\displaystyle+U_{PMNS}^{\dagger}N_{NU}^{\dagger}\left(\begin{array}[]{ccc}V_{nc}(x)+V_{cc}(x)&0&0\\ 0&V_{nc}(x)&0\\ 0&0&V_{nc}(x)\end{array}\right)N_{NU}U_{PMNS},$$

$$N_{NU}$$

$$\frac{\Gamma ({\pi }^{\pm }\to e^{\pm }\nu )}{\Gamma ({\pi }^{\pm }\to {\mu }^{\pm }\nu )}$$

$$\frac{\Gamma ({\pi }^{\pm }\to e^{\pm }\nu )}{\Gamma ({\pi }^{\pm }\to {\mu }^{\pm }\nu )}{|}_{\mathrm{experimental}}$$

$$\frac{m_e^2}{m_{\mu }^2}\frac{(1-m_e^2/m_{{\pi }^{\pm }}^2{)}^2}{(1-m_{\mu }^2/m_{{\pi }^{\pm }}^2{)}^2}(1+{ϵ}_{\mathrm{SM}})$$

$$\frac{{\alpha }_{11}^2}{{\alpha }_{22}^2+|{\alpha }_{21}{|}^2}$$

$$\frac{1}{{\alpha }_{22}^2+|{\alpha }_{21}{|}^2}$$

$$|V_{ud}|$$

$$|V_{us}|f_{+}(0)$$

$$f_{+}(0)$$

$$\frac{1}{{\alpha }_{22}^2+|{\alpha }_{21}{|}^2}$$

$$4(1-{\alpha }_{11})(1-{\alpha }_{22})$$

$${\alpha }_{11}$$

Full text

**Heavy neutrino mixing in the T2HK, the T2HKK and

an extension of the T2HK with a detector at Oki Islands **

Yugo Abea,b,  Yusuke Asanoa,  Naoyuki Habaa  and  Toshifumi Yamadaa

a* Graduate School of Science and Engineering, Shimane University, Matsue 690-8504, Japan

b National Institute of Technology, Miyakonojo College, Miyakonojo-shi Miyazaki 885-8567, Japan *

Abstract

We study the discovery potential for the mixing of heavy isospin-singlet neutrinos in extensions of the Tokai-to-Kamioka (T2K) experiment, the Tokai-to-Hyper-Kamiokande (T2HK), the Tokai-to-Hyper-Kamiokande-to-Korea (T2HKK) with a Korea detector with $\simeq 1000$ km baseline length and 1∘ off-axis angle, and a plan of adding a small detector at Oki Islands to the T2HK. We further pursue the possibility of measuring the neutrino mass hierarchy and the standard $CP$-violating phase $\delta_{CP}$ in the presence of heavy neutrino mixing by fitting data with the standard oscillation parameters only. We show that the sensitivity to heavy neutrino mixing is highly dependent on $\delta_{CP}$ and new $CP$-violating phases in the heavy neutrino mixing matrix, and deteriorates considerably when these phases conspire to suppress interference between the standard oscillation amplitude and an amplitude arising from heavy neutrino mixing, at the first oscillation peak. Although this suppression is avoided by the use of a beam with smaller off-axis angle, the T2HKK and the T2HK+small Oki detector do not show improvement over the T2HK. As for the mass hierarchy measurement, the wrong mass hierarchy is possibly favored in the T2HK because heavy neutrino mixing can mimic matter effects. In contrast, the T2HKK and the T2HK+small Oki detector are capable of correctly measuring the mass hierarchy despite heavy neutrino mixing, as measurements with different baselines resolve degeneracy between heavy neutrino mixing and matter effects. Notably, adding a small detector at Oki to the T2HK drastically ameliorates the sensitivity, which is the central appeal of this paper. As for the $\delta_{CP}$ measurement, there can be a sizable discrepancy between the true $\delta_{CP}$ and the value obtained by fitting data with the standard oscillation parameters only, as large as $1\sigma$ resolution of the $\delta_{CP}$ measurement. Hence, if a hint of heavy neutrino mixing is discovered, it is necessary to incorporate the effects of heavy neutrino mixing to measure $\delta_{CP}$ accurately.

1 Introduction

It is a viable possibility that isospin-singlet neutrinos mix with active flavors $\nu_{e},\nu_{\mu},\nu_{\tau}$ after electroweak symmetry breaking. Notably, such a mixing is inherent in the Type-I seesaw model [1], and although the mixing angle is undetectably small in most of the parameter space, it can be sizable when there is a cancellation in the seesaw mass. To see this, suppose, for example, there are three isospin-singlet neutrinos with a Majorana mass $M_{N}$, having a Dirac mass term $m_{D}$ with active neutrinos. Since the tiny neutrino mass is derived as $-m_{D}M_{N}^{-1}m_{D}^{T}$, one can express $m_{D}$ as [2]

[TABLE]

where $U_{PMNS}$ denotes the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [3, 4], $m_{1},m_{2},m_{3}$ are the tiny neutrino mass, and $R_{3}$ is a 3-dimensional rotation matrix with complex-valued rotation angles $\alpha,\beta,\gamma$. Unless $m_{1},m_{2},m_{3}$ are degenerate, one can arbitrarily enhance $m_{D}$ and the mixing angles between isospin-singlet and active neutrinos by taking large imaginary values for $\alpha,\beta,\gamma$. In this case, a cancellation is taking place in the components of $-m_{D}^{T}M_{N}^{-1}m_{D}$. Large mixing between isospin-singlet and active neutrinos is also realized in extensions of the Type-I seesaw model, such as the inverse seesaw model [5].

In the presence of isospin-singlet neutrino mixing, transitions among active flavors require for their consistent description extra parameters beyond the standard oscillation parameters, the three mixing angles $\theta_{12},\theta_{23},\theta_{13}$, one $CP$-violating phase $\delta_{CP}$ and two mass differences $\Delta m_{12}^{2},\Delta m_{13}^{2}$. Tension in fitting neutrino oscillation data with the standard parameters would thus be evidence for isospin-singlet neutrino mixing. If isospin-singlet neutrinos are ”heavy”, by which we signify that their mass is much above the neutrino beam energy, they do not contribute to oscillations and simply alter the mixing matrix $U_{PMNS}$ into a non-unitary matrix in the following fashion:

[TABLE]

where $\alpha_{ii}$’s are real and $\alpha_{ij}$’s $(i\neq j)$ are complex numbers.

One purpose of this paper is to estimate the largest possible significance of heavy isospin-singlet neutrino mixing in the Tokai-to-Hyper-Kamiokande (T2HK) [6] and the Tokai-to-Hyper-Kamiokande-to-Korea (T2HKK) [7] (for early proposals, see Refs. [8, 9]) experiments, as well as an experiment in which a new detector at Oki Islands is added to the T2HK (for early proposals on a detector at Oki, see Refs. [10, 9]). Emphasis is on how the significance varies with $CP$-violating phases in the PMNS and the non-unitary mixing matrices and to what extent the sensitivity deteriorates with these phases. The secondary goal of this paper is to pursue the possibility that the neutrino mass hierarchy and the standard $CP$-violating phase in the PMNS matrix $\delta_{CP}$ are measured in the presence of heavy neutrino mixing by fitting neutrino data with the standard oscillation parameters only, namely, without assuming heavy neutrino mixing. This is achievable if data allow us to discern contributions from the standard oscillation and those from the non-unitary mixing matrix. Such measurement is of practical importance when a hint of heavy neutrino mixing were discovered but were not conclusive. The mass hierarchy and $\delta_{CP}$ measurement in the presence of heavy neutrino mixing has been investigated in a general context [11] and for the T2K, T2HK, NO$\nu$A and DUNE experiments [12, 13, 14, 15].

The T2HK and the T2HKK are proposed extensions of the on-going T2K experiment [17], utilizing an off-axis $\nu_{\mu}$+$\bar{\nu}_{\mu}$ beam delivered from J-PARC. In the T2HK, there are two water Cerenkov detectors, each with 187 kton fiducial volume, located solely at Kamioka, where the baseline length is $L=295$ km and the off-axis-angle is 2.5∘. In one plan of the T2HKK, one water Cerenkov detector with 187 kton fiducial volume is placed at Kamioka and the same detector is also at a site in Korea where the baseline length is $L\simeq 1000$ km and the off-axis-angle is 1∘. We also consider a plan of placing, in addition to the T2HK, a small 10 kton water Cerenkov detector at Oki Islands, where the baseline length is $L=653$ km and the off-axis-angle is 1.0∘.

Merits of the above long baseline experiments for the search for heavy neutrino mixing are as follows: (i) In disappearance experiments, the effects of heavy neutrino mixing as parametrized by Eq. (8) cancel between measurements at the near and far detectors, while the T2K and its extensions are primarily appearance experiments. (ii) $\nu_{\mu}\to\nu_{e}$ ($\bar{\nu}_{\mu}\to\bar{\nu}_{e}$) transition amplitude arising from the non-unitary mixing matrix Eq. (8) interferes with the standard oscillation amplitude, enhancing the impact of non-unitary mixing on the transition probability. The merit (ii), however, bears two subtleties:

• •

$CP$-violating phases in the PMNS and the non-unitary mixing matrices possibly lead to suppression of the interference.

• •

Non-unitary mixing $\alpha_{21}$ in Eq. (8) can mimic matter effects and the effect of the standard phase $\delta_{CP}$, which distorts the measurement of the mass hierarchy and $\delta_{CP}$.

These issues become transparent by expressing the transition amplitude in the leading order of $\sin\theta_{13}$, $\Delta m_{21}^{2}$, $\alpha_{21}$ and matter effects as

[TABLE]

where $E$ denotes the neutrino energy and $L$ the baseline length, and $V_{cc}$ is the potential due to charged current interaction inside matter.

We notice that the following interference terms,

[TABLE]

are simultaneously suppressed at the first oscillation peak, $\Delta m_{31}^{2}L/(4E)=\pi/2$, if Arg($\alpha_{21})\simeq\delta_{CP}\pm\pi/2$ holds accidentally (note the relation $\Delta m_{21}^{2}L/(4E)\simeq 0.1\cdot 2\sin\theta_{13}$ at the first oscillation peak). This suppression is mitigated if the energy distribution is not concentrated at the first oscillation peak. We thus infer that when Arg($\alpha_{21})\simeq\delta_{CP}\pm\pi/2$ holds, the T2HKK with a Korea detector placed at $1^{\circ}$ off-axis angle exhibits better sensitivity to $\alpha_{21}$ than the T2HK, and that the extension of the T2HK with an Oki detector gives qualitative improvement. Of course, these anticipations depend on statistics, namely, competition with the loss of statistics in the T2HKK, or the size of an Oki detector.

We observe in Eqs. (9), (10) that when $\delta_{CP}\simeq 0,\pi$, the real part of non-unitary mixing Re($\alpha_{21}$) plays exactly the same role as matter effects, which afflicts the measurement of the mass hierarchy because it is probed through interference with matter effects. However, when measurements with different baselines and hence different $V_{cc}L$’s are combined, one resolves degeneracy between the non-unitary mixing and matter effects, thereby resurrecting sensitivity to the mass hierarchy. We thus expect that the mass hierarchy measurement in the T2HK is subject to strong influence from heavy neutrino mixing, whereas the T2HKK maintains sensitivity to the mass hierarchy even with heavy neutrino mixing. Also, the addition of a small detector at Oki to the T2HK leads to qualitative improvement in the mass hierarchy measurement.

Eqs. (9), (10) tell us that when $\delta_{CP}\simeq$Arg$(\alpha_{21})$, non-unitary mixing $\alpha_{21}$ and the first term in the first line of Eqs. (9),(10), which is proportional to $e^{i\delta_{CP}}$, interfere maximally. In such a case, the presence of $\alpha_{21}$ affects the measurement of $\delta_{CP}$ and causes a discrepancy between the true $\delta_{CP}$ and the value measured by fitting data with the standard oscillation parameters only. This effect is mitigated if the baseline is longer, because more information about the energy dependence of the standard oscillation amplitude helps distinguish it from non-unitary mixing $\alpha_{21}$. We thus expect that the discrepancy is smaller in the T2HKK than in the T2HK.

In the main body of this paper, we demonstrate, through a simulation, the following:

(a) The sensitivity to heavy neutrino mixing is suppressed for Arg($\alpha_{21})\simeq\delta_{CP}\pm\pi/2$. Contrary to our anticipations, the T2HK shows better sensitivity than the T2HKK for any values of $CP$-violating phases, simply because of higher statistics. Also, adding a small detector at Oki to the T2HK does not make qualitative difference.

(b) In the presence of heavy neutrino mixing, the wrong mass hierarchy can be favored in the T2HK. The T2HKK maintains sensitivity to the mass hierarchy and measures it correctly despite heavy neutrino mixing. The extension of the T2HK with a small detector at Oki exhibits revived sensitivity to the mass hierarchy, even if statistics at Oki is quite subdominant compared to that at Kamioka.

(c) In the presence of heavy neutrino mixing, the true $\delta_{CP}$ and the value obtained by fitting data with the standard oscillation parameters can have a sizable discrepancy that is of the same order as $1\sigma$ resolution of the $\delta_{CP}$ measurement. The discrepancy tends to be larger in the T2HK than in the T2HKK. Adding a small Oki detector to the T2HK does not make any difference due to its limited statistics.

Our most prominent finding is related to (b), which is that adding a small 10 kton detector at Oki to the T2HK leads to revived sensitivity to the mass hierarchy.

This paper is organized as follows: In Section 2, we present the formalism for incorporating the effects of heavy isospin-singlet neutrino mixing with a non-unitary mixing matrix. In Section 3, we estimate the largest possible significance of heavy neutrino mixing in the T2HK, the T2HKK and the extension of the T2HK with a small detector at OKi, and further study the mass hierarchy and $\delta_{CP}$ measurement in the presence of heavy neutrino mixing. Here, we describe our benchmark model, procedures for simulating the long baseline experiments, and analysis with a $\chi^{2}$ fitting. Section 4 summarizes the paper.

2 Formalism for incorporating heavy isospin-singlet neutrino mixing

Suppose we have $N-3$ isospin-singlet neutrinos $\nu_{s_{f}}\,(f=4,5,...,N)$ in addition to the three active flavors $\nu_{e},\nu_{\mu},\nu_{\tau}$. Accordingly, there are $N$ mass eigenstates $\nu_{i}\,(i=1,2,...,N)$ belonging to mass eigenvalues $m_{i}^{2}$ and the neutrino mixing is expressed as

[TABLE]

where $U_{\alpha i}$ is a $N\times N$ unitary matrix. In this paper, the $N-3$ isospin-singlet neutrinos are assumed to be heavier than about 10 GeV.

In long baseline neutrino oscillation experiments, only the light mass eigenstates $\nu_{1},\nu_{2},\nu_{3}$ are produced and propagate. Therefore, the probability for a neutrino of active flavor $\alpha$ with energy $E$ to transition to active flavor $\beta$ after traveling a distance $L$ inside matter is

[TABLE]

where $U_{3\times 3}$ denotes the $\alpha=e,\mu,\tau$ and $i=1,2,3$ part of the mixing matrix $U$ in Eq. (13), and $(U^{\dagger})_{3\times 3}$ denotes the same part of the matrix $U^{\dagger}$. Note that $U_{3\times 3}\neq(U^{\dagger})_{3\times 3}$ in general. The factor $-\sqrt{E^{2}-m_{i}^{2}}$ $(i=1,2,3)$ is because a phase shift due to time evolution, $E{\rm d}t$, is equivalent to $-\sqrt{E^{2}-m_{i}^{2}}{\rm d}x$ for the mass eigenstate of mass $m_{i}$. $V_{nc}(x)$ is the potential of the neutral current interaction inside matter divided by the neutrino velocity, and $V_{cc}(x)$ is that of the charged current interaction, which are evaluated to be

[TABLE]

where $n_{n}$ and $n_{e}$ respectively denote the neutron and electron number density, $Y_{n}$ and $Y_{e}$ are respectively the number of neutrons and electrons per one nucleon, which are typically $1/2$, and $\rho(x)$ is the matter density. $[...]_{3\times 3}$ in Eq. (16) signifies that one extracts the $i,j=1,2,3$ part of the $N\times N$ matrix.

It is convenient to parametrize $U$ in terms of mixing angles, $\theta_{ij}$, and $CP$-violating phases, $\phi_{ij}$ 111 Only $(N-1)(N-2)/2$ linear combinations out of $N(N-1)/2$ phases $\phi_{ij}$’s are physical.

, $(i<j,\,i,j=1,2,...,N)$ as [19, 16] (for the case with $N=6$, see also Ref. [20])

[TABLE]

Note that the PMNS mixing matrix $U_{PMNS}$ corresponds to $(\Omega_{2,3}\Omega_{1,3}\Omega_{1,2})$ with phases $\phi_{12}$, $\phi_{23}$ set to 0 by the phase redefinition of $|\nu_{1}\rangle,\ |\nu_{2}\rangle$. We can thus express $U$ in terms of the PMNS mixing matrix as

[TABLE]

where $S_{3\times N-3}$, $W_{N-3\times 3}$ and $T_{N-3\times N-3}$ denote $3\times(N-3)$, $(N-3)\times 3$ and $(N-3)\times(N-3)$ matrices, respectively, $\alpha_{11}$, $\alpha_{22}$ and $\alpha_{33}$ are real numbers and $\alpha_{21}$, $\alpha_{31}$ and $\alpha_{32}$ are complex ones.

Since the mixing of isospin-singlet and active neutrinos is tiny, we make the approximation that we ignore terms in the second order of the mixing angles of active and isospin-singlet neutrinos. This gives $(U^{\dagger})_{3\times 3}=(U_{3\times 3})^{\dagger}$, and also

[TABLE]

We thus arrive at the following formula for the transition probabilities.

[TABLE]

where we have made a further approximation with $E^{2}\gg m_{1}^{2},m_{2}^{2},m_{3}^{2}$ and discarded terms proportional to the unit matrix in $H(x)$. The transition probabilities for antineutrinos, $P(\bar{\nu}_{\alpha}\to\bar{\nu}_{\beta})$, are obtained by changing $\delta_{CP}\to-\delta_{CP}$ in $U_{PMNS}$, $\alpha_{ij}\to\alpha_{ij}^{*}$ in $N_{NU}$, and $V_{cc}(x)\to-V_{cc}(x),\ V_{nc}(x)\to-V_{nc}(x)$ in Eq. (28).

3 Significance of heavy neutrino mixing and its effects on the mass hierarchy and $\delta_{CP}$ measurement

We estimate the largest possible significance of heavy neutrino mixing and its effects on the mass hierarchy and $\delta_{CP}$ measurement, through the following simulation: First we introduce a benchmark model where $\alpha_{11},\alpha_{22},|\alpha_{21}|$ in Eq. (8) are set at the current experimental bounds, in order to evaluate maximal impact of heavy neutrino mixing. Neutrino detection in the T2HK, the T2HKK, and the extension of the T2HK with a small detector at Oki is simulated based on the above benchmark, with Particle Data Group values employed for the standard parameters $\theta_{12},\theta_{23},\theta_{13},\Delta m^{2}_{21},|\Delta m^{2}_{31}|$. Since no conclusive data on the mass hierarchy and $\delta_{CP}$ exist, and the phase of $\alpha_{21}$ is undetermined in the benchmark, we repeat the simulation for both mass hierarchies and for various values of $\delta_{CP}$ and Arg$(\alpha_{21})$. To assess merits of the extension of the T2HK with an Oki detector, we additionally consider, for comparative study, an experiment where the detector at Oki is instead placed at Kamioka. Finally, we perform a $\chi^{2}$ fit of the simulation results under the assumption that no heavy neutrino mixing is present, namley, $N_{NU}=I_{3}$ in Eq. (24). The minimum of $\chi^{2}$ is the square of the significance of heavy neutrino mixing. The mass hierarchy and the value of $\delta_{CP}$ with which $\chi^{2}$ is minimized correspond to their measured values obtained without assuming heavy neutrino mixing.

3.1 Benchmark model

Our analysis is based on the following benchmark model: We employ the parametrization of Eq. (24). Since $\nu_{\mu}\to\nu_{e}$ process is most affected by $\alpha_{11},\,\alpha_{22}$ and $\alpha_{21}$, we further set

[TABLE]

At present, the most severe constraint on $\alpha_{11},\,\alpha_{22}$ and $\alpha_{21}$ derives from experimental tests of lepton flavor universality and a mathematical inequality among $\alpha_{11},\,\alpha_{22}$ and $\alpha_{21}$. These are enumerated below:

• •

One test of lepton flavor universality has been done in $\pi^{\pm}$ decays. When the non-unitary matrix Eq. (32) enters into the neutrino mixing Eq. (13), the ratio of the branching fractions becomes

[TABLE]

where $\epsilon_{{\rm SM}}$ indicates a Standard Model correction taking into account final state radiation. The experimental value and the Standard Model prediction are respectively given by [21, 22, 23, 27, 28]

[TABLE]

from which we obtain the following bound:

[TABLE]

• •

Another test of lepton flavor universality comes from unitarity test of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In the presence of the non-unitary matrix Eq. (32), the Fermi constant measured in the muon decay becomes $G_{\mu}=G_{F}\,\alpha_{11}\sqrt{\alpha_{22}^{2}+|\alpha_{21}|^{2}}$, whereas the Fermi constant measured in decays of hadrons involving $e\nu$ becomes $G_{\beta}=G_{F}\,\alpha_{11}$. Deviation of $G_{\beta}/G_{\mu}$ from 1 mimics non-unitarity in the first row of the CKM matrix if $|V_{ud}|$ is measured in $\beta$-decays of nuclei and $|V_{us}|$ is measured in $K\to\pi e\nu$ decays, since $G_{\mu}$ is regarded as the true Fermi constant in these measurements while it is $G_{\beta}$ that is involved in the processes, and $|V_{ub}|^{2}$ is known to be smaller than experimental error of $|V_{us}|^{2}$. We therefore have

[TABLE]

Measurements of superallowed $\beta$-decays of nuclei yield the following experimental value [24]:

[TABLE]

Combining measurements of the branching ratios of only three processes $K^{\pm}\to\pi^{0}e\nu$, $K_{L}\to\pi^{\pm}e\nu$ and $K_{S}\to\pi^{\pm}e\nu$ [24], we obtain the following experimental value for $|V_{us}|$ times the kaon form factor, $f_{+}(0)$:

[TABLE]

The average of lattice calculations of the kaon form factor $f_{+}(0)$ in Refs. [25, 26] is given by

[TABLE]

From Eqs. (37), (38), (39), we obtain the following bound:

[TABLE]

• •

The mathematical inequality among $\alpha_{11},\,\alpha_{22}$ and $\alpha_{21}$ stems from the fact that the non-unitary mixing matrix $N_{NU}$ is derived from the product of mixing matrices as in Eq. (18). The inequality reads [27, 28, 29]

[TABLE]

We combine the 3$\sigma$ experimental bounds in Eqs. (35), (40) and the mathematical inequality Eq. (41), and obtain the following 3$\sigma$ bound 222 A different bound on $\alpha_{11},\,\alpha_{22}$ and $|\alpha_{21}|$ has been derived in Ref. [16] and used to determine a benchmark in Refs. [12, 13, 14, 15], which takes into account experimental bounds on lepton flavor universality similar to Eqs.(35), (40) and the NOMAD bound Eq. (48), but does not accommodate the mathematical inequality Eq. (41). :

[TABLE]

To estimate the maximum significance of heavy neutrino mixing, we set the values of $\alpha_{11},\,\alpha_{22}$ and $|\alpha_{21}|$ at the above 3$\sigma$ bounds and therefore employ a benchmark model with the following non-unitary mixing matrix $N_{NU}$:

[TABLE]

Since no experimental bound is reported on the phase of $\alpha_{21}$, we leave it as a free parameter.

We comment in passing on constraints from past neutrino oscillation experiments. Among those experiments, the NOMAD experiment [30], a short baseline experiment searching for $\nu_{\mu}\to\nu_{e}$ and $\bar{\nu}_{\mu}\to\bar{\nu}_{e}$ transitions, gives the most severe bound, which is translated into the following:

[TABLE]

Our benchmark Eq. (47) satisfies the above bound.

Finally, we mention constraints from lepton flavor violating processes. These constraints can always be evaded by tuning heavy neutrino mass $m_{4},m_{5},m_{6},...,m_{N}$ around the $W$ boson mass, because the amplitude of a lepton flavor violating process involving charged leptons $\ell_{\alpha}$ and $\ell_{\beta}$, $A_{\alpha\beta}$, depends on $N_{NU}$ through the following combination:

[TABLE]

where $m_{i},m_{j}$ are neutrino mass eigenvalues, $M_{W}$ is the $W$ boson mass, and $F(x)$ is a function depending on the process, and in the second line, we make the approximation with $m_{1},m_{2},m_{3}\ll M_{W}$ as $m_{1},m_{2},m_{3}$ correspond to the tiny neutrino mass. It is evident that one can negate the term $(N_{NU}N_{NU}^{\dagger})_{\beta\alpha}F(0)$ by taking appropriate $O(1)$ values for $m_{j}/M_{W}$ $(j=4,5,6,...,)$. We therefore neglect any constraints from lepton flavor violating processes.

In Appendix C, we consider another benchmark that satisfies the experimental bounds of Eqs. (35), (40) at 2$\sigma$ level as well as the mathematical inequality Eq. (41).

3.2 Simulation

We simulate the neutrino propagation and detection based on the benchmark model Eq. (47). The oscillation probability is calculated numerically with the formula Eq. (28) without making any approximation. For the standard oscillation parameters, we employ the central values of $\theta_{12},\,\theta_{13},\,\theta_{23},\,\Delta m^{2}_{21},\,|\Delta m^{2}_{32}|$ found in the Particle Data Group [24]. As the $CP$-violating phase $\delta_{CP}$ and the sign of $\Delta m_{32}^{2}$ have not been measured conclusively, and the phase of $\alpha_{21}$, Arg$(\alpha_{21})$, is undetermined in the benchmark, we repeat the simulation for various values of $\delta_{CP}$ and Arg$(\alpha_{21})$ and for both cases with $\Delta m_{32}^{2}>0$ and $\Delta m_{32}^{2}<0$. The values of parameters used in the simulation are shown in Table 1.

We assume the baseline length $L$, fiducial volume, matter density along the baseline and neutrino beam off-axis angle [9, 18] for the water Cerenkov detectors used in the T2HK, the T2HKK, the extension of the T2HK with a detector at Oki (T2HK+Oki), and the extension with the same detector at Kamioka (T2HK+Kami) as Table 2. The matter density is approximated to be uniform [31, 18].

We assume that J-PARC operates with 1.3 MW beam power for 10 years in the neutrino mode and for another 10 years in the antineutrino mode, delivering $27\times 10^{21}$ proton-on-target (POT) flux of neutrino-focusing beam and $27\times 10^{21}$ POT flux of antineutrino-focusing beam, as Table 3.

The energy distributions of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ in neutrino-focusing and antineutrino-focusing beams at each site are calculated based on Ref. [32]. $\nu_{e}$ and $\bar{\nu}_{e}$ components in the beam are ignored in this analysis. In Appendix A, we present the energy distributions.

The cross sections for charged current quasi-elastic scattering between a neutrino and a proton, $\nu_{\ell}\,n\to\ell^{-}\,p$, and that between an antineutrino and a neutron, $\bar{\nu}_{\ell}\,p\to\ell^{+}\,n$, ($p$ and $n$ denote proton and neutron, respectively, and $\ell$ denotes $e$ or $\mu$) are quoted from Ref. [33]. In Appendix B, we present the cross sections.

We assume the ideal measurement of $\nu_{e}$ and $\bar{\nu}_{e}$ and do not take into account acceptance and efficiency factors and finite energy resolution. Signals coming from processes other than charged current quasi-elastic scattering, as well as backgrounds due to neutral current quasi-elastic scatterings and other processes are not considered in the analysis.

The total numbers of neutrino events at the detectors at Kamioka, Oki and in Korea in the T2HK, T2HKK, T2HK+Oki, and T2HK+Kami with the detector size of Table 2 and the beam flux of Table 3 are tabulated in Table 4.

We calculate the numbers of $\nu_{e}$ and $\bar{\nu}_{e}$ events in 0.05 GeV bins of the neutrino energy $E$ in the range 0.4 GeV$\leq E\leq$3 GeV for various values of $\delta_{CP}$ and Arg$(\alpha_{21})$ and for both cases with $\Delta m_{31}^{2}>0$ and $\Delta m_{31}^{2}<0$. We remind that in real experiments, the neutrino flux is determined by the measurement of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ at the near detector, for which the standard oscillation is negligible but the non-unitary mixing reduces the $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ flux by the following amount:

[TABLE]

In the analysis, we consider a virtual neutrino flux that is larger than the one found in Ref. [32] by $1/(|\alpha_{21}|^{2}+\alpha_{22}^{2})$, to compensate the reduction of Eq. (50). We therefore compute the number of events with a neutrino-focusing beam, $N_{e,i}$, and that with an antineutrino-focusing beam, $\tilde{N}_{e,i}$, as

[TABLE]

where $\Phi_{\nu_{\mu}}(E)$ and $\Phi_{\bar{\nu}_{\mu}}(E)$ respectively denote the neutrino and antineutrino flux per energy in a neutrino-focusing beam, and $\tilde{\Phi}_{\nu_{\mu}}(E)$ and $\tilde{\Phi}_{\bar{\nu}_{\mu}}(E)$ are those in an antineutrino-focusing beam. $P$ denotes the transition probability Eq. (28), with $E$ dependence made explicit. $N_{n}$ and $N_{p}$ are respectively the number of neutrons and protons in a water Cerenkov detector. $\sigma$ denotes the cross sections for $\nu_{\ell}n\to\ell^{-}p$ and $\bar{\nu}_{\ell}p\to\ell^{+}n$ processes.

3.3 $\chi^{2}$ analysis

We fit the numbers of events in the bins with neutrino-focusing and antineutrino-focusing beams, detected at Kamioka/Korea/Oki, $N_{e,i,{\rm site}}$ and $\tilde{N}_{e,i,{\rm site}}$ $(i=1,2,...,52)$ (site=Kamioka, Korea, Oki), under the assumption that no heavy neutrino mixing is present, that is, $N_{NU}=I_{3}$. This is performed by minimizing $\chi^{2}(\Pi)$,

[TABLE]

with respect to the standard oscillation parameters $\Pi=(\theta_{12},\theta_{13},\theta_{23},\Delta m^{2}_{21},\Delta m^{2}_{32},\delta_{CP})$. Here, $\chi^{2}_{{\rm statistical}}$ represents statistical uncertainty, $N_{e,i,{\rm site}}^{{\rm unitary}}(\Pi)$ is the number of $\nu_{e}$ and $\bar{\nu}_{e}$ events in a bin with a neutrino-focusing beam when $N_{NU}=I_{3}$, calculated as a function of the standard oscillation parameters $\Pi$, and $\tilde{N}_{e,i,{\rm site}}^{{\rm unitary}}(\Pi)$ is the corresponding number with an antineutrino-focusing beam.

The following simplification is made for $\chi^{2}_{{\rm parametrical}}(\Pi)$: We note that $\theta_{12},\theta_{13},\theta_{23},\Delta m^{2}_{21},|\Delta m^{2}_{32}|$ can be accurately measured in disappearance experiments with $\bar{\nu}_{e}\to\bar{\nu}_{e}$, $\nu_{\mu}\to\nu_{\mu}$ and $\bar{\nu}_{\mu}\to\bar{\nu}_{\mu}$, which are not influenced by heavy neutrino mixing as parametrized by $N_{NU}$, because $N_{NU}$ equally affects the neutrino flux of the same flavor at the near and far detectors. ($\cos\delta_{CP}$ and the sign of $\Delta m^{2}_{32}$ are also possibly measured in disappearance experiments, but the sensitivity is limited.) We therefore assume that the true values of $\theta_{12},\theta_{13},\theta_{23},\Delta m^{2}_{21},|\Delta m^{2}_{32}|$ are known prior to the experiments considered in this paper, and accordingly, we fix $\theta_{12},\theta_{13},\theta_{23},\Delta m^{2}_{21},|\Delta m^{2}_{32}|$ in $\chi^{2}(\Pi)$ at the values used in the simulation and set $\chi^{2}_{{\rm parametrical}}(\Pi)=0$.

Finally, we have $\chi^{2}_{{\rm systematic}}=0$ in the current analysis, because no acceptance, efficiency, energy resolution and background events are considered.

To summarize, $\chi^{2}(\Pi)$ is approximated as

[TABLE]

We numerically evaluate the minimum of $\chi^{2}(\Pi)$ Eq. (53), which corresponds to the square of the significance for heavy neutrino mixing. In Figure 1, we show contour plots of the minimum of $\chi^{2}(\Pi)$,

[TABLE]

on the plane spanned by the $CP$-violating phase in the PMNS matrix $\delta_{CP}$ and the phase in the non-unitary mixing matrix Arg($\alpha_{21}$), when the true mass hierarchy is normal. The four subplots correspond to the T2HK, the T2HKK, the extension of the T2HK with a detector at OKi, and that with the same detector at Kamioka.

We study the impact of heavy neutrino mixing on the mass hierarchy measurement performed by fitting data without incorporating heavy neutrino mixing. For this purpose, we evaluate the difference between the minimum (with respect to $\delta_{CP}$ only) of $\chi^{2}(\Pi)$ for the wrong hierarchy and that for the true mass hierarchy, which corresponds to the square of the significance of rejecting the wrong mass hierarchy. We show in Figure 2 contour plots of the difference,

[TABLE]

on the same plane when the true mass hierarchy is normal. The black-filled regions in the left subplots are where the wrong mass hierarchy is favored over the true mass hierarchy, i.e., $\min\chi^{2}({\rm wrong\ H})-\min\chi^{2}({\rm true\ H})<0$.

We study the impact of heavy neutrino mixing on the $\delta_{CP}$ measurement performed by fitting data without incorporating heavy neutrino mixing. To assess the impact, we compute the difference between the true value of $\delta_{CP}$, and the value that minimizes $\chi^{2}(\Pi)$ of Eq. (53) with the true mass hierarchy. We show in Figure 3 contour plots of the difference,

[TABLE]

on the same plane when the true mass hierarchy is normal.

Figures 4, 5, 6 are the corresponding figures when the true mass hierarchy is inverted.

The following observations are made:

(A) In Figures 1,4, we see that the significance of heavy neutrino mixing can be above $3\sigma$ for Arg($\alpha_{21})\sim\delta_{CP}$, but decreases considerably for Arg($\alpha_{21})\simeq\delta_{CP}\pm\pi/2$ in all the experiments, confirming that the sensitivity deteriorates due to suppression of interference between the standard oscillation amplitude and non-unitary mixing $\alpha_{21}$ at the first oscillation peak, as can be read from Eqs. (11), (12). The regions of depleted sensitivity is slightly off the lines Arg($\alpha_{21})=\delta_{CP}\pm\pi/2$, because of the influence of the second, subleading interference term in Eqs. (11), (12). We also find that the significance does not ameliorate in the T2HKK and in the extension of the T2HK with a 10 kton detector at Oki. Rather, the T2HKK shows inferior performance due to the loss of statistics, which is seen in Table 4.

(B) In Figures 2,5, we observe that in the T2HK, the wrong mass hierarchy is favored when $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(0,\,0),\,(\pi,\,\pi)$ holds in the normal hierarchy case, and when $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(\pi,\,0)$ holds in the inverted hierarchy case. To estimate the impact of non-unitary mixing on the mass hierarchy measurement in the T2HK, we show in Figure 7 the values of $\{\min\chi^{2}({\rm wrong\ H})-\min\chi^{2}({\rm true\ H})\}$ when non-unitarity is absent, i.e. $N_{NU}=I$, calculated with the same procedure as Figures 2, 5.

Comparison of Figures 2,5 with Figure 7 indicates that for $\delta_{CP}\sim 0,\pi$, the T2HK has limited sensitivity to the mass hierarchy in the unitary case, and the presence of non-unitary mixing may aggravate it to the extent that the wrong mass hierarchy is favored. Such impact of non-unitary mixing is because $\alpha_{21}$ cancels matter effects in the neutrino transition probability Eq. (9) when $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(0,\,0),\,(\pi,\,\pi)$ holds in the normal hierarchy case, and when $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(\pi,\,0)$ holds in the inverted hierarchy case. Although $\alpha_{21}$ and matter effects add positively in the antineutrino transition probability Eq. (10) in the same cases, it does not change the situation because neutrino flux is larger than antineutrino flux in our setup (see Table 3 and Figure 9).

The T2HKK maintains strong sensitivity to the mass hierarchy in the presence of heavy neutrino mixing. Also, adding a 10 kton detector at Oki to the T2HK drastically ameliorates sensitivity to the mass hierarchy, in spite of the small size of the Oki detector. The advantage of the T2HKK and T2HK+Oki is because measurements with different matter effects resolve degeneracy between non-unitary mixing $\alpha_{21}$ and matter effects.

We stress that even a small 10 kton detector at Oki negates the impact of heavy neutrino mixing on the mass hierarchy measurement and allows correct measurement.

(C) In Figures 3,6, we find that the value of $\delta_{CP}$ that minimizes $\chi^{2}$ under the assumption of standard unitary oscillation can deviate from the true $\delta_{CP}$ by more than 0.15 radian in the T2HK, and by more than 0.1 radian in the T2HKK. The result should be compared to the resolution of $\delta_{CP}$. To this end, we present in Figure 8 the difference between $\delta_{CP}$ that gives $\chi^{2}=1$ and the true $\delta_{CP}$ when non-unitarity is absent ($N_{NU}=I$),

[TABLE]

which quantifies $1\sigma$ resolution of the $\delta_{CP}$ measurement in the standard oscillation case. Eq. (57) is calculated through a simulation performed with the same procedure.

The $1\sigma$ resolution for $\delta_{CP}$ is better than 0.1 radian in the T2HK, and better than 0.08 radian in the T2HKK, which is of the same order as the deviation of $\delta_{CP}$ due to non-unitarity. This indicates that non-unitarity can have a non-negligible effect on the measurement of $\delta_{CP}$, and hence, if a hint of non-unitarity has been discovered, it must be incorporated in the data fitting in order to measure $\delta_{CP}$ correctly.

The deviation due to non-unitarity is enhanced when $({\rm true}\ \delta_{CP})\sim{\rm Arg}(\alpha_{21})$ holds. This is because interference between the part of the standard oscillation amplitude proportional to $e^{i\delta_{CP}}$ and the amplitude involving non-unitary mixing $\alpha_{21}$ is maximized, and so is the influence of $\alpha_{21}$ on the value of $\delta_{CP}$ that minimizes $\chi^{2}$.

The deviation is smaller in the T2HKK than in the T2HK. This is due to the fact that the energy dependence of the part of the standard oscillation amplitude proportional to $e^{i\delta_{CP}}$, which helps distinguish the standard and non-unitary amplitudes, is more accurately measured with a longer baseline. Hence, the $\delta_{CP}$ measurement is less affected by non-unitarity in the T2HKK.

4 Summary

We have studied the discovery potential for the mixing of heavy isospin-singlet neutrinos in the Tokai-to-Hyper-Kamiokande (T2HK), the Tokai-to-Hyper-Kamiokande-to-Korea (T2HKK), and a plan of adding a small detector at Oki Islands to the T2HK, and further examined the feasibility of measuring the mass hierarchy and the standard $CP$-violating phase $\delta_{CP}$ in the presence of heavy neutrino mixing by fitting data without assuming heavy neutrino mixing. The mixing of heavy neutrinos is parametrized with a non-unitary mixing matrix for active flavors. A benchmark that maximizes the non-unitary mixing and is consistent with the current experimental and theoretical bounds is employed to estimate the largest possible significance of heavy neutrino mixing and its impact on the mass hierarchy and $\delta_{CP}$ measurement.

Through a simulation, we have revealed that the significance of heavy neutrino mixing can be above $3\sigma$ in all the experiments when the standard phase $\delta_{CP}$ and a new $CP$-violating phase originating from heavy neutrino mixing, Arg($\alpha_{21}$), satisfy Arg($\alpha_{21})\simeq\delta_{CP}$. The significance decreases considerably for Arg($\alpha_{21})\simeq\delta_{CP}\pm\pi/2$, due to suppression of interference between the standard oscillation amplitude and a non-unitary mixing term.

In the T2HK, the mass hierarchy measurement is so affected by heavy neutrino mixing that the wrong mass hierarchy is favored for $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(0,\,0),\,(\pi,\,\pi)$ in the normal hierarchy case, and for $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))\sim(\pi,\,0)$ in the inverted hierarchy case. This is because a non-unitary mixing term possibly cancels matter effects. In contrast, the T2HKK and the extension of the T2HK with a 10 kton detector at Oki show strong sensitivity to the mass hierarchy even with heavy neutrino mixing, because measurements with different baseline lengths at Kamioka and in Korea or at Oki resolve degeneracy between a non-unitary mixing term and matter effects. We thus conclude that (i) the mass hierarchy measurement in the T2HKK is highly stable against any effects of heavy neutrino mixing, and that (ii) although the mass hierarchy measurement in the T2HK is easily affected, the addition of a small detector at Oki negates effects of heavy neutrino mixing.

In the presence of non-unitarity, the value of $\delta_{CP}$ measured without assuming non-unitarity deviates from the true $\delta_{CP}$. The deviation is enhanced when $\delta_{CP}$ and ${\rm Arg}(\alpha_{21})$ are similar, because in such a case, the part of the standard oscillation amplitude proportional to $e^{i\delta_{CP}}$ interferes maximally with a non-unitary mixing term at oscillation peaks. The deviation can be of the same order as $1\sigma$ resolution of the $\delta_{CP}$ measurement, which suggests that if a hint of non-unitary mixing has been discovered, one must employ an appropriate ansatz including the non-unitary mixing terms to fit neutrino oscillation data and measure $\delta_{CP}$ correctly.

Acknowledgement

The authors would like to thank Kaoru Hagiwara (KEK), Hiroyuki Ishida (NCTS) and Yuya Yamaguchi (Hokkaido University). This work is partially supported by Scientific Grants by the Ministry of Education, Culture, Sports, Science and Technology of Japan (Nos. 24540272, 26247038, 15H01037, 16H00871, and 16H02189).

Appendix A

In Figure 9, we show the flux of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ in neutrino-focusing and antineutrino-focusing beams from J-PARC, detected at a water Cerenkov detector at Kamioka, Oki and in Korea if the neutrino oscillation were absent. The baseline length and beam off-axis angle assumed are given in Table 2. In each plot, the blue lines correspond to a neutrino-focusing beam and the red lines correspond to an antineutrino-focusing beam. The solid lines denote $\nu_{\mu}$ flux and the dashed lines denote $\bar{\nu}_{\mu}$ flux.

Appendix B

In Figure 10, we show the cross sections for charged current quasi-elastic scatterings $\nu_{\ell}n\to\ell^{-}p$ and $\bar{\nu}_{\ell}p\to\ell^{+}n$ ($\ell=e,\,\mu$). The solid line corresponds to the cross section for $\nu_{\ell}$, and the dashed line to that for $\bar{\nu}_{\ell}$.

Appendix C

We consider an alternative benchmark with a smaller non-unitary mixing that satisfies the experimental bounds Eqs. (35), (40) at 2$\sigma$ level as well as the mathematical inequality Eq. (41). The non-unitary mixing matrix in the new benchmark is given by

[TABLE]

Based on the benchmark with Eq. (62), we repeat the same simulation and calculation as subsection 3.2 and 3.3, and obtain the results in Figures 11, 12, 13, 14, 15, 16.

The plots for min$\chi^{2}$, Figures 11, 14, show that the significance of non-unitary mixing is reduced below $3\sigma$ in most parameter regions in the T2HK and in the entire region in the T2HKK. The dependence on $(\delta_{CP},\,{\rm Arg}(\alpha_{21}))$ is similar to the plots for the original benchmark, Figures 1, 4.

The plot for $\{\min\chi^{2}({\rm wrong\ H})-\min\chi^{2}({\rm true\ H})\}$ in the normal hierarchy case, Figure 12, tells us that in the T2HK, the wrong mass hierarchy can still be favored if the true hierarchy is normal. In contrast, Figure 15 gives that this is no longer the case if the true hierarchy is inverted.

The plots for $\Delta_{CP}$, Figures 13, 16, show that the value of $\delta_{CP}$ that minimizes $\chi^{2}(\Pi)$ of Eq. (53) still possibly deviates from the true value by more than 0.1 radian in the T2HK, and more than 0.05 radian in the T2HKK, even with the alternative benchmark with the smaller non-unitary mixing. This augments our argument that the non-unitary mixing terms must be incorporated in the $\delta_{CP}$ measurement when a hint of non-unitarity has been discovered.