\section{Introduction\label{introduction}}

    Dialogue state tracking (DST) is a core component in modular task-oriented dialogue systems \cite{hosseini2020simple, yang2021ubar, sun-etal-2022-bort,sun-etal-2022-mars}. It extracts users' intents from the dialogue history and converts them into structured dialogue states, i.e., sets of slot-value pairs.

	An accurate dialogue state is crucial for generating correct dialogue action and suitable natural language responses, which are the main tasks of dialogue management and natural language generation components \citep{williams2007partially, thomson2010bayesian, young2010hidden}.

	Earlier DST approaches predict the state directly from the dialogue history (natural language utterances) \citep{mrkvsic2017neural, xu2018end, wu2019transferable, chen2020parallel}.
	Since the dialogue state is condensed structural information memorizing all historical information, recent methods incorporate the previously predicted state as the input besides the dialogue history.

    \citep{ouyang2020dialogue, kim2020efficient,ye2021slot}.

    Conventional DST models taking the previous state as the input usually show the characteristic that the previously predicted slot values tend to be kept unchanged when predicting the current state, defined as \textit{state momentum} in this paper.
    The state momentum makes DST models struggle to modify the previous prediction, which affects the performance when the values of some slots need to be \textit{updated} as the user's intent changes, and there exist wrongly predicted slot values that need to be \textit{corrected}.

    Figure \ref{diagexm} shows an example of a dialogue involving three turns involving two state momentum issues.

    State momentum (1): the predicted state ``hotel-book day-\textcolor{red}{Saturday}'' keeps unchanged from the previous state ``hotel-book day-\textcolor{blue}{Saturday}'', which should be \textit{updated} into ``hotel-book day-Sunday'' in the second turn.
    State momentum (2): the predicted state ``hotel-book day-\textcolor{red}{Saturday}'' keeps unchanged from the previous state ``hotel-book day-\textcolor{red}{Saturday}'', which should be \textit{corrected} into ``hotel-book day-Sunday''

    in the third turn.

    By using the general DST model (the baseline model in Section \ref{base_arc}) to make predictions on the dataset MultiWOZ 2.4, there are around 32.4\% errors belonging to state momentum which accounts for a large proportion of bad cases, indicating that the solving state momentum issues is crucial to the DST task.

	In this paper, we propose \textbf{MoNET} to tackle the state
    \underline{mo}mentum issue via a \underline{n}oise-\underline{e}nhanced \underline{t}raining strategy.
    First, the previous state of each turn in the training data is noised via replacing some of its slot values.

    Specifically, for each active slot (with a non-$none$ value), we replace its value with a certain probability.

    Then, the noised previous state concatenated with the dialogue history is used as the input to learn to predict the current state, improving the model's ability to \textit{update} and \textit{correct} slot values.
    Furthermore, a contrastive context matching framework is designed to narrow the representation distance between a state and its corresponding noised variant, which reduces the impact of noised state and helps the model better understand the dialogue history.

    Such approaches make the model ``less insensitive'' to the noise, enhancing the model's ability to modify the slot values of previous states in current predictions.

	Experiments on the multi-domain dialogue datasets MultiWOZ 2.0/2.1/2.4 show that our MoNET outperforms previous DST models.
	Ablations and analysis verify the effectiveness of the proposed noised DST training and the contrastive context matching framework in alleviating state momentum and improving anti-noise ability.

	The contributions are summarized as follows:

	(1) We define the state momentum issue in DST, where DST models tend to keep the predicted slot values unchanged, namely, struggling to \textit{update} and \textit{correct} them in the last turn.
	(2) We propose MoNET to tackle the state momentum issue via noised DST training and the contrastive context matching framework.
	(3) We conduct comprehensive experiments on three datasets, i.e., MultiWOZ 2.0/2.1/2.4, and the results of MoNET outperform previous DST methods. We will release our code to further the DST research.

	\begin{figure*}[ht]
			\centering
			\includegraphics[width=15.5cm]{fig/MoNET8}
            \caption{The model description and noised input construction example. The left part (a) shows the architecture of the MoNET model. A context input representation $H^n_t$ in an $N$-size batch is shown in the contrastive context matching framework as an example, where $H^{n+}_t$ is the context representation of its noised variant.
            The right part (b) gives an example of constructing noised context input. For each active slot in state $B_{t-1}$, given a noise threshold $p$, a random number $a$ is selected. If $a < p$, then its value is replaced with another one randomly selected from the ontology (suppose the pair ``train-day-Saturday'' is replaced with ``train-day-\textcolor{red}{Friday}''); otherwise the value will be kept unchanged (suppose the pairs ``train-departure-Birmingham'' and ``train-destination-Cambridge'' are unchanged).}

			\label{fig:STC}
		\end{figure*}

	\section{Related Work}
	\subsection{Dialogue State Tracking}
	Traditional DST approaches focus on single-domain dialogue state tracking \citep{williams2007partially,thomson2010bayesian,lee2016dialog}.

	Recent researches pay more attention to multi-domain DST using distributed representation learning \citep{wen2017network,mrkvsic2017neural}.
	Previous works implement Seq2seq frameworks to encode the dialogue history, then predict the dialogue state from scratch at every turn \citep{rastogi2017scalable,ren2018towards,lee2019sumbt, wu2019transferable, chen2020parallel}. Utilizing dialogue history is limited for larger turns, since the state of each turn is accumulated from all previous turns, while it's hard to retrieve state information from a long history.

	Current works mainly incorporate the previous state as the model input, which is regarded as an explicit fixed-sized memory \citep{ouyang2020dialogue, ye2022assist, wang-etal-2022-luna}.
    \citet{kim2020efficient} propose a state operation sub-task, where the model is trained to first predict the operation of each slot-value pair, such as UPDATE, CARRYOVER, etc., then only the value of a minimal subset of slots will be newly modified

    \citep{zeng2020jointly,zhu2020efficient}. These methods enhance model prediction efficiency and the ability to \textit{update} slot-value pairs.

    \citet{tian2021amendable} deal with the error propagation problem that mistakes are prone to be carried over to the next turn, and design a two-pass generation process, where a temporary state is first predicted then used to predict the final state, enhancing the ability to \textit{correct} wrong predictions.

    In this paper, we use ``state momentum'' to define the issue where the wrong dialogue state is predicted due to that the previous prediction keeps unchanged, either it should be \textit{updated} or \textit{corrected}.
    To the best of our knowledge, this is the first time to systematically tackle the issue caused by continuous unchanged predictions in the multi-turn DST task.

	\subsection{Contrastive Learning}
	Contrastive learning aims to generate good-quality representations by constructing a set of paired examples and learning the similarity \citep{hadsell2006dimensionality,mnih2012fast,baltescu2015pragmatic,peters2018deep}.
	The goal is to help the model semantically divide the augmented data samples into different groups.
	During training, the neighbors with similar semantic representations (\textbf{positive pairs}) will be gathered, while the non-neighbors (\textbf{negative pairs}) will be pushed apart, and then the model will learn better representations.
    In NLP area, semantic representations can be learned through self-supervised methods, such as center word prediction in Word2Vec, next sentence prediction in BERT, sentence permutation in BART, etc \citep{mikolov2013efficient, devlin2019bert, lewis2020bart}.
	Recent approaches build augmented data samples through token shuffle, word deletion, dropout, and other methods \citep{cai2020group,klein2020contrastive,yan-etal-2021-consert,wang2021cline,gao2021simcse}.
	In this paper, we construct augmented samples based on the noised and original dialogue state. Given context inputs with the same dialogue history and different states, the model is trained to gather them into similar objects, aiming to learn better representations, reduce the impact of noise state and better understand the dialogue history.

	\section{Methodology}
	\subsection{Problem Formulation}

	In this paper, we focus on building a dialogue state tracking (DST) model which predicts a current dialog state based on the dialogue history and the previous state during a multi-turn dialogue interaction.
	A dialogue state consists of domain-slot-value tuples, typically corresponding to the dialogue topic, the user's goal, and the user's intent.
	Following previous studies, in the rest of this paper, we omit ``domain'' and use ``slot'' to refer to a ``domain-slot'' pair.
	All slot-value pairs are from a pre-defined ontology.

	Formally, let's define $D_t = [U_t, R_t]$ as a pair of system utterance $U_t$ and user response $R_t$ in the $t$-th turn of a multi-turn dialogue, and $B_t$ as the corresponding dialogue state.
	Each state $B_t$ contains a set of slot-value pairs, i.e., $B_t=\{(S_j,V^{i}_{j}) | j \!\in\! [1\!:\!J]\}$, where $J$ is the number of all slots and $V^{i}_{j} \in \mathcal{V}_j$ is one of the values in $\mathcal{V}_j$ for the $j$-th slot $S_j$ in the ontology.
	Given the dialogue history $\{D_1,..., D_t\}$ and previous state $B_{t-1}$, the goal of DST task is to predict the current state $B_t$.

	\subsection{MoNET}
	As introduced in Section~\ref{introduction}, solving the \textit{state momentum} issue is crucial to the DST task.
	Therefore, in this paper, we propose MoNET to tackle the state
    \underline{mo}mentum issue via a \underline{n}oise-\underline{e}nhanced \underline{t}raining strategy to enhance the model's ability to \textit{update} and \textit{correct} slot values.
    The architecture of MoNET is shown in Figure \ref{fig:STC}(a), which consists of context BERT encoders, slot and value BERT encoders, the slot-context attention module, the slot-value matching module, and the contrastive context matching framework. Each of them will be elaborated on in this section.

	\subsubsection{Base Architecture}
	We first introduce the base architecture of our MoNET.
	A model trained only with the base architecture of MoNET is noted as ``Baseline'' and evaluated in Section \ref{5} to compare the difference in performance with the whole MoNET model.
	\label{base_arc}

	\paragraph{Context Encoder.} A BERT encoder encodes the context input, which is the concatenation of dialogue history and state in the last turn:
	\begin{equation}
	\resizebox{0.98\hsize}{!}{
	$
	\begin{aligned}
	&X_t = f(M_t, B_{t-1}, D_t) \\
	&= [CLS] \oplus M_t \oplus B_{t-1} \oplus [SEP] \oplus D_t \oplus [SEP],
	\end{aligned}
	$}
	\end{equation}
	where $M_t=D_1 \oplus ,...,\oplus D_{t-1}$ contains previous utterances, $B_{t-1}$ is the state in the last turn (\textcolor{black}{contains the active slots only}), $[CLS]$ and $[SEP]$ are the special tokens of the BERT encoder. Then the representations of the context input are derived:
	\begin{equation}
		H_{t} = BERT(X_t) \in \mathbb{R}^{|X_t| \times d},
	\end{equation}
    where $|X_t|$ is the total number of tokens in $X_t$, and $d$ is the encoded hidden size.
	\paragraph{Slot and Value Encoders.}
	The BERT encoders with fixed parameters are used to get the slot and value representations:
	\begin{equation}
	\begin{aligned}
		&h_{S_j}=BERT_{fixed}(S_j)_{[CLS]}, \\
		&h_{V^{i}_{j}}=BERT_{fixed}(V^{i}_{j})_{[CLS]},
	\end{aligned}
	\end{equation}
	where $h_{S_j}, h_{V^{i}_j} \in \mathbb{R}^{1 \times d}$ are the $[CLS]$ representations of the slot and value.

	\vskip 0.5em \noindent
	\paragraph{Slot-Context Attention.}
	For each slot $S_j$, its slot-context specific feature is extracted by the multi-head attention mechanism \citep{vaswani2017attention}:
	\begin{equation}
	    r^t_{S_j}=LN(MultiHead(h_{S_j}, H_t, H_t)) \in \mathbb{R}^{1 \times d},
	\end{equation}
	where $LN$ is the normalization layer.

	\paragraph{Slot-Value Matching.}
	The probability of predicting the value $V^{i}_{j}$ of the slot $S_j$ is derived by calculating the $L2$-distance between the value representation $h_{V^{i}_{j}}$ and the slot-context representation $r^{t}_{S_j}$, which is denoted as:
	\begin{equation}
	\resizebox{0.95\hsize}{!}{
		$P_{\theta}(V^{i}_{j}|X_t, S_j) = \frac{exp(-||r^{t}_{S_j}-h_{V^{i}_{j}} ||_2)}{\sum_{k \in [1:|\mathcal{V}_j|]}exp(-||r^t_{S_j}-h_{V^{k}_{j}}||_2)}$,}
	\end{equation}
	where $\theta$ are trainable parameters of the model.
	\paragraph{Training and Inference.}
	During training, the oracle dialogue state is used to form the context input $X_t$ (teacher forcing). For the $t$-th turn, the loss is the sum of the negative log-likelihood among all $J$ slots as follows:
	\begin{equation}
		L_{d_{ori}} = \sum_{j=1}^J -log(P_{\theta}(V^{i^{*}}_{j}|X_t, S_j)),
	\end{equation}
	where $V^{i^{*}}_{j}$ is the ground truth value of the slot $S_j$ at turn $t$.
	During inference, the previously predicted state is used to form the context input $X_t$, and the value of the slot $S_j$ is predicted by selecting the one with the smallest distance, i.e., the largest probability:
	\begin{equation}
	    {V^{\widehat{i}}_{j}} = \mathop{argmax}_{i \in [1:|\mathcal{V}_j|]} P_{\theta}(V^{i}_{j}|X_t, S_j).
	\end{equation}

    \subsubsection{Noised Data Construction}
    \label{nos_thresh}
    By using the Baseline model to make predictions on the dataset MultiWOZ 2.4, there exists around $32.4\%$ errors belonging to state momentum which accounts for a large proportion of bad cases, indicating that solving state momentum issues is crucial to DST.
    We think the reason behind this is that most slots' values in the $(t-1)$-th turn

    are the same as the values in the $t$-th turn in the training data, which makes conventional DST models lack the ability to \textit{update} or \textit{correct} slot values during inference.
    An intuitive idea is to augment training instances with a higher ratio of slots whose values in the $(t-1)$-th turn are different from those in the $t$-th turn.

    Based on this point, in this paper, we propose to tackle the issue by utilizing noised data to train the DST model.

    Generally, for each active slot (with a non-$none$ value) in the previous dialogue state, we involve noise by replacing its value with another value with a probability $p$ (used as the noise threshold), e.g., the example shown in Figure \ref{fig:STC}(b).

	Formally, at each training step, given a batch of training instances, a noised context $X^+_{t}$ is constructed for each instance based on its original context input $X_t=f(M_t, B_{t-1}, D_t)$ as follows:
	\begin{equation}
		\begin{aligned}
		     X^+_{t}=& f(M_t, B^+_{t-1}, D_t), \\

		  B^+_{t-1}=&\{(S_j,V^{i^+}_{j}) | j \!\in\! [1\!:\!J]\}.
		\end{aligned}
		\label{lable_pos_sample}
	\end{equation}
	For each active slot $S_j$ in $B_{t-1}=\{(S_j,V^{i}_{j})\}$, a real number $a\in[0,1]$ is sampled to determine whether the original $V^{i}_j$ is replaced with a randomly selected value $V^{k}_{j} \in \mathcal{V}_j \setminus \{V_{j}^{i}\}$ from the ontology:
    \begin{equation}
    \begin{aligned}

            V^{i^+}_{j}= \left \{
                \begin{array}{ll}
                    V^{k}_{j}, &if \quad a < p \\
                    V^{i}_{j}, &if \quad a \geq p.
                \end{array}
                    \right.

    \end{aligned}
    \end{equation}

    \subsubsection{Noised State Tracking}
	\label{noised_DST}
    Like $X_{t}$, the noised context input $X_{t}^+$ is also used as the input to predict the state $B_{t}$ as the training target, aiming to improve the model's ability to \textit{update} and \textit{correct} slot values.
	Specifically, the representation $H_t^+$ of $X^+_t$ is first derived by the BERT context encoder mentioned in Section \ref{base_arc}:
	 \begin{equation}
	    H_t^+ = BERT(X^+_{t})  \in \mathbb{R}^{|X^+_{t}| \times d}.
	 \end{equation}
	Then, similar to the previous process, for each slot $S_j$, $X^+_t$ is used to predict its value based on the distribution $P_{\theta}(V^{i}_{j}|X^+_t, S_j)$.
	Eventually, the loss of the noised state tracking can be denoted as:
	\begin{equation}
		L_{d_{nos}} = \sum_{j=1}^J -log(P_{\theta}(V^{i^*}_{j}|X^+_t, S_j)).
	\end{equation}

	\subsubsection{Contrastive Context Matching}

    Inspired by contrastive learning approaches which group similar samples closer and diverse samples far from each other, a contrastive context matching framework is designed to narrow the representation distance between $X_t$ and its noised variant $X_t^+$, aiming to reduce the impact of noised state $B_{t-1}^+$ and help the model better understand the dialogue history.
	Specifically, in a batch of $N$ instances with the original context input $\mathbf{X}_t = \{X_t^n\}_{n=1}^{N}$, we construct $N$ corresponding noised instances with the context input $\mathbf{X}_t^+ = \{X_t^{n+}\}_{n=1}^{N}$.
	Note that, in this section, we involve $n$ into $X_t$ (or $H_t$) as $X_t^n$ (or $H_t^n$) to indicate the index in a batch.
	For each context input $X_t^n$, the noised sample $X_t^{n+}$ is regarded as its positive pair, and the rest $2N\!-\!2$ instances in the same batch with different dialogue histories are considered to be its negative pairs.
	Then the model is trained to narrow the distance of the positive pair and enlarge the distance of negative pairs in representation space with the following training objective:
	\begin{equation}
	    \resizebox{0.98\hsize}{!}{
	    $L_c = -log( \frac{exp(sim(H^{n[cls]}_{t}, H_t^{n+[cls]})/\tau)}{\sum_{k=1}^{2N}\mathbbm{1}_{[k\neq n]}exp(sim(H_t^{n[cls]}, H_t^{k[cls]})/\tau)})$,}
	\end{equation}
	where $H_t^{n[cls]}$ and $H_t^{n+[cls]}$ are the $[CLS]$ representations of $H_t^n$ and $H_t^{n+}$, $\tau$ is the temperature parameter, and $sim(\cdot)$ indicates the cosine similarity function \citep{chen2020simple}.
    \subsubsection{Optimization}
    The total training loss of an instance is the sum of losses from the slot-value matching for DST and the contrastive context matching for representation, where the former is the average of the losses using the original or the noised context input mentioned in Section
	\ref{base_arc} and \ref{noised_DST}:
	\begin{equation}
	    L_{tot} = (L_{d_{ori}} + L_{d_{nos}})/2 + L_c.
	\end{equation}

	\begin{table*}[!t]
    \centering

    \small
    \resizebox{\textwidth}{!}{
        \scalebox{0.2}{\begin{tabular}{@{}llllllll@{}}
        \toprule
        \multirow{2}{*}{Baseline} & \multirow{2}{*}{Pre-trained Model} & \multicolumn{2}{c}{MultiWOZ 2.0} & \multicolumn{2}{c}{MultiWOZ 2.1} & \multicolumn{2}{c}{MultiWOZ 2.4} \\ \cmidrule(l){3-4} \cmidrule(l){5-6} \cmidrule(l){7-8}

        &  & Joint & Slot & Joint & Slot & Joint & Slot \\ \midrule
        TRADE \citep{wu2019transferable}  & - & 48.62 & 96.92 & 45.60 & 96.55 & 55.05 &  97.62 \\
        SUMBT \citep{lee2019sumbt}  & BERT-base & 42.40 & - & 49.01 & 96.76 & 61.86 & 97.90            \\
        PIN  \citep{chen2020parallel}   & - & 52.44 & 97.28      & 48.40 & 97.02      & 58.92 & 98.02        \\
        SOM-DST \citep{kim2020efficient}  & BERT-base & 51.72 & - & 53.01 & - & 66.78 & 98.38 \\
        CSFN-DST \citep{zhu2020efficient} & BERT-base & 52.23 & - & 53.19 & - & - & - \\
        DST-Picklist \citep{zhang2020find}  & BERT-base &  54.39 & - &  53.30 & 97.40 & - & -            \\
        SAVN \citep{wang2020slot} & BERT-base & 54.52 & 97.42 & 54.86 & 97.55 & 60.55 & 98.05 \\
        SST \citep{chen2020schema} & BERT-base & 51.17 & - & 55.23 & - & - & - \\
        SimpleTOD  \citep{hosseini2020simple}  & DistilGPT2 & - & - & 55.26 & - & - &-    \\
        TripPy \citep{heck2020trippy} & BERT-base & - & - & 55.29 & - & - & -\\

        Seq2SeqDU \citep{feng2021sequence}  & BERT-base & - & - & 56.10 & - & - & -  \\
        STAR \citep{ye2021slot}      & BERT-base & 54.53 & - & 56.36 & 97.59 & 73.62 & 98.85   \\
        SDP-DST \citep{lee2021dialogue} & T5-base & - & - & 56.66 & - & - & - \\
        DS-Graph \citep{lin2021knowledge} & GPT2 & 54.86 & 97.47 & - & - & - & - \\
        DSGFNet \citep{feng2022dynamic} & BERT-base & - & - & 56.70 & - & - & - \\
        PPTOD \citep{su-etal-2022-multi} & T5-large & 53.89 & - & 57.45 & - & - & - \\ \midrule

        Baseline   & BERT-base &  54.38 & 97.47  & 55.82  & 97.51 & 73.81 & 98.82  \\
    	MoNET   & BERT-base & \textbf{55.48} ($\uparrow$ 1.10)  &  97.55 & \textbf{57.71}  ($\uparrow$ 1.89) & 97.71      & \textbf{76.02}  ($\uparrow$ 2.21) & 98.99       \\

        \bottomrule[1.3pt]
        \end{tabular}}
    }
    \caption{Joint and slot goal accuracy of our MoNET and several previous methods on three MultiWOZ test sets.}
	\label{jgamain}
    \end{table*}

	\begin{table}[!t]
		\centering
		\small
		\scalebox{0.95}{
		\begin{tabular}{@{}lccc@{}}
			\toprule[1.3pt]
			Model            &NoisedCM  &NoisedST  & Accuracy \\ \midrule
			Baseline w/o state   & $\times$ & $\times$  & 64.94 \\
			Baseline & $\times$ & $\times$ & 73.81 \\ \midrule
			MoNET-ST & $\times$& \checkmark  & 75.54 \\
		    MoNET-CM &\checkmark &$\times$ & 75.76 \\ \midrule
			MoNET & \checkmark & \checkmark  & 76.02 \\ \bottomrule[1.3pt]
		\end{tabular}}
		\caption{Joint goal accuracy on MultiWOZ 2.4 test set of MoNET and four ablated modifications.}
		\label{bttest}
	\end{table}

	\section{Experiment Setting}
	\subsection{Datasets}
	We choose MultiWOZ, 2.0, 2.1, and 2.4 versions as our datasets.
	MultiWOZ 2.0 \citep{budzianowski2018multiwoz} is a standard human-human conversational dialogue corpus with seven domains.
	MultiWOZ 2.1 \citep{eric2020multiwoz} has the same dialogues as the 2.0 version, where some incorrect state labels are re-annotated.
	Both of them are widely used in previous DST approaches, and MultiWOZ 2.4 \citep{ye2021multiwoz} is the latest refined version correcting all the incorrect state labels in validation and test sets. All three datasets contain the same number of dialogues, which are  8438/1000/1000 in train/validation/test sets.
    For the three datasets, we follow the previous work \citep{wu2019transferable} to use five domains (attraction, hotel, restaurant, taxi, train) with 30 domain-slot pairs in experiments, since the dialogues in the remaining domains are not in the validation and test sets.

	\subsection{Evaluation Metrics}
	The joint and slot goal accuracy are used as the evaluation metrics. Joint goal accuracy is the ratio of dialogue turns where the values of all slots are correctly predicted. Slot goal accuracy is the ratio of domain-slot pairs whose values are correctly predicted. Both of them include correctly predicting those inactive slots with the value $none$.

    \subsection{Baselines}
	We compare the performance of our MoNET with several existing methods, i.e., TRADE, SUMBT, PIN, SOM-DST, CSFN-DST, DST-Picklist, SAVN, SST, SimpleTOD, TripPy, STAR, SDP-DST, DS-Graph, DSGFNet, PPTOD shown in Table~\ref{jgamain}, and our base architecture mentioned in Section \ref{base_arc}, denoted as Baseline.

	\subsection{Training Details}
	The BERT-base-uncased model is used as the context, slot and value encoders, with 12 attention layers and a hidden size of 768.
	During training, only the parameters of the context BERT encoder are updated, while the parameters of the slot and value BERT encoders are not.
	The batch size is set to 8. The AdamW optimizer is applied to optimize the model with the learning rate 4e-5 and 1e-4 for the context encoders and the remaining modules, respectively \citep{loshchilov2018decoupled}.
    The temperature parameter $\tau$ is set to 0.1.

    The noise threshold $p$ defined in Section~\ref{nos_thresh} is set to 0.3,

    and its impact on model performance is discussed in Section \ref{5}. All models are trained on a P40 GPU device for 6-8 hours.

	\section{Results and Analysis}
	\label{5}
	\subsection{Main Results}
	Table \ref{jgamain} shows the performance of MoNET and baselines on MultiWOZ 2.0, 2.1 and 2.4.
	As can be observed, MoNET achieves the joint goal accuracy of 55.48\%, 57.71\%, 76.02\% in three datasets, which are impressive results compared with previous methods, and has 1.10\%, 1.89\%, and 2.21\% improvement on the Baseline model, indicating that our proposed noise-enhanced training helps the model make better predictions.

    \subsection{Ablation Study}

	To explore the performance of each part of our model, we remove the previous dialogue state from the context input of the Baseline model,
	i.e., $X_t = [CLS] \oplus M_t \oplus [SEP] \oplus D_t \oplus [SEP]$, denoted as Baseline w/o state;
	besides, the two noise-enhanced methods are removed from the MoNET respectively, denoted as MoNET-CM (contrastive context matching only) and MoNET-ST (noised state tracking only).
	Table \ref{bttest} shows the performance of the full MoNET model and its four modifications in joint goal accuracy on MultiWOZ 2.4.
	Baseline w/o state gets the lowest accuracy, demonstrating that explicitly using the previous dialogue state as input is beneficial to make predictions.
	Besides, both MoNET-CM and MoNET-ST perform better than Baseline, demonstrating functionalities of the noised state tracking in modifying slot-value pairs that need to be \textit{updated} or \textit{corrected} in further turns, and the context matching framework in learning better semantic representations. Moreover, MoNET gets the best performance, demonstrating the effectiveness of integrating the two parts into a unified noised-enhance training strategy.

	\begin{figure}[!t]
		\centering
		\includegraphics[width=3.0in]{fig/turnlevel3}
		\caption{Turn-level joint goal accuracy and accuracy difference between MoNET and Baseline on MultiWOZ 2.4 test set.}
		\label{turnacc}
	\end{figure}

	\subsection{Turn-Level Evaluation}
	Figure \ref{turnacc} shows the turn-level joint goal accuracy of MoNET and Baseline models, and the accuracy difference percentage (difference between two models divided by the accuracy of Baseline) on MultiWOZ 2.4. The state momentum issue becomes more apparent in dialogues with larger turns, since they always contain more slot-value pairs, and any one of the wrong pairs kept unchanged will affect the further prediction accuracy. With the increase of turns, the accuracy of Baseline harshly degrades, while MoNET gets a relatively lower degradation, and the accuracy difference percentage gradually increases, demonstrating its superiority in avoiding accuracy decrease caused by the state momentum issue, especially for dialogues with larger than 6 turns.

	\begin{figure}[!t]
		\centering
	 	\includegraphics[width=3.0in]{fig/nosrat24}
		\caption{Performance on MultiWOZ 2.4 validation set w.r.t the noise threshold of adding noise.}
		\label{nosexp}
	\end{figure}

	\begin{figure}[!t]
		\centering
	    \includegraphics[width=3.0in]{fig/accsimdegnew1}
		\caption{Performance on MultiWOZ 2.4 test set w.r.t. the noise threshold (corresponding to the noised slot-value pair ratio in the dialogue state input). The left one is the joint goal accuracy, and the right one is the active slot-context features similarity.}
		\label{accsimdeg}
	\end{figure}

	\subsection{Noise Threshold Selection for Training}
	To test the impact of different probabilities of adding noise into the context input for training, we vary the noise threshold $p$ from 0 to 0.5 to train our MoNET. The results on MultiWOZ 2.4 validation set are shown in Figure \ref{nosexp}, where MoNET gets the best performance when the probability $p$ is set to 0.3.
	Intuitively, a small $p$ makes the noised context input contain less noised slot-value pairs (hard to learn semantics from noised data);
	conversely, a large $p$ makes the noised context input far from the original context input in representation space (hard to group them closer).
	Both make it difficult to learn effective representations from the noised context input, leading to a lower prediction accuracy. Hence a proper probability of adding noise is important to derive the best performance.

    \begin{figure*}[!t]
		\centering
		\includegraphics[width=16cm]{fig/nosheat2}
		\caption{Attention visualizations of the two dialogue examples mentioned in Table \ref{cases}.}
		\label{nosheat}
	\end{figure*}

		\begin{table}[!t]
		\centering
		\scriptsize
		\begin{tabular}{@{}l@{}}
			\toprule[1.3pt]
			\begin{tabular}[c]{@{}l@{}}Sys: I have many trains available. What day would you like to leave? \\ Usr: We will be leaving on Sunday afternoon. \\
			GT: train-day-Sunday \\
			Baseline: train-day-\textcolor{blue}{Sunday} \qquad
			MoNET: train-day-\textcolor{blue}{Sunday} \\
			\end{tabular} \\ \midrule
			\begin{tabular}[c]{@{}l@{}}Sys: There are still many trains to pick from, can you narrow down a \\ \qquad departure and arrival time frame? \\ Usr: Yes, it should leave after 20:15 and leave on Monday, not Sunday. \\ \qquad Please give me a price, travel time, and arrival time for any train.\\
			GT: train-day-Monday \\
			Baseline: train-day-\textcolor{red}{Sunday} \qquad
			MoNET: train-day-\textcolor{blue}{Monday} (\textit{update}) \\
			 \end{tabular} \\ \midrule[1.1pt] \midrule[1.1pt]
			\begin{tabular}[c]{@{}l@{}}Sys: Great! You are booked at Autumn House for 1 night. Your reference \\ \qquad number is n4tvfkgs. Would you like more information?\\
			Usr: I would also like a taxi between the places if possible. \\
            GT: taxi-destination-Autumn House \\

			Baseline:
			taxi-destination-\textcolor{red}{Gonville and Caius College} \\
			MoNET:
			taxi-destination-\textcolor{red}{Gonville and Caius College} \\
			\end{tabular} \\ \midrule
			\begin{tabular}[c]{@{}l@{}}Sys: I can help with that. When would you like to either leave or arrive?\\ \qquad And do you want the taxi from the college to the hotel or the other \\ \qquad way around? \\
			Usr: I need to go from the college to the hotel, and I want to leave the \\ \qquad college by 20:45, please. \\
			GT:
			taxi-destination-Autumn House \\
			Baseline: taxi-destination-\textcolor{red}{Gonville and Caius College} \\
            MoNET: taxi-destination-\textcolor{blue}{Autumn House} (\textit{correct})

             \end{tabular} \\
			\bottomrule[1.3pt]
		\end{tabular}
		\caption{Predictions of two dialogue examples on MultiWOZ 2.4 separated by the double solid line, corresponding to two state momentum cases. Wrong and correct predicted values are marked in \textcolor{red}{red} and \textcolor{blue}{blue}.}
		\label{cases}
	\end{table}

    \subsection{Anti-noise Probing with Noise Testing}

	In this section, we conduct \textit{noise testing} to explore the impact of anti-noise ability on DST models.
	We first evaluate the DST performances of MoNET and Baseline by introducing different ratios of noise (with $p$ from $0$ to $1$) into the oracle previous dialogue states as the model input.
    Figure \ref{accsimdeg} shows the performance of MoNET and Baseline on MultiWOZ 2.4.
    Both of them get high accuracy when the noise ratio is $0$, which is because we use the oracle previous dialogue state as the model input;
    with the increase of noise ratio, the joint goal accuracy of Baseline gets a sharp decline, while MoNET degrades much more smoothly.
    This verifies the effectiveness of the noise-enhanced training design of our MoNET which achieves better DST performance with noised input than Baseline.
    Furthermore, for each dialogue turn, we also show the $L2$-distance between the original and noised context representations, i.e., the mean pooling of all token representations $H_t$ and $H_t^+$.
    As can be observed, along with the increase of noise ratio, the distance between the two representations of MoNET is much lower than those of Baseline.
    This phenomenon shows that MoNET achieves higher anti-noise ability than Baseline by giving relatively more similar representations to the original and noised contexts, which helps the model keep an acceptable performance even with a large ratio of noise as the input.

	\subsection{Case Study}
	\label{casestudy}
	Table \ref{cases} shows two prediction examples using MoNET and Baseline on MultiWOZ 2.4. In the first one, they correctly predict the slot-value pair ``train-day-\textcolor{blue}{Sunday}'', while only MoNET updates it in the next turn along with the ground truth changing into ``train-day-\textcolor{blue}{Monday}''. In the second one, they make a wrong prediction ``taxi-destination-\textcolor{red}{Gonville and Caius College}'', while Baseline keeps it unchanged in the next turn, and MoNET corrects it. These indicate the effectiveness of our MoNET in alleviating the two state momentum issues, no matter consistently \textit{updating} the state with the label changes or \textit{correcting} the wrong predictions in further turns.

	\subsection{Attention Visualization}
	We further explore the two examples mentioned in Table \ref{cases} by calculating and visualizing the overall attention scores, which is shown in Figure \ref{nosheat}.
    Given a slot, its overall attention score for each token is the weighted sum of the self-attended scores by all the tokens in $X_t$. The weights come from the slot-context attention, and the self-attended scores are the average of attention scores over multiple layers in BERT. As can be observed, Baseline pays more attention to the values in the previously predicted state, and fails to solve the state momentum issues; MoNET pays relatively higher attention to the correct tokens (``monday'' in the first case and ``autumn house'' in the second one), consequently, successfully \textit{updates} \textcolor{blue}{Sunday} into \textcolor{blue}{Monday} and \textit{corrects} \textcolor{red}{Gonville and Caius College} into \textcolor{blue}{Autumn House}.

	\section{Conclusion}
	State momentum issue in DST task is found and defined in this paper.
	To tackle this issue, we propose MoNET, a method equipped with noised DST training and the contrastive context matching framework.
	Extensive experiments on MultiWOZ 2.0/2.1/2.4 verify the effectiveness of MoNET.
	Supplementary studies and analyses demonstrate that MoNET has a stronger anti-noise ability which helps alleviate the state momentum issues.


\section{Introduction}

The Glashow resonance (GR)\cite{1} at $E_{\nu}\simeq6.3$ PeV in the rest frame of electrons is a  probe of high-energy astrophysical antineutrinos $\bar{\nu}_{e}$\cite{2,3}.
The observation of  GR events can provide useful information on the energy spectra, production source, and flavor transition of high-energy astrophysical neutrinos (HANs).
In recent years, several PeV events of HANs have been reported by the IceCube neutrino observatory\cite{4}. In particular, the first GR event was detected at the $2.3\sigma$ level\cite{5}.
Although there is just one event, the source where no antineutrinos are produced may be excluded as the sole origin of HANs.
Furthermore, the detection of the GR event can also set constraints on new physics effects which become significant at the PeV energy-scale.
Many novel effects in the field of HANs, such as  those from pseudo-Dirac neutrinos\cite{6,7,8,9}, neutrino-decay\cite{10,11,12,13,14,15}, neutrino secret interactions\cite{16,17,18},
etc\cite{19,20,21,22}, have been proposed. In this paper, we take into account the impacts of the Lorentz violation (LV)\cite{23,24,25,26,27,28,29,30,31,32} on
the flux of $\bar{\nu}_{e}$ around the resonance energy and show an interesting observation that the GR events may be damped by the LV.

As is known, the $\bar{\nu}_{e}$ flux controlling the GR event rate is affected by neutrino oscillations over cosmological distances.
The effects of the LV on neutrino oscillations are determined by the Hamiltonian $H=H_{0}+H_{LV}$, where the first and the second term are expressed respectively as following\cite{33}:
\begin{equation}
\label{eq:1}
H_{0}=
U\left(
    \begin{array}{ccc}
      0 & 0 & 0 \\
      0 & \frac{\Delta m^{2}_{21}}{2E_{\nu}} & 0 \\
      0 & 0 & \frac{\Delta m^{2}_{31}}{2E_{\nu}} \\
    \end{array}
  \right)U^{+},
\end{equation}
\begin{equation}
\label{eq:2}
H_{LV}=\pm\left(
                   \begin{array}{ccc}
                     0 & a^{T}_{e\mu} & a^{T}_{e\tau} \\
                     (a^{T}_{e\mu})^{*} & 0 & a^{T}_{e\mu\tau} \\
                      (a^{T}_{e\tau})^{*} &  (a^{T}_{\mu\tau})^{*} & 0 \\
                   \end{array}
                 \right)-\frac{4E_{\nu}}{3}\left(
                   \begin{array}{ccc}
                     0 & c^{TT}_{e\mu} & c^{TT}_{e\tau} \\
                     (c^{TT}_{e\mu})^{*} & 0 & c^{TT}_{e\mu\tau} \\
                      (a^{TT}_{e\tau})^{*} &  (a^{TT}_{\mu\tau})^{*} & 0 \\
                   \end{array}
                 \right).
\end{equation}
Here $U$ is the leptonic mixing matrix in vacuum \cite{34} and matter effects are not considered. For antineutrinos, all terms should be replaced by their complex conjugates and the signal of the $a$ matrix is minus. According to the
recent constraints on the LV parameters from neutrino oscillation experiments\cite{33,35}, the energy-scale of the $a$ matrix can take $10^{-24}$GeV at the $95\%$ confidence limit.
Hence, employing the global fit data on the squared mass differences of neutrinos \cite{36}, we can see that $H_{LV}$  becomes dominant in $H$ when the neutrino energy is of order TeV.
Moreover, even if the $a$ ($c$) parameters are as tiny as $10^{-28}$ GeV ($10^{-31}$), $H_{LV}$ is still significant when $E_{\nu}$ reaches PeV.
Accordingly, the $\bar{\nu}_{e}$ flux leading to the GR events is sensitive to the LV effect which cannot be detected by the neutrino oscillation experiments at low energies.

In the light of the realistic situation of neutrino astronomy, to obtain information on LV from the GR events, two troublesome issues should be addressed.
First, although the first promising point-source of HANs, namely the blazar TXS0506+056, was identified at the $3\sigma$ level\cite{37}, the genuine discovery of the origin of HANs is still in lack at present epoch.
Thus, the original flavor composition of the HANs at the source is totally unknown. For the sake of undermining this problem, two typical sources, namely the pion ($\pi^{\pm}$) decay source, the muon ($\mu^{\pm}$) damped
source, are considered  in this paper. Second, the fitted energy spectra of HANs by now have large uncertainties.
Several models, such as the single power-law, the double power-law,  the power-law with cut-off, etc, can coordinate the observed events of HANs\cite{38,39}.
Furthermore, the constraints on the parameters of the spectra are loose. Thus, moderate LV effects on GR events may be concealed by the uncertainties of the energy spectra.
To overcome this issue, we propose the matrix $H_{LV}$ of a special texture which can  bring the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$).
The GR events from the typical sources are damped by the decoupling effect.
Consequently, stringent constraints on the LV Hamiltonian and the production source of HANs may be obtained.

In the following section, the specific impacts of the LV on the flavor transition of HANs, the flux of $\bar{\nu}_{e}$, and the strength of the GR event rate are shown. Finally, conclusions are given.

\section{Glashow resonance impacted by Lorentz violation}

\subsection{The flavor transition of HANs and the texture of  $H_{LV}$}
The coherence of neutrino mass eigenstates  has disappeared after HANs travelled a cosmological distance. Hence, the averaged neutrino oscillation probability is written as
 \begin{equation}
\label{eq:3}
\overline{P}_{\alpha\beta}^{s}=\sum_{i}|U_{\alpha i}|^{2}|U_{\beta i}|^{2},
\end{equation}
where $\alpha, \beta = e, \mu, \tau$, $i = 1, 2, 3$.
This is the so-called  standard flavor transition probability.
When the LV is taken into account, $\overline{P}_{\alpha\beta}^{s}$ is changed into $\overline{P}_{\alpha\beta}$ which reads
\begin{equation}
\label{eq:4}
\overline{P}_{\alpha\beta}=\sum_{i}|U^{N}_{\alpha i}|^{2}|U^{N}_{\beta i}|^{2}.
\end{equation}
The new mixing matrix $U^{N}$ is obtained by the diagonalization of the total Hamiltonian $H=H_{0}+H_{LV}$.
The eigenvalues of $H$ are written as\cite{33}
\begin{equation}
\label{eq:5}
E_{i}=-2\sqrt{Q}\cos(\frac{\theta_{i}}{3})-\frac{a}{3},~~i=1,~2,~3.
\end{equation}
The components are expressed as follows\cite{33}:
\begin{equation}
\label{eq:6}
Q=\frac{a^{2}-3b}{9},
\end{equation}
\begin{equation}
\label{eq:7}
\theta_{1}=\arccos(RQ^{-\frac{3}{2}}),~~\theta_{2}=\theta_{1}+2\pi,~~\theta_{3}=\theta_{1}-2\pi,
\end{equation}
\begin{equation}
\label{eq:8}
a=-Tr(H),~~b=\frac{Tr(H)^{2}-Tr(H^{2})}{2},~~ c=det(H),
\end{equation}
\begin{equation}
\label{eq:9}
R=\frac{2a^{3}-9ab+27c}{54},
\end{equation}
where the notations Tr and det denote the trace and determinant respectively.
Accordingly, the elements of $U^{N}$ are written as \cite{33}
\begin{equation}
\label{eq:10}
U^{N}_{ei}=\frac{B_{i}^{*}C_{i}}{N_{i}},~~U^{N}_{\mu i}=\frac{A_{i}C_{i}}{N_{i}},~~U^{N}_{\tau i}=\frac{A_{i}B_{i}}{N_{i}},
\end{equation}
in which
\begin{equation}
\label{eq:11}
A_{i}=H_{\mu\tau}(H_{ee}-E_{i})-H_{\mu e}H_{e\tau},
\end{equation}
\begin{equation}
\label{eq:12}
B_{i}=H_{\tau e}(H_{\mu\mu}-E_{i})-H_{\tau\mu }H_{\mu e},
\end{equation}
\begin{equation}
\label{eq:13}
C_{i}=H_{\mu e}(H_{\tau\tau}-E_{i})-H_{\mu\tau}H_{\tau e},
\end{equation}
\begin{equation}
\label{eq:14}
N_{i}^{2}=|A_{i}B_{i}|^{2}+|A_{i}C_{i}|^{2}+|B_{i}C_{i}|^{2}.
\end{equation}

Given the analytical expression of $U^{N}$ and the specific LV parameters, we can predict how the flavor ratio at the source is changed into that at Earth.
However, the issue of the flavor transition including the LV effects is more subtle. On the one hand, since the probability $\overline{P}_{\alpha\beta}$ is dependent on both the magnitude
and the texture of the matrix $H_{LV}$, too many patterns of $\overline{P}$ can be obtained from  the LV parameters in the allowed ranges shown in Tab.\ref{tab:1}.
On the other hand, the experimental constraints on the flavor ratio of HANs at the PeV energy-scale are loose at present. In order to give robust observations on the impacts of the LV on the GR events,
we propose a special pattern with the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$), i.e.,
\begin{equation}
\label{eq:15}
\overline{P}\sim\left(
                  \begin{array}{ccc}
                    \frac{1}{2} & 0 & \frac{1}{2} \\
                    0 & 1 & 0 \\
                     \frac{1}{2}  & 0 &  \frac{1}{2}  \\
                  \end{array}
                \right).
\end{equation}

The decoupling pattern can be obtained from the Hamiltonian $H_{LV}$ dominated by the parameter $(H_{LV})_{e\tau}$, namely $|H_{LV}|_{e\tau} \gg |H_{LV}|_{\alpha\beta}$ with $\alpha(\beta)\neq e, \tau$.
According to the expression of $H_{LV}$ shown in Eq.\ref{eq:2}, we know that the texture of $H_{LV}$ is dependent on the neutrino energy $E_{\nu}$.
For the antineutrinos (neutrinos) with $E_{\nu}\geq1\rm TeV$, the decoupling texture can be realised with the LV parameters in the allowed ranges in Tab.\ref{tab:1}.
The specific decoupling energy-scale $E_{d}$ is determined by the  magnitudes of the LV parameters.
Here we consider two representative decoupling energy-scales for antineutrinos (neutrinos). For the LV parameters $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $C^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$,
we have $E_{d1}\sim5\rm Tev$. For the parameters  $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $C^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$, we obtain $E_{d2}\sim1\rm Pev$.
The dependence of $\overline{P}_{\alpha\beta}$ on $E_{\nu}$ is shown in Fig.\ref{fig:1}.
In the first case, the behavior of $\overline{P}_{\alpha\beta}$ is simple. When $E_{\nu}> E_{d1}$, $\overline{P}_{\alpha\beta}$ shows the decoupling pattern. In the energy range [1TeV, 5TeV],
a small variation is added to the  pattern.
In the second case, the energy $E_{\nu}$ can be classified into 3 areas.
In the range [1TeV, 10TeV], $H_{0}$ is dominant in $H$, namely  $\overline{P}_{\alpha\beta} \sim\overline{P^{s}}_{\alpha\beta}$.
In the range [10TeV, 1PeV], the impact of $H_{LV}$ is noticeable. The variation of  $\overline{P}_{\alpha\beta}$ with $E_{\nu}$ is steep, especially in the range [100TeV, 1PeV].
In the range $E_{\nu}>E_{d2}$, $\overline{P}_{\alpha\beta}$ is also converged to the decoupling pattern.
General speaking, we can see that tiny LV parameters may bring an interesting phenomenology of the flavor transition of HANs, which hence impacts the flux of $\bar{\nu}_{e}$.

\begin{table}
\caption{\label{tab:1} Constraints on the Lorentz violation parameters \cite{33}.}
  \centering
  \begin{tabular}{c c c}
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     LV Parameters ~~&~~ $95\%$ confidence limit  ~~&~~ Best fit \\
     \noalign{\smallskip}\hline
     Re$(a^{T}_{e\mu})$ ~~~~&~~~~ $1.8\times10^{-23}$ GeV ~~~~&~~~~ $1.0\times10^{-23}$ GeV \\
    Im$(a^{T}_{e\mu})$  ~~~~&~~~~$1.8\times10^{-23}$ GeV ~~~~&~~~~ $4.6\times10^{-24}$ GeV \\
    Re$(c^{TT}_{e\mu})$ ~~~~&~~~~ $8.0\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$ \\
   Im$(c^{TT}_{e\mu})$   ~~~~&~~~~$8.0\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$  \\
   \noalign{\smallskip}\hline
      Re$(a^{T}_{e\tau})$  ~~~~&~~~~$4.1\times10^{-23}$ GeV ~~~~&~~~~ $2.2\times10^{-24}$ GeV\\
     Im$(a^{T}_{e\tau})$~~~~&~~~~ $2.8\times10^{-23}$ GeV ~~~~&~~~~ $1.0\times10^{-28}$ GeV \\
     Re$(c^{TT}_{e\tau})$~~~~&~~~~$9.3\times10^{-25}$ ~~~~&~~~~$1.0\times10^{-28}$  \\
    Im$(c^{TT}_{e\tau})$  ~~~~&~~~~ $1.0\times10^{-24}$ ~~~~&~~~~$3.5\times10^{-25}$  \\
    \noalign{\smallskip}\hline
       Re$(a^{T}_{\mu\tau})$~~~~&~~~~ $6.5\times10^{-24}$ GeV ~~~~&~~~~ $3.2\times10^{-24}$ GeV \\
      Im$(a^{T}_{\mu\tau})$ ~~~~&~~~~$5.1\times10^{-24}$ GeV ~~~~&~~~~ $1.0\times10^{-28}$ GeV \\
       Re$(c^{TT}_{\mu\tau})$ ~~~~&~~~~ ~$4.4\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$ \\
     Im$(c^{TT}_{\mu\tau})$   ~~~~&~~~~ ~$4.2\times10^{-27}$ ~~~~&~~~~$7.5\times10^{-28}$ \\
     \hline
   \end{tabular}
\end{table}

\begin{figure}
  \centering

  \includegraphics[width=.49\textwidth]{fig1a}
  \hfill
  \includegraphics[width=.49\textwidth]{fig1b}
  \caption{\label{fig:1} The flavor transition probability of antineutrinos.
$\overline{P}_{\alpha\beta}$ is obtained from the leptonic mixing parameters at the $3\sigma$ level of the global fit data with normal
  mass ordering(NO)\cite{36} and the LV parameters.
 For blue lines, $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$. For red lines, $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$.  $(a^{T}_{\alpha\beta})_{bf}$, $(c^{TT}_{\alpha\beta})_{bf}$  are the  best fit values listed in Tab.\ref{tab:1}.
 Black line: $\overline{P^{s}}_{\alpha\beta}$ from the best fit values of the global fit data with NO\cite{36}.
 }
\end{figure}

\subsection{The energy spectrum of $\bar{\nu}_{e}$ modified by the LV}
As has been shown, the LV resulting in the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling  has significant impacts on the flavor transition of HANs.  Thus, the energy spectrum of $\bar{\nu}(e)$
may be modified by the LV effect. In addition to $H_{LV}$, the flux of $\bar{\nu}_{e}$ at Earth is also dependent on the production source of HANS.
Since the source of $\bar{\nu}_{e}$ is uncertain at present, here we consider two typical sources to set the original flavor composition.
They are listed as follows: the pion ($\pi^{\pm}$)  decay source with the original flavor ratio
\begin{equation}
\label{eq:16}
R^{s(\pi^{\pm})}=( 1/3, ~2/3,~ 0 ),
\end{equation}
and the muon ($\mu^{\pm}$) damped source with the ratio
\begin{equation}
\label{eq:17}
R^{s(\mu^{\pm})}=( 0,~ 1,~ 0 ).
\end{equation}
For the both sources, the flux of neutrinos is equal to that of antineutrinos.
Employing the flavor ratio at the source and the transition matrix $\overline{P}$, we can derive the ratio of $\bar{\nu}_{e}$ at Earth, namely
\begin{equation}
\label{eq:18}
r_{e}(E_{\nu})=\frac{\phi_{\bar{\nu}_{e}}(E_{\nu})}{\phi_{\nu+\bar{\nu}}(E_{\nu})}=\frac{1}{2}\sum_{\alpha=e,\mu,\tau}\overline{P}_{e\alpha}(E_{\nu})R^{s}_{\alpha},
\end{equation}
in which $\phi_{\bar{\nu}_{e}}$ and $\phi_{\nu+\bar{\nu}}$ are the  differentia flux  of $\bar{\nu}_{e}$ and that of total HANs respectively,
$R^{s}_{\alpha}$ is the component of the flavor ratio vector at the source.
Accordingly, the differentia flux of $\bar{\nu}_{e}$ is predicted to be
\begin{equation}
\label{eq:19}
\phi_{\bar{\nu}_{e}}=\frac{1}{2}\phi_{\nu+\bar{\nu}}(E_{\nu})\times(\sum_{\alpha=e,\mu,\tau}\overline{P}_{e\alpha}(E_{\nu})R^{s}_{\alpha}~).
\end{equation}

At present, the uncertainty of $\phi_{\nu+\bar{\nu}}$ is large. Several models of the energy spectra  can work, e.g., the single power-law,
double power-law, single power law with spectral cutoff, log-parabola, segmented power-law.
For the sake of illustration, we consider the single power-law and the double power-law spectrum in this paper.
They are expressed respectively as\cite{38}
\begin{equation}
\label{eq:20}
\phi_{\nu+\bar{\nu}}^{spl}=\phi_{astro}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{astro}}\cdot10^{-18}\rm GeV^{-1}cm^{-2}s^{-1}sr^{-1},
\end{equation}
\begin{equation}
\label{eq:21}
\phi_{\nu+\bar{\nu}}^{dpl}=(\phi_{hard}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{hard}}+\phi_{soft}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{soft}})\cdot10^{-18}\rm GeV^{-1}cm^{-2}s^{-1}sr^{-1}.
\end{equation}
The ranges of the normalization and the spectral index at the 68.3\% confidence level are listed in Tab.\ref{tab:2}.
\begin{table}
\caption{\label{tab:2} Energy-spectrum parameters of HANs at the 68.3\% confidence level\cite{38}. }
  \centering
  \begin{tabular}{c c c c c }
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     Model~~&~~~~~~ normalization  ~~&~~~~&~ ~ spectral index  ~&~~ \\
     \noalign{\smallskip}\hline
     single power-law ~~~& $\phi_{astro}$ ~~~&~&~~~ $\gamma_{astro}$ ~ \\
  ~~~~~&$5.68^{+1.56}_{-1.55}$ ~~~&~&~~~~ $2.89^{+0.23}_{-0.20}$ \\
   \noalign{\smallskip}\hline
       double power-law ~~~&$\phi_{hard}$ ~~&~ $\phi_{soft}$ ~~&~~ $\gamma_{hard}$ &~$\gamma_{soft}$~ \\
    ~~~~&$3.26^{+1.88}_{-2.44}$ ~~&~ $0.08^{+3.64}_{-0.08}$~&~~~~~ $2.78^{+0.24}_{-0.32}$~~&~~~~ $3.12^{+1.01}_{-0.30}$\\
     \hline
   \end{tabular}
\end{table}
Using these parameters, we show the LV effect on the energy spectrum of $\bar{\nu}_{e}$ in Fig.\ref{fig:2} and Fig.\ref{fig:3}.

For the  muon ($\mu^{\pm}$) damped source (see Fig.\ref{fig:3}), the LV can bring a noticeable damping of the flux of $\bar{\nu}_{e}$ when $E_{\nu}\geq E_{d}$.
The damping rate $\phi_{\bar{\nu}_{e}}^{LV}$/ $\phi_{\bar{\nu}_{e}}^{S }$ is of order [$10^{-6}$, $10^{-4}$] at the resonance energy. Here $\phi_{\bar{\nu}_{e}}^{LV}$ and $\phi_{\bar{\nu}_{e}}^{S }$ denote the flux including the LV and that without the LV respectively. The specific  damping value is dependent on the form of $\phi_{\nu+\bar{\nu}}$.
For the  pion ($\pi^{\pm}$) decay source (see Fig.\ref{fig:2}), the damping effect is moderate and concealed by the large uncertainty of $\phi_{\nu+\bar{\nu}}$, irrespective of the magnitude of the decoupling energy.
In general, we can expect that  the LV effect can lead to a noticeable decrease of  the $\bar{\nu}_{e}$ flux around the resonance energy for a source dominated by $\bar{\nu}_{\mu}$ ($\nu_{\mu}$).
Thus, the detection of GR events may set stringent constraints on the texture of $H_{LV}$  and the original flavor ratio.

\begin{figure}
\label{fig:2}
  \centering

  \includegraphics[width=.49\textwidth]{fig2a}
   \hfill
  \includegraphics[width=.49\textwidth]{fig2b}
   \includegraphics[width=.49\textwidth]{fig2c}
   \hfill
  \includegraphics[width=.49\textwidth]{fig2d}
  \caption{\label{fig:2} The energy spectrum of $\bar{\nu}_{e}$ from the pion ($\pi^{\pm}$) decay source. Left plots: from the total flux of the single power-law. Right plots: from the total flux of the double power-law.
   The ranges of the parameters of the energy spectrum $\phi_{\nu+\bar{\nu}}$ are taken as those listed in Tab.\ref{tab:2}.
   Blue lines: arising from $\overline{P^{s}}$ without the LV effects.  Red lines:  from $\overline{P}$ with  the LV effects. For $\overline{P^{s}}$ and $\overline{P}$,
   the leptonic mixing parameters are in the $3\sigma$ allowed ranges of the global fit
   data with NO\cite{36}. The plots in the first row: with the LV parameters  taken as $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$.
   The plots in the second row: $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$.
  The black line: arising from $\overline{P}$ with the best fit values of the energy spectrum $\phi_{\nu+\bar{\nu}}$ and leptonic mixing parameters. The dashed black line: from $\overline{P^{s}}$ with the best fit values of the leptonic mixing parameters.}
\end{figure}
\begin{figure}
\label{fig:3}
  \centering

  \includegraphics[width=.49\textwidth]{fig3a}
   \hfill
  \includegraphics[width=.49\textwidth]{fig3b}
  \includegraphics[width=.49\textwidth]{fig3c}
   \hfill
  \includegraphics[width=.49\textwidth]{fig3d}
  \caption{\label{fig:3} The energy spectrum of $\bar{\nu}_{e}$ from the muon ($\mu^{\pm}$) damped source. The conventions of the parameters, colors and the lines are the same as those in Fig.\ref{fig:2}.}
\end{figure}

\subsection{The strength of the GR event rate damped by the LV}
The GR events provide a window to explore the texture of $H_{LV}$ and the flavor composition of the HANs at the PeV energy-scale.
Now let us examine the impacts of the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) on the GR event rate.
To weaken the influence of the uncertain energy spectrum  $\phi_{\nu+\bar{\nu}}$,  for down-coming events, the ratio of the resonance events to the nonresonant continuum events is proposed
to describe the strength of the GR event rate, i.e.\cite{2},
\begin{equation}
\label{eq:22}
\begin{aligned}
  \frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})}&=\frac{10\pi}{18}(\frac{\Gamma_{W}}{M_{W}})(\frac{\sigma^{peak}_{Res}}{\sigma^{CC}_{\nu N}(E_{\nu}=6.3\rm PeV)})\frac{(\alpha-1.4)(\frac{E_{\nu}^{min}}{6.3\rm PeV})^{\alpha-1.4}}{[1-(\frac{E_{\nu}^{min}}{E_{\nu}^{max}})^{(\alpha-1.4)}]}\times [r_{e}]_{E_{\nu}=6.3\rm PeV} \\
  &=11\times\frac{(\alpha-1.4)(\frac{E_{\nu}^{min}}{6.3\rm PeV})^{\alpha-1.4}}{[1-(\frac{E_{\nu}^{min}}{E_{\nu}^{max}})^{(\alpha-1.4)}]}\times [r_{e}]_{E_{\nu}=6.3\rm PeV}.
\end{aligned}
\end{equation}
Here the spectral index $\alpha$ is dependent on the acceleration mechanism of neutrinos. In absence of GR events, $\alpha\geq2.3$\cite{40,41}. In this paper, we take $\alpha=2$.
We note that a moderate variation of $\alpha$ will not lead to a noticeable change of the observation on the LV effect.
Employing the best fit values of the leptonic mixing parameters with NO\cite{36} and the LV parameters $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$ with $(a^{T}_{\alpha\beta})_{bf}$, $(c^{TT}_{\alpha\beta})_{bf}$ listed in Tab.\ref{tab:1}, we calculate the ratio of resonance event rate around the resonance energy.
The specific values are listed in Tab.\ref{tab:3}. For the sake of comparison, the parenthetic values obtained from the standard flavor transition matrix $\overline{P^{s}}$ are also shown.
\begin{table}
\caption{\label{tab:3} The ratio of resonance event rate to nonresonant event rate with and without the LV effect for the  pion ($\pi^{\pm}$) decay source, muon ($\mu^{\pm}$) damped source, and a source rich in $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) with $\phi_{\nu}=\phi_{\bar{\nu}}$. Here $\alpha=2$, $E_{\nu}^{min}$=1, 2, 3, 4, 5 PeV, $E_{\nu}^{max}=\infty$, $r_{e}$ takes the value $ [r_{e}]_{E_{\nu}=6.3\rm PeV}$.
The parenthetic values are obtained from the standard flavor transition matrix $\overline{P^{s}}$. }
  \centering
  \begin{tabular}{c c c c c c c}
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     $E_{\nu}^{min}(\rm PeV)$ ~&~~1  ~&~2~~&~~3~~&~~ 4~~&~~ 5&~~ \\
     \noalign{\smallskip}\hline
$R^{s(\pi^{\pm})}=( 1/3, ~2/3,~ 0 )$ ~&~~0.18~~&~~~0.28~ &~~~0.35 &~~~0.42&~~~0.48&~$r_{e}$=0.083\\
~~&~~(0.32)~~&~~~(0.50)~ &~~~(0.63) &~~~(0.75)&~~~(0.86)&~($r_{e}$=0.149)\\ \hline
$R^{s(\mu^{\pm})}=( 0,~ 1,~ 0 )$  ~~~~~&~~$7.8\times10^{-6}$~~~~&~~$1.2\times10^{-5}$~~&~$1.5\times10^{-5}$&~$1.8\times10^{-5}$ &~$2.0\times10^{-5}$&~$r_{e}=3.5\times10^{-6}$\\
~~~~&~~~~(0.19)~~~~&~~~(0.28)~~&~(0.36)&~(0.43) &~(0.49)&~($r_{e}$=0.086)\\ \hline
$R^{s(x)}=( 0.1,~ 0.9,~ 0 )$ ~~&~~0.05~~&~~~0.08~~ &~~0.10~~&~~0.12~~&~~0.14~~&~$r_{e}$=0.025\\
~~~~&~~~~~(0.23)~~~~&~~~(0.35)~~&~~(0.44)&~~(0.53) &~~(0.60)&~~($r_{e}$=0.105)\\
     \hline
   \end{tabular}
\end{table}

From Tab.\ref{tab:3}, we can see that the ratio of  resonance event rate to the nonresonant event rate
is mainly determined by the ratio of $\bar{\nu}_{e}$ at the resonance energy-scale. The damping rate of the GR event rate approximates that of $ [r_{e}]_{E_{\nu}=6.3\rm PeV}$,
namely
\begin{equation}
\label{eq:23}
  (\frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})})^{LV}/(\frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})})^{S}\sim ([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}/([r_{e}]_{E_{\nu}=6.3\rm PeV})^{S},
\end{equation}
in which the superscript 'LV' ('S') denotes the value with (without) the LV effect.
For the original flavor ratios $R^{s(\pi^{\pm})}$, $R^{s(\mu^{\pm})}$, and $R^{s(x)}$, $([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}/([r_{e}]_{E_{\nu}=6.3\rm PeV})^{S}$ takes 0.557, $4.07\times10^{-5}$, 0.238, respectively.
The damping rate is sensitive to the increase of the original ratio of $\bar{\nu}_{\mu}$ around the resonance energy when $\phi_{\bar{\nu}_{\mu}}/\phi_{\bar{\nu}}>0.9$.
Furthermore, when the decoupling energy is around $1\rm PeV$, i.e., $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$,
the values of  $([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}$ are similar to those listed in Tab.\ref{tab:3}. Thus, the damping rate of the GR event rate is mainly
sensitive to the texture of $H_{LV}$, irrespective the magnitudes of the LV parameters.
However, we should keep in mind  that the magnitudes of the  parameters determine the decoupling energy which impacts the evolution of $\overline{P}$ with $E_{\nu}$.
In the range $E_{\nu}<E_{d}$, the flavor ratio of HANs at Earth is dependent on the magnitudes of the LV parameters.

\subsection{Comments on the realistic source and the $\nu_{\tau}$ event at the PeV energy-scale }

At present, the production source and the acceleration mechanism of HANs are unknown. However, it is widely believed that the strong magnetic field plays an important role in the
acceleration of HANs\cite{42}, e.g. acceleration by the magnetic reconnection\cite{43,44,45}. As is known, the synchrotron radiation of the charged particle in the strong magnetic field becomes significant at high energies. \
We can expect that the muon ($\mu^{\pm}$) damped source may be dominant at the PeV energy-scale.  Accordingly, when $E_{d}<6.3\rm PeV$, the decoupling pattern can show a noticeable impact on the GR events.
Furthermore, we note that the decoupling pattern can bring the following observation:
\begin{equation}
\label{eq:24}
\phi_{\nu_{e}}\sim\phi_{\nu_{\tau}},~~~~\phi_{\bar{\nu}_{e}}\sim\phi_{\bar{\nu}_{\tau}}.
\end{equation}
Hence, for the muon ($\mu^{\pm}$) damped source, the proposed LV texture can lead to a notable decrease of both the GR event rate and the double bang events from $\nu_{\tau}$ at the PeV energy-scale.
In other words, we can obtain the double constraints on the texture of $H_{LV}$ from the GR event and the double bang event.

\section{Conclusions }

The Lorentz invariance is a fundamental symmetry in nature. Strong constraints on the LV parameters have been obtained at low energies.
To connect the LV to the GR event at the PeV energy-scale, we proposed the LV Hamiltonian of the special texture which can  bring the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling  at high energies.
On the base of the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling, we analysed the impacts of the LV parameters on the flavor transition of HANs.
For the typical sources of HANs, the damping effect from the decoupling pattern on the energy spectrum of $\bar{\nu}_{e}$ was shown.
For the  muon ($\mu^{\pm}$) damped source, the decrease of the $\bar{\nu}_{e}$ flux at the PeV energy-scale is notable, irrespective of the form of the total flux $\phi_{\nu+\bar{\nu}}$.
The LV effect can lead to a noticeable damping of the GR event rate if the HANs are mainly from the source dominated by $\bar{\nu}_{\mu}$ ($\nu_{\mu}$).
Furthermore, the decoupling pattern can also decrease the double bang events from $\nu_{\tau}$ at the PeV energy-scale.
Therefore, the detection of GR events and other high-energy events may bring extra constraints on the LV Hamiltonian and the production source of HANs.

\acknowledgments
 This work is supported by the National Natural Science Foundation of China under grant No. 12065007, the Guangxi Scientific Programm Foundation under grant No. Guike AD19110045.\\

\section{Introduction}

The Glashow resonance (GR)\cite{1} at $E_{\nu}\simeq6.3$ PeV in the rest frame of electrons is a  probe of high-energy astrophysical antineutrinos $\bar{\nu}_{e}$\cite{2,3}.
The observation of  GR events can provide useful information on the energy spectra, production source, and flavor transition of high-energy astrophysical neutrinos (HANs).
In recent years, several PeV events of HANs have been reported by the IceCube neutrino observatory\cite{4}. In particular, the first GR event was detected at the $2.3\sigma$ level\cite{5}.
Although there is just one event, the source where no antineutrinos are produced may be excluded as the sole origin of HANs.
Furthermore, the detection of the GR event can also set constraints on new physics effects which become significant at the PeV energy-scale.
Many novel effects in the field of HANs, such as  those from pseudo-Dirac neutrinos\cite{6,7,8,9}, neutrino-decay\cite{10,11,12,13,14,15}, neutrino secret interactions\cite{16,17,18},
etc\cite{19,20,21,22}, have been proposed. In this paper, we take into account the impacts of the Lorentz violation (LV)\cite{23,24,25,26,27,28,29,30,31,32} on
the flux of $\bar{\nu}_{e}$ around the resonance energy and show an interesting observation that the GR events may be damped by the LV.

As is known, the $\bar{\nu}_{e}$ flux controlling the GR event rate is affected by neutrino oscillations over cosmological distances.
The effects of the LV on neutrino oscillations are determined by the Hamiltonian $H=H_{0}+H_{LV}$, where the first and the second term are expressed respectively as following\cite{33}:
\begin{equation}
\label{eq:1}
H_{0}=
U\left(
    \begin{array}{ccc}
      0 & 0 & 0 \\
      0 & \frac{\Delta m^{2}_{21}}{2E_{\nu}} & 0 \\
      0 & 0 & \frac{\Delta m^{2}_{31}}{2E_{\nu}} \\
    \end{array}
  \right)U^{+},
\end{equation}
\begin{equation}
\label{eq:2}
H_{LV}=\pm\left(
                   \begin{array}{ccc}
                     0 & a^{T}_{e\mu} & a^{T}_{e\tau} \\
                     (a^{T}_{e\mu})^{*} & 0 & a^{T}_{e\mu\tau} \\
                      (a^{T}_{e\tau})^{*} &  (a^{T}_{\mu\tau})^{*} & 0 \\
                   \end{array}
                 \right)-\frac{4E_{\nu}}{3}\left(
                   \begin{array}{ccc}
                     0 & c^{TT}_{e\mu} & c^{TT}_{e\tau} \\
                     (c^{TT}_{e\mu})^{*} & 0 & c^{TT}_{e\mu\tau} \\
                      (a^{TT}_{e\tau})^{*} &  (a^{TT}_{\mu\tau})^{*} & 0 \\
                   \end{array}
                 \right).
\end{equation}
Here $U$ is the leptonic mixing matrix in vacuum \cite{34} and matter effects are not considered. For antineutrinos, all terms should be replaced by their complex conjugates and the signal of the $a$ matrix is minus. According to the
recent constraints on the LV parameters from neutrino oscillation experiments\cite{33,35}, the energy-scale of the $a$ matrix can take $10^{-24}$GeV at the $95\%$ confidence limit.
Hence, employing the global fit data on the squared mass differences of neutrinos \cite{36}, we can see that $H_{LV}$  becomes dominant in $H$ when the neutrino energy is of order TeV.
Moreover, even if the $a$ ($c$) parameters are as tiny as $10^{-28}$ GeV ($10^{-31}$), $H_{LV}$ is still significant when $E_{\nu}$ reaches PeV.
Accordingly, the $\bar{\nu}_{e}$ flux leading to the GR events is sensitive to the LV effect which cannot be detected by the neutrino oscillation experiments at low energies.

In the light of the realistic situation of neutrino astronomy, to obtain information on LV from the GR events, two troublesome issues should be addressed.
First, although the first promising point-source of HANs, namely the blazar TXS0506+056, was identified at the $3\sigma$ level\cite{37}, the genuine discovery of the origin of HANs is still in lack at present epoch.
Thus, the original flavor composition of the HANs at the source is totally unknown. For the sake of undermining this problem, two typical sources, namely the pion ($\pi^{\pm}$) decay source, the muon ($\mu^{\pm}$) damped
source, are considered  in this paper. Second, the fitted energy spectra of HANs by now have large uncertainties.
Several models, such as the single power-law, the double power-law,  the power-law with cut-off, etc, can coordinate the observed events of HANs\cite{38,39}.
Furthermore, the constraints on the parameters of the spectra are loose. Thus, moderate LV effects on GR events may be concealed by the uncertainties of the energy spectra.
To overcome this issue, we propose the matrix $H_{LV}$ of a special texture which can  bring the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$).
The GR events from the typical sources are damped by the decoupling effect.
Consequently, stringent constraints on the LV Hamiltonian and the production source of HANs may be obtained.

In the following section, the specific impacts of the LV on the flavor transition of HANs, the flux of $\bar{\nu}_{e}$, and the strength of the GR event rate are shown. Finally, conclusions are given.

\section{Glashow resonance impacted by Lorentz violation}

\subsection{The flavor transition of HANs and the texture of  $H_{LV}$}
The coherence of neutrino mass eigenstates  has disappeared after HANs travelled a cosmological distance. Hence, the averaged neutrino oscillation probability is written as
 \begin{equation}
\label{eq:3}
\overline{P}_{\alpha\beta}^{s}=\sum_{i}|U_{\alpha i}|^{2}|U_{\beta i}|^{2},
\end{equation}
where $\alpha, \beta = e, \mu, \tau$, $i = 1, 2, 3$.
This is the so-called  standard flavor transition probability.
When the LV is taken into account, $\overline{P}_{\alpha\beta}^{s}$ is changed into $\overline{P}_{\alpha\beta}$ which reads
\begin{equation}
\label{eq:4}
\overline{P}_{\alpha\beta}=\sum_{i}|U^{N}_{\alpha i}|^{2}|U^{N}_{\beta i}|^{2}.
\end{equation}
The new mixing matrix $U^{N}$ is obtained by the diagonalization of the total Hamiltonian $H=H_{0}+H_{LV}$.
The eigenvalues of $H$ are written as\cite{33}
\begin{equation}
\label{eq:5}
E_{i}=-2\sqrt{Q}\cos(\frac{\theta_{i}}{3})-\frac{a}{3},~~i=1,~2,~3.
\end{equation}
The components are expressed as follows\cite{33}:
\begin{equation}
\label{eq:6}
Q=\frac{a^{2}-3b}{9},
\end{equation}
\begin{equation}
\label{eq:7}
\theta_{1}=\arccos(RQ^{-\frac{3}{2}}),~~\theta_{2}=\theta_{1}+2\pi,~~\theta_{3}=\theta_{1}-2\pi,
\end{equation}
\begin{equation}
\label{eq:8}
a=-Tr(H),~~b=\frac{Tr(H)^{2}-Tr(H^{2})}{2},~~ c=det(H),
\end{equation}
\begin{equation}
\label{eq:9}
R=\frac{2a^{3}-9ab+27c}{54},
\end{equation}
where the notations Tr and det denote the trace and determinant respectively.
Accordingly, the elements of $U^{N}$ are written as \cite{33}
\begin{equation}
\label{eq:10}
U^{N}_{ei}=\frac{B_{i}^{*}C_{i}}{N_{i}},~~U^{N}_{\mu i}=\frac{A_{i}C_{i}}{N_{i}},~~U^{N}_{\tau i}=\frac{A_{i}B_{i}}{N_{i}},
\end{equation}
in which
\begin{equation}
\label{eq:11}
A_{i}=H_{\mu\tau}(H_{ee}-E_{i})-H_{\mu e}H_{e\tau},
\end{equation}
\begin{equation}
\label{eq:12}
B_{i}=H_{\tau e}(H_{\mu\mu}-E_{i})-H_{\tau\mu }H_{\mu e},
\end{equation}
\begin{equation}
\label{eq:13}
C_{i}=H_{\mu e}(H_{\tau\tau}-E_{i})-H_{\mu\tau}H_{\tau e},
\end{equation}
\begin{equation}
\label{eq:14}
N_{i}^{2}=|A_{i}B_{i}|^{2}+|A_{i}C_{i}|^{2}+|B_{i}C_{i}|^{2}.
\end{equation}

Given the analytical expression of $U^{N}$ and the specific LV parameters, we can predict how the flavor ratio at the source is changed into that at Earth.
However, the issue of the flavor transition including the LV effects is more subtle. On the one hand, since the probability $\overline{P}_{\alpha\beta}$ is dependent on both the magnitude
and the texture of the matrix $H_{LV}$, too many patterns of $\overline{P}$ can be obtained from  the LV parameters in the allowed ranges shown in Tab.\ref{tab:1}.
On the other hand, the experimental constraints on the flavor ratio of HANs at the PeV energy-scale are loose at present. In order to give robust observations on the impacts of the LV on the GR events,
we propose a special pattern with the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$), i.e.,
\begin{equation}
\label{eq:15}
\overline{P}\sim\left(
                  \begin{array}{ccc}
                    \frac{1}{2} & 0 & \frac{1}{2} \\
                    0 & 1 & 0 \\
                     \frac{1}{2}  & 0 &  \frac{1}{2}  \\
                  \end{array}
                \right).
\end{equation}

The decoupling pattern can be obtained from the Hamiltonian $H_{LV}$ dominated by the parameter $(H_{LV})_{e\tau}$, namely $|H_{LV}|_{e\tau} \gg |H_{LV}|_{\alpha\beta}$ with $\alpha(\beta)\neq e, \tau$.
According to the expression of $H_{LV}$ shown in Eq.\ref{eq:2}, we know that the texture of $H_{LV}$ is dependent on the neutrino energy $E_{\nu}$.
For the antineutrinos (neutrinos) with $E_{\nu}\geq1\rm TeV$, the decoupling texture can be realised with the LV parameters in the allowed ranges in Tab.\ref{tab:1}.
The specific decoupling energy-scale $E_{d}$ is determined by the  magnitudes of the LV parameters.
Here we consider two representative decoupling energy-scales for antineutrinos (neutrinos). For the LV parameters $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $C^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$,
we have $E_{d1}\sim5\rm Tev$. For the parameters  $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $C^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$, we obtain $E_{d2}\sim1\rm Pev$.
The dependence of $\overline{P}_{\alpha\beta}$ on $E_{\nu}$ is shown in Fig.\ref{fig:1}.
In the first case, the behavior of $\overline{P}_{\alpha\beta}$ is simple. When $E_{\nu}> E_{d1}$, $\overline{P}_{\alpha\beta}$ shows the decoupling pattern. In the energy range [1TeV, 5TeV],
a small variation is added to the  pattern.
In the second case, the energy $E_{\nu}$ can be classified into 3 areas.
In the range [1TeV, 10TeV], $H_{0}$ is dominant in $H$, namely  $\overline{P}_{\alpha\beta} \sim\overline{P^{s}}_{\alpha\beta}$.
In the range [10TeV, 1PeV], the impact of $H_{LV}$ is noticeable. The variation of  $\overline{P}_{\alpha\beta}$ with $E_{\nu}$ is steep, especially in the range [100TeV, 1PeV].
In the range $E_{\nu}>E_{d2}$, $\overline{P}_{\alpha\beta}$ is also converged to the decoupling pattern.
General speaking, we can see that tiny LV parameters may bring an interesting phenomenology of the flavor transition of HANs, which hence impacts the flux of $\bar{\nu}_{e}$.

\begin{table}
\caption{\label{tab:1} Constraints on the Lorentz violation parameters \cite{33}.}
  \centering
  \begin{tabular}{c c c}
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     LV Parameters ~~&~~ $95\%$ confidence limit  ~~&~~ Best fit \\
     \noalign{\smallskip}\hline
     Re$(a^{T}_{e\mu})$ ~~~~&~~~~ $1.8\times10^{-23}$ GeV ~~~~&~~~~ $1.0\times10^{-23}$ GeV \\
    Im$(a^{T}_{e\mu})$  ~~~~&~~~~$1.8\times10^{-23}$ GeV ~~~~&~~~~ $4.6\times10^{-24}$ GeV \\
    Re$(c^{TT}_{e\mu})$ ~~~~&~~~~ $8.0\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$ \\
   Im$(c^{TT}_{e\mu})$   ~~~~&~~~~$8.0\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$  \\
   \noalign{\smallskip}\hline
      Re$(a^{T}_{e\tau})$  ~~~~&~~~~$4.1\times10^{-23}$ GeV ~~~~&~~~~ $2.2\times10^{-24}$ GeV\\
     Im$(a^{T}_{e\tau})$~~~~&~~~~ $2.8\times10^{-23}$ GeV ~~~~&~~~~ $1.0\times10^{-28}$ GeV \\
     Re$(c^{TT}_{e\tau})$~~~~&~~~~$9.3\times10^{-25}$ ~~~~&~~~~$1.0\times10^{-28}$  \\
    Im$(c^{TT}_{e\tau})$  ~~~~&~~~~ $1.0\times10^{-24}$ ~~~~&~~~~$3.5\times10^{-25}$  \\
    \noalign{\smallskip}\hline
       Re$(a^{T}_{\mu\tau})$~~~~&~~~~ $6.5\times10^{-24}$ GeV ~~~~&~~~~ $3.2\times10^{-24}$ GeV \\
      Im$(a^{T}_{\mu\tau})$ ~~~~&~~~~$5.1\times10^{-24}$ GeV ~~~~&~~~~ $1.0\times10^{-28}$ GeV \\
       Re$(c^{TT}_{\mu\tau})$ ~~~~&~~~~ ~$4.4\times10^{-27}$ ~~~~&~~~~$1.0\times10^{-28}$ \\
     Im$(c^{TT}_{\mu\tau})$   ~~~~&~~~~ ~$4.2\times10^{-27}$ ~~~~&~~~~$7.5\times10^{-28}$ \\
     \hline
   \end{tabular}
\end{table}

\begin{figure}
  \centering

  \includegraphics[width=.49\textwidth]{fig1a}
  \hfill
  \includegraphics[width=.49\textwidth]{fig1b}
  \caption{\label{fig:1} The flavor transition probability of antineutrinos.
$\overline{P}_{\alpha\beta}$ is obtained from the leptonic mixing parameters at the $3\sigma$ level of the global fit data with normal
  mass ordering(NO)\cite{36} and the LV parameters.
 For blue lines, $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$. For red lines, $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$.  $(a^{T}_{\alpha\beta})_{bf}$, $(c^{TT}_{\alpha\beta})_{bf}$  are the  best fit values listed in Tab.\ref{tab:1}.
 Black line: $\overline{P^{s}}_{\alpha\beta}$ from the best fit values of the global fit data with NO\cite{36}.
 }
\end{figure}

\subsection{The energy spectrum of $\bar{\nu}_{e}$ modified by the LV}
As has been shown, the LV resulting in the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling  has significant impacts on the flavor transition of HANs.  Thus, the energy spectrum of $\bar{\nu}(e)$
may be modified by the LV effect. In addition to $H_{LV}$, the flux of $\bar{\nu}_{e}$ at Earth is also dependent on the production source of HANS.
Since the source of $\bar{\nu}_{e}$ is uncertain at present, here we consider two typical sources to set the original flavor composition.
They are listed as follows: the pion ($\pi^{\pm}$)  decay source with the original flavor ratio
\begin{equation}
\label{eq:16}
R^{s(\pi^{\pm})}=( 1/3, ~2/3,~ 0 ),
\end{equation}
and the muon ($\mu^{\pm}$) damped source with the ratio
\begin{equation}
\label{eq:17}
R^{s(\mu^{\pm})}=( 0,~ 1,~ 0 ).
\end{equation}
For the both sources, the flux of neutrinos is equal to that of antineutrinos.
Employing the flavor ratio at the source and the transition matrix $\overline{P}$, we can derive the ratio of $\bar{\nu}_{e}$ at Earth, namely
\begin{equation}
\label{eq:18}
r_{e}(E_{\nu})=\frac{\phi_{\bar{\nu}_{e}}(E_{\nu})}{\phi_{\nu+\bar{\nu}}(E_{\nu})}=\frac{1}{2}\sum_{\alpha=e,\mu,\tau}\overline{P}_{e\alpha}(E_{\nu})R^{s}_{\alpha},
\end{equation}
in which $\phi_{\bar{\nu}_{e}}$ and $\phi_{\nu+\bar{\nu}}$ are the  differentia flux  of $\bar{\nu}_{e}$ and that of total HANs respectively,
$R^{s}_{\alpha}$ is the component of the flavor ratio vector at the source.
Accordingly, the differentia flux of $\bar{\nu}_{e}$ is predicted to be
\begin{equation}
\label{eq:19}
\phi_{\bar{\nu}_{e}}=\frac{1}{2}\phi_{\nu+\bar{\nu}}(E_{\nu})\times(\sum_{\alpha=e,\mu,\tau}\overline{P}_{e\alpha}(E_{\nu})R^{s}_{\alpha}~).
\end{equation}

At present, the uncertainty of $\phi_{\nu+\bar{\nu}}$ is large. Several models of the energy spectra  can work, e.g., the single power-law,
double power-law, single power law with spectral cutoff, log-parabola, segmented power-law.
For the sake of illustration, we consider the single power-law and the double power-law spectrum in this paper.
They are expressed respectively as\cite{38}
\begin{equation}
\label{eq:20}
\phi_{\nu+\bar{\nu}}^{spl}=\phi_{astro}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{astro}}\cdot10^{-18}\rm GeV^{-1}cm^{-2}s^{-1}sr^{-1},
\end{equation}
\begin{equation}
\label{eq:21}
\phi_{\nu+\bar{\nu}}^{dpl}=(\phi_{hard}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{hard}}+\phi_{soft}\times(\frac{E_{\nu}}{100\rm TeV})^{-\gamma_{soft}})\cdot10^{-18}\rm GeV^{-1}cm^{-2}s^{-1}sr^{-1}.
\end{equation}
The ranges of the normalization and the spectral index at the 68.3\% confidence level are listed in Tab.\ref{tab:2}.
\begin{table}
\caption{\label{tab:2} Energy-spectrum parameters of HANs at the 68.3\% confidence level\cite{38}. }
  \centering
  \begin{tabular}{c c c c c }
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     Model~~&~~~~~~ normalization  ~~&~~~~&~ ~ spectral index  ~&~~ \\
     \noalign{\smallskip}\hline
     single power-law ~~~& $\phi_{astro}$ ~~~&~&~~~ $\gamma_{astro}$ ~ \\
  ~~~~~&$5.68^{+1.56}_{-1.55}$ ~~~&~&~~~~ $2.89^{+0.23}_{-0.20}$ \\
   \noalign{\smallskip}\hline
       double power-law ~~~&$\phi_{hard}$ ~~&~ $\phi_{soft}$ ~~&~~ $\gamma_{hard}$ &~$\gamma_{soft}$~ \\
    ~~~~&$3.26^{+1.88}_{-2.44}$ ~~&~ $0.08^{+3.64}_{-0.08}$~&~~~~~ $2.78^{+0.24}_{-0.32}$~~&~~~~ $3.12^{+1.01}_{-0.30}$\\
     \hline
   \end{tabular}
\end{table}
Using these parameters, we show the LV effect on the energy spectrum of $\bar{\nu}_{e}$ in Fig.\ref{fig:2} and Fig.\ref{fig:3}.

For the  muon ($\mu^{\pm}$) damped source (see Fig.\ref{fig:3}), the LV can bring a noticeable damping of the flux of $\bar{\nu}_{e}$ when $E_{\nu}\geq E_{d}$.
The damping rate $\phi_{\bar{\nu}_{e}}^{LV}$/ $\phi_{\bar{\nu}_{e}}^{S }$ is of order [$10^{-6}$, $10^{-4}$] at the resonance energy. Here $\phi_{\bar{\nu}_{e}}^{LV}$ and $\phi_{\bar{\nu}_{e}}^{S }$ denote the flux including the LV and that without the LV respectively. The specific  damping value is dependent on the form of $\phi_{\nu+\bar{\nu}}$.
For the  pion ($\pi^{\pm}$) decay source (see Fig.\ref{fig:2}), the damping effect is moderate and concealed by the large uncertainty of $\phi_{\nu+\bar{\nu}}$, irrespective of the magnitude of the decoupling energy.
In general, we can expect that  the LV effect can lead to a noticeable decrease of  the $\bar{\nu}_{e}$ flux around the resonance energy for a source dominated by $\bar{\nu}_{\mu}$ ($\nu_{\mu}$).
Thus, the detection of GR events may set stringent constraints on the texture of $H_{LV}$  and the original flavor ratio.

\begin{figure}
\label{fig:2}
  \centering

  \includegraphics[width=.49\textwidth]{fig2a}
   \hfill
  \includegraphics[width=.49\textwidth]{fig2b}
   \includegraphics[width=.49\textwidth]{fig2c}
   \hfill
  \includegraphics[width=.49\textwidth]{fig2d}
  \caption{\label{fig:2} The energy spectrum of $\bar{\nu}_{e}$ from the pion ($\pi^{\pm}$) decay source. Left plots: from the total flux of the single power-law. Right plots: from the total flux of the double power-law.
   The ranges of the parameters of the energy spectrum $\phi_{\nu+\bar{\nu}}$ are taken as those listed in Tab.\ref{tab:2}.
   Blue lines: arising from $\overline{P^{s}}$ without the LV effects.  Red lines:  from $\overline{P}$ with  the LV effects. For $\overline{P^{s}}$ and $\overline{P}$,
   the leptonic mixing parameters are in the $3\sigma$ allowed ranges of the global fit
   data with NO\cite{36}. The plots in the first row: with the LV parameters  taken as $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$.
   The plots in the second row: $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$.
  The black line: arising from $\overline{P}$ with the best fit values of the energy spectrum $\phi_{\nu+\bar{\nu}}$ and leptonic mixing parameters. The dashed black line: from $\overline{P^{s}}$ with the best fit values of the leptonic mixing parameters.}
\end{figure}
\begin{figure}
\label{fig:3}
  \centering

  \includegraphics[width=.49\textwidth]{fig3a}
   \hfill
  \includegraphics[width=.49\textwidth]{fig3b}
  \includegraphics[width=.49\textwidth]{fig3c}
   \hfill
  \includegraphics[width=.49\textwidth]{fig3d}
  \caption{\label{fig:3} The energy spectrum of $\bar{\nu}_{e}$ from the muon ($\mu^{\pm}$) damped source. The conventions of the parameters, colors and the lines are the same as those in Fig.\ref{fig:2}.}
\end{figure}

\subsection{The strength of the GR event rate damped by the LV}
The GR events provide a window to explore the texture of $H_{LV}$ and the flavor composition of the HANs at the PeV energy-scale.
Now let us examine the impacts of the decoupling of $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) on the GR event rate.
To weaken the influence of the uncertain energy spectrum  $\phi_{\nu+\bar{\nu}}$,  for down-coming events, the ratio of the resonance events to the nonresonant continuum events is proposed
to describe the strength of the GR event rate, i.e.\cite{2},
\begin{equation}
\label{eq:22}
\begin{aligned}
  \frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})}&=\frac{10\pi}{18}(\frac{\Gamma_{W}}{M_{W}})(\frac{\sigma^{peak}_{Res}}{\sigma^{CC}_{\nu N}(E_{\nu}=6.3\rm PeV)})\frac{(\alpha-1.4)(\frac{E_{\nu}^{min}}{6.3\rm PeV})^{\alpha-1.4}}{[1-(\frac{E_{\nu}^{min}}{E_{\nu}^{max}})^{(\alpha-1.4)}]}\times [r_{e}]_{E_{\nu}=6.3\rm PeV} \\
  &=11\times\frac{(\alpha-1.4)(\frac{E_{\nu}^{min}}{6.3\rm PeV})^{\alpha-1.4}}{[1-(\frac{E_{\nu}^{min}}{E_{\nu}^{max}})^{(\alpha-1.4)}]}\times [r_{e}]_{E_{\nu}=6.3\rm PeV}.
\end{aligned}
\end{equation}
Here the spectral index $\alpha$ is dependent on the acceleration mechanism of neutrinos. In absence of GR events, $\alpha\geq2.3$\cite{40,41}. In this paper, we take $\alpha=2$.
We note that a moderate variation of $\alpha$ will not lead to a noticeable change of the observation on the LV effect.
Employing the best fit values of the leptonic mixing parameters with NO\cite{36} and the LV parameters $a^{T}_{\alpha\beta}=10^{-3}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-2}(c^{TT}_{\alpha\beta})_{bf}$ with $(a^{T}_{\alpha\beta})_{bf}$, $(c^{TT}_{\alpha\beta})_{bf}$ listed in Tab.\ref{tab:1}, we calculate the ratio of resonance event rate around the resonance energy.
The specific values are listed in Tab.\ref{tab:3}. For the sake of comparison, the parenthetic values obtained from the standard flavor transition matrix $\overline{P^{s}}$ are also shown.
\begin{table}
\caption{\label{tab:3} The ratio of resonance event rate to nonresonant event rate with and without the LV effect for the  pion ($\pi^{\pm}$) decay source, muon ($\mu^{\pm}$) damped source, and a source rich in $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) with $\phi_{\nu}=\phi_{\bar{\nu}}$. Here $\alpha=2$, $E_{\nu}^{min}$=1, 2, 3, 4, 5 PeV, $E_{\nu}^{max}=\infty$, $r_{e}$ takes the value $ [r_{e}]_{E_{\nu}=6.3\rm PeV}$.
The parenthetic values are obtained from the standard flavor transition matrix $\overline{P^{s}}$. }
  \centering
  \begin{tabular}{c c c c c c c}
     \noalign{\smallskip}\hline
     \noalign{\smallskip}\hline

     $E_{\nu}^{min}(\rm PeV)$ ~&~~1  ~&~2~~&~~3~~&~~ 4~~&~~ 5&~~ \\
     \noalign{\smallskip}\hline
$R^{s(\pi^{\pm})}=( 1/3, ~2/3,~ 0 )$ ~&~~0.18~~&~~~0.28~ &~~~0.35 &~~~0.42&~~~0.48&~$r_{e}$=0.083\\
~~&~~(0.32)~~&~~~(0.50)~ &~~~(0.63) &~~~(0.75)&~~~(0.86)&~($r_{e}$=0.149)\\ \hline
$R^{s(\mu^{\pm})}=( 0,~ 1,~ 0 )$  ~~~~~&~~$7.8\times10^{-6}$~~~~&~~$1.2\times10^{-5}$~~&~$1.5\times10^{-5}$&~$1.8\times10^{-5}$ &~$2.0\times10^{-5}$&~$r_{e}=3.5\times10^{-6}$\\
~~~~&~~~~(0.19)~~~~&~~~(0.28)~~&~(0.36)&~(0.43) &~(0.49)&~($r_{e}$=0.086)\\ \hline
$R^{s(x)}=( 0.1,~ 0.9,~ 0 )$ ~~&~~0.05~~&~~~0.08~~ &~~0.10~~&~~0.12~~&~~0.14~~&~$r_{e}$=0.025\\
~~~~&~~~~~(0.23)~~~~&~~~(0.35)~~&~~(0.44)&~~(0.53) &~~(0.60)&~~($r_{e}$=0.105)\\
     \hline
   \end{tabular}
\end{table}

From Tab.\ref{tab:3}, we can see that the ratio of  resonance event rate to the nonresonant event rate
is mainly determined by the ratio of $\bar{\nu}_{e}$ at the resonance energy-scale. The damping rate of the GR event rate approximates that of $ [r_{e}]_{E_{\nu}=6.3\rm PeV}$,
namely
\begin{equation}
\label{eq:23}
  (\frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})})^{LV}/(\frac{N_{Res}}{N_{non-Res}(E_{\nu}>E_{\nu}^{min})})^{S}\sim ([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}/([r_{e}]_{E_{\nu}=6.3\rm PeV})^{S},
\end{equation}
in which the superscript 'LV' ('S') denotes the value with (without) the LV effect.
For the original flavor ratios $R^{s(\pi^{\pm})}$, $R^{s(\mu^{\pm})}$, and $R^{s(x)}$, $([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}/([r_{e}]_{E_{\nu}=6.3\rm PeV})^{S}$ takes 0.557, $4.07\times10^{-5}$, 0.238, respectively.
The damping rate is sensitive to the increase of the original ratio of $\bar{\nu}_{\mu}$ around the resonance energy when $\phi_{\bar{\nu}_{\mu}}/\phi_{\bar{\nu}}>0.9$.
Furthermore, when the decoupling energy is around $1\rm PeV$, i.e., $a^{T}_{\alpha\beta}=10^{-7}(a^{T}_{\alpha\beta})_{bf}$, $c^{TT}_{\alpha\beta}=10^{-7}(c^{TT}_{\alpha\beta})_{bf}$,
the values of  $([r_{e}]_{E_{\nu}=6.3\rm PeV})^{LV}$ are similar to those listed in Tab.\ref{tab:3}. Thus, the damping rate of the GR event rate is mainly
sensitive to the texture of $H_{LV}$, irrespective the magnitudes of the LV parameters.
However, we should keep in mind  that the magnitudes of the  parameters determine the decoupling energy which impacts the evolution of $\overline{P}$ with $E_{\nu}$.
In the range $E_{\nu}<E_{d}$, the flavor ratio of HANs at Earth is dependent on the magnitudes of the LV parameters.

\subsection{Comments on the realistic source and the $\nu_{\tau}$ event at the PeV energy-scale }

At present, the production source and the acceleration mechanism of HANs are unknown. However, it is widely believed that the strong magnetic field plays an important role in the
acceleration of HANs\cite{42}, e.g. acceleration by the magnetic reconnection\cite{43,44,45}. As is known, the synchrotron radiation of the charged particle in the strong magnetic field becomes significant at high energies. \
We can expect that the muon ($\mu^{\pm}$) damped source may be dominant at the PeV energy-scale.  Accordingly, when $E_{d}<6.3\rm PeV$, the decoupling pattern can show a noticeable impact on the GR events.
Furthermore, we note that the decoupling pattern can bring the following observation:
\begin{equation}
\label{eq:24}
\phi_{\nu_{e}}\sim\phi_{\nu_{\tau}},~~~~\phi_{\bar{\nu}_{e}}\sim\phi_{\bar{\nu}_{\tau}}.
\end{equation}
Hence, for the muon ($\mu^{\pm}$) damped source, the proposed LV texture can lead to a notable decrease of both the GR event rate and the double bang events from $\nu_{\tau}$ at the PeV energy-scale.
In other words, we can obtain the double constraints on the texture of $H_{LV}$ from the GR event and the double bang event.

\section{Conclusions }

The Lorentz invariance is a fundamental symmetry in nature. Strong constraints on the LV parameters have been obtained at low energies.
To connect the LV to the GR event at the PeV energy-scale, we proposed the LV Hamiltonian of the special texture which can  bring the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling  at high energies.
On the base of the $\nu_{\mu}$ ($\bar{\nu}_{\mu}$) decoupling, we analysed the impacts of the LV parameters on the flavor transition of HANs.
For the typical sources of HANs, the damping effect from the decoupling pattern on the energy spectrum of $\bar{\nu}_{e}$ was shown.
For the  muon ($\mu^{\pm}$) damped source, the decrease of the $\bar{\nu}_{e}$ flux at the PeV energy-scale is notable, irrespective of the form of the total flux $\phi_{\nu+\bar{\nu}}$.
The LV effect can lead to a noticeable damping of the GR event rate if the HANs are mainly from the source dominated by $\bar{\nu}_{\mu}$ ($\nu_{\mu}$).
Furthermore, the decoupling pattern can also decrease the double bang events from $\nu_{\tau}$ at the PeV energy-scale.
Therefore, the detection of GR events and other high-energy events may bring extra constraints on the LV Hamiltonian and the production source of HANs.

\acknowledgments
 This work is supported by the National Natural Science Foundation of China under grant No. 12065007, the Guangxi Scientific Programm Foundation under grant No. Guike AD19110045.\\