\section{Introduction}
\label{sec:intro}

Eigenvalues of graphs have diverse applications in combinatorics. We refer to \cite{Bi93, CDS95, GR, BH, Sp12} and references therein.
On the other hand, results about eigenvalues of digraphs are sparse. One reason is that it is not clear which matrix associated to a digraph $D$ would best reflect interesting combinatorial properties in its spectrum. One candidate is the \emph{adjacency matrix} $A=A(D)$ whose $(u,v)$-entry is 1 if there is an arc from the vertex $u$ to $v$, and 0 otherwise. A well-known theorem of Wilf \cite{Wilf} bounding the chromatic number in terms of its largest eigenvalue extends to this setting as shown in \cite{M10}. However, the disadvantage of this matrix is that it is not symmetric and we lose the property that eigenvalues are real. Moreover, algebraic and geometric multiplicities of eigenvalues may be different. Another candidate is the \emph{skew-symmetric adjacency matrix} $S(D)$, where the $(u,v)$-entry is 1 if there is an arc from $u$ to $v$, and $-1$ if there is an arc from $v$ to $u$ (and $0$ otherwise). This choice is quite natural but it only works for \emph{oriented graphs} (i.e.\ when we have no digons). We refer to a survey by Cavers et al.\ \cite{CaCietal12}.

A Hermitian adjacency matrix for digraphs and mixed graphs was introduced recently by Guo and Mohar \cite{GuMo}, who provided many basic properties of this matrix. The same matrix was independently used by Liu and Li \cite{LiuLi15}, who considered an application in mathematical chemistry. The $uv$-entry of this Hermitian matrix is equal to 1 if there is an unoriented edge between $u$ and $v$, it is equal to $i=\sqrt{-1}$ if $uv$ is an arc and is equal to $-i$ if $vu$ is an arc. Of course, one may ask whether taking complex numbers $i$ and $-i$ is the most natural choice. In this note we offer another choice and argue why that one may be more natural to be used when investigating relationship between eigenvalues and combinatorial properties of digraphs and, more generally, of mixed multigraphs.

With this new Hermitian matrix, each arc directed from $v$ to $u$ contributes the sixth root of unity $\omega = \bigl(1+i\sqrt{3}\,\bigr)/2$ to the $vu$-entry in the matrix and contributes $\overline{\omega} = \bigl(1-i\sqrt{3}\,\bigr)/2$ to the $uv$-entry. Note that two oppositely directed arcs between $u$ and $v$ together contribute 1 to each of the entries. In this way, digons have the same effect on the adjacency matrix as undirected edges.

The main reason why the sixth root of unity is natural in relation to combinatorial properties is that $\omega \cdot \overline{\omega} = 1$ and $\omega + \overline{\omega} = 1$. The product of ``weights'' of edges is natural with the counting of closed walks (which is related to the entries of the powers of the adjacency matrix). The sum is needed when dealing with multiple arcs or with weighted edges.

The sixth root of unity appears naturally across applications. It appears in the definition Eisenstein integers; in relation to abstract linear independence (matroids), the sixth root matroids play a special role next to regular and binary matroids \cite{OVW98,Wh97}. They also seem to have connections to Quantum Field Theory \cite{KaKn10}.

In this paper we introduce this new Hermitian matrix and prove a few basic results about its eigenvalues.

\section{Basic properties}
\label{sec:defn}

\subsection{Digraphs and their Hermitian matrices}

A \emph{directed graph} (or \emph{digraph}) $D$ consists of a finite set $V=V(D)$ of vertices together with a set $E=E(D)$ of \emph{arcs} or \emph{directed edges}. Each arc $e\in E$ joins two vertices $u,v\in V$ and one of these, say $v$, is its \emph{initial vertex}, while the other one, $u$, is its \emph{terminal vertex}. To denote this incidence we write $e=vu$ for short. If $vu \in E$ and $uv \in E$, we say that the pair $\{uv, vu\}$ of these oppositely directed arcs is a \emph{digon} of $D$.

Multiple arcs and multiple edges between the same pair of vertices are allowed. Although we sometimes write $uv\in E(D)$, we mean by this any arc from $u$ to $v$. Although we will not discuss loops, they may be present, but they should be considered as undirected (or as ``two oppositely oriented" loops at the same vertex).

Let $\a = a+bi$ be a complex number with absolute value 1, $|\a|=1$, where $a\ge0$, and let $\bar\a = a-bi$ be its conjugate. For a digraph $D$ with vertex set $V = V(D)$ and arc set $E = E(D)$, we consider the \emph{Hermitian adjacency matrix} $N^\a = N^\a(D) \in \mathbb{C}^{V\times V}$, whose entries $N_{uv}^\a$ are given by
\[
N_{uv}^\a = e(u,v)\,\a + e(v,u)\, \bar\a,
\]
where $e(x,y)$ denotes the number of arcs from $x$ to $y$. Observe that $N_{uv}^\a$ and $N_{vu}^\a$ are conjugate to each other, and therefore $N^\a$ is a Hermitian matrix. If $D$ has loops, each loop contributes $\a+\bar\a = 2a$ to the corresponding diagonal entry.

\subsection{About the most natural choice of $\a$}

Let
$$
   \omega = \frac{1+i\sqrt{3}}{2}
$$
be the primitive sixth root of unity and let $\overline{\omega} = \tfrac{1-i\sqrt{3}}{2}$ be its conjugate.

For a digraph $D$, we consider the corresponding Hermitian adjacency matrix $M(D) = N^\omega(D)$ and we refer to it as the \emph{Hermitian matrix of the second kind}.

If every edge of $D$ lies in a digon, then $M(D) = A(D)$, because $\omega + \overline{\omega}=1$. This reflects that $D$ is, essentially, equivalent to an undirected graph in such a case. More generally, a mixed graph is a graph where directed and undirected edges may coexist. Formally, a \emph{mixed graph} is an ordered triple $(V, E, A)$ where $V$ is the vertex-set, $E$ is a set of undirected edges, and $A$ is set of arcs, or directed edges. The Hermitian adjacency matrix of the second kind is defined in such a way that all undirected edges may be replaced by digons and, from this perspective, mixed graphs are equivalent to the class of digraphs that we consider here.

The number $\omega$ used in defining the matrix $M(D)$ satisfies the following identities:
$$
    \omega \cdot \overline{\omega} = 1 \quad \textrm{and} \quad \omega + \overline{\omega} = 1.
$$
While the first condition holds for any $\a$ of absolute value 1, the second identity gives an important \emph{additivity} property and enables us to view oppositely oriented arcs as an unoriented edge and works naturally with multiple arcs and multiple edges. The first one is multiplicative and plays the role in expressing powers of the matrix $M^k$ which correspond, in the same way as the powers of the usual adjacency matrices, to counting ``walks'' of length $k$, the second one allows this property to extend when multiple edges are present, and also works in the setting where the edges have arbitrary positive weights.

Besides the additivity property being implied by our choice of the sixth root of unity, there are several other reasons why this choice may be the most natural among all possible choices for the entries of Hermitian adjacency matrices. The sixth root of unity appears naturally across applications. One such instance is that the sixth-root matroids play a special role in formal theory of linear independence, see, for example, \cite{OVW98,Wh97,Re09,CDN15}. Another natural setup comes from theoretical physics, see \cite{KaKn10} for an example.

\subsection{Eigenvalues}

Observe that $N^\a(D)$ is a Hermitian matrix and so is diagonalizable with real eigenvalues. The following proposition contains properties that are true for adjacency matrices which also carry over to the Hermitian case.

\begin{proposition}\label{prop:real}
Every Hermitian adjacency matrix $N^\a = N^\a(D)$ of a digraph $D$ with vertex-set $V$ has the following properties:
\begin{enumerate}[\rm (i)]
\item All eigenvalues of $N^\a$ are real.
\item The matrix $N^\a$ has $|V|$ pairwise orthogonal eigenvectors in $\mathbb{C}^V$ and is unitarily similar to a diagonal matrix.
\item The numerical range of $N^\a$, defined as the set $R = \{z^*N^\a z \mid z\in \mathbb{C}^V, \Vert z \Vert = 1\}$ is an interval of real numbers and $\min R = \lambda_{\rm min}(N^\a)$ is the smallest eigenvalue of $N^\a$ and $\max R = \lambda_{\rm max}(N^\a)$ is the largest eigenvalue of $N^\a$.
\end{enumerate}
\end{proposition}

The eigenvalues of $N^\a(D)$ are the \emph{$N^\a$-eigenvalues of $D$} and the spectrum of $N^\a(D)$ (i.e.\ the multiset of eigenvalues, counting their multiplicities) is the \emph{$N^\a$-spectrum of $D$}. Similar terminology is used for the particular case of the Hermitian matrix $M(D)$  of the second kind. The $N^\a$-eigenvalues of a digraph $D$ will be ordered in the decreasing order, the $j$th largest eigenvalue will be denoted by $\nu_j^{(\a)}(D)$, so that $\nu_1^{(\a)}(D)\ge \nu_2^{(\a)}(D)\ge \cdots \ge \nu_n^{(\a)}(D)$ ($n=|V|$). For the special case of the matrix $M(D)$, the eigenvalues are $\mu_1(D)\ge \mu_2(D)\ge \cdots \ge \mu_n(D)$.

A direct consequence of Proposition \ref{prop:real} is the min-max formula for $\nu_j^{(\a)}(D)$:
\begin{equation}
  \nu_j^{(\a)}(D) ~= \max_{\dim U=j} ~ \min_{\substack{z\in U\\ \Vert z \Vert = 1}} ~ z^*N^\a z ~= \min_{\dim U=n-j+1} ~ \max_{\substack{z\in U\\ \Vert z \Vert = 1}} ~ z^*N^\a z
\label{eq:minmax jth largest}
\end{equation}
where the outer maximum (minimum) is taken over all subspaces $U$ of $\mathbb{C}^V$ of dimension $j$ ($n-j+1$, respectively) and the inner minimum (maximum) is taken over all unit vectors $z\in U$.

For $z\in \mathbb{C}^V$, the following expression of the quadratic form $z^*N^\a z$ shows how individual arcs contribute to it. In this expression we use the summation over all arcs $vu$, which we indicate as $\sumvtou{f(v,u)}$. The sum runs over all arcs $e=vu$, where each arc is considered as many times as its multiplicity. Let $z=x+iy$, where $x,y\in \mathbb{R}^V$. With $\a = a+bi$ we have:
\begin{eqnarray}
   z^*N^\a z & = & \sum_{v\in V} \overline{z_v} \sum_{u\in V} N_{vu}^\a z_u \nonumber \\
         & = & \sum_{v\in V} (x_v-iy_v) \sum_{u\in V} N_{vu}^\a (x_u+iy_u) \nonumber \\
         & = & \sumvtou{((x_v-iy_v)(x_u+iy_u)\,\a+(x_u-iy_u)(x_v+iy_v)\,\bar\a)} \nonumber \\
         & = & \sumvtou{(2a\,x_vx_u + 2a\,y_vy_u - 2b\,x_vy_u + 2b\,y_vx_u)}.
\label{eq:zNz}
\end{eqnarray}

This expression has a clear combinatorial meaning when combined together with the min-max formula (\ref{eq:minmax jth largest}).

Another property from linear algebra that is used in many combinatorial applications of graph eigenvalues is interlacing. Since $N^\a(D)$ is Hermitian, it has the interlacing property which we describe next.

Suppose that $\nu_1\ge \nu_2\ge \cdots \ge \nu_n$ and $\kappa_1\ge \kappa_2\ge \cdots \ge \kappa_{n-t}$ (where $t\ge1$ is an integer) are sequences of real numbers. We say that the sequences $\nu_l$ ($1\le l\le n$) and $\kappa_j$ ($1\le j\le n-t$) \emph{interlace} if for every $s=1,\dots,n-t$, we have
$$
    \nu_s\ge \kappa_s\ge \nu_{s+t}.
$$

The usual version of the eigenvalue interlacing property states that the eigenvalues of any principal submatrix of a Hermitian matrix interlace those of the whole matrix (see \cite[Theorems 4.3.8 and 4.3.15]{HoJo13}).
This implies that the eigenvalues of any induced subdigraph interlace those of the digraph itself.

\begin{corollary}
\label{cor:interlacing digraph}
The $N^\a$-eigenvalues of an induced subdigraph interlace the $N^\a$-eigenvalues of the digraph.
\end{corollary}

To see a simple example how useful the interlacing theorem is, let us consider the following notion.
Let $\eta^+(D)$ denote the number of non-negative $N^\a$-eigenvalues of a digraph $D$ and $\eta^{-}(D)$ denote the number of non-positive $N^\a$-eigenvalues. The spectral bound of Cvetkovi\'{c} (see \cite{CDS95}) for the largest independent set of a graph extends to digraphs.
Here we will say that a vertex-set $S\subseteq V(D)$ is \emph{independent} if no two vertices in $S$ are joined by an arc in $D$.

\begin{proposition}
\label{prop:indep}
If $D$ has an independent set of size $k$, then for every $\alpha$, we have that $\eta^+(D) \geq k$ and $\eta^-(D) \geq k$.
\end{proposition}

\begin{proof} Let $\nu_1 \geq \cdots \geq \nu_n$ be the eigenvalues of $N^\a(D)$. By interlacing, we see that $\nu_k \geq 0$ and so $N^\a(D)$ has at least $k$ non-negative eigenvalues. Applying the same argument to $-N^\a(D)$ shows that there are at least $k$ non-positive eigenvalues as well.
\end{proof}

Since the independent sets in $D$ correspond to independent sets in the underlying undirected graph $G$, one can optimize the bound of Proposition \ref{prop:indep} over all orientations of $G$ and over different choices of $\alpha$. One can even add any real weights on the edges (which gives rise to a ``complex version" of the Lov\'asz $\vartheta$-function).

\subsection{Symmetry of the spectrum}

There is an essential difference when considering the symmetry of the spectra of the two kinds of Hermitian adjacency matrices. For the Hermitian matrix of the first kind from \cite{GuMo}, the spectrum of every digraph without loops or digons is symmetric about 0. This is no longer true for the Hermitian matrix of the second kind (see the example in Figure \ref{fig:digraph large negative}(a)). However, we still have the behaviour from the undirected case. Recall that a digraph $D$ is \emph{bipartite} if there is a \emph{bipartition} $V(D) = X \cup Y$ such that every arc in $D$ has one end in $X$ and the other end in $Y$.

\begin{theorem}
\label{thm:bipartite}
If D is bipartite, then the $M$-spectrum of $D$ is symmetric about 0, i.e., if $\nu$ is an eigenvalue of $M(D)$, then $-\nu$ is an eigenvalue of $M(D)$ with the same multiplicity.
\end{theorem}

\begin{proof}
The theorem is a simple consequence of the following fact: If $Mz=\nu z$, let $z'$ by the vector which agrees with $z$ on $X$ and agrees with $-z$ on $Y$, where $X\cup Y$ is the bipartition of $D$. Then $Mz' = -\nu z'$. The details are left to the reader.
\end{proof}

An alternative proof of Theorem \ref{thm:bipartite} can be given by looking at the coefficients of the characteristic polynomial $\phi(M(D),\nu) = \sum_{k=0}^n c_k \nu^{n-k}$. These coefficients have combinatorial interpretation via expansion of the determinant as the sum over all permutations of $V(D)$. If $D$ is bipartite, there are no cycles of odd length, thus all coefficients $c_k$ with $k$ odd are zero, and hence the eigenvalues are symmetric with respect to 0.

\begin{figure}
\centering
\includegraphics[scale=1.0]{SymmetricSpectrum.pdf}
\caption{A nonbipartite graph whose characteristic polynomial is $\phi(N,\nu)=\nu^5-7\nu^3+6\nu$ and its spectrum is symmetric about 0.}\label{fig:symmetricspectrum}
\end{figure}

There are nonbipartite digraphs whose $M$-spectrum is symmetric about 0. An example is given in Figure \ref{fig:symmetricspectrum}.

\section{Spectral radius}
\label{sec:rho}

Recall that the \emph{spectral radius} $\rho(Q)$ of a matrix $Q$ is defined as
\[ \rho(Q) = \max\{|\lambda| \mid \lambda \text{ an eigenvalue of }Q \}
\]
and we also define $\rho(D) = \rho(M(D))$ as the \emph{spectral radius} of the digraph $D$.

Let us start with a simple observation.

\begin{proposition}
\label{prop:average degree}
Suppose that a mixed graph $G$ with $n$ vertices has $s$ arcs and $t$ undirected edges. Let $d = \tfrac{s+2t}{n}$. Then
$\mu_1(G) \ge d$.
\end{proposition}

\begin{proof}
To prove the inequality, we just take the constant real vector $x$ with coordinates $x_v=n^{-1/2}$. Clearly, $\Vert x\Vert = 1$. By (\ref{eq:zNz}), it is easy to see that $x^*Mx = d$ (since $y=0$). This implies that $\mu_1(G)\ge d$.
\end{proof}

\begin{figure}
\centering
\includegraphics[scale=1.2]{Circulant5.pdf}
\caption{(a) A simple circulant $Z_5$ with $\mu_1(Z_5) \approx 2.165$ and $\mu_5(Z_5) \approx -3.165$. (b) A circulant $C$ with multiple arcs has $\mu_1(C) \approx 4.0418$ and $\mu_5(C) \approx -6.8195$.}\label{fig:digraph large negative}
\end{figure}

For a digraph $D$, let the eigenvalues of $N^\a = N^\a(D)$ be $\nu_1 \geq \cdots \geq \nu_n$. Since $N^\a$ is not a matrix with non-negative real entries, there is no analogue of the Perron value of the adjacency matrix and the properties of $\nu_1$ may be a bit unintuitive. In particular, it may happen that $\nu_1 < |\nu_n| = \rho(N^\a)$. For example, the digraph shown in Figure \ref{fig:digraph large negative}(b) has $|\mu_5|/\mu_1 > 1.68$. The purpose of this section is to show that $|\nu_n|$ is still bounded in terms of $\nu_1$.

We will first give a general upper bound that holds for arbitrary Hermitian matrices $N^\a(D)$ as long as the real part of $\a$ is non-negative.

\begin{theorem}
\label{thm:spectral_radius}
Let $D$ be a digraph and $\a = a+bi$, where $a,b\in \mathbb{R}$ and $a\ge0$. Then
$$
   \tfrac{1}{3}\, \rho(N^\a(D)) \le \nu_1(N^\a(D)) \le \rho(N^\a(D)).
$$
\end{theorem}

\begin{proof}
While the second inequality follows from the definition of the spectral radius, the first one needs a proof. We may assume that $|\a|=1$. If $\rho(N^\a)\ne \nu_1$, then $\rho = \rho(N^\a) = |\nu_n|$. Let $z=x+iy$ be a unit eigenvector for $\nu_n$. We may assume w.l.o.g. that $x\ne 0$ and $y\ne 0$. By (\ref{eq:zNz}), we have
\begin{equation}\label{eq:XYZ}
   -\rho = \nu_n = z^*N^\a z = X + Y + Z,
\end{equation}
where
$$
   X = 2a \sumvtou{x_vx_u}, \quad Y = 2a \sumvtou{y_vy_u}, \quad \textrm{and} \quad Z = 2b \sumvtou{(y_vx_u - x_vy_u)}.
$$

\newcommand\hatx{\hat x}
\newcommand\haty{\hat y}

Suppose first that $|X|+|Y|\ge \rho/3$. Let ${\hatx},{\haty}\in\mathbb{R}^V$ be real vectors whose coordinates are $|x_v|$ and $|y_v|$, respectively. Then ${\hatx}^* N^\a {\hatx} \ge |X|$ and ${\haty}^* N^\a {\haty} \ge |Y|$. Since $\Vert {\hatx}\Vert^2 + \Vert {\haty}\Vert^2 = \Vert z\Vert^2 = 1$, it follows that
${\hatx}^* N^\a {\hatx} / \Vert {\hatx}\Vert^2 \ge |X|+|Y|$ or ${\haty}^* N^\a {\haty} / \Vert {\haty}\Vert^2 \ge |X|+|Y|$. In either case, it follows that $\nu_1\ge |X|+|Y|\ge \rho/3.$

Suppose now that $|X|+|Y|\le \rho/3$. Since $\rho = -X-Y-Z \le |X|+|Y| - Z$, we conclude that $\tfrac{2}{3}\,\rho \le -Z$.
Let us consider $\overline z = x -iy$. It follows that
$$
   \nu_1 \ge {\overline z}^* N^\a {\overline z} = X + Y - Z.
$$
By adding (\ref{eq:XYZ}) to this inequality, we obtain that
$$
   \rho = -X-Y-Z \le -X-Y + \nu_1 -X - Y \le 2(|X|+|Y|) + \nu_1 \le \tfrac{2}{3}\,\rho + \nu_1,
$$
which implies that $\rho \le 3\nu_1$.
\end{proof}

Theorem \ref{thm:spectral_radius} is very similar to the corresponding theorem in \cite[Theorem 5.6]{GuMo} and gives the same constant $\tfrac{1}{3}$ in the lower bound. However, the proof here is essentially different and gives a new proof of \cite[Theorem 5.6]{GuMo}. On the other hand, the proof of \cite[Theorem 5.6]{GuMo} does not work for our Theorem \ref{thm:spectral_radius}.

While the factor $\tfrac{1}{3}$ in Theorem \ref{thm:spectral_radius} is tight for Hermitian matrices of the first kind \cite{GuMo}, it is not tight for the second kind, where it can be strengthened as follows.

\begin{theorem}
\label{thm:spectral_radius_new}
If $D$ is a digraph and $M=M(D)$ is its Hermitian matrix of the second kind, then
$$
   \tfrac{1}{2} \rho(M) \le \mu_1(D) \le \rho(M).
$$
\end{theorem}

\begin{proof}
We may assume that $\rho = \rho(M) = -\mu_n > \mu_1$.
Let $z\in \mathbb{C}^V$ be a unit eigenvector for $\mu_n$ and let $\widehat{z}\in \mathbb{R}^V$ be the vector with entries $\widehat{z}_v = |z_v|$, $v\in V$. Let us observe that $\Vert \widehat{z} \Vert = 1$ and that the following holds:
\begin{equation}\label{eq:half}
   \widehat{z}_v \widehat{z}_u = |\omega|\, |z_v|\, |z_u| = |\omega \bar{z}_v z_u| =
   \tfrac{1}{2}\,|\omega \bar{z}_v z_u| + \tfrac{1}{2}\,|\bar\omega \bar{z}_u z_v| \ge
   \tfrac{1}{2}\,|\omega \bar{z}_v z_u + \bar\omega \bar{z}_u z_v|.
\end{equation}
By using this inequality and (\ref{eq:zNz}), we obtain the following:
$$
  \tfrac{1}{2}\,\rho = \tfrac{1}{2}\,\left| z^* M z \right| = \tfrac{1}{2}\,\left| \sumvtou (\omega\bar{z}_v z_u + \bar\omega \bar{z}_u z_v)\right| \le
  \sumvtou \widehat{z}_v \widehat{z}_u = \widehat{z}^* M \widehat{z} \le \mu_1.
$$
This completes the proof.
\end{proof}

The factor $\tfrac{1}{2}$ in Theorem \ref{thm:spectral_radius_new} is best possible. This is justified with the example of the directed cycle of length 3, whose $M$-eigenvalues are $1$ (with multiplicity 2) and $-2$. By using the Cartesian product operation (see \cite{GuMo}), one can boost this example to obtain bigger graphs and larger spectral radii. The example of the 3-cycle shows another ``anomaly" that the largest eigenvalue of a (strongly) connected digraph need not be simple.

The above proof gives the following improvement for any $\alpha=a+bi$, where $|\a|=1$ with $a > \tfrac{1}{3}$:
$$
   a \, \rho(N^\a(D)) \le \nu_1(N^\a(D)) \le \rho(N^\a(D)).
$$

\subsection*{Acknowledgement}
The author is grateful to an anonymous referee whose insight provided a strengthening of Theorem \ref{thm:spectral_radius} that is now given in Theorem \ref{thm:spectral_radius_new}.

\bibliographystyle{abbrv}


\section{Introduction}

Let $\Gamma$ be a triangulation of a simplex $\Delta$ of dimension $d-1$.
The $h$-polynomial $h(\Gamma; x) = h_0 + h_1 x + \cdots + h_d x^d$ is a common and
convenient way of encoding the number of faces of $\Gamma$ in each dimension.
It is characterized by the equation
\[
\sum_{i=0}^{d} h_i (x+1)^{d-i} = \sum_{i=0}^d f_{i-1} x^{d-i},
\]
where $f_{-1} = 1$ and $f_i$ is the number of $i$-dimensional faces of $\Gamma$,
for $i \geq 0$.  The coefficients $h_i$ are non-negative integers.
One powerful tool for studying $h(\Gamma;x)$ is the local $h$-polynomial
$\ell(\Gamma;x) = \ell_0 + \ell_1x + \cdots + \ell_d x^d$,
introduced by Stanley in his seminal paper \cite{Stanley92}.
It is characterized via M\"obius inversion by the equation
\begin{equation} \label{eq:h=sum-of-l}
h(\Gamma;x) = \sum_{F \leq \Delta} \ell(\Gamma_F;x),
\end{equation}
where $\Gamma_F$ denotes the restriction of the triangulation $\Gamma$ to a face $F$
(which may be empty or all of $\Delta$), together with the condition $\ell(\varnothing; x) = 1$.

The local $h$-polynomial has remarkable properties.
In particular, the coefficients $\ell_i$ are nonnegative and satisfy $ \ell_i = \ell_{d-i}$.
Moreover, if the subdivision is regular, then these coefficients are unimodal.
Among other applications, Stanley used local $h$-polynomials to prove that $h$-polynomials
increase coefficientwise under refinement.

As discussed in Section 2, the local $h$-polynomial also behaves predictably with respect
to basic operations on subdivisions.  It is additive for refinements that nontrivially
subdivide only one facet, multiplicative for joins, and vanishes on the trivial subdivision.
In particular, if $\Gamma'$ is a refinement of $\Gamma$ that nontrivially subdivides only
one facet, and that subdivision is the cone over a subdivision of a codimension 1 face,
then $\ell(\Gamma';x) = \ell(\Gamma;x)$.  We call such subdivisions
\emph{conical facet refinements}.

Our aim is to study geometric triangulations with vanishing local $h$-polynomial.
One might hope that all such subdivisions are obtained from the trivial subdivision by a
sequence of conical facet refinements.
However, not all triangulations $\Gamma$ with vanishing local $h$-polynomial can be obtained from
the trivial subdivision in this way.  One notable example
in dimension 2 is the following subdivision, which we call the
\emph{triforce}\footnote{This name reflects the subdivision's realization in a sacred
golden relic that is the ultimate source of power in the action-adventure video game
series \emph{The Legend of Zelda}.}.

\medskip

\begin{center}
\begin{tikzpicture}[scale=2]
\draw (0.5,0.8)  -- (1,0) -- (0,0) -- (0.5,0.8);
\draw (0.75,0.4) -- (0.5,0) -- (0.25,0.4) -- (0.75,0.4);
\end{tikzpicture}
\end{center}

All of our theorems are for \emph{geometric triangulations}, and all triangulations in this paper are assumed to be geometric

\begin{theorem} \label{thm:dim2}
In dimension 2, any triangulation with vanishing local $h$-polynomial is obtained from either the trivial subdivision or the triforce by a sequence of conical facet refinements.
\end{theorem}

\noindent The two cases in the theorem are distinguished by the \emph{internal edge graph}, i.e., the union of the edges that meet the interior of the subdivided simplex, which also figures prominently in the proof.  For an iterated conical facet refinement of the triforce subdivision, the internal edge graph has Euler characteristic zero and contains no vertices of the original triangle.  For an iterated conical facet refinement of the trivial subdivision in dimension 2, the internal edge graph is a tree that contains exactly one of the vertices of the original triangle.

In dimension three, the structural classification is even simpler.

\begin{theorem} \label{thm:dim3}
In dimension 3, any triangulation with vanishing local $h$-polynomial is obtained from the trivial subdivision  by a sequence of conical facet refinements.
\end{theorem}

\noindent The proof again relies on an analysis of the internal edge graph, which in dimension three is a union of trees, each of which contains exactly one vertex supported on a face of codimension at least 2 in the original simplex.

This structure of the internal edge graph is similar in higher dimensions.  Note that $\ell_0 = 0$ for any triangulation of a nonempty simplex.

\begin{theorem}  \label{thm:internaledgegraph}
Let $\Gamma$ be a triangulation of a simplex of dimension at least 3 such that $\ell_1 = \ell_2 = 0$.  Then the internal edge graph of $\Gamma$ is a union of trees each of which contains exactly one vertex supported on a face of codimension at least 2.
\end{theorem}

\noindent However, the pattern of obtaining all triangulations whose local $h$-polynomials vanish from a finite collection of examples by iterated conical facet refinements does not continue in higher dimensions.  See Section~\ref{sec:higherdim}.

Our investigation into the structure of triangulations with vanishing local $h$-polynomial
is motivated by recent connections to algebraic and arithmetic geometry.
Local $h$-polynomials appear prominently in formulas for dimensions of homology groups of
intersection complexes for toric morphisms \cite{deCataldoMiglioriniMustata18} and
multiplicities of eigenvalues of monodromy \cite{KatzStapledon16, Stapledon17}.
In Igusa's $p$-adic monodromy conjecture \cite{Igusa75}, and the motivic generalization of
Denef and Loeser \cite{DenefLoeser98}, the essential question is understanding whether or
not these multiplicities vanish.  This problem is also natural and interesting from a
purely combinatorial viewpoint; Stanley specifically asked for a nice characterization of
such triangulations in his original paper \cite[Problem~4.13]{Stanley92}.

\medskip

\noindent \textbf{Acknowledgments.}  This project was carried out as part of the SUMRY undergraduate research program. We thank the other participants and mentors in this program for their support and encouragement.  We are also grateful to Matt Larson for helpful conversation related to this project.  Our work was partially supported by NSF grant DMS-1702428.

\section{Preliminaries}

We consider only geometric subdivisions $\Gamma$ of a $(d-1)$-simplex $\Delta$, except where explicitly stated otherwise.  In particular, all of the triangulations that we consider are realized by subdividing a linearly embedded simplex into subsimplices.  For further details and background on local $h$-polynomials, we refer the reader to the recent survey article of Athanasiadis \cite{Athanasiadis16}, as well as \cite{Stanley92}.

\subsection{Formulas for the local $h$-polynomial}
We recall two useful formulas for the local $h$-polynomial $\ell(\Gamma;x)$.  First, by applying M\"{o}bius inversion to \eqref{eq:h=sum-of-l}, we can express $\ell(\Gamma;x)$ as an alternating sum of $h$-polynomials:
\begin{equation} \label{eq:alternating sum}
\ell(\Gamma;x) = \sum_{F \leq \Delta}(-1)^{\codim (F)} h(\Gamma_F;x).
\end{equation}
Here, the codimension is $\codim(F) = d - 1 - \dim(F)$, since $\dim(\Delta) = d-1$.

Let $\sigma: \Gamma \rightarrow \Delta$ be the map taking a face $G \in \Gamma$ to the smallest face $F \leq \Delta$ that contains it.  One says that $G$ is \emph{carried by} $\sigma(G)$, and the \emph{excess} of $G$ is
\[
e(G) := \dim (\sigma(G)) - \dim (G).
\]
As in \cite[Prop. 2.2]{Stanley92}, we can then express the local $h$-polynomial as
\begin{equation}\label{eq:excess formula}
\ell(\Gamma; x) = \sum_{G \in \Gamma} (-1)^{\codim(G)} x^{d- e(G)} (x-1)^{e(G)}.
\end{equation}

\begin{definition}
We say that $G \in \Gamma$ is an \emph{interior face} if $\sigma(G) = \Delta$.
\end{definition}

\subsection{Elementary properties of the local $h$-polynomial}

The following basic properties of local $h$-polynomials were mentioned in the introduction and will be used in our main arguments.

Given triangulations $\Gamma$ and $\Gamma'$ of simplices $\Delta$ and $\Delta'$, respectively, the join $\Gamma * \Gamma'$ is naturally a triangulation of the simplex $\Delta * \Delta'$.  It is well-known that $h$-polynomials are multiplicative for joins, i.e.,
\[
h(\Gamma*\Gamma';x) = h(\Gamma;x) \cdot h(\Gamma';x).
\]
It follows that the local $h$-polynomial has the analogous property \cite[Lemma~2.2]{AthanasiadisSavvidou12}, i.e.,
\begin{equation} \label{eq:join}
\ell(\Gamma*\Gamma';x) = \ell(\Gamma;x) \cdot \ell(\Gamma';x).
\end{equation}

A triangulation $\Gamma'$ of $\Delta$ is a \emph{refinement} of $\Gamma$ if every $G \in \Gamma$ is a union of faces of $\Gamma'$.  We then write $\Gamma'_G$ for the restriction of $\Gamma'$ to $G$. We are particularly interested in the restriction of $\Gamma'$ to maximal faces, or \emph{facets}.

\begin{definition}
Let $\Gamma'$ be a refinement of $\Gamma$, and let $G$ be a facet of $\Gamma$.  We say that $\Gamma'$ is a \emph{facet refinement of $\Gamma$ along $G$} if, for any facet $H \neq G$, the restriction $\Gamma'_H$ is trivial.
\end{definition}

\begin{figure}[h!] \label{fig:facet-refinement}
{
\centering
\begin{tabular}{cc}
\begin{tikzpicture}[scale=2]
\draw[fill=white] (0,0) -- (1,0) -- (.5, .8) -- cycle;
\filldraw (0,0) circle (.6pt) -- (1,0) circle (.6pt) -- (.5, .8) circle (.6pt) -- (0,0);
\filldraw (.5,0) circle (.6pt) -- (.75, .4) circle (.6pt) -- (.25, .4) circle (.6pt) -- (.5,0);
\end{tikzpicture}
&
\begin{tikzpicture}[scale=2]
\draw[fill=white] (0,0) -- (1,0) -- (.5, .8) -- cycle;
\filldraw (0,0) circle (.6pt) -- (1,0) circle (.6pt) -- (.5, .8) circle (.6pt) -- (0,0);
\filldraw (.5,0) circle (.6pt) -- (.75, .4) circle (.6pt) -- (.25, .4) circle (.6pt) -- (.5,0);
\filldraw (.25, .4) -- (.25, 0) circle (.6pt);
\filldraw (.125, .2) circle (.6pt) -- (.5, 0);
\filldraw (0,0) -- (.25, .133) circle (.6pt);
\end{tikzpicture}
\\
$\Gamma$
&
$\Gamma'$
\end{tabular}\\
}
\caption{A facet refinement of the triforce along its lower left facet.}
\end{figure}

Local $h$-polynomials behave additively with respect to facet refinements.

\begin{proposition}\label{additivity of local h}
Let $\Gamma'$ be a facet refinement of $\Gamma$ along $G$, and let $\Gamma'_{G}$ denote the triangulation of $G$ induced by restricting $\Gamma'$.  Then
\[
\ell(\Gamma';x) = \ell(\Gamma;x)  + \ell(\Gamma'_G; x).
\]

\end{proposition}

\noindent This follows from \cite[Corollary~4.7]{KatzStapledon16}, a much more general result about local $h$-polynomials for compositions of strong formal subdivisions, applied to $\Gamma' \to \Gamma \to \Delta$.  For the reader's convenience, we include a direct proof in our setting.

\begin{proof}
The assertion is vacuously true in dimension zero.  We proceed by induction on dimension.  First, we observe that
\begin{equation}
\label{eq:additiveh} h(\Gamma';x) = h(\Gamma;x) + h(\Gamma'_G;x) - h(G;x).
\end{equation}
This follows from the formula $\sum_{i = 0}^d f_i (x-1)^{d-i} = \sum_{i=0}^d h_i x^{d-i}$, by inclusion-exclusion.

Next, write each $h$-polynomial in \eqref{eq:additiveh} as a sum of local $h$-polynomials, using \eqref{eq:h=sum-of-l}. The contributions of the empty faces cancel, and each nonempty face of $G$ contributes 0, since it is trivially subdivided. Therefore, we have
\begin{equation} \label{eq:sumofls}
\sum_{\{\varnothing\} \neq H \leq \Delta} \ell(\Gamma'_H;x) = \sum_{\{\varnothing\} \neq H \leq \Delta} \ell(\Gamma_H;x) + \sum_{\{\varnothing\} \neq H' \leq G} \ell(\Gamma'_{H'};x).
\end{equation}
For each nonempty face $H \leq \Delta$, either $\Gamma'_H = \Gamma_H$ or $\Gamma'_H$ is a facet refinement of $\Gamma_H$ along $G \cap H$.  If $H$ is a proper face of $\Delta$ then, by induction on dimension, we may assume $\ell(\Gamma'_H;x) = \ell(\Gamma_H;x) + \ell(\Gamma'_{G \cap H};x)$.  Similarly, for any nonempty face $H' \leq G$, either $\Gamma'_{H'}$ is the trivial subdivision, or $H'$ is carried by a face $H \leq \Delta$ such that  $H' = G \cap H$ and $\Gamma'_H$ is a facet refinement of $\Gamma_H$ along $H'$ . Thus, the contributions of proper faces in \eqref{eq:sumofls} cancel, and we conclude that $\ell(\Gamma';x) = \ell(\Gamma;x) + \ell(\Gamma'_G;x)$.
\end{proof}

\begin{definition}
A facet refinement $\Gamma'$ of $\Gamma$ along $G$ is a \emph{conical facet refinement} if $\Gamma'_G$ is the cone over $\Gamma'_H$ for some codimension 1 face $H < G$.
\end{definition}

\begin{corollary}  \label{cor:nfacetjoin}
Let $\Gamma'$ be a facet refinement of $\Gamma$ along $G$, and suppose that $\Gamma'_G$ is the join of two triangulations of faces of $G$, one of which has vanishing local $h$-polynomial.  Then $\ell(\Gamma';x) = \ell(\Gamma;x)$.  In particular, if $\Gamma'$ is a conical facet refinement of $\Gamma$, then $\ell(\Gamma';x) = \ell(\Gamma;x)$.
\end{corollary}

\begin{proof}
This follows immediately from \eqref{eq:join} and Proposition~\ref{additivity of local h}.
\end{proof}

\begin{figure}[h!] \label{fig:conical-facet-refinement}
{
\centering
\begin{tikzpicture}[scale=3.5]
\draw[fill=white] (0,0) -- (1,0) -- (.5, .8) -- cycle;
\filldraw (0,0) circle (.4pt) -- (1,0) circle (.4pt) -- (.5, .8) circle (.4pt) -- (0,0);
\filldraw (.5,0) circle (.4pt) -- (.75, .4) circle (.4pt) -- (.25, .4) circle (.4pt) -- (.5,0);

\draw [fill= black](.1,0) circle (.4pt) -- (.25,.4) circle (.4pt) -- (.2,0) circle (.4pt) -- (.25,.4) circle (.4pt) -- (.3,0) circle (.4pt) -- (.25,.4) circle (.4pt) -- (.4,0) circle (.4pt);
\draw [fill= black] (.45,.72) circle (.4pt) -- (.6,.64) circle (.4pt) -- (.4,.64) circle (.4pt) -- (.75,.4) circle (.4pt) -- (.3,.48) circle (.4pt);
\draw [fill= black] (.75,0) circle (.4pt) -- (.75,.4) circle (.4pt);
\draw [fill= black] (.75,0) -- (.85,.24) circle (.4pt);
\draw [fill= black] (.75,0) -- (.9,.16) circle (.4pt);
\end{tikzpicture}
\caption{A triangulation obtained from the triforce by a series of conical facet refinements.}
}
\end{figure}

\section{The internal edge graph}

The proofs of the main results of this paper all involve the \emph{internal edge graph}, formed by the edges of a triangulation $\Gamma$ that meet the interior of the simplex $\Delta$.  In this section, we study the properties of this graph when $d \geq 3$, and $\ell_1 = \ell_2 = 0$.

Let $f_i^j$ denote the number of $i$-simplices in $\Gamma$ that are carried by $j$-faces of $\Delta$.  From equation \eqref{eq:excess formula}, we see that $\ell_0 = 0$, $\ell_1 = f_0^{d-1}$, and
\begin{equation}\label{eqn:formula for l2}
\ell_2 = f_1^{d-1} -f_0^{d-2} - (d-1)f_0^{d-1}.
\end{equation}
In particular, if $\ell_1 = 0$ then $\ell_2 = f_1^{d-1} - f_0^{d-2}$.

Stanley proved non-negativity of all $\ell_i$, using methods from commutative algebra.
It follows that if there are no interior vertices (i.e., if $\ell_1 = 0$) then
$f_1^{d-1} \geq f_0^{d-2}$.  Note that every vertex carried by a $(d-2)$-face is contained
in an interior edge, so this inequality is a statement about the internal edge graph.
We now give a combinatorial proof of this inequality in a stronger form, showing in particular
that it holds separately on each connected component of the internal edge graph.

\begin{definition}
The \emph{internal edge graph} of a subdivision $\Gamma$ of $\Delta$ is the graph
$\Sigma(\Gamma)$ whose edges are the edges of $\Gamma$ carried by the improper face $\Delta$,
and whose vertices are the vertices incident to those edges.
\end{definition}

\noindent When $\Gamma$ is clear from context, we will write $\Sigma$ rather than $\Sigma(\Gamma)$.

\begin{proposition}\label{prop:each component has a low excess or low chi}
Let $\Gamma$ be a triangulation of a $(d - 1)$-simplex $\Delta$ with $d \geq 3$ and
$\ell_1 = 0$.
Then each connected component $C$ of $\Sigma(\Gamma)$ contains either a vertex of excess
less than $d-2$ or a simple $3$-cycle.

Furthermore, if $d \geq 4$, and $C$ has no vertices of excess less than $d-2$ then
$C$ contains at least two distinct simple $3$-cycles.
\end{proposition}

\begin{proof}
We may assume that $\Delta$ is the standard $(d-1)$-simplex in $\mathbb{R}^{d}$, i.e.,
the convex hull of the standard basis vectors $e_1,\ldots,e_{d}$. Let
$F_i = \Conv{e_1,\ldots,\widehat{e_i},\ldots,e_{d}}$, and let
$\pi_i : \mathbb{R}^d \to \mathbb{R}$ be the projection onto the $i^{\mathrm{th}}$ coordinate axis.

Suppose $C$ is a connected component of $\Sigma$ all of whose vertices have excess $d-2$.
We now show that $C$ contains a simple 3-cycle.

Consider the two codimension 1 faces carrying the endpoints of an edge of $C$. Without
loss of generality, we may assume these to be $F_1$ and $F_2$.  Let $L$ be the linear
functional $ \pi_1 + \pi_2$.  Let $e= \Conv{v_1, v_2}$ be an edge of $C$ such that
$v_1 \in F_1$, and $v_2 \in F_2$, such that $L(v_1 + v_2)$ is maximal among all such
edges.
For $(i,j) \in \{(1,2),(2,1)\}$, let $P_i : \mathbb{R}^{d}\to\mathbb{R}^d$ be the projection to the
linear hyperplane spanned by $F_i$ along $v_j-v_i$.
More explicitly, for all $w \in \mathbb{R}^{d}$,
\[
P_i(w) = w - \frac{\pi_i(w)}{\pi_i(v_j)}(v_j - v_i),
\]
Substituting in the definitions above, we observe that
\begin{align*}
\pi_i(v_j)L(P_i(w) - v_i) & = \pi_i(v_j)\pi_j(P_i(w)) - \pi_i(v_j)\pi_j(v_i) \\
& = \pi_i(v_j)\pi_j(w) + \pi_j(v_i)\pi_i(w) - \pi_i(v_j) \pi_j(v_i).
\end{align*}
This is invariant under switching $i$ and $j$, and hence
\begin{equation}\label{e:projection equality}
\pi_1(v_2)L(P_1(w) - v_1) = \pi_2(v_1)L(P_2(w) - v_2).
\end{equation}

Applying $P_1$ to the closed star $\Cl(\Star(e))$, we see that $e$ projects to $v_1$, and
2-faces project to segments ending at $v_1$. Since $e$ is an interior edge, $v_1$ is in
the relative interior of $P_1(\Cl(\Star(e)))$.  Therefore,
for any hyperplane $H$ in the image of $P_1$ that contains $v_1$, there is a vertex of
the link of $e$ that projects to either side of that hyperplane. In particular,
there is a vertex $w$ in the link of $e$ such that $L(P_1(w) - v_1) > 0$.
Recalling that all vertices of $C$ have excess $d-2$, we note that
$\pi_1(v_2)>0$ and $\pi_2(v_1)>0$, otherwise they would be contained in
$F_1\cap F_2$.
Therefore, by \eqref{e:projection equality}, $L(P_2(w) - v_2) > 0$.

We claim that $w$ is not carried by $F_1$ or $F_2$.
Indeed, if $w$ is carried by $F_1$, we would have $P_1(w) = w$ and
$L(v_2 + w) > L(v_2 + v_1)$, contradicting our maximality assumption on $e$. Thus
$w \notin F_1$. An identical argument shows that $w$ is not carried by $F_2$. This proves
the claim.  Then without loss of generality, we may assume that $w$ is carried by $F_3$.
The three interior edges of $C$ connecting $v_1$, $v_2,$ and $w$ form a simple 3-cycle.
If $d = 3$, then we are done.

Suppose $d \geq 4$.  We must show that $C$ contains another simple $3$-cycle.
Since $w$ is in $\lk e$, $\Gamma$ has an interior $2$-face $G = \Conv{v_1,v_2, w}$.
Let $p$ denote the projection along $G$ onto the codimension 2 linear subspace spanned by
$F_1 \cap F_2$. Then $p(G)$ is in the relative interior of $p(\Cl(\Star(G))$, and hence
there exists a vertex $w' \notin F_1 \cap F_2 \cap F_3$ in the link of $G$. Without loss
of generality, $w' \notin F_1$. Then $\Conv{v_1, w'}$ is an interior edge in $C$, and
hence $w'$ is carried by a codimension 1 face of $\Delta$. Since $w' \notin F_2 \cap F_3$,
we may assume, without loss of generality, that $w' \notin F_2$. Then the edges of $G$
and $\Conv{v_1, v_2, w'}$ form two distinct simple $3$-cycles in $C$.
\end{proof}

Let $f_i^j(C)$ denote the number of $i$-faces of $C$ that are carried by $j$-faces of
$\Delta$.

\begin{corollary} \label{cor:non-negative}
Suppose $d \geq 3$ and $\ell_1 = 0$.  Then each connected component $C$ of $\Sigma$ has
\[
f_0^{d-2}(C) \leq f_1^{d-1}(C).
\]
\end{corollary}

\begin{proof}
Recall that the Euler characteristic of a connected graph $C$ is
\[
\chi(C) = \# \{ \mbox{vertices of $C$}\} - \#\{ \mbox{edges of $C$} \}.
\]
This is also equal to $1 - h^1(C)$, where $h^1$ denotes the first Betti number.  Hence $\chi(C) \leq 1$, with equality if and only if $C$ is a tree.  Moreover, $\chi(C) = 0$ if and only if $C$ contains a unique simple cycle.

Let $C$ be a connected component of the internal edge graph.  Then $$f_0^{d-2}(C) - f_1^{d-1}(C) = \chi(C) - f_0^{< d-2}(C).$$  By Proposition~\ref{prop:each component has a low excess or low chi}, either $\chi(C) = 1$ and $f_0^{< d-2}$ is positive, or $\chi(C) \leq 0$.
\end{proof}

In Section~\ref{sec:dim23}, we will repeatedly use the following structural description of the internal edge graph.

\begin{proposition}\label{prop:components are trees or have cycles}
Suppose $d \geq 3$, $\ell_1=\ell_2 = 0$, and $C$ is a connected component of $\Sigma$.  Then either
\begin{enumerate}
\item $C$ is a tree with a unique vertex of excess less than $d-2$, or
\item $d=3$, every vertex of $C$ has excess $d-2$, and $C$ has a unique simple 3-cycle.
\end{enumerate}
\end{proposition}

\begin{proof}
Since $\ell_1 = 0$, we can write $\ell_2$ as a sum over connected components of $\Sigma$
\begin{align*}
\ell_2 & = \sum_C f_1^{d-1}(C) - f_0^{d-2}(C),\\
& = \sum_C f_0^{< d-2}(C) - \chi(C).
\end{align*}
By Proposition~\ref{prop:each component has a low excess or low chi} and Corollary~\ref{cor:non-negative}, each summand is nonnegative, and the summand corresponding to $C$ is zero if and only if either $C$ is a tree with a unique vertex of excess less than $d-2$, or $d=3$, every vertex of $C$ has excess $d-2$, and $C$ contains a unique simple 3-cycle.
\end{proof}

Note that Theorem \ref{thm:internaledgegraph} follows immediately from Proposition \ref{prop:components are trees or have cycles}.

\begin{remark} \label{rem:connected}
Note that the internal edge graph of any triangulation of a $2$-simplex is connected. To see this, consider a polyhedral subdivision of a triangle with disconnected internal edge graph.  It must contain a polygon that meets two different components of the internal edge graph, and such a polygon has at least four edges: at least two segments of the boundary, and at least one edge from each component of the internal edge graph that it meets. In particular, it is not a triangulation.
\end{remark}

\begin{definition}
Let $F$ be a facet in a triangulation of the simplex $\Delta$.
We say that $F$ is a \emph{pyramid} if there is some proper face of $\Delta$ that
contains every vertex of $F$ except one.
\end{definition}

\begin{remark} \label{rem:pyramids}
Understanding when a facet of a triangulation is a pyramid is helpful for potential
applications, e.g., to the monodromy conjecture for nondegenerate hypersurfaces, as
mentioned in the introduction.  Note that, if $\Gamma'$ is a conical facet refinement of
$\Gamma$ along $G$, then both $G$ and every facet of $\Gamma'_G$ is such a pyramid. Hence,
it follows from Theorem~\ref{thm:dim3} that, if $d = 4$ and $\ell(\Gamma;x) = 0$, then every facet is
a pyramid.  Similarly, it follows from Theorem~\ref{thm:dim2} that, if $d = 3$ and
$\ell(\Gamma;x) = 0$, then at most one facet is not a pyramid.
\end{remark}

The following proposition is not used in the remainder of the paper, but illustrates how
the structure of the internal edge graph gives useful statements about which facets of a
triangulation with vanishing local $h$-polynomial are pyramids.

\begin{proposition} \label{prop:pyramid}
Let $\Gamma$ be a triangulation of a $(d - 1)$-simplex $\Delta$ with $d \ge 4$. If $\ell_1= \ell_2 = 0$,
then every facet $G$ of $\Gamma$ that contains an interior edge is a pyramid.
\end{proposition}

\begin{proof}
Let $\Gamma$ be a facet of a triangulation of a $(d - 1)$-simplex $\Delta$ with $d \ge 4$, such that $\ell_1= \ell_2 = 0$.  Let $G$ be a facet of $\Gamma$ that contains an interior edge.  We assume that $G$ is not a pyramid, and will derive a contradiction. Let $e_1 = \Conv{v_1, v_2}$ be an interior edge in $G$. Let $C$ be the connected component of the interior edge graph $\Sigma(\Gamma)$ containing $e_1$.

By Proposition~\ref{prop:components are trees or have cycles}, either $v_1$ or $v_2$ is carried by a codimension 1 face.
Without
loss of generality, $v_1$ is carried by a codimension 1 face $F_1$. Note that $v_2 \notin F_1$. Since $G$ is not a pyramid, it has a vertex $w_2 \notin F_1$ distinct from $v_2$. Then $e_2 = \Conv{v_1, w_2}$ lies in $C$.
By Proposition~\ref{prop:components are trees or have cycles}, either $v_2$ or $w_2$ is carried by a codimension 1 face. Without
loss of generality, $v_2$ is carried by a codimension 1 face $F_2$. If $w_2 \notin F_2$, then $\Conv{v_2, w_2}$
is an interior edge and $C$ contains a cycle, contradicting Proposition~\ref{prop:components are trees or have cycles}.
Hence, $w_2 \in F_2$. Since $G$ is not a pyramid, it has a vertex $w_1 \notin F_2$ distinct from $v_1$.
Then $e_3 = \Conv{w_1, v_2}$ lies in $C$. If $w_1 \notin F_1$, then $\Conv{v_1, w_1}$
is an interior edge and  $C$ contains a cycle, contradicting Proposition~\ref{prop:components are trees or have cycles}.
Hence $w_1 \in F_1$. By Proposition~\ref{prop:components are trees or have cycles},
either $w_1$ or $w_2$ is carried by a codimension 1 face, meaning that they cannot both
be contained in some third codimension 1 face $F_3$. Therefore,
$e_4 = \Conv{w_1, w_2}$ lies in $C$, and $e_1,e_3,e_4,e_2$ forms a
cycle in $C$, contradicting Proposition~\ref{prop:components are trees or have cycles}.
\end{proof}

\section{Dimensions 2 and 3} \label{sec:dim23}

In this section we prove our main structural results, that all subdivisions with vanishing local $h$ can be obtained by a sequence of conical facet refinements from the trivial subdivision and the triforce subdivison in dimension 2, and from the trivial subdivision in dimension 3.  For both results our proof is by induction on the number of vertices or interior edges in the subdivision.  In the induction step, we identify a subcomplex that arises from a conical subdivision of a facet in a coarser subdivision.

Recall that the support of a subcomplex $\Gamma' \subset \Gamma$ is the union of the faces in $\Gamma'$.  It is a closed subset of $\Delta$.
The relative boundary of any closed subset $F \subset \Delta$ is the intersection of $F$ with the closure of its complement $\Delta \smallsetminus F$.  Note, in particular, that if $F$ is the support of a subcomplex $\Gamma' \subset \Gamma$, then the relative boundary of $F$ is the support of a subcomplex of $\Gamma'$.  When no confusion seems possible, we will refer to this subcomplex as the relative boundary of $\Gamma'$.

\begin{lemma}\label{coarsening-condition}
Let $\Gamma$ be a triangulation of a simplex $\Delta$ and let $\Gamma'\subseteq \Gamma$ be a subcomplex whose support is a simplex $F$ of dimension $d-1$.  Suppose that $\Gamma'$ induces the trivial subdivision on the relative boundary of $F$ in $\Delta$.
Then $\ell(\Gamma; x) - \ell(\Gamma'; x)$ has nonnegative coefficients.
In particular, if $\ell(\Gamma;x)=0$, then $\ell(\Gamma';x)=0$.
\end{lemma}

\begin{proof}
Since $\Gamma'$ does not nontrivially subdivide the relative boundary of $F$ in $\Delta$, we get a subdivision $\Gamma''$ of $\Delta$ by replacing $\Gamma'$ with $F$ itself. Note that $\Gamma$ is then a facet refinement of $\Gamma''$ along $F$, so by Proposition~\ref{additivity of local h},
we have $\ell(\Gamma;x)=\ell(\Gamma'';x)+\ell(\Gamma';x)$, and the lemma follows.
\end{proof}

\begin{remark}
Note that, in the proof of the lemma, $\Gamma$ is a conical facet refinement of $\Gamma''$  (resp. obtained from $\Gamma''$ by a sequence of conical facet refinements) if and only $\Gamma'$ is a conical facet refinement of $F$ (resp. obtained from the trivial subdivision of $F$ by a sequence of conical facet refinements).
\end{remark}

\subsection{Dimension 2}
We now apply the lemma to prove Theorem~\ref{thm:dim2}, which says that any triangulation of a 2-dimensional simplex with vanishing local $h$-polynomial is obtained from either the trivial subdivision or the triforce by a sequence of conical facet refinements.

\begin{proof}[Proof of Theorem~\ref{thm:dim2}]
If $\Gamma$ has no interior edges, then it is the trivial subdivision; this case of the theorem is obvious.  We proceed by induction on the number of interior edges, and consider two cases, according to the possible structures of the internal edge graph given by  Proposition~\ref{prop:components are trees or have cycles}.  (Recall also that the internal edge graph is connected, by Remark~\ref{rem:connected}.)

\proofsection{Case 1:}{$\Sigma(\Gamma)$ is a tree that contains exactly one vertex of excess zero}  We claim that $\Gamma$ is obtained from the trivial subdivision of $\Delta$ by a sequence of conical facet refinements.  To see this, note that any interior edge that contains a vertex of excess zero divides $\Gamma$ into two triangular subcomplexes whose relative boundary is not subdivided, as shown.
\begin{center}
\begin{tikzpicture}
\filldraw (0,0) circle (.6pt) -- (.5,0) circle (.6pt)-- (1,0)circle (.6pt) -- (.5,.8)circle (.6pt) -- (0,0);
\filldraw (.5,0)circle (.6pt) -- (.5,.8) ;
\node at (.5,-.2) {};
\node at (.5, 1) {};
\end{tikzpicture}
\end{center}
These two subcomplexes have vanishing local $h$-polynomials, by Lemma~\ref{coarsening-condition}. By induction, each is obtained from the trivial subdivision by a sequence of conical facet refinements.  Hence $\Gamma$ is obtained from the subdivision into two triangles by a sequence of conical facet refinements.  The subdivision into two triangles is itself a conical facet refinement of the trivial subdivision, and we conclude that $\Gamma$ is obtained from the trivial subdivision by a sequence of conical facet refinements, as claimed.

\proofsection{Case 2:}{$\Sigma(\Gamma)$ contains a simple 3-cycle and all of its vertices have excess 1} We claim that $\Gamma$ is obtained from the triforce by a sequence of conical facet refinements.  To see this, note that the three vertices of the simple 3-cycle must be carried by the three sides of $\Delta$.  This splits $\Gamma$ into four
triangular subcomplexes, each with unsubdivided relative boundary, i.e., $\Gamma$ is a refinement of the triforce, and the induced subdivision on the relative boundary of each of the four triangles is trivial.  Hence, by Lemma \ref{coarsening-condition}, each of the four induced subdivisions has vanishing local $h$-polynomial.  In fact, the induced subdivision  of the interior triangle must be trivial (because $\Gamma$ has no interior vertices).  In each of the other three triangles, the induced subdivision cannot contain a cycle of interior edges, and hence the induced subdivisions are obtained from the trivial subdivision by a sequence of conical facet refinements.  We conclude that $\Gamma$ is obtained from the triforce by a sequence of conical facet refinements, as claimed.
\end{proof}

\subsection{Dimension 3}

We now apply similar inductive arguments to prove Theorem~\ref{thm:dim3}, which says that any subdivision $\Gamma$ of a 3-dimensional simplex $\Delta$ with vanishing local $h$-polynomial is obtained from the trivial subdivision by a sequence of conical facet refinements.

\begin{proof}[Proof of Theorem~\ref{thm:dim3}]
If $\Gamma$ has only 4 vertices, then it is the trivial subdivision and the conclusion is obvious.  We proceed by induction on the number of vertices.  By Proposition~\ref{prop:components are trees or have cycles}, each connected component of the internal edge graph $\Sigma(\Gamma)$ is a tree with exactly one vertex of excess less than 2.

\proofsection{Case 1:}{$\Sigma(\Gamma)$ is empty}

Let $F$ be a 2-face of $\Delta$.  We claim that $\ell(\Gamma_F;x) = 0$.  To see this, note that any vertex carried by $F$ must be contained in an interior edge. Hence $\Gamma_F$ has no interior vertices, and $\ell_1(\Gamma_F;x) = 0$.  Now $\Gamma_F$ is a triangulation of the 2-simplex, so the symmetry of the local $h$-polynomial implies that $\ell_i(\Gamma_G;x) = \ell_{3-i}(\Gamma_F;x)$. Hence $\ell_2(\Gamma_F;x) = \ell_1(\Gamma_F;x) = 0$, and we conclude that $\ell(\Gamma_F;x) = 0$, as claimed.

\medskip

We proceed to consider three subcases, according to the internal edge graphs of the restrictions of $\Gamma$ to the 2-faces of $\Delta$.

\proofsection{Subcase 1.1:}{All 2-faces have empty internal edge graph}

Any interior 2-face of $\Gamma$ contains an edge that is carried by $\Gamma_F$ for some 2-face $F \leq \Delta$.  So $\Gamma$ has no interior 2-faces, and hence is the trivial subdivision.

\proofsection{Subcase 1.2:}{Some 2-face has a 3-cycle in its internal edge graph}

Let $F_1$ be a 2-face of $\Delta$ such that $\Gamma_{F_1}$ contains a 3-cycle of interior edges.
The vertices of this 3-cycle, $v_2$, $v_3$, and
$v_4$, are carried by distinct edges of $F_1$.
Label the other 2-faces of $\Delta$ as $F_2,F_3,F_4$ so that $v_i$ is carried by $F_1\cap F_i$.

The 2-face $\Conv{v_2,v_3,v_4}$ is contained in some 3-face $T=\Conv{v_2,v_3,v_4,x}$ of $\Gamma$. Since $\Sigma(\Gamma)$
is empty, the edges $\Conv{x,v_i}$ must be contained in the boundary of $\Delta$. Hence $x$ must lie in the intersection $F_2\cap F_3\cap F_4$, which is the vertex of $\Delta$ opposite $F_1$. Then the 2-faces $\Conv{v_1,v_2,x}$,
$\Conv{v_2,v_3,x}$ and $\Conv{v_1,v_3,x}$ cut $\Delta$ into four tetrahedral regions (the cone over a triforce subdivision of $F_1$).  By Lemma \ref{coarsening-condition}, the restriction of $\Gamma$ to each of these regions has vanishing local $h$-polynomial.  By induction, the induced subdivision of each of these regions is obtained from the trivial subdivision by a sequence of conical facet refinements, and we are done.

\proofsection{Subcase 1.3:}{No 2-face has a 3-cycle, and some 2-face has an interior edge}

Let $F$ be a $2$-face of $\Delta$ such that $\Gamma_F$ contains an interior edge.  By Proposition~\ref{prop:components are trees or have cycles}, there is an edge $e$ carried by $F$ with a vertex of excess 0.  Label the vertices of $\Delta$ as $A$, $B$, $C$, and $D$, so that $F = \Conv{B,C,D}$ and $e$ contains $D$.  Let $E \in \Conv{B,C}$ be the second vertex of $e$.

Let $G$ be an interior 2-face of $\Gamma$ containing $e$. Then $G=\Conv{e,v}$ for some vertex $v$ of $\Gamma$. This vertex $v$ must lie on $\Conv{A,B}$ or $\Conv{A,C}$, so that the edges $\Conv{E,v}$ and $\Conv{D,v}$ are not interior.  Furthermore, $v\ne B,C$, since $G$ is interior.  Then either $v=A$, or, without loss of generality, we may assume $\sigma(v)=\Conv{A,C}$.  These possibilities are illustrated in Figure~\ref{fig:twocases}.

\begin{figure}[h!]
\centering
\hfill
\begin{subfigure}{0.3\textwidth}
\centering
\begin{tikzpicture}[scale=3]
\coordinate (A) at (0,0.5); \node[below left=0.5ex of A] (Alab) {$A$};
\coordinate (B) at (1,0.4); \node[below right=0.5ex of B] (Blab) {$B$};
\coordinate (C) at (0.3,0); \node[below right=0.5ex of C] (Clab) {$C$};
\coordinate (D) at (0.4,1); \node[above=0.5ex of D] (Dlab) {$D$};
\path (B) -- (C) coordinate[midway] (E);
\node[below right=0.5ex of E] (Elab) {$E$};
\draw[dotted,color=white!50!black]
(A)
-- (B);
\draw
(A)
-- (C)
-- (B);
\draw
(B)
-- (D)
-- (C);
\fill[color=white!30!black,opacity=0.4]
(A)
-- (E)
-- (D);
\draw
(A)
-- (D);
\draw
(E)
-- (D) coordinate[pos=0.4] (ePos);
\draw[dashed,color=white!30!black]
(A)
-- (E);
\node[right=0.3ex of ePos] (e) {$e$};
\end{tikzpicture}
\caption{$v=A$}
\label{fig:v-equals-a}
\end{subfigure}
\hfill
\begin{subfigure}{0.3\textwidth}
\centering
\begin{tikzpicture}[scale=3]
\coordinate (A) at (0,0.5); \node[below left=0.5ex of A] (Alab) {$A$};
\coordinate (B) at (1,0.4); \node[below right=0.5ex of B] (Blab) {$B$};
\coordinate (C) at (0.3,0); \node[below right=0.5ex of C] (Clab) {$C$};
\coordinate (D) at (0.4,1); \node[above=0.5ex of D] (Dlab) {$D$};
\path (B) -- (C) coordinate[midway] (E);
\node[below right=0.5ex of E] (Elab) {$E$};
\path (A) -- (C) coordinate[pos=0.5] (V1);
\node[below left=0.5ex of V1] (V1lab) {$v$};
\draw[dotted,color=white!50!black]
(A)
-- (B);
\draw
(A)
-- (C)
-- (B);
\draw
(B)
-- (D)
-- (C);
\fill[color=white!30!black,opacity=0.4]
(V1)
-- (E)
-- (D);
\draw
(V1)
--(D);
\draw
(A)
-- (D);
\draw
(E)
-- (D) coordinate[pos=0.4] (ePos);
\draw[dashed,color=white!30!black]
(V1)
-- (E);
\node[right=0.3ex of ePos] (e) {$e$};
\end{tikzpicture}
\caption{$\sigma(v)=\Conv{A,C}$}
\label{fig:v-not-equal-a}
\end{subfigure}
\hfill {}
\caption{Possible locations for $v$ in Subcase 1.3}
\label{fig:loc-for-v}
\label{fig:twocases}
\end{figure}

If $v = A$ then $G$, marked in grey in Figure \ref{fig:v-equals-a}, divides $\Delta$ into two distinct tetrahedra, that satisfy the conditions of the Lemma \ref{coarsening-condition}, and by the usual induction argument we are done.

We may therefore assume $\sigma(v) = \Conv{A,C}$ as in Figure \ref{fig:v-not-equal-a}. In this case, $G$ cuts off a tetrahedron $\Conv{v,C,D,E}$, which satisfies the conditions of Lemma \ref{coarsening-condition}. By induction, it suffices to consider the case where this tetrahedron is trivially subdivided. Also since $G$ is an interior 2-face, it must be contained in another tetrahedron $\Conv{v,D,E,w}$ of $\Gamma$. We must have $w\in \Conv{v,A}\cup\Conv{A,B}\cup \Conv{E,B}\setminus\set{v,E}$, in order to prevent the edges $\Conv{w,v}$ and $\Conv{w,E}$ from being interior, and so that $w$ is on the other side of $G$ from $C$.

Now if $\sigma(w)=\Conv{A,B}$, then $\Gamma_{\Conv{A,B,C}}$ will contain a cycle of interior edges, with vertices $v$, $E$, $w$. But we assumed that no 2-face has a cycle in its internal edge graph. Thus, $w\in\Conv{v,A}\setminus \set{v}$ or $w\in\Conv{E,B}\setminus\set{E}$. Without loss of generality, we may assume $w\in\Conv{E,B}\setminus\set{E}$. Then regardless of whether $w=B$ or is instead
in the interior of the segment, the tetrahedron $\Conv{v,D,w,C}$ is one to which Lemma \ref{coarsening-condition} applies. Moreover, $\Gamma_{\Conv{v,D,w,C}}$ is nontrivial, because $G$ lies in the interior. Hence $\Gamma_{\Conv{v,D,w,C}}$ is obtained from the trivial subdivision by a series of conical facet refinements, and the conclusion follows by induction.

\bigskip

\proofsection{Case 2:}{$\Sigma(\Gamma)$ is nonempty}

By Proposition~\ref{prop:components are trees or have cycles}, we know
that each connected component of $\Sigma(\Gamma)$ is a tree rooted at a unique vertex of excess less than $2$. We consider one of these components $C$ and let $v_1$ be one of its leaves at maximal distance from the root. Let $w$ be the unique vertex adjacent to $v_1$ in $C$, and let $F$ be the 2-face carrying $v_1$. Let $v_2, \cdots, v_r$ be the other leaves that are adjacent to $w$ and carried by $F$, as in Figure \ref{fig:possible-sigma}.

\begin{figure}[h!]
\centering
\begin{tikzpicture}
\filldraw (0,0) circle (.6pt)-- (0,-1)circle (.6pt) -- (-1,-2) circle (.6pt)-- (-2, -3)circle (.6pt);
\filldraw (-1.5, -3) circle (.6pt) -- (-1, -2) circle (.6pt)-- (-1, -3)circle (.6pt);
\filldraw (-.5, -3) circle (.6pt)--( -1, -2);
\filldraw (0,-1) -- (1, -2) circle (.6pt) -- (1, -3)circle (.6pt);
\filldraw (-1,-2) -- (0,-3) circle (.6 pt);
\node at (0, .3) {root};
\node at (-1.3, -2) {$w$};
\node at (-2, -3.2) {$v_1$};
\node at (-1.25, -3.2) {$\cdots$};
\node at (-.5, -3.2) {$v_r$};
\draw[decoration={brace},decorate](-.4,-3.4) -- (-2.1,-3.4) ;
\node at (-1.25,-3.8) {carried by $F$};
\end{tikzpicture}
\caption{Structure of $C$.}
\label{fig:possible-sigma}
\end{figure}

Let $X$ be the subcomplex of $\Gamma_F$ that is the union of the closed stars of $v_1, \ldots, v_r$.

\proofsection{Claim:} {$X \subseteq \Cl(\Star(w))$.}

In other words, we claim that $\Gamma$ contains, as a subcomplex, the cone over $X$ with vertex $w$. To prove this, it suffices to show that if $G$ is a $k$-face of $\Gamma_F$ containing some vertex $v_i$, then $G$ is contained in a face of $\Star(w)$. Any such $G$ must be contained in an interior $(k+1)$-face of $\Gamma$.  Let $w'$ be a vertex such that $\Conv{G,w'}$ is an interior $(k+1)$-face. Note that $\Conv{v_i,w'}$ is an interior edge. Since $v_i$ is a leaf of the internal edge graph, we have $w'=w$.  Hence $\Conv{G,w}$ is a face of $\Star(w)$ that contains $G$.  This proves the claim.

\proofsection{Subcase 2.1:} {There is a vertex of $X$ carried by $F$ that is not in $\{ v_1, \ldots, v_r \}$}

Let $v$ be a vertex of $X$ that is carried by $F$ and not contained in $\{v_1, \ldots, v_r\}$. Since $v\in X$, there is an edge $\Conv{v,v_i}$ in $X$ for some index $i$. By the claim proved above, $\Conv{v,v_i,w}$ is a face of $\Gamma$, and $\Conv{v,w}$ is an edge of $C$. Since $v \not \in \{v_1, \ldots, v_r\}$, it must be the parent of $w$.

Consider the link $\lk_{\Gamma_F}(v)$.  This is a cycle of edges in $\Gamma_F$ that contains $v_i$.  If every edge in $\lk_{\Gamma_F}(v)$ contains some $v_j \in \{ v_1, \ldots, v_r \}$, then $\Star(v)$ is contained in $\Cl(\Star(w))$ and hence $\Conv{v,w}$ is the unique interior edge containing $v$.  This contradicts the fact that $v$ is the parent of $w$.  Hence $\lk_{\Gamma_F}(v)$ contains multiple vertices not in $\{v_1, \ldots, v_r \}$.

Choose vertices $w_1$ and $w_2$ in $\lk_{\Gamma_F}(v)$ that are not in $\{v_1, \ldots, v_r\}$ so that every vertex in the path between $w_1$ and $w_2$ that contains $v_i$ is in $\{v_1, \ldots, v_r\}$. Note that, since $w_1$ and $w_2$ are adjacent to $w$ and not in $\{v_1, \ldots, v_r \}$, and since $v$ is the parent of $w$ in $C$, the edges $\Conv{w,w_i}$ must not be interior.  Similarly, $w$ cannot be the root of $C$. Hence $w$ is carried by a $2$-face that contains both $w_1$ and $w_2$.

It follows that $w_1$ and $w_2$ lie on an edge of $F$. Thus the subcomplex with support $\Conv{v,w_1,w_2,w}$ satisfies the conditions of Lemma \ref{coarsening-condition}.

\proofsection{Subcase 2.1.1:}{$\Conv{v,w_1,w_2,w}$ is trivially subdivided}
We claim that this is impossible.  Indeed, if $\Conv{v,w_1,w_2,w}$ is trivially subdivided, then $\lk_{\Gamma_F}(v)$ consists of the edge $\Conv{w_1, w_2}$ together with a sequence of edges connecting $w_1$ to $w_2$ in which all other vertices are contained in $\{ v_1, \ldots, v_r \}$.  Then every facet containing $v$ is a pyramid with base on $F$ and vertex at $w$.  Hence, $v$ is a leaf of $C$.  This is impossible, since we have already shown that $v$ is the parent of $w$. This proves the claim.

\proofsection{Subcase 2.1.2:}{$\Conv{v,w_1,w_2,w}$ is nontrivially subdivided}
In this case, by the induction hypothesis, $\Gamma_{\Conv{v,w_1,w_2,w}}$ is obtained from the trivial subdivision by a series of conical facet refinements, and the result follows by induction.

\proofsection{Subcase 2.2:}{Every vertex of $X$ carried by $F$ is in $\{v_1, \ldots, v_r\}$}

In this case, every edge in the boundary of $X$ connects two vertices that are in the boundary of $F$.  Then the closure $X'$ of any connected component of the relative interior of $X$ is a subcomplex with support a polygon whose relative boundary in $F$ is trivially subdivided. If the support of $X'$ is a triangle, then we can apply Lemma~\ref{coarsening-condition} for the subcomplex of $\Gamma$ given by the cone over $X'$ with vertex $w$, and be done by induction.

If $X'$ is not a triangle, then we show that $X'$ is contained in a larger subcomplex $Y \subset \Gamma_F$ such that
\begin{enumerate}
\item The relative boundary of $Y$ in $F$ is trivially subdivided.
\item The cone over $Y$ with vertex $w$ is a subcomplex of $\Gamma$.
\item The support of $Y$ is a triangle.
\end{enumerate}
Once we find such a $Y$, then the conclusion of the theorem follows by Lemma~\ref{coarsening-condition} and induction on the number of vertices.

In order to find $Y$ as above, we begin with $Y_0=X'$, and construct an increasing sequence of subcomplexes
\[
Y_0 \subset Y_1 \subset \cdots \subset Y_n = Y
\]
such that each $Y_i$ has properties $(1)$ and $(2)$, and the support of $Y$ is a triangle.

We describe the construction of $Y_{i+1}$ from $Y_{i}$ in terms of corner cutting edges.

\begin{definition}
A boundary edge $e$ of $Y_i$ is \emph{corner cutting} if its vertices are carried by edges of $F$, and $Y_i$ lies on the side of $e$ opposite from the corner that it cuts off.
\end{definition}

In our construction, if $Y_i$ is not a triangle then $Y_{i+1}$ has strictly fewer corner cutting edges than $Y_i$.  Since $Y_0$ has at most $3$ corner cutting edges, the procedure will terminate for some $n \leq 3$.

\begin{figure}[h!]
\begin{center}
\begin{tikzpicture}[scale = 3]
\filldraw [color=white!70!black] [decorate, decoration = {snake}] (.4,0) to [out = 110, in = 220] (.65,.56);
\filldraw[color=white!70!black] (.4,0) -- (.6,0) -- (.8, 0.32) -- (.65, .56) -- (.4, 0) ;
\node at (.6, .24) {$Y_i$};
\draw (0,0)  -- (1,0) circle (.6pt) -- (.5,.8) -- (0,0);
\filldraw (1,0) circle (.6pt);
\filldraw (.6,0) circle (.6pt) -- (.8, 0.32) circle (.6pt) ;
\node at (1.1,-.2) {$v$};
\node at (.6, -.2) {$w_2$};
\node at (.9, .4) {$w_1$};
\node at (.76, .1) {$e$};
\end{tikzpicture}
\end{center}
\caption{A corner cutting edge $e$ with vertices $w_1$ and $w_2$}
\label{fig:ccedge}
\end{figure}

\proofsection{Subcase 2.2.1:}{$Y_i$ has no corner cutting edges} In this case, we claim that the support of $Y_i$ is triangular.  First, if $Y_i$ has no boundary edges, then its support is $F$. Suppose $Y_i$ has some boundary edge $e$. If both vertices of $e$ are carried by edges, then $Y_i$ is the on the same side as the vertex of $F$ cut off by $e$.  Any other boundary edge of $Y_i$ would have to be corner cutting.  Hence, the support of $Y_i$ is the triangle cut off by $e$, and we are done.

Otherwise, $e$ contains a vertex of $F$, and hence divides $F$ into two triangles.  Once again, any other boundary edge would have to be corner cutting, so the support of $Y_i$ is a triangle.  This proves the claim, and completes this subcase.

\bigskip

It remains to consider the situation where $Y_i$ has a corner cutting edge $e$, with vertices $w_1$ and $w_2$. Let $v$ be the vertex of $F$ cut off by $e$, as in Figure \ref{fig:ccedge}.  We consider two further subcases, according to whether or not $w$ is the root of $C$, as illustrated in Figure~\ref{fig:locations-for-w}.

\begin{figure}[h!]
\begin{center}
\hfill
\begin{subfigure}{0.3\textwidth}
 \begin{tikzpicture}[scale=2]
 \coordinate (v_0) at (1.5,0);
 \coordinate (v_a) at (.25,.5);
 \coordinate (v_b) at (-.5,-.5);
 \coordinate (top) at (.5, 1.5);
 \coordinate (w) at (0.25,0.25);

\path  (v_0) -- (v_a) coordinate[pos=.4] (w_1);
\path (v_0) -- (v_b) coordinate[pos=.4] (w_2);
\path (v_0) --(v_b) coordinate [pos=.8] (w_3);

\fill [white!70!black] (w_2)--(w_3)--(w_1) --(w_2);

\draw[dotted,color=white!30!black] (v_0) -- (v_a) ;
\draw  [thick] (v_0) -- (v_b) ;
\draw[dashed,color=white!30!black] (w_1) -- (w_2);
\draw[dotted,color=white!30!black] (v_a) -- (v_b);
\draw [thick] (top) -- (v_0);
\draw  [thick] (top) -- (v_b);
\draw[dotted,color=white!30!black] (top) -- (v_a);
\draw[dashed,color=white!30!black] (w) -- (w_1);
\draw (w) -- (w_3);
\draw[dashed,color=white!30!black] (w_1) -- (w_3);
\draw (w) -- (w_2);

\filldraw (v_0) circle (.6pt);
\filldraw (w_1) circle (.6pt);
\filldraw (w_2) circle (.6pt);
\filldraw (w) circle (.6pt);
\filldraw (w_3) circle (.6pt);

\node[below right=0.5ex of v_0]  {$v$};
\node[above right  =0.2ex of w_1]  {$w_1$};
\node[below =0.5ex of w_2]  {$w_2$};
\node[below =0.5ex of w]  {$w$};
\node[below =0.5ex of w_3]  {$w_3$};
\end{tikzpicture}
\caption{$w$ carried by a $2$-face}
\label{subfig:w3triangle}
\end{subfigure}
\hfill
\begin{subfigure}{0.3\textwidth}
 \begin{tikzpicture}[scale=2]
 \coordinate (v_0) at (1.5,0);
 \coordinate (v_a) at (.25,.5);
 \coordinate (v_b) at (-.5,-.5);
 \coordinate (top) at (.5, 1.5);

\path  (v_0) -- (v_a) coordinate[pos=.4] (w_1);
\path  (v_0) -- (v_a) coordinate[pos=.5] (z_1);
\path (v_0) -- (v_b) coordinate[pos=.4] (w_2);
\path (v_0) -- (v_b) coordinate[pos=.5] (z_2);
\path (v_0) -- (top) coordinate[pos=.6] (w);

\filldraw [color=white!70!black] [decorate, decoration = {snake}] (z_2) to [out = 110, in = 220] (z_1);
\filldraw [color=white!70!black] (w_1) to (z_1) to  (z_2) to (w_2);

\draw[dotted,color=white!30!black] (v_0) -- (v_a) ;
\draw  [thick] (v_0) -- (v_b) ;
\draw[dashed,color=white!30!black] (w_1) -- (w_2);
\draw[dotted,color=white!30!black] (v_a) -- (v_b);
\draw [thick] (top) -- (v_0);
\draw  [thick] (top) -- (v_b);
\draw[dotted,color=white!30!black] (top) -- (v_a);
\draw[dashed,color=white!30!black] (w) -- (w_1);
\draw (w) -- (w_2);

\filldraw (v_0) circle (.6pt);
\filldraw (w_1) circle (.6pt);
\filldraw (w_2) circle (.6pt);
\filldraw (w) circle (.6pt);

\node[below right=0.5ex of v_0]  {$v$};
\node[above right  =0.1ex of w_1]  {$w_1$};
\node[below =0.5ex of w_2]  {$w_2$};
\node[right =0.5ex of w]  {$w$};

\end{tikzpicture}
\caption{$w$ contained in $F_1\cap F_2$}
\label{subfig:w-on-edge}
\end{subfigure}
\hfill{}
\caption{Possible locations for $w$, with the region $Y_i$ shaded}
\label{fig:locations-for-w}
\end{center}
\end{figure}

\proofsection{Subcase 2.2.2:}{$Y_i$ has a corner cutting edge, and $w$ is not the root of $C$}

If $w$ is not the root, then it is carried by some 2-face of $\Delta$ other than $F$. It cannot be carried by the face opposite $v$, since then both of the edges $\Conv{w_1,w}$ and $\Conv{w_2,w}$ would be interior, meaning that $w_1$ and $w_2$ would both be low excess vertices on the same component of the internal edge graph. Therefore without loss of generality, we can assume that $w$ is carried by the 2-face $F' \neq F$ containing $w_2$.

Note that the support of $Y_i$ is a polygon, and if we now think of its boundary not in $F$, but in the plane containing $F$, there are exactly two vertices of the boundary adjacent to $w_1$; one of these is $w_2$, and the other is some vertex, which we call $w_3$.

The edge $\Conv{w,w_3}$ cannot be interior, since $w_3$ is contained in the boundary of $F$, and hence has low excess, but $w_1$ is the root of $C$. Thus $w_3 \in F\cap F'$, as shown in Figure \ref{subfig:w3triangle}. It follows that the polygon $Y_i$ both contains and is contained in the triangle bounded by $\Conv{w_1, w_2}$,  $\Conv{w_2, w_3}$, and $\Conv{w_1, w_3}$.  Therefore, $Y_i$ is the triangle $\Conv{w_1, w_2, w_3}$, and we are done.

\proofsection{Subcase 2.2.3:}{$Y_i$ has a corner cutting edge, and $w$ is the root of $C$}

In this case, $w$ cannot have an interior edge to either $w_1$ or $w_2$, since $w_1$ and $w_2$ both have low excess, and $w$ is the root of $C$. Therefore, if $F_1$ and $F_2$ are the faces other than $F$ containing $w_1$ and $w_2$ respectively, then $w$ must lie on the edge $F_1\cap F_2$, possibly at the vertex, though this will not matter. This situation is illustrated in Figure \ref{subfig:w-on-edge}. Since $\Conv{w,w_1,w_2}$ is in fact a 2-face of $\Gamma$, the subcomplex with support $\Conv{w,w_1,w_2,v}$ satisfies the hypotheses of Lemma \ref{coarsening-condition}. Thus we may assume that it is trivially subdivided. We then set $Y_{i+1}$ to be the union of $Y_i$ with the triangle $\Conv{w_1, w_2, v}$.  Then $Y_{i+1}$ satisfies properties (1) and (2), and has fewer corner cutting edges than $Y_i$.  Repeating this procedure, as $Y_0$ has at most three corner cutting edges, we eventually arrive at $Y_n$ satisfying properties (1)-(3), and the theorem follows.
\end{proof}

\section{Higher dimensions} \label{sec:higherdim}

While Theorem~\ref{thm:internaledgegraph}, which is an immediate consequence of Proposition~\ref{prop:components are trees or have cycles}, says that the structure of the internal edge graph is essentially the same in higher dimensions, it seems that this is not sufficient for a useful classification of triangulations with vanishing local $h$-polynomial.  Indeed, there are infinitely many different triangulations of the 4-simplex with empty internal edge graph.  For instance, the join of the triforce with an arbitrary subdivision of the 1-simplex has this property.  These triangulations give rise to the following negative result.

\begin{proposition}
There is no finite collection of triangulations of the $4$-simplex with vanishing local $h$-polynomials from which all others can be obtained by a series of conical facet refinements.
\end{proposition}

\begin{proof}
Let $\Gamma_n$ be the subdivision of the $1$-simplex with $n$ interior vertices, and let $\Gamma'$ be the triforce.  If $e$ is an edge of $\Gamma_n$ such that both endpoints are interior, then the join of $e$ with the center facet of $\Gamma'$ is not a pyramid.  Hence $\Gamma_n * \Gamma'$ is a triangulation with vanishing local $h$-polynomial that has $n-1$ non-pyramid facets.  The proposition follows, since, as discussed in Remark~\ref{rem:pyramids}, all of the new facets introduced by a conical facet refinement are pyramids.
\end{proof}

\noindent We leave the problem of classifying triangulations with vanishing local $h$-polynomial in dimensions 4 and higher open, for future research.


\section{Introduction}
\label{sec:introduction}

Identifying the dark matter is the issue of greatest concern both in particle physics and cosmology.
Amongst various models,
cosmological constant and cold dark matter ($\Lambda$CDM) model is the most favored by many observations~\cite{Komatsu:2010fb, Aghanim:2018eyx}.

$\Lambda$CDM model significantly succeeds in explaining large-scale structure of the universe though,
it still has astronomical problems in small scales.
One of the famous problems is the core-cusp problem,
in which the tension between the simulated cusp profile~\cite{Navarro_1997, Fukushige:1996nr, Ishiyama:2011af}
and observed cored profile of the galactic center~\cite{Borriello:2000rv, Gilmore:2007fy, Oh:2008ww} is reported.
This problem is solved, for example,
when dark matter is ultra-light ($m \sim 10^{-22}\ {\rm eV}$)
and coherently oscillating~\cite{Hu:2000ke, Hui:2016ltb}.
This is because the large de Broglie wavelength of the ultra-light particle, approximately the same scale as the galactic center ($\sim {\rm kpc}$),
smears out the central excessive density.
This particle is often called ultra-light axion-like particle (ULAP)
because it is ultra-light and coherently oscillating like QCD axion~\cite{Weinberg:1977ma, Wilczek:1977pj, Dine:1981rt}.

In this case,
clumps of ULAP can be formed by the gravity interaction like axion star~\cite{Kolb:1993zz}
(we call it "ULA star" with an analogy to QCD axion).
ULA star is a kind of soliton
that the attractive forces from gravity and repulsive forces from the kinetic pressure are balanced.
Some studies show
that ULAP forms ULA star by gravity interaction at the galactic center~\cite{Schive:2014dra, Schive:2014hza, Veltmaat:2018dfz}.
The existence of the star solution is also proved in the contest of axion~\cite{Braaten:2015eeu, Visinelli:2017ooc}
(the axion star formation is also confirmed in Ref.~\cite{Levkov:2018kau}).

Before the ULA star formation, however,
self-interaction of ULAP can lead to the oscillon~\cite{Bogolyubsky:1976yu, Gleiser:1993pt, Copeland:1995fq} formation by parametric resonance~\cite{Shtanov:1994ce, Kofman:1995fi, Kofman:1997yn,Yoshimura:1995gc} in the radiation dominated universe.
Oscillon is a spatially localized solitonic state of a real scalar field,
which often survives for a long time because of the conservation of the adiabatic invariance~\cite{Kasuya:2002zs,Kawasaki:2015vga}.
The existence of such dense objects in the universe
affects the current observational constraint~\cite{Li:2013nal, Nori:2018pka, Porayko:2018sfa}
and proposed experiments~\cite{Graham:2017ivz, Nagano:2019rbw}.

The possibility of long lifetime oscillons in a ULAP potential is pointed out in Ref.~\cite{Olle:2019kbo},
but their formation and precise lifetime is unclear
because analytical estimation of the fluctuation growth and lifetime is difficult.
In this paper,
we confirm the oscillon formation by classical lattice simulation in a particular parameter region.
We also calculate the oscillon lifetime by the analytical method derived in Ref.~\cite{Ibe:2019vyo},
which shows that produced ULAP oscillons live long, maybe up to the present day.

In Sec.~\ref{sec:oscillon},
we show that the condition for the oscillon formation is satisfied in a typical ULAP potential.
In Sec.~\ref{sec:formation},
we numerically confirm the oscillon formation by the classical lattice simulation.
In Sec.~\ref{sec:lifetime},
we analytically derive the lifetime of produced oscillon and briefly discuss the cosmological consequence.
Finally, in Sec.~\ref{sec:conclusion},
we conclude our results.

\section{Oscillon of Ultra-Light Axion Particle}
\label{sec:oscillon}

In the context of ULAP, we often use the cosine potential with an analogy to axion.
In that case, however, large initial fluctuations are necessary for the oscillon formation
because the instability generated by the cosine potential is too weak to enhance the fluctuations.
In addition,
the resultant oscillon is unstable
(the lifetime is $\mathcal{O}(100) m$)~\cite{Piette:1997hf}
and hardly leave interesting cosmological effects~\cite{Vaquero:2018tib, Buschmann:2019icd}.

Still many other potentials of ULAP are suggested~\cite{Dong:2010in, Kallosh:2013hoa, Kallosh:2013tua, Kallosh:2013yoa},
instead of the cosine potential we use the following potential~\cite{Silverstein:2008sg, McAllister:2008hb, Nomura:2017ehb}
\begin{equation}
    V(\phi) = \frac{m^2F^2}{2p} \left[1 - \left(1+\frac{\phi^2}{F^2}\right)^{-p} \right],
    \hspace{5mm} (p > -1)
    \label{eq:potential}
\end{equation}
where $p > -1$ is necessary to guarantee the existence of the oscillon solution
(see \ref{sec:solution} and Eq.~(\ref{eq:oscillon-condition})).
In the following subsections,
we will show that the two general conditions necessary for the oscillon formation are satisfied;
\begin{enumerate}
   \item The existence of the oscillon solution,
   \item The large fluctuations $\delta\phi /\phi_0 \sim \mathcal{O}(1)$.
\end{enumerate}
Note that these conditions are not sufficient but just necessary.
To confirm the oscillon formation,
we should perform classical lattice simulation.

\subsection{Oscillon solution}
\label{sec:solution}

First, let us show that the oscillon configuration can be realized in the potential Eq.~(\ref{eq:potential}).

Oscillon is the non-topological pseudo soliton
approximately conserving the adiabatic invariant $I$~\cite{Kasuya:2002zs} defined as
\begin{equation}
    I \equiv \frac{1}{\omega} \int d^{3}x \overline{\dot{\phi}^{2}},
\end{equation}
where the overline represents the time average over periodic motion of $\phi$
and $\omega$ does its oscillation frequency almost the same as the mass of $\phi$.

The oscillon configuration is obtained by minimizing the time-averaged energy $\overline{E}$
for a given $I$.
With the use of the Lagrangian multiplier method,
\begin{align}
    E_{\lambda}
    &= \overline{E}
    + \lambda \left( I - \frac{1}{\omega} \int dx^{3} \overline{\dot{\phi}^2} \right), \\
    &= \int d^{3}x \left[ \left(1-\frac{2\lambda}{\omega}\right)\frac{1}{2}\overline{\dot{\phi}^2}
    + \frac{1}{2} \overline{(\nabla\phi)^2} + \overline{V(\phi)} \right] + \lambda I.
    \label{eq:LM-energy}
\end{align}
Assuming that the periodic motion of $\phi$ is written as
$\phi(x) \simeq \Phi(\bm{x}) \cos \omega t$
and defining $V(\Phi)\equiv2\overline{V(\phi)}$, Eq.~(\ref{eq:LM-energy}) becomes
\begin{align}
    E_{\omega}
    &= \frac{1}{2} \int d^{3}x
    \left[\omega \left(\omega-2\lambda\right) \frac{1}{2}\Phi^{2}
    +\frac{1}{2} (\nabla\Phi)^{2} + V(\Phi) \right] + \lambda I, \\
    & =\frac{1}{2} \int d^{3}x
    \left[\frac{1}{2}(\nabla\Phi)^{2}+V(\Phi)-\frac{1}{2}\omega^2\Phi^{2} \right]
    + \omega I.
\end{align}
where we use the relation $\lambda = \omega$~\cite{Ibe:2019vyo} in the last line.
Because the lowest energy configuration is realized in the spherically symmetric configuration,
we set $\Phi(\bm{x})=\Phi(r)$ and impose the boundary condition as
\begin{equation}
\left.\Phi\right|_{r\rightarrow\infty}=0,\ \ \ \left.\frac{d\Phi}{dr}\right|_{r\rightarrow0}=0.\label{eq:boundary condition}
\end{equation}
Under this condition,
by differentiating $E_{\omega}$ by $\Phi$ to get the extremum we obtain
\begin{equation}
    \frac{d^{2}\Phi}{dr^{2}}
    + \frac{2}{r} \frac{d\Phi}{dr}
    + \frac{d}{d\Phi} \left( \frac{1}{2}\omega^2\Phi^{2}-V(\Phi) \right)
    = 0.
    \label{eq:oscillon-eq}
\end{equation}
This equation is considered as the equations of motion of $\Phi$
moving in the potential $\omega^2\Phi^{2}/2-V(\Phi)$ with friction.
Thus the condition for existing the solution satisfying the boundary
condition Eq.~(\ref{eq:boundary condition}) is
\begin{equation}
    \min \left[ \frac{V(\Phi)}{\Phi^{2}}\right] < \omega^2 < m^2.
    \label{eq:oscillon-condition}
\end{equation}
Because the ULAP potential Eq.~(\ref{eq:potential}) of $p > -1$ satisfies Eq.~(\ref{eq:oscillon-condition}),
the oscillon configuration is realized when the ULAP fluctuations are large enough.

\subsection{Instability growth}
\label{sec:Floquet}

Not only the existence of the oscillon solution but also the large spatial inhomogeneity is necessary for the oscillon formation.
In the ULAP case, however, the field value is almost uniform at the beginning of the oscillation
because the homogeneous initial value is set by inflation.
Therefore, we discuss the growth of the small initial fluctuations by Floquet analysis in this subsection.

For analytical purposes,
we expand the potential Eq.~(\ref{eq:potential}) up to a quartic order under $\phi < F$,
\begin{equation}
    V(\phi) \simeq
    \frac{m^2F^2}{2} \left( \frac{\phi^2}{F^2}
    - \frac{p+1}{2} \frac{\phi^4}{F^4} \right).
\end{equation}

Assuming that the background is dominated by the harmonic oscillation $\phi_0 \simeq \Phi \cos (mt)$,
the equation of motion of the $k$-mode fluctuations without cosmic expansion leads to the Mathieu equation,
\begin{align}
    \ddot{\delta \phi_k} + \left[k^2 + V''(\phi_0) \right] \delta \phi_k
    &= 0,\\
    \Leftrightarrow\ \ \
    \ddot{\delta \phi_k}
    + \left[ k^2 + m^2
    - \frac{3(p+1)}{2} m^2 \left(\frac{\Phi}{F}\right)^2
    \left( 1 + \cos 2mt \right)
    \right] \delta \phi_k
    &\simeq 0.
\end{align}
When $\Phi \lesssim F$,
the narrow instability band is induced with the approximate growth rate (Floquet exponent)
\begin{equation}
    \frac{\mu_{\rm max}}{m} \simeq \frac{3(p+1)}{8} \left(\frac{\Phi}{F}\right)^2,
    \hspace{0.5cm} {\rm at} \hspace{0.4cm}
    \frac{k}{m} \simeq \sqrt{\frac{3(p+1)}{2}} \frac{\Phi}{F},
    \label{eq:Floquet_exponent}
\end{equation}
which shows that the larger $p$ leads to the stronger instability
\footnote{
    See Ref.~\cite{Olle:2019kbo} for Floquet exponent of other $p$.
}.
In an expanding universe, when $\mu_{\rm max}/H \gtrsim 1$,
the fluctuation growth beats the cosmic expansion
and it can lead to large enough fluctuations for the oscillon formation.

\section{Oscillon Formation}
\label{sec:formation}

In this section,
we verify the oscillon production by the classical lattice simulation
after a brief explanation of the simulation setup.
The oscillon formation in the potential Eq.~(\ref{eq:potential}) is confirmed in many papers in the inflationary context~\cite{Amin:2011hj, Lozanov:2017hjm, Kitajima:2018zco, Hong:2017ooe}.
Here, we test the formation in the radiation dominated universe.

\subsection{Simulation setup}
\label{sec:setup}

In the simulation,
the units of the field, time, space, etc. are taken to be $F$ and $m^{-1}$, that is,
\begin{equation}
    \bar{\phi} \equiv \frac{\phi}{F},\ \
    \bar{\tau} \equiv m\tau,\ \
    \bar{x} \equiv mx,\ \ \dots\ {\rm etc}.
\end{equation}
where the overline denotes the dimensionless program variables and $\tau$ is the conformal time.

The equation of motion of $\bar{\phi}$ is represented by
\begin{equation}
    \bar{\phi}'' + 2 \frac{a'}{a}\bar{\phi}'
    - \bar{\Delta} \bar{\phi}
    + a^2\frac{\partial \overline{V}}{\partial \bar{\phi}} = 0,
    \hspace{0.6cm}
    \overline{V}
    = \frac{1}{2p} \left[ 1 - \left(1+\bar{\phi}^2\right)^{-p} \right].
    \label{eq:num_eom}
\end{equation}
where dash denotes the derivative of $\overline{\tau}$.
As the initial condition,
\begin{itemize}
    \item Hubble parameter:
        ULAP starts to oscillate in the radiation dominated universe, so we take $H_i = 1/2t = m$.
    \item Scale factor:
        the initial scale factor is set to be unity $a_i = 1$.
        (in the radiation dominated universe $a = \bar{\tau}$).
    \item Field values:
        the initial field value and its derivative are set as

        \begin{eqnarray}
            \bar{\phi}_i = \frac{2}{3}\pi (1 + \zeta),\ \ \bar{\phi}'_i = 0.
        \end{eqnarray}

        where $\zeta$ is the uniform random fluctuations of $\mathcal{O}(10^{-5})$
        \footnote{
            The initial field average $\langle \bar{\phi}_i \rangle = 2\pi/3$ is just the benchmark.
        }.
\end{itemize}
The other simulation parameters are shown in Table~\ref{Ta:params}.

\begin{table}[t]
    \centering
    \begin{tabular}{cc}
        \hline \hline
        $p$ & varying \tabularnewline
        Initial field value $\langle \bar{\phi}_i \rangle$ & $2\pi/3$ \tabularnewline
        Box size $L$  & $8$ \tabularnewline
        Grid size $N$ & $256^3$ \tabularnewline
        Final time    & $50$ \tabularnewline
        Time step     & $0.01$ \tabularnewline
        \hline \hline
    \end{tabular}
    \caption{
        Simulation parameters.
        $p$ is changed in every simulation to set the appropriate potential.
    }
    \label{Ta:params}
\end{table}

We utilize our lattice simulation code used in Refs.~\cite{Ibe:2019vyo, Ibe:2019lzv},
in which the time evolution is calculated by the fourth-order symplectic integration scheme
and the spatial derivatives are calculated by the fourth-order central difference scheme.
We impose the periodic boundary condition on the boundary.

\subsection{Simulation results}
\label{results}

\begin{figure}[t]
    \begin{minipage}{.48\linewidth}
        \begin{center}
            \includegraphics[width=\linewidth]{figs/result_p=25.pdf}
        \end{center}
    \end{minipage}
    \begin{minipage}{.48\linewidth}
        \begin{center}
            \includegraphics[width=\linewidth]{figs/result_p=4.pdf}
        \end{center}
    \end{minipage}
    \caption{
        Simulation results.
        We plot the ratio of the energy density $\rho$ to the averaged one $\bar{\rho}$ in the simulation box at $a = 30$.
        Left and right figures show the case of $p = 2.5$ and $4$ respectively.
        From these figures,
        We can find that a lot of oscillons are produced by self-resonance.
    }
    \label{fig:sim_results}
\end{figure}

The simulation results are shown in Fig.~\ref{fig:sim_results} for $p = 2.5$ (left) and $p = 4$ (right).

In both cases, the high energy local objects are produced.
However, we have confirmed by simulations that the oscillon formation does not occur
when $p \lesssim 2$.
There are two main reasons for this.
The first is that smaller $p$ leads to the weaker resonance as shown in Eq.~(\ref{eq:Floquet_exponent}).
The twice larger $\mu$ after few oscillations makes fluctuations $e^{10} \sim 10^4$  times larger.
The second is that
the effective mass at the beginning of the oscillation $m_{\rm eff}^2 = V''(\phi_i)$ becomes smaller in larger $p$.
This causes the later onset of the oscillation,
which leads to the increase of the number of oscillations in each Hubble time
and instability enhancement.

We also comment on the initial field value.
In this simulation, we set $\langle \bar{\phi}_i \rangle = 2\pi/3$ as the benchmark of the $\mathcal{O}(1)$ initial field value.
However, if we assume the highly fine-tuned initial condition as
$\langle \bar{\phi}_i \rangle = \mathcal{O}(10^{-1})$ or so,
fluctuations may not be grown enough because of the opposite reasons to the above two.
On the other hand,
if we assume the larger initial value,
even smaller $p$ may cause the oscillon formation.

Although we focused only on the potential Eq.~(\ref{eq:potential}) in this paper,
other ULAP potentials satisfying the condition Eq.~(\ref{eq:oscillon-condition}) also have the possibility of the oscillon formation.

\section{Oscillon Lifetime}
\label{sec:lifetime}

In this section, we briefly estimate the oscillon lifetime based on Ref.~\cite{Ibe:2019vyo}
and discuss the effect of the resultant oscillon.

In Ref.~\cite{Ibe:2019vyo},
we derive the method to calculate the classical decay rate of oscillon
by solving the equation of motion of fluctuations around the theoretical oscillon profile.
Let us decompose a scalar field $\phi(x)$ into the oscillon profile and fluctuations as
\begin{equation}
    \phi(x) = \psi(r) \cos \omega t + \xi(x).
\end{equation}
where $\psi(r)$ obeys the oscillon profile equation Eq.~(\ref{eq:oscillon-eq})
\footnote{
    We impose the spherical symmetry on the system
    because we focus only on the lowest energy state.
}.
With the use of the equation of motion of $\phi$ and Eq.~(\ref{eq:oscillon-eq}),
the equation of motion of $\xi$ is solved and
the energy loss rate of oscillon
\begin{equation}
    \frac{dE}{dt} = 4\pi r^2 T^{0r},
\end{equation}
where $T_{0r} = \partial_0 \xi \partial_r \xi$ denotes the Poynting vector
can be calculated.

\begin{figure}[t]
    \begin{minipage}{.45\linewidth}
        \begin{center}
            \includegraphics[width=\linewidth]{figs/lifetime_p=25.pdf}
        \end{center}
    \end{minipage}
    \hspace{.3cm}
    \begin{minipage}{.45\linewidth}
        \begin{center}
            \includegraphics[width=\linewidth]{figs/lifetime_p=4.pdf}
        \end{center}
    \end{minipage}
    \caption{
        The decay rate of oscillon.
        The left and right figures show the case of $p = 2.5, 4$ respectively.
        For $\omega \gtrsim 0.98$ the oscillon does not satisfy the stability condition $\frac{\partial \omega}{\partial I} > 0$ \cite{Mukaida:2016hwd}.
        Larger $\omega$ means the smaller oscillon by definition,
        so the minimum lifetime is $\mathcal{O}(10^7)m$.
        In addition, if larger osillons (e.g., $\omega \lesssim 0.96$) are formed,
        the lifetime may be much longer
        because the decay rate becomes much smaller there.
    }
    \label{fig:lifetime}
\end{figure}

The derived decay rate $\Gamma = |\dot{E}|/E$ is shown in Fig.~\ref{fig:lifetime}.
Since the larger $\omega$ corresponds to the smaller oscillon (see Sec.~\ref{sec:solution} for the definition of $\omega$), the decay rate of a large oscillon evolves following the curve in Fig.~\ref{fig:lifetime} from left to right.
From Fig.~\ref{fig:lifetime} it is found that the produced oscillon is stable at least for $\mathcal{O}(10^7)m$.
For the typical ULAP mass $m \sim 10^{-22}\ {\rm eV}$,
the minimum lifetime is about $10^8$ years.
Thus, they may affect the cosmic evolution, particularly structure formation.

In addition,
the lifetime of oscillon may become longer
because its lifetime depends on its size.
As shown in the figures, the decay rate has a lot of poles
where the decay rate is extremely small.
At these points, oscillon hardly evolves to the smaller one,
which results in a longer lifetime ($\gtrsim 10^{10}$ years).
Therefore, if the large oscillons are formed,
ULAP oscillons may still live in the present universe.

\section{Conclusions}
\label{sec:conclusion}

In this paper, we have examined the oscillon formation in a ULAP potential, Eq.~(\ref{eq:potential}).
We have shown that
oscillons are really produced when the potential index is $p > 2$ with $\langle \phi_i \rangle/F = 2\pi/3$.
Their lifetime is at least $\mathcal{O}(10^7)m$, which equals to $10^8$ years with $m \sim 10^{-22}\ {\rm eV}$.
Because such long lifetime objects survive until the structure formation,
they may affect dynamical history of the universe.

Moreover, the lifetime could be much longer if large oscillons are formed
because the decay rate of the oscillon extremely decreases in the specific profiles.
In either case,
a deep understanding of the evolution of ULAP oscillons is indispensable to consider effects on cosmology.
Detailed study should be discussed in the future work.

\begin{acknowledgments}

This work was supported by JSPS KAKENHI Grant Nos. 17H01131 (M.K.) and 17K05434 (M.K.), MEXT KAKENHI Grant Nos. 15H05889 (M.K.), World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan, and JSPS Research Fellowships for Young Scientists Grant No. 19J12936 (E.S.).

\end{acknowledgments}

\bibliographystyle{JHEP}