\section{Introduction}

Let $(Y,\omega)$ be a compact K\"ahler manifold, and $d\mu$ be a measure on $Y$. We say $(Y,\omega, d\mu)$ satisfies the \emph{Skoda type inequality}, if for any K\"ahler potential $u \in PSH(Y,\omega)$ normalised to $\sup u=0$,
\begin{equation}\label{Skodaassumption}
\int_Y e^{-\alpha u} d\mu \leq A,
\end{equation}
where $\alpha, A$ are independent of  $u$. A prototype theorem is

\begin{thm}\label{Skodastandardversion}\cite{Tian}
	On a fixed compact K\"ahler $(Y, \omega)$, the Skoda type inequality holds for $d\mu=\omega^n$.
\end{thm}

\begin{rmk}
Here the supremum of all such $\alpha$
 is known as Tian's alpha invariant, important for existence questions of K\"ahler-Einstein metrics.
\end{rmk}

We are interested  in keeping track of these constants $\alpha, A$ as  $(Y,\omega,d\mu)$ varies. The main theme of this paper is that oftentimes the Skoda constants can be chosen uniformly for quite flexible choices of probability measures $d\mu$, even when the complex structure degenerates severely. 	In the literature $\alpha$ is much studied (\cf \cite{Tian}\cite{GZ}), and a very recent preprint \cite{Eleonora} made aware to the author after the completion of this work contains a  uniform estimate for both $\alpha, A$ in the related context of K\"ahler-Einstein manifolds.

Our main application is to algebraic degenerations of Calabi-Yau manifolds.
We work over $\C$. Let $S$ be a smooth affine algebraic curve, with a point $0\in S$.
An \emph{algebraic degeneration family} is given by a submersive projective morphism $\pi: X\to S\setminus \{0\}$ with smooth connected $n$-dimensional fibres $X_t$ for $t\in S\setminus \{0\}$. A \emph{polarisation} is given by an ample line bundle $L$ over $X$; the sections of a sufficiently high power of
$L$ induces an embedding $X\to \mathbb{CP}^N$, hence a Fubini-Study metric $\omega_X$ on $(X, L)$. For $0<|t|\ll 1$,
a fixed choice of $\omega_X$ induces \emph{rescaled} background metrics $\omega_t= \frac{1}{ |\log |t| |  } \omega_X|_{ X_t }$ on $X_t$ in the class
$\frac{1}{|\log |t| |   }c_1(L)  $.

 A \emph{model} of $X$ is a normal flat projective $S$-scheme $\mathcal{X}$ which agrees with $\pi: X\to S\setminus \{0\}$ over the punctured curve. It is called a \emph{semistable snc model} if $\mathcal{X}$ is smooth,
 the central fibre over $0\in S$ is \emph{reduced} and is a simple normal crossing divisor in $\mathcal{X}$. By the semistable reduction theorem \cite[chapter 2]{ToroidalembeddingsI}, after finite base change to another smooth algebraic curve $S'$, we can always find some semistable snc model for the degeneration family $X\times_S (S'\setminus \{0\})$. Everything here is quasi-projective.

 We say the degeneration family is \emph{Calabi-Yau} if there is a trivialising section $\Omega$ of the canonical bundle $K_X$. Over a small disc $\mathbb{D}_t$ around $0\in S$, this induces holomorphic volume forms $\Omega_t$ on $X_t$ via $\Omega= dt\wedge \Omega_t$. The \emph{normalised Calabi-Yau measure} on $X_t$ is the probability measure
 \begin{equation}\label{CalabiYaumeasureeqn}
 d\mu_t= \frac{   \Omega_t \wedge \overline{\Omega}_t }{   \int_{X_t} \Omega_t \wedge \overline{\Omega}_t   }.
 \end{equation}

Our main result is

\begin{thm}(Uniform Skoda estimate)\label{UniformSkodamaintheorem}
	Given a polarised algebraic Calabi-Yau degeneration family $\pi: X\to S\setminus \{0\}$ as above. Then there are uniform positive constants $\alpha, A$ independent of $t$ for $0<|t|\ll 1$, such that for the normalised Calabi-Yau measures $d\mu_t$,
	\[
	\int_{X_t} e^{-\alpha u}  d\mu_t \leq A, \quad \forall u\in PSH(X_t,\omega_t) \text{ with } \sup_{X_t} u=0.
	\]

\end{thm}

This is proved by reducing to the semistable snc model case, and prove a general Skoda type estimate there (\cf Theorem \ref{UniformSkodasemistable}). A major consequence, readily reaped using Kolodziej's estimate (\cf Theorem \ref{pluripotentialthm1}), is

\begin{thm}\label{UniformLinfty}
(Uniform $L^\infty$-estimate) Let $\phi_t$ be the K\"ahler potential of the Calabi-Yau metric in the class $(X_t, [\omega_t])$, namely
\[
\frac{(\omega_t+ \sqrt{-1} \partial\bar{\partial}\phi)^n }{ \int_{X_t} \omega_t^n  }= d\mu_t,\quad \sup_{X_t} \phi_t=0.
\]
Then $\norm{ \phi_t}_{L^\infty} \leq C$ independent of $t$ for $0<|t|\ll 1$.

\end{thm}

\begin{rmk}
Applications of pluripotential theory to Calabi-Yau metrics when the K\"ahler class is degenerating can be found in \cite{EGZ2}, which is used further in \cite{Tosatti}. Our main results generalize certain aspects of \cite{LiSYZ} which focuses on  degenerating projective hypersurfaces near the large complex structure limit.
\end{rmk}

\begin{Acknowledgement}
The author is a 2020 Clay Research Fellow, currently based at the Institute for Advanced Study. He thanks Song Sun and Simon Donaldson for discussions, and Sebastien Boucksom, Eleonora Di Nezza, and Valentino Tossati for comments.
\end{Acknowledgement}

\section{Uniform Skoda inequality}\label{UniformSkodachapter}

We work in the context of semistable simple normal crossing (snc) models. Concretely,
let $\pi: \mathcal{X}\to \mathbb{D}_t$ be a flat projective family of $n$-dimensional varieties over a small disc $\mathbb{D}_t$, such that the total space $\mathcal{X}$ is smooth, $\pi$ is a submersion over the punctured disc with connected fibres, the central fibre $X_0$ is \emph{reduced} and is an snc divisor in $\mathcal{X}$. Denote the components of $X_0$ as $E_i$ with $i\in I$. We equip $\mathcal{X}$ with a fixed background K\"ahler metric $\omega_{\mathcal{X}}$, inducing a distance function $d_{\omega_{\mathcal{X}}}$. This induces a family of \emph{rescaled} K\"ahler metrics $\omega_t= \frac{1}{ |\log |t||  } \omega_{ \mathcal{X}  }|_{X_t  }$. We shall derive a uniform Skoda type estimate (\ref{Skodaassumption}) for $(X_t, \omega_t, d\mu_t)$, where $d\mu_t$ belongs to a natural class of measures. The main result is Theorem \ref{UniformSkodasemistable}.

\subsection{Quantitative stratification and good test functions}

There is a \emph{quantitative stratification} on any smooth fibre $X_t$ induced by the intersection pattern of $E_i$: for $J\subset I$ such that $E_J=\cap_{i\in J} E_i\neq \emptyset$, the corresponding statum is
\[
E_J^0= \{   x\in X_t | d_{\omega_{\mathcal{X}}}( x, E_J )\lesssim \epsilon     \} \setminus \{    x\in X_t | d_{\omega_{\mathcal{X}}}( x, E_{J'} )\lesssim \epsilon ,  \text{ some } J'\supsetneq J     \},
\]
namely a small `$\epsilon$-tubular neighbourhood' of $E_J$ minus the deeper strata. For $J=\{ i\}$ we write $E_i^0= E_{\{i\}  }^0$. Here the disc $\mathbb{D}_t$ and the small parameter $\epsilon\ll 1$ can be shrinked for convenience; the essential thing is that all parameters should be independent of the coordinate $t$.

It is useful to introduce local coordinates $\{ z_i \}_0^{n}$ around $E_J\subset \mathcal{X}$, such that $z_0, \ldots, z_p$ with $p=|J|-1$ are the local defining equations of $E_j$ for $j\in J$, and locally the  fibration map is $t=z_0\ldots z_p$.  Then up to uniform equivalence, locally
\begin{equation*
\omega_{ \mathcal{X} }\sim \sum_0^{n} \sqrt{-1}dz_i\wedge d\bar{z}_i.
\end{equation*}

The rest of this section is devoted to the construction of good test functions.
Given any of these divisors $E_0$, we can find a nonnegative function $h=h_{E_0}$ on $\mathcal{X}$, such that
\begin{itemize}
\item

In the local charts near $E_0$ with $z_0$ being the defining function for $E_0$,
\[
h= |z_0|^2 \tilde{h}(z_0, \ldots z_n)
\]
for some positive smooth function $\tilde{h}$;
\item Away from $E_0$ the function $h$ is comparable to 1.
\end{itemize}

We observe
\begin{itemize}
	\item
The form $\partial \bar{\partial} \log h= \partial \bar{\partial} \log \tilde{h}$ extends smoothly;
\item
For $|t|^2\ll h\lesssim \delta\ll 1$ inside $X_t$, so that $|z_0|\gg |t|$, by a local calculation  near $E_J$ with $0\in J$,
\[
\begin{split}
 &\sqrt{-1}   \partial  \log h   \wedge \bar{\partial } \log h   \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1} \geq \frac{ \sqrt{-1} }{2|z_0|^2} dz_0\wedge d\bar{z}_0\wedge\omega_{ \mathcal{X}  }|_{X_t}^{n-1}
 \\
 &\gtrsim   \min\{  \frac{1}{|z_0|^2}, \max_{1\leq i\leq p} |z_i|^{-2}  \} \omega_{ \mathcal{X}  }|_{X_t}^n
 \\
 &
 \gtrsim
  \min\{  \frac{1}{h}, h^{1/p} |t|^{-2/p}  \}\omega_{ \mathcal{X}  }|_{X_t}^n
  \\
  &
 \gtrsim
 \min\{  \frac{1}{h}, h^{1/n} |t|^{-2/n}  \}\omega_{ \mathcal{X}  }|_{X_t}^n .
\end{split}
\]
Here in the first line we need to fix $\delta\ll 1$ so that the effect of $\partial \log h$ is dominated by $d\log z_0$.  The second line uses that for $1\leq k\leq p$, the volume forms on $X_t$
\[
\frac{1}{|z_0|^2} dz_0\wedge d\bar{z}_0\wedge
\prod_{j\neq k, 1\leq j\leq n} \sqrt{-1}dz_j\wedge d\bar{z}_j
\sim
\frac{1}{|z_k|^2} \prod_{ 1\leq j\leq n} \sqrt{-1}dz_j\wedge d\bar{z}_j,
\]
and the third line uses $|t|=|z_0\ldots z_p|\sim h^{1/2}|z_1\ldots z_p|$.

\item
On $X_t$ the function $h\gtrsim |t|^2$.
The region $\{ |t|^2\sim h  \}\subset X_t$ can be identified as $E_0^0$, namely the vicinity of $E_0$ away from deeper strata. Here
\[
\begin{split}
&\sqrt{-1}   \partial  \log h   \wedge \bar{\partial } \log h   \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1} \geq 0,
\\
&
\sqrt{-1}   \partial  \bar{\partial } \log h   \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1} \gtrsim -\omega_{ \mathcal{X}  }|_{X_t}^{n}  .
\end{split}
\]

\end{itemize}

\begin{lem}\label{goodtestfunction}
	(Good test function) Given the divisor $E_0$, we can choose a $C^2$ test function $v$ on $X_t$ such that the following hold uniformly for small $t\neq 0$:
	\begin{itemize}
		\item $v$ is zero for $h\geq \delta$.

		\item Globally $0\leq v\leq -\log |t|$.

		\item
		For any divisor $E_j$ intersecting $E_0$, there is a subset of $E_j^0$ with measure at least $C_2$ on which $\sqrt{-1} \partial \bar{\partial}v \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1}  \geq C_3\omega_{ \mathcal{X}  }|_{X_t}^n $.

		\item
		For $C_4 |t|^2 \leq h\leq \delta$, the form $\sqrt{-1}\partial \bar{\partial} v\wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1} \geq 0$.

		\item
		For $h\leq C_4 |t|^2$, the form $\sqrt{-1}\partial \bar{\partial} v \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1}\geq -C_5\omega_{ \mathcal{X}  }|_{X_t}^n $.

	\end{itemize}
\end{lem}

\begin{proof}
We seek the test function in the form $v=\Phi\circ \log h$ for some convex, non-increasing, non-negative $C^2$-function $\Phi$. Compute
\[
\partial \bar{\partial} v=  \Phi'' \partial  \log h   \wedge \bar{\partial } \log h   + \Phi' ( \partial \bar{\partial } \log \tilde{h}),
\]
so using the properties of $h$ above,
\[
\sqrt{-1} \partial \bar{\partial}v \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1}
\geq
\begin{cases}
\left( \Phi'' C_1' \min\{  \frac{1}{h}, h^{1/n} |t|^{-2/n}  \} +\Phi'C_2'  \right)    \omega_{ \mathcal{X}  }|_{X_t}^n     , \quad &|t|^2\lesssim h\leq \delta,
\\
C_3' \Phi' \omega_{ \mathcal{X}  }|_{X_t}^n, \quad & h\lesssim |t|^2.
\end{cases}
\]

To satisfy our conditions on $v$, it is enough to have
\begin{itemize}
\item
$\Phi(x)=0$ for $x\geq \log \delta$.
\item
$|\Phi'(x)|\lesssim 1$ for $2\log |t|\lesssim x\leq \log \delta$.
\item
$ -\frac{d}{dx} \log |\Phi'|=  \frac{\Phi''}{|\Phi'|  }  \geq C_4' \max\{ h, h^{-1/n} |t|^{2/n}  \}$ for $ h=e^x\leq \delta,$ where
$C_4'> C_2'/C_1'$. Morever, for $x<\delta$, we need $\Phi'<0$ so that
$\sqrt{-1} \partial \bar{\partial}v \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1}$ has some strict positivity for $\delta/2<h<\delta$. Notice convexity of $\Phi$ is a consequence of these conditions.

\end{itemize}

To construct such $\Phi$, we can prescribe the behaviour near $x=\log \delta$ by
$
\Phi'(x)= -e^{1/(x-\log \delta  )}
$ for $x< \log \delta$,
and match this with a solution to
\[
-\frac{d}{dx} \log |\Phi'|  = C_4' \max\{e^x, e^{-x/n} |t|^{2/n}  \}, \quad x<\log \delta
\]
for some large enough $C_4'$, such that $\Phi'$ remains $C^1$ at the matching point. Integration shows that $ |\Phi'|$ remains uniformly bounded at $h\sim |t|^2$, or equivalently $x\sim 2\log |t|$.
\end{proof}

\subsection{Convexity}

Consider $u\in PSH(X_t, \omega_t)$ normalised to $\sup_{X_t} u=0$. Equivalently, we can cover $X_t$ by a bounded number of charts as before, and use the
 local potentials of $\omega_{\mathcal{X} }$ to represent $u$  as a collection of local plurisubharmonic (psh) functions $\{ u_\beta\}$ with $|u_\beta-u|\leq C$.

\begin{lem}(Convexity)\label{pshconvexity}
	Let $\phi$ be any psh function on the open subset of $  \{   1<|z_i|<\Lambda, i=1,\ldots p, |z_k|<1, k=p+1, \ldots n \}\subset (\C^*)^p\times \C^{n-p}$. Then the function \[
	\bar{\phi}(x_1, \ldots x_n)= \frac{1}{(2\pi)^n} \int_{D(1)^{n-p} } \prod_{p+1}^n \sqrt{-1} dz_k \wedge d\bar{z}_k
		\int_{T^p} \phi( e^{x_1+i\theta_1}, \ldots e^{x_p+i\theta_p}) d\theta_1\ldots d\theta_p
	\]
	is convex.
\end{lem}

\begin{proof}
For any choice of $\theta_i$  the function  $\phi(z_1 e^{i\theta_1}, \ldots z_p e^{i\theta_n}, z_{p+1},\ldots, z_n)$ is psh, since the $T^p$-action on $(\C^*)^p$ is holomorphic. Thus the average function $\bar{\phi}$ is also psh as a function of $z_1,\ldots z_p$. Any $T^p$-invariant psh function must be convex in the log coordinates, because for $x_i=\log |z_i|$,
	\[
	\sqrt{-1}\partial \bar{\partial} \bar{\phi}= \frac{1}{4} \sum \frac{\partial^2 \bar{\phi}}{  \partial x_i \partial x_j  } \sqrt{-1} d\log z_i \wedge d\overline{\log z_j} \geq 0.
	\]
\end{proof}

\subsection{Harnack type inequality}

\begin{lem}
	(Almost maximum on top strata)
	For $u\in PSH(X_t, \omega_t)$ normalised to $\sup_{X_t} u=0$, there is some $i\in I$, such that
	\[
	\sup_{E_i^0} u \geq -C, \quad \int_{E_i^0} u \omega_{ \mathcal{X} }|_{X_t}^n \geq -C'.
	\]

\end{lem}

\begin{proof}
	Let the global maximum of $u$ be achieved at $q_0\in E_J^0$, and denote the local potential of $u$ as $u_\beta$. Without loss of generality $u_\beta\leq 0$. We have $u_\beta(q_0) \geq -C$ since $|u-u_\beta|\leq C$. Applying the mean value inequality around $q_0$, we find that the local average function $\bar{u}_\beta$ produced in Lemma \ref{pshconvexity} satisfies $\sup \bar{u}_\beta\geq -C$ for another uniform constant $C$. By the convexity of $\bar{u}_\beta$ its sup is almost achieved at the boundary of the chart, which is contained in a union of less deep strata $E_{J'}^0$ with $J'\subsetneq J$. Thus we can find a point $q'$ with $u(q')\geq -C$ that belongs to a less deep stratum; an induction shows that there is some $i\in I$, such that $\sup _{E_i^0} u \geq -C$.

For the $L^1$-bound we recall the following Harnack inequality argument. Suppose a coordinate ball $B(q, 3R)$ is contained in a local chart in a small neighbourhood of $E_i^0$.
Applying the mean value inequality to the local psh function associated to $u$, we see for $y\in B(q,R)$ that
\[
u(y) \leq  C+ \dashint_{B(y,2R)} u \lesssim 1+ \dashint_{B(q,R)} u.
\]
hence the Harnack inequality
\[
\dashint_{B(q,R)} |u |\lesssim 1+ \inf_{B(q, R)} (-u).
\]
Applying this to a chain of balls connecting any two points in $E_i^0$ gives the $L^1$-bound $\int_{E_i^0} u \omega_{ \mathcal{X} }|_{X_t}^n \geq -C'$; the bound is uniform because the number of balls involved in the chain can be controlled independent of $t$.
\end{proof}

\begin{prop}(Almost maximum on top strata II)\label{AlmostmaximumII}
	There is a uniform lower bound for all $|t|\ll 1$ and all $i\in I$:
\begin{equation}
\sup_{E_i^0} u \geq -C, \quad \int_{E_i^0} u \omega_{ \mathcal{X} }|_{X_t}^n \geq -C'.
\end{equation}
\end{prop}

\begin{proof}
The $L^1$-estimate follows from the sup estimate as above, so the real problem is to transfer bounds between different $E_i^0$. This is nontrivial because the necks connecting $E_i^0$ with each other are highly degenerate.

Given one divisor $E_0$ such that $
 \int_{E_0^0} u \omega_{ \mathcal{X} }|_{X_t}^n \geq -C,
$
we produce a good test function $v$ by Lemma \ref{goodtestfunction}. Integrating by parts,
\[
\int_{X_t} v \sqrt{-1}\partial \bar{\partial} u \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1} = \int_{X_t} u \sqrt{-1}\partial \bar{\partial} v \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}.
\]
The LHS is the difference of $\int_{X_t} v(\omega_t+ \sqrt{-1}\partial \bar{\partial} u) \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}$ and $\int_{X_t} v\omega_t \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}$, and since $-\log |t|\gtrsim v\geq 0$ both terms are bounded between $0$ and $C$. Thus
\[
|\int_{X_t} u \sqrt{-1}\partial \bar{\partial} v \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}|\leq C.
\]

Now the form
$\sqrt{-1}\partial \bar{\partial} v \wedge \omega_{ \mathcal{X}  }|_{X_t}^{n-1}$ can only be negative on $\{ h\sim |t|^2  \}=E_0^0$, and is bounded below by $-C  \omega_{ \mathcal{X}  }|_{X_t}^{n}$. Thus the positive part of the signed measure $u \sqrt{-1}\partial \bar{\partial} v \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}$ has total mass controlled by $ \int_{E_0^0} |u| \omega_{ \mathcal{X} }|_{X_t}^n \leq C$.
Consequently, the negative part of the signed measure must also have total mass $\leq C$.

By construction, for any divisor $E_j$ intersecting $E_0$ there is a nontrivial amount of $\sqrt{-1}\partial \bar{\partial} v \wedge \omega_{\mathcal{X} }|_{X_t}^{n-1}$-measure inside $E_j^0$. This forces $\sup_{E_j^0} u \geq -C$. To summarize, we have transferred the sup bound from $E_0^0$ to any $E_j^0$ with $E_j\cap E_0\neq \emptyset$. Since the central fibre $X_0$ is connected, in at most $|I|$ steps
this sup bound is transferred to all $E_i^0$ with $i\in I$.
\end{proof}

\begin{rmk}
This proof is inspired by the intersection theoretic argument of \cite[section 6.1]{Boucksomsemipositive}, which can be viewed as a non-archimedean analogue.
\end{rmk}

\subsection{Local $L^1$ estimate}

Given a local chart on $E_J^0$ with $\C^*$-coordinates $z_1,\ldots z_p$ and $\C$-coordinates $z_{p+1}, \ldots, z_n$, and a point $q$ therein, we shall refer to the subregion
\[
	\{  \frac{1}{2} |z_i(q)|  \lesssim|z_i| \lesssim 2|z_i(q)|  , \quad 1\leq i\leq p    \}
\]
as a \emph{log scale}.

\begin{lem}\label{LocalL1estimate}
	(Local $L^1$-estimate)
Within every log scale there is a uniform bound on the $L^1$-average integral
\[
\dashint_{loc} |u| \prod_1^p \sqrt{-1}d\log z_i \wedge d\log \bar{z}_i \wedge \prod_{p+1}^n \sqrt{-1}d z_k \wedge d \bar{z}_k \leq C.
\]
\end{lem}

\begin{proof}
We induct on the depth of the strata. For $p=0$ this follows from Prop. \ref{AlmostmaximumII}. So let us assume the bound is achieved for depth $<p$. For a given chart, we consider the local psh function $u_\beta$ associated to $u$ and produce the convex average function $\bar{u}_\beta$ as in Lemma \ref{pshconvexity}. Since a definite neighbourhood of the boundary of the chart lies inside less deep strata, we know that near the boundary $|\bar{u}_\beta|\leq C$ by the induction hypothesis and the convexity condition. Using convexity again in the interior of the chart we see $|\bar{u}_\beta|\leq C$ in the whole chart.

Within any log scale, by construction the local average
$
\dashint_{loc} (u_\beta- \bar{u}_\beta)=0.
$
But $u_\beta\leq C$ by $u\leq 0$, hence
\[\dashint_{loc} |u_\beta-\bar{u}_\beta| \lesssim \dashint_{loc} (u_\beta- \bar{u}_\beta )_+ \leq C.
\]
Using $|u-u_\beta|\leq C$ we conclude the local $L^1$-estimate on $u$.
\end{proof}

\subsection{Local Skoda  estimate}\label{Skodainequalitybackgroundsection}

We recall a basic version of the \emph{Skoda inequality}:

\begin{prop}\label{Skodabasicversion}
	(\cf \cite[Thm 3.1]{Zeriahi}) If $\phi$ is psh on $B_2\subset \C^n$, with $\int_{B_2} |\phi| \omega_E^n \leq 1$ with respect to the standard Euclidean metric $\omega_E$, then
	there are dimensional constants $\alpha$, $C$, such that
	\[
 \int_{B_1} e^{-\alpha \phi} \omega_E^n \leq  C.
	\]

\end{prop}

Applying this with Lemma \ref{LocalL1estimate},

\begin{cor}
(Local Skoda estimate) Within every log scale, there are uniform positive constants $\alpha$ and $C$, such that
\[
\dashint_{loc} e^{-\alpha u} \prod_1^p \sqrt{-1}d\log z_i \wedge d\log \bar{z}_i \wedge \prod_{p+1}^n \sqrt{-1}d z_k \wedge d \bar{z}_k \leq C.
\]
\end{cor}

\subsection{Uniform global Skoda estimate}

We are interested in the following class of measures, motivated by Calabi-Yau measures (\cf section \ref{CalabiYaumeasure}). Let $a_i$ be non-negative real numbers assigned to $i\in I$, with $\min a_i=0$. Let
\[
m= \max\{ |J|-1:  E_J\neq \emptyset, a_i=0 \text{ for } i\in J     \}.
\]
We say the measures $d\mu_t$ on $X_t$ satisfy a \emph{uniform upper bound of class $(a_i)$}, if on the local charts of each $E_J^0$,
\begin{equation}\label{measureupperbound}
d\mu_t \leq \frac{ C}{ |\log |t||^m  } |z_0|^{2a_0}\cdots |z_p|^{2a_p} \prod_1^p  \sqrt{-1}d\log z_i \wedge d\log \bar{z}_i \wedge \prod_{p+1}^n \sqrt{-1}d z_k \wedge d \bar{z}_k.
\end{equation}
 The normalisation factor ensures $\int_{X_t}d\mu_t \leq C$ independent of $t$, by a straightforward local calculation.

\begin{thm}\label{UniformSkodasemistable}
(Uniform Skoda estimate) Suppose the measures $d\mu_t$ on $X_t$ satisfy a uniform upper bound of class $(a_i)$. Then there are uniform positive constants $\alpha$ and $A$, such that
\[
\int_{X_t} e^{-\alpha u}  d\mu_t \leq A, \quad \forall u\in PSH(X_t,\omega_t) \text{ with } \sup_{X_t} u=0.
\]

\end{thm}

\begin{proof}
We choose the charts so that each point on $X_t$ is covered by $\leq C$ log scales. Summing over the local Skoda estimates from all log scales, $\int_{X_t} e^{-\alpha u}d\mu_t$ is bounded by
\[
\begin{split}
 &\frac{C}{ |\log |t||^m  }  \sum_{\text{log scales}} \int_{loc} |z_0|^{2a_0}\ldots |z_p|^{2a_p}
\prod_1^p \sqrt{-1}d\log z_i \wedge d\log \bar{z}_i \wedge \prod_{p+1}^n \sqrt{-1}d z_k \wedge d \bar{z}_k
\\
&\leq C.
\end{split}
\]
\end{proof}

\section{Application to Calabi-Yau degeneration}

We work in the setting of polarised algebraic degeneration of Calabi-Yau manifolds, as in the Introduction.

\subsection{Calabi-Yau measure}\label{CalabiYaumeasure}

The Calabi-Yau measure (\ref{CalabiYaumeasureeqn}) is studied thoroughly in \cite{Boucksom1}, but it is illustrative to recall it explicitly on a semistable snc model $\mathcal{X}$. The discussion is local on the base, and we will follow the notations of section \ref{UniformSkodachapter}, e.g. the components of the central fibre are denoted as $E_i$ for $i\in I$.

The canonical divisor $K_{ \mathcal{X} }=\sum_i a_i E_i$ is supported on the central fibre, since $K_X$ is trivialised. Multiplying $\Omega$ by a power of $t$, which does not change $d\mu_t$, we may assume $\min a_i=0$.
In the local coordinates around $E_J$ away from the deeper strata,
\[
\Omega= f_J \prod_0^p z_i^{a_i} dz_i \wedge \prod_{p+1}^n dz_j
\]
for some nowhere vanishing local holomorphic function $f_J$. Since $t= z_0\ldots z_p$,
\[
\Omega_t=  f_J z_0^{a_0}\ldots z_p^{a_p} \prod_1^p  d\log z_i \wedge \prod_{p+1}^n dz_j,
\]
hence
\[
\sqrt{-1}^{n^2}\Omega_t \wedge \overline{\Omega}_t= |f_J|^2 |z_0|^{2a_0}\ldots |z_p|^{2a_p} \prod_1^p \sqrt{-1}d\log z_i \wedge d\overline{\log z_i} \wedge \prod_{p+1} ^n\sqrt{-1}d z_j \wedge d\overline{z}_j.
\]
The total measure $\int_{X_t} \sqrt{-1}^{n^2}\Omega_t \wedge \overline{\Omega}_t$ is of the order $O(|\log |t| |^m  )$ where
\[
m= \max\{ |J|-1:  E_J\neq \emptyset, a_i=0 \text{ for } i\in J     \}.
\]
Thus $d\mu_t $ satisfies a uniform upper bound of class $(a_i)$ (\cf (\ref{measureupperbound})).

\subsection{Uniform Skoda estimate}

We now prove the main theorem \ref{UniformSkodamaintheorem}.

\begin{proof}
First we observe that the choice of the Fubini-Study metric $\omega_X$ is immaterial. Given any two choices, the relative K\"ahler potential between them is bounded by $O(|\log |t| |)$ for $0<|t|\ll 1$, because the pole order of a section near $t=0$ must be finite. Thus the relative K\"ahler potential between two choices of $\omega_t$ is bounded by $O(1)$ independent of $t$, which affects the Skoda constant $A$ but not its uniform nature.

We now pass to a finite base change and find a semistable reduction. The Calabi-Yau measure $d\mu_t$ on $X_t$ is independent of the parametrisation of the base, and is preserved under finite base change. Thus it is enough to prove it assuming $\omega_X$ agrees with a smooth K\"ahler metric on a semistable snc model $\mathcal{X}$; this is a special case of Theorem \ref{UniformSkodasemistable}.
\end{proof}

\subsection{Uniform  $L^\infty$-estimate}\label{Toolsfrompsh}

We recall the following result proved using Kolodziej's pluripotential theoretic methods (\cf \cite[section 2.2]{LiSYZ} for an exposition based on \cite{EGZ}\cite{EGZ2}):

\begin{thm}\label{pluripotentialthm1}
	Let $(Y, \omega)$ be a compact K\"ahler manifold, and the K\"ahler potential $\phi$ solves the complex Monge-Amp\`ere equation
	\[
\frac{	(\omega+\sqrt{-1} \partial \bar{\partial} \phi )^n}{  \int_Y \omega^n }=d\mu, \quad \sup \phi=0.
	\]
 Assume there are positive constants $\alpha, A$, such that the Skoda type estimate (\ref{Skodaassumption}) holds
 for $(Y,\omega, d\mu)$:
	\begin{equation*}
	\int_Y e^{-\alpha u}  d\mu \leq A, \quad \forall u\in PSH(Y,\omega) \text{ with } \sup_Y u=0.
	\end{equation*}
Then  $\norm{\phi}_{C^0} \leq C(n,\alpha,A)$.

\end{thm}

The uniform $L^\infty$-estimate for the Calabi-Yau potentials in Theorem \ref{UniformLinfty} is an immediate consequence.


\section{Introduction}
The COVID-19 pandemic \emph{started} in Switzerland on February 25\textsuperscript{th} 2020, when the first infection was officially reported in the Italian-speaking canton of Ticino~\cite{FirstCase:EN,FirstCase:DE}. Soon the pandemic spread around the country on all cantons, and Switzerland became one of the most infected countries on a per-capita basis~\cite{Vaud:DE,Spread:EN}.

The Swiss government started putting in place various measures to control and suppress the pandemic. Gatherings were limited and later totally banned, following by closure of all except essential business, and finally closing land borders with neighbouring countries~\cite{Measures:DE}.

These measures helped control the spread of the virus and significantly decreased the number of active and daily new cases in Switzerland. With the success confirmed, government started gradually lifting the established restrictions from late April~\cite{Lifting}. Finally, the June 15\textsuperscript{th} re-opening land borders with the neighbouring countries marked the \emph{"end"} of the pandemic in Switzerland, at least for the first wave~\cite{BorderOpening}.

Since the first recorded case of the COVID-19 virus in Switzerland and far before as it was gaining attention around the world, The Swiss media started covering the topic from various aspects, including the everyday news about the state of the country, the immediate effects, and longer-term consequences of the pandemic. Given the multi-lingual and multi-cultural nature of Switzerland, interesting analyses can be accomplished to see how the media coverage of the pandemic has been managed and what topics in respect to the pandemic have been important to the Swiss media, and hopefully, by proxy to the Swiss public.

In order to help the research community and the public be able to analyse and seek answers to the above questions, we at ELCA decided to release our COVID-19 news dataset, containing more than 10,000 links to news articles related to the Coronavirus pandemic published in the Swiss media since early January 2020.

We hope this dataset helps researchers make insightful analyses on the reaction of the Swiss public to the pandemic and deliver results that help shape a better response in the prospective future cases.

\section{The Data}
We tried to cover the most popular Swiss newspapers and news websites. Therefore, we chose a total of 10 news sources in German, five in French, three in Italian, and also two English-speaking Swiss news websites in order to make the data more accessible to to researchers outside Switzerland. Table \ref{tab11} depicts the list of venues and some information about them.

\begin{table}[htbp]
\caption{Selected News Sources}
\begin{center}
\begin{tabular}{|c|c|c|}
\hline
\textbf{News Source} & \textbf{Language}& \textbf{Number of Articles} \\
\hline
Blick     & German &                2080 \\ \hline
Neue Zürcher Zeitung & German &       1961 \\ \hline
Tages Anzeiger  & German &            1109 \\ \hline
Aargauer Zeitung & German &           1099 \\ \hline
Basler Zeitung & German &             859  \\ \hline
SRF & German &           718   \\ \hline
20 Minuten & German &            553 \\\hline
Berner Zeitung  & German &            281 \\\hline
20 Minutes & French &            225 \\\hline
Tribune de Genève & French &          223 \\\hline
Corriere del Ticino & Italian &        199 \\\hline
24 heures  & French &                 189 \\\hline
World Radio Switzerland & English &     176 \\\hline
Il portale del Ticino  & Italian &     156 \\\hline
RTS & French &             188 \\\hline
SWI swissinfo.ch & English &      112 \\\hline
Le Temps & French &     89 \\\hline
SWI swissinfo.ch  & French &  71 \\\hline
SWI swissinfo.ch  & German &       55 \\\hline
SWI swissinfo.ch  & Italian &       52 \\\hline
\end{tabular}
\label{tab11}
\end{center}
\end{table}

\subsection{Selection of Relevant Articles}
We backdated our data collection to late 2019, and started scanning front pages of the selected news sources in consecutive days, extracting headlines of the articles. Initially, articles with any of the following \emph{keywords} in the title were deemed \emph{"Coronavirus-related"}:
\begin{itemize}
\item Coronavirus,
\item Covid,
\item Lockdown,
\item Pandem* (To account for different spellings of the concept in different languages).
\end{itemize}
This inevitably leads to false negatives. In order to reduce such false negatives, we read at a later stage the synopsis of the article and searched for the keywords also in the body, yielding more positive results. The distribution of the languages in the dataset is depicted in Figure \ref{fig:my_label}.
\begin{figure}
    \centering
    \begin{tikzpicture}
    [
        pie chart,
        slice type={comet}{blu},
        slice type={legno}{rosso},
        slice type={coltello}{giallo},
        slice type={sedia}{viola},
        slice type={caffe}{verde},
        pie values/.style={font={\small}},
        scale=2
    ]

        \pie{}{85.19/comet,8/legno,3.9/sedia,2.8/coltello}

        \legend[shift={(0cm,-1.2cm)}]{{German}/comet, {French}/legno,  {Italian}/coltello,  {English}/sedia}

\end{tikzpicture}
    \caption{Language Distribution of Articles in the Dataset}
    \label{fig:my_label}
\end{figure}

The first article we could find in the Swiss media has been published on January 8\textsuperscript{th} in the French-speaking news portal \emph{20 Minutes}, titled \emph{"A new Coronavirus appears in China"}~\cite{FirstArticle}. We have made a web application to simplify exploring and browsing the data, and reading the collected news articles. The web application is available at \url{https://covidnewsdataset.herokuapp.com/}.
\section*{Summary}

We explained in this article our Swiss COVID-19 dataset and how it has been collected. We publish the dataset hereby for public use, along with an online visualisation application to help explore and look at the news articles of Swiss media during the pandemic in Switzerland. We hope this dataset proves useful in analysis of the pandemic era and the public response to it in Switzerland.

\def\do\/\do-{\do\/\do-}
\Urlmuskip=0mu plus 1mu\relax
\bibliographystyle{./bibliography/IEEEtran}


\section{Introduction}

Recently XENON 1T Collaboration reported that they found electron recoil excess around 2-3 keV
with 3.5 $\sigma$ significance analyzing data for an exposer of 0.65 ton-year \cite{Aprile:2020tmw}.
There is an issue on the tritium contamination to be resolved.
This energy region is sensitive to solar axion search, but the interpretation of this excess in terms
of solar axion is in conflict with astrophysical bounds on the axion coupling to electron.
XENON1T Collaboration also interprets the excess in the context of magnetic moment of solar
neutrino and absorption of light bosonic dark matter \cite{Aprile:2020tmw}.
After the accouncement of XENON1T Collaboration, there appeared a number of papers that
address various issues related with this excess
\cite{Takahashi:2020bpq,OHare:2020wum,Kannike:2020agf,Alonso-Alvarez:2020cdv,Fornal:2020npv,Amaral:2020tga,Boehm:2020ltd,Harigaya:2020ckz,Bally:2020yid,Su:2020zny,Du:2020ybt,DiLuzio:2020jjp,Dey:2020sai,Chen:2020gcl,Bell:2020bes,AristizabalSierra:2020edu,Buch:2020mrg,Choi:2020udy,Paz:2020pbc,Lee:2020wmh,Cao:2020bwd,Robinson:2020gfu,Khan:2020vaf,Primulando:2020rdk,Nakayama:2020ikz,Gelmini:2020xir,Jho:2020sku,Bramante:2020zos,Baryakhtar:2020rwy,An:2020bxd,Zu:2020idx,Gao:2020wer,Budnik:2020nwz,Lindner:2020kko,Bloch:2020uzh,Chala:2020pbn,DeRocco:2020xdt,Dent:2020jhf,McKeen:2020vpf,
Coloma:2020voz,An:2020tcg,DelleRose:2020pbh,Dessert:2020vxy,Bhattacherjee:2020qmv,Ge:2020jfn,Chao:2020yro,Gao:2020wfr,Ko:2020gdg,Hryczuk:2020jhi}.

In this paper, we interpret this electron recoil excess in terms of exothermic DM scattering on atomic
electron bound to Xe in the inelastic DM models.   We shall consider both complex scalar
\cite{Baek:2014kna} and Dirac fermion DM models \cite{Ko:2019wxq} with local $U(1)$ dark gauge
symmetry which is spontaneously broken into its $Z_2$ subgroup by Krauss-Wilczek mechanism
\cite{Krauss:1988zc}.  In this framework, the mass difference ($\delta$)
between the DM and the excited DM (XDM) is generated by dark Higgs mechanism,
and there is no explicit violation of local gauge symmetry related with the presence of dark photon.
On the other hand,  in a number of literature, the mass difference $\delta$ is often introduced
by hand in terms of dim-2 (3) operators for scalar (fermion) DM.
Then local gauge symmetry is broken explicitly and softly.
Introducing dark gauge boson (or dark photon) in such a case would be theoretically inconsistent,
since the current dark gauge fields couple is not a conserved current.
There will appear some channels where high energy behavior of the scattering amplitudes violate
perturbative unitarity, in a similar way with the $W_L W_L \rightarrow W_L W_L$ scattering violates
unitarity if $W$ boson mass is put in by hand.  One of the present authors  pointed out this issue
in fermionic DM model (see Appendix A of Ref.~\cite{Ko:2019wxq}).

Local $Z_2$ scalar \cite{Baek:2014kna}  and fermion DM models \cite{Ko:2019wxq}
have been studied by authors for the $O(100)$ GeV -- O(1) TeV WIMP scenarios.
In this paper we explore the same models for light DM mass  $\lesssim O(1)$ GeV in order to
evade the strong bounds from direct detections experiments searching for signals of DM scattering
on various nuclei. In particular we will emphasize that the DM thermal relic density and the XENON1T
electron recoil excess with a few keV could be simultaneously accommodated if dark Higgs boson
is light enough that ${\rm DM} + {\rm XDM} \rightarrow Z^{' *} \rightarrow Z' \phi$ is kinematically
open.  This channel will play an important role when $m_{DM} < m_{Z'}$, as we  shall demonstrate
in the following.  In order to explain the XENON1T excess in terms of
${\rm XDM} + e_{atomic} \rightarrow {\rm DM} + e_{free}$ with a kinetic mixing, both dark photon and
(X)DM mass should be sub-GeV, more specifically $\sim O(100)$ MeV, in order to avoid the stringent
bounds on the kinetic mixing parameter. For such a light DM, one has to consider the DM annihilation
should be mainly in $p$-wave, and not in $s$-wave,  in order to avoid strong constraints from CMB
(see \cite{Slatyer:2015jla,Leane:2018kjk} and references therein).

For this purpose it is crucial to have dark Higgs $(\phi)$, since they can play a key roles in the $p$-wave
annihilations of DM at freeze-out epoch:
\begin{eqnarray*}
X X^\dagger & \rightarrow & Z^{' *} \rightarrow Z^{'} \phi  , \\
\chi \overline{\chi} & \rightarrow & \phi \phi ,
\end{eqnarray*}
where $X$ and $\chi$ are complex scalar and Dirac fermion DM, respectively.
At freeze-out epoch, the mass gap is too small ($\Delta m \ll T$) and we can consider DM as
complex scalar or Dirac fermion.  In the present Universe,  we have $T \ll \Delta m$ and so we have
to work in the two component DM picture for XENON1T electron recoil.
It can not be emphasized enough that these channels would not be possible without dark Higgs
$\phi$, and  it would be difficult to make the DM pair annihilation be dominated by the $p$-wave
annihilation.

\section{Models for (excited) DM}
\subsection{Scalar DM model}
The dark sector  has a gauged $U(1)_X$ symmetry. There are two scalar particles in the dark sector $X$ and $\phi$ with
$U(1)_X$ charges 1 and 2, respectively. They are neutral under the SM gauge group.
After $\phi$ gets VEV, $\langle \phi \rangle = v_\phi/\sqrt{2}$, the gauge symmetry is spontaneously broken down to discrete $Z_2$.
The  $Z_2$-odd $X$  becomes the DM candidate.
The model Lagrangian is in the form~\cite{Baek:2014kna}
\begin{eqnarray}
{\cal L} & = & {\cal L}_{\rm SM}
- \frac{1}{4} \hat{X}_{\mu\nu} \hat{X}^{\mu\nu} - \frac{1}{2} \sin \epsilon \hat{X}_{\mu\nu} \hat{B}^{\mu\nu}
+ D^\mu \phi^\dagger D_\mu \phi + D^\mu X^\dagger D_\mu X
 - m_X^2 X^\dagger X + m_\phi^2 \phi^\dag \phi
 \nonumber \\
 &&
- \lambda_\phi \left( \phi^\dagger \phi \right)^2 - \lambda_X \left( X^\dagger X \right)^2\
- \lambda_{\phi X} X^\dagger X \phi^\dagger \phi
- \lambda_{\phi H} \phi^\dagger \phi H^\dagger H
- \lambda_{HX} X^\dagger X H^\dagger H
\nonumber\\
&&-  \mu \left( X^2 \phi^\dagger +  H.c. \right),
\label{eq:model}
\end{eqnarray}
where $\hat{X}_{\mu\nu}$ ($B_{\mu\nu}$) is the field strength tensors of $U(1)_X$ ($U(1)_Y$) gauge boson in the interaction basis.

We decompose the $X$ as
\begin{align}
X ={1 \over \sqrt{2}} (X_R + i X_I),
\end{align}
and $H$ and $\phi$ as
\begin{align}
H =
\begin{pmatrix}
0 \\
{1 \over \sqrt{2}} (v_H + h_H)
\end{pmatrix}, \quad
\phi = {1 \over \sqrt{2}} (v_\phi + h_\phi),
\label{eq:ScalarDecomp}
\end{align}
in the unitary gauge.

The dark photon mass is given by
\begin{align}
m_{Z'}^2 \simeq (2 g_X v_\phi)^2,
\end{align}
where we neglected the corrections from the kinetic mixing, which is second order in $\epsilon$ parameter.
The masses of $X_R$ and $X_I$ are obtained to be
\begin{align}
m_R^2 &= m_X^2+{1 \over 2} \lambda_{HX} v_H^2 + {1 \over 2} \lambda_{\phi X} v_\phi^2 + {\mu \over \sqrt{2}} v_\phi, \nonumber\\
m_I^2 &= m_X^2+{1 \over 2} \lambda_{HX} v_H^2 + {1 \over 2} \lambda_{\phi X} v_\phi^2 - {\mu \over \sqrt{2}} v_\phi,
\end{align}
and the mass difference, $\delta \equiv m_R - m_I \simeq \mu v_\phi/\sqrt{2} m_X$.
Since the original $U(1)_X$ symmetry is restored by taking $\mu=0$, small $\mu$ does not give rise to fine-tuning problem.
The mass spectrum of the scalar Higgs sector can be calculated by diagonalising the mass-squared matrix
\begin{align}
\begin{pmatrix}
2 \lambda_H v_H^2 & \lambda_{\phi H} v_H v_\phi \\
\lambda_{\phi H} v_H v_\phi &2 \lambda_\phi v_\phi^2
\end{pmatrix},
\label{eq:H-Mass}
\end{align}
which is obtained in the $(h_H, h_\phi)$ basis.
We denote the mixing angle to be $\alpha_H$ and the mass eigenstates to be $(H_1,H_2)$,
where $H_1$ is the SM Higgs-like state and $H_2 (\equiv \phi)$ is mostly dark Higgs boson.
Since we work in the small $\alpha_H$ in this paper, the VEV of $\phi$ is approximated to be,
$v_\phi \simeq {m_{H_2} / \sqrt{2 \lambda_\phi}}$, while $\alpha_H \simeq  \lambda_{\phi H} v_\phi / 2 \lambda_H v_H$.

The mass eigenstates $Z_\mu$ and $Z'_\mu$ of the neutral gauge bosons can be obtained using the procedure shown in Ref.~\cite{Babu:1997st}.
In the linear order approximation in $\epsilon$ we can write the covariant derivative as
\begin{align}
D_\mu \simeq \partial_\mu + i e Q_{\rm em} A_\mu + i \Big( g_Z (T^3 -Q_{\rm em} s_W^2) + \epsilon g_X Q_X s_W \Big) Z_\mu
+ i \Big( g_X Q_X - \epsilon e Q_{\rm em} c_W \Big) Z'_\mu,
\end{align}
where $Q_{\rm em}$ ($Q_X$) is the electric ($U(1)_X$) charge and $A_\mu$ is the photon field.
We note that $Z'$ couples  to the electric charge but not to the weak isospin component $T^3$.
For example, $Z'$ does not couple to neutrinos at this order of $\epsilon$.

To evade the bound from the DM scattering off the nuclei we are considering sub GeV scale DM.
To calculate the relic abundance of the DM we take the $X$ as the physical state  with mass $m_X$ instead of $X_I$ and $X_R$,
{\it i.e.} $m_X \simeq m_R \simeq m_I$,
because the mass difference $\delta$ is much smaller than the freezeout temperature
$T_f \sim m_X/10$~\cite{Harigaya:2020ckz}.
For this light DM the CMB constraint rules out the $s$-wave annihilation of the DM pair.
So the contribution to the DM annihilation should start from $p$-wave.
We suppress the $X X^\dagger \to Z' Z'$ by choosing $m_{Z'} > m_X$.
We also make $X X^\dagger \to H_2 H_2$ subdominant.
To achieve this we suppress the direct coupling of the DM to $\phi$ and $H$ by taking
small $\lambda_{\phi X}$ and $\lambda_{HX}$.
We also take small $\lambda_{\phi H}$ to evade the bound from the Higgs invisible decay.
However, the coupling $\lambda_{\phi H}$ should not be too small to make the DM in thermal contact with the SM plasma.
Too small $\lambda_{\phi H}$ also makes the $H_2$ lifetime too long, causing cosmological problems.
For example, $H_2$ with mass $\sim 1$ GeV decays dominantly into muon (or strange quark) pairs through mixing $\lambda_{\phi H}$,
whose decay with is given by

\begin{align}
\Gamma(H_2 \to \mu^+ \mu^-) &\simeq \frac{\alpha_H^2}{8 \pi} m_{H_2} \left(m_\mu \over v_H \right)^2 \left(1 -{4 m_\mu^2 \over m_{H_2}^2} \right)^{3/2} \nonumber\\
&\approx \left( {1.1 \times 10^{16}  \, {\rm sec}^{-1}} \right) \alpha_H^2 \left(m_{H_2} \over 1\, {\rm GeV}\right).
\label{eq:Gmm}
\end{align}
We require that $H_2$ lives shorter than 1 sec to evade the constraints from Big Bang nucleosynthesis  (BBN)
 and the mixing angle is small. For $m_{H_2}=1$ GeV it is translated into $\alpha_H > 9.5 \times 10^{-9}$.
When the muon channel is kinematically forbidden but $m_{H_2} > 2 m_e$, it decays into electron-positron pair. The decay rate
is obtained by replacing $m_\mu$ by $m_e$ in (\ref{eq:Gmm}). For example, for $m_{H_2} = 10$ MeV, we
need $\alpha_H > 2.0 \times 10^{-5}$.
When $m_{H_2} < 2 m_e$, the scalar particle decays into two photons.
In this paper we consider $m_{H_2} \gtrsim 2 m_e$.
The small mixing parameters are also technically natural by extending the Poincar\'e symmetry~\cite{Foot:2013hna,Baek:2019wdn}.

The $Z'$ can also decay into charged SM particles through mixing $\epsilon$.
The decay width for $Z' \to e^+ e^-$ is given by
\begin{align}
\Gamma(Z' \to e^+ e^-) &=\frac{Q_f^2 \epsilon^2 e^2 c_W^2}{12 \pi} m_{Z'} \left(1+ {2 m_e^2 \over m_{Z'}^2}\right) \left(1- {4 m_e^2 \over m_{Z'}^2}\right)^{1/2} \nonumber\\
& \approx 1.87 \times 10^{-11} \, {\rm GeV} \left(\epsilon \over 10^{-4} \right)^2 \left(m_{Z'} \over 1 \, {\rm GeV} \right)^2.
\end{align}
Its lifetime is much shorter than 1 sec in the parameter space of our interest.

The light $\phi$ and/or $Z'$ may contribute to the effective neutrino number $N_{\rm eff}$, which is another possible constraint in the model.
Since the mediators $Z'$ and $\phi$ decay before 1 sec, there is no relativistic extra degree freedom which mimics the neutrino at
the recombination era ($T_{\rm CMB} \approx$ 4 eV).
Another potential source for $\Delta N_{\rm eff}$ is light mediator with mass below 1 MeV.
It mainly decays into $e^\pm$ or $\gamma$ not into $\nu$,  making the difference between the temperatures $T_\gamma$ and $T_\nu$ larger
than the one given by the standard cosmology by imparting its entropy only to
$\gamma$~\cite{Matsumoto:2018acr}.
This also causes $\Delta N_{\rm eff} \not= 0$.
We evade this problem by taking their masses greater than 1 MeV.

In this restricted region of parameter space the main channel for the relic density is
$X X^\dagger \to Z^{\prime *}\to Z' \phi$ (see the left panel in Fig. \ref{SDM} )
\footnote{From now on we call $H_2$ as $\phi$.}.
The leading contribution to cross section is $p$-wave with the cross section
\begin{align}
\sigma v & \simeq \frac{g_X^4 v^2}{384  \pi \,m_X^4 (4 m_X^2-m_{Z'}^2 )^2}\left(16 m_X^4 + m_{Z'}^4 + m_\phi^4+40 m_X^2 m_{Z'}^2 -8 m_X^2 m_\phi^2 -2 m_{Z'}^2 m_\phi^2\right) \nonumber\\
&\times \Big[\left\{4 m_X^2-(m_{Z'}+m_\phi)^2\right\}\left\{4 m_X^2-(m_{Z'}-m_\phi)^2\right\}\Big]^{1/2}+{\cal O}(v^4),
\end{align}
where $m_\phi  \equiv \sqrt{2 \lambda_\phi} v_\phi \simeq m_{H_2}$.
The resulting dark matter density is obtained by~\cite{Griest:1990kh}
\begin{align}
\Omega_X h^2 = \frac{2 \times 8.77 \times 10^{-11} \, {\rm GeV}^{-2} \,x_f}{g_*^{1/2} (a+3 b/x_f)},
\label{eq:Omega}
\end{align}
where $a$ and $b$ are defined by $\sigma v = a + b v^2$, and the additional factor 2
comes from the fact that $X$ is a complex scalar instead of real scalar.
For example, with $m_X =1$ GeV, $m_{Z'} =1.2$ GeV, $m_\phi=0.2$ GeV, $\lambda_\phi = 4.5 \times 10^{-5}$,
 $g_* \approx 10$, and $x_f \approx 10$, we can explain the current DM relic abundance: $\Omega_X h^2 \approx 0.12$.
Other $p$-wave contributions include $X X^\dagger \to Z^{\prime *} \to f \bar{f}$ where $f$ is an SM fermion.
But these contributions are suppressed by $\epsilon^2$ compared to the above annihilation, and we neglect them.
The SM $Z$ boson contributions are further suppressed by both small mixing angle $\epsilon^2$ and small mass ratio $m_{Z'}^4/m_{Z}^4$.

To calculate the inelastic down-scattering cross section for the XENON1T anomaly, instead of $X$ and $X^\dagger$ we now consider
two real scalars $X_R, X_I$ with mass difference $\delta$. With the kinetic mixing term given in (\ref{eq:model}) we get the dark-gauge interactions
with the DM and the electron~\cite{Babu:1997st}
\begin{align}
{\cal L} \supset g_X Z^{\prime\mu} (X_R \partial_\mu X_I - X_I \partial_\mu X_R) - \epsilon\, e c_W Z^\prime_\mu \overline{e} \gamma^\mu e,
\end{align}
where $c_W$ is the cosine of the Weinberg angle,  $Z$ and $Z'$ are mass eigenstates, and we assumed that $\epsilon(\sim 10^{-4})$ is small.
The cross section for the inelastic scattering $X_R e \to X_I e$ for $m_X \gg m_e$ and small momentum transfer is given by
\begin{align}
\sigma_e = \frac{16 \pi \epsilon^2 \alpha_{\rm em} \alpha_X c_W^2  m_e^2}{m_{Z'}^4},
\label{eq:sigma_e}
\end{align}
where $\alpha_{\rm em}\simeq 1/137$ is the fine structure constant and $\alpha_X \equiv g_X^2/4 \pi$.
This can be used to predict the differential cross section of the dark matter scattering off the xenon atom for the DM velocity $v$, which reads
\begin{align}
\frac{d \sigma v}{d E_R}=\frac{\sigma_e}{2 m_e v} \int_{q_-}^{q_+} a_0^2 q dq K(E_R,q),
\end{align}
where $E_R$ is the recoil energy, $q$ is the momentum transfer, $K(E_R,q)$ is the atomic excitation factor.
From energy conservation we obtain the relation~\cite{Harigaya:2020ckz},
\begin{align}
E_R = \delta + v q \cos\theta -\frac{q^2}{2 m_R},
\end{align}
where $\theta$ is the angle between the incoming $X_R$ and the momentum transfer $\boldsymbol{q}=\boldsymbol{p'}_e -\boldsymbol{p}_e$.
The integration limits are~\cite{Harigaya:2020ckz},
\begin{align}
q_\pm &\simeq m_R v \pm \sqrt{m_R^2 v^2 - 2 m_R (E_R -\delta)}, \quad \text{for} \; E_R \ge \delta, \nonumber\\
q_\pm &\simeq \pm m_R v +\sqrt{m_R^2 v^2 - 2 m_R (E_R -\delta)}, \quad \text{for} \; E_R \le \delta.
\end{align}
Then we can obtain the differential event rate for the inelastic scattering of DM with electrons in the xenon atoms given by
\begin{align}
\frac{d R}{d E_R} = n_T {n_R} \frac{d \sigma v}{d E_R},
\end{align}
where $n_T \approx 4 \times 10^{27}/$ton is the number density of xenon atoms and $n_R \approx 0.15 \, {\rm GeV}/m_R/{\rm cm}^3$ is
the number density of the heavier DM component $X_R$, assuming $n_R =n_I$. Integrating over $E_R$, we get the event rate
\begin{align}
R \approx 3.69 \times 10^9 \, \epsilon^2 \, g_X^2 \left( 1 {\rm GeV} \over m_R \right) \left( 1 {\rm GeV} \over m_{Z'} \right)^4 /\text{ton}/\text{year}.
\end{align}

Since $X_R$ is a dark matter component in our model with the same abundance with $X_I$, its lifetime should be much longer than the age of the universe.
It can decay via $X_R \to X_I \gamma \gamma \gamma$ as shown in~\cite{Harigaya:2020ckz}.
Its decay into three-body final state, $X_R \to X_I \nu \overline{\nu}$, is also possible in our model. The relevant interactions are
\begin{align}
{\cal L} \supset \epsilon g_X s_W Z^\mu (X_R \partial_\mu X_I -X_I \partial_\mu X_R) -{g_Z \over 2} Z_\mu \overline{\nu}_L \gamma^\mu \nu_L.
\end{align}
The decay width is given by
\begin{align}
\Gamma \simeq \frac{\epsilon^2 \alpha_X s_W^2}{5\sqrt{2} \pi^2} \; \frac{G_F \delta^5 }{m_Z^2} \simeq 1.9 \times 10^{-49} \,\text{GeV}
\left(\epsilon \over 10^{-4}\right)^2 \left(\alpha_X \over 0.078\right)\left(\delta \over 2\,\text{keV}\right)^5.
\label{eq:GammaX}
\end{align}
Although this channel is much more effective than $X_R \to X_I \gamma\gamma\gamma$ considered in~\cite{Harigaya:2020ckz},
the lifetime of $X_R$ is still much longer than the age of the universe.

In the right panel of Fig. \ref{SDM} , we show the allowed region in the $(m_{Z'}, \epsilon)$ plane
where we can explain the XENON1T excess with correct thermal relic density of DM within the
standard freeze-out scenario.  For illustration, we chose the DM mass to be $m_R = 0.1$ GeV,
and varied the dark Higgs mass $m_\phi = 20, 40, 60, 80$ MeV denoted with different colors.
The sharp drops on the right allowed region is from the kinematic boundary,
$m_{Z'} + m_\phi < 2 m_R$. It is nontrivial that we could explain the XENON1T excess with
inelastic DM models with spontaneously broken  $U(1)_X \rightarrow Z_2$ gauge symmetry.
In particular it is important to include light dark Higgs for this explanation.
It would be straightfoward to scan over all the parameters to get the whole allowed region.
\begin{figure}
\includegraphics[width=0.40\linewidth]{Relic_ScalarDM.pdf}
\includegraphics[width=0.50\linewidth]{mZp_VS_eps_ScalarDM100MeV.pdf}
\caption{ \label{SDM}
({\it left}) Feynman diagrams relevant for thermal relic density of DM:
 $X X^\dagger \rightarrow Z' \phi$ and ({\it right}) the region in the $( m_{Z'}, \epsilon)$ plane
that is allowed for the XENON1T electron recoil excess and the correct thermal relic density
for scalar DM case for $\delta = 2$ keV : (a) $m_{\rm DM} = 0.1$ GeV. Different colors represents
$m_\phi = 20,40,60,80$ MeV.    The gray areas are excluded by various experiments, from BaBar \cite{Lees:2014xha}, E774 \cite{Bross:1989mp}, E141 \cite{Riordan:1987aw}, Orasay \cite{Davier:1989wz}, and E137 \cite{Batell:2014mga}, assuming $Z' \rightarrow X_R X_I$ is kinematically forbidden.
}
\end{figure}

\subsection{Fermion DM model }

We start from a dark $U(1)$ model, with a Dirac fermion dark matter (DM) $\chi$ appointed with a nonzero dark $U(1)$ charge $Q_\chi$ and dark photon.   We also introduce a complex dark Higgs field $\phi$, which takes a nonzero vacuum expectation value, generating nonzero mass for the dark photon. We shall consider a special case where $\phi$ breaks the dark $U(1)$ symmetry into a dark $Z_2$ symmetry with a judicious choice of its dark charge $Q_\phi$.

Then the gauge invariant and renormalizable Lagrangian for this system is given by
\begin{eqnarray}
\mathcal{L} & = & -\frac{1}{4} \hat{X}^{\mu \nu} \hat{X}_{\mu \nu}
- \frac{1}{2} \sin \epsilon \hat{X}_{\mu\nu}B^{\mu\nu}
+ \overline{\chi} \left( i\slashed{D} - m_{\chi} \right) \chi + D_{\mu} \phi^\dagger D^{\mu} \phi
\label{Lag_fermion}
\\
& - & \mu^2 \phi^\dagger  \phi - \lambda_{\phi} |\phi|^4
-  \frac{1}{\sqrt{2}} \left( y \phi^\dagger \overline{\chi^C} \chi + \text{h.c.} \right)
- \lambda_{\phi H} \phi^\dagger \phi H^\dagger H
\nonumber
\end{eqnarray}
where $\hat{X}_{\mu \nu} = \partial_{\mu} \hat{X}_{\nu} - \partial_{\nu} \hat{X}_{\mu}$.
$D_{\mu} = \partial_{\mu} + i g_X Q_X \hat{X}_{\mu}$ is the covariant derivative,
where $g_X$ is the dark coupling constant, and $Q_X$ denotes the dark charge of $\phi$ and $\chi$:
$Q_{\phi}  = 2,  Q_{\chi} = 1$, respectively.
Then $U(1)_X$ dark gauge symmetry is spontaneously broken into its $Z_2$ subgroup, and
the Dirac DM $\chi$ is split into two Majorana DM $\chi_R$ and $\chi_I$ defined as
\begin{eqnarray}
\chi &  = & \frac{1}{\sqrt{2}} ( \chi_R + i \chi_I ) ,
\\
\chi^c &  = & \frac{1}{\sqrt{2}} ( \chi_R - i \chi_I ) ,
\\
\chi_R^c & = & \chi_R , \ \ \  \chi_I^c = \chi_I ,
\end{eqnarray}
with
\begin{equation}
m_{R,I} = m_\chi \pm  y v_\phi = m_\chi \pm \frac{1}{2} \delta .
\end{equation}

We assume  $y>0$ so that $\delta \equiv m_R - m_I = 2 y v_\phi > 0$.
Then the above Lagrangian is written as
\begin{eqnarray}
{\cal L} & = & \frac{1}{2} \sum_{i=R,I} \overline{\chi_i} \left( i \slashed{\partial} - m_i \right) \chi_i
- i \frac{g_X}{2}  ( Z'_\mu+\epsilon s_W Z_\mu) \left( \overline{\chi_R} \gamma^\mu \chi_I
- \overline{\chi_I} \gamma^\mu \chi_R \right)
\\
& - & {1 \over 2} y h_\phi \left( \overline{\chi_R} \chi_R  - \overline{\chi_I} \chi_I \right),
\end{eqnarray}
where $h_\phi$ is neutral CP-even component of $\phi$ as defined in (\ref{eq:ScalarDecomp}).

When we calculate the DM relic density, we can assume the mass difference is small compared
to the DM mass as in the case of the scalar DM, {\it i.e.} $m_\chi \simeq m_R \simeq m_I$.
For the fermionic DM the annihilation processes into the scalar pair, $\stackrel{(-)}{\chi} \stackrel{(-)}\chi \to \phi \phi$,  are $p$-wave.
To evade the CMB constraint we suppress the $s$-wave annihilation by assuming $2 m_\chi<2m_{Z'}, m_{Z'}+m_\phi$.
The calculation of the annihilation process (the top panel of Fig. \ref{FDM}) yields
\begin{align}
\sigma v = \frac{y^2 v^2 \, \sqrt{m_\chi^2-m_\phi^2}}{96\pi m_\chi}
\Bigg[\frac{27 \lambda_\phi^2 v_\phi^2}{(4 m_\chi^2 -m_\phi^2)^2}
+\frac{4 y^2 m_\chi^2 (9 m_\chi^4 -8 m_\chi^2 m_\phi^2 +2 m_\phi^4)}{(2 m_\chi^2 -m_\phi^2)^4}
 \Bigg] +{\cal O}(v^4),
\end{align}
where $m_\phi  \equiv \sqrt{2 \lambda_\phi} v_\phi \simeq m_{H_2}$.
The current DM relic abundance is obtained by (\ref{eq:Omega}).  Since the annihilation cross
section for the fermion DM case if proportional to $y^2 \propto (\delta/v_\phi )^2$ and
$\delta \sim 2$ keV, the $v_\phi$  should be not too large. If we ignore $m_\phi$ in the above
equation, we find that $m_\chi \sim O(1-10)$ MeV will be required to get the correct thermal relic
density, and dark Higgs $\phi$ should be even lighter. Therefore dark sector particles in this
case should be lighter than the scalar DM case.
Later for illustration, we will consider $m_\chi \sim O(10)$ MeV. and $m_\phi \sim O(1)$ MeV
to get the correct DM relic density and explain the XENON1T excess.

Now let's consider the inelastic scattering of DM  with the electron in the xenon atom to explain the XENON1T anomaly.
The scattering occurs through the interactions
\begin{align}
{\cal L} \supset
- i \frac{g_X}{2}   Z'_\mu \left( \overline{\chi_R} \gamma^\mu \chi_I
- \overline{\chi_I} \gamma^\mu \chi_R \right)
- \epsilon\, e c_W Z^\prime_\mu \overline{e} \gamma^\mu e,
\end{align}
It turns out that in the limit $m_\chi \gg m_e$, $\sigma_e$ has exactly the same form with (\ref{eq:sigma_e}) of the scalar DM case.

We require the $\chi_R$ to be long-lived so that it is also a main component of the dark matter.
It decays mainly via the SM $Z$-mediating $\chi_R \to \chi_I \nu \overline{\nu}$, using the interactions
\begin{align}
{\cal L} \supset -{i \over 2} \epsilon s_W g_X Z_\mu (\overline{\chi}_R \gamma^\mu \chi_I -\overline{\chi}_I \gamma^\mu \chi_R)
-{1 \over 2} g_Z Z_\mu \overline{\nu}_L \gamma^\mu \nu_L.
\end{align}
The expression for the decay with, $\Gamma(\chi_R \to \chi_I \nu \overline{\nu})$, also agrees exactly with (\ref{eq:GammaX}).
As shown in (\ref{eq:GammaX}), the lifetime of $\chi_R$ is much longer than the age of the universe, which guarantees the $\chi_R$ is as good a dark matter as $\chi_I$.

\begin{figure}
\includegraphics[width=0.40\linewidth]{Relic_FermionDM1.pdf}
\hspace{0.5cm}
\includegraphics[width=0.40\linewidth]{Relic_FermionDM2.pdf}
\\
\vspace{0.5cm}
\includegraphics[width=0.50\linewidth]{mZp_VS_eps_FermionDM10MeV.pdf}
\caption{ \label{FDM}
({\it top}) Feyman diagrams for $\chi \bar{\chi} \rightarrow \phi\phi$.
({\it bottom}) the region in the $( m_{Z'}, \epsilon)$ plane
that is allowed for the XENON1T electron recoil excess and the correct thermal relic density
for fermion DM case for $\delta = 2$ keV and the fermion DM mass to be $m_R = 10$ MeV.
Different colors represents $m_\phi = 2, 4, 6, 8$ MeV.    The gray areas are excluded by various
experiments, assuming $Z' \rightarrow \chi_R \chi_I$ is kinematically allowed, and the experimental
constraint is weaker in the $\epsilon$ we are interested in, compared with the scalar DM case
in Fig. 1 (right). We also show the current experimental bounds by NA64 \cite{NA64:2019imj}.
}
\end{figure}

In the bottom panel of Fig. \ref{FDM}, we show the allowed region in the $(m_{Z'}, \epsilon)$ plane
where we can explain the XENON1T excess with correct DM thermal relic density within the
standard freeze-out scenario.  For illustration, we choose the fermion DM mass to be
$m_\chi = m_R = 10$ MeV,   and varied the dark Higgs mass $m_\phi = 2, 4, 6, 8$ MeV denoted
with different colors.    Note that the kinetic mixing $\epsilon \sim 10^{-7 \pm 1}$, which is much
smaller than the scalar DM case.
We have checked if the gauge coupling $g_X$ and the quartic coupling of dark Higgs
($\lambda_\phi$) remain in the perturbative regime. The solid (dashed) lines denote the region where
$g_X$ satisfy (violate) perturbativity condition, depending $\alpha_X < 1$ or not.
Within this allowed region, $\lambda_\phi$ remain perturbative.
Again it is nontrivial that we could explain the XENON1T excess with inelastic fermion DM
models with spontaneously broken  $U(1)_X \rightarrow Z_2$ gauge symmetry.
In particular it is important to include light dark Higgs for this explanation as in the scalar DM case.

\section{Conclusion}

In this paper, we showed that the electron recoil excess reported by XENON1T Collaboration could
be accounted for by exothermic DM scattering on atomic electron in Xe, with sub-GeV light
DM: $m_X \sim O(100)$ MeV for the scalar and $m_\chi \sim O(10)$ MeV for the fermion DM,
and dark Higgs $\phi$ neing even lighter that DM particle for both cases.
Dark photon should be heavier than DM in order that we can forbid the DM pair  annihilation
into the $Z' Z'$ channels.
This scenario could be described by DM models with dark $U(1)$ gauge symmetry broken
into its $Z_2$ subgroup by Krauss-Wilczek mechanism. And dark photon $Z'$ and dark Higgs $\phi$
in such dark gauge models play important roles in DM phenomenology.
In particular in the calculation of thermal relic density,  new channels involving a dark Higgs can
open $X X^\dagger \rightarrow \phi Z'$ and  $\chi \bar{\chi} \rightarrow \phi \phi$,
which are $p$-wave annihilations. Then one could evade  the stringent constraints from CMB
for such light DM.   Other dangerous $s$-channel annihilations can be kinematically forbidden
by suitable choice of parameters.
Thus the exothermic scattering in inelastic $Z_2$ DM models within standard freeze-out scenario
can explain the XENON1T excess without modifying early universe cosmology.
We emphasize again that the existence of dark Higgs $\phi$ is crucial for us to get the desired DM
phenomenology to explain the XENON1T excess with the correct thermal relic density in case of
both  scalar and fermion DM models within the standard freeze-out scenario.

\section*{Note Added}
While we were preparing this manuscript, there appeared a few papers which explain
the XENON1T excess in terms of scalar or fermion exothermic DM \cite{Harigaya:2020ckz,Su:2020zny,Lee:2020wmh,Bramante:2020zos,Baryakhtar:2020rwy,An:2020tcg}.
Our paper is different from these previous works in that we consider dark $U(1)$ gauge symmetry
broken to its $Z_2$ subgroup by dark Higgs mechanism, and include the light dark Higgs in the
calculations of thermal relic density for the two component DM in the standard freeze-out scenario:
$X X^\dagger \rightarrow Z' \phi$ for scalar DM and $\chi \overline{\chi} \rightarrow \phi\phi$.
Other $s$-wave channels are kinematically forbidden by suitable choice of mass parameters.
This possibility of dark Higgs in the final state in the (co)annihilation channels are not included
in other works.  By including these new channels, we could achieve the correct thermal relic
density and the desired heavier DM fluxes on the Xe targets simultaneously, without conflict with
strong constraints on light DM annihilations from CMB.
And the models considered in this paper is renormalizable and DM stability is guaranteed by
underlying $U(1)$ dark gauge symmetry and its unbroken $Z_2$  subgroup.

\begin{acknowledgments}
The work is supported in part by KIAS Individual Grants, Grant No. PG021403 (PK) and
Grant No. PG074201 (JK) at Korea Institute for Advanced Study, and
by National Research Foundation of Korea (NRF) Grant No. NRF-2018R1A2A3075605 (SB)
and No. NRF-2019R1A2C3005009 (PK), funded by the Korea government (MSIT).
\end{acknowledgments}

\bibliographystyle{utphys}