1 Introduction

Neutrino is a Standard Model neutral weakly interacting fermion [1]. It is the second most abundant particle in the Universe after photons of light. Neutrino oscillation (NO) implies that the neutrino has a non-zero mass. In the framework of the simplest standard model of three-neutrino mixing, three different flavors of neutrino are electron e, muon \(\mu \), and tau \(\tau \) leptons, in which the three flavor states are unitary linear combinations of three mass eigenstates [2, 3]. NO shows that a given flavor may change into another flavor in the neutrino propagation. The probability of measuring a particular flavor for a neutrino varies periodically as it propagates through space, and can be measured at the arbitrary time. The values of the oscillations parameters have been measured and analyzed in both theory and experiment in recent years [4,5,6,7,8]. Remarkably, oscillation probabilities of neutrino can be used to study the different properties from classical to quantum mechanical prediction of such an interesting system.

As an analog of Bell’s inequality in temporal interpretation, LGI study the correlations of a single system measured at different times, based on two assumptions of macroscopic realism and non-invasive measurability [9,10,11]. It is shown that experimentally-observed neutrino oscillations can violate the classical limits imposed by the LGI [12,13,14,15], provide an evidence that quantum coherence can apply broadly to microscopic systems over a macroscopic distance. However, the violation of LGI can be taken as an effective indicator for quantifying the amount of quantumness in the framework of quantum resource theories (QRTs). Recently, the quantification of quantumness in terms of flavor oscillation probabilities of NOs was investigated by quantities in QRTs, including entanglement [16,17,18], Svetlichny inequalities [19, 20], entropic uncertainty relation [21, 22], quantum coherence [23,24,25], and quantum correlations [26,27,28,29].

Among them, entanglement is the most fundamental concept that can be rigorously quantified and characterized by the tools of QRTs. It has many potential applications in quantum information processing, including quantum cryptography [30, 31], quantum teleportation [32], entanglement swapping [33], and so on. In particular, multipartite entanglement turns out to be a versatile tool in quantum error correction [34]. There are many measures to quantify the entanglement in multiqubit quantum systems, such as the generalized geometric measure (GGM) [35, 36], three-\(\pi \) entanglement measure [37], the genuinely multipartite concurrence (GMC) [38], and concurrence fill [39]. GGM is identified as an optimized distance of the given state from the set of all states that are closest biseparable states. The three-\(\pi \) entanglement measure consider all the bipartite residual entangle quantify three-qubit entanglement, based on the negativity. GMC is an intuitive measure of multipartite entanglement, based on the concurrence. Its lower bound can be obtained by the powerful detection criteria. Concurrence fill is introduced as a genuine tripartite entanglement measure in terms of the area of the concurrence triangle. It represents a global entanglement, which fulfills the following properties: non-negativity, monotonicity, discriminance, normalization, smoothness, and convexity [40]. For the GHZ state, the concurrence fill can reach its maximum of 1, since the lengths of the three edges of concurrence triangle are all equal to 1. For the W state, its maximum is \(F_{123}=8/9\approx 0.89\). In particular, \(F_{123}^{max}=0.89\) for three-flavor electron NOs considered in this work. Investigations of genuine multipartite entanglement of neutrino oscillation might pave the way for the further quantum computation using oscillating neutrinos.

Here, we quantify the quantumness in three-flavor NOs by several tripartite entanglement measures, including GGM, three-\(\pi \), GMC, and concurrence fill in QRTs, and compare their performances with respect to the ratio between the distance and neutrino energy L/E. The definitions of concurrence fill and three-\(\pi \) take all bipartite entanglements into account, and thus their change are always smooth. For the three-flavor electron and muon neutrino oscillations, concurrence fill is always larger than three-\(\pi \), suggesting that concurrence fill contains more quantum resource. The definitions of GGM and GMC are given by the minimal entanglement measures of bipartite entanglement, which lead to several nonanalytical sharp peaks. Moreover, we use the concurrence triangles for three-qubit systems to visualize the difference between GMC and concurrence fill. It is shown that concurrence fill is always larger than or equal to the GMC in both sources of NOs. The results show that concurrence fill is a genuine tripartite measure of tripartite entanglement in three-flavor NOs.

The paper is arranged as follows. In Sect. 2, we introduce three-flavor neutrino model. In Sect. 3, we briefly review several tripartite entanglement measures, including GGM, three-\(\pi \), GMC, and concurrence fill. In Sects. 4 and 5, we compare the performances of the four tripartite entanglement measures in the three-flavor electron and muon antineutrino oscillations, respectively. Finally, a summary is concluded in Sect. 6.

2 The three-flavor neutrino model

In the standard three-flavor neutrino oscillation model, a linear superpositions of the mass eigenstates \(|\nu _1\rangle \), \(|\nu _2\rangle \), \(|\nu _3\rangle \), constitutes the three flavor states \(|\nu _e\rangle \), \(|\nu _\mu \rangle \), and \(|\nu _\tau \rangle \):

$$\begin{aligned} |\nu _{\alpha }\rangle =\sum _{k}U^*_{\alpha k}|\nu _{k} \rangle , \end{aligned}$$
(1)

where \(k=1,2,3\) and \(\alpha =e,\mu ,\tau \), which represents the neutrino flavor state. The related \(U_{\alpha k}\) is a \(3\times 3\) unitary matrix, is called as the Pontecorvo–Maki–Nakagawa–Sakata matrix [41]. It can be characterized by three mixing angles and a charge conjugation and parity (CP) violating phase,

$$\begin{aligned} U\!\!=\!\left( \!\! \begin{array}{ccc} c_{12}c_{13} &{}\quad \! s_{12}c_{13} &{}\quad \! s_{13}e^{-i\delta _{CP}} \\ -\!s_{12}c_{23}\!-\!c_{12}s_{13}s_{23}e^{i\delta _{CP}} &{}\quad \!c_{12}c_{23}\!-\!s_{12}s_{13}s_{23}e^{i\delta _{CP}} &{}\quad \! c_{13}s_{23}\\ s_{12}s_{23}\!-\!c_{12}s_{13}c_{23}e^{i\delta _{CP}} &{}\quad \! -\!c_{12}s_{23}\!-\!s_{12}s_{13}c_{23}e^{i\delta _{CP}} &{}\quad \! c_{13}c_{23} \\ \end{array} \!\!\right) , \end{aligned}$$
(2)

where \(c_{ij}\equiv \cos \theta _{ij}\) and \(s_{ij}\equiv \sin \theta _{ij}\) \((i,j=1,2,3)\). The CP violating phase has not yet been observed, so we neglect it for simplicity. The state \(|\nu _{k}(t)\rangle \) is the mass eigenstate of the free Dirac Hamiltonian H with an positive energy \(E_k\), and its time evolution of the mass eigenstates \(|\nu _{k}\rangle \) is expressed as

$$\begin{aligned} |\nu _{k}(t)\rangle =e^{-iE_kt/\hbar }|\nu _{k}(0)\rangle . \end{aligned}$$
(3)

Substituting the Eqs. (1) and (2) into Eq. (3), we can get the time evolution of flavor neutrino states \(|\nu _{\alpha }(t)\rangle \) as

$$\begin{aligned} |\nu _{\alpha }(t)\rangle =a_{\alpha e}(t)|\nu _{e}\rangle +a_{\alpha \mu }(t)|\nu _{\mu }\rangle +a_{\alpha \tau }(t)|\nu _{\tau }\rangle , \end{aligned}$$
(4)

where \(a_{\alpha \beta }(t)= \sum _kU^*_{\alpha k}e^{-iE_kt/\hbar }U_{\beta k}\). Finally, the transition probability for detecting \(\beta \) neutrino in the original \(\alpha \) neutrino state is given by [42],

$$\begin{aligned} P_{\alpha \beta }= & {} \delta _{\alpha \beta }-4\sum _{k>l}\text {Re}(U_{\alpha k}^*U_{\beta k}U_{\alpha l}U_{\beta l}^*)\sin ^2\left( \varDelta m^2_{kl}\frac{Lc^3}{4\hbar E}\right) \nonumber \\&+2\sum _{k>l}\text {Im}(U_{\alpha k}^*U_{\beta k}U_{\alpha l}U_{\beta l}^*)\sin \left( \varDelta m^2_{kl}\frac{Lc^3}{2\hbar E}\right) , \end{aligned}$$
(5)

where \(\varDelta m^2_{kl}=m^2_{k}-m^2_{l}\), E is the energy of the neutrino in neutrino experiments, and \(L\thickapprox ct\) (c is the speed of light in free space) is the distance between the source and the detector.

For convenience, we can write the oscillatory quantity \(\sin ^2\left( \varDelta m^2_{kl}\frac{Lc^3}{4\hbar E}\right) \) that appears in Eq. (5) as [1],

$$\begin{aligned} \sin ^2\left( \varDelta m^2_{kl}\frac{Lc^3}{4\hbar E}\right) =\sin ^2\left( 1.27\varDelta m^2_{kl}[eV^2]\frac{L[km]}{E[GeV]}\right) . \end{aligned}$$
(6)

Meantime, the oscillation parameters in normal ordering of the neutrino mass spectrum (\(m_1< m_2 < m_3\)) are chosen as

$$\begin{aligned} \varDelta m_{21}^2= & {} 7.50 \times {10^{ - 5}}e{V^2},\nonumber \\ \varDelta m_{31}^2= & {} 2.457 \times {10^{ - 3}}e{V^2},\nonumber \\ \varDelta m_{32}^2= & {} 2.382 \times {10^{ - 3}}e{V^2},\nonumber \\ {\theta _{12}}= & {} {33.48^ \circ }, {\theta _{23}} = {42.3^ \circ }, {\theta _{13}} = {8.50^ \circ }. \end{aligned}$$
(7)

The occupation number in the neutrinos used here can be established in the following correspondence [16]

$$\begin{aligned} \left| {{\nu _e }} \right\rangle&\equiv {\left| 1 \right\rangle _e} \otimes {\left| 0 \right\rangle _\mu } \otimes {\left| 0 \right\rangle _\tau }&\equiv \left| {100} \right\rangle , \end{aligned}$$
(8)
$$\begin{aligned} \left| {{\nu _\mu }} \right\rangle&\equiv {\left| 0 \right\rangle _e} \otimes {\left| 1 \right\rangle _\mu } \otimes {\left| 0 \right\rangle _\tau }&\equiv \left| {010} \right\rangle ,\end{aligned}$$
(9)
$$\begin{aligned} \left| {\nu _\tau } \right\rangle&\equiv {\left| 0 \right\rangle _e} \otimes {\left| 0 \right\rangle _\mu } \otimes {\left| 1 \right\rangle _\tau }&\equiv \left| {001} \right\rangle . \end{aligned}$$
(10)

Thus, the flavor oscillations of neutrino can be seen as the time evolution of a tripartite quantum state. Then, from Eq. (4), we have

$$\begin{aligned} {\left| {\psi _\alpha (t)} \right\rangle } = {a_{\alpha e}}(t)\left| {100} \right\rangle + {a_{\alpha \mu }}(t)\left| {010} \right\rangle + {a_{\alpha \tau }}(t)\left| {001} \right\rangle . \end{aligned}$$
(11)

As a result, we can use the framework of quantum resource theory to study the flavor oscillations of neutrino. In the following, we focus on the performances of the tripartite entanglement measures, including GGM, three-\(\pi \), GMC, and concurrence fill, and give a detailed comparison among them in this three-flavor neutrino systems.

3 Measure of tripartite entanglement

3.1 Generalized geometric measure

The GGM [35, 36] of an N-partite pure state \(\left| \psi _{N} \right\rangle \), can measure the entanglement of pure quantum states of an arbitrary number of parties. It is defined as an optimized distance of the given state from the set of all states that are not genuinely multiparty entangled, given by

$$\begin{aligned} \varepsilon (\left| \psi _{N} \right\rangle )=1- \varLambda ^2_{max}(\left| \psi _{N} \right\rangle ), \end{aligned}$$
(12)

where \(\varLambda _{\max }(\left| \, \psi _{N} \,\right\rangle )=\max |\langle \chi \,|\,\psi _{N}\rangle |= \mathop {\max } \limits _{\left| \chi \right\rangle } F (\left| \psi _N\right\rangle , \left| \chi \right\rangle )\). The maximization is performed over all \(\left| \chi \right\rangle \) that are not multiparty entangled. \(F(\left| \psi _N\right\rangle , \left| \chi \right\rangle )\) denotes the fidelity between two pure states \(\left| \psi _N\right\rangle \) and \(\left| \chi \right\rangle \). Based on the Hilbert-Schmidt distances, we can obtain an equivalent mathematical expression of Eq. (12)

$$\begin{aligned} G(\left| \psi _{n}\right\rangle )&=1-max\{\lambda ^2_{I:L}|I\cup L=\{A_1, \ldots , A_N\},\nonumber \\ I \cap L&= \emptyset \}, \end{aligned}$$
(13)

where \(\lambda _{I:L}\) is the maximal Schmidt coefficient in the bipartite split I : L of \(\left| \psi _{N} \right\rangle \).

3.2 The three-\(\pi \) entanglement

For a three-partite pure state \(\left| \psi \right\rangle _{ABC}\), the Coffman–Kundu–Wootters-inequality-like [43] monogamy inequality by the negativity quantified the entanglement is given by

$$\begin{aligned} N^2_{AB}+N^2_{AC}\le N^2_{A(BC)}, \end{aligned}$$
(14)

where \(N_{AB}\) and \(N_{AC}\) are the sum of the negative eigenvalues of the partial transpose of the states \(\rho _{AB}=\text {Tr}_C(\rho _{ABC})\) and \(\rho _{AC}=\text {Tr}_B(\rho _{ABC})\), respectively. Following the idea of Ref [37], we have

$$\begin{aligned} N_{A(BC)}= C_{A(BC)}=\sqrt{2[1-\text {Tr}(\rho ^2_A)]}, \end{aligned}$$
(15)

where \(\rho _A=\text {Tr}_{BC}(\rho _{ABC})\). The \(C_{A(BC)}\) is defined as a one-to-other bipartite entanglement, i.e., the entanglement between qubit A and the remaining two (B and C) taken together as an “other” single party, where a subscript i \((i=A, B, C)\) refers to the system’s ith qubit. The residual entanglement is defined as the difference between the two sides of Eq. (14)

$$\begin{aligned} \pi _A=N^2_{A(BC)}-N^2_{AB}-N^2_{AC}. \end{aligned}$$
(16)

Similarly, if one takes the different focus B and C, the other two residual entanglement are created

$$\begin{aligned}&\pi _B=N^2_{B(AC)}-N^2_{BA}-N^2_{BC}, \end{aligned}$$
(17)
$$\begin{aligned}&\pi _C=N^2_{C(AB)}-N^2_{CA}-N^2_{CB}, \end{aligned}$$
(18)

respectively. Note that the residual entanglement for the different focus changes under transformations of the qubits. Finally, the three-\(\pi \) entanglement of the tripartite systems are obtained as the average of \(\pi _A\), \(\pi _B\), \(\pi _C\)

$$\begin{aligned} \pi _{ABC}=\frac{1}{3}(\pi _A+\pi _B+\pi _C). \end{aligned}$$
(19)

3.3 Genuinely multipartite concurrence

The GMC is a computable measure quantifying the amount of multipartite entanglement based on the well-known concurrence [38]. For n-partite pure states \(\left| \varPsi \right\rangle \in H_1 \otimes H_2 \otimes \cdot \cdot \cdot \otimes H_n\), where \(\mathrm{dim}(H_i)=d_i\), \(i=1,2,\ldots ,n\), the GMC is defined as

$$\begin{aligned} {C_{GMC}}(\left| \varPsi \right\rangle )= \mathop {\min }\limits _{{\gamma _i} \in \gamma }\sqrt{2[1-\text {Tr}(\rho ^2_{A_{\gamma _i}})]} \end{aligned}$$
(20)

where \(\gamma =\{\gamma _i\}\) represents the set of all possible bipartitions \(\{A_i|B_i\}\) of \(\{1,2,\ldots n.\}\). The GMC can also be generalized to the case of mixed states via the convex roof construction

$$\begin{aligned} {C_{GMC}}(\rho )=\mathop {\inf }\limits _{\{p_i, \left| \psi _i \right\rangle \} }\mathop {\sum }\limits _{i }p_iC_{GMC}(\left| \psi _i \right\rangle ), \end{aligned}$$
(21)

where the infimum is taken over all possible decompositions \(\rho =\sum _i\left| \psi _i \right\rangle \left\langle \psi _i \right| \). For three-qubit systems, GMC is exactly the square root of the length of the shortest edge of the concurrence triangle. Following the idea of Ref. [39], hereinafter we also ignore the square root and treat \(C_{GMC}\) as the length of the shortest edge for simplicity.

3.4 The concurrence fill

In 2021, Xie et al. introduced the method of concurrence fill to genuinely capture triangle measure of tripartite entanglement, based on Heron’s formula. Compared to other tripartite entanglement measures, concurrences fill has the two advantages: (i) it contains more information; (ii) it is always smooth, while other measures contain a minimum argument, which will leads to a nonanalytical sharp peaks. In their paper, a concurrence triangle is presented where the three squared one-to-other concurrences are the lengths of the three edges of a triangle for a three-qubit system. Then concurrence fill is defined as the square root of the area of the concurrence triangle. This gives the following expression of genuine tripartite entanglement measure [39]

$$\begin{aligned} F_{ABC}\equiv \left[ \frac{16}{3}Q(Q-C^2_{A(BC)})(Q-C^2_{B(AC)})(Q-C^2_{C(AB)})\right] ^{1/4}, \end{aligned}$$
(22)

where

$$\begin{aligned} Q=\frac{1}{2}(C^2_{A(BC)}+C^2_{B(AC)}+C^2_{C(AB)}), \end{aligned}$$
(23)

with \(C^2_{A(BC)}\), \(C^2_{B(AC)}\), and \(C^2_{C(AB)}\) being \(2[1-\text {Tr}(\rho ^2_A)]\), \(2[1-\text {Tr}(\rho ^2_B)]\), and \(2[1-\text {Tr}(\rho ^2_C)]\), respectively. Q is the half-perimeter, the prefactor 16/3 ensures the normalization \(0\le F_{ABC} \le 1\), and the extra square root beyond Heron’s formula guarantees local monotonicity under local quantum operations assisted with classical communications.

4 Entanglement in electron antineutrino oscillations

When the electron neutrino is generated at initial time \(t=0\), the evolutive states for three-flavor NOs is given by

$$\begin{aligned} {\left| {\psi _e (t)} \right\rangle } = {a_{ee}}(t)\left| {100} \right\rangle + {a_{e\mu }}(t)\left| {010} \right\rangle + {a_{e\tau }}(t)\left| {001} \right\rangle . \end{aligned}$$
(24)

To quantify the tripartite entanglement, we shall focus on the the density matrix \(\rho _{ABC}^e(t)={\left| {\psi _e (t)} \right\rangle }\left\langle \psi _e (t)\right| \), that is,

$$\begin{aligned} \rho _{ABC}^e(t) = \left( {\begin{array}{*{20}{c}} 0 &{}\quad 0 &{}\quad 0 &{}\quad 0 &{}\quad 0 &{}\quad 0 &{}\quad 0 &{}\quad 0\\ 0 &{}\quad {\rho _{22}^e} &{}\quad {\rho _{23}^e} &{}\quad 0 &{}\quad {\rho _{25}^e} &{}\quad 0 &{}\quad 0 &{}\quad 0\\ 0&{}\quad {\rho _{32}^e}&{}\quad {\rho _{33}^e}&{}\quad 0&{}\quad {\rho _{35}^e}&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad {\rho _{52}^e}&{}\quad {\rho _{53}^e}&{}\quad 0&{}\quad {\rho _{55}^e}&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0 \end{array}} \right) , \end{aligned}$$
(25)

where the matrix elements are written as \(\rho _{22}^e\! =\!{\left| {{a_{e\tau }}(t)} \right| ^2};\) \(\rho _{23}^e(t)\! =\! {a_{e\tau }}(t)a_{e\mu }^ * (t);\) \(\rho _{25}^e(t) \!= \!{a_{e\tau }}(t)a_{ee}^ * (t);\) \(\rho _{32}^e \!= \!{a_{e\mu }}(t)a_{e\tau }^ * (t);\) \(\rho _{33}^e(t) \!= \!{\left| {{a_{e\mu }}(t)} \right| ^2};\) \(\rho _{35}^e(t) \!= \!{a_{e\mu }}(t)a_{ee}^ * (t);\) \(\rho _{52}^e \!=\! {a_{ee}}(t)a_{e\tau }^*(t);\) \(\rho _{53}^e(t) \!= \!{a_{ee}}(t)a_{e\mu }^*(t);\) \(\rho _{55}^e(t)\! = \!{\left| {{a_{ee}}(t)} \right| ^2}.\)

And the corresponding oscillation probabilities, \(P_{ee}(t)=|a_{ee}(t)|^2\), \(P_{e\mu }(t)=|a_{e\mu }(t)|^2\) and \(P_{e\tau }(t)=|a_{e\tau }(t)|^2\), in electron NOs as a function of L/E are plotted in Fig. 1. When the initial neutrino is in electron flavor neutrino, the survive probability for maintaining in this flavor is always higher than 0.1, and the other two transition probabilities are smaller than 0.7 in a range [0, 40] of L/E with dimension \(\mathrm{km/MeV}\). The differences between survive and transition probabilities are large over a wide range with respect to L/E.

Fig. 1
figure 1

Oscillation probabilities as a function of the ratio between the traveled distance L and the energy E. The figure plots the oscillation probability \({v_e}(0) \rightarrow {v_e}(t)\) (red, solid line), \({v_e}(0) \rightarrow {v_\mu }(t)\) (purple, dashed line), \({v_e}(0) \rightarrow {v_\tau }(t)\) (blue, dashed-dotted line) when the initial neutrino flavor is electron flavor

For the three-flavour neutrino state \({\left| {\psi _e (t)} \right\rangle }\) , using the Eq. (13), the GGM can be calculated as

$$\begin{aligned} G(\rho ^e)=1-max\{{\lambda ^e _A},\, {\lambda ^e _B},\, {\lambda ^e _C}\}, \end{aligned}$$
(26)

where \(\lambda ^e_A\), \(\lambda ^e_B\) and \(\lambda ^e_C\) are the largest eigenvalues of the reduced density matrix \(\rho ^e_A\), \(\rho ^e_B\), \(\rho ^e_C\), respectively. It can found that \(\lambda ^e_A=max\{P_{e\tau }+P_{e\mu }, P_{ee} \}\), \(\lambda ^e_B=max\{P_{ee}+P_{e\tau }, P_{e\mu }\}\), \(\lambda ^e_C=max\{P_{ee}+P_{e\mu }, P_{e\tau } \}\), respectively.

To get the three-\(\pi \) entanglement of the three-qubit electron neutrino system, we should firstly trace over the qubit C and B of \(\rho _{ABC}\) to get the reduced density matrices: \(\rho _{AB}\) and \(\rho _{AC}\), respectively. The partial transpose \(\rho _{AB}^{T}\) and \(\rho _{AC}^{T}\) are obtained by focus on A. By solving the corresponding eigen equation, we can find the negativity of \(N_{AB}\) and \(N_{AC}\) with their negative eigenvalues. Using Eqs.(15) and Eqs.(16), we obtain the residual entangle \(\pi _A\). In a similar way, we have \(\pi _B\) and \(\pi _C\). Finally, the systemic three-\(\pi \) entanglement can be calculated as

$$\begin{aligned}&\pi _{ABC}(\rho ^e)\nonumber \\&\quad =\frac{4}{3}\Bigg [-P^2_{ee}-P^2_{e\mu }-P^2_{e\tau }+P_{ee}\sqrt{P^2_{ee}+4P_{e\mu }P_{e\tau }}\nonumber \\&\qquad +P_{e\mu }\sqrt{P^2_{e\mu }+4P_{ee}P_{e\tau }}+P_{e\tau }\sqrt{P^2_{e\tau }+4P_{ee}P_{e\mu }}\,\Bigg ]. \end{aligned}$$
(27)

For GMC in the three-flavour electron neutrino system, from the Eq. (20), the GMC is expressed as \(\text {min}\big \{2[1-\text {Tr}(\rho ^2_A)]\), \(2[1-\text {Tr}(\rho ^2_B)]\), \(2[1-\text {Tr}(\rho ^2_C)]\big \}\). This gives

$$\begin{aligned} C_{GMC}(\rho ^e)&=\mathrm{min}\,\big \{4P_{ee}\big (P_{e\mu }+P_{e\tau }\big ),\, 4P_{e\mu }\big (P_{ee}+P_{e\tau }\big ), \nonumber \\&\quad 4P_{e\tau }\big (P_{ee}+P_{e\mu } \big )\big \}. \end{aligned}$$
(28)
Fig. 2
figure 2

The four kinds of multipartite entanglement measure, including GGM (black dotted line), three-\(\pi \) (blue, dashed-dotted line), GMC (purple, dashed line), and concurrence fill (red, solid line) for the electron antineutrino oscillation system

Finally, by calculating the area of the concurrence triangle with the Eq. (22), we can get concurrence fill

$$\begin{aligned} F_{123}(\rho ^e)=8\Bigg \{\frac{P^2_{ee}P^2_{e\mu }P^2_{e\tau }[P_{e\tau }P_{e\mu }+P_{ee}(P_{e\mu }+P_{e\tau })]}{3}\Bigg \}^{1/4}. \end{aligned}$$
(29)

We can find that the amount of tripartite quantum entanglement depends on the three oscillation probabilities, which are subjected to the normalization constraint: \(\sum _{\alpha }P_{\alpha \beta }=\sum _{\beta }P_{\alpha \beta } =1\) (\(\alpha ,\beta =e,\mu ,\tau \)). The maximal value of four tripartite entanglement measures: GGM, three-\(\pi \), GMC, and concurrence fill, in three-flavored electron NOs is 1/3, \(4(\sqrt{5}-1)/9\), 8/9, and 8/9, respectively, when all the oscillation probabilities are equal to 1/3.

In Fig. 2, we plot the tripartite entanglement measures for three-flavored electron NOs as a function of ratio L/E, including GGM, three-\(\pi \), GMC, and concurrence fill. At the point \(L/E=0\), all the measure of multipartite entanglement are 0. In particular, each edge of the concurrence triangle has a zero length, so the area of the concurrence triangle is also zero. Similarity, the four measures repeat a same change trend in one period of L/E that increase firstly and then decrease. However, they take different amount of information in the neutrino propagation. Moreover, the three oscillation probabilities get close to the critical condition that they take the same value, i.e., \(P_{\alpha \beta }=1/3\), at around \(L/E = 10.83~\mathrm{km/MeV}\). In this sense, the maximal value of entanglement for these measures are approximatively 0.32, 0.55, 0.88, and 0.89, respectively. One can see that GMC and GGM have nonanalytical sharp peaks due to the nonanalytic minimum argument in its expression, concurrence fill is always smooth and contains more information. Although the three-\(\pi \) presents a smooth curve, concurrence fill is always larger than it in the whole range of L/E. These results show that concurrence fill has advantages over other measures in the three-flavored electron NOs.

In addition, to visualize the difference between GMC and concurrence fill for the electron NO systems, we make the use of the concurrence triangles to characterize the amount of information they carried. The concurrence triangles for three cases: \(P_{e\mu }=P_{e\tau }=0.115\), \(P_{ee}=0.77\) at around \(L/E = 4.61~\mathrm{km/MeV}\), \(P_{e\mu }=P_{e\tau }=0.2\), \(P_{ee}=0.6\) at around \(L/E = 8.10~\mathrm{km/MeV}\), and \(P_{ee}=P_{e\mu }=0.41\), \(P_{e\tau }=0.18\) at around \(L/E = 10.31~\mathrm{km/MeV}\), respectively, are plotted in Fig. 3. They are all isosceles triangles, in which \(F_{123}\) is the square root of the area of the concurrence triangle, and \(C_{GMC}\) is the corresponding shortest edge. One can see that concurrence fill is always larger than or equal to the GMC, meaning that concurrence fill possesses more information than GMC.

Fig. 3
figure 3

The lengths of the edges and the areas for the concurrence triangles with \(L/E = 4.61~\mathrm{km/MeV}\), \(L/E = 8.10~\mathrm{km/MeV}\), and \(L/E = 10.31~\mathrm{km/MeV}\), respectively, are shown in the three-flavored electron NOs

5 Entanglement in muon antineutrino oscillations

If the muon flavor state is prepared at initial time \(t=0\), the evolution of states for three-flavored NOs can be expressed as

$$\begin{aligned} {\left| {\psi _\mu (t)} \right\rangle } = {a_{\mu e}}(t)\left| {100} \right\rangle + {a_{\mu \mu }}(t)\left| {010} \right\rangle + {a_{\mu \tau }}(t)\left| {001} \right\rangle . \end{aligned}$$
(30)

To quantify the entanglement measure, we shall focus on the the density matrix \(\rho _{ABC}^\mu (t)={\left| {\psi _\mu (t)} \right\rangle }\left\langle \psi _\mu (t)\right| \). This gives

$$\begin{aligned} \rho _{ABC}^\mu (t) = \left( {\begin{array}{*{20}{c}} 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad {\rho _{22}^\mu }&{}\quad {\rho _{23}^\mu }&{}\quad 0&{}\quad {\rho _{25}^\mu }&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad {\rho _{32}^\mu }&{}\quad {\rho _{33}^\mu }&{}\quad 0&{}\quad {\rho _{35}^\mu }&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad {\rho _{52}^\mu }&{}\quad {\rho _{53}^\mu }&{}\quad 0&{}\quad {\rho _{55}^\mu }&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0&{}\quad 0 \end{array}} \right) , \end{aligned}$$
(31)

where the matrix elements are given by

\(\rho _{22}^\mu \!=\! {\left| {{a_{\mu \tau }}(t)} \right| ^2}\); \(\rho _{23}^\mu (t) \!=\! {a_{\mu \tau }}(t)a_{\mu \mu }^ * (t)\); \(\rho _{25}^\mu \!=\!{a_{\mu \tau }}(t)a_{\mu e}^ * (t);\) \(\rho _{32}^\mu \!\!=\!{a_{\mu \mu }}(t)a_{\mu \tau }^ * (t)\); \(\rho _{33}^\mu (t)\!\!= \!\!{\left| {{a_{\mu \mu }}(t)} \right| ^2}\); \(\rho _{35}^\mu (t)\!=\!{a_{\mu \mu }}(t)a_{\mu e}^ * (t);\) \(\rho _{52}^\mu \!=\!{a_{\mu e}}(t)a_{\mu \tau }^*(t)\); \(\rho _{53}^\mu (t)\! =\! {a_{\mu e}}(t)a_{\mu \mu }^*(t)\); \(\rho _{55}^\mu (t)\!=\!{\left| {{a_{\mu e}}(t)}\right| ^2}.\)

Fig. 4
figure 4

Oscillation probabilities as a function of the ratio between the traveled distance L and the energy E. The figure plots the oscillation probability \({v_\mu }(0) \rightarrow {v_\mu }(t)\) (red, solid line), \({v_\mu }(0) \rightarrow {v_e}(t)\) (purple, dashed line), \({v_\mu }(0) \rightarrow {v_\tau }(t)\) (blue, dashed-dotted line) when the initial neutrino flavor is muon flavor

The probabilities for finding the neutrino in state \(|\nu _{e}\rangle \), \(|\nu _{\mu }\rangle \) and \(|\nu _{\tau }\rangle \) are, respectively, \(P_{\mu e}(t)=|a_{\mu e}(t)|^2\), \(P_{\mu \mu }(t)=|a_{\mu \mu }(t)|^2\) and \(P_{\mu \tau }(t)=|a_{\mu \tau }(t)|^2\), for the initial muon flavor neutrino. In Fig. 4, we show the oscillation probabilities of the muon neutrino oscillation as a function of L/E with logarithmic scale. While the transition probability of \(P_{\mu \tau }(t)\) shows a relatively large range of variation, the transition probability of \(P_{\mu e}(t)\) takes a trivial value with respect to L/E with dimension \({\mathrm{km}}/{\mathrm{GeV}}\).

For the three-flavour state \({\left| {\psi _ \mu (t)} \right\rangle }\) , the GGM is

$$\begin{aligned} G(\rho ^ \mu )=1-max\{\lambda ^\mu _A,\, \lambda ^\mu _B,\, \lambda ^\mu _C\}, \end{aligned}$$
(32)

where \(\lambda ^\mu _A\), \(\lambda ^\mu _B\) and \(\lambda ^\mu _C\) are the largest eigenvalues of the reduced density matrix \(\rho ^\mu _A\), \(\rho ^\mu _B\), \(\rho ^\mu _C\), respectively. By calculation, \(\lambda ^\mu _A=max\{P_{\mu \tau }+P_{\mu \mu }, P_{\mu e} \}\), \(\lambda ^\mu _B=max\{P_{\mu e}+P_{\mu \tau }, P_{\mu \mu }\}\), \(\lambda ^ \mu _C=max\{P_{\mu e}+P_{\mu \mu }, P_{\mu \tau } \}\).

Similar to the case of electron neutrino oscillation system, we can calculate that the corresponding three-\(\pi \) entanglement as

$$\begin{aligned}&\pi _{ABC}(\rho ^ \mu )\nonumber \\&\quad =\frac{4}{3}\Bigg [-P^2_{\mu e}-P^2_{\mu \mu }-P^2_{\mu \tau }+P_{\mu e}\sqrt{P^2_{\mu e}+4P_{\mu \mu }P_{\mu \tau }}\nonumber \\&\qquad +P_{\mu \mu }\sqrt{P^2_{\mu \mu }+4P_{\mu e}P_{\mu \tau }}+P_{\mu \tau }\sqrt{P^2_{\mu \tau }+4P_{\mu e}P_{\mu \mu }}\Bigg ]. \end{aligned}$$
(33)

The GMC, as the smallest of three one-to-other bipartite entanglements of the three-flavour muon neutrino system, is expressed as

$$\begin{aligned} C_{GMC}(\rho ^\mu )&=\mathrm {min}\,\big \{4P_{\mu e}\big (P_{\mu \mu }+P_{\mu \tau }\big ),\, 4P_{\mu \mu }\big (P_{\mu e}+P_{\mu \tau }\big ) \nonumber \\&\quad 4P_{\mu \tau }\big (P_{\mu e}+P_{\mu \mu } \big )\big \}. \end{aligned}$$
(34)

Concurrence fill, as a new method of measure the genuine multipartite entanglement, is obtained

$$\begin{aligned}&F_{123}(\rho ^ \mu )\nonumber \\&\quad =8\Bigg \{\frac{(P_{\mu e}P_{\mu \mu }P_{\mu \tau })^2[P_{\mu \tau }P_{\mu \mu }+P_{\mu e}(P_{\mu \mu }+P_{\mu \tau })]}{3}\Bigg \}^{1/4}. \end{aligned}$$
(35)
Fig. 5
figure 5

In the muon antineutrino oscillation system, the multipartite entanglement measured by GGM (black dotted line), three-\(\pi \) (blue, dashed-dotted line), GMC (purple, dashed line), and concurrence fill (red, solid line)

Fig. 6
figure 6

The lengths of the edges and the areas for the concurrence triangles with \(L/E = 262.2~{\mathrm{km}}/{\mathrm{GeV}}\), \(L/E = 479.9~{\mathrm{km}}/{\mathrm{GeV}}\), and \(L/E = 1130~{\mathrm{km}}/{\mathrm{GeV}}\), respectively, are shown in the three-flavored muon NOs

In Fig. 5, we plot GGM, three-\(\pi \), GMC, and concurrence fill in three-flavored muon neutrino oscillations as a function of ratio L/E with dimension \({\mathrm{km}}/{\mathrm{GeV}}\). The minimum of the first concave interval of four entanglement measures are approximatively given by 0.013, 0.018, 0.051, and 0.063, respectively, when at around \(L/E = 513.4~{\mathrm{km}}/{\mathrm{GeV}}\) for three-flavored muon neutrino oscillations. At the point \(L/E=0\), only the survive probability is 1, and the other two transition probabilities are 0, and the four measures of multipartite entanglement are 0. While GMC and GGM have six nonanalytical sharp peaks, the change of concurrence fill are smooth in the range of [10, 1600] of L/E. Also, the concurrence fill is always larger than GGM and three-\(\pi \) in the whole change range of parameter. Again, we can find that concurrence fill is a more genuine tripartite measure of tripartite entanglement possessing the advantage of smoothing change and more information simultaneously, compared to other three tripartite entanglement measures.

In Fig. 6, we plot the concurrence triangles for the three-flavored muon neutrino systems in the case that \(P_{\mu \mu }=P_{\mu \tau }=0.488\), \(P_{\mu e}=0.024\) at around \(L/E = 262.2~{\mathrm{km}}/{\mathrm{GeV}}\) , \(P_{\mu \mu }=0.022\), \(P_{\mu \tau }=0.937\), \(P_{\mu e}=0.041\) at around \(L/E = 479.9~\mathrm{km}/\mathrm{GeV}\), and \(P_{\mu \mu }=0.894, P_{\mu \tau }=0.087\), \(P_{\mu e}=0.019\) at around \(L/E = 1130~{\mathrm{km}}/{\mathrm{GeV}}\), respectively. The \(F_{123}\) and \(C_{GMC}\) are the square root of the area and shortest edge of the concurrence triangle, respectively. For the \(L/E = 262.2~{\mathrm{km}}/{\mathrm{GeV}}\), \(L/E = 479.9~{\mathrm{km}}/{\mathrm{GeV}}\), and \(L/E = 1130~{\mathrm{km}}/{\mathrm{GeV}}\), the square root of the area of the concurrence triangle are 0.33, 0.09, and 0.13, respectively, and the corresponding shortest edges are 0.09, 0.09, and 0.08. The shortest edges is always not more than the square root of the area, suggesting that concurrence fill contains more information than GMC.

6 Summary

In summary, we show that concurrence fill is a genuine tripartite measure to quantify the quantumness of three-flavor electron and muon neutrino oscillations. For initial electron-neutrino oscillations, concurrence fill can reach the maximum 0.89, while the maximum of other three tripartite entanglement measures, including GGM, three-\(\pi \), and GMC are 0.32, 0.55, and 0.88, respectively. It has been proven that concurrence fill is a more natural measure compared to other three measures. The reasons are two-manifold: (i) Concurrence fill shows a smoothing change with respect to L/E, since it considers all the bipartite entanglement measures. GGM and GMC have nonanalytical sharp peaks due to the nonanalytic minimum argument in its expression, which lose sight of the global distribution of entanglement among the parties. (ii) Concurrence fill contains more quantum resource in the electron and muon neutrino propagation. Concurrence fill is always larger than or equal to the other three entanglement measures. All these features make concurrence fill a genuine information-theoretic quantification of tripartite quantum entanglement of quantum resource theories in three-flavor NOs.