\section{Introduction}
One of the unresolved puzzles in neutrino physics is whether neutrinos are Dirac or Majorana fermions.
In the Dirac description, neutrinos are different from their antiparticles and the lepton number (L) is conserved.
Whereas in the Majorana picture, neutrinos and anti-neutrinos are not physically distinguishable. If neutrinos are Majorana fermions, neutrinoless double beta decay ($0\nu 2 \beta$) process ($X^Z_A \rightarrow Y^A_{Z + 2} + 2 e^-$) could occur with the violation of lepton number ($\Delta L = 2$)~\cite{vergados2012theory,bilenky2015neutrinoless}. To this day, there is no experimental evidence for $0\nu 2 \beta$ process.
The most recent combined analysis of KamLAND-Zen 400 (2011-2015) and KamLAND-Zen 800 (started 2019) data predicts the half-life of $T_{0\nu 2 \beta}^{1/2}>2.3 \times 10^{26}$ yr at 90$\%$ C.L. and this corresponds to an effective neutrino mass ($m_{\beta \beta}$) less than 36-156 meV~\cite{zen2022first}. The smallness of neutrino mass has been successfully explained in beyond the standard model (BSM) physics, where neutrinos are assumed to be Majorana particles. Therefore, to resolve this puzzle, holds a strong theoretical motivation for BSM physics.

Experimental data from nearly two decades has established neutrino oscillations as a leading mechanism for neutrino flavour transitions. However, sub-leading effects like non-standard interactions, quantum decoherence, neutrino decay, still hold a torch for new physics scenarios beyond the standard model.
In this precision era of neutrino physics experiments, it is the need of the hour to investigate the implications of these sub-leading effects. In this context, we study the effect of quantum decoherence on the neutrino oscillation probabilities in an open quantum system framework. The open quantum system is modelled by considering the interaction of the neutrino subsystem with the environment~\cite{farzan2008reconciling,oliveira2010quantum,bakhti2015revisiting,guzzo2016quantum,gomes2017parameter,capolupo2019decoherence}. The interactions of this kind could originate from the effects of quantum gravity \cite{Hawking_Unpredictability, Hawking_Virtual_black_holes, Hawking_Wormholes, Hawking_Gravitational_entropy}, strings and branes  \cite{Ellis_String_theory, Ellis_CPT_violation, Benatti_Damped_harmonic} at the Plank scale. Consequently, they manifest as dissipation effects and modify the neutrino oscillation probabilities. Moreover, in ref.~\cite{benatti2001massless,richter2017leggett,capolupo2019decoherence,buoninfante2020revealing,carrasco2022uncovering}, it has been shown that one can probe the nature of neutrinos when the neutrino system interacts with the environment. The authors in ref~\cite{richter2017leggett} have presented that Leggett-Garg $K_3$ quantity takes different values for Majorana and Dirac neutrinos in a dissipative environment.

In the present work, we analyse the effect of a dissipative environment on the neutrino oscillation probabilities in a two flavour picture including the matter effect. The probabilities depend on the Majorana phase and thus provide a window to probe the nature of neutrinos for different baselines. In this context, we study the transition probabilities at different baselines of T2K, NOvA and the upcoming experiments ESS$\nu$SB, T2HKK, DUNE experiments and their dependency on the Majorana phase. In addition, we verify how the measures of quantum correlation like Bell-CHSH inequality, Steering inequality and Non-local advantage of quantum coherence (NAQC) vary with respect to the neutrino beam energy in the case of all the five experiments T2K, NOvA, ESS$\nu$SB, T2HKK, DUNE. We show, that NAQC serves as a stronger and legitimate quantifier among all three to study the quantumness of a system. We further define a new parameter $\Delta$NAQC and analyse how it varies with respect to the Majorana phase for different baselines.

This paper is structured as follows:  In sec. II, we present a basic formalism to determine the neutrino oscillation probabilities in matter while assuming decoherence among neutrino mass eigen states.  We discuss in sec. III, the short description of the neutrino oscillation experiments that have been used in this study. In sec. IV, we show the oscillation probabilities relevant to these experiments and discuss their implications by considering two cases of decoherence matrix. In sec. V, we present a brief account of several non-classical measures of quantum coherence and further study probe the nature of neutrinoa using a new parameter $\Delta NAQC$. We finally summarise our findings in sec. VI.

\section{Formalism of neutrino oscillations in matter assuming decoherence}
In an open quantum system framework, the neutrino subsystem interacts weakly with the environment leading to a loss of coherence in the subsystem.
The decoherence phenomenon in Markovian systems is described by the Lindblad-Kossakowski master equation~\cite{gorini1978properties,manzano2020short}
\begin{align}\label{lindblad}
    \mathcal{L} \rho(t) = \frac{\partial\rho}{\partial t} = -i[H_{eff} , \rho(t)] + \mathcal{D}[\rho(t)]~.
 \end{align}
 Here, the infinitesimal generator $\mathcal{L}$ and its action on the density matrix $\rho(t)$ depend on effective Hamiltonian $H_{eff}$ and the dissipative term $\mathcal{D}$.
The dissipative factor $\mathcal{D}$ has the following form,
\begin{align}\label{dissipative}
    \mathcal{D}[\rho] = - \frac{1}{2}\sum_{j=0}^{N^2 - 1}(A_j^\dagger A_j \rho + \rho A_j^\dagger A_j) + \sum_{j=0}^{N^2 - 1} A_j \rho A_j^\dagger~.
\end{align}
Here $A_j$ are the operators that depend on the system dimensions. In an N-level system, $A_j$ has $N \times N$ dimension and they form $N^2 - 1$ linearly independent basis without including identity. For the two flavor mixing, $A_j$ are the Pauli matrices $\sigma_i$ and that for the three flavor are represented by Gell-Mann matrices $\lambda_i$. In the former case the density matrix $\rho$ and operator $A_j$ in eq.~(\ref{dissipative}), can be written as, $\rho = \frac{1}{2}\rho_\mu \sigma_\mu$ , $A_j = a_\mu^j\sigma_\mu$ where, $\mu = 0, 1, 2, 3$, $\sigma_0$ is $2 \times 2$ identity matrix and $\sigma_{i =1,2,3}$ are the Pauli matrices.

In this work, we consider ultrarelativistic electron neutrinos ($\nu_e$) and muon neutrinos ($\nu_{\mu}$) in two dimensional Hilbert space. In two generation model of neutrinos, the flavor states ($\nu_e$, $\nu_{\mu}$) are related to the mass states ($\nu_1$, $\nu_2$) by a $2 \times 2$ unitary mixing matrix U
\begin{align}\label{mixing}
    \begin{pmatrix}
    \nu_e\\
    \nu_\mu
    \end{pmatrix}
    = U
    \begin{pmatrix}
    \nu _1\\
    \nu_2
    \end{pmatrix},
    \end{align}
    where
    \begin{align}\label{mixing_matrix}
       U = \begin{pmatrix}
    \cos\theta  & \sin\theta \hspace{0.2cm} e^{-i \phi}\\
    - \sin\theta\hspace{0.2cm} e^{i \phi} & \cos\theta
    \end{pmatrix}~,
    \end{align}

$\theta$ is the mixing angle and the phase $\phi$ represents the Majorana phase.

 In eq.~(\ref{lindblad}), one can note that the evolution of $\rho$ depends on a time independent effective Hamiltonian  $H_{eff}$, which is combination of vacuum Hamiltonian ($H_{vac}$) and matter interaction Hamiltonian ($H_{mat}$)
\begin{align}
    H_{eff} = H_{vac} + H_{mat} = h_\mu \sigma_\mu~ ,
\end{align}
here $h_\mu$ are the components of $H_{eff}$.
 If neutrinos propagate in vacuum, the free Hamiltonian will be in the following form,
\begin{align}
   H_{vac} =  \begin{pmatrix}
   - \omega & 0\\
   0 &  \omega
   \end{pmatrix}
\end{align}\\
where $\omega = \frac{\Delta m^2}{4E}$ contains the square mass difference of two mass eigenstates ($\Delta m^2$) and $E$ represents the neutrino energy. The interaction of the neutrinos with the electrons in the medium is given by the interaction Hamiltonian,
   \begin{align}
        H_{mat} = U^\dagger \begin{pmatrix}
        A & 0\\
        0 & 0
        \end{pmatrix} U
   \end{align}
   \begin{align}\label{Hmat}
     = A
    \begin{pmatrix}
    \cos^2\theta  & \cos\theta \sin\theta \hspace{0.1cm} e^{-i \phi}\\
     \cos\theta \sin\theta \hspace{0.1cm} e^{i \phi} & \sin^2\theta
    \end{pmatrix}.
\end{align}

where $A = \sqrt{2}G_Fn_e$. Here $G_F$ is called the Fermi constant and the electron number density in the medium is denoted by $n_e$.

The dissipative matrix $\mathcal{D}_{mn}$ (based on positivity and trace-preserving conditions) depends on six real and independent parameters $\alpha,~\beta,~\gamma,~a,~b,~c$. $\mathcal{D}_{mn}$ is given by,

   \begin{align} \label{Dmn}
    \mathcal{D}_{mn} = -2\begin{pmatrix}
    0&0&0&0\\
    0&\alpha&\beta&\gamma\\
    0&\beta&a&b\\
    0&\gamma&b&c
    \end{pmatrix}~,
\end{align}
where $m,n=0,1,2,3$.

Now, the time evolution of the density matrix $\rho$ in eq.~(\ref{lindblad}) can be expressed through Schr$\Ddot{o}$dinger like equation
\begin{eqnarray}\label{sch}
\frac{d}{dt}\ket{\rho(t)} = -2 \mathcal{H}\ket{\rho(t)}~,
\end{eqnarray}
with
    \begin{align}
    \mathcal{H} = 2\begin{pmatrix}
    0 & 0 & 0 & 0\\
    0 & \alpha & \beta + \xi & \gamma - \lambda \sin\phi\\
    0 & \beta - \xi & a & b + \lambda \cos\phi\\
    0 & \gamma + \lambda \sin\phi & b - \lambda \cos\phi & c
    \end{pmatrix}
\end{align}

and \begin{align}
    \xi = \frac{A}{2}\cos2\theta - \omega, \hspace{0.5cm} \lambda = \frac{A}{2}\sin2\theta
\end{align}

After imposing trace preserving condition $\Dot{\rho_0}(t) = 0$, eq.~(\ref{sch}) will take the following form
\begin{align}
    \begin{pmatrix}
     \Dot{\rho_1}(t) \\ \Dot{\rho_2}(t) \\ \Dot{\rho_3}(t)
    \end{pmatrix} = -2 \mathcal{H} \begin{pmatrix}
     {\rho_1}(t) \\ {\rho_2}(t) \\ {\rho_3}(t)
    \end{pmatrix}.
\end{align}

 Hence the evolved state at time t is written as

\begin{align}\label{evolveddensity}
    \begin{pmatrix}
     \rho_1(t) \\ \rho_2(t) \\ \rho_3(t)
    \end{pmatrix} = \mathcal{M}(t) \begin{pmatrix}
     {\rho_1(0)} \\ {\rho_2(0)} \\ {\rho_3(0)}
    \end{pmatrix}
\end{align}
where $\mathcal{M}(t) = e^{-2\mathcal{H}t} = \mathcal{S}e^{-2\mathcal{H'}t}\mathcal{S}^{-1}$ is a $3 \times 3$ matrix.
Here, $\mathcal{S}$ is the similarity transformation matrix and $\mathcal{H'}$ is the diagonal of the eigenvalues of $\mathcal{H}$.

The density matrix at arbitrary time t is
\begin{align}\label{density}
 \rho(t) = \frac{1}{2}
    \begin{pmatrix}
    \rho_0(t) + \rho_3(t) & \rho_1(t) - i \rho_2(t)\\
    \rho_1(t) + i \rho_2(t) & \rho_0(t) - \rho_3(t)
    \end{pmatrix}~.
\end{align}

Using eq.~(\ref{mixing_matrix}), the density matrix at initial time t = 0 can be obtained as
\begin{align}
    \rho_{\nu_e}(0) = \begin{pmatrix}
    \cos^2\theta & \frac{1}{2}\sin2\theta \hspace{0.2cm} e^{i\phi}\\
    \frac{1}{2}\sin2\theta \hspace{0.2cm} e^{-i\phi}
& \sin^2\theta
\end{pmatrix}
\end{align}
and
\begin{align}\label{muondensity}
    \rho_{\nu_\mu}(0) = \begin{pmatrix}
    \sin^2\theta & - \frac{1}{2}\sin2\theta \hspace{0.2cm} e^{i\phi}\\
    - \frac{1}{2}\sin2\theta \hspace{0.2cm} e^{-i\phi}
& \cos^2\theta
\end{pmatrix}~.
\end{align}

Further, the evolved state $\rho_{\nu_\mu}(t)$ can be obtained using eq.~(\ref{evolveddensity})-(\ref{muondensity}). The probability of transition of an initial state $\nu_\alpha$ to a final state $\nu_\beta$ can be evaluated using
\begin{align}
    P_{\nu_\alpha \rightarrow \nu_\beta}(t) = Tr[\rho_{\nu_\alpha}(t) \rho_{\nu_\beta}(0)]~.
\end{align}
Upon substitution, the appearance $ P_{\nu_\mu \rightarrow \nu_e}(t)$ and disappearance $P_{\nu_\mu \rightarrow \nu_\mu}(t)$ probabilities are obtained to be
 \begin{equation}\label{nue-app prob}
 \begin{split}
     P_{\nu_\mu \rightarrow \nu_e}(t) = \frac{1}{2} \Bigl[1 - M_{33} \cos^2(2\theta) + \sin(2 \theta) \cos(2\theta) \bigl\{(M_{23} + M_{32}) \sin\phi\\
       - (M_{13} + M_{31})\cos\phi \bigr\} +\sin^2(2\theta) \bigl\{(M_{12} + M_{21})\sin\phi \cos\phi\\ -
       (M_{11}\cos^2\phi + M_{22}\sin^2\phi)\bigr\}\Bigr]~,
 \end{split}
 \end{equation}

  \begin{equation}\label{numu-dis prob}
 \begin{split}
     P_{\nu_\mu \rightarrow \nu_\mu}(t) = \frac{1}{2} \Bigl[1 + M_{33} \cos^2(2\theta) - \sin(2 \theta) \cos(2\theta) \bigl\{(M_{23} + M_{32}) \sin\phi\\
       - (M_{13} + M_{31})\cos\phi \bigr\} - \sin^2(2\theta) \bigl\{(M_{12} + M_{21})\sin\phi \cos\phi\\ -
       (M_{11}\cos^2\phi + M_{22}\sin^2\phi)\bigr\}\Bigr] = 1 -  P_{\nu_\mu \rightarrow \nu_e}(t)~,
 \end{split}
 \end{equation}

where $M_{ij}$ are the matrix elements of $\mathcal{M}$ and $i,j = 1,2,3$.

The probability of antineutrinos under the same conditions can be obtained from the probability of neutrinos by replacing the $A \rightarrow -A$ and $U \rightarrow U^*$ in eq.~(\ref{Hmat}). Then $P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_e}(t)$ and $P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_\mu}(t)$ are deduced as

\begin{equation}\label{antinue app prob}
 \begin{split}
     P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_e}(t) = \frac{1}{2} \Bigl[1 - M_{33} \cos^2(2\theta) - \sin(2 \theta) \cos(2\theta) \bigl\{(M_{13} + M_{31})\cos\phi \\
        + (M_{23} + M_{32}) \sin\phi\bigr\} - \sin^2(2\theta) \bigl\{(M_{12} + M_{21})\sin\phi \cos\phi\\ +
       (M_{11}\cos^2\phi + M_{22}\sin^2\phi)\bigr\}\Bigr]~,
 \end{split}
 \end{equation}

 and

 \begin{equation}\label{antinumu dis prob}
 \begin{split}
     P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_\mu}(t)  = 1 -  P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_e}(t)
 \end{split}~.
 \end{equation}

The eqns.(\ref{nue-app prob} - \ref{antinumu dis prob}) show that in the presence of decoherence the neutrino and anti-neutrino oscillation probabilities  depend on the Majorana phase $\phi$. This modified oscillation probabilities open a window to explore the nature of neutrinos in the current and the upcoming neutrino oscillation experiments. In this context, we analyse the probabilities of various long baseline neutrino oscillation experiments and their dependency on the Majorana phase $\phi$ in the presence of decoherence. In the following section, we give a brief account of the five experiments considered in this study.

\section{Experimental details}
\begin{table}[b]
      \centering
      \begin{tabular}{l r r}
        \toprule[1.5pt]
       \textbf{Experiment}~&~\textbf{L(km)}~&~\textbf{E(GeV)} \\
        \midrule[0.5pt]
T2K & 295 & 0.6\\
   ESS & 540 & 0.2 \\
   NOvA & 810 & 2.0 \\
   T2HKK & 1100 & 0.7 \\
   DUNE & 1300 & 2.8  \\
        \bottomrule[1.5pt]
      \end{tabular}
      \caption{The baseline length and the peak neutrino beam energy for different experiments are considered in this work.}
      \label{table:exp-details}
    \end{table}

    \begin{table}[b]
      \centering
      \begin{tabular}{l r}
        \toprule[1.5pt]
        \textbf{Parameter} & \textbf{Value}  \\
        \midrule[0.5pt]
        $\Delta m^2_{32}$ & $\pm$ 2.54 $\times 10^{-3}$ $eV^2$  \\
        $\theta_{23}$ & 49.2$^{\circ}$  \\
        $\alpha$ = $a$ = $c$ & 7.8 $\times 10^{-24}$ GeV \\
        $\beta$ & 3.8 $\times 10^{-24}$ GeV  \\
        $b$ & 3.0 $\times 10^{-24}$ GeV \\
        \bottomrule[1.5pt]
      \end{tabular}
      \caption{Standard neutrino oscillation parameters and decoherence parameters. Plus and minus sign in $\Delta m^2_{32}$ refers to the normal and inverted ordering respectively.}
      \label{table:osc-deco-params}
    \end{table}

T2K and T2HKK : Tokai to Kamiokande (T2K) experiment~\cite{abe2013t2k} is an off-axis (off-axis angle (OAA) of $2.5^{\degree}$) oscillation experiment with Japan Proton Accelerator Research Complex (J-PARC) based $\nu_{\mu}$ beam facility. The far detector of volume 22.5 kt is placed at Kamiokande with a baseline of 295 km. The motivation of this experiment is to precisely measure oscillation parameters, $\theta_{13}$, $\theta_{23}$ and $\Delta m^2_{32}$.
The T2HKK experiment~\cite{hyper2018physics} is proposed to have an off-axis $\nu_{\mu}$ beam (OAA ranging from 1-3$\degree$) from J-PARC facility, to travel a distance of 1100 km before it reaches the Water Cherenkov detector of 187 kt based in Korea.

NOvA : The NuMI Off-axis $\nu_e$ Appearance Experiment~\cite{itow2001jhf} is an ongoing neutrino oscillation experiment with a baseline of 810 km. It has an off-axis (OAA 0.8$\degree$) muon neutrino beam with a peak energy of $\sim$2 GeV. NOvA has a near detector at the Fermilab site and NuMI beam focussed towards a far detector of volume 14 kt, placed at Minnesota. The main goal is to understand the atmospheric neutrino flavor transition and also measure atmospheric mass square difference to the higher precision level ($10^{-4}$ $eV^2$).

ESS$\nu$SB : The future neutrino experiment European Spallation Source Neutrino Super Beam (ESS$\nu$SB) experiment~\cite{baussan2012use} is a long baseline experiment (baseline 540 km). The peak energy of the neutrino beam is around 0.2 GeV which is the energy corresponding to the second oscillation maxima. A megaton Water Cherenkov far detector is placed underground in ESS site at Lund. This experiment is sensitive to observe leptonic CP-violation phase at 5$\sigma$ confidence level.

DUNE : The Deep Underground Neutrino Experiment (DUNE)~\cite{acciarri2016long} is an on-axis accelerator based long baseline neutrino experiment. The DUNE mainly consists of a beamline, a near detector at Fermilab and a far detector at Sanford Underground Research Facility (SURF) which is 1300 km away from the near detector. It will measure the neutrino events having broad range of energy (1-8 GeV). The neutrino flux is peaked around 2.8 GeV corresponding to the energy at the first oscillation maxima. This experiment will help to determine charge-parity (CP) violation phase and the neutrino mass ordering with very high precision.

The baselines and the corresponding peak neutrino beam $\nu_{\mu}$ energy values of the above experiments are tabulated in table~\ref{table:exp-details}. The details of the oscillation parameters and the decoherence parameters considered in this work are listed in table~\ref{table:osc-deco-params}.

\section{Numerical Analysis}
 We consider the two-flavour neutrino oscillation analysis of $\nu_{\mu}$-$\nu_e$ oscillations in a dissipative medium presented in section-II.
 We provide the probability plots corresponding to the eqns.~(\ref{nue-app prob} - \ref{antinumu dis prob}) as these are the major oscillation channels to be studied in the experiments T2K, NOvA, ESS$\nu$SB, T2HKK and DUNE. We further substitute $t=L$ in natural units, where L is the distance travelled by the neutrino beam.
We assume a constant matter density profile with $\rho = 2.8~ g/cm^3$ corresponding to the matter density potential $A = 1.01 \times 10^{-13}~eV$. In this work, we analyse two cases of decoherence matrix. Firstly, for simplicity,  we assume all the off-diagonal elements of the decoherence matrix in eq.~(\ref{Dmn}) to be zero and secondly we consider non-zero diagonal and two non-zero off-diagonal elements. We assume the neutrino mass ordering as normal ordering throughout the work unless otherwise mentioned.

\subsection{\textit{Non-zero diagonal elements in $D_{mn}$:}}

Firstly, we assume that the decoherence matrix $D_{mn}$ in eq.~{\ref{Dmn}} has only non-zero diagonal elements and impose that the off-diagonal elements are zero. Under this assumption $D_{mn}$ in eq.~{\ref{Dmn}} takes a simple form
   \begin{align} \label{diag_Dmn}
    \mathcal{D}_{mn} = -2\begin{pmatrix}
    0&0&0&0\\
    0&\alpha&0&0\\
    0&0&a&0\\
    0&0&0&c
    \end{pmatrix}~,
\end{align}
where $\alpha = a = c$. Now, using the density matrix formalism presented in sec.II, we numerically obtain the oscillation probabilities $P_{\nu_\alpha \rightarrow \nu_\beta}$ and $P_{\bar{\nu}_\alpha \rightarrow \bar{\nu}_\beta}$ where $(\alpha,\beta)=(e,\mu)$.

\begin{figure}[!htbp]
\begin{center}
\includegraphics[height=10cm,width=\textwidth]{PvsE_non-zero_diag_0.1-5GeV.png}
\end{center}
\vspace{-0.8cm}
\caption{\footnotesize Appearance and disappearance probabilities of neutrinos (upper panel) and anti-neutrinos (lower panel) with respect to energy E (assuming normal mass ordering). The orange, green, red, cyan and black curves correspond to the probability versus energy E for the baselines 295 km, 540 km, 810 km, 1100 km and 1300 km respectively.}
\label{fig:diag-Prob-NH-phi-90}
\end{figure}

\begin{figure}[!htbp]
\begin{center}
\includegraphics[height=6cm,width=\textwidth]{DeltaPvsphi_non-zero_diag_Pme_Pbar_me.png}
\end{center}
\vspace{-0.8cm}
\caption{\footnotesize $\Delta P_{\alpha \beta}$ versus $\phi$ in the left panel and $\Delta P_{\bar{\alpha} \bar{\beta}}$ versus $\phi$ in the right panel with $(\alpha, \beta) = (e, \mu)$.}
\label{fig:diag-del-Prob-NH-phi-90}
\end{figure}

In fig~\ref{fig:diag-Prob-NH-phi-90}, we plot the transition probabilities with respect to the neutrino beam energy E. We plot the appearance probability $P_{\mu e}=P_{\nu_\mu \rightarrow \nu_e}$ and the disappearance probability $P_{\mu \mu} = P_{\nu_\mu \rightarrow \nu_\mu}$ of the $\nu_{\mu}$ beam with respect to E, in the top left and right panels respectively. Additionally, in the bottom panel we show $P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_e}$ and $P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_\mu}$. The solid lines correspond to $\phi \neq 0$ which represents the non-zero Majorana phase while dashed lines correspond to $\phi=0$ which is the case with Dirac neutrinos. We consider maximal phase  $\phi = \pi/2$ as a representative value in these plots.
 The orange, green, red, cyan and black curves correspond to the probability versus energy E for the baselines 295 km, 540 km, 810 km, 1100 km and 1300 km respectively. As can be seen from the orange curves, in the case of T2K experiment where baseline length is 295 km the dashed curve and the solid curve are almost overlapping with each other, indicating minimal effect of decoherence at this L and E. As the baseline increases, we can see from the green (540 km), red (810 km), cyan (1100 km) and black (1300km) curves that the impact of decoherence increases and reaches maximum for DUNE experiment.
 However, the ESS$\nu$SB experiment (540 km) and T2HKK (1100 km) are designed to work at the second oscillation maxima i.e. 0.2 GeV and 0.7 GeV respectively. At these energy values one can see negligible difference between the Dirac (solid curve) and Majorana (dashed curve) cases from the green and cyan curves.
 While the solid and dashed black curves corresponding to DUNE baseline 1300 km show maximum sensitivity to discriminate between Dirac and Majorana phase.

 We further study the effect of decoherence on the determination of neutrino nature at the five experiments considered, by defining two quantities $\Delta P_{\alpha \beta}$ and $\Delta {P_{\bar{\alpha} \bar{\beta}}}$.
The two quantities $\Delta P_{\alpha \beta}$ and $\Delta {P_{\bar{\alpha} \bar{\beta}}}$ in the case of neutrinos and anti-neutrinos are defined as below
   \begin{align}\label{DelP-nu}
        \Delta P_{\alpha \beta} = |P_{\nu_\alpha \rightarrow \nu_\beta} (Dirac) - P_{\nu_\alpha \rightarrow \nu_\beta} (Majorana)|,
    \end{align}
   \begin{align}\label{DelP-nubar}
        \Delta  {P_{\bar{\alpha} \bar{\beta}}} = |P_{\bar{\nu}_\alpha \rightarrow \bar{\nu}_\beta} (Dirac) - P_{\bar{\nu}_\alpha \rightarrow \bar{\nu}_\beta} (Majorana)|.
    \end{align}
The first term in eq.~(\ref{DelP-nu}) and eq.~(\ref{DelP-nubar}) is obtained by assuming $\phi = 0$ while the second is obtained by taking $\phi \subset [0,2\pi]$. When $\alpha=\beta=\mu$ in eq.~(\ref{DelP-nu}), one can derive that $\Delta P_{\mu \mu}=\Delta P_{\mu e}$ since our analysis is done for two flavour neutrino oscillation picture. Similarly, from eq.~(\ref{DelP-nubar}), we obtain $\Delta P_{\bar{\mu} \bar{\mu}}=\Delta P_{\bar{\mu} \bar{e}}$.
Therefore, we numerically obtain two quantities $\Delta P_{\mu e}$ and $\Delta P_{\bar{\mu} \bar{e}}$ with respect to $\phi$ from eq.~({\ref{DelP-nu}) and eq.~(\ref{DelP-nubar}).

Fig.~\ref{fig:diag-del-Prob-NH-phi-90}, shows $\Delta P_{\mu e}$ vs $\phi$ in the left panel and $\Delta P_{\bar{\mu} \bar{e}}$ vs $\phi$
in the right panel. To obtain this figure, we fix the peak neutrino beam energy values as per the experimental specifications in table-\ref{table:exp-details} and vary the Majorana phase $\phi \subset [0,2\pi]$. For all non-zero values of $\phi$, the green curve corresponding to the 540 km baseline (ESS$\nu$SB) and the black curve corresponding to the 1300 km baseline (DUNE) show minimum and maximum potential to determine the nature of neutrinos. The green (540 km) and cyan (1100 km) curves show lesser sensitivity to discriminate between Dirac and Majorana neutrinos for all values of $\phi \neq 0$. This can be explained from fig~\ref{fig:diag-Prob-NH-phi-90}, for $\phi = \pi/2$ we have the green, cyan solid (Majorana case) and the dashed (Dirac) curves almost overlapping at the second oscillation maximum. Since, ESS$nu$SB and T2HKK are planned at second oscillation maximum, the effect of decoherence is lesser than that at the first oscillation maxima. For this reason NOvA (810 km) will be more sensitive to determining the nature of neutrinos, next to DUNE as can be seen from the red curve.

 \subsection{\textit{Non-zero off diagonal elements in $D_{mn}$}}
 In this case, we assume that the decoherence matrix $D_{mn}$ in eq.~{\ref{Dmn}} has two non-zero off-diagonal elements and the diagonal elements are equal i.e. $\alpha = a = c$, $\gamma=0$, $\beta$ and $b$ non-zero. Similar to the previous case, we further obtain the oscillation probabilities under this approximation.

 In fig~\ref{fig:Prob-NH-phi-90}, we plot the transition probabilities with respect to the beam energy E. We show the appearance probability $P_{\nu_\mu \rightarrow \nu_e}$ and the disappearance probability $P_{\nu_\mu \rightarrow \nu_\mu}$ of the $\nu_{\mu}$ beam with respect to E, in the top left and right panels respectively. In the bottom left(right) panel, we present $P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_e}$ ($P_{\bar{\nu}_\mu \rightarrow \bar{\nu}_\mu}$) with respect to anti-neutrino $\bar{\nu_{\mu}}$ beam energy. We assume the neutrino mass ordering as normal ordering. The solid lines correspond to $\phi \neq 0$ which represents the case of non-zero Majorana phase while the dashed lines correspond to $\phi=0$, the Dirac neutrinos. We consider maximal phase  $\phi = \pi/2$ as a representative value in these plots. The orange, green, red, cyan and black curves correspond to the probability versus energy E for the baselines 295 km, 540 km, 810 km, 1100 km and 1300 km respectively. In comparison between fig~\ref{fig:diag-Prob-NH-phi-90} and fig~\ref{fig:Prob-NH-phi-90}, we can see that the additional non-zero off-diagonal elements do not change the probabilities vs E significantly. Accordingly, the conclusions drawn for fig~\ref{fig:diag-Prob-NH-phi-90}, also apply in this case. That is there is maximum difference between the Dirac and the Majorana case for 1300 km baseline when one assumes decoherence in a two flavour neutrino framework. The green and cyan curves corresponding to ESS and T2HKK baselines, show that at second oscillation maxima 0.2 GeV and 0.7 GeV the difference between probabilities for Majorana (solid curve) and Dirac (dashed curve) case is almost negligible.

\begin{figure}[!htbp]
\begin{center}
\includegraphics[height=10cm,width=\textwidth]{PvsE_non-zero_b_bita_0.1-5GeV.png}
 \end{center}
\vspace{-0.8cm}
\caption{\footnotesize Appearance and disappearance probabilities of neutrinos (upper panel) and anti-neutrinos (lower panel) with respect to energy E (assuming normal mass ordering). The orange, green, red, cyan and black curves correspond to the probability versus energy E for the baselines 295 km, 540 km, 810 km, 1100 km and 1300 km respectively.}
\label{fig:Prob-NH-phi-90}
\end{figure}

\begin{figure}[!htbp]
\begin{center}
\includegraphics[height=6cm,width=\textwidth]{DeltaPvsphi_non-zero_b_bita_Pme_Pbar_me.png}
\end{center}
\vspace{-0.8cm}
\caption{\footnotesize $\Delta P_{\alpha \beta}$ versus $\phi$ in the left panel and $\Delta P_{\bar{\alpha} \bar{\beta}}$ versus $\phi$ in the right panel with $(\alpha, \beta) = (e, \mu)$.}
\label{fig:del-Prob-NH-phi-b-and-beta-non-90}
\end{figure}
We further obtain $\Delta P_{\mu e}$ and $\Delta {P_{\bar{\mu} \bar{e}}}$ using eq.~(\ref{DelP-nu}) and eq.~(\ref{DelP-nubar}) respectively, for this case of non-zero off-diagonal elements. In fig~\ref{fig:del-Prob-NH-phi-b-and-beta-non-90}, we plot $\Delta P_{\mu e}$ vs $\phi$ in the left panel and $\Delta {P_{\bar{\mu} \bar{e}}}$ vs $\phi$ in the right panel. We fix the L and E values as listed in table-\ref{table:exp-details}. From the black curve, one can see that for all non-conserving values of $\phi \subset [0,2\pi]$, DUNE shows maximum $\Delta P_{\mu e}$ and $\Delta {P_{\bar{\mu} \bar{e}}}$ values and is followed by NOvA experiment. T2K and T2HKK have same $\Delta P_{\mu e}$ and $\Delta {P_{\bar{\mu} \bar{e}}}$ irrespective of their baselines because T2HKK is planned at second oscillation maxima. From the left plot we note that in the case of neutrino beam, the sensitivity to the Majorana phase is seen for all non-zero values. However for the anti-neutrino beam, the $\Delta {P_{\bar{\mu} \bar{e}}}$ values are non-negligible for smaller range of $\phi$ values ([$0.6\pi$ to $2.6\pi$], [$4\pi$ to $5.6\pi$]).

\section{Quantum quantifiers}
To study a quantum system, it is essential to explore the quantumness of the system. It has been shown that various measures of quantum coherence such as, Bell nolocality~\cite{clauser1969proposed,horodecki1995violating,brunner2014bell}, Leggett-Garg inequality (LG-I)~\cite{leggett1985quantum,emary2013leggett}, entanglement~\cite{yin2013lower,blasone2009entanglement}, steering inequality~\cite{wiseman2007steering,costa2016quantification,uola2020quantum}, non-local advantage of quantum coherence (NAQC)~\cite{mondal2017nonlocal} can be expressed in terms of neutrino oscillation probabilities~\cite{alok2016quantum,ming2020quantification} and are useful to shed some light on the unknown neutrino oscillation parameters. Recently, some progress has been made where non-classical correlations are defined to analyze charge-parity (CP) violating phase and the neutrino mass hierarchy problem~\cite{naikoo2020quantum}.

In this section, we briefly discuss the various well established measures of quantum correlations and their relations with neutrino oscillation probability. The two flavour neutrino oscillation framework, considered in this work is ideal to study the quantum quantifiers like Bell-CHSH inequality, Steering inequality and non-local advantage of quantum coherence (NAQC), as they are defined for a two-qubit system. We show that these nonclassical entities take different forms for Majorana and Dirac neutrinos. Therefore, one can consider them as potential candidates to probe the nature of neutrinos.

Firstly, we give a brief account of the measures of quantum correlations Bell-CHSH inequality, Steering inequality and NAQC. Later, we numerically obtain these quantities in terms of neutrino oscillation parameters and plot them with respect to the neutrino beam energy E. For analysis purposes we show the results for three baselines, smallest - 295 km, intermediate - 540 km and largest - 1300 km baselines corresponding to T2K, ESS$\nu$SB, DUNE.

\begin{figure}
\begin{adjustwidth}{-0.2cm}{}
\setlength{\tabcolsep}{5pt}

\begin{tabular}{c c c}
\includegraphics[width =0.32\textwidth]{BellCHSH-E-295.png} & \includegraphics[width =0.32\textwidth]{BellCHSH-E-540.png} & \includegraphics[width =0.32\textwidth]{BellCHSH-E-1300.png}\\[0.5cm]

\includegraphics[width =0.32\textwidth]{Steering-E-295.png} & \includegraphics[width =0.32\textwidth]{Steering-E-540.png} & \includegraphics[width =0.32\textwidth]{Steering-E-1300.png}\\[0.5cm]

\includegraphics[width =0.32\textwidth]{NAQC-E-295.png} & \includegraphics[width =0.32\textwidth]{NAQC-E-540.png} & \includegraphics[width =0.32\textwidth]{NAQC-E-1300.png}\\[0.5cm]
 \end{tabular}
\caption{\footnotesize{Bell-CHSH inequality versus energy (E) in first row, steering inequality versus energy (E) in second row and NAQC versus energy (E) in third row for T2K (left), ESS$\nu$SB (middle), DUNE (right). Solid line, dashed line and dot-dashed line represents decoherence including matter effect, only matter effect (without decoherence) and vacuum respectively}.}
\label{fig:Bell-steering-NAQC-vs-E}
\end{adjustwidth}
\end{figure}

\textit{Bell-CHSH Inequality :}
The 2022 Nobel Prize in Physics has been awarded to Alain Aspect, John F. Clauser and Anton Zeilinger for establishing the violation of Bell's inequalities in experiments with entangled photons. In 1963, it was shown by John Bell that, quantum mechanics is not compatible with local hidden variable theory~\cite{bell1964einstein}. The generalisation of Bell’s theorem presented by Clauser, Horne, Shimony and Holt, to apply on realistic experiment is popularly known as Bell-CHSH inequality~\cite{clauser1969proposed}.
This inequality is used to test the quantum correlation between experimentally measurable quantities which are operated on two spatially separated systems.
The CHSH inequality in terms of the Bell operator $B_{CHSH}$ is given by
\begin{eqnarray} \label{CHSH1}
  |\langle B_{CHSH}\rangle| \leq 2.
  \end{eqnarray}

The violation of the inequality in eq.~({\ref{CHSH1}}), can be rewritten in terms of a quantity $M(\rho)=\max_{i\neq j}(u_i + u_j)$ where $u_{i = 1,2,3}$ are the eigenvalues of a real matrix $T^\dagger T$ (T - correlation matrix) ~\cite{horodecki1995violating}, as
\begin{eqnarray} \label{CHSH2}
|\langle B_{CHSH}\rangle| = 2\sqrt{M(\rho)} > 2
\end{eqnarray}
In the context of two flavour neutrino oscillations, A and B can be considered as neutrinos with two different flavours. The authors of ~\cite{alok2016quantum} have derived $M(\rho)$ in terms of the neutrino oscillation probability $P_{\rm osc}$ and the survival probability $P_{\rm sur}$ as

\begin{eqnarray}
M(\rho) &=& 1 + 4 P_{\rm sur} P_{\rm osc}~.
\end{eqnarray}
Therefore, from eq.~(\ref{CHSH2})
\begin{align}\label{bell}
|\langle B_{CHSH}\rangle| = 2 \sqrt{1 + 4 P_{osc} P_{sur}}~ > 2.
  \end{align}

\textit{Steering Inequality :}
The notion of quantum steering was first proposed by Schrodinger in 1935~\cite{schrodinger1935discussion} to investigate the Einstein-Podolsky-Rosen (EPR) paradox~\cite{einstein1935can}. In a system with two parties, quantum steering is a phenomenon where one party influences/steers the outcome of a measurement on another party. The steering criteria of a bipartite state $\Gamma_{AB}$ shared between two parties (Alice and Bob) can be diagnosed by an inequality that was developed in a seminal paper by Cavalcanti-Jones-Wiseman-Reid~\cite{cavalcanti2009experimental}. If two parties, both are permitted to measure $n$ observables, then the inequality is formulated as~\cite{costa2016quantification}
\begin{eqnarray}
S_n(\Gamma_{AB},\varsigma) = \frac{1}{\sqrt{n}}\left|\sum^n_{i = 1}\langle A_i \otimes B_i\rangle\right| \leq 1
\end{eqnarray}
with $\langle A_i \otimes B_i\rangle$ = Tr($\Gamma_{AB} A_i \otimes B_i$),
where $A_i = {\hat{x}_i\cdot\hat{\sigma}}$, $B_i = {\hat{y}_i\cdot\hat{\sigma}}$, $\hat{x}_i \in \mathbb{R}^3$ are unit vectors, $\hat{y}_i \in \mathbb{R}^3$ are orthonormal vectors, $\varsigma = \{\hat{x}_i,...,\hat{x}_n,\hat{y}_i,...,\hat{y}_n\}$ represents set of measurement directions.

In the case of two flavor neutrino oscillation the steering inequality can be represented as~\cite{ming2020quantification}
\begin{eqnarray}\label{steering}
    S_n(\Gamma_{AB},\varsigma) = \sqrt{1 + 8 P_{osc} P_{sur}}~ \leq 1.
\end{eqnarray}

\textit{Non-local Advantage of Quantum Coherence (NAQC)} :
Recently, authors in ref.~\cite{mondal2017nonlocal} have proposed a quantifier called non-local advantage of quantum coherence (NAQC). This NAQC parameter is based on $l_1$-norm of coherence measure and it represents the quantumness of a system in terms of entanglement and steerability.
In \cite{ming2020quantification} the authors have shown that NAQC acts as a stronger quantifier when compared to Bell-CHSH and steering inequality.

The NAQC-parameter $N_{l_1}$ measures the quantumness of a system  using steerability and quantum coherence criteria. For a two flavor neutrino system, this parameter can be derived in terms of oscillation ($P_{osc}$) and survival probability ($P_{sur}$) as~\cite{ming2020quantification}
\begin{align}\label{NAQC}
   N_{l_1} = 2 + 2\sqrt{P_{osc}P_{sur}}~.
\end{align}
Thus using the above definition of $N_{l_1}$, the condition to achieve NAQC for two flavor neutrino system is $N_{l_1} \geq 2.45$.
\begin{figure}[!htbp]
\begin{center}
\includegraphics[height=10cm,width=\textwidth]{Delta_NAQC_vs_phi.png}
 \end{center}
\vspace{-0.8cm}
\caption{\footnotesize  variation of $\Delta NAQC$ as a function of $\phi$ using approximations in fig. 6a. for a diagonal dissipator ($\alpha = a = c$, $\beta = \gamma = b = 0$), fig. 6b. one off-diagonal element $\beta$ non zero ($\alpha = a = c$, $\beta = \gamma = 0$, $\beta \neq 0$), fig. 6c. one off-diagonal element $b$ non zero ($\alpha = a = c$, $b = \gamma = 0$, $b \neq 0$) and in fig. 6d. $\alpha = a = c$, $\gamma=0$, $b$ and $\beta$ non zero}
\label{fig:Delta_NAQC_vs_phi}
\end{figure}
In fig.~\ref{fig:Bell-steering-NAQC-vs-E} we present the variation of the Bell-CHSH (top panel), steering inequality (middle panel) and NAQC inequality (bottom panel) with respect to neutrino energy E using eqn.~(\ref{bell}), eqn.~(\ref{steering}) and eqn.~(\ref{NAQC}) to investigate the quantumness of neutrino oscillations. The oscillation and the decoherence parameters considered are listed in table~\ref{table:osc-deco-params}.
In each row, we show the plots for three different baselines 295 km (left), 540 km (middle) and 1300 km (right). The solid lines in each plot are obtained by considering the neutrino propagation over a distance L, assuming decoherence and matter interactions while the dashed and the dot-dashed lines correspond to the cases of neutrino propagation with matter effect (no decoherence) and propagation in vacuum respectively.
From the first row of fig~\ref{fig:Bell-steering-NAQC-vs-E}, it can be seen that the Bell-CHSH inequality is violated i.e. $B_{CHSH}\geq 2$, for all E values in the case of all the three baselines. Second row of fig~\ref{fig:Bell-steering-NAQC-vs-E} shows that the flavor states violate steering inequality i.e. $S_3 \geq 1$ for all the energies in the case of T2K, ESS$\nu$SB and DUNE baselines. However, in the third row it is illustrated that NAQC i.e. $N_{l_1} \geq 2.45$ is not achieved for some energy ranges by the neutrino flavor states in the case of all the three baselines. Hence, neutrino flavour states attaining NAQC can also achieve Bell nonlocality and quantum steering. NAQC can be considered as a stronger quantifier than steering and Bell-CHSH quantity to reveal the quantumness of the neutrino system. Similar conclusions have been reported by Fei Ming et al~\cite{ming2020quantification} examining the neutrino oscillations in vacuum. Basing on this conclusion, we consider the NAQC quantifier as a better candidate to throw some light on the nature of neutrinos.

From eq.~(\ref{nue-app prob}) and eq.~(\ref{NAQC}) it is obvious that ${N}_{l_1}$ is a function of the Majorana phase ($\phi$). When we consider $\phi=0$ we can obtain $P_{osc}$ and thus ${N}_{l_1}$ for the case of Dirac neutrinos.
Here we propose a candidate $\Delta NAQC $, to determine the nature of neutrinos
    \begin{align}
        \Delta NAQC = | N_{l_1}(Dirac) - N_{l_1}(Majorana) |.
    \end{align}
We use the same values for the oscillation and decoherence parameters to calculate $\Delta NAQC$, as given in table II. In fig~\ref{fig:Delta_NAQC_vs_phi} we plot $\Delta NAQC$ as a function of $\phi$ imposing some extra conditions on decoherence parameters. We consider the following combinations for the sake of completeness: in fig \ref{fig:Delta_NAQC_vs_phi}a both $\beta$ and $b$ are zero (only non-zero diagonal elements); in fig~\ref{fig:Delta_NAQC_vs_phi}b, $\beta$ is non-zero and $b$ is zero (one non-zero off-diagonal element); in fig~\ref{fig:Delta_NAQC_vs_phi}c $b$ is non-zero and $\beta$ is zero (one non-zero off-diagonal element); and in fig~\ref{fig:Delta_NAQC_vs_phi}d both $\beta$ and $b$ are non-zero (two non-zero off-diagonal elements).
The orange, green, red, cyan and black lines are assigned for T2K, ESS$\nu$SB, NOvA, T2HKK and DUNE experiments respectively. In fig~\ref{fig:Delta_NAQC_vs_phi}, the $\Delta NAQC$ values in the presence of decoherence indicate that one can probe the nature of neutrinos using DUNE, NOvA and T2K but the possibility is low for the baselines T2HKK and ESS$\nu$SB. Additionally, it is visible from fig~\ref{fig:Delta_NAQC_vs_phi} that the $\Delta NAQC$ for 295 km baseline (orange curve) shows more sensitivity to the nature of neutrinos when compared to $\Delta P_{\mu e}$ plotted in fig~\ref{fig:diag-del-Prob-NH-phi-90} and fig~\ref{fig:del-Prob-NH-phi-b-and-beta-non-90}. Thus reinforcing the conclusion that NAQC, quantum quantifiers in general, can be used as a measure to understand the nature of neutrinos.
\section{Conclusions}
We analyse the two flavour neutrino oscillations in matter in a dissipative environment, by assuming two scenarios for the decoherence matrix. Firstly, for simplicity, we consider the case where we have non-zero diagonal elements and zero off-diagonal elements in the decoherent matrix. Whereas in the second case, we assume non-zero diagonal and non-zero off-diagonal elements. In these two scenarios, we see that the transition probabilities depend explicitly on the Majorana phase $\phi$. We further study the oscillations of Dirac and Majorana neutrinos at L (baseline) and E (neutrino beam energy) values corresponding to five long-baseline neutrino oscillation experiments T2K, NOvA, ESS$\nu$SB, T2HKK and DUNE. We find that, in principle, one can probe the nature of neutrinos at these experiments under the assumption that neutrino subsystem interacts with the environment.
Moreover, we extend our analysis to estimate the quantum correlations in the system using quantifiers - Bell-CHSH, steering and NAQC for these L, E values. We note that the neutrino flavour states achieving NAQC attain Bell-CHSH, steering inequalities for all E.
Thus, we consider NAQC as a stronger quantifier for a quantum system and analyse its sensitivity to the nature of neutrinos.
The NAQC quantifier carries different values for Majorana and Dirac neutrinos. Hence, the non-vanishing difference $\Delta NAQC$ can be used to discriminate between Dirac and Majorana neutrinos. In the presence of matter effect, a dissipative matrix presents non-vanishing $\Delta NAQC$, which in principle allows us to differentiate between the Dirac and Majorana neutrinos at neutrino oscillation experiments.

{\bf Acknowledgements:} K.~N.~Deepthi would like to thank DST-SERB SIRE program for the financial support to visit Indiana University, Bloomington, Indiana, USA.

\bibliographystyle{elsarticle-num}


\section{Introduction}
 Since by definition a gauge symmetry reduces to the identity on the observables, its physical meaning has been questioned and debated from a foundational/philosophical point of view.(For  the conceptual and philosophical discussions on the empirical meaning of gauge symmetries, which appear to be relevant for our  understanding of the  physical world,  see \cite{BC}).  Clearly, the final description of a physical system should be written in a gauge invariant language and therefore gauge symmetry are bounded to play only an intermediate role.

The following physical principles help to properly set the problem of understanding the strategic role of gauge symmetries. The {\em operational/experimental description of a physical system}
(not necessarily quantum!) makes reference to the following elements:

\vspace{1mm}
\noindent {\bf A}) the set  $\A$ of its measurable quantities, briefly
called {\bf observables};

\vspace{1mm}
\noindent {\bf B}) the measurement of the {\bf time evolution} of the observables;

\vspace{1mm}
\noindent {\bf C}) the set $\Sigma$ of configurations or {\bf states}  in which the
system may be prepared, according to well defined protocols of experimental preparations;
 the measurement of an observable $A$, when the system is in the state
$\om$, is  operationally defined by the {\bf experimental expectation} $< A  >_\omega$.

Given a state $\om$, the set of states  which can be prepared starting from $\om$, through physically realizable operations, is denoted by $\Gamma_\om$ and called the {\bf phase} of $\om$. This means that the protocols of preparations of states belonging to different phases are not related by physically realizable operations.
Clearly, by definition, different phases describe disjoint realizations of the system (or disjoint ''worlds'') which cannot communicate through realizable  observable operations.

The above \textit{operational framework} has an essentially unique \textit{mathematical  transcription}.
General physical considerations indicate that the set of the observables generate a normed algebra, actually a $C^*$-algebra with identity, called the {\bf algebra of observables}, for simplicity still denoted by $\A$. The time evolution defines a one-parameter  group of transformations (technically automorphisms) of $\A$: $t: \A \ra \alpha_t(\A)$.

 By  its experimental expectations, each state $\om$, defines a linear positive functional  $\om(A) \eqq <A>_\om, \, \forall A \in \A$, and by the GNS theorem a representation $\pi_\om$ of $\A$ (as operators $\pi_\om(A)$) in a Hilbert space $\H_\om$; $\om$ is represented by a vector $\Psi_\om \in \H_\om$ and the experimental expectations are represented by the matrix elements $<A>_\om = (\Psi_\om, \pi_\om(A) \, \Psi_\om)$. (For a  detailed discussion of the general description of a physical system, see \cite{FS12}).

The  states obtained by acting with   (the representatives of)  the  observables  on  (the representative of ) the state $\om$,  $\pi_\omega(\A) \, \Psi_\om$, are dense in $\H_\om$, which may be considered as the mathematical description of $\Gamma_\om$. In view of such a dual structure the family of  physical states on $\A$ is denoted by $\Sigma(\A)$.

In the following, we shall consider infinitely extended systems, which generically display  the occurrence of different phases and the related role of gauge symmetries (see below).  For this kind of systems a very relevant  property  is the {\bf local structure} of $\A$: each  observable   is identified by the experimental apparatus
used for its measurement and, since the physically realizable operations, as well as the corresponding experimental apparatuses, are inevitably localized in space, the algebra of observables is generated  by  localized observables .  This means that the algebras $\A(V)$ of observables localized in bounded regions  $V$, generate $\A$,  $\A = \cup_V \,\A(V)$, (hereafter called local algebra).

Another relevant property is the stability  of $\A$ under the group of {\bf space translations}  $\alpha_{\bf a}, \,\,{\bf a } \in \Rbf^3$, with  $\alpha_{\bf a}(A)$ denoting the ${\bf a} $-translated of $A$. (For lattice systems the space translations are replaced by the lattice translations).

  For a reasonable physical interpretation,
 the  measurement  of a localized observable $A$ should not  be influenced by measurements of  observables at infinite space separations, i.e. the following condition (called  {\bf asymptotic abelianess} must hold, ($B$ a localized observable):
  $\lim_{|\x| \ra \infty}\, [\,A, \,\alpha_\x(B)\,] = 0.$

A distinguished role is played by the \textit {pure homogeneous states} $\om_0$, characterized  by being \textit{invariant under a subgroup $\T$ of translations}
\be{ \om_0( \alpha_{\bf a} (A)) = \om_0(A), \,\,\,\forall A \in \A}\ee
and by satisfying the {\em cluster property}\index{cluster property}
\be{ \lim_{n \ra \infty}\,\om_0(A\,\alpha_{n \,{\bf a}}(B)) = \om_0(A)\,\om_0(B),}\ee
where $n \,{\bf a}$ denote the group parameters of $\T$. Typically, but not necessarily, $\om_0$ is a ground state. The set of such states $\om_0$ on $\A$ is denoted by $\Sigma_0(\A)$.

A very important consequence  is that:

\noindent $\om_0$ {\em is the unique $\T$ invariant state in the pure homogeneous phase $\Gamma_{\om_0}$}.

\vspace{1mm}
The interest of considering the phases  defined by such states  is twofold.

In relativistic quantum field theory, the vacuum state is  invariant under space translations and the validity of the cluster property is a necessary condition for the possibility of defining the scattering matrix, which requires the factorization of expectations of infinitely (space) separated clusters (describing scattering processes). Thus, the phase defined by a vacuum state in quantum field theory is a homogeneous pure phase in the above sense.

Quite generally,  in a homogeneous pure phase  the macroscopic observables, defined by space averages of local observables, take sharp (classical) values in agreement with the characteristic property of the standard pure phases in thermodynamics.

In conclusion, general physical principles characterize the mathematical description of a physical system, in terms of the local algebra of observables $\A$, its time evolution $\alpha_t(\A), t \in \Rbf$,  and the family of states $\Sigma(\A)$ in which the system may be prepared. The task of theoretical physics is to devise concrete effective strategies for implementing  and controlling such a general structure by  determining    the representations of the observable algebra corresponding to the various (physically realizable) states.

\def \Fo  {\F_{obs}}

In quantum field theory, as advocated by Wightman, the standard and well established  strategy is to study the vacuum representations of the relevant  field algebra, provided by the vacuum correlation functions. However, already in the case of a system of free particles,  physically important  particle states are not present in the vacuum representation of the observable algebra, called the vacuum sector.

A simple example  is provided by the system of free massive Dirac fermions, since the observable algebra and therefore its vacuum representation have zero fermionic charge ;  a complete  QFT description is obtained by introducing the  local field algebra $\F$ generated by the fermionic fields, whose (unique) vacuum representation contains  the representations of the  observable sub-algebra  $\Fo$, labeled by the quantum number of the fermionic charge.

This is  the standard strategy  adopted in the  formulation  and control of QFT models. One studies the vacuum representation of a local field algebra $\F$ generated by local fields, which describe the degrees of freedom inferred and borrowed  from the limit of vanishing interaction. Such a local field algebra contains a (suitably characterized) observable field sub-algebra $\Fo$ and  one looks  for the representations of $\Fo$ contained in a vacuum representation of   $\F$.

Therefore, the assumed role of the field algebra $\F$ is to provide the building stones of the theory: the fields of $\F$ describe the relevant states,  generate the observable algebra and allow for a direct and simple definition of the dynamics (through a Hamiltonian or Lagrangian which is a polynomial function of the fields of $\F$).

 In conclusion,  the difficult problem of determining the representations of the local observable algebra and its time evolution,  beyond the vacuum sectors, i.e. the solution of the problem $(\A, \alpha_t(\A), \Sigma(\A))$, is attacked  by finding the vacuum representations of a local field algebra $\F \supset \A$, i.e. by solving the relatively easier problem $(\F, \alpha_t(\F), \Sigma_0(\F))$.

Behind  such a strategic  choice there is  the implicit assumption (extrapolated from the non-interacting case) of considering  the states outside the vacuum sector which  are local states, i.e. may be obtained by applying the local fields of $\F$ to the vacuum; this property allows for their localizability in the DHR sense (for general discussion of such an important property, see \cite{Ha}):  the state $\om$ is localized in the double cone  $\O$ if
\be{ \om(A) = \om_0(A), \,\,\,\,\, \forall A \in \A(\O'), }\ee
where $\O'$ denote the set of points which are spacelike  w.r.t   $\,O$ and $\A(\O')$ the algebra of observables localized in $\O'$.

The quantum numbers which characterize the  inequivalent representations of $\A$ contained  in the Hilbert space $\H_{\om_0(\F)}$  of a vacuum representation of $\F$, define superselection rules; they   identify  a global gauge group $G$, which  defines non-trivial symmetries  of  $\F$ and   it is represented by the states of $\H_{\om_0(\F)}$ (if $\om_0(\F)$ is invariant under $G$).

Under general assumptions (namely the absence of infinite degeneracy  of particle types with equal mass and the completeness of the asymptotic states) the so derived global gauge group is compact. These deep results, which explain the {\textit {origin of compact gauge groups  from general physical principles}} were obtained in a series of papers by Doplicher, Haag and Roberts. (For a comprehensive  account see \cite{Ha}, esp. Section IV.4, and references therein.)

	In conclusion, the experimentally detectable existence of \textit{localizable states labeled by superselected quantum numbers}  imply that a local field algebra which generates them from the vacuum  has the symmetry of a global compact gauge group; this explains why such gauge groups  arise in the theory of infinitely extended systems, typically in quantum field theory, where locality plays a crucial role and vacuum representations of  local field algebras are the object of the standard approach.

More subtle is the problem of justifying  the role of local gauge symmetries on the basis of general physical principles, since they reduce to the identity both on the observables and on the (physical) states. There is no doubt that local gauge symmetries proved to be  useful  for the formulation of the quantum field theory of elementary particles, but an interesting question is whether their \textit{raison d'ètre}  may be traced back to general a priori arguments beyond their a posteriori successful use.

The standard justification of local gauge symmetries, in particular   in the prototypical case of quantum electrodynamics,  is tight to the introduction of redundant degrees of freedom through the vector potential, used either  for describing a massless spin one particle (the photon) or for  defining  the electron-photon interaction.\cite{We} \goodbreak

These reasons do not appear to be  rooted in general principles better than the origin of local gauge symmetries they should explain.
Actually, the photon is equally well described by the electromagnetic (quantum) field $F_{\mu\, \nu}$, whose Lorentz transformation does not require a gauge transformation.
As a matter of fact,  already for the solution of the  free field equations $\partial^\mu F_{\mu\, \nu} = 0$, the introduction of the quantum field $A_\mu$ as a Hilbert space operator is not free of problematic choices, since it cannot be done    without giving up covariance and/or locality (see \cite{FS16}, Chapter 7, Appendix 8).

   Moreover, a Lagrangian invariant under the group $\G$ of  local gauge transformations, labeled by the set of infinitely differentiable  gauge functions \textit{localized in space}, does not  define a deterministic time evolution of the field algebra and therefore $\G$ must be broken by a gauge fixing; such a breaking need not to be down to the identity, as might be inferred from the functional integral argument,  but only up to the extent of restoring  deterministic evolution  (see  \cite{FS19}).

In the standard  approach to gauge theories (particularly in the perturbative treatment) the identification and construction of the physical states is usually obtained by starting with a Lagrangian invariant under a group  $\G$ of local gauge transformations and  by adding a gauge fixing which preserves the invariance  under the related  global gauge group $G$ and allows to use a local field algebra $\F$. Then, the (infinitesimal) transformations of $G$ on $\F$ may be  generated by local currents  $J^a_\mu$, ($a = 1, ...n,  \,n  =$ the dimension of $G$).

As a consequence of the  gauge fixing, the second Noether theorem does not apply  and such currents do not satisfy a Local Gauss Law, i.e. they are not  the divergence of  antisymmetric tensors $ F^a_{\nu\,\mu}$. Nevertheless,  one may show that a characteristic property of the physical states constructed in a vacuum representation of $\F$ is that the currents $J^a_\mu$ satisfy Local Gauss laws  on them \cite{FS16, FS19}; i.e.
for any physical state $\Psi \in \H_{\om_0(\F)}$
\be{ (\Psi, (J_\mu^a - \partial^\nu F^a_{\nu\,\mu}\, \Psi) = 0,}\ee
(with $F^a_{\mu\,\nu}$ actually the field  strength) or, in a manifestly gauge invariant way,
\be{(\Psi, \sum_a (G^a_\mu)^*\, G^a_\nu \,\Psi) = 0, \,\,\,\,\,\,\,\,G^a_\mu \eqq J_\mu^a - \partial^\nu F^a_{\nu\,\mu}.}\ee
Such a form of the Local Gauss Law  is obviously satisfied in QED and may be easily checked to hold in the (local) BRST quantization of Yang-Mills theories, thanks to the nillpotency of the BRST charge. Actually, for the same reasons, any monomial  of  $G^a_\mu$ has vanishing expectation on the physical states. For brevity, such physical states  shall be called \textit{LGL states}. Clearly a prototypical  example is provided by the charged states in QED.

Thus, one may argue that even if local gauge symmetries reduce to thee identity on the physical states, they play an intermediate role in the construction of the representations of the observable algebra by guaranteeing that they are defined by  physical states obeying LGL.

Conversely, given the existence of LGL states,  one may investigate which gauge symmetries emerge for a local field algebra $\F$ which allows their construction through its vacuum representation; we shall argue that in this way one obtains local gauge symmetries.  This means that, in order to allow for the construction of LGL states, the local algebra  $\F$ must contains fields with non-trivial transformation under  a local gauge symmetry.

This offers a possible explanation of the emerge and strategic role of local gauge symmetries,  \textit{without  ever mentioning the vector potential}, which, being generically determined up to a (local) gauge transformation, automatically and obviously brings with it the freedom of  local gauge symmetries.

In this perspective local gauge symmetries are not introduced by an \textit{a priori} ansatz or Local Gauge principle, nor as a consequence of the introduction of  redundant degrees of freedom through the vector potential, but automatically arise as symmetries  of a local field algebra which  allows for  the realization of states satisfying Local Gauss Laws.

\section{Local Gauss Laws and local gauge symmetries}

The aim of this Section is to argue that  physical LGL states (e.g. the charged states in QED or  unconfined quark states in QCD) lead to local gauge symmetries for the  local field algebra used for constructing them.

In the abelian case (QED)  the characterization of the  physical LGL states is simply provided by the electrically charged  states, labeled by the superselected electric charge (giving rise to a $U(1)$ global gauge group) and satisfying the  LGL  given by the Maxwell equations.

 Less obvious  is the non-abelian case and we adopt the following characteristic properties of LGL states, extracted from their actual construction:

\noindent i) they carry superselected quantum numbers corresponding to a compact global gauge group $G$, as in the DHR case,

\noindent ii) in contrast with the DHR states, LGL states cannot  be described by local states in a vacuum representation of an auxiliary  local field algebra $\F$ which contains the field strengths $F_{\mu\, \nu}^a$, as well as the  currents $J^a_\mu$, which generate the infinitesimal transformations of $\F$ under $G$; for brevity,  the corresponding charges shall  be called \textit{Gauss charges}, to emphasize  that as a consequence of the LGL they are not localizable charges,

\noindent iii) in the abelian case of Quantum Electrodynamics (QED) they do not define local states on the algebra of the observables and therefore they  are not localizable in the DHR sense.

\vspace{1mm}
\noindent {\bf Remark 1}.\,\,
Property ii) follows from the fulfillment of LGL, if the vacuum representation of $\F$ satisfies semipositivity or the relativistic spectral condition.
In fact, eq. (1.5) with $\Psi = \Psio$, implies that $ G^a_\nu \,\Psio$ is a null vector and therefore, if  semipositivity holds,  $(\Psi, G^a_\nu\,\Psio) = 0$, $\forall \Psi \in  \H_{\om_0(\F)}$. As a consequence of this last equation, given  a local state $\Psi = F\,\Psio$, where $F \in \F$ transforms non-trivially under $G$,
\be{ \delta^a F = i  \limR [ \,Q^a_R, \, F\,],   \,\,\,\,\,Q^a_R = J^a_0(f_R\, \alpha),}\ee
(with the standard notation for the smearing of $J^a_0$,  see e.g. \cite{FS16}, Chapter 7, Section 2), by the locality of $F$ one has
$$ 0 \neq (\delta^a F\,\Psio, \delta^a F \Psio) = i \limR (\delta^a F \Psio,  [\,J^a_0(f_R \alpha) - \partial^i F_{i\, 0}(f_R \alpha), \,F\,] \Psio).$$ Then, if $F \Psio$, and therefore $\delta^a F \Psio$, satisfies LGL, the r.h.s. reduces to    $$ i \limR ( F^*\,\delta^a F \Psio,  G^a_0(f_R \alpha)\, \Psio) = 0, $$
leading to  a contradiction.

The same conclusion is reached even if the vacuum correlation functions of the auxiliary local field algebra $\F$ do not satisfy semipositivity, but the relativistic spectral condition   holds. In fact, as above, by locality for a local state $\Psi = F\,\Psio$ one has
$$ \limR < \Psi, [\, Q^a_R, F\,] \Psio> =  \limR < \Psi, [\,G^a_R, F\,] \Psio>,  $$
the limit is reached for finite $R$ and, if $\Psi $ satisfies   the LGL, the r.h.s reduces to $< \Psio, F^* \,F \,G^a_R \,\Psio >$, $R$ large enough.

Now, by the relativistic spectral condition  $$W(y-x, z-x) \eqq < \Psio, F^*(x) \,F(y) \,G^a_R(z) \,\Psio >,$$ (where $F(y) \eqq U(y) F  U(y)^{-1} $ denotes the $y$-translated of $F$, and similarly for the other operators), is the boundary value of a function $W(\zeta_1, \zeta_2)$ analytic in the tube $\T_2$. Such  a function vanishes for  $\zeta_1, \,\zeta_2$ real, $\zeta_2 = (\z_2, z_{0, 2})$, $|\z_2|$ sufficiently large, since then   $[ \,F(y), \,G^a_R(z)\,] = 0$ and LGL applies. Then, by the edge of the wedge theorem $W(\zeta_1, \zeta_2) = 0$ and $F \,\Psio$ is chargeless, contradicting i).

\vspace{1mm}
\noindent {\bf Remark 2}.\,\,In the QED case, the current, which generates the global  gauge group $U(1)$, and the field strength $F_{\mu\, \nu}$ are observable fields; then, if  a state  $\om$ is localized in the DHR sense, say in a double cone $\O$, for sufficiently large $R$, $\partial^i F_{0 \, i}(f_R \alpha) \in \F_{obs}(\O')$, so that, according to the DHR criterion,    $\om(\partial^i F_{0 \, i}(f_R \alpha) ) = \om_0(\partial^i F_{0 \, i}(f_R \alpha)) =0$ and by the LGL $\om$ is chargeless.

\vspace{2mm}
The embedding of the observable algebra into a larger algebra $\F$ for the description of LGL states, through  a vacuum representation of $\F$, is not unique and the local gauge symmetries of $\F$ depend on $\F$. As we shall see, the generators of the infinitesimal  local gauge transformations on $\F$ are  provided by the Gauss operators $G^a_0$.    A limiting case in QED is the choice of the Coulomb gauge field algebra $\F_C$, since in $\F_C$ the Gauss operator $G_0$ vanishes, $\F_C$ is non-local  and  the Coulomb  gauge fixing excludes any local gauge symmetry of $\F_C$, (for the general structure of the Coulomb gauge see \cite{Di, Sy}. For the necessary ultraviolet regularization, see, for a perturbative control, \cite{Ste} and, for a general control which exploits the properties of the Feynman-Gupta-Bleuler (FGB) gauge, \cite{BDMRS}; see also the discussion in \cite{FS19} Section 2.3).

Actually,  a crucial property for the strategy outlined above, leading to local gauge symmetries,  is the locality of the field algebra; this choice, motivated by the no-interaction limit,  is also suggested by technical reasons, since locality helps for the control of the dynamics (it is well known that the renormalizable gauges are local gauges).
Also from a constructive point of view, the infinite volume limit is better handled  for a local field algebra with a local dynamics.\footnote{For the problems arising in the case of  a non-local dynamics, see e.g. \cite{FS21}, Appendix A.}

For these reasons,  the auxiliary local field algebra $\F$  should not satisfy the LGL; otherwise, by locality, $\F$ would be pointwise invariant under the global gauge group $G$.
In fact, if the LGL hold in the local algebra $\F$ one would have
$$ \delta^a F = \limR [\, J^a_0(f_R \alpha), \, F\,] = \limR [\, \partial^i F^a_{i \,0}(f_R \alpha), \, F\,] = 0,\,\,\,\,\,\forall F \in \F, $$
 and, consequently,  $\H_{\om_0(\F)} = \overline{\F \, \Psio}$ would not contain states carrying  non-trivial charges of $G$.

   The different choices of the local algebra $\F$ are characterized by the way LGL fail, i.e. by the non-vanishing Gauss operators $G^a_0$. To this purpose, it is useful to remark that LGL correspond to a combination of hyperbolic evolution equations and constraint/elliptic  equations for the field strengths $ F^a_{\mu \,\nu}$:
\be{ \square F^a_{\mu \,\nu} = \partial_\mu J^a_\nu - \partial_\nu J^a_\mu, \,\,\,\,\,\,\,\,\,\,\,\,\,\,\partial^i F^a_{i \,0} = J^a_0.}\ee

If, as required, LGL do not hold in $\F$, the above equations get modified by the non-vanishing Gauss operators   $G^a_\mu$,\,\, $ \partial^\mu G^a_\mu =0$.

\vspace{1mm}
\noindent a) {\bf{ \textit{Time independent local gauge symmetries}}}

As remarked above, a non-trivial representation of the global gauge group $G$ by a local field algebra $\F$  requires that $G^a_0 \neq 0$, since, by locality  the infinitesimal transformations of $\F$  by $G$ are given by
$$\delta^a F = \limR [\, J^a_0(f_R \alpha), \, F\,] = \limR [\,J^a_0(f_R \alpha) -  \partial^i F^a_{i \,0}(f_R \alpha), \, F\,] =$$
\be{ =  \limR [\, G^a_0(f_R \alpha), \, F\,].}\ee
 A  possible realization of such a framework,  which, in a certain sense, minimizes the needed violation of the LGL in $\F$, even at the expense of loosing manifest covariance, is given by  $G^a_0 \neq 0$,
$ G^a_i  = 0$, which plays the role of a gauge fixing; then, the continuity equation for $G^a_\mu$ requires  that $G^a_0(\x, t)$  is time independent.
As a consequence, eqs.\,(2.2) are replaced by
\be{ \square F^a_{i \,j} = \partial_i J^a_j - \partial_j J^a_i,\,\,\,\,\,\,\,\,\,\,\partial^0 F^a_{0\, i} =  - \partial^j F^a_{j \, i} + J^a_i,}\ee
with the equal time constraint
\be{\partial^i F^a_{i\, 0 } =   J^a_0 - G^a_0.}\ee
This is the choice adopted by the temporal gauge; in fact, eqs.\,(2.4, 2.5) are the evolution equations of the fields $F^a_{\mu \,\nu}$ in the temporal gauge.

Equations (1.5) imply that  for any physical state $\Psi$, $(G^a_0 \,\Psi, \,G^a_0\, \Psi) = 0$, i.e., if positivity holds,
$G^a_0\, \Psi = 0$.

One should remark that a mathematical subtlety occurs in the standard local and positive   temporal gauge, namely    the local field algebra $\F$ represented by a vacuum state  is generated by the exponentials  of the standard (non-observable) fields, with algebraic relations corresponding to those of their formal generators. Then, the condition which selects  the physical states   should rather read $$V^a(\Lambda) \Psi = \Psi,\,\,\,\,\,	\Lambda  \in \D(\Rbf^3),$$
 where the unitary field operator  $ V^a(\Lambda) $ is formally the unitary  exponential of $G^a_0(\Lambda)$.

Furthermore, since $\F$  is not covariant under  relativistic transformations the relativistic spectral condition is not satisfied by the vacuum correlation function of $\F$, the Reeh-Schlieder theorem does not apply and $G^a_0 \Psio = 0$, or better $(V^a(\Lambda) - \id) \, \Psio = 0$ does not imply that the local operator $G^a_0$ or better $(V^a(\Lambda) - \id)$ vanishes.\footnote{ For a more detailed discussion see  \cite{FS16} Chapter 8, Section 2.1.}

 For simplicity, in the following discussion, we shall sometimes use  the formal generators of the exponential fields, the more accurate mathematical discussion being easy to obtain.

The operator $G^a_0$ plays the role of a static charge density which is not seen by the physical states and compensates the vanishing flux of $ F^a_{i \,0}$  at  (space) infinity on the local states of $\H_{\om_0(\F)}$.

The non-vanishing  $G^a_0$ modifies the equal time constraint implied by the LGL on the local states; the  vanishing of $G^a_0$ on the physical states requires that the physical states of $\H_{\om_0(\F)}$, carrying a non-trivial charge of $G$, are non-local limits of local states.

Now, for any test function $\Lambda(\x) \in \D(\Rbf^3)$, the operator $G^a_0(\Lambda, t) $ generates a \textit{time independent derivation} on $\F$
\be{ \delta^{a, \Lambda} F \eqq i [\, G^a_0(\Lambda, t), \,F\,] = i [\, G^a_0(\Lambda, h), \,F\,], \,\,\,\,\,\,h \in \D(\Rbf), \,\,\,\,\,\int d t\,h(t) = 1. }\ee
This is the basic property of the derivations generated by local (covariant) currents, corresponding to infinitesimal symmetry transformations, where the time independence is a consequence of current conservation and locality.

If the auxiliary local field algebra $\F$ contains the formal exponentials of $G^a_0$, (as in the standard temporal gauge), the following local transformations are defined on $\F$:
 \be{ \b^\Lambda(F) = V^a(\Lambda)\, F \,V^a(\Lambda)^*.}\ee

Moreover, since  the subspace of physical state vectors  must be pointwise invariant under the application of the observable operators, these operators  should commute  with $V^a(\Lambda)$, so that the above  transformation defines a \textit{local time independent gauge transformation}.

As a matter of fact, in the temporal gauge, where  the local algebra $\F$ is generated by canonical fields, the derivation (2.6) corresponds to  the standard time  independent gauge transformations, with local gauge parameter $\Lambda(\x)$.

The time independence of the generators $G^a_0$ implies that the time evolution of $\F$ commutes with such local gauge transformations; thus, the Hamiltonian, as a function of the fields of $\F$, must  be gauge invariant, a property which in the standard approach corresponds to the requirement of minimal coupling.

Thus, by the above arguments,  the strategic role  of local gauge symmetries may be traced back  to the realization of LGL states through a vacuum representation of a local field algebra $\F$, where the necessarily \textit{non-vanishing  Gauss operators  generate  local gauge symmetries}.\goodbreak

\vspace{10mm}
\noindent b) {\bf \textit{Local gauge symmetries in QED}}

Another distinguished example of local gauge symmetries, arising according to the pattern discussed above, is provided by Quantum Electrodynamics.\footnote{For a general discussion of the occurrence  of local gauge symmetries in QED, in a  C*-algebra setting, see F. Ciolli, G. Rizzi, E. Vasselli, QED Representation for the Net of Causal Loops, Rev. Math. Phys., {\bf{27}}, 1550012 (2013);   D. Buchholz, F. Ciolli, G. Rizzi, E. Vasselli, The universal algebra of the electromagnetic field. III. Static charges and emergence of gauge fields, arXiv: 2111.01538 [math-ph]. }

This example is obtained by requiring that the time evolution of the observable  electromagnetic field as operator in $\F$ is not modified by the non-vanishing Gauss operator $G_\mu$, i.e. that
\be{  \square F_{\mu \,\nu} = \partial_\mu J_\nu - \partial_\nu J_\mu,}\ee
the only effect of $G_\mu \neq 0$ being a modification of the equal time constraint
\be{\partial^i F_{i \,0} = J_0 - G_0.}\ee
This implies that $\partial_\mu G_\nu - \partial_\nu G_\mu = 0$ and in a \textit{Lorentz-covariant local field algebra} this is obtained by  $G_\mu  = \partial_\mu \L$, with $\L(x)$ a scalar field of $\F$; then, the continuity equation obeyed by $G_\mu$ implies  $\square \L(x) = 0$.

Such a choice of the local field algebra corresponds to the  Feynman-Gupta-Bleuler (FGB) gauge, albeit in a more general context, since no reference is made to the vector potential; in fact,   eqs.(2.7), (2.8) coincides with the equations for $F_{\mu\, \nu}$ in  that gauge.

The relativistic covariance of $\F$ lead to  the validity of the Reeh-Schlieder theorem and therefore positivity cannot hold, since otherwise eq.\,(1.5) implies $G_\mu \,\Psio = 0$; hence, by the Reeh-Schlieder theorem the local operator  $G_\mu$ vanishes, $\F$ commutes with the charge and $\H_{\omega_0(F)}$ does not contain charged states.

 A way out of this difficulty is to guarantee the validity of eq.\,(1.5) on the vacuum state by a non-local condition; since $\L(x)$ is a free field, its negative energy part $\L(x)^-$ is well defined, it is a non-local operator and  the equation
\be{ \partial_\mu \L(x)^- \, \Psi   = 0}\ee
implies that $\Psi$ satisfies eq.\,(1.5); then, such a condition may  be chosen for  selecting the physical states.

As in the temporal gauge discussed above for the general case of a compact gauge group $G$, one may show that in QED the Gauss operator generates a time independent derivation on the local field algebra $\F$.

To this purpose, one considers  infinitely differentiable functions $\Lambda(x)$ satisfying
\be{ \square \Lambda(x) = 0,  \,\,\,\,\,\, \,\,\,\,\,\,\,\Lambda(\x, 0),\,\,\partial_0 \Lambda (\x, 0) \in \D(\Rbf^3), }\ee
and  the operator
\be{ G(\Lambda)_{R}(x_0) \eqq   \int d^3 x \,f_R(\x)\,( \Lambda(x) \stackrel{\leftrightarrow}{ \partial _0} \L(x)),}\ee
 with $(A \stackrel{\leftrightarrow}{ \partial _0} B) \eqq (A \partial_0 B - \partial_0 A \,B)$.

Since  $\Lambda(x) $ and  $\L(x)$ satisfy the free wave equation, by  locality the following commutator   is independent of time, i.e.
\be{ i \partial_0 \limR [ G(\Lambda)_{R}(x_0), \, F] = 0,\,\,\,\,F \in \F.}\ee
In fact,
$$\partial_0 \limR [ G(\Lambda)_{R}(x_0) =   \int d^3 x\, f_R(\x) (\Lambda(x) \Delta \L(x) - \Delta \Lambda(x)  \L(x)) = $$
$$ =  \int d^3 x\, \partial_i f_R(\x) (\Lambda(x) \partial_i \L(x) - \partial_i \Lambda(x)  \L(x)).$$
Then, since supp $\partial_i f_R(\x) \subset (R \leq |\x| \leq R(1 + \eps))$, by locality the commutator (2.11) vanishes.
In conclusion, one has  a time independent derivation on $\F$, labeled by the  infinitely differentiable functions $\Lambda(x)$ of compact support in space
\be{ \delta^\Lambda \,F \eqq \int d x_0\, \alpha(x_0) i \limR \,[\, G(\Lambda)_R(x_0), \, F\,].}\ee

The stability of the subspace of physical states under application of observable operators is guaranteed if the observables commute with $\L(x)^-$, (and therefore  with $\L(x)$). Then the above derivation has the meaning of  an infinitesimal local gauge transformation with gauge function $\Lambda(x)$.

Indeed, if the algebra $\F$ is generated by local canonical fields, as it is the case of the FGB realization, one gets the standard local gauge symmetries   of the FGB gauge.

As discussed above in point a), the time independence  of the derivation (2.14) implies that, as a   function of the local fields of $\F$, the Hamiltonian should be invariant under local gauge transformations with gauge parameter $\Lambda(x)$, (satisfying $\square \Lambda = 0$).

\section{Conclusion}

In conclusion, a possible physical/empirical explanation of gauge quantum field theories is the existence of (particle) \textit{states which carry the quantum numbers of a  (global) compact  gauge group $G$  and satisfy Local Gauss Laws}, for the conserved currents associated to $G$. The prototype is clearly  Quantum Electrodynamics where the electric charge (the generator of the global compact $U(1)$ gauge group)  labels  the (physical) charged states, which  satisfy the  Local Gauss Law corresponding to the Maxwell equations.  Quite generally, the  description of such  LGL states through a vacuum representation of a {\bf \textit{local}} \textit{field algebra} $\F$  leads to the emergence  of {\bf \textit{local}} \textit{ gauge symmetries} for $\F$, which commute with the time evolution of $\F$.

 Thus, whereas   states carrying localizable superselected charges lead to global compact  gauge  groups  for  the local field algebra which obtains them from the vacuum,  the realization of  states carrying  Gauss charges through a  vacuum representation of a local field algebra $\F$ implies a  local gauge symmetry for $\F$,  i.e. such a field algebra  must contains fields which transform non-trivially under local gauge transformations. Therefore, the local extension of the global gauge group  $G$ (related to the superselected quantum numbers) needs not to be \textit{ a priori} assumed with no compelling physical motivations (as in the standard   characterization/definition of  gauge theories),  but it is required by the \textit{physical} existence  of states carrying Gauss charges and their description/construction through the vacuum representation of a local field algebra.

\vspace{5mm}
\noindent
Acknowlegements.
I am indebted to Francesco Serra for a stimulating discussion.

 \newpage