Testing violation of the Leggett–Garg-type inequality in neutrino oscillations of the Daya Bay experiment

The Leggett–Garg inequality (LGI), derived under the assumption of realism, acts as the temporal Bell inequality. It is studied in electromagnetic and strong interaction like photonics, superconducting qubits and nuclear spin. The weak interaction two-state oscillations of neutrinos affirmed the violation of Leggett–Garg-type inequalities (LGtI). We make an empirical test for the deviation of experimental results with the classical limits by analyzing the survival probability data of reactor neutrinos at a distinct range of baseline dividing energies, as an analog to a single neutrino detected at different times. A study of the updated data of the Daya Bay experiment unambiguously depicts an obvious cluster of data over the classical bound of LGtI and shows a 6.1σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6.1\sigma $$\end{document} significance of the violation of them.


Introduction
Nonclassical features of the quantum system have experienced extensive study since the inception of quantum mechanics. After a long debate between the believers of local realism and quantum mechanics, a breakthrough study, Bell's inequality (BI), was provided by Bell [1]. The unique feature of BI is its testable formula from the consequence of the famous hypothesis called local realism (LR). The LR believers assume that any observable value of an object, even if not detected, must have a definite value and that results of any individual measurement of the observables remain unaffected if they have a space-like separation. Extensive experimental investigations [2][3][4]  tested the violation of BI. These studies conclude that any local realism view of a microscopic object needs to be nonlocal. Based on these studies of BI, Leggett and Garg further derived a new series of inequalities [5] on the assumption of macrorealism (MR), now known as the Leggett-Garg inequalities (LGIs), that any system behaving as a macroscopic realism must obey. From the structure of the LGIs, we can see them as an analog of Bell's inequalities in temporal interpretation, which also makes it possible to implement a rigorous test of quantum mechanics on a macroscopic level, which is usually very difficult in designing experiments in space-like separation condition. By testing the LGtI, we can also perform a rigorously loophole-free test of quantum mechanics [2,[6][7][8].
Besides the nonlocal behavior and quantum correlation between different particles, for single particle states there can also exist entanglement by the flavor transition [9]. Due to the oscillation of neutrino flavors, it offers an ideal source to test quantum mechanics in the case of the weak interaction and from a macroscopic point of view. For two-flavor neutrino oscillation, a two-level state's matrix can be expressed in the form Here the elements of the matrix T are For a two-qubit situation, many kinds of quantum correlations [10][11][12] like entanglement, fidelity, quantum discord and geometric discord have been derived to get their general expressions. Mermin and Svetlichny [13,14] inequalities were derived for three or even many body system that have two macroscopic distinct states. Using these results of two-level states, a series of quantum-information theory calculations [15][16][17][18] have been applied on neutrino flavor oscillation in the last several years. These quantum correlations can be directly linked to the probabilities of flavor oscillation, which leads to violation of the classical boundary limits when assuming the neutrino oscillation mixing angle is not vanishing. Neutrino flavor oscillation is such a special process in that it is merely affected by the neutrino's own properties like mass square differences, mixing angles and the energies. Neutrinos just interact with matter by weak interaction with a quite low cross section. The influence of the environment on the neutrinos' propagation is much more negligible comparing with an optical or an electrical system, which makes neutrinos an ideal particle on testing the LGIs. As the mass eigenstate of a neutrino is not the same as its flavor eigenstate, during propagation, neutrinos undergo flavor mixing as regards the three flavored eigenstates. The MINOS experiment has been studied in Ref. [19], which observed the violation of Leggett-Garg-type inequalities, K3 and K4 terms, with a significance greater than 6σ [19]. The MINOS experiment is an accelerator neutrino experiment using decay in flight neutrinos with a fixed baseline distance of 735 km and a large range of ν μ energy from 0.5 to 50 GeV, which happens to cover the largest violation of LGIs K3 and K4. The Daya Bay Collaboration reported an updated data analysis of the electron anti-neutrino disappearance channel [20], which gives a best fit of sin 2 2θ 13 = 0.084 ± 0.005. We will investigate whether the Daya Bay reactor neutrino experiment can observe the violation of LGIs.

The Leggett-Garg-type inequalities
We focus on the simplest L-G-type inequality, which is constructed as follows. Consider a system with two absolutely distinguishable states corresponding to an observable quantity Q(t) which can have two different values, +1 or −1. Assume that whenever the system was being measured, the observable quantity occupies a value of either +1 or −1 for being in state 1 or 2, respectively. Then we can define a macroscopic observable Q(t) for the macroscopic system. We have C t i ,t j = Q(t i )Q(t j ) as its two-time correlation function, where Q(t i ) and Q(t j ) are the observable quantity's values when being measured at time t i and time t j . In this work we consider the two states as the survival of electron anti-neutrino and the disappearance of the electron anti-neutrino. As the neutrinos are being created in the reactor by beta decay process, they are totally in the state of the flavor eigenstate. Since the PMNS matrix does not change with time, the two-flavor neutrino oscillation obeys the same survival probability. We shall introduce this stationarity assumption [21], which requires that the evolution of the neutrino for different ordered time intervals is the same. Then Next, consider a sequence of times t 1 , t 2 , t 3 and t 4 (here, t 1 < t 2 < t 3 < t 4 ). If we take a series of measurements for Q(t) in these four times, it is straightforward to determine four time correlations (C 12 , C 23 , C 34 and C 14 ). Then it is possible to adopt the stationary condition on the standard LGI procedure leading to K4 LG-type inequality involving four correlation functions. For any sequence of measurements, any Q(t i ) has a definite observable value, regardless of the choice of the pair lies always between -2 and +2. Similarly, the K 3 inequality lies between -1 and +1. If all the terms in the above formula are replaced by time correlations (average), the Leggett-Garg-type inequalities are in the form (2) The above inequalities impose a constraint on macroscopic realism as regards the temporal separated joint probabilities in any two-state system.

Three-flavoured neutrino oscillations
It has extensively been verified that the flavor component of a neutrino oscillates during its propagation. The oscillation properties of different neutrino flavors are determined by their mixing angles (θ 12 , θ 23 , and θ 13 ), a CP phase of the Pontecorvo-Maki-Nakagawa-Sakata matrix and their masssquared differences ( m 2 32 , m 2 21 ) [22,23]. Here, we will treat the Leggett-Garg-type inequalities using the updated measurement of ν e survival channel in the results of the Daya Bay experiment, where give the latest best fit of the mixing angle θ 13 with large significance [20]. The measurement was updated later with a full detector configuration [20]. For the θ 13 measurement, one used the baseline length dividing the ν e energy as the variable to depict the survival probability of ν e as P ν e →ν e = 1 − cos 4 θ 13 sin 2 2θ 12 sin 2 where E is the energy of ν e in MeV, L is the propagation distance between the near and far point detector, θ 12 is the solar neutrino mixing angle and m 2 21 is their mass-squared difference in eV 2 . Notice that m 2 ee is an effective mass-squared difference [24] in electron anti-neutrino disappearance with the form m 2 ee = cos 2 θ 12 m 2 31 + sin 2 θ 12 m 2 32 .
Since m 2 12 = (7.50±0.20)×10 −5 eV 2 , while m 2 ee = (2.42± 0.11) × 10 −3 eV 2 according to Ref. [20], we can choose an appropriate value of the ratio L/E to make one of the terms of sin 2 1.267 m 2 L E vanishing. For the Daya Bay experiment, the effect of the parameter θ 12 becomes far less, sufficiently so that, compared with θ 13 , θ 12 can be regarded as negligibly small. Given that there is an initial pure electron anti-neutrino source, after propagation for a time t, the survival probability of ν e will be P ν e →ν e = 1 − sin 2 2θ 13 sin 2 1.267 m 2 ee ct E .
However, there is the MSW (Mikheyev-Smirnov-Wolfenstein) effect (usually called the matter effect) during the propagation of neutrino in matter. The effect is only significant for high energy neutrinos and a long range of matter, like the solar neutrino experiment. The KamLAND and Super-K P ee day-night discrepancies are only obvious for larger than 6 MeV neutrinos [25,26]. Furthermore, the solar neutrino experiments involve the matter effect caused by the electron in the solar medium, which electron density is much larger than that in the Earth. Generally speaking, a neutrino vector of state in flavor basis |ν(t) = (ν e (t) ν μ (t) ν τ (t)) T obeys the Schrödinger equation: where the Hamiltonian can be replaced by an effective one as where V is the effective charged potential contribution to ν e [27], given in the form where ρ(x) is the matter density along the track path of the neutrino,Y e (x) (for the Earth 0.5) is the number of electrons normalized to the number of nucleons. For the matter of constant density, the series expansion for three-flavor neutrino oscillation probabilities can be derived from the Hamiltonian Eq. (7) [28]. For ν e survival, the survival probability expansion to second order is where α = [GeV] , For this Daya Bay analysis, we calculate the discrepancy of the P ee probability of a 6 MeV neutrino in the range of 0 to 10 km covering the range of the experiment, about 2 km. From Fig. 1, we can draw the conclusion that the matter effect is too small to be included in the "short" baseline. According to the expansion of the P ee with matter effect mentioned above, we will use the oscillation probability just in the vacuum. Using the stationary assumption, one can derive the four joint probabilities P ν α ,ν β (t i , t j ), here α and β stands for ν e and another neutrino flavor, i and j are from 1 to 4 defined above. The two-time correlation function C 12 in this P ν e →ν e is given by Similarly, the correlation functions C 23 , C 34 , and C 14 can be calculated. Using the Eq. (5), the quantity K Q n can be evaluated as defined in Eq. (2). By choosing the time intervals in a particular way, we can achieve a maximum value of K 4 when t 4 − t 3 = t 3 − t 2 = t 2 − t 1 = δt. Under this condition, the correlation functions depend on the baseline length L and the neutrino energy E ν . We select the neutrinos' measured L e f f /E to make the oscillation phase ψ a = 1.267 m 2 E c(δt) obey the sum rule: ψ 12 + ψ 23 + ψ 34 = ψ 14 . We have an experimental arrangement in which measurements occur at some fixed distance from the neutrino sources. Assuming the neutrino begins in the pure |ν e state, Here n can be 3 or 4 in this paper, which corresponds to K 3 or K 4 LGI. In quantum mechanics, the commutators of operators can be nonvanishing. However, in a classical system, operators with observable values must commute; then the macrorealism derived K n will become The Daya Bay Collaboration released updated oscillation results as a function of the effective baseline distance L eff over the average energy E ν in bins [20]. For their six antineutrino detectors (ADs) placed in three separate experimental halls (EHs) and three nuclear reactors neutrino sources, the effective baseline varies for each detected anti-neutrino. The Daya Bay experiment covers an energy between 1 and 8 MeV. The ranges of effective baseline and energy correspond to a phase range of (0, 3/4π ), within which the violations of LGI will be observed near the minimum point of the antineutrino survival probability.
To test the violations of the K 3 and K 4 inequalities, we address the data from the Daya Bay neutrino experiment. The Daya Bay experiment extracted the survival probabilities of neutrinos using Daya Bay and Ling-Ao nuclear power stations' reactors. We use all the measurement positions including EH1, EH2 and EH3 of the Daya Bay. The reactors provide different sources of neutrinos with several fixed baselines and an energy spectrum with peaks. We make a θ 13 fit over the Daya Bay updated data and get the fit error band and center value of P ee shown in Fig. 2. With the best fit of sin 2 2θ 13 and the 1σ error band of it, we generate a large set of pseudodata. Then we select all sets of data points in Fig. 3 which obey the sum rule of phase with the precision of 0.5% (ψ 1 + ψ 2 ∈ ψ 3 ) and 0.1% (ψ 1 + ψ 2 + ψ 3 ∈ ψ 4 ) for the K 3 and K 4 respectively. For the K 3 (K 4 ) situation, 48 (56) correlation triples (quadruples) satisfy the sum rule. Meanwhile the updated measurement only includes the static errors and we simply assumed that the errors at small phase of the oscillation probability are the fitting error.
The violation of Leggett-Garg-type inequalities has been tested and confirmed by the MINOS experiment, with the K 3 and K 4 being inconsistent with the realism prediction over 5σ [19]. Since the violation of Leggett-Garg-type inequalities happens when the mixing angle of two flavors is not zero, we suppose that the violation could be observed in the ν e survival channel at Daya Bay. In order to estimate the significance from events as regards the number of violations, we simulated the statistical quantity by creating a large sample of pseudodata based on the fitting result of the observed P ee Fig. 2 The electron anti-neutrino survival probability versus effective propagation distance L eff over anti-neutrino energy E ν . The dashed blue curve indicates the prediction fitting directly to the measured Daya Bay values of P ee . The red band indicates a 1σ confidence interval around the fitted prediction. The blue triangles, green rectangles and black dots are the binned data of the Daya Bay EH1, EH2 and EH3 from Ref. [20], respectively values. The pseudodata is generated by a Gaussian distribution model with the means and variances matched to the center values and deviations of the best fit. Each set of simulated data gives an artificial number of LGI violations for K 3 and K 4 , from which we can calculate the level of inconsistency of the predictions between quantum and classical K n .
To estimate the confidence level of these results being inconsistent with expectation and prediction from realism, we make a fit of the histogram filled by predicted LGI violations number under the realism model of Eq. (7) to a betabinomial distribution, thus to estimate the deviation of classical predictions from the actually observed number of LGI violations. For the actual number of LGI violations (41 in 48 data points), there exists a 6.1σ deviation from the expected distribution of the classical prediction.
A similar statistical test is made for LGI K 4 . Using the filter of the phase sum rule described above, we get a number of 30 (in total of 56 data points) exceeding the classical limits. As Fig. 4 shows, there are obvious clusters of points over the classical bound of K 3 and K 4 . The discrepancy between the observed events number and the classical predicted events originating from the fluctuation is very clear. Our K 4 data also possesses a 6σ deviation from the classical prediction.

Discussions
The results mentioned above clearly constrain the validity of quantum mechanics in such a macroscopic area. Values of LGtI K 3 and K 4 are violated with the QM prediction at a confidence level of over 6σ compared with the classical bound Fig. 3 The histograms of the number of K3 (upper) and K4 (lower) values that violate the LGI bound. The left curves with red filling indicate the expected classical distributions, while the right curves with blue filling indicate the corresponding quantum quantity for the neutrino θ 13 mixing in our estimation. Anti-electron neutrino oscillations also violate the limits of Leggett-Garg inequality. The detected violations act as a new affirmation of quantum nonlocality existing in the neutrino system during its long-range propagation. These violations were observed over the near and far detectors placed at three experimental halls (EHs) with the baseline long enough to make the test not being a Bell-like inequality test. Besides, it should be worthwhile to make a detailed data analysis on the Daya Bay experiment involving three-flavor neutrino oscillation, in order to achieve more data points of LGI K 3 and K 4 . It could be worth to test the quantum mechanics in such a weak interaction context. Although tests of incompatible of LGI and QM have been achieved in photonics and electronic experiments [29,30], nuclear spin qubits [31] and even condensed states [32], there are few reports of LGI violation in particle physics. Even though the MINOS and Daya Bay experimental setup show the LGtI violations, these two experiments are all in the context of two-flavor neutrino oscillation, which can not reveal the CP violation. Since entanglement exists between a pair of neutral meson and anti-meson, which will violate the Bell inequality [33], three-flavor oscillation analysis involving neutrinos and anti-neutrinos may shed light on the study of CP-violating phase. [34] Fig. 4 The distribution of K3 (upper) and K4 (lower) versus the effective propagation length divided by neutrino energy reconstructed from P ee . The black dot data show a cluster over the LGI bound. We also show the expected distributions of classical (red circles) and quantum (blue circle) predictions. Note that the K 3 and K 4 can be multiple values, since there are many triples and quadruples satisfying the phase sum rule