Confronting Tri-direct CP-symmetry models to neutrino oscillation experiments

Tri-direct CP symmetry is an economical neutrino model building paradigm, and it allows to describe neutrino masses, mixing angles and CP violation phases in terms of four free parameters. Viability of a class of tri-direct CP models is examined with a comprehensive simulation of current and future neutrino oscillation experiments. Two benchmark models as well as the full parameter space are carefully scanned, and the parameter degeneracy problem is observed from the constraint of one group of neutrino oscillation experiments. Complementary roles from the accelerator neutrino experiments (e.g., T2HK and DUNE) and reactor neutrino experiments (e.g., JUNO) are crucial to break the degeneracy and nail down fundamental neutrino mixing parameters of the underlying theory.


I. INTRODUCTION
Neutrinos in the standard model (SM) of particle physics are strictly massless. Neutrino oscillation requires mass-squared differences and non-zero neutrino masses, which is striking new physics beyond the SM and calls for new degrees of freedom. In the framework of the three-generation neutrino oscillation paradigm, we have two mass-squared differences (∆m 2 21 , ∆m 2 31 ), three mixing angles (θ 12 , θ 13 and θ 23 ) and Dirac CP phase δ CP [1]. The precision of measuring θ 13 is dominated by reactor neutrino experiments [2][3][4], θ 12 and ∆m 2 21 are dominated by solar and reactor neutrino experiment KamLAND [5][6][7][8], as well as θ 23 and |∆m 2 31 | are dominated by atmospheric neutrino experiments [4,9,10]. A global analysis of different experiments provides the precise values of mixing parameters at the percentage level [1]. However, the mass ordering ∆m 2 31 > 0 or ∆m 2 31 < 0 and the value of δ CP remains unclear, although hints exist from experiments which are currently running.
Many models have been proposed to accommodate massive neutrinos without violating overwhelming constraints from previous experimental results. The origin of neutrino masses, flavor mixing and CP violation is a longstanding open question in particle physics. It turns out that a broken flavour symmetry based on a discrete group is particularly suitable to explain the structure of the leptonic mixing matrix, see Refs. [11][12][13][14][15] for review. If the discrete flavor symmetry is extended to also involve CP as a symmetry, the CP violation phases in the quark sector (observed) and lepton sector can be predicted [16][17][18]. Recently a new discrete flavor symmetry model building approach called tri-direct CP was proposed [19,20], and it is dictated by residual symmetries such that it is quite predictive. The light neutrino mass matrix only depends on four real free parameters to describe the entire neutrino sector (three neutrino masses as well as the lepton mixing matrix). Moreover the CP violations in neutrino oscillations and leptogenesis generally arise from the same phase in the tri-direct CP model, consequently they are closely related to each other.
Precision measurements of neutrino oscillation parameters will guide us to the new physics domain. While in the quark sector, the precision is at the sub-percentage level [21], in the neutrino sector, the parameter uncertainties remain at the percentage level [1]. New physics might be hidden in the uncertainties of measured neutrino mixing parameters, as neutrino oscillations bridge neutrino mixings and other factors which affect the propagation of coherent states. The new models to accommodate massive neutrinos intend to bring new fundamental symmetries, new particles and their new interactions beyond the standard model. It is promising to conduct precision measurements in accelerator neutrino experiments to search for new physics, including nonstandard interactions, and neutrino decays (e.g. Ref. [22][23][24][25][26][27][28][29][30][31][32]). It is a question of whether we are able to test different flavor and CP-violation models directly in running accelerator neutrino experiments like T2K [4,33]and NOνA [9,34,35], and future neutrino oscillation experiments like DUNE [36] and T2HK [37]. In this paper, we shall determine the potential of current and upcoming neutrino facilities to test the tri-direct CP approach, and the sensitivity regions of oscillation parameters will be presented.
The paper is organized as follows: we firstly review the tri-direct CP symmetry in Sec. II. In Sec. III, we investigate the precision measurements of oscillation parameters, either represented by standard three neutrino mixing parameters or denoted by the benchmark model parameters, in running experiments such as T2K and NOνA. We expect better sensitivities in future neutrino experiments, such as T2HK, DUNE and Jiangmen Underground Neutrino Observatory (JUNO). In Sec. IV, we show our simulation results. We study the precision of model parameters for different experimental configurations, and then discuss a degeneracy problem, which can be resolved by including JUNO data. We further study how to constrain oscillation parameters with the restriction of the tri-direct CP model, before we discuss two benchmark models. Finally, we summarize our results in Sec. V.

II. REVIEW OF TRI-DIRECT CP SYMMETRY MODELS
Let us firstly recapitulate a benchmark tri-direct CP model proposed in [19]. This model is based on the S 4 flavor symmetry and CP symmetry. The flavor group S 4 and CP are broken to the subgroups Z T 3 , Z T ST 2 2 × X atm and Z U 2 × X sol in the charged lepton, atmospheric neutrino and solar neutrino sectors, respectively, where S, T , U are the generators of S 4 and X atm = SU and X sol = U denote the residual CP symmetry. In the generic tri-direct CP paradigm, the structure of the neutrino and charged lepton mass matrices essentially arise from the vacuum alignment of flavon fields which are fixed by the residual symmetry. In the working basis of [19], the residual flavor symmetry Z T 3 enforces that the charged lepton mass matrix is diagonal [19]. The atmospheric and solar flavon vacuum alignments are determined to be φ atm ∝ 1, ω 2 , ω T and φ sol ∝ (1, x, x) T , where ω = e 2πi/3 is a cube root of unity and the parameter x is real because of the imposed CP symmetry. As a result, the Dirac neutrino mass matrix reads as The right-handed neutrino Majorana mass matrix is diagonal The light effective left-handed Majorana neutrino mass matrix is given by the seesaw formula: where m a = |y 2 a /M atm |, m s = |y 2 s /M sol |, and the only physically important phase η depends on the relative phase between y 2 a /M atm and y 2 s /M sol . It is noteworthy that only four parameters m a , m s , η and x are involved to describe both neutrino masses and lepton mixing parameters. As a consequence, this model is quite predictive. The low-energy phenomenology of this model has been studied both numerically and analytically in [19]. Although x and the relative phase η are free parameters in the general setup of tri-direct CP, they can be fixed to some particular values through the vacuum alignment technique in a discrete flavor symmetry model. It is found that a quite good fit to the experimental data can be obtained for certain choices of x and η, two benchmark examples are x = −7/2, η = π and x = −4, η = 5π/4. In these benchmark models, the corresponding vacuum alignments take a simple form such that they can be easily realized in concrete models [19,20]. Moreover, the neutrino mass matrix as well as neutrino masses and mixing parameters only depend on two free parameters m a and m s in the benchmark models and the experimental data can be accommodated very well.
The neutrino mass spectrum is predicted to be normal ordering in this model, and the lightest neutrino is massless m 1 = 0. The other two non-vanishing neutrino masses m 2 and m 3 are expressed in terms of the input parameters as follows. where As regards the predictions for lepton flavor mixing, the first column of the mixing matrix is determined to be proportional to ( The angles θ and ψ are specified by sin ψ = (y * z + wz * ) |y * z + wz * | , cos ψ = (y * z + wz * ) |y * z + wz * | , As a consequence, we find the exact expressions for the mixing angles are We see that the solar and reactor mixing angles satisfy the following sum rule cos 2 θ 12 cos 2 θ 13 = 3x 2 5x 2 + 2x + 2 .
If the atmospheric mixing angle is maximal, this sum rule implies that the Dirac CP phase would be maximal (i.e. δ CP = ±π/2) as well. Furthermore, the result for the Jarlskog invariant is from which we can extract the value of sin δ CP , sin δ CP = ± csc 2θ 23 1 + (x 2 + x + 1) 2 cot 2 θ 13 cos 2 2θ 23 with "+" for x cos ψ > 0 and "−" for x cos ψ < 0. The above results for cos δ CP and sin δ CP allow us to fix the value of δ CP . Comprehensive numerical analyses show that the allowed region of the parameters x, η, r = m a /m s and m a are −5.475 ≤ x ≤ −3.370, 0.455π ≤ η ≤ 1.545π, 0.204 ≤ r ≤ 0.606 and 3.343 meV ≤ m a ≤ 4.597 meV respectively in order to accommodate the experimental data on neutrino masses and lepton mixing angles [38]. It is remarkable that both solar mixing angle and Dirac CP phase are predicted to lie in a narrow range 0.329 ≤ sin 2 θ 12 ≤ 0.346 and 1.371π ≤ δ CP ≤ 1.629π in this model. For the benchmark values of the vacuum parameters x = −7/2 and η = π, the effective light neutrino mass matrix in Eq. (3) only depends on two free parameters m a and m s . Using the general results presented above, we find that the lepton mixing matrix is of the form where sin 2θ = 10|14r − 1| 13 √ 4 + 32r + 289r 2 , cos 2θ = 3 (8 + 57r) 13 with r = m s /m a . The lepton mixing angles read sin 2 θ 13 = 26 75 sin 2 θ, sin 2 θ 12 = 26 cos 2 θ 62 + 13 cos 2θ , and the Jarlskog invariant is which implies that the Dirac CP phase is exactly maximal, i.e.
Notice that both θ 23 and δ CP are predicted to be maximal and they are favored by the latest data from T2K [39] and NOνA [9,40], the reason is the neutrino mass matrix of Eq. (3) fulfills the µ − τ reflection symmetry in this case. The lightest neutrino masses as functions of m a and r are In order to accommodate the experimental values of the mixing angles and neutrino mass splittings ∆m 2 21 and ∆m 2 31 [38], we find that m a and r are constrained to lie in rather narrows regions 3.560 meV ≤ m a ≤ 3.859 meV and 0.5282 ≤ r ≤ 0.5904. Accordingly the allowed regions of the reactor and solar mixing angles are strongly constrained 0.02206 ≤ sin 2 θ 13 ≤ 0.02349 and 0.3310 ≤ sin 2 θ 12 ≤ 0.3319.
Then we proceed to discuss the second representative values of the vacuum parameters x = −4 and η = 5π/4, the lepton mixing matrix reads as where a Majorana phase matrix is omitted, and the parameters θ and ψ are functions of the mass ratio r, The expressions of the mixing angles are sin 2 θ 13 = 13 37 sin 2 θ, sin 2 θ 12 = 26 cos 2 θ 61 + 13 cos 2θ , The Jarlskog CP invariant takes the form The sum rules for the Dirac CP phase in terms of lepton mixing angles are given by It noteworthy that all the lepton mixing angles as well as δ CP only depend on the parameter r through θ and ψ in this case. Moreover, the results for the light neutrino masses are In order to describe the experimentally measured values of both lepton mixing angles and neutrino mass squared differences, we find the allowed ranges of the input parameters are 3.568 meV ≤ m a ≤ 3.871 meV and 0.3983 ≤ r ≤ 0.4473. As a consequence, the solar and reactor mixing angles are constrained to lie in the narrow intervals 0.02254 ≤ sin 2 θ 13 ≤ 0.02280 and 0.3362 ≤ sin 2 θ 12 ≤ 0.3364, and the atmospheric mixing angle is predicted to be in the second octant 0.5559 ≤ sin 2 θ 23 ≤ 0.5636. The predicted values of δ CP are distributed around 3π/2, namely 1.582π ≤ δ CP ≤ 1.594π. Since the model is very predictive and the mixing angles as well as Dirac CP phase are constrained to lie in rather narrow regions, in particular we have 0.329 ≤ sin 2 θ 12 ≤ 0.346 for the most general case, we expect the benchmark tri-direct model could be excluded in future neutrino experiments. If θ 23 and δ CP are measured precisely enough, the two values x = −7/2, η = π and x = −4, η = 5π/4 may be distinguished from each other. It will be nice to probe these features in detailed simulations of current and future neutrino oscillation experiments.

III. IMPLEMENTATION OF NEUTRINO EXPERIMENTS IN SIMULATION
In this section, we will introduce the current and future experiments -T2K, NOνA, T2HK, DUNE and JUNO. All sensitivities in experiments are simulated in a state-of-the-art tool GLoBES [41,42] where the experimental details can be very nicely implemented by an Abstract Experimental Design Language (AEDL) file. As soon as the publicly available signal and background spectra are reproduced, we can safely claim the expected sensitivities in the precision measurements. In the simulation, input values of neutrino mixing parameters are taken as the best fit values of the latest NuFit4.0 [1]: sin 2 θ 12 = 0.310, sin 2 θ 13 = 0.0224, sin 2 θ 23 = 0.580, δ CP = 215 • , ∆m 2 21 = 7.39 × 10 −5 eV 2 , ∆m 2 31 = 2.525 × 10 −3 eV 2 . In the current study, we will choose a normal mass hierarchy as a demonstration. In the meantime, the Preliminary Reference Earth Model (PREM) density profile is considered in the numerical calculations [43]. We are using two methods to present our results: • Standard three neutrino oscillations expressed by θ 12 , θ 13 , θ 23 , δ CP , ∆m 2 21 and ∆m 2 31 are taken as the truth in Nature, we expect that precision measurements of mixing parameters are correlated, and uncertainties of current global fit results are taken into account. For given oscillation parameters, we define a set of parameters: and predict the expected event rate in the bin i µ i ( − → O ). We suppose a given experiment reconstructs neutrino spectra in N bins sequentially. The event rate in the bin i is recorded as n i . We can build a χ 2 ( − → O ) to quantify the sensitivity: The final results come from a minimization of the summation of the χ 2 ( − → O ) in every oscillation channel of all experiments over a set of parameters, or the so-called marginalization.
• Once we fit the model parameters, the number of degrees of freedom is reduced from six to four, as shown in the previous section. We consider the following parameters from the tri-direct CP symmetry models: x, η, m a and r. In this case, we have to change the oscillation parameters predicted by the specific model: Other steps in the likelihood analysis will follow the same strategy as the above method, but replace the equation Eq. (26) by with the PMNS parameters as functions of model parameters We can expect better measurements of input parameters after a combination of experimental results and symmetry-induced constraints from the theory.

A. T2K
T2K stands for Tokai to Kamioka a long-baseline experiment in Japan. In Tokai, muon neutrinos or antineutrinos are produced by bombarding a 30 GeV proton beam onto a graphite target station in the J-PARC accelerator center. The neutrino beams are detected firstly at the near detectors which are 280 meters away from the target station. The far detector which reconstructs oscillated neutrino/antineutrino signals is Super-Kamiokande which has a fiducial mass of 22.5 kilotons and is 295 kilometers away with an off-axis angle of 2.5 • from the beam direction. With the carefully chosen off-axis angle, the neutrino beam energy is peaked at about 0.6 GeV and matches the first maximum in the neutrino oscillation channels: P(ν µ → ν e ) and P(ν µ →ν e ).
In 2011, the T2K collaboration published their first result on P(ν µ → ν e ) with 1.43 × 10 20 Protons On Target (POT). It is the first hint of non-zero θ 13 at 2.5σ confidence level (C.L.) [33]. In 2012, they presented an analysis of neutrino oscillation for P(ν µ → ν µ ) based on the same POT data, where we have best-fit values of ∆m 2 32 = 2.63 × 10 −3 eV 2 and sin 2 2θ 23 = 0.98 in the three-neutrino mixing framework [44]. In 2016, the first antineutrino result was published based on 4.01 × 10 20 POT, where we have best-fit values of ∆m 2 32 = 2.51 × 10 −3 eV 2 and sin 2 θ 23 = 0.45. The latest results of searching for CP violation in neutrino and antineutrino oscillations by T2K are based on 2.2 × 10 21 POT [4]. In our simulation, we equally split 7.8 × 10 21 POT into two modes for T2K as the final total POT number.

B. NOνA
NOνA is a long-baseline neutrino oscillation experiment in the United States. Muon neutrinos or antineutrinos are produced by the NuMI beam at Fermilab. The experiment also adopts an off-axis angle of 14.6 mrad to reach the first neutrino oscillation maximum at a peak energy of 2 GeV, since the far detector using 14 kt active scintillator is 810 km away from the target station. The far detector is on the surface. An identical detector with a mass of 290 ton scintillator is 100 meter deep at a distance of 1 km in order to monitor the neutrino flux and cancel the systematic uncertainties.
In 2016, the NOνA collaboration published their first result in the ν µ → ν e channel [34] and in the ν µ → ν µ channel [35] based on 2.74×10 20 POT. In 2017, they updated results on the electron neutrino appearance channel based on 6.05 × 10 20 POT [45]. The degeneracy of θ 23 shows up at 2.6σ C.L. The latest results in a combination of neutrino and antineutrino runs are given in Ref. [9]. In our simulation, we assume total 36 × 10 20 POT for ν andν modes until 2024 for NOνA.

C. T2HK
An evolution of Water Cherenkov detectors from Kamiokande to Hyper-Kamiokande makes it possible to conduct an upgrade of T2K to T2HK [37]. The HyperK detector will have 560 kt fiducial mass to reconstruct neutrino oscillation spectra. T2HK shares the same baseline of 295 km as T2K while the offaxis beam remains in the same direction with an upgraded proton beam at 1.3 MW. We assume T2HK is running in the neutrino mode in 2.5 years and in the antineutrino mode in 7.5 years. The second far detector in Korea is actively under consideration. In our simulation, we will keep the conservative option without the second far detector.

D. DUNE
DUNE is the next-generation accelerator neutrino oscillation experiment with a baseline of 1300 km from FNAL to the underground laboratory in South Dakota. The experiment will search for CP violation in the leptonic sector and conduct precision measurements using appearance and disappearance channels by ν µ andν µ beams. DUNE is going to reconstruct oscillated neutrino spectra with a detector complex of four 10-kt Liquid Argon Time Projection Chamber (LArTPC). We adopt an AEDL file provided by Ref. [36]. We assume the experiment will be running in the neutrino/antineutrino mode in 3.5 years, and adopt the 3-horn optimised beam design, which consists of the 62.5 GeV proton beam with a power of 1.83 × 10 21 POT per year [46,47].

E. JUNO
JUNO is a multi-purpose underground neutrino experiment, which will build a 20 kt liquid scintillator detector in South China and is planned to be online in 2021 [48]. The primary goal of JUNO [49] is to determine the neutrino mass ordering and precision measurement of oscillation parameters using the reactor electron neutrino disappearance channel thanks to the unprecedented energy resolution of 3%/ √ E. Regarding the precision measurement of |∆m 2 31 |, ∆m 2 21 and sin 2 θ 12 , JUNO can reach the levels of 0.44%, 0.59% and 0.67% respectively, after six years of data taking. Moreover, the determination of the neutrino mass ordering at reactors is free from the contamination of matter effects [50] and possible new physics [51]. JUNO will be rather robust when combined with accelerator neutrino experiments. The sub-percent level precision for three of all six standard oscillation parameters will certainly be powerful for selecting the flavor-symmetry models. In our simulation, we use the standard precision levels as our input priors to combine with accelerator neutrino experiments using the state-of-the-art GLoBES tool.

IV. MODEL TESTING WITH NEUTRINO OSCILLATION EXPERIMENTS
In this section, we show our simulation results with the experiments introduced in Sec. III. The configurations we considered are the synergy of T2K and NOνA, DUNE, T2HK, the combination of all LBLs, and the interplay of LBLs and the reactor experiment JUNO. In Sec. IV A, we will firstly investigate the precision of four model parameters -x, η, r and m a . One will see in those results, there is a degeneracy problem. In Sec. IV B, we will explain how this degeneracy problem appears, and propose a way to resolve this problem. In the following section, we will study how the uncertainties are changed for oscillation parameters by tri-direct CP models. The above subsections are based on the general tri-direct CP model. In Tab. I, we show two benchmark models. Finally, we will predict how these two benchmark models can be tested in future experiments.

A. Precision measurement of model parameters
In Fig. 1, we show the ∆χ 2 values for each model parameter. We use the true values for the model parameters; (x, η, r, m a ) = (−3.65, 1.1π, 0.5, 3.7meV), which is the best fit for the NuFit4.0 result. We also include the prior according to the NuFit4.0 result. We consider the configurations: DUNE (dashed-blue), T2HK (dashedgreen), the combination of NOνA and T2K (dashed-brown), the synergy of all LBLs (dashed-grey), and the optimised configuration by combing all LBLs and the JUNO experiment (solid-red). Except for the m a result, we see a great improvement in the DUNE result compared to the T2K-and-NOνA combination. T2HK further improves the measurements, and its performance is similar to the combination of all LBLs. This demonstrates the fact that T2HK dominates the contributions. The feature that the performance of T2HK is better than that of DUNE reflects the well-known result that T2HK works better than DUNE with fixed mass ordering, which is naturally imposed by the tri-direct CP model. In more detail, for , respectively. The 3σ uncertainty for combining all LBLs is almost the same as that for T2HK. The relative symmetry is seen in the result for η. The size of the 3σ uncertainty for the T2K-and-NOνA combination, DUNE and T2HK are about 0.3π, 0.2π, 0.15π, respectively. The correlation between η and r worsens the sensitivity for η smaller than the assumed true value. Thanks to the high precision of T2HK and combing all LBLs, the degeneracy problem can be resolved when η is very close to the true value. Therefore, we see a twist around η = π for these two configurations. Details about this degeneracy will be introduced in Sec. IV B.
For the above three parameters x, η and r, it is hard to see the improvement by including the data of JUNO to those of all LBLs. Data from JUNO is important for the m a measurement. We see the overlapping of all curves for all LBL configurations (dashed-blue, dashed-green, dashed-brown, black curves) in the m a result. The uncertainty is mainly contributed from ∆m 2 21 , which is not measured well by LBLs. As a result, we see a great improvement by including data from JUNO, which well measures ∆m 2 21 . We also show the 3σ (∆χ 2 = 11.83) contour between any two model parameters in Fig. 2. We see some correlations among x, η and r for all configurations on x − η, x − r and η − r planes. This correlation is consistent with what we see in Eq. (3), in which m a is less dependent on the other three parameters. We discover a degeneracy problem related to this correlation for all possible LBL configurations-the combination of NOνA and T2K, DUNE, T2HK, and all-LBL synergy. This degeneracy is mainly caused by the poor measurement of θ 12 . More details about this degeneracy can be seen in Sec. IV B. These correlations are not removed even if we include JUNO, but combing LBLs and JUNO data can resolve the degeneracy problem.

B. Breaking degeneracies
The degeneracy in Fig. 2 can be understood by the equal-oscillation-parameter-value contours on different planes as shown in Fig. 3. In Fig. 3, we show these contours on the x − η, x − r, and η − r planes. In the upper panels, we set the model paramters at the true values (x, η, r, m a ) = (−3.65, 1.1π, 0.5, 3.7 meV), which predicts the value for oscillation parameters θ 12 ∼ 35. (short-dashed-black). The contours are shown with these conditions. Therefore, the intersection of all contours is at the assumed true values. In the lower panels, we focus on the degeneracy region: for the left, middle, and the right panels, we set r ∼ 0.1, η ∼ 1.84π, x ∼ −8 and m a ∼ 3.81meV. We see that the local minimum of the degeneracy region (magenta triangles) takes place where the green, blue, red, orange and black curves meet together or go very close. LBL experiments are not sensitive to θ 12 (grey curve) and ∆m 2 21 (orange curve). As a , and ∆m 2 31 ∼ 2.52 × 10 −3 eV 2 (black) on the x − η, x − r, and η − r planes. In the upper panels, we let model parameters be the best fit values, except for those which are varied. In the lower panels we focus on the degeneracy regions, where r ∼ 0.1, η ∼ 1.84π, x ∼ −8 and ma ∼ 3.81meV. The gray curve for θ12 is below r = 0.07 so that it is hardly visible. We show the true values and the local minimum using black dots and magenta triangles respectively. result, these LBL experiments cannot exclude this region by improving precision. This also explains why once we include reactor data that is sensitive to θ 12 , the degeneracy region is excluded. One may notice that the different curves do not intersect at the magenta triangle in the x − r plane and there is the gray curve of θ 12 in the last panel for η − r. The reason is that the triangle presents a local minimum at which 11.83 > ∆χ 2 > 9. Though that is a local minimum, it does not need to cross all curves. The grey curve of θ 12 is below r=0.07, where the bottom of the panel is. We did not show it in this panel because the main feature here crossing or going close to θ 13 , θ 23 , δ and ∆m 2 31 results in the degeneracy issue for LBLs.

C. Standard oscillation parameters under tri-direct CP symmetry model
In Fig. 4, we show ∆χ 2 values against all oscillation parameters for the combination of NOνA and T2K (brown), DUNE (dahsed-blue), T2HK (dashed-green), the synergy of these four LBLs (grey), and including all LBLs and JUNO (red), assuming the tri-direct model. We see that under the tri-direct assumption, the combination of NOνA and T2K performs the worst, while DUNE performs much better, except for θ 13 and ∆m 2 21 . T2HK works slightly better than DUNE, and dominates the performance of the combination of all LBLs. The asymmetry for ∆m 2 31 comes from the asymmetry behaviour of x and r in Fig. 1   for η is passed to those for θ 13 , θ 23 and δ by the tri-direct model. We note that even LBL experiments are not sensitive to θ 12 , the uncertainty can be improved by precisely measuring other oscillation parameters within the tri-direct model. We further point out a great improvement by including JUNO data which can be seen in the result for ∆m 2 21 . (lower-left), θ13-δ(lower-right) for the synergy of these four LBLs (grey), and including all LBLs and JUNO (red). We also compare these results to those without the restrictions from tri-direct models for LBLs synergy (grey contour) and combing all experiments (red contour).
In Fig. 5, we show the points at the 3σ surface projected on θ 23 -∆m 2 31 (upper-left), θ 13 -∆m 2 31 (upper-right), θ 12 -∆m 2 21 (lower-left), θ 13 -δ(lower-right) for the synergy of these four LBLs (grey), and including all LBLs and JUNO (red). Because of the nonlinear relations between model parameters and standard parameters, the data do not spread uniformly. We also compare them with those without the restriction from the tri-direct CP model: the grey curve is for including all LBLs, while the dashed black curve is for a combination of LBLs and JUNO. There is a discontinuity when θ 13 is larger than ∼ 8.8 • , because of the degeneracy with θ 23 . These results show that assuming tri-direct CP improves the key measurements for future experiments. We see that based on one model, the better way to exclude the other one is by precision measurement of r. The experimental configuration does not affect the uncertainty for r. This uncertainty is ∼ 0.2 at 3σ under both models. Two models predict very similar values for m a . As a result, it is impossible to exclude the wrong model by measuring this parameter alone. Moreover, the uncertainty of m a depends on the model and the experimental configuration. The precision under Model A is generally better than that for Model B, expect for the combination of LBLs and JUNO data. The rank of precision from the worst to the best experimental configuration is the combination of NOνA and T2K, DUNE, T2HK, the synergy of all LBLs, and combining all LBLs and JUNO. Both m a and m s will be determined precisely by experiments. As given in the definition, the model discriminator r ≡ m a /m s points to a requirement to measure both mass squared differences in neutrino experiments as precisely as possible.

V. SUMMARY
The tri-direct CP symmetry model offers fruitful features to accommodate neutrino masses and explain neutrino mixing and oscillations. The more powerful aspect is the model predicted correlations of standard neutrino mixing parameters preserved by an underlying symmetry. We looked into a probe of the tri-direct CP symmetry model by simulating the current and future neutrino oscillation experiments, including T2K, NOνA, T2HK, DUNE and JUNO. We found that the degeneracy problem cannot be avoided at a single long baseline experiment in the precision measurement of model parameters while a combination of long-baseline and reactor experiments will resolve the problem. This fact highlights the complementarity of different neutrino oscillation experiments. In addition, we scanned the standard neutrino mixing parameters expressed by the underlying model "true" values in order to determine how powerful precision measurements in the traditional analysis will be. It seems that shape of contours in the projected parameter space can give us hints of the underlying theory but the information remains limited by a multiple-channel analysis in a single experiment. This limitation points to a combined analysis by multiple experiments with different beams and baseline configurations. Finally, we can discriminate benchmark models after a discovery of CP violation in the leptonic sector by any one of these experiments.

VI. ACKNOWLEDGEMENT
This work is supported in part by the National Natural Science Foundation of China under Grant Nos. 11505301, 11881240247, 11522546 and 11835013. JT appreciates ICTP's hospitality and nice discussions with participants during the workshop PANE2018. The work was initiated and expanded at the Chinese High-Energy Physics Conference and the MOMENT&EMuS meeting in 2018. We also thank Dr. Nick W. Prouse to kindly provide the simulation package for T2HK. Finally, we appreciate Dr. Neil Drouard Raper's help to improve the readability of our paper. In this section, we demonstrate the experimental potential for different configurations by showing the 1σ, 2σ and 3σ contours on θ 23 -δ and θ 23 -∆m 2 31 planes. These measurements are the main goals for current and future LBLs. For current running experiments NOνA (upper) and T2K (middle), we show their expected final performance in Fig. 7. For NOνA, we assume total 36 × 10 20 POT for ν andν modes until 2024, while for T2K we equally split 7.8 × 10 21 POT into two modes. We also show the combination of these two experiments in the lower panels of Fig. 7. In Fig. 8, we show the performance of DUNE (upper) and T2HK (lower). For DUNE, we consider the 3-horn-optimised design with 1.83 × 10 21 POT per year, and we adopt 3.5 years for each mode. For T2HK, we assume a 1.3 MW proton beam for the neutrino source, and run ν andν modes for 2.5 and 7.5 years respectively. More details about these experiments can be seen in Secs. III. 3σ 2σ 1σ 40