\section{Introduction}

Mass, charge radius, lifetime, electric (magnetic) transition probability, and deformation are among the most fundamental observables for the many-body nuclear system.
A systematic analysis of these data have been successful in bringing forth a global picture of the atomic nuclei.
For example, experimental binding energies or charge radii for neighboring nuclei do not differ much
except at several specified regions like closed shells or onset of shape.

On the other hand, comparing these observables of atomic nuclei, which differ by one or a few neutrons or protons,
have yielded many empirical relations or filters for special interaction strengths between the valence nucleons.
Of them, the Garvey-Kelson relations for nuclear binding energies are probably one of the best known
examples~\cite{GK1966PRL16,GK-np2001,GK2008PhysRevC.77.041304,GK2011PhysRevC.83.054309,GK2013-PhysRevC.87.057304,GK2013PhysRevC.87.014313,GK2013PhysRevC.88.064325,GK2014PhysRevC.89.061304}.
The validity of this relation in the nuclear charge radius has  also been examined recently~\cite{Piekar2010EPJA46}.

The radius of the charge (proton) distribution can be assumed to be equal to that  of the nuclear mass distribution,
considering the nucleus as a liquid drop with the protons homogeneously distributed over the sphere of the nucleus.
Although not accurate in prediction, this simple liquid-drop model can serve as a guide, and  the interesting physics
can be found in local deviations from the global behavior.
In our recent work~\cite{Sun2014PhysRevC.90.054318,Sun2015PhysRevC.91.019902}
we proposed a set of nuclear charge radius relations $\delta R_{ip-jn}(Z,N)$,
\begin{eqnarray}
   \delta R_{ip-jn}(Z,N)&= &R(Z,N)-R(Z,N-j)     \nonumber \\
                        & &-[R(Z-i,N)-R(Z-i,N-j)]    \nonumber \\
                        &\simeq & 0  \;,
   \label{eq1}
\end{eqnarray}
where $R(Z,N)$ is the root-mean-square (rms) charge radius of the nucleus with $N$ neutrons and $Z$ protons. $i$ and $j$ are integers.
The validity of such relation  is a consequence of the smooth transition in the nuclear structure that is often found when
going from a nucleus to its neighboring nuclei.  Eq.~\ref{eq1} holds precisely over almost the whole nuclear chart
except at a few regions characterized by shape transition and shape coexistence at, e.g., $N \sim 60$, $N \sim 90$ and $Z \sim 80$.
These exceptions raise the possibilities that more accurate local systematics may be developed from the experimental data.

One simple case connecting only even-even nuclei is
\begin{eqnarray}
   \delta R_{2p-2n}(Z,N)&= &R(Z,N)-R(Z,N-2)    \nonumber \\
                        & &-[R(Z-2,N)-R(Z-2,N-2)]   \nonumber \\
                        &\simeq & 0  \;.
   \label{eq2}
\end{eqnarray}
The term $R(Z,N)-R(Z,N-2)$, the so-called isotope shift,
involves the variation of the charge distribution when only two neutrons are added to the system.
In this sense, $\delta R_{2p-2n}(Z,N)$ is nothing but
the difference of isotope shifts for neighboring two isotopic chains. Hereafter we simply
rewrite $\delta R_{2p-2n}(Z,N)$ as $\delta R(Z,N)$.

In this work, we aim to examine and quantify the correlation between the local charge radius relations  $\delta R(Z,N)$
and those of deformation data.  The correlation is made from cases that both charge radius and quadrupole deformation data are experimentally availbale.
This then leads us to an improved relation by correcting the contribution from quadrupole deformation effect in atomic nuclei.
Moreover, this new relation can be naturally extended to the reduced electric quadrupole transition probability $B(E2)$ between the first $2^+$ state and the $0^+$ ground state,
and the mean lifetime $\tau$ of the first 2$^+$ state.

\section{Shape effect on charge radii}\label{sec1}

For the system of spherical nuclei, the rms charge radii can be empirically described by
\begin{equation}\label{eq:Alaw}
R(Z,N)=\sqrt{3/5} r_0 A^{1/3} \;,
\end{equation}
where $A$ is the mass number and $r_0$ is fixed to 1.2 fm throughout this paper.
Thus $\delta R(Z,N)$ for the even-even isotopes is
\begin{equation}
   \delta R(Z,N)   \equiv  \delta R(Z,N)_{\mbox{sph}} = \sqrt{\frac{3}{5}}r_0\delta(A^{1/3}) \;.
   \label{eq4}
\end{equation}
Numerically it is easy to see that  $\delta R(Z,N)_{\mbox{sph}}$ goes down to a few times $10^{-4}$ fm with increasing mass number.

For a deformed nucleus, the charge  radius can be simply decoupled to the spherical and deformation part.
For the important case of an axially symmetric shape,
by neglecting the high-order corrections one can express this approximately as:
\begin{equation}\label{eq.5}
 R(Z,N)=R_0(Z,N)\left[1+\frac{5}{8\pi}\beta_2^2(Z,N)\right] \;
\end{equation}
where $\beta_2(Z,N)$ is the rms quadrupole deformation for the nuclide ($Z,N$). As will be described later,  it is derived experimentally from the  reduced electric quadrupole transition probability $B(E2)$.
 $R_0(Z,N)$ corresponds to the charge radius of a spherical nucleus that the nuclear volume is conserved, and is defined by Eq.~\ref{eq:Alaw}.
Eq.~(\ref{eq.5}) may be generalized to include higher order multipoles or triaxial shape~\cite{Bohr1969book,Greiner1996book}.

Accordingly, the experimental four-point relation $\delta R(Z,N)$ can be expressed in terms of two variables representing the spherical equivalent radius and the deformation.
This corresponds to a two-parameter model of describing $\delta R(Z,N)$,
\begin{eqnarray}
   \delta R(Z,N)   & = &  \delta R(Z,N)_{\mbox{sph}}   +  \delta R(Z,N)_{\mbox{def}}  \nonumber \\
   & = & \delta R(Z,N)_{\mbox{sph}} + \frac{5}{8\pi}\delta(R_0\beta_2^2)  \;,
   \label{eq6}
\end{eqnarray}
where $ \delta R(Z,N)_{\mbox{sph}}$ is defined in Eq.~(\ref{eq4}), and $\delta R(Z,N)_{\mbox{def}}$ comes from the variance of deformation in the relevant nuclei,
\begin{eqnarray}
   \delta R(Z,N)_{\mbox{def}}  & \equiv& \frac{5}{8\pi}\delta(R_0\beta_2^2) \nonumber \\
   &=&  \sqrt{\frac{3}{5}}\frac{5}{8\pi}r_0[A^{1/3}\beta_2^2(Z,N)   \nonumber \\
   &&  + (A-4)^{1/3}\beta_2^2(Z-2,N-2)  \nonumber \\
   && -(A-2)^{1/3}\beta_2^2(Z-2,N)   \nonumber \\
   && - (A-2)^{1/3}\beta_2^2(Z,N-2)]   \nonumber \\
   &\approx&  \sqrt{\frac{3}{5}} \frac{5}{8\pi}r_0A^{1/3}\delta\beta_2^2 \;.
   \label{eq7}
\end{eqnarray}
The approximation of the last term is valid especially for heavier system.
Because of the negligible contribution of $\delta R(Z,N)_{\mbox{sph}}$,  the resulting  $\delta R(Z,N)$ is mostly determined by the terms relevant to nuclear deformation.
Although dynamic deformations and higher-order multipoles are not included in this equation,
they can be subsumed in principle under the deformation term in which $\delta\beta_2^2$ is replaced by $\sum_i\delta\beta_i^2$.

\section{Correlation between charge radius and quadrupole deformation}

\subsection{Experimental data}
\begin{figure}[htbp]\noindent
\centering
\includegraphics[width=0.45\textwidth]{fig1-2016.eps}
 \caption{(Color online) $\delta R(Z,N)$ as a function of $\delta R(Z,N)(\beta_2^2)$ for all experimentally known cases. The linear fit is indicated by the dash line.}\label{fig.1}
\end{figure}

We can now examine the correlation between the experimental $\delta R(Z,N)$ and $\delta R(Z,N)_{\mbox{def}}$. The resulting correlation plot is shown in Fig.~\ref{fig.1}.
The experimental charge radius and deformation data are from the latest evaluations~\cite{Angeli2013,beta2016}.  There are in total 149 even-even nuclei from Ne to Cm.
It is seen that almost all the data follow a linear trend within 1 standard deviation ($\sigma$).  A coefficiency of 0.29(6)  is determined with a reduced $\chi^2$ of 0.8.
This indicates that experimental data of charge radii and quadrupole deformation are in well consistency.

Specifically, this correlation can be seen  for the Sr isotopes in Fig.~\ref{fig.2}.
The sudden onset of the shape transition at $N=60$  is reflected distinctly by both $\delta R(Z,N)_{\mbox{def}}$ and $\delta R(Z,N)$.
The deformation parameters for the relevant 4 nuclei, $^{98}_{38}$Sr, $^{96}_{38}$Sr, $^{96}_{36}$Kr and $^{94}_{36}$Kr are 0.40(1), 0.175(6), 0.25(3), 0.19(1),  resulting in
 $\delta R(Z,N)_{\mbox{def}}$ of  0.088(14)  fm for $^{98}$Sr.
Therefore, once considering the deformation correlation, the large $\delta R(Z,N)$ value at $N=60$, the well-known region of phase transitions, can be largely diminished.
Similar correlation has been observed at $N\sim 90$ for the Nd isotopes. Unfortunately, it is not possible yet to test the region at $Z\sim80$ due to missing deformation data experimentally.

\begin{figure}[htbp]\noindent
\centering
\includegraphics[width=0.45\textwidth]{Z=38.eps}
 \caption{(Color online) Experimentally known $\delta R(Z,N)$ (filled circles) and  $\delta R(Z,N)_{\mbox{def}}$ (open squares) for the Sr isotopic chain.
 }\label{fig.2}
\end{figure}

5 cases at $^{20}_{10}$Ne, $^{46}_{20}$Ca, $^{46}_{22}$Ti,   $^{76}_{34}$Se, and  $^{78}_{36}$Kr,  show deviations from the linear trend by more than 2$\sigma$.
Of them, $^{20}_{10}$Ne is the lightest nuclide with available charge radius and deformation data. It is known already
that the precision of the charge radius formula deteriorates with decreasing mass number (in particularly for $A<60$)~\cite{Sun2014PhysRevC.90.054318}.
This may be understood in the fact that the collective ``deformation'' property is more suitable for heavy nuclei in comparision with lighter nuclei.
Further check on the mass dependence confirms this argument, where 8 of 10 cases  with $A<60$ have deviations from the linear trend by more than 1 $\sigma$.
It should be noted that 1 $\sigma$ discussed here is even as small as 0.0040 fm.

These few cases may be (partially) related to the various
sources of charge radii and deformation, given the different methods of systematic errors associated with different techniques.
We noted that a recent analysis~\cite{BE2-Birch2016145} of $B(E2)$ measurements, from which $\beta_2$ is derived, concluded that most prevalent methods of measuring $B(E2)$
values are equivalent.
Such comparison is not available yet for charge radius data across the entire chart of nuclides.
Anyway, a consistent and equivalent set of nuclear charge radii and $B(E2)$ data are definitely crucial.
A combined analysis of the cases for $^{46}_{20}$Ca and $^{46}_{22}$Ti shows that
increasing the charge radius by 0.3\% or decreasing the $\beta_2$ by about 40\% (unlikely) for $^{44}_{20}$Ca
 will result in their $\delta R(Z,N)$ values in better agreements with the linear trend.
Similar arguments hold also for the cases of $^{76}_{34}$Se and $^{78}_{36}$Kr.
Therefore, the correlation identified here may provide us a very accurate way
to investigate the consistence of both the charge radius and deformation surface.

In Ref.~\cite{Sun2014PhysRevC.90.054318}, it was found that Eq.~(\ref{eq1}) is remarkably successful even at nuclei with magic neutron and/or proton numbers.
This can be easily understood with the correlation identified. Nuclei with magic number of neutrons and/or protons is mostly in the spherical shape, i.e.,
their $\beta_2$ values of less than 0.1, thus leading to naturally a net  $\delta R(Z,N)$.
This is very different from the counterpart in nuclear mass, the valence proton-neutron interactions $\delta V_{pn}$~\cite{Brenner2006PhysRevC.73.034315,Cakirli2005PhysRevLett.94.092501,Chen2009Phys.Rev.Lett.122503}. It depends strongly on the spatial overlap of the valence orbits and presents a dramatic variation when crossing neutron shell closures.

\subsection{Correlations in nuclear models}

The knowledge of both experimental data on charge radii and deformation are still very limited.
It will be very useful if nuclear models can provide these data either in absolute values or in differential values at a reasonable precision. For example,
when one experimental deformation parameter is missing, one can resort to the theoretical predictions in nuclear models.
Care should be taken that quantities like $B(E2)$ refer to the charge (proton) distribution in the nucleus and that in particular $\beta$ is
the charge deformation related to this charge distribution. This should be kept in mind when comparing the experimental results to nuclear models.

\begin{figure}[htbp!]\noindent
\centering
\includegraphics[width=0.45\textwidth]{fig1-RMFbeta2.eps}
\includegraphics[width=0.45\textwidth]{fig1-RMFbeta.eps}
 \caption{(Color online) Same as Fig.~\ref{fig.1} but using the predicted $\delta R(Z,N)$ and $\delta R(Z,N)_{\mbox{def}}$ in the RMF model (a), and
 using experimental $\delta R(Z,N)$ and predicted $\delta R(Z,N)_{\mbox{def}}$ in the RMF model (b). Errorbars are not shown for experimental
 data.}\label{fig.3}
\end{figure}

Global nuclear models can provide both charge radius and deformation data self-consistently. We choose Hartree-Fock-Bogoliubov (HFB-24)
model~\cite{Goriely2013PhysRevC.88.024308} and the relativistic mean-field (RMF) model~\cite{Geng2005PTP} to check the same correlation.  As shown in Fig.~\ref{fig.3}(a),
a linear correlation is predicted between $\delta R(Z,N)$ and $\delta R(Z,N)_{\mbox{def}}$ in the RMF.
The slope is determined to be 0.60, about a factor two larger than that of experimental data.
The difference from the experimental trend should be related to the fact that all nuclei are treated with axially symmetric shapes in the RMF approach.
The same correlation has also been found in the HFB-24 but with a different coefficience parameter (0.85).

However, such correlation vanish once we use the theoretical $\beta_2$ values instead of experimental data.
As an example, $\delta R(Z,N)$ vs. $\delta R(Z,N)_{\mbox{def}}$ calculated using the $\beta_2$ in the RMF is shown in Fig.~\ref{fig.3}(b).
This indicates that current nuclear models are not accurate enough yet to reproduce the correlation seen in the experiment.

\section{Discussion}
\subsection{Improved charge radius formula}

As verified in the previous section, the $\delta R(Z,N)$ can be quantitatively reproduced with the $\delta R(Z,N)_{\mbox{def}}$ for the existing data.
This leads to the improved charge radius formula as follows:
\begin{equation}\label{eq8}
\delta R(Z.N)^{\mbox{corr}}=\delta R(Z,N)-\mbox{C}\cdot \delta R(Z,N)_{\mbox{def}} \approx 0 \;
\end{equation}
where $\mbox{C}$ is 0.29(6) determined from Fig.~\ref{fig.1}. $\delta R(Z,N)$ and $\delta R(Z,N)_{\mbox{def}}$ are given in Eq.~(\ref{eq2}) and  Eq.~(\ref{eq7}), respectively.

The cases with experimental data are used to the check the accuracy of the relation without/with deformation correction.
The weighted mean values of $\delta R(Z,N)^{\mbox{corr}}$ amount to only -8$\times 10^{-4}$ fm with the weighted standard deviation of 5$\times 10^{-3}$ fm.
This is about 15\% improvement in precision comparing with that without correction (i.e., $\delta R(Z,N)$).
The significance is that Eq.~(\ref{eq8}) can be extended to cases even when sudden variances occur in nuclear shapes.

\subsection{Correlation of charge radius with $B(E2)$ and $\tau$}

Experimental quadrupole deformation values are derived from the model-independent
experimental values of the reduced electric quadrupole transition probability $B(E2)$, between
the $0^+$ ground state and the first $2^+$ state in even-even nuclides, using the semi-empirical approach,
\begin{equation}
\beta_2 = \frac{4\pi}{3ZR_0^2}[\frac{B(E2)}{\mbox{e}^2}]^{1/2} \;.
\label{eq:be2}
\end{equation}
Here $R_0=r_0A^{1/3}=1.2A^{1/3}$ fm and $B(E2)$ is in units of e$^2$b$^2$.
It is assuming a uniform charge distribution out to the distance $R$ and zero charge beyond~\cite{Greiner1996book,RAMAN2001ADNDT78}.

The $B(E2)$ values are fundamentally important quantities for determining the collectivity in nuclei.
Moreover, $B(E2)$ are related to  the  mean lifetime $\tau$ of the first $2^+$ state  through
\begin{equation}
  \tau(1+\alpha) = 40.81\times 10^{13}E^{-5}[\frac{B(E2)}{\mbox{e}^2\mbox{b}^2}]^{-1} \;,
\end{equation}
where $E$ is the excitation energy of the first $2^+$ state (in units of keV), and $\tau$ in ps.
The total internal conversion coefficient $\alpha$ for a specific $E$ is needed for correction.

$\delta R(Z,N)_{\mbox{def}}$ can be accordingly rewritten in terms of $B(E2)$ and $\tau$,
\begin{eqnarray}
   \delta R(Z,N)_{\mbox{def}} &=&  4.35\times 10^3 \delta(\frac{B(E2)/ \mbox{e}^2\mbox{b}^2}{Z^2A})\mbox{b}^{1/2} \nonumber \\
   &=& 1.77\times 10^{21}\delta\left(\frac{E^{-5}Z^{-2}A^{-1}}{\tau(1+\alpha)}\right) {\mbox{b}^{1/2}} \;,
   \label{eq9}
\end{eqnarray}
where all the quantities $E$, $B(E2)$, $Z$, $A$ and $\alpha$ in the $\delta$ term are for the four neighboring even-even nuclei.
In case of no abrupt shape transition,   $\delta R(Z,N)_{\mbox{def}} \simeq 0$,  and the following relation involving again four neighboring doubly even nuclei should hold well,
\begin{eqnarray}
 \delta(A^{1/3}\beta_2^2 ) \simeq \delta(\frac{B(E2)}{Z^2A}) \simeq  \delta\left(\frac{E^{-5}Z^{-2}A^{-1}}{\tau(1+\alpha)}\right) \simeq 0 \;.
\end{eqnarray}
For heavy nuclear system, where the difference in  ${Z^2A}$ can be safely neglected, we can then get the relation
\begin{eqnarray}
   \delta B &= &B(E2)(Z,N)-B(E2)(Z,N-2)     \nonumber \\
                        & &-B(2)(Z-2,N)  +B(E2)(Z-2,N-2))    \nonumber \\
                        &\simeq & 0  \;,
   \label{eq11}
\end{eqnarray}
The same relation was proposed independently in Ref.~\cite{BE2relation-PhysRevC.12.2038}.
To examine the validity of  eq.(\ref{eq11}), we use the same data set as in Fig.~1 were used to examine the validity of .

The ``theoretical'' $B(E2)$ value of a given nucleus $(Z,N)$,  $B(E2)_{\mbox{pred}}$,  is calculated  in
terms of experimental $B(E2)$ of its three neighboring nuclei $(Z-2,N)$, $(Z,N-2)$, $(Z-2,N-2)$.  Fig.~\ref{fig.4} shows the
relative differences of the predictions defined as $[B(E2)_{\mbox{pred}}-B(E2)_{\mbox{exp}}]/B(E2)_{\mbox{exp}}$.
It is seen that $B(E2)$  values can be calculated within an accuracy of $\pm$25\%, and often better.
The large deviations are shown at $N\sim 60$ and 90. It should be aware that the precisions of experimental data are getting worse as well at these regions.

\begin{figure}[htbp]\noindent
\centering
\includegraphics[width=0.45\textwidth]{fig-BE2.eps}
 \caption{(Color online)  Relative differences in $B(E2)$ between the predictions of Eq. (12) and experimental data. The $\pm$ 25\% accuracy is indicated by two dotted lines.}\label{fig.4}
\end{figure}

Inserting Eq.~(\ref{eq9}) into Eq.~(\ref{eq8})  we get a correlation between four-point charge radius relation  with that of $B(E2)$ or $\tau$.
This new relation, in principle,  should be more accurate for predictions of unknown  $B(E2)$ values than Eq.~(\ref{eq11}), because
possible shape transitional effect can be at least partially compensated by the relevant charge radius data.
When seven quantities out of  eight , e.g., four charge radii and three $B(E2)$ data, are known, then the last $B(E2)$ can be calculated.
Unfortunately, this will give mostly a too large uncertainty to make meanful predictions.
The uncertainty is mainly propagated from the charge radius data and is typically about 1 order of magnitude higher than that from Eq.~(\ref{eq11}).

\section{Conclusion}

With the available experimental data,  a linear correlation has been found between the charge radius relation $\delta R(Z,N)$ and
the according quadrupole deformation (and thus $B(E2)$) relation $\delta R(Z,N)_{\mbox{def}}$. This correlation can provide a consistence check or analysis of experimental data on charge radii and deformation data.
In the near future it is also interesting to see whether the linear coefficiency 0.29(6) remains for more exotic nuclei, especially at shape transitional regions.

The large deviation of four-point charge radius relation $\delta R(Z,N)$ at shape transitional regions, can be quantitatively reproduced
with the $\delta R(Z,N)_{\mbox{def}}$ when experimental data are available.
This in turn gives an improved charge radius formula, and therefore is very useful to make reliable short-range extrapolations of
charge radii over the nuclear chart. Same correlation has been found in globle nuclear models, but so far the model itself is not accurate enough to reproduce the experimental data. Moreover, the relation can be generalized to a new relation between charge radius and $B(E2)$ or $\tau$. A simple four-point $B(E2)$ relation can reproduce experimental $B(E2)$ values within an accuracy of about $\pm$25\%.

Finally, we would like to mention that a consistent description~\cite{Wood1999Nucl.Phys.A323,E0-PhysRevLett.101.022502,E0-PhysRevC.85.034331,PhysRevC.79.054301,PhysRevC.80.061301,PhysRevC.85.034321,Li2013866,Zhao2014Phys.Rev.C11301}
of radius and transition probabilities of atomic nuclei are important to understand their correlation and thus for a better
interpretation of experimental results. A recent example is shown in $^{111-129}$Cd~\cite{Yordanov2016PhysRevLett.116.032501},
in which the parabolic behavior of the charge radii  is found due to the linear tendency of the quadruple deformation.

\section{Acknowledgments}
This work has been supported by the NSFC
under No. 11235002, 11475014 and National Program on Key Basic
Research Project (2016YFA0400502). The authors thank Z. P. Li, Z. M. Niu, P.W. Zhao, and L. H. Zhu for useful comments.

\input{v1.bbl}
\end{document}


\section{Introduction}
News media are publishing ideas, events and opinions in an increasingly wide range of data formats such as news articles, headlines, videos, tweets, hashtags and others. The explosion of Big news data has sparked the text and data mining research communities to focus on developing systems for news data exploration and analysis. Broadly, two types of news data exploration systems are developed till date: \textit{Event centric} (GDELT \cite{D:8}, EventRegistry \cite{D:7}) and \textit{Content centric} (STICS \cite{D:9}, EMM \cite{D:10}). In Event centric system, input query maps to real world events, whereas, the content centric system outputs related news articles of a given query.

Although, both types of systems provide up-to-date news information in real time, but they overload the user with the large amounts of results. For instance, given input query ``2014 FIFA World Cup'', event centric EventRegistry suggested 11,504 news events, and the content centric STICS suggested 1,286,369 news articles having multiple organizations, people and places mentioned. Clearly, there is a need of a system that enables readers to get a broad overview of the news data generated in response of user query.

In this paper, we propose a keyphrase based news exploration engine E$^3$ to summarize high volume news data. In our context, keyphrase is a short and meaningful chunk of text that describes an important news concepts, news entities, etc. For instance, ``Bihar election'', ``Bihar bjp'' are examples of keyphrases. Our proposed work is keyphrase centric as recent literature has shown that Keyphrase mining is able to generate a wide range of informational and important phrases from large documents (KEA \cite{D:1}, Micro-ngram \cite{D:2}, ToPMine \cite{D:3}, and SegPhrase \cite{D:4}). As shown in Figure~\ref{fig:architecture}, engine E$^3$ works in two phases: keyphrase extraction followed by keyphrase enrichment. Keyphrase extracation is performed on multi-form news (like article, keywords, and others). As existing keyphrase mining approaches performs poorly on news data, we propose a novel keyphrase extraction technique that leverages linguistic-syntactic feature, and performs very well on short texts like news headlines and video captions. The keyphrase enrichment phase finds important and interesting information related to the keyphrase, such as connected entities, novel (emerging), and active news concepts, and the role played by most frequent entities present in set of keyphrases.

In summary, given an input query $q$, engine E$^3$ generates a keyphrase template, as shown in Figure~\ref{fig:example}, where keyphrases are organized into three sections: Type Discovery, Keyphrase Ranking and InfoBox Miner. The keyphrase of type Person, Location and Organization are kept in Type section. The Novel, Active and Frequent keyphrases are annotated in Keyphrase Ranking section. Information like Role played, top most associated Types and Keyphrase for a selected keyphrase is stored in InfoBox section.

\section{SYSTEM OVERVIEW}

Figure~\ref{fig:architecture} gives an overview of E$^3$ architecture.

\begin{figure}[!ht]
\centering
\includegraphics[scale=0.17]{overview}
\vspace{-0.2in}
\caption{$E^3$ System Architecture}
\vspace{-0.2in}
\label{fig:architecture}
\end{figure}

\begin{figure*}[ht!]
\centering
\includegraphics[scale=0.5]{biharelection-1}
\vspace{-0.1in}
\caption{$E^3$ System Working Example for $q$: ``Bihar Election''}
\vspace{-0.2in}
\label{fig:example}
\end{figure*}

\subsection{Data Collection}
For input query $q$, our data collection module prepares related news data published by news media. The module can either use news search engines like Google News, Yahoo News or online news structured repositories like GDELT \cite{D:8}, EMM \cite{D:10}, iMM \cite{D:6}. We use our in-house iMM system which periodically extracts news headlines (video title) along with their URL, publication date and meta-keywords. To prepare news data related to $q$, we select URL, if URL's headlines (video title) or URL's meta-keywords contains $q$.

In summary, for given query $q$, we prepare a dataset $R$ containing several news records $R_q$, where each record is described by $Quintuple$ \{Headline, Keywords, Meta-description, Article, Publication Date\}. For ``Bihar Election'' query, we retrieved 216 records from iMM for further processing.

\subsection{Keyphrase Extraction}
Next, we extract keyphrases from all the records of news data $R$. As meta-keywords are small group of meaningful words, a naive solution is to output the meta-keywords as keyphrases. However, not all records in $R$ contain meta-keyword. As a result, meta-keywords are not sufficient enough to describe the news data completely. For example, around 20\% news headlines, obtained for ``Bihar Election'', do not have a meta-keyword. Hence, we require an efficient selection of the relevant phrases.

We observed that news headlines are short in length and contains special tokens such as colon, apostrophe, quotes, hash, dash to emphasize important information. On the other hand, meta-descriptions and news articles are long passage texts and are governed by grammatical rules. Thus, we propose two different keyphrase extractors to handle both kinds of writing styles.
\begin{itemize}[leftmargin=*]
\item \texttt{Syntactic Extractor} utilises special characters such as colon (:), apostrophe ($'$), quotes (", $'$), hash (\#), dash (-) for keyphrase extraction. In case of colon (dash), the news headline is tokenized into two parts using colon (dash) and both parts are declared as keyphrases. In case of quotes (hash) the part of text enclosed inside the quotes and hashtag containing the hash is declared as keyphrase respectively.

\item \texttt{Linguistic Extractor} applies language specific part of speech (POS) tagging on meta-description and article and then annotate collocated nouns, adjectives, noun apostrophe ($'$) connector and numbers present in the input text. These annotated tokens are further used as keyphrases.

\end{itemize}

At the end, when all the news records in $R$ are processed, we obtain a set of keyphrases $R_k$, along with the number of times they are generated. For ``Bihar Election'' query, we obtained around 6000 keyphrases for further processing.

\subsection{Keyphrase Enrichment}
The size of generated keyphrases $R_k$ may be large and noisy. To resolve this problem, our Keyphrase Enrichment module helps in extracting valuable and actionable information by filtering and ranking the extracted keyphrases. The keyphrases are filtered using news media specific stopwords such as update, video, photo, pti and others. Next, we apply case normalization and remove duplicate keyphrases. At this point noisy keyphrases are removed. The remaining keyphrases are passed through the Type discovery, Keyphrase ranking and InfoBox mining modules.

\begin{itemize}[leftmargin=*]
\item In \texttt{Type\ Discovery} module, NER tagger is used to classify keyphrase into three types: $Person$, $Location$ and $Organization$. As existing NER taggers do not perform well on Indian named entities, we use a separate list\footnote{https://github.com/NikkiJain09/Transliteration} for Indian named entities, prepared through in-house research work. A keyphrase without any above NER type, are termed as a $News\ Concept$. For ``Bihar Election'' query, a sample keyphrases for each type is shown in Figure~\ref{fig:example}.

\item Using \texttt{Keyphrase\ Ranking} module, keyphrases are organized according to the value of frequency, novelty and activeness. The frequency of the keyphrase is already computed during keyphrase extraction process. To compute the value of novelty and activeness, we first extracts the time intervals $q_t$ of an input query $q$, during which the $q$ was highly popular in news headlines. Next, a keyphrase is \textit{novel} (denoted by \# in Figure~\ref{fig:example}) if its frequency is very high in news headlines only during $q_t$. Similarly, a keyphrase is \textit{active} (denoted by + in Figure~\ref{fig:example}) if its frequency is very high around $q_t$. For instance, ``Grand Alliance'' and ``NDA'' are discovered as novel and active keyphrases respectively for ``Bihar Election'' query.

\item \texttt{InfoBox\ Miner} discovers personalized information for selected keyphrase. The InfoBox displays role of keyphrase $k$ that it played with respect to query $q$. The type-wise top most connections $k$ have, determined with help of co-occurrence value of the keyphrases in $R_k$. A phrase frequently located near $k$ in the collected news corpus is labeled as the $k$'s role. Generally, keyphrases with type person and organization are preferred for InfoBox mining. Figure~\ref{fig:example} shows InfoBox for entity ``Narendra Modi'' for query ``Bihar Election''.

\end{itemize}

\section{Summary}
The engine is tested for varying the input query ranging from general topics (e.g., Election, ISIS) to specific topics (e.g., Paris Attack, Gravitational Wave) and compared the results\footnote{goo.gl/yoLXTh} with KEA, ToPMine and Micro-ngram. We found that our system outperforms existing approaches in terms of quality and quantity of keyphrases generated. As our engine is online, we can demonstrate the working to the conference participants.


\section{Introduction}
\setcounter{equation}{0}

In our quest for the unifying understanding of microscopic and macroscopic physics, we are puzzled by the existence of the large universe itself.
Let us denote the Planck scale by $M_P$, the Hubble constant by $H$, and the characteristic scale of Dark Energy $\Lambda$:
\begin{align}
H^2 M_P^2 &\sim \Lambda^4,\hspace{1em}\Lambda \sim m_{\nu}, \notag\\
\big(\frac{\Lambda}{M_P}\big)^4&=\big(\frac{H}{M_P}\big)^2\sim 10^{-120}.
\end{align}
Why is there such a large hierarchy?
Does Dark Energy change with time?
How does inflation theory fit in this hierarchy?
Apparently they are incomprehensible in our present knowledge of field theory and string theory.
Quantum gravity in de Sitter spacetime is relevant to these questions.
There still may be large unknown territory in nonperturbative and nonequilibrium physics.

In de Sitter spacetime, scale independent metric fluctuations are universally generated.
The two-point function of the tensor mode of the metric $\delta g^T$ exhibits secular growth as more super-horizon modes accumulate with cosmic time $t$:
\begin{align}
\langle\delta g^T \delta g^T\rangle\simeq \frac{H^2}{M_P^2}\int_{p_\text{min}/a(t)}^H \frac{dP}{P} = \frac{H^2}{M_P^2}\log \big(a(t)/a_0\big) = \frac{H^2}{M_P^2}H(t-t_0),
\label{RW}\end{align}
where an initial time is introduced as $a_0=e^{Ht_0}=p_\text{min}/H$.
The metric is jittered randomly by the horizon crossing modes like a dust particle  is jittered by molecules.
It performs Brownian motion with fractal dimension two.
This IR quantum effect breaks de Sitter symmetry as the correlators are not invariant under a constant translation of cosmic time.
Here we have fixed the minimum comoving momentum cut-off $p_\text{min}$.
In such a case, the physical minimum momentum scales as $P_\text{min}(t)=p_\text{min}/a(t)\propto 1/a(t)$.
Physically speaking, we consider a situation that a universe with a finite spatial extension starts de Sitter expansion at an initial time $t_0$.

The secular growth was found in the propagator for a massless and minimally coupled scalar field \cite{VF1982,L1982,S1982}.
Scale invariant fluctuations are common between the tensor mode of gravity and the light scalar.
Although we work with the Poincar\'e coordinate, we believe it does not miss important degrees of freedom as the large $t$ behavior (\ref{RW}) also holds in the global coordinate \cite{AF}.

A possible way to eliminate this secular growth is to fix the physical momentum cut-off $P_\text{min}$.
This procedure may not be consistent with unitarity in a more generic context.
In inflation theory super-horizon modes may eventually come back into the cosmic horizon.
It is argued that de Sitter symmetry breaking effects are not observable as the super-horizon modes are pure gauge \cite{HMM,Urakawa2009}.
A caveat of this argument is that the relevant gauge transformation diverges at spatial infinity, namely it is a large gauge transformation.
We do not throw away super-horizon modes in inflation theory.
We believe it is a consistent strategy to investigate the entire evolution of a universe of a finite spatial extension when it started de Sitter expansion.

We have investigated IR logarithmic effects due to the quantum fluctuation of the metric in Schwinger--Keldysh formalism where the both metric and matter are quantized.
We evaluated the quantum equations of motion, namely the first derivative of the effective action.
The matter fields do not give rise to such effects unless it is a minimally coupled massless scalar field.
It thus appears to be enough to integrate metric fluctuations only leaving matter fields as classical\footnote{The IR logarithmic effects of a minimally coupled massless scalar field can be included separately and they are manifestly Lorentz invariant.}.
In the case of the kinetic term, we have considered two types of diagrams: Type A and Type B.
Type A is the quantum average over the metric fluctuation in the  local vertex.
Type B corresponds to the intermediate states of the matter and the soft graviton.
Physically it is difficult to distinguish a particle state from the same state with a soft metric fluctuation.
In this sense Type B diagram is hard to distinguish from a one particle reducible diagram.

This fact may be reflected in the quantization ambiguity of the matter field \cite{KitamotoPa}.
As it turned out, the IR effects depend on the parametrization of the matter field with respect to the conformal metric fluctuation.
Let us consider a scalar field $\phi$ for example.
We could rescale the field by the conformal mode of the metric $\omega$ as $e^{\alpha\omega}\phi$.
The problem is that IR effects are found to be sensitive to this parameter $\alpha$.
Physically speaking, the matter field is accompanied by soft metric fluctuations and IR effects depend on how to specify them.
We have proposed that the ambiguity is fixed by requiring that the IR logarithmic effect preserves Lorentz symmetry at a sub-horizon scale.

We recall here that the propagators of the metric are not Lorentz covariant at the super-horizon scale in de Sitter spacetime.
We have adopted a non-covariant gauge which provides us the simplest and complete description of the super-horizon modes.
In this gauge, the conformal mode is not a singlet under full Lorentz transformation but only under spatial rotation.

Lorentz invariance is one of the fundamental principles of the quantum field theory and microscopic physics.
In general relativity, Lorentz symmetry must hold at least locally when the spacetime curvature can be ignored.
In the case of the de Sitter spacetime, it corresponds to the sub-horizon scale.
In fact the requirement of Lorentz symmetry at this limit follows from the fundamental principle of general relativity.
From microscopic quantum field theory point of view, general covariance is a necessary consequence of Lorentz invariance for massless spin-2 particles \cite{SW}.
As the physical degrees of freedom: $\pm 2$ helicity states being the small representation of the Poincar\'{e} group, the general covariance is required to ensure the unitarity of the theory.

We argue that IR logarithmic effects cannot spoil this important symmetry.
Although it is suppressed by a factor $H^2/M_P^2$, it could become large at late times due to the logarithmic factor $\log a(t)$.
In the literature, different IR logarithmic corrections are reported with different parametrization of the matter field \cite{KitamotoSD,KitamotoG}, \cite{Kahya2007,Miao2005}.
Nevertheless the unique result is obtained by requiring Lorentz symmetry after reparametrizing the matter field with respect to the conformal mode \cite{KitamotoPa}.

Our requirement of Lorentz invariance at the sub-horizon scale effectively minimizes the IR logarithmic effect.
In fact we have shown that we can eliminate IR logarithmic effect in the free field theories after time dependent wave function renormalization.
Nevertheless it is not possible to do so in the interacting theories as couplings acquire time dependence.

In this paper, we construct an effective Lagrangian in de Sitter spacetime by integrating super-horizon mode of the metric fluctuation.
In this approach, we consider the contribution from Type A diagram only.
The Type B diagram contains both local and non-local contributions.

We postulate that the effect of Type B diagram can be reproduced by a local redefinition of fields as far as IR effects are concerned.
We represent this effect by rescaling the matter field by the conformal mode $\phi\rightarrow e^{\alpha\omega}\phi$.
We require that the IR logarithmic effects respect the Lorentz symmetry at the sub-horizon scale.
This requirement uniquely fixes the rescaling freedom of the matter field.
We show that this strategy leads to the identical physical predictions with the Schwinger--Keldysh approach under the same requirement of Lorentz symmetry at the sub-horizon scale.
In this approach, locality of IR logarithms is manifest.
The validity of our fundamental postulate is underscored by the consistency with the existing results.

In section 2, we show IR logarithmic effects are local by investing Schwinger--Dyson equations.
In section 3, we construct an effective Lagrangian.
In section 4, we examine gauge dependence of physical observables.
We conclude in section 5.
We point out a similarity of our strategy to construct an effective Lagrangian in de Sitter spacetime to 2d quantum gravity \cite{DDK}.
In 2d gravity, there exits a conformal mode dressing ambiguity with respect to an operator $\int d^2x\ e^{2\omega}O(x)\rightarrow \int d^2x\ e^{\beta\omega}O(x)$.
The dressing parameter $\beta$ is fixed by requiring that this operator becomes conformally invariant.
This requirement follows from the consistency with the general covariance.
Here we require Lorentz invariance of the effective Lagrangian at the sub-horizon scale.
It also follows from the consistency with the general covariance.
This condition fixes the rescaling ambiguity of the matter fields by the conformal mode in our construction of the effective action.

\section{Locality of IR logarithms}
\setcounter{equation}{0}

In this section, we show that IR logarithmic effects are local as they cancel in non-local contributions.
We have investigated the two-point function of the conformally coupled scalar field $G^{-+}(x_1,x_2)=\langle\phi (x_1)\phi(x_2)\rangle$ to examine the one-loop metric fluctuation effect \cite{Kitamoto2013}.
It obeys the following Schwinger--Dyson equation:
\begin{align}
&\ G_0^{-1}|_{x_1}G^{-+}(x_1,x_2)-G_0^{-1}|_{x_2}G^{-+}(x_1,x_2) \notag\\
=&\ \int d^4x'\ \Sigma_\text{4-pt}(x_1,x')G^{-+}_0(x',x_2) \notag\\
&-\int d^4x'\ G^{-+}_0(x_1,x')\Sigma_\text{4-pt}(x',x_2) \notag\\
&+\int d^4x'\ \Sigma^R_\text{3-pt}(x_1,x')G^{-+}_0(x',x_2)+\int d^4x'\ \Sigma^{-+}_\text{3-pt}(x_1,x')G^A_0(x',x_2) \notag\\
&-\int d^4x'\ G^{-+}_0(x_1,x')\Sigma^A_\text{3-pt}(x',x_2)-\int d^4x'\ G^R_0(x_1,x')\Sigma^{-+}_\text{3-pt}(x',x_2),
\label{B0}\end{align}
where $G_0^{-1}=i(\partial_0^2-\partial_i^2)$ is the Laplacian and $\Sigma$ denotes the self-energy of the scalar field due to the metric fluctuation.
$R$ and $A$ denote retarded and advanced Green functions respectively.

In momentum space,
\begin{align}
&2i\partial_{\tau_c}\partial_{\Delta \tau}G^{-+}(\tau_1,\tau_2,\vec{p}) \notag\\
=&\ \Sigma_\text{4-pt}(\tau_1)G^{-+}_0(\tau_1,\tau_2,\vec{p})
- G^{-+}_0(\tau_1,\tau_2,\vec{p})\Sigma_\text{4-pt}(\tau_2) \notag\\
&+\int d\tau'\ \Sigma^R_\text{3-pt}(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
+\int d\tau'\ \Sigma^{-+}_\text{3-pt}(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p}) \notag\\
&-\int d\tau'\ G^R_0(\tau_1,\tau',\vec{p})\Sigma^{-+}_\text{3-pt}(\tau',\tau_2,\vec{p})
-\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})\Sigma^A_\text{3-pt}(\tau',\tau_2,\vec{p}).
\label{Bmn}\end{align}
where $\tau_c=(\tau_1+\tau_2)/2$ and $\Delta \tau=\tau_1-\tau_2$.
This is a standard procedure to derive Boltzmann type equations.
The left-hand side picks up $\tau_c$ dependence of the propagator.
It is caused by the ``collision terms'' on the right-hand side.

They are integrated over the conformal time of the interaction vertices.
The self-energy part $\Sigma$ contains the IR singularity $\int d^3q/q^3$ due to super-horizon mode of the metric.
We focus on the non-local contributions in the ``collision terms''.
The IR singularity of them can be estimated as\footnote{In our gauge, IR singularity in $\Sigma_\text{non-local}$ cancels as  $\langle h^{\mu\nu}h^{\rho\sigma}(q)\rangle p_\mu p_\nu p_\rho p_\sigma \sim 0$.
It is no longer the case in a different gauge.}
\begin{align}
\Sigma_\text{non-local}\sim \int d^3q\ \langle h^{\mu\nu}h^{\rho\sigma}(q)\rangle p_{\mu}p_{\nu}p_{\rho}p_{\sigma}G_0(\tau_1,\tau_2,\vec{p}).
\label{NLSE}\end{align}

First we observe that such IR singularities cancel in the last two lines of (\ref{Bmn})
\begin{align}
&\int d^3q\ \langle h^{\mu\nu}h^{\rho\sigma}(q)\rangle p_{\mu}p_{\nu}p_{\rho}p_{\sigma}\times \notag\\
&\Big\{\int d\tau'\ G_0^R(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
+\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p}) \notag\\
&-\int d\tau'\ G^R_0(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
-\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p})\Big\}.
\end{align}
They represent a particle $\rightarrow$ a particle $+$ a soft graviton and the inverse process.
The signs are opposite since the process and the inverse process decrease and increase the weight of a particle state.
The cancellation takes place as we cannot distinguish a particle and a particle with a soft graviton.

Second, the integration over the conformal time could still produce IR divergences from large negative region.
They also cancel between real and virtual contributions since such a cancellation can be manifestly demonstrated in the last two lines of (\ref{Bmn}) as
\begin{align}
\int_{\tau_2}^{\tau_1} d\tau'\ \Sigma^{-+}_\text{3-pt}(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
-\int_{\tau_2}^{\tau_1}  d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})\Sigma^{-+}_\text{3-pt}(\tau',\tau_2,\vec{p}).
\label{RVC}\end{align}
On the other hand it has been explicitly shown that the local contributions lead to IR logarithmic effects:
\begin{align}
&\partial_{\tau_c}Z=\frac{3}{8}\frac{\kappa^2H^2}{4\pi^2}\frac{1}{\tau_c}, \notag\\
&Z=1-\frac{3}{8}\frac{\kappa^2H^2}{4\pi^2}\log a(\tau_c),
\label{Z}\end{align}
where $Z$ is the overall normalization (wave function renormalization) factor of the propagator.

Let us consider the perturbation by a constant metric $\bar{h}^{\mu\nu}$.
Namely we replace the super-horizon mode by a constant metric.
The cubic term is
\begin{align}
\frac{1}{2}\int d^4x\ \bar{h}^{\mu\nu}\partial_{\mu}\phi\partial_{\nu}\phi.
\end{align}
The second order self-energy is
\begin{align}
\Sigma \sim \bar{h}^{\mu\nu}\bar{h}^{\rho\sigma}p_{\mu}p_{\nu}p_{\rho}p_{\sigma}G_0(\tau_1,\tau_2,\vec{p}).
\end{align}
The time dependence of the two-point function is estimated as
\begin{align}
&\bar{h}^{\mu\nu}\bar{h}^{\rho\sigma}p_{\mu}p_{\nu}p_{\rho}p_{\sigma}\times \notag\\
&\Big\{\int d\tau'\ G_0^R(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
+\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p}) \notag\\
&-\int d\tau'\ G^R_0(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
-\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p})\Big\}.
\end{align}
Therefore the cancellation of IR logarithm: $\log a(\tau_c)$ in non-local contributions holds if $\tau_c$ dependence cancels with a constant metric perturbation.
The latter cancellation follows from the time translation invariance with a constant metric perturbation.
It is clear that this argument can be generalized beyond the one-loop level.

We can make a correspondence between non-local contributions of the IR metric fluctuations and
those of the constant metric perturbation to all orders through identifying super-horizon modes with a constant background.
So the cancellation of the IR singularity in the non-local contributions follows from the absence of $\tau_c$ dependence due to a constant background.
Hence the former is suppressed by a factor $\Delta \tau/\tau_c\ll 1$ in the sub-horizon scale where the IR logarithmic factor  $\log a$ changes slowly.
It is because super-horizon modes are not quite constant but jolted randomly with the characteristic scale of $H$.
We conclude that IR logarithms must cancel in the non-local contributions in generic situations.

In order to understand the cancellation mechanism of IR logarithms in non-local contributions, we may consider a simpler model.
We replace metric fluctuations by that of the minimally coupled scalar field $\varphi$ with a cubic interaction to $\phi$
\begin{align}
\lambda_3\int d^4x\ a^2(\tau ) \varphi\phi^2.
\end{align}
The Schwinger--Dyson equation is
\begin{align}
&2i\partial_{\tau_c}\partial_{\Delta \tau}G^{-+}(\tau_1,\tau_2,\vec{p}) \notag\\
=&\ \int d\tau'\ \Sigma^R(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
+\int d\tau'\ \Sigma^{-+}(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p}) \notag\\
&-\int d\tau'\ G^R_0(\tau_1,\tau',\vec{p})\Sigma^{-+}(\tau',\tau_2,\vec{p})
-\int d\tau'\ G^{-+}_0(\tau_1,\tau',\vec{p})\Sigma^A(\tau',\tau_2,\vec{p}),
\label{BYK}\end{align}
where $\Sigma$ has the following IR singularity
\begin{align}
\Sigma(\tau_1,\tau_2,\vec{p})\sim \lambda_3^2\int \frac{d^3q}{q^3}\ \frac{1}{\tau_1^2\tau_2^2}G_0(\tau_1,\tau_2,\vec{p}).
\end{align}
We observe again that the IR singularity in the integrand cancel on the right-hand side of (\ref{BYK}) at sub-horizon scale where $\Delta \tau/\tau_c \ll1$
\begin{align}
&\lambda_3^2\int \frac{d^3q}{q^3}\times \notag\\
&\Big\{\int d\tau'\ \frac{1}{\tau_1^2\tau'^2}G^R_0(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
+\int d\tau'\ \frac{1}{\tau_1^2\tau'^2}\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p}) \notag\\
&-\int d\tau'\ \frac{1}{ \tau'^2\tau_2^2}\ G^R_0(\tau_1,\tau',\vec{p})G^{-+}_0(\tau',\tau_2,\vec{p})
-\int d\tau'\ \frac{1}{\tau'^2\tau_2^2}\ G^{-+}_0(\tau_1,\tau',\vec{p})G^A_0(\tau',\tau_2,\vec{p})\Big\}.
\end{align}
The integration of the conformal time does not give rise to IR singularity as an analogous formula with (\ref{RVC}) holds.

Note that $\Sigma$ is identical with the second order self-energy due to a mass perturbation under the identification $m^4=\lambda_3^2\int d^3q/q^3$.
Therefore the cancellation of IR logarithms holds in this model if $\tau_c$ dependence cancels with a mass perturbation.
The latter cancellation is a consequence of de Sitter symmetry of a mass perturbation.
de Sitter symmetry implies no time dependence of microscopic physics as $\tau_c$ can be scaled out.
We conclude that the cancellation of IR logarithms with the cubic interaction follows from this correspondence.

Note that the mass term can be identified with the constant expectation value of the scalar field in the cubic interaction $m^2\sim \lambda_3\bar{\varphi}$.
We can make a correspondence between diagrams of the cubic interaction and those of the mass perturbation to all orders through this identification
\begin{align}
m^{2n}\sim  \lambda_3^{n}\bar{\varphi}^{n}.
\end{align}
So the cancellation of the IR singularity in the non-local
contributions follows from the absence of $\tau_c$ dependence due to the constant background.
Since a constant background preserves de Sitter symmetry, it must be the case in general.

Let us come back to the IR logarithms due to the metric fluctuations.
It remains the same that the cancellation of the IR singularity in the non-local contributions follows from the absence of $\tau_c$ dependence due to a constant background metric perturbation.
It also remains true that no $\tau_c$ dependence arise for a constant background metric perturbation since it preserves de Sitter symmetry.
We conclude that IR logarithms must cancel in the non-local contributions with metric fluctuations in generic situations.

It is not possible to distinguish one particle state from the same state with a soft graviton.
The situation is similar to one particle state with soft photons in QED \cite{Kinoshita1962,Lee1964}.
It is the physical reason behind the cancellation of IR logarithms in the non-local contributions.
The cancellation of non-local IR logarithms follows from the absence of $\tau_c$ dependence with a constant background.
We believe that the cancellation of IR divergences at large negative conformal time region is a universal phenomena in unitary field theories.

On the other hand, there is no such cancellation mechanism for local contributions.
The self-energy for the local contributions can be estimated as
\begin{align}
\Sigma_\text{local}(x,x')\sim i\delta^{(4)}(x-x')\times\big\{\log \big(a(\tau')/a_0\big) \partial'^2-\partial'_0\log \big(a(\tau')/a_0\big)\cdot\partial'_0\big\},
\end{align}
where the IR logarithm comes from the gravitational propagator at the coincident point.
Since the IR logarithm is not quite constant, the time derivative of it is nonzero.
Considering partial integrations, we can confirm that the local contributions do not cancel out:
\begin{align}
&\ \int d^4x'\ \Sigma_\text{local}(x_1,x')G^{-+}(x',x_2)-\int d^4x'\ G^{-+}(x_1,x')\Sigma_\text{local}(x',x_2) \notag\\
\sim&-2iHa(\tau_c)\partial_0|_{x_1}G(x_1,x_2).
\end{align}
That is how the local contributions lead to the time dependent normalization factor (\ref{Z}).
The IR logarithmic behavior of the propagator is a consequence of the accumulation of the super-horizon modes with cosmic evolution.
We therefore argue that the locality of IR logarithmic effect holds in a generic situation.

\section{Effective Lagrangian}\label{Effective}
\setcounter{equation}{0}

In this section, we construct an effective Lagrangian in de Sitter spacetime which incorporates
IR logarithmic effects.
The de Sitter metric is
\begin{align}
ds^2&=-dt^2+a^2(t)dx^idx_i \notag\\
&=a^2(\tau)(-d\tau^2+dx^idx_i).
\end{align}
We parametrize the metric including the quantum fluctuation as
\begin{align}
g_{\mu\nu}&=a^2(\tau )e^{2\omega}\tilde{g}_{\mu\nu}, \notag\\
\text{det}\ \tilde{g}_{\mu\nu}&=-1,\hspace{1em}\tilde{g}_{\mu\nu}= (e^h)_{\mu\nu},
\end{align}
where $h_{\mu\nu}$ and $\omega$ represent traceless and conformal modes of the metric respectively.
In this paper, the Lorentz indices are raised and lowered by the flat metric $\eta^{\mu\nu}$ and $\eta_{\mu\nu}$ respectively.

The Lagrangian of Einstein gravity on the $4$-dimensional de Sitter background is
\begin{align}
\mathcal{L}_\text{gravity}&=\frac{1}{\kappa^2}\sqrt{-g}\big[R-6H^2\big] \notag\\
&=\frac{1}{\kappa^2}\big[\Omega^2\tilde{R}
+6\tilde{g}^{\mu\nu}\partial_\mu\Omega\partial_\nu\Omega-6H^2\Omega^4\big],
\label{gravity}\end{align}
where $\tilde{R}$ denotes the Ricci scalar constructed from $\tilde{g}_{\mu\nu}$.
The gravitational coupling $\kappa$ is related to Newton's constant as $\kappa^2=16\pi G_N$.

In order to fix the gauge with respect to general coordinate invariance, we adopt the following gauge fixing term \cite{Tsamis1992}:
\begin{align}
\mathcal{L}_\text{GF}&=-\frac{1}{2}a^{2}F_\mu F^\mu, \notag\\
F_\mu&=\partial_\rho h_\mu^{\ \rho}-2\partial_\mu \omega+2h_\mu^{\ \rho}\partial_\rho\log a+4\omega\partial_\mu\log a.
\label{GF}\end{align}

We decompose the spatial part of the metric as
\begin{align}
h^{ij}=\tilde{h}^{ij}+\frac{1}{3}h^{kk}\delta^{ij}=\tilde{h}^{ij}+\frac{1}{3}h^{00}\delta^{ij}.
\end{align}
After diagonalizing the quadratic action in terms of
\begin{align}
X=2\sqrt{3}\omega-\frac{1}{\sqrt{3}}h^{00},\hspace{1em}Y=h^{00}-2\omega,
\label{diagonalize}\end{align}
we find that some metric modes behave as the massless and minimally coupled scalar field $\varphi$ and the other modes behaves as the massless and conformally coupled mode $\phi$:
\begin{align}
\langle X(x)X(x')\rangle&=-\langle\varphi(x)\varphi(x')\rangle, \notag\\
\langle\tilde{h}^i_{\ j}(x)\tilde{h}^k_{\ l}(x')\rangle&=(\delta^{ik}\delta_{jl}+\delta^i_{\ l}\delta_j^{\ k}-\frac{2}{3}\delta^i_{\ j}\delta^k_{\ l})\langle\varphi(x)\varphi(x')\rangle, \notag\\
\langle b^i(x)\bar{b}^j(x')\rangle&=\delta^{ij}\langle\varphi(x)\varphi(x')\rangle,
\label{minimally}\end{align}
\begin{align}
\langle h^{0i}(x)h^{0j}(x')\rangle&=-\delta^{ij}\langle\phi(x)\phi(x')\rangle, \notag\\
\langle Y(x)Y(x')\rangle&=\langle\phi(x)\phi(x')\rangle, \notag\\
\langle b^0(x)\bar{b}^0(x')\rangle&=-\langle\phi(x)\phi(x')\rangle,
\label{conformally}\end{align}
where $b$, $\bar{b}$ denote the ghost and anti-ghost fields.

Since we focus on the de Sitter symmetry breaking effects, we may introduce an approximation.
We can neglect the conformally coupled modes of gravity (\ref{conformally}) since they do not induce the IR logarithm.
In such an approximation, the following identity holds
\begin{align}
h^{00}\simeq 2\omega\simeq\frac{\sqrt{3}}{2}X.
\label{diagonalize1}\end{align}
Note that we are still left with scalar $X$ and spin-2 modes $\tilde{h}^{ij}$ which contain the tensor modes.
Their total contribution is found to preserve Lorentz symmetry at sub-horizon scale in this gauge.

We consider a generic renormalizable Lagrangian just like the standard model:
\begin{align}
\int \sqrt{-g}d^4x\ \big[&-g^{\mu\nu}D_{\mu}\phi (D_{\nu}\phi)^*-(\frac{1}{6}R+m^2)\phi\phi^*+i\bar{\psi}e^\mu_{\ a}\gamma^aD_{\mu}\psi - m_f\bar{\psi}\psi \notag\\
&-\frac{1}{2}\lambda_4(\phi\phi^*)^2 -\lambda_Y\phi\bar{\psi}\psi+\text{(h.c.)}-\frac{1}{4g^2}g^{\mu\rho}g^{\nu\sigma}F^a_{\mu\nu}F^a_{\rho\sigma}\big].
\end{align}
We assume that there are no massless minimally coupled scalar field.

We propose to construct the effective Lagrangian as
\begin{align}
\mathcal{L}_\text{eff}=\langle\mathcal{L}\rangle_\text{metric},
\label{EL}\end{align}
where the average is taken over the super-horizon metric fluctuation only.
The effective Lagrangian $\mathcal{L}_\text{eff}$ summarizes IR logarithmic effects in de Sitter spacetime.
The justification of this proposal is that IR logarithmic effects are contained in the metric fluctuations only.
The effect of the matter field fluctuation should be neglected except possible local redefinition of fields.
When we construct the effective action, we rescale the fields as
\begin{align}
\phi &\rightarrow a^{-1}e^{(\alpha-1)\omega}\phi, \notag\\
\psi &\rightarrow a^{-\frac{3}{2}}e^{(\beta-\frac{3}{2})\omega}\psi, \notag\\
A_{\mu} &\rightarrow A_{\mu}.
\label{rescale}\end{align}
We point out that the gauge field is protected against the rescaling ambiguity due to gauge symmetry.

After the field redefinition, we obtain
\begin{align}
\int d^4x\ \big[&-e^{2\alpha\omega}\tilde{g}^{\mu\nu}D_{\mu}\phi (D_{\nu}\phi)^*-m^2a^2e^{2(1+\alpha)\omega}\phi\phi^*
+ie^{2\beta\omega}\bar{\psi}\tilde{e}^\mu_{\ a}\gamma^aD_{\mu}\psi + m_fae^{(1+2\beta)\omega}\bar{\psi}\psi \notag\\
&-\frac{1}{2}e^{4\alpha\omega}\lambda_4(\phi\phi^*)^2 -\lambda_Ye^{(\alpha+2\beta)\omega}\phi\bar{\psi}\psi+\text{(h.c.)}
-\frac{1}{4g^2}\tilde{g}^{\mu\rho}\tilde{g}^{\nu\sigma}F^a_{\mu\nu}F^a_{\rho\sigma}\big].
\label{rescaled}\end{align}
We have neglected derivatives of the conformal mode as we focus on the IR effects.
Our strategy is to adjust $\alpha$ and $\beta$ to preserve Lorentz invariance at the sub-horizon scale.
The metric propagators are suppressed by $H^2/M_P^2$ factors.
So the leading IR effects appear as $(H^2/M_P^2\cdot\log a)^n$ at the $n$-th loop.
We investigate the leading one-loop effect.

As an illustration, let us consider the kinetic term of the scalar field in (\ref{rescaled})
\begin{align}
&-\langle e^{2\alpha\omega}\tilde{g}^{\mu\nu}\rangle D_{\mu}\phi (D_{\nu}\phi)^* \notag\\
\sim &-\langle e^{2\alpha\omega}\tilde{g}^{00}\rangle D_0\phi (D_0\phi)^*
-\langle e^{2\alpha\omega}\tilde{g}^{ij}\rangle D_i\phi (D_j\phi)^*.
\end{align}
The $\alpha=0$ contribution is
\begin{align}
-\frac{3}{8}\langle\varphi^2\rangle D_0\phi (D_0\phi)^*-\frac{13}{8}\langle\varphi^2\rangle D_i\phi (D^i\phi)^*,
\end{align}
where $\langle\varphi^2\rangle=\frac{\kappa^2H^2}{4\pi ^2}\log a(\tau )$.
The linear term in $\alpha$ is
\begin{align}
-\frac{3\alpha}{4}\langle\varphi^2\rangle D_0\phi (D_0\phi)^*-\frac{\alpha}{4}\langle \varphi^2\rangle D_i\phi (D^i\phi)^*.
\end{align}
We find that the requirement of Lorentz invariance fixes $\alpha=-2$.
The total result including $\alpha^2$ effect is
\begin{align}
-\frac{3}{8}\langle\varphi^2\rangle D_0\phi (D_0\phi)^*+\frac{3}{8}\langle\varphi^2\rangle D_i\phi (D^i\phi)^*.
\end{align}
We thus find that the IR logarithmic effect can be canceled by the time dependent wave function renormalization of
$\phi\rightarrow Z\phi$ where $Z^2=(1+\frac{3}{8}\langle\varphi^2\rangle)$.

As we have seen, the IR logarithmic effects can be eliminated for free fields in this way.
However it is not so for interacting fields.
Let us consider the scalar quartic coupling in our parametrization with canonically normalized kinetic term
\begin{align}
\lambda_4 (\phi\phi^*)^2 \rightarrow \lambda_4 Z^4e^{-8\omega}(\phi\phi^*)^2.
\end{align}
We find that the coupling decreases with time
\begin{align}
\lambda_4 Z^4\langle e^{-8\omega}\rangle=\lambda_4(1-\frac{21}{4}\langle\varphi^2\rangle),
\end{align}
in agreement with \cite{KitamotoSD}.

We next consider the mass term
\begin{align}
m^2a^2e^{2\omega}\phi\phi^* \rightarrow m^2a^2 e^{-2\omega}Z^2\phi\phi^*.
\end{align}
Since $\langle e^{-2\omega}\rangle Z^2\sim 1$, the mass term is not renormalized after the wave function renormalization.

Subsequently we consider the kinetic term of the Dirac field in (\ref{rescaled})
\begin{align}
S=\int d^4x\ ie^{2\beta\omega}\bar{\psi}\tilde{e}^\mu_{\ a}\gamma^a D_\mu\psi.
\label{Dirac}\end{align}
The analogous considerations with the the scalar field fixes $\beta=-1$.
The IR logarithmic contribution to the kinetic term is
\begin{align}
-\frac{3}{32}\langle \varphi^2\rangle i\bar{\psi}\gamma^{\mu} D_\mu\psi.
\end{align}
This effect can be canceled after the time dependent wave function renormalization of $\psi\rightarrow Z_{\psi}\psi$ as $Z_{\psi}^2=(1+\frac{3}{32}\langle\varphi^2\rangle)$.

We consider the Yukawa coupling in our parametrization with the canonically normalized kinetic terms for the both scalar and Dirac fields
\begin{align}
\lambda_Y \phi\bar{\psi}\psi\rightarrow \lambda_Y Z Z_\psi^2e^{-4\omega}\phi\bar{\psi}\psi.
\end{align}
We find that the Yukawa coupling decreases with time
\begin{align}
\lambda_Y Z Z_\psi^2\langle e^{-4\omega}\rangle=\lambda_Y (1-\frac{39}{32}\langle\varphi^2\rangle),
\end{align}
in agreement with \cite{KitamotoSD}.

We next consider the mass term
\begin{align}
m_f ae^{\omega}\bar{\psi}\psi \rightarrow m_f ae^{-\omega}Z_{\psi}^2\bar{\psi}\psi.
\end{align}
Since $\langle e^{-\omega}\rangle Z_{\psi}^2\sim 1$, the mass term is not renormalized after the wave function renormalization.

These considerations parallel our physical interpretations of the IR logarithmic effects in terms of the conformal mode dynamics \cite{review}.
This work may be regarded as a more precise statement of such an idea.

We consider the kinetic term of the gauge field in (\ref{rescaled})
\begin{align}
-\frac{1}{4g^2}\langle\tilde{g}^{\mu\rho}\tilde{g}^{\nu\sigma}\rangle F^a_{\mu\nu}F^a_{\rho\sigma}=-\frac{1}{4g^2}(1+\frac{3\kappa H^2}{8\pi^2}\log a(\tau))F^a_{\mu\nu}F^{a\mu\nu}.
\end{align}
Thus IR effects screen the gauge coupling with time.
\begin{align}
g^2(\tau) =g^2(1-\frac{3\kappa H^2}{8\pi^2}\log a(\tau)).
\end{align}
The one-loop IR effect is to make the gauge coupling time dependent
\begin{align}
-\frac{1}{4g^2(\tau )}F^a_{\mu\nu}F^{a\mu\nu}-D_{\mu}\phi (D^{\mu}\phi)^*+i\bar{\psi }\gamma^{\mu}D_{\mu}\psi.
\end{align}
It is because the gauge field appears just like the derivatives in the covariant derivatives.
It shows the consistency of our result with the gauge invariance.
This conclusion is in agreement with \cite{KitamotoG} given in a background gauge.
It clearly shows that the IR logarithmic effect is gauge invariant and local.

As we have shown, the effective Lagrangian approach universally applies to scalar, spinor and vector fields.
The advantage is its simplicity and robustness of the conclusion.
It clearly shows that IR logarithmic effect is local.
The requirement of Lorentz invariance at sub-horizon scale minimizes the IR logarithmic effects in such a way that it disappears in the free field theories.
Nevertheless it makes couplings to evolve with time  in interacting theories.

\section{Gauge dependence}
\setcounter{equation}{0}

So far, we have worked in a particular gauge.
Although the two-point functions of the tensor mode are gauge invariant, it is not the case for the other modes of the metric.
Nevertheless the contributions from tensor mode alone cannot lead to the Lorentz invariant result.
In other words, Lorentz invariance requires contributions from both the gauge invariant and gauge dependent degrees of freedom.
So there is no Lorentz invariant quantity per se which does not depend on the gauge parameter.
In order to show that the time evolution of the couplings is physical, we need to identify relations among them in which the gauge dependence cancels out.

In order to investigate this question, we consider a more general gauge with a parameter $\beta$:
\begin{align}
\mathcal{L}_\text{GF}&=-\frac{1}{2}a^{2}F^{\beta}_\mu F^{\beta\mu}, \notag\\
F^{\beta}_\mu&=\beta(\partial_\rho h_\mu^{\ \rho}-2\partial_\mu \omega)+\frac{1}{\beta} (2h_\mu^{\ \rho}\partial_\rho\log a+4\omega\partial_\mu\log a).
\label{generalGF}\end{align}

As it turns out, the evolution speed of each coupling depends on a gauge parameter $\beta$.
However the gauge dependence cancels out in the ratio of the evolution speed of the couplings at least for an infinitesimal change of the gauge parameter  $\beta \sim 1$.
To the first order of $\delta =\beta^2-1$,
\begin{align}
&\langle h^{00}h^{00}\rangle = (1-\delta)\langle h^{00}h^{00}\rangle_0, \notag\\
&\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle=(1-\delta)\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle_0,
\label{P-Scaling}\end{align}
where $\langle h^{00}h^{00}\rangle_0$ denotes the correlator in the original gauge $\delta=0$.
To evaluate the effective Lagrangian (\ref{EL}), these propagators are investigated at a coincident point.

From (\ref{P-Scaling}), the time evolutions of the effective couplings are evaluated as
\begin{align}
\lambda_4(\tau)/\lambda_4 = f^{\Delta_4}(\tau),\hspace{1em}\lambda_Y(\tau)/\lambda_Y=f^{\Delta_Y}(\tau),\hspace{1em}g(\tau)/g=f^{\Delta_g}(\tau),
\label{Scaling}\end{align}
\begin{align}
\Delta_4=\frac{21}{4},\hspace{1em}\Delta_Y=\frac{39}{32},\hspace{1em}\Delta_g=\frac{3}{4},
\label{Exponents}\end{align}
where $f(\tau)$ is a gauge dependent decaying function:
\begin{align}
f(\tau)=1-(1-\delta)\frac{\kappa^2 H^2}{4\pi^2}\log a(\tau).
\end{align}
It is no surprise that the time dependence of each effective coupling is gauge dependent.
That is because time is observer-dependent as we may reparametrize it.
Our proposal is that the ratio of the evolution speed of the couplings are the observables in de Sitter spacetime \cite{KitamotoSD,KitamotoG}.
Specifically, we can construct gauge invariant quantities as\footnote{The scaling exponents $\Delta_i$ themselves are gauge dependent.
We can set $f(\tau)$ as $f(\tau)=1-\frac{\kappa^2 H^2}{4\pi^2}\log a(\tau)$, and then they become $\delta$-dependent.}
\begin{align}
\Delta_Y/\Delta_4=\frac{13}{56},\hspace{1em}\Delta_g/\Delta_4=\frac{1}{7}.
\label{Ratios}\end{align}
The gauge invariance of the relative scaling exponents can be interpreted as follows:
It is sensible to pick a particular coupling and use its time evolution as a physical time.
We assigned the role to the coupling of the quartic interaction in (\ref{Ratios}).
In this setting, the relative scaling exponents measure the time evolution of the couplings in terms of a physical time.

Such a relation is reminiscent of the following phenomena in non-equilibrium physics.
In quantum quench when a system is suddenly brought into criticality,
the one-point functions behave as \cite{Calabrese}
\begin{align}
\langle O_i\rangle=\tilde{f}^{\tilde{\Delta}_i}(\tau),
\end{align}
where $\tilde{\Delta}_i$ is the conformal dimension of $O_i$ and $\tilde{f}(\tau)$ is a non-universal decaying function.
The ratios of relaxation times of different operators are universal in $1+1$ dimension.

It is important to investigate whether this idea works to all orders of gauge parameter shift $\delta$.
To the second order of $\delta$, we find in the Appendix
\begin{align}
&\langle h^{00}h^{00}\rangle_{\delta} = (1-\delta-\frac{5}{2}\delta^2)\langle h^{00}h^{00}\rangle_0, \notag\\
&\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle_{\delta}=(1-\delta+\frac{5}{4}\delta^2)\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle_0.
\end{align}
The gauge parameter no longer cancels in the ratio of two independent correlators.
We need a further consideration to circumvent this problem.

Under this circumstance, we may reparametrize the matter field such that
\begin{align}
\phi \rightarrow e^{\gamma\tau\partial_0}\phi,\hspace{1em}
\psi \rightarrow e^{\gamma\tau\partial_0}\psi,\hspace{1em}
A_0^a \rightarrow e^\gamma e^{\gamma\tau\partial_0}A_0^a,\hspace{1em}
A_i^a \rightarrow e^{\gamma\tau\partial_0}A_i^a.
\end{align}
This may be interpreted as a reparametrization of the conformal time:
\begin{align}
\tau \rightarrow e^{\gamma} \tau.
\label{fieldrepara}\end{align}
Using this freedom, we may modify the way $h^{00}$, $\omega$ couple to the matter fields as\footnote{This kind of counter term is necessary to renormalize UV divergences in a non-covariant gauge \cite{Kitamoto2013}.}
\begin{align}
h^{00} \rightarrow h^{00} +\frac{3}{2}\gamma,
\hspace{1em} \omega \rightarrow \omega+\frac{3}{4}\gamma.
\label{metrictrans}
\end{align}
Since $\partial_{\mu}(e^\gamma \tau )=e^\gamma\delta_\mu^{\ 0}+O(k\tau)$, it is a small transformation for super-horizon mode.
The same cannot be said for a reparametrization of spatial coordinates $x^i\rightarrow (e^h)^i_{\ j} x^j$.
Here $\partial_\mu((e^h)^i_{\ j} x^j)=(e^h)^i_{\ j}\delta_\mu^{\ j}+O(k x^i)$, the $O(kx^i)$ term cannot be neglected.
Thus we cannot transform tensor modes in an analogous way.

We can choose $\gamma=\frac{5}{4}\delta^2h^{00}$ in such a way that the gauge dependence cancels in the ration of two independent correlators: $h^{00}_\text{eff}=2\omega_\text{eff}=(1+\frac{15}{8}\delta^2)h^{00}$ and $\tilde{h}^{ij}$.
We also need to modify the rescaling of the fields as
\begin{align}
\phi &\rightarrow a^{-1}e^{(\alpha-1)\omega_\text{eff}}\phi, \notag\\
\psi &\rightarrow a^{-\frac{3}{2}}e^{(\beta-\frac{3}{2})\omega_\text{eff}}\psi,  \notag\\
A_{\mu} &\rightarrow A_{\mu}.
\label{rescale'}\end{align}
As far as the IR logarithmic effects to the matter fields are concerned, the Lorentz invariance is preserved at sub-horizon scale in this procedure.
Furthermore the scaling relation (\ref{Scaling})-(\ref{Exponents}) holds to the second order of gauge parameter change $\delta$ with $f(\tau)=1-(1-\delta+\frac{5}{4}\delta^2)\frac{\kappa^2 H^2}{4\pi^2}\log a(\tau)$.
Since we have one parameter freedom $\gamma$ for the gauge parameter change $\delta$, we can repeat the same procedure to all orders of $\delta$.
This argument provides a further evidence that the scaling relation (\ref{Scaling})-(\ref{Exponents}) is observable in de Sitter spacetime.

The essential point here is that there are two independent correlators of metric: scalar and spin-2 modes.
The gauge dependence of the relative magnitude can be canceled by a reparametrization of the matter field.
Such a procedure is required by the Lorentz invariance of matter dynamics at sub-horizon scale.
So Lorentz invariance leads to gauge independence of the observables.

Let us consider the graviton propagator in a different class of gauge: namely covariant de Donder gauge \cite{deDondergauge}.
In this gauge the two-point function of spin-2 mode $\tilde{h}^{ij}$ is identical to that in the original gauge with $\delta=0$
\begin{align}
\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle_\text{de Donder}=\langle\tilde{h}^{ij}\tilde{h}^{kl}\rangle_0.
\end{align}
However the remaining (scalar) component of the propagator in de Donder gauge does not agree with that in the original gauge $\delta=0$.
In fact it exhibits more singular IR behavior.
We focus on the traceless part of the propagator as it is relevant to possible Lorentz symmetry breaking.
The trace part decouples when the matter field is conformally coupled.
The important point is that the IR singular part of this component is pure gauge.
We can effectively eliminate that part of the graviton propagator by a judicious reparametrization of the matter field.
We subsequently put $h^{00}_\text{eff}$ identical to that of the original gauge by using the transformation (\ref{fieldrepara}) with $\gamma=\frac{2}{3}h^{00}$
\begin{align}
\langle h^{00}_\text{eff}h^{00}_\text{eff}\rangle_\text{de Donder}=\langle h^{00}h^{00}\rangle_0.
\end{align}
Recall that the correlators in the original gauge are singled out by requiring Lorentz invariance of matter dynamics at the sub-horizon scale.
We thus conclude that the scaling relation in (\ref{Scaling}) also holds in de Donder gauge once we impose effective Lorentz symmetry at sub-horizon scale.

\section{Conclusion}
\setcounter{equation}{0}

We construct effective Lagrangian in de Sitter spacetime which summarizes IR logarithmic effects on local matter field dynamics.
We integrate super-horizon mode of the metric fluctuations.
We require Lorentz symmetry at sub-horizon scale.
This strategy utilizes rescaling ambiguity of the matter fields with respect the conformal mode.

The Lorentz symmetry at sub-horizon scale is a natural consequence of general covariance.
We argue that IR logarithmic effects cannot spoil this fundamental symmetry.
From microscopic physics point of view, general covariance is required from Lorentz symmetry and unitarity for quantum gravity.
We thus argue that our requirement must be equivalent to the unitarity requirement.

In de Sitter space, Lorentz symmetry is violated at horizon scale.
However Lorentz symmetry holds at the sub-horizon scale in accordance with general covariance.

The important consequence of our strategy is that we obtain unique physical predictions.
In particular, we find that the couplings of $\phi^4$, Yukawa and gauge interactions are decaying with time.
The screenings of these couplings can be expressed as logarithmic dependences of the scale factor and then they are slowly-progressing with cosmic evolution.
Furthermore the relative evolution speed of the couplings are gauge independent once we impose Lorentz invariance of matter dynamics at the sub-horizon scale.

The gauge change of metric propagator is equivalent to the coordinate change of the observer.
However the observer can choose a new coordinate in which Lorentz symmetry holds at sub-horizon scale.
This is realized by a reparametrization of matter fields.
This procedure cancels the gauge change of the metric propagator as the matter-metric interaction is invariant under the reparametrization of both metric and matter fields.
It is clear that a local quantity in quantum gravity depends on the observer.
We need to specify the coordinate system of the observer to obtain gauge independent quantity.
The requirement of Lorentz symmetry specifies such a system.

\section*{Acknowledgment}

This work is supported by Grant-in-Aid for Scientific Research (B) No. 26287044 and (C) No. 16K05336.
We thank S. Iso, H. Kawai and T. Tanaka for discussions.