\section{Introduction and statement of the results}

Given an open bounded subset $\Omega$ of $\mathbb{R}^n$,
we consider the following eigenvalue problem
\begin{equation}\label{f:pb1} \begin{cases}
 - \Delta^N_p u = \lambda_p u & \text{ in } \Omega
 \\  \noalign{\medskip}
 u = 0 & \text{ on } \partial \Omega\,,
 \end{cases}
 \end{equation}
where $\Delta^N_p$ denotes the normalized or game-theoretic $p$-Laplacian, defined for any $p \in (1, + \infty)$ by
\[
\begin{array}{ll}
\Delta^N_p u & \displaystyle := \frac{1}{p} |\nabla u| ^ {p-2}  {\rm div} \big (|\nabla u| ^ {p-2} \nabla u \big )
\\ \noalign{\medskip}
& \displaystyle
= \frac{p-2}{p} |\nabla u|^{-2} \pscal{\nabla ^2 u\, \nabla u}{\nabla u} + \frac{1}{p}\, \Trace(\nabla ^2 u)\,,
\end{array}
\]
 where $\nabla^2u$  stands for the Hessian of $u$.
Equivalently, see \cite{K0}, it can be defined as a convex combination of the limit operators
as $p \to 1$ and $p \to + \infty$, since
\begin{equation}\label{convexcomb}\Delta^N_p u =\frac{p-1}{p} \Delta _\infty ^N u +  \frac{1}{p}  \Delta _1 ^N u  \,,
\end{equation}
with
$$
\Delta _\infty ^N u = \frac{1}{|\nabla u| ^ 2} \pscal{\nabla ^ 2 u \nabla u} {\nabla u}\
\hbox{ and } \,
\Delta _1 ^N u := |\nabla u|\, {\rm div} \Big ( \frac{\nabla u} {|\nabla u|} \Big )\, . $$
Let us point out that solutions to \eqref{f:pb1} are in general not classical, i.e. of class $C^2$, but have to be understood as viscosity solutions and these are defined in Section \ref{sec:proofs}.

The normalized $p$-Laplacian has recently received increasing attention, partly because of its application in image processing \cite{K0,Does} and in the description of tug-of-war games (see \cite{PSSW1, PSSW2}).
Without claiming to be complete we list \cite{BK18,CFd,CFe,CFf,CF7,EKNT, JK, Juut07, KH, K11, kuhn , MPR1, MPR2} for some related works.

 Following Berestycki, Nirenberg, and Varadhan \cite{BNV},
in the paper \cite{BiDe2006} (where actually they deal with a wider class of operators),
Birindelli and Demengel introduced  the  {\it first eigenvalue of $\Delta^N_p$ in $\Omega$} as
 $$\overline \lambda _p (\Omega) := \sup \Big \{ \lambda_p \in \mathbb{R} \ :\  \exists u >0 \text{ such that } \
 \Delta^N_p u + \lambda _p u \leq 0\ \text{ in the viscosity sense}
 \Big \}\,.$$
They proved that calling it first eigenvalue is justified, see \cite[Theorems 1.3 and 1.4]{BiDe2006}.
In particular they showed that there exists a positive eigenfunction associated with
$\overline \lambda _p (\Omega)$. In other words for $\lambda _p = \overline \lambda _p (\Omega)$
problem \eqref{f:pb1} admits a positive viscosity solution.
They also posed the open problem to determine whether $\overline \lambda _p (\Omega)$ is  simple.
We show that the answer is affirmative.  More precisely, we prove:

\begin{theorem}\label{t:simple}
Let $\Omega$ be an open bounded domain in $\mathbb{R} ^n$, with $\partial \Omega$ smooth and connected.
If $u$ and $v$ are two positive eigenfunctions associated with
$\overline \lambda _p (\Omega)$, then $u$ and $v$ are proportional, that is there exists $t \in \mathbb{R} _+$ such that
$u = tv $ in $\Omega$.
 \end{theorem}

Here and in the following, $\partial\Omega$ smooth means that
it is of class $C^{2,\alpha}$.
Theorem \ref{t:simple} has the following  immediate consequence:

  \begin{corollary}\label{symmetry}
 Let $\Omega$ be an open bounded domain in $\mathbb{R} ^n$, with $\partial \Omega$ smooth and connected. If $\Omega$ is invariant under elements from a symmetry group such as reflections or rotations, then so are the first eigenfunctions of the normalized $p$-Laplace operator.
 \end{corollary}

In order to obtain Theorem \ref{t:simple} we follow the approach used by Sakaguchi in \cite{Sak}.
In particular, it will be clear by inspection of the proof that this method does not work if one drops the assumption that
$\partial \Omega$ is connected. It is conceivable that the result continues to be true for more general domains, as it is known in the literature for other kinds of operators at least in dimension two  (see for instance \cite[Theorem 4.1]{BirDem2010}).

 As a fundamental preliminary tool, our proof of Theorem \ref{t:simple} exploits a Hopf type lemma (see Lemma \ref{t:hopf})  and, incidentally, it requires also the strict positivity of the eigenvalue. The latter can be easily established by comparison with the behaviour on balls (see Lemma \ref{l:ball} and Lemma \ref{l:positivity}).
  In fact, the observation that $\overline\lambda_p(\Omega_1)\geq\overline\lambda_p(\Omega_2)$ for $\Omega_2\subset\Omega_1$ leads to the bounds
 \begin{equation}\label{ballcomp}
\overline\lambda_p(B_R)\leq \overline\lambda_p(\Omega)\leq\overline\lambda_p(B_\rho),
\end{equation}
where $\rho$ and $R$ denote inradius and outer radius of $\Omega$, see the recent papers \cite{blanc, KH}.
These bounds are sharp if $\Omega$ is a ball, but they are far from optimal if $R-r$ becomes large, e.g. for slender ellipsoids. On the other hand, the problem of finding more accurate bounds for the eigenvalue seems to be an interesting and mostly unexplored question.   In this respect \eqref{ballcomp}  is complemented by the following lower estimate for $\overline \lambda _p (\Omega)$ in terms of the Lebesgue measure of $\Omega$.

\begin{theorem}\label{t:FK}
 For every open bounded  domain $\Omega$ in $\mathbb{R} ^n$   we have the lower bound
$$\overline \lambda _p (\Omega) \geq K _{n, p} |\Omega| ^ { - 2/n} \, , $$
with
\begin{equation}\label{f:Knp}
K_{n,p}:= \frac{ \big ( n[(p-1) \wedge 1 ] ) ^ 2 } {p(p-1) }   \,
4^{-1+1/n}\, \pi^{1+1/n}\, \Gamma\left(\frac{n+1}{2}\right)^{-2/n}\,.
\end{equation}

\end{theorem}

The proof of Theorem \ref{t:FK}  will be obtained by the Alexandrov--Bakelman--Pucci method, as addressed by Cabr\'e in \cite{C15} (see also \cite{CDDM}). Unfortunately, it seems to be an intrinsic drawback of this  approach to provide  a non-optimal estimate. Actually it is natural to conjecture that, as in case of the well-known Faber-Krahn inequality for the $p$-Laplacian, the product  $\overline \lambda _p (\Omega) |\Omega| ^ { 2/n}$ should be minimal on balls.
 In other words, the optimal lower bound expected for the product $\overline \lambda _p (\Omega) |\Omega| ^ {  2/n}$ is the constant
 $K^*_{n,p}:= \overline \lambda _p (B)  |B| ^ {2/n} $.
Notice that due to the scaling invariance $B$ can be an arbitrary ball here.
To prove such an optimal bound seems to be a very interesting and delicate problem. The symmetrization technique usually employed to prove the Faber-Krahn inequality for the $p$-Laplacian does not work here because the normalized $p$-Laplacian operator does not have a variational nature.

 To demonstrate that \eqref{f:Knp} is not optimal for balls let us sketch a quick comparison between the values of $K _{n,p}$ and $K ^* _{n,p}$.
  Clearly, by Theorem \ref{t:FK}, the quotient $K^*_{n,p} / K _{n, p}$ is larger than or equal to $1$.
 In order to evaluate the presumed accuracy of our estimate, one can evaluate how far it is from $1$.
 As shown in Lemma \ref{l:ball} below, we have
\begin{equation}\label{f:Knp2}
K^*_{n,p} =  \frac{\pi(p-1)}{p} \Gamma \Big ( 1 + \frac{n}{2} \Big ) ^ { -2/n} \big ( \mu _1 ^ { (- \alpha ) }  \big ) ^ 2 \,,
\end{equation}
where $\mu_1^{(-\alpha)}$ denotes the first zero of the Bessel function
$J _{-\alpha}$, with $\alpha = \tfrac{p-n}{2(p-1)}$.
The plots in Figure \ref{fig:comparison} left and right, obtained with Mathematica,  represent this ratio in two and three dimensions as a function of $p$.
Observe that both maps
 $$p \mapsto  g_2 (p):= \frac{K^*_{2,p}}{K _{2, p}} \, , \qquad p \mapsto g_3 (p):= \frac{K^*_{3,p}}{K _{3, p}} \,$$
 turn out to be minimal at $p = 2$, with
 $$g_2 (2) \approx 1.446\, , \qquad g_3 (2)\approx 1.561\,.$$
This shows that the constant $K_{n,p}$  in Theorem \ref{t:FK} is not optimal, not even in the linear case $p=2$.
\begin{figure}[ht]
\begin{minipage}{0.5\linewidth}
\centering
\includegraphics[height=4cm]{fig1.pdf}
\end{minipage}%
\begin{minipage}{0.5\linewidth}
\centering
\includegraphics[height=4cm]{fig2.pdf}
\end{minipage}
\caption{Plots of $g_2(p)$ and $g_3(p)$}
\label{fig:comparison}
\end{figure}

The proofs of Theorems \ref{t:simple} and \ref{t:FK} are given in Section \ref{sec:proofs} below,
after recalling the definition of viscosity solution to problem \eqref{f:pb1} and providing some preliminary results.

\section{Proofs}\label{sec:proofs}

In the notation of viscosity theory, the equation $- \Delta^N_p u = \lambda _p u $  can be rewritten as
\begin{equation}\label{F_p}
F_p ^N(\nabla u, \nabla ^ 2 u ) = \lambda_p u\, ,
\end{equation}
 where $F_p^N$  is defined on $(\mathbb{R}^n\setminus \{ 0 \})\times S(n)$ and $S(n)$ denotes the space of $n\times n$ symmetric matrices, with
\begin{equation}\label{f:Fp}
F_p^N(\xi, X) := - \frac{p-2}{p} |\xi|^{-2} \pscal{X\xi}{\xi} - \frac{1}{p}\, \Trace(X)\quad \forall \xi \in \mathbb{R} ^n \setminus \{ 0 \}\,,\ X \in S(n).
\end{equation}
At $\xi=0$ the function $F_p^N$ is discontinuous.  In this case, following \cite{CIL} we request from a viscosity solution of \eqref{F_p}
that it is a viscosity subsolution of $(F^{N}_{p})_*(Du,D^2u)=\lambda_pu$ and a viscosity supersolution of $(F^{N}_p)^*(Du,D^2u)=\lambda_pu$.
 Here $(F^{N}_p)^*$ is the upper semicontinuous hull and $(F^N_{p})_*$ is the lower semicontinuous hull of $F^N_p$.

Now since $F_p^N$ is given by
\begin{equation*}
F_p^N(\xi,X)=
-\frac{1}{p}\left(\delta_{ij}+(p-2)\frac{\xi_i\xi_j}{|\xi|^2}\right)X_{ij}
 \hbox{ for }\xi\not=0
\end{equation*}
we have to compute its semicontinuous limits as $\xi\to 0$. Each
symmetric matrix $X$ has real eigenvalues, and we order them according to magnitude as $\lambda_1(X)\leq\lambda_2(X)\leq\cdots\leq\lambda_n(X)$.
Then a simple calculation shows that
\begin{equation}\label{F_p^N_*}
{(F_{p}^N)_*}(0,X)=\begin{cases}
\,-\, \frac{1}{p}\sum_{i=1}^{n-1}\lambda_i-\frac{p-1}{p}\lambda_n & \hbox{ if } p\in[2,\infty],\\
-\frac{1}{p}\sum_{i=2}^{n}\lambda_i-\frac{p-1}{p}\lambda_1 &\hbox{ if } p\in[1,2],\\
\end{cases}
\end{equation}
and
\begin{equation}\label{F_p^N^*}
(F_p^N)^*(0,X)=\begin{cases}
-\frac{1}{p}\sum_{i=2}^{n}\lambda_i\,-\frac{p-1}{p}\lambda_1 &\hbox{ if } p\in[2,\infty],\\
-\frac{1}{p}\sum_{i=1}^{n-1}\lambda_i-\frac{p-1}{p}\lambda_n & \hbox{ if } p\in[1,2].
\end{cases}
\end{equation}
In \cite{Bru} these bounds for the normalized $p$-Laplacian are called dominative and submissive $p$-Laplacians and studied in more detail. Anyway, the above considerations serve as a motivation for the following

\begin{definition}\label{def:vs}
Given a symmetric matrix $A\in S(n)$, we denote by $M (A)$ and $m (A)$ its greatest and smallest eigenvalue.

\smallskip
-- An upper semicontinuous function $u :\Omega \to \mathbb{R}$ is  a viscosity subsolution of $- \Delta^N_p u = \lambda_p u$
in $\Omega$ if, for every point $x$ in $\Omega$ and every smooth function $\varphi$  which touches $u$ from above at $x$ (and for which $u - \varphi$ attains a local maximum at $x$) it holds
$$
\begin{cases}
- \Delta^N_p \varphi (x) \leq \lambda_p \varphi  (x)  & \text{ if } \nabla \varphi (x) \neq 0
\\  \medskip
- \frac{1}{p} \Delta \varphi (x) - \frac{(p-2) }{p} M (D ^ 2 \varphi (x) ) \leq \lambda_p \varphi (x) & \text { if  } \nabla \varphi (x) = 0 \text{ and } p \geq 2
\\  \medskip
- \frac{1}{p} \Delta \varphi (x) - \frac{(p-2) }{p} m (D ^ 2 \varphi (x) ) \leq \lambda_p \varphi (x) & \text { if  } \nabla \varphi (x) = 0 \text{ and } p \leq 2 .
\end{cases}
$$
-- A lower semicontinuous function $u :\Omega \to \mathbb{R}$ is  a viscosity supersolution of $- \Delta^N_p u = \lambda_p u$
in $\Omega$ if, for every point $x$ in $\Omega$ and every smooth function $\varphi$  which touches $u$ from below at $x$ (and for which $u - \varphi$ attains a local minimum at $x$) it holds
$$
\begin{cases}
- \Delta^N_p \varphi (x) \geq \lambda_p \varphi  (x)  & \text{ if } \nabla \varphi (x) \neq 0
\\  \medskip
- \frac{1}{p} \Delta \varphi (x) - \frac{(p-2) }{p} m (D ^ 2 \varphi (x) ) \geq \lambda_p \varphi (x) & \text { if  } \nabla \varphi (x) = 0 \text{ and } p \geq 2
\\  \medskip
- \frac{1}{p} \Delta \varphi (x) - \frac{(p-2) }{p} M (D ^ 2 \varphi (x) ) \geq \lambda_p \varphi (x) & \text { if  } \nabla \varphi (x) = 0 \text{ and } p \leq 2 .
\end{cases}
$$
-- A continuous function $u :\Omega \to \mathbb{R}$ is a  viscosity supersolution to $- \Delta^N_p u = \lambda_p u$  if it is both a viscosity supersolution and a viscosity subsolution.
\end{definition}

\begin{remark}\label{r:propF}
For later use we mention that the function $F_p^N$  satisfies the following identities:
\begin{itemize}
\item[(i)] $-F_p^N (t \xi, \mu X) =  - \mu F_p^N ( \xi, X)\quad \forall t \in \mathbb{R} \setminus \{ 0 \}$,  $\xi \in \mathbb{R} ^n\setminus \{ 0 \} $, $\mu \in \mathbb{R}$, and $X\in S(n)$.
\item[(ii)]  $-F_p^N (\xi, X ) \leq 0$ for any $\xi \in \mathbb{R} ^n \setminus \{ 0 \}$ and $X\in S(n)$ with $X\leq 0$.
This follows from  \eqref{f:Fp}, since the eigenvalues $\lambda_i(X)$ are assumed nonpositive.

\item[(iii)]  As a consequence of \eqref{f:Fp}, \eqref{F_p^N_*} and \eqref{F_p^N^*},
for every $\xi\in\mathbb{R}^n$ and $X\in S(n)$ we have that
\begin{gather*}
\frac{(p-1)\wedge p}{p}\,
\Trace(X)
\leq
-F_p^N(\xi,X)
\leq
\frac{(p-1)\vee p}{p}\,
\Trace(X)\,,
\qquad
\text{if}\  X\geq 0,
\\
\frac{(p-1)\vee p}{p}\,
\Trace(X)
\leq
-F_p^N(\xi,X)
\leq
\frac{(p-1)\wedge p}{p}\,
\Trace(X)\,,
\qquad
\text{if}\ X\leq 0.
\end{gather*}
\end{itemize}
\end{remark}

\bigskip
 For $x \in \mathbb{R} ^n$ and $R>0$, we denote by $B _R (x)$ the open ball of radius $R$ centred at $x$. We also set for brevity $B _ R := B _ R (0)$.

\begin{lemma}[First eigenvalue of the ball]\label{l:ball}
For any $p \in (1, + \infty)$, we have
$$\overline \lambda _p ( B _R) = \frac{p-1}{p} \Big ( \frac{\mu_1^{(-\alpha)}} {R} \Big ) ^ 2   = K ^ *_ {n,p} |B _R| ^ {-2/n}\, ,  $$
where $\mu_1^{(-\alpha)}$ denotes the first zero of the Bessel function
$J _{-\alpha}$, for $\alpha = \tfrac{p-n}{2(p-1)}$  (and the constant $K ^ * _{n,p} $ is defined in \eqref{f:Knp2}).
 \end{lemma}

\begin{proof}
Set $\lambda_p (R):= \frac{p-1}{p} \Big ( \frac{\mu_1^{(-\alpha)}} {R} \Big ) ^ 2$. We first prove that $\overline \lambda _p ( B _R) \geq \lambda_p(R)$.  By definition, this amounts to show that problem \eqref{f:pb1} admits a positive viscosity subsolution when $\lambda_p = \lambda_p (R)$. We search for a radial solution and make the ansatz $u (x) = g (|x|)$. In terms of the function $g = g (r)$, problem \eqref{f:pb1} can be written as (see \cite{KKK})
\begin{equation}\label{f:pbradial}
\begin{cases}
\displaystyle -  g''(r) - \Big ( \frac{n-1}{p-1} \Big )  \frac{g'(r)}{r} =\Big (  \frac{p}{p-1}  \Big ) \lambda _p \, g (r) & \text{ on } (0, R)
\\ \noalign{\medskip}
g (R) = 0  &
\\  \noalign{\medskip}
g' (0) = 0 \,. &
\end{cases}
\end{equation}
For $p=2$ the left hand side in the differential equation is just the classical Laplacian, evaluated in polar coordinates for $g(|x|)$. For other $p$ it can be interpreted as a linear Laplacian in a fractional dimension. This was done in \cite{KKK}, and a full spectrum and orthonormal system of radial eigenfunctions was derived. The first eigenfunction is a (positive) multiple of $r^\alpha J_{-\alpha}(\mu_1^{(-\alpha )}\frac{r}{R})$. This function is positive in $B_R$.

Finally, let us show that the equality $\overline \lambda _p ( B _R)=  \lambda _p(R)$ holds. For this we use an idea from \cite{MPR2}, there given for $p>n$. Assume by contradiction that
$\overline \lambda _p ( B _R)>  \lambda _p(R)$. Choose $\rho \in (0, R)$ such that
$\overline \lambda _p ( B _R)> \lambda_p (\rho) > \lambda_p(R)$,  and
let $g _\rho$ be a positive solution to problem
\begin{equation}\label{f:pbradial2}
\begin{cases}
\displaystyle -  g''(r) - \Big ( \frac{n-1}{p-1} \Big )  \frac{g'(r)}{r} =\Big (  \frac{p}{p-1}  \Big ) \kappa _p (\rho) \, g (r) & \text{ on } (0, \rho)
\\ \noalign{\medskip}
g (R) = 0  &
\\  \noalign{\medskip}
g' (0) = 0 \,. &
\end{cases}
\end{equation}
Then
the function $w$ defined on $B _R$ by
$w (x) = g _\rho (|x| )$ if $|x| \leq \rho$ and $0$ otherwise turns out to satisfy
$- \Delta^N_p w \leq \lambda _p (\rho) w$ in $B _R$  and $w \leq 0$ on $\partial B _R$.
In view of Remark \ref{r:propF} (i) and (ii), the  operator $\Delta^N_p$
 satisfies the assumptions of the
comparison result stated in \cite[Theorem 1.1]{BiDe2006}.
We infer that $w\leq 0$ in $B_R$, a contradiction.
\end{proof}

 \bigskip
 \begin{lemma}[Positivity of the eigenvalue] \label{l:positivity} For every open bounded domain  $\Omega \subset \mathbb{R} ^n$, we have $\overline \lambda _p (\Omega) >0$.
 \end{lemma}

\begin{proof}
From its definition, it readily follows that $\overline \lambda _p$ is monotone decreasing under domain inclusion, {\it i.e.} $\overline \lambda _p (\Omega_1 ) \geq \overline \lambda _p (\Omega_2 )$ if $\Omega _1 \subseteq \Omega _2$. In particular, for every open bounded domain $\Omega$, we have
$\overline \lambda _p (\Omega) \geq \overline \lambda _ p  ( B _ R)$, where $R = R (\Omega) = \inf \big\{ r >0 \, :\, \Omega \subset B _ r (x)\text{ for some } x \big \}$.   Invoking Lemma \ref{l:ball}, we obtain the positivity of $\overline \lambda _p (\Omega)$.
\end{proof}

\bigskip

In the following Lemma we do not assume differentiability of $u$ on the boundary. Nevertheless we can bound the difference quotient in interior normal direction from below.

\begin{lemma}[Hopf type Lemma]\label{t:hopf}
Assume that $\Omega\subset\mathbb{R}^n$ satisfies a uniform interior sphere
condition,
and let $u\in C(\overline{\Omega})$ be a positive viscosity supersolution
of $-\Delta^N_p u = 0$  in $\Omega$ such that $u=0$ on $\partial\Omega$.
Then there exists a constant $\kappa > 0$ such that for any $y\in \partial\Omega$
\begin{equation}\label{f:hopf}
\liminf_{t\to 0+} \frac{u(y-t\nu(y))}{t} \geq \kappa.
\end{equation}
Here $\nu$ denotes the unit outer normal to $\partial \Omega$,
\end{lemma}

\begin{proof}
This follows from realizing that the normalized $p$-Laplacian satisfies the assumptions in \cite[Theorem 1]{BDL}.
\end{proof}

\bigskip

{\bf Proof of Theorem \ref{t:simple}}.
Let $u$ and $v$ be two positive eigenfunctions associated with $\overline \lambda _p (\Omega)$.
Inspired by the appendix in \cite{Sak} we set
\[
\begin{array}{ll}
& a:= \sup \Big \{ t \in \mathbb{R} \ : \ u - tv >0 \text{ in } \Omega \Big \}
\\
\noalign{\medskip}
& b:= \sup \Big \{ t \in \mathbb{R} \ : \ v - tu >0 \text{ in } \Omega \Big \}\,.
\end{array}
\]
 Clearly, we have
 \begin{equation}\label{f:pos}
 u - a v \geq 0 \quad  \text{ and } \quad v-bu \geq 0 \qquad \text{ in } \Omega \,.
 \end{equation}
We claim that $a$ and $b$ are strictly positive.
Indeed, the functions $u$ and $v$ are of class $C ^ {1, \alpha}$ up to the boundary (see \cite[Proposition 3.5]{BirDem2010} or  \cite[Theorem 1.1]{APR17}). Then, applying Lemma \ref{t:hopf}
to $u$ and $v$, we see that
\begin{equation}\label{f:hopf2}
\frac{\partial u }{\partial \nu} <0 \quad \text{ and } \quad \frac{\partial v }{\partial \nu} <0 \qquad \text{ on } \partial \Omega\,.
\end{equation}
Hence, for $t \in \mathbb{R} _+$ small enough,  $\frac{\partial}{\partial \nu} (u - t v)$ is strictly negative on $\partial \Omega$, so that  there exists $\overline t>0$   and a neighbourhood $\mathcal U$ of $\partial \Omega$ such that
$u- \overline tv >0$ in $\mathcal U$.  It follows that
$$u- m  v >0 \text{  in } \Omega \quad \text{ for }
m<\min \Big \{ \overline t, \frac{\min _{\Omega \setminus \overline { \mathcal U } } u } {\max_{\Omega \setminus \overline {\mathcal U } } v} \Big \}\,.$$
Thus $a \geq m >0$. Arguing in the same way with $u$ and $v$ interchanged we obtain $b>0$, and our claim is proved.

Now, to obtain the result, we are going to show that there exists a neighbourhood $\mathcal V$ of $\partial \Omega$ such that
\begin{equation}\label{f:neig}
u - a v = 0 \quad \text{ and } \quad v - b u = 0 \quad \text{ in } \mathcal V\,.
\end{equation}
This implies $u - (ab) u = 0$ in $\mathcal V$ and, in view of the condition $u >0 $ in $\Omega$, $b = a ^ {-1}$.
 The latter equality, combined with \eqref{f:pos}, implies $u - av = 0$  in $\Omega$ as required.

Let us show how to obtain the first equality in \eqref{f:neig}, the derivation of the second one is  completely analogous.

By the regularity of $\partial \Omega$, its unit outer normal $\nu$ can be extended to a smooth unit vector field, still denoted by $\nu$, defined in an open connected neighbourhood of $\partial \Omega$.
Then, by \eqref{f:hopf2} and the $C ^ 1$ regularity of $u$ and $v$ on $\overline \Omega$, we infer that there exist $\delta >0$ and an open connected neighbourhood $\mathcal V$ of $\partial \Omega$ such that
\begin{equation}\label{f:nondeg}
\frac{\partial u }{\partial \nu} <-\delta \quad \text{ and } \quad \frac{\partial v }{\partial \nu} <- \delta \qquad \text{ in } \overline {\mathcal V}\,.
\end{equation}
This implies first of all that the PDE solved by $u$ and $v$ is nondegenerate in $\mathcal V$, which in turn, by standard elliptic regularity (see \cite{GT}) yields that $u$ and $v$ are of class $C ^ \infty$ in $\mathcal V$.  Moreover, from the inequality
$$0 \leq \overline \lambda _p u - \overline \lambda _p (av) \qquad \text{ in } \Omega$$
we infer that
$$0 \leq - \Delta^N_p u -  \big ( - \Delta^N_p (av) \big )  = L_p (u - av)  \qquad \text{ in } \Omega\, , $$
 where $L_p w = \sum _{i,j= 1}^n c_{ij} w _{x_i x_j} + \sum _{i = 1} ^ n d _i w _{x_i}$ is the linear operator defined by
 $$\begin{array}{ll}
 & \displaystyle  c _{ij}:= \int _0 ^ 1 \frac{\partial F_p^N }{\partial X _{ij} } \big ( s  \nabla u+ ( 1-s ) \nabla v,  s  \nabla^2 u+ ( 1-s ) \nabla^2  v  \big ) \, ds
 \\
 \noalign{\bigskip}
 & \displaystyle  d _{i}:= \int _0 ^ 1 \frac{\partial F_p^N }{\partial \xi _{i}} \big ( s  \nabla u+ ( 1-s ) \nabla v,  s  \nabla^2 u+ ( 1-s ) \nabla^2  v  \big ) \, ds \,.
 \end{array}
 $$
 In particular, since
 $$ \frac{\partial F_p^N }{\partial X _{ij} } = - \frac{p-2}{p}\frac{1}{ |\xi| ^ 2} \xi _i \xi _j - \frac{1}{p}
  \delta _{ij}\,$$
  and, from \eqref{f:nondeg},
  $$\forall s \in [0,1], \quad s  \frac{\partial u}{\partial \nu} +  (1-s) \frac{\partial(a v) }{\partial \nu}  \leq  - \min \{ \delta, a \delta \} <0 \quad \text{ in } \mathcal V\, ,$$
  we see that $L_p$ is {\it uniformly elliptic} in
the connected set $\mathcal V$.
 Then, to achieve our proof, it is enough to show that there exists some point $x^* \in \mathcal V$ where the function $
 u-av$ vanishes. Indeed, if this is the case, we have:
 $$
 \begin{cases}
L_p (u-av) \geq 0   & \text{ in } \mathcal V
\\
\noalign{\medskip}
u - av \geq 0 & \text{ in } \mathcal V
\\
\noalign{\medskip}
(u-av) (x^*) = 0\,.
\end{cases}
 $$
 By the strong maximum principle for uniformly elliptic operators  \cite[Theorem 3.5]{GT}, it will follow that
 $u-av \equiv 0$ in $\mathcal V$ as required.
 We point out that, without the connectedness of $\partial \Omega$ (and hence of $\mathcal V$), the two equalities in \eqref{f:neigh} might be obtained in two, a priori distinct, connected components of $\mathcal V$, and this would not be sufficient to infer that $u$ and $v$ are proportional.

 To conclude, let us now show that $u-av$ vanishes at some point $x^*$ in $\mathcal V$. As an intermediate step we notice that the function $u-av$ must vanish at some point $\overline x$ in $\Omega$. Otherwise, we would have:
 $$
 \begin{cases}
L_p (u-av) \geq 0   & \text{ in } \mathcal V
\\
\noalign{\medskip}
u - av > 0 & \text{ in } \mathcal V
\\
\noalign{\medskip}
u-av \equiv  0 & \text{ on } \partial \Omega\,.
\end{cases}
 $$
By applying Hopf's boundary point lemma for uniformly elliptic operators \cite[Lemma 3.4]{GT}, we infer that $\frac{\partial}{\partial \nu} (u - av) < 0$ on $\partial \Omega$. By continuity, this inequality, combined with the strict one $u-av >0 $ in $\Omega$ that we are assuming by contradiction, implies that $u - (a + \eta) v >0$  in $\Omega$ for some $\eta>0$. But this contradicts the definition of $a$.

Now, we choose an open bounded set $\omega$ with smooth boundary such that
\[
\overline \omega \subset\Omega\, , \quad \overline x \in \omega\, , \quad \partial \omega \subset \mathcal V\,.
\]
We assert that there is a point $x ^* \in \partial \omega $ where $u-av$ vanishes (and this point does the job since
$\partial \omega \subset \mathcal V$).  Assume the contrary. Then by continuity we have
$u - av \geq \varepsilon >0$ on $\partial \omega$ for some
$\varepsilon >0$.  Then the two functions $u$ and $w:= av +\varepsilon$  satisfy
\[
\begin{cases}
- \Delta^N_p u  = \overline \lambda _p u  \geq \overline \lambda _p (av)= - \Delta^N_p w & \text{ in } \omega
\\
u \geq w & \text{ on } \partial \omega\,.
\end{cases}
\]
 In view of Lemma \ref{l:positivity}, the continuous function $f:=\overline \lambda _p u$ is strictly positive in $\omega$.  Now we can apply the comparison principle proved in \cite[Thm.~2.4]{LuWang2008},
and we infer that
\[
u \geq w \quad \text{ in } \omega\,.
\]
In particular, since $\omega$ contains the point $\overline x$, we have
$$u (\overline x ) \geq w (\overline x) = av (\overline x) +\varepsilon\,,$$
which gives a contradiction since $u (\overline x ) = av (\overline x)$.
\qed

\bigskip

 \bigskip
 In order to prove Theorem \ref{t:FK}, we need some preliminary results.

 Let $u$ be  a positive eigenfunction associated with $\overline \lambda _p (\Omega)$.
The approximations of
$u$ via supremal convolution are defined for $\e >0$ by
\begin{equation}\label{f:ue}
u ^ \e (x) := \sup _{y \in \Omega} \Big \{ u (y)  - \frac { |x-y| ^ 2 }{2 \e}  \Big \} \qquad \forall x \in \Omega\,.
\end{equation}

Let us start with a preliminary lemma in which we recall
some basic well-known properties of the functions $u ^ \e$.
To fix our setting let us define
\[
\rho(\e) := 2 \sqrt{\e \, \|u\|_\infty},
\qquad
 \Omega^{\rho(\e)} := \{x\in\Omega:\ d_{\partial\Omega}(x) > \rho(\e)\}\,,
\]
then for every \(x \in \Omega^{\rho(\e)}\) the supremum in \eqref{f:ue}
is attained at a point \(y_\e(x) \in \overline{B}_{\rho(\e)}(x)\subset\Omega\).
Thus, setting
\begin{equation}\label{f:Ue}
U _\e:= \big \{ x \in \Omega \ :\ u (x) > \e \big \} \, , \qquad
A_\e := \big \{ x \in U _\e \ :\ d_{\partial U _\e}(x) > \rho(\e) \big \}\, ,
\end{equation}
so that by definition
\begin{equation}\label{d:ue2}
u ^ \e (x) =
u(y_\e(x)) - \frac{|x-y_\e(x)|^2}{2\e} =
\sup _{y \in U _\e} \Big \{ u (y)  - \frac { |x-y| ^ 2 }{2 \e}  \Big \}
\qquad \forall x \in \overline A_\e\,.
\end{equation}
In what follows, we shall always assume that $\e \in (0, 1)$ is small enough
to have $A_\e \neq \emptyset$.
Moreover, let us define
\begin{equation}\label{f:omegae}
m_\e := \max_{\partial A_\e} u^\e,
\qquad
\Omega_\e := \{x\in A_\e : \ u^\e (x) > m_\e \}\,.
\end{equation}

\begin{lemma}\label{l:approx1}
Let
$u$ be a positive eigenfunction associated with $\overline \lambda _p (\Omega)$, let
   $u ^\e$ be its supremal convolutions according to \eqref{f:ue}, and let  $\Omega_\e$ be the domains defined in \eqref{f:omegae}. Then:
\begin{itemize}
\item[(i)] $u ^ \e$ is semiconvex in $\Omega _\e$;

\smallskip
\item[(ii)] $u ^\e$ is a viscosity sub-solution to $-\Delta^N_p  u -\overline \lambda _p (\Omega) u = 0$ in $\Omega_\e$;

\smallskip
\item[(iii)] as $\e \to 0 ^+$,  $u ^ \e$ converge to $u$ uniformly in $\overline \Omega$.
Hence $m_\e \to 0$ and $\Omega_\e$ converges to $\Omega$
in Hausdorff distance;

\smallskip
\item[(iv)]
as $\e \to 0 ^+$,  $\nabla u ^ \e \to \nabla u$ locally uniformly in $\Omega$.
\end{itemize}
\end{lemma}

\begin{proof}
(i)  We have $u ^ \e = - ( -u ) _\e$, where $( -u ) _\e$ is the infimal convolution defined by
\[
( -u ) _\e (x) := \inf _{y \in U _\e} \Big \{-  u (y)  + \frac { |x-y| ^ 2 }{2 \e}  \Big \} \qquad \forall x \in \Omega _\e\,.
\]
{}From \cite[Proposition 2.1.5]{CaSi}, it readily follows that $(-u) _\e$ is semiconcave on $\Omega _\e$, and hence that
$u ^ \e$ is semiconvex on $\Omega _\e$.

 (ii) The notion of of viscosity subsolution according to Definition \ref{def:vs} can be reformulated by asking that, for every $x \in \Omega$ and every $(\xi, X)$ in the second order
superjet $J ^ {2, +} _\Omega u (x)$   (classically defined as in \cite{CIL}), it holds
$$
\begin{cases}
F_p^N(\xi, X) \leq \lambda_p u  (x)  & \text{ if } \xi \neq 0
\\  \medskip
- \frac{1}{p} \Trace(X)  - \frac{(p-2) }{p} M (X) \leq \lambda_p u (x) & \text { if  } \xi = 0 \text{ and } p \geq 2
\\  \medskip
- \frac{1}{p} \Trace(X)  - \frac{(p-2) }{p} m (X) \leq \lambda_p u (x) & \text { if  }\xi = 0 \text{ and } p \leq 2 .
\end{cases}
$$
Then, in order to prove (ii), it is enough to show that, for every fixed point $x \in \Omega _\e$, any pair  $(p, X) \in  J  ^ {2, +} _{\Omega _\e} u ^ \e (x)$ belongs to $J ^ {2, +} _{\Omega} u (y)$ for some other point $y \in \Omega _\e$. In fact, the so-called magic properties of supremal convolution ({\it cf.} \cite[Lemma A.5]{CIL}) assert precisely that any
$(p, X) \in  J  ^ {2, +} _{\Omega _\e} u ^ \e (x)$ belongs to
$J ^ {2, +} _{\Omega _\e} u (y)$, where $y$ is a point at which the supremum which defines
$u ^ \e (x)$ is attained.  Since $y \in U _\e \subset \Omega _\e$, it holds $J ^ {2, +} _{\Omega} u (y)= J  ^ {2, +} _{\Omega _\e} u ^ \e(x)$.

(iii) For these convergence properties we refer to \cite[Thm.\ 3.5.8]{CaSi},
\cite[Lemma 4]{CFd}.

(iv)
Since $u\in C^1(\Omega)$, this property follows
from \cite[Lemma 10]{CF7}.
\end{proof}

\bigskip\bigskip

\begin{lemma}\label{l:approx2}
Let
$u$ be a positive eigenfunction associated with $\overline \lambda _p (\Omega)$,  let
   $u ^\e$ be its supremal convolutions according to \eqref{f:ue}, and let $\Omega_\e$ be the domains defined   in \eqref{f:omegae}.   Let $v ^ \e$ be the continuous functions defined by
\begin{equation}\label{f:ve}
v ^ \e (x) := \begin{cases}
\log ( u ^ \e ) & \text{ if } x \in \Omega _\e
\\
\log (m _ \e) & \text{ if } x \in \mathbb{R} ^n \setminus \Omega _\e
 \end{cases}
 \end{equation}
and, for $\sigma >0$, let $\Gamma _\sigma (v ^ \e)$ be the concave envelope of $v ^ \e$ on the set
   \begin{equation}\label{f:neigh}
   (\Omega _\e ^* ) _\sigma:= \Big \{ x \in \mathbb{R} ^n \ :\ {\rm dist} (x, \Omega _\e ^ *) \leq \sigma \Big \}\,,\end{equation}
 $\Omega_\e^*$ being the convex envelope of $\Omega _\e$.
  Then:
\begin{itemize}
\item[(i)]  $\Gamma_\sigma ( v ^ \e)$ is locally $C ^ {1, 1}$ in $(\Omega _\e ^* ) _\sigma$;
\smallskip
\item[(ii)]  at any $x \in (\Omega _\e ^*)_\sigma$ such that $\det   D ^ 2 (\Gamma _\sigma  ( v ^ \e)(x) \neq 0$,  it holds
$v^\e(x) = \Gamma_\sigma ( v ^ \e)(x)$;
\smallskip
\item[(iii)]  $v ^ \e$ is a viscosity sub-solution to $- \Delta^N_p v = \overline \lambda _p (\Omega ) + \frac{p-1}{p} |\nabla v  | ^ 2$ in $\Omega _\e$.
 \end{itemize}
\end{lemma}

\begin{proof}
We observe that by \cite[Prop.2.1.12]{CaSi} and Lemma \ref{l:approx1}(i) also $v^\e$ is semiconvex.
Statements (i) and (ii) follow now from \cite[Lemma 5]{CDDM} since, for every fixed $\e>0$, the function $v ^ \e - \log (m _ \e) $ satisfies the assumptions of such result on the convex domain $\Omega_\e^*$.

Statement (iii) follows from part (iii) in Lemma \ref{l:approx1} above, combined with the fact that, if a smooth function $\varphi$
touches $v ^ \e$ from above at $x$, the smooth function $e ^ \varphi$ touches $u ^ \e$ from above at $x$.
\end{proof}

\bigskip
{\bf Proof of Theorem \ref{t:FK}}.
Throughout the proof we write for brevity $ \lambda _p$ in place of  $\overline \lambda _p (\Omega)$.
 Set
\begin{equation}\label{f:defg}
g (s): = \frac{1} {\Big ( \lambda _p + \frac{p-1}{p} s ^ 2  \Big ) ^n }\, , \qquad s \geq 0  \,,
\end{equation}
 and
 $$I _g := \int _{\mathbb{R} ^n} g (|\xi | ) \, d \xi \,.$$
 By direct computation in polar coordinates, the value of $I _g$ is given by
 \begin{equation}\label{f:Ig}
\begin{split}
I _ g &  = \frac{\omega _n}{ \lambda _p ^n } \int _0 ^ { + \infty}
  \frac{\rho ^ {n-1}} {\Big (1 + \frac{p-1}{p  \lambda _p} \rho ^ 2  \Big ) ^n } \, d \rho   =
  \frac{\omega _n}{ \lambda _p ^n } \Big (  \frac{p  \lambda _p} {p-1} \Big ) ^ {n/2} \int _0 ^ { + \infty}
  \frac{t ^ {n-1}} {\big (1 + t^ 2  \big ) ^n } \, d t
\\
  &  = 2^{1-n} \, \pi^{(n+1)/2}   \Big (  \frac{p } {p-1} \Big ) ^ {n/2}
\lambda _p ^{-n/2}\Gamma \Big ( \frac{n+1}{2} \Big )^{-1}\,,
\end{split}
\end{equation}
where $\omega _n := \mathcal H ^ {n-1} (S ^ {n-1}) = 2 \pi^{n/2} \Gamma(n/2)^{-1}$.

 On the other hand, a natural idea in order to estimate $I _g$ (and hence $\lambda _p$) in terms of the measure of $\Omega$, is to apply the change of variables formula to the map $\xi = -\nabla v (x)$, with $v (x) = \log u (x)$ and
$u$ being a positive eigenfunction associated with $\lambda _p $.

This is suggested by the fact that, as one can easily check, $v$ is a viscosity solution to
\begin{equation}\label{f:eqv}\begin{cases}
-\Delta^N_p v = \lambda _p + \frac{p-1}{p} |\nabla v| ^ 2 & \text{ in } \Omega
\\  \noalign{\medskip}
v = - \infty & \text{ on } \partial \Omega\,,
\end{cases}
 \end{equation}
combined with the observation that $-\nabla v$ maps $\Omega$ onto $\mathbb{R} ^n$, namely
 \begin{equation}\label{f:image}
 -\nabla v (\Omega) = \mathbb{R} ^n \,.
 \end{equation}
Indeed, for every $p \in \mathbb{R} ^n$,   the minimum over $\overline \Omega$ of the function $ -v (y) - p \cdot y$  is necessarily attained a point $x$ lying in the interior  of $\Omega$ (since $v = - \infty$ on $\partial \Omega$), and at such point $x$ we have  $p =- \nabla v (x)$.

 In view of \eqref{f:image}, we have
 $$I _ g  = \int _{  -\nabla v (\Omega ) } g (|\xi | ) \, d \xi   \,,$$

but unfortunately the map $\xi = -\nabla v (x)$
  is a priori not regular enough to apply directly the area formula. Therefore, we need to proceed by approximation.

Let $u ^\e$ be the supremal convolutions of $u$ according to \eqref{f:ue}, and let  $\Omega_\e$ be the domains defined in \eqref{f:omegae}. Then consider the functions $v ^ \e$ and the sets $(\Omega _\e  ^*)_\sigma$ defined  as in \eqref{f:ve} and \eqref{f:neigh}, and let $\Gamma _\sigma (v ^ \e)$ be the concave envelope of $v ^ \e$ on    $(\Omega _\e ^* ) _\sigma$.

By Lemma \ref{l:approx2} (i), we are in a position to
apply the area formula on $(\Omega_{\e} ^* ) _\sigma $ (see \cite[Section 3.1.5]{GMS1}) to the map $\xi =-\nabla \Gamma _\sigma( v^\e)$, and we obtain
\[
\begin{split}
\int _{ - \nabla \Gamma _\sigma( v^\e)  (  (\Omega _\e ^*)_\sigma)  } g (|\xi | ) \, d \xi
&
\leq \int _{  -\nabla \Gamma _\sigma( v^\e)  (  (\Omega _\e ^*)_\sigma) } g (|\xi | )\, {\rm card} ((-\nabla \Gamma _\sigma ( v ^\e) )^ {-1} (\xi) \cap (\Omega_{\e} ^* ) _\sigma )  \, d \xi
\\
&
= \int _{  (\Omega _\e ^*)_\sigma   } g (|\nabla \Gamma _\sigma(v ^\e) (x)| ) \, \det ( -  D ^ 2 \Gamma _\sigma( v^\e))  (x)  \, d x \,.
\end{split}
\]
 Now, we introduce  the {\it contact set}
 $$C_{\e, \sigma}:= \Big \{ x \in (\Omega _\e  ^*)_\sigma \ :\ v^\e(x) = \Gamma_\sigma ( v ^ \e)(x)\Big \} \,. $$
 Thanks to Lemma \ref{l:approx2} (ii), we have
  $$
 \int _{  (\Omega _\e ^*)_\sigma   } g (|\nabla \Gamma _\sigma(v ^\e) (x)| ) \, \det(- D ^ 2 \Gamma _\sigma( v^\e)  )(x) \, d x  =
 \int _{ C _{\e, \sigma}   } g (|\nabla v ^ \e(x)|)  \, \det (-D ^ 2  v ^ \e ) (x) \, d x  \,.   $$

Then we use the following pointwise estimates on $C_{\e, \sigma}$:
\begin{gather}
\det (-  D ^ 2 v^\e  ) \leq \Big (-\frac{1}{n} \Delta v ^ \e \Big ) ^ n   \label{f1}
\\
- \Delta v ^ \e  \leq - \frac {p}{(p-1) \wedge 1 } \Delta^N_p v ^ \e
\label{f2}
\\
- \Delta^N_p v ^\e \leq \lambda _p + \frac{p-1}{p} |\nabla v ^\e| ^ 2\,.
\label{f3}
\end{gather}
 Indeed, \eqref{f1} is consequence of the arithmetic-geometric inequality
 observing that by construction $- D ^ 2 v_\e$ is non-negative definite on $C_{\e, \sigma}$,
 \eqref{f2} holds by  Remark \ref{r:propF} (iii), and
 \eqref{f3} holds  thanks to Lemma \ref{l:approx2}  (i) and (iii), at every point of $C _{\e, \sigma}$ where $v ^ \e$ is twice differentiable (hence a.e.\ on $C _{\e, \sigma}$).

In this way we arrive at
 $$\begin{array}{ll} \displaystyle \int _{- \nabla \Gamma _\sigma( v^\e)  (  (\Omega _\e ^*)_\sigma) } g (|\xi | ) \, d \xi   & \displaystyle  \leq \int_{C _{\e, \sigma}  }  g (|\nabla v ^\e (x)|) \Big ( \frac {p}{n[(p-1) \wedge 1] } (- \Delta^N_p v  ^\e))  \Big )   ^n \, dx
 \\ \noalign{\medskip}
  & \displaystyle\leq     \Big ( \frac {p}{n[(p-1) \wedge 1] }  \Big )  ^n  |C _{\e, \sigma} | \, ,
  \end{array}
  $$
where in the last inequality we have exploited the choice of the function $g$ in \eqref{f:defg}.

So far, we have obtained the upper bound
$$ \int _{-\nabla \Gamma _\sigma( v^\e)  (  (\Omega _\e ^*)_\sigma)  } g (|\xi | ) \, d \xi  \leq  \Big ( \frac {p}{n[(p-1) \wedge 1] }  \Big )  ^n  |C _{\e, \sigma}|
\,.$$

Now we pass to the limit in the above inequality, first as $\sigma \to 0 ^+$, and then as $\e \to 0 ^+$.
In view of Lemma \ref{l:approx1} (iii) and \eqref{f:image}, we obtain
\[
\lim _{\e \to 0 ^+} \lim _{\sigma \to 0 ^+} \big ( -\nabla \Gamma _\sigma( v^\e)  (  (\Omega _\e ^*)_\sigma)  )  = \mathbb{R} ^n \qquad \text{ and } \qquad  \lim _{\e \to 0 ^+} \lim _{\sigma \to 0 ^+}  |C_{\e, \sigma}| \leq |\Omega | \,.
\]
We conclude that
\[
I _g \leq \Big ( \frac {p}{n[(p-1) \wedge 1] }  \Big )  ^n  |\Omega|
\,.
\]
The statement follows by inserting into the above inequality the explicit value of $I _g$ as given by \eqref{f:Ig}.
\qed

\bigskip \bigskip

{\bf Acknowledgments.}
G.C.\ and I.F.\
have been supported by the Gruppo Nazionale per l'Analisi Matematica,
la Probabilit\`a e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM).

\def$'${$'$}
\begin{bibdiv}
\begin{biblist}

\bib{APR17}{article}{
  author={Attouchi, A.},
  author={Parviainen, M.},
  author={Ruosteenoja, E.},
  title={$C^{1,\alpha }$ regularity for the normalized $p$-Poisson problem},
  language={English, with English and French summaries},
  journal={J. Math. Pures Appl. (9)},
  volume={108},
  date={2017},
  number={4},
  pages={553--591},
  issn={0021-7824},
  review={\MR {3698169}},
  doi={10.1016/j.matpur.2017.05.003},
}

\bib{BK18}{article}{
  author={Banerjee, {A.}},
  author={Kawohl, {B.}},
  title={Overdetermined problems for the normalized {$p$}-{L}aplacian},
  date={2018},
  journal={Proc. Amer. Math. Soc. Ser. B},
  volume={5},
  pages={18\ndash 24},
}

\bib{BDL}{article}{
  author={Bardi, M.},
  author={Da Lio, F.},
  title={On the strong maximum principle for fully nonlinear degenerate elliptic equations},
  journal={Arch. Math. (Basel)},
  volume={73},
  date={1999},
  number={4},
  pages={276--285},
  issn={0003-889X},
  review={\MR {1710100}},
  doi={10.1007/s000130050399},
}

\bib{BNV}{article}{
  author={Berestycki, {H.}},
  author={Nirenberg, {L.}},
  author={Varadhan, {S. R. S.}},
  title={The principal eigenvalue and maximum principle for second-order elliptic operators in general domains},
  date={1994},
  journal={Comm. Pure Appl. Math.},
  volume={47},
  number={1},
  pages={47\ndash 92},
}

\bib{BiDe2006}{article}{
  author={Birindelli, {I.}},
  author={Demengel, {F.}},
  title={First eigenvalue and maximum principle for fully nonlinear singular operators},
  date={2006},
  journal={Adv. Differential Equations},
  volume={11},
  number={1},
  pages={91\ndash 119},
}

\bib{BirDem2010}{article}{
  author={Birindelli, {I.}},
  author={Demengel, {F.}},
  title={Regularity and uniqueness of the first eigenfunction for singular fully nonlinear operators},
  date={2010},
  issn={0022-0396},
  journal={J. Differential Equations},
  volume={249},
  number={5},
  pages={1089\ndash 1110},
  url={https://doi.org/10.1016/j.jde.2010.03.015},
  review={\MR {2652165}},
}

\bib{blanc}{misc}{
  author={Blanc, {P.}},
  title={A lower bound for the principal eigenvalue of fully nonlinear elliptic operators},
  date={2017},
  note={preprint arXiv:1709.02455},
}

\bib{Bru}{article}{
  author={Brustad, {K.\ K.}},
  title={Superposition in the {$p$}-{L}aplace equation},
  date={2017},
  issn={0362-546X},
  journal={Nonlinear Anal.},
  volume={158},
  pages={23\ndash 31},
  url={https://doi.org/10.1016/j.na.2017.04.004},
  review={\MR {3661428}},
}

\bib{C15}{article}{
  author={Cabr\'{e}, {X.}},
  title={Isoperimetric, {S}obolev, and eigenvalue inequalities via the {A}lexandroff-{B}akelman-{P}ucci method: a survey},
  date={2017},
  journal={Chin. Ann. Math. Ser. B},
  volume={38},
  number={1},
  pages={201\ndash 214},
}

\bib{CaSi}{book}{
  author={Cannarsa, {P.}},
  author={Sinestrari, {C.}},
  title={Semiconcave functions, {H}amilton-{J}acobi equations and optimal control},
  series={Progress in Nonlinear Differential Equations and their Applications},
  publisher={Birkh\"auser},
  address={Boston},
  date={2004},
  volume={58},
}

\bib{CDDM}{article}{
  author={Charro, {F.}},
  author={De~Philippis, {G.}},
  author={Di~Castro, {A.}},
  author={M\'{a}ximo, {D.}},
  title={On the {A}leksandrov-{B}akelman-{P}ucci estimate for the infinity {L}aplacian},
  date={2013},
  journal={Calc. Var. Partial Differential Equations},
  volume={48},
  number={3-4},
  pages={667\ndash 693},
}

\bib{CIL}{article}{
  author={Crandall, {M.G.}},
  author={Ishii, {H.}},
  author={Lions, {P.L.}},
  title={User's guide to viscosity solutions of second order partial differential equations},
  date={1992},
  journal={Bull. Amer. Math. Soc. (N.S.)},
  volume={27},
  pages={1\ndash 67},
}

\bib{CFd}{article}{
  author={Crasta, {G.}},
  author={Fragal{\`a}, {I.}},
  title={On the {D}irichlet and {S}errin problems for the inhomogeneous infinity {L}aplacian in convex domains: regularity and geometric results},
  date={2015},
  issn={0003-9527},
  journal={Arch. Ration. Mech. Anal.},
  volume={218},
  number={3},
  pages={1577\ndash 1607},
  url={http://dx.doi.org/10.1007/s00205-015-0888-4},
  review={\MR {3401015}},
}

\bib{CFe}{article}{
  author={Crasta, {G.}},
  author={Fragal{\`a}, {I.}},
  title={A {$C^1$} regularity result for the inhomogeneous normalized infinity {L}aplacian},
  date={2016},
  issn={0002-9939},
  journal={Proc. Amer. Math. Soc.},
  volume={144},
  number={6},
  pages={2547\ndash 2558},
  url={http://dx.doi.org/10.1090/proc/12916},
  review={\MR {3477071}},
}

\bib{CFf}{article}{
  author={Crasta, {G.}},
  author={Fragal{\`a}, {I.}},
  title={Characterization of stadium-like domains via boundary value problems for the infinity {L}aplacian},
  date={2016},
  issn={0362-546X},
  journal={Nonlinear Anal.},
  volume={133},
  pages={228\ndash 249},
  url={http://dx.doi.org/10.1016/j.na.2015.12.007},
  review={\MR {3449756}},
}

\bib{CF7}{misc}{
  author={Crasta, {G.}},
  author={Fragal{\`a}, {I.}},
  title={Rigidity results for variational infinity ground states},
  date={2017},
  note={To appear in Indiana Univ. Math. J.},
}

\bib{Does}{article}{
  author={Does, K.},
  title={An evolution equation involving the normalized {$p$}-{L}aplacian},
  date={2011},
  issn={1534-0392},
  journal={Commun. Pure Appl. Anal.},
  volume={10},
  number={1},
  pages={361\ndash 396},
  url={https://doi.org/10.3934/cpaa.2011.10.361},
  review={\MR {2746543}},
}

\bib{EKNT}{article}{
  author={Esposito, L.},
  author={Kawohl, B.},
  author={Nitsch, C.},
  author={Trombetti, C.},
  title={The {N}eumann eigenvalue problem for the {$\infty $}-{L}aplacian},
  date={2015},
  issn={1120-6330},
  journal={Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl.},
  volume={26},
  number={2},
  pages={119\ndash 134},
  url={https://doi.org/10.4171/RLM/697},
  review={\MR {3341101}},
}

\bib{GMS1}{book}{
  author={Giaquinta, {M.}},
  author={Modica, {G.}},
  author={Sou{\v {c}}ek, {J.}},
  title={Cartesian currents in the calculus of variations. {I}},
  series={Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge. A Series of Modern Surveys in Mathematics [Results in Mathematics and Related Areas. 3rd Series. A Series of Modern Surveys in Mathematics]},
  publisher={Springer-Verlag},
  address={Berlin},
  date={1998},
  volume={37},
  isbn={3-540-64009-6},
  note={Cartesian currents},
  review={\MR {1645086 (2000b:49001a)}},
}

\bib{GT}{book}{
  author={Gilbarg, {D.}},
  author={Trudinger, {N.S.}},
  title={Elliptic partial differential equations of second order},
  publisher={Springer-Verlag},
  address={Berlin},
  date={1977},
}

\bib{Juut07}{article}{
  author={Juutinen, {P.}},
  title={Principal eigenvalue of a very badly degenerate operator and applications},
  date={2007},
  journal={J. Differential Equations},
  volume={236},
  number={2},
  pages={532\ndash 550},
}

\bib{JK}{article}{
  author={Juutinen, {P.}},
  author={Kawohl, {B.}},
  title={On the evolution governed by the infinity {L}aplacian},
  date={2006},
  issn={0025-5831},
  journal={Math. Ann.},
  volume={335},
  number={4},
  pages={819\ndash 851},
  url={https://doi.org/10.1007/s00208-006-0766-3},
  review={\MR {2232018}},
}

\bib{K0}{article}{
  author={Kawohl, B.},
  title={Variational versus PDE-based approaches in mathematical image processing},
  conference={ title={Singularities in PDE and the calculus of variations}, },
  book={ series={CRM Proc. Lecture Notes}, volume={44}, publisher={Amer. Math. Soc., Providence, RI}, },
  date={2008},
  pages={113--126},
  review={\MR {2528737}},
}

\bib{K11}{article}{
  author={Kawohl, B.},
  title={Variations on the $p$-Laplacian},
  conference={ title={Nonlinear elliptic partial differential equations}, },
  book={ series={Contemp. Math.}, volume={540}, publisher={Amer. Math. Soc., Providence, RI}, },
  date={2011},
  pages={35--46},
  review={\MR {2807407}},
  doi={10.1090/conm/540/10657},
}

\bib{KH}{article}{
  author={Kawohl, {B.}},
  author={Hor\'{a}k, {J.}},
  title={On the geometry of the {$p$}-{L}aplacian operator},
  date={2017},
  issn={1937-1632},
  journal={Discrete Contin. Dyn. Syst. Ser. S},
  volume={10},
  number={4},
  pages={799\ndash 813},
  url={https://doi.org/10.3934/dcdss.2017040},
  review={\MR {3640538}},
}

\bib{KKK}{article}{
  author={Kawohl, {B.}},
  author={Kr\"{o}mer, {S.}},
  author={Kurtz, {J.}},
  title={Radial eigenfunctions for the game-theoretic {$p$}-{L}aplacian on a ball},
  date={2014},
  journal={Differential Integral Equations},
  volume={27},
  number={7-8},
  pages={659\ndash 670},
}

\bib{kuhn}{article}{
  author={K\"{u}hn, {M.}},
  title={Power- and log-concavity of viscosity solutions to some elliptic {D}irichlet problems},
  date={2018},
  journal={Commun. Pure Appl. Anal.},
  volume={17},
  number={6},
  pages={2773\ndash 2788},
}

\bib{LuWang2008}{incollection}{
  author={Lu, {G.}},
  author={Wang, {P.}},
  title={A uniqueness theorem for degenerate elliptic equations},
  date={2008},
  booktitle={Geometric methods in {PDE}'s},
  series={Lect. Notes Semin. Interdiscip. Mat.},
  volume={7},
  publisher={Semin. Interdiscip. Mat. (S.I.M.), Potenza},
  pages={207\ndash 222},
  review={\MR {2605157}},
}

\bib{MPR1}{article}{
  author={Mart\'{i}nez-Aparicio, {P.J.}},
  author={P\'{e}rez-Llanos, {M.}},
  author={Rossi, {J.D.}},
  title={The limit as {$p\rightarrow \infty $} for the eigenvalue problem of the 1-homogeneous {$p$}-{L}aplacian},
  date={2014},
  journal={Rev. Mat. Complut.},
  volume={27},
  number={1},
  pages={241\ndash 258},
}

\bib{MPR2}{article}{
  author={Mart\'{i}nez-Aparicio, {P.J.}},
  author={P\'{e}rez-Llanos, {M.}},
  author={Rossi, {J.D.}},
  title={The sublinear problem for the 1-homogeneous {$p$}-{L}aplacian},
  date={2014},
  journal={Proc. Amer. Math. Soc.},
  volume={142},
  number={8},
  pages={2641\ndash 2648},
}

\bib{PSSW1}{article}{
  author={Peres, {Y.}},
  author={Schramm, {O.}},
  author={Sheffield, {S.}},
  author={Wilson, {D. B.}},
  title={Tug-of-war and the infinity {L}aplacian},
  date={2009},
  journal={J. Amer. Math. Soc.},
  volume={22},
  number={1},
  pages={167\ndash 210},
}

\bib{PSSW2}{article}{
  author={Peres, {Y.}},
  author={Sheffield, {S.}},
  title={Tug-of-war with noise: a game-theoretic view of the {$p$}-{L}aplacian},
  date={2008},
  journal={Duke Math. J.},
  volume={145},
  number={1},
  pages={91\ndash 120},
}

\bib{Sak}{article}{
  author={Sakaguchi, {S.}},
  title={Concavity properties of solutions to some degenerate quasilinear elliptic {D}irichlet problems},
  date={1987},
  issn={0391-173X},
  journal={Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4)},
  volume={14},
  number={3},
  pages={403\ndash 421 (1988)},
  url={http://www.numdam.org/item?id=ASNSP_1987_4_14_3_403_0},
  review={\MR {951227 (89h:35133)}},
}

\end{biblist}
\end{bibdiv}

\end{document}

\bib{BiDe2004}{article}{
  author={Birindelli, {I.}},
  author={Demengel, {F.}},
  title={Comparison principle and {L}iouville type results for singular fully nonlinear operators},
  date={2004},
  issn={0240-2963},
  journal={Ann. Fac. Sci. Toulouse Math. (6)},
  volume={13},
  number={2},
  pages={261\ndash 287},
  url={http://afst.cedram.org/item?id=AFST_2004_6_13_2_261_0},
  review={\MR {2126744}},
}


\section{Introduction \label{intro}}
Hot, dense stellar matter appears in  explosive astrophysical phenomena
 such  as core-collapse supernovae and neutron-star mergers,
 but its equation of State (EOS) and weak-interaction rates are not well understood.
The EOS determines the thermodynamic quantities and composition
in terms of nucleons and various nuclei, as functions of temperature, density, and charge fraction.
In the near future, gravitational waves from neutron-star mergers may provide us
 with information about thermodynamic quantities such as pressure at supranuclear density \cite{abott17,shibata17}.
In contrast,
neutrinos from core-collapse supernova may enable investigations of
the nuclear composition
of  the EOS and the weak interaction rates at subnuclear densities.

Nuclear electron captures and coherent neutrino-nucleus scattering
during the collapse of a  massive star play important roles in supernova dynamics and neutrino observations.
In particular, the evolution of the charge fraction and entropy in the collapsing core
is very sensitive to weak interactions  involving  nuclei that appear at densities of  $\rho_B \approx10^{11}\textrm{--}10^{12}$~g/cm$^3$ immediately before the onset of neutrino trapping.
As a result, the nuclear weak interactions potentially affect the total mass of a newborn neutron star, the peak neutrino luminosity, and the strength of a shock wave \cite{sullivan16}.
Recently, some groups have suggested that medium-mass, neutron-rich nuclei may dominate
 the deleptonization process \cite{sullivan16,furusawa17b,raduta16,raduta17,titus18},
although the approach used to determine the nuclear composition or the nuclear model in the EOS  differs from group to group.

Several investigations have compared different  EOS models \cite{buyukcizmeci13,hempel15, oertel17}
but have not been systematic in their approaches; therefore, it is difficult to obtain quantitative conclusions.
Many research groups have  already compared modern statistical EOSs that include full ensemble of nuclei, and classical single-nucleus EOSs in that a representative nucleus is used. The latter have been shown to exhibit defects  \cite{burrows84,souza09,gulminelli15,furusawa17b,furusawa17c}; accordingly,
this study focuses only on the former in this work.
The main inputs for the statistical EOSs are experimental data
and theoretical estimations for uniform nuclear matter, nuclear masses, and nuclear excitations, which are related to each other.

Bulk properties, such as incompressibility, characterize uniform nuclear matter,
and their roles in supernovae have been  discussed previously, with the primary focus being the stiffness
 of nuclear matter around nuclear densities.
It is known that softer bulk properties lead to faster contractions of proto-neutron stars,
 higher neutrino luminosities, and the generation of more acoustic waves,
 resulting in more energetic explosions (see, e.g., Refs. \cite{suwa13,nagakura18a}).
 Most supernova EOSs employ relativistic mean-field calculations  \cite{hempel10,typel10,furusawa11,furusawa13a,furusawa17a,sheng11a,sheng11b,steiner13} or Skrme-type interactions \cite{schneider17}, using fixed parameter sets.
Recently, some realistic nuclear EOSs  based on a variational method that employs the bare nuclear potentials have also been presented \cite{togashi17,furusawa17d}.

Nuclear masses are another essential ingredient for the statistical EOSs as the nuclear composition
 at low densities and temperatures  is sensitive  to the nuclear binding energies.
 During  core collapse, not only do nuclei close to the stability line  appear but so do extremely heavy and/or neutron-rich nuclei;
we require theoretical mass data for such nuclei.
Blinnikov {\it et al.} \cite{blinnikov11} and the FYSS EOSs \cite{furusawa11,furusawa13a,furusawa17a} utilize the KTUY mass data \cite{koura05},
whereas the
SHO EOS \cite{sheng11a,sheng11b} adopts the old version of FRDM \cite{moller95}.
In addition, the authors of the HS EOSs \cite{hempel10,steiner13} have prepared various EOSs based on different theoretical mass data; see, e.g., Refs. \cite{geng05, moller95,roca08,lalazissis99}.

At finite temperatures, nuclear excitations must be  included in the EOS.
They increase the internal degrees of freedom, especially for large and non-magic nuclei, which have nucleons in open-shells.
Ideally, in  the statistics,  we should precisely count all the states with individual excitation energies, for all nuclei; however,
 experimental data for the nuclear level densities are insufficient.
In addition, theoretical calculations, such as shell-model calculation  \cite{tsunoda17},   for all states would be quite difficult.
For now, we employ two approaches to take account of the finite-temperature effects.
One introduces internal partition functions as a function of temperature
in determining the  number densities of nuclei.
The other approach phenomenologically represents the ground state and all excited states by a hot nucleus
so that the nuclear internal free energy depends on the temperature and is identical  to nuclear mass at zero temperature.

The former approach allows us to partly  use  the experimental data on level  densities directly
if the temperature is not too high and  the excitations are weak.
Ishizuka {\it et al.} \cite{ishizuka03} and the HS EOS employ a semi-empirical expression for internal partition functions provided by  Fai and Randrup \cite{fai82}, which depends only on the temperature and mass number.
The SRO EOS  \cite{schneider17} utilizes tabulated partition functions provided by Raucher \cite{rauscher97,rauscher00,rauscher03}, which is  based on  a Fermi-gas approach, using experimental data for level densities and the FRDM mass data \cite{moller95} to evaluate the separation energies of nucleons.
The latter approach employs a temperature-dependent internal free energy, such as in  the SMSM (statistical model for supernova matter) EOS  \cite{botvina04, botvina10,buyukcizmeci14}.
The SMSM EOS is an extension of the SMM (statistical multifragmentation model) EOS \cite{bondorf95},
 which reproduces  fragmentations in low-energy heavy-ion collisions.
The SHO EOS also assumes similar excitation effects for the bulk free energy of the nuclei.
This approach may be valid for reproducing the ensemble of hot nuclei when the temperature is sufficiently high to make the excited states dominant.
Ground-state properties, such as neutron magic numbers and pair energies, are considered to be almost washed out at $T\approx$ 2--3~MeV \cite{brack74, bohr87, sandulescu97,nishimura14}, which are temperatures typical of the late phases of core collapse.
The FYSS EOS \cite{furusawa17a,furusawa17d} employs
 temperature
interpolation between  the internal free energies of cold and hot nuclei.

In the present work, I investigate the sensitivity of the nuclear composition to several nuclear uncertainties, the bulk properties of nuclear matter, the nuclear mass data,
and the finite-temperature modeling in order to clarify ambiguities in the supernova EOS at subnuclear densities
and to organize the target properties and nuclides  for future theoretical and experimental studies.
The organization of this paper is as follows:
In Sec.~\ref{sec:model}, I describe the formulations and inputs for the EOSs.
Comparisons among  EOSs under thermodynamics conditions during core collapse of massive 15M$_{\odot}$ and 25M$_{\odot}$ stars are discussed in Sec.~\ref{sec:res}.
Finally, a summary is given in Sec.~\ref{sec:conc}

\section{Nuclear physics inputs for statistical equations of state \label{sec:model}}
The free-energy density is minimized with respect to the number densities of all particles,
subject to the constrains of baryon and charge conservation,
to obtain the EOS as a function of the volume-averaged values of the baryon number density $n_B$,
  temperature $T$, and charge fraction $Y_p$.
I systematically change three ingredients and
 Table~\ref{tab_model} summarizes the specific cases.
Below, I explain the 14 models---
1B, 1D, 1E, 2FB, 2FD, 2FE, 2HE, 2KE, 2WE, 3FE, 3HE, 3KE, 3WE, and 4FE
---for which the number and abbreviations denote, in order,
the treatment of finite-temperature effects for nuclear excitations,
the choice of mass data for models 2--4, and
the parameter set for bulk nuclear matter.

The constituents of nuclear matter under nuclear statistical equilibrium are the dripped nucleons and all nuclei.
The free-energy density can be given by
\begin{eqnarray}
f  = \xi \left[ n'_p m_p+n'_n m_n +(n'_p+n'_n) \omega \left(n'_p+n'_n,T,\frac{n'_p}{n'_p+n'_n}\right) \right] \nonumber \\
 +  \sum_{AZ}  n_{AZ} \left\{ T \left[ \ln \left(\frac{ n_{AZ}/\kappa}{g_{AZ} (M_{AZ} T/2\pi \hbar ^2 )^{3/2}  }\right)- 1 \right]  + M_{AZ} \right\},
\label{total}
\end{eqnarray}
where $n'_p$ and $n'_n$  are the local
number densities of protons and neutrons, respectively, and $m_p$ and $m_n$, are the corresponding masses
\cite{furusawa18a}.
The volume factor for free nucleons is defined as $\xi=V'/V$ with the vapor volume $V'$ and the total volume $V$.
The local free energy density of the dripped nucleons
is obtained from a nuclear-matter calculation of $\omega(n_B,T,x)$
 as a function of the local values of  $n_B$, $T$, and the charge fraction $x$,
 as in Furusawa and Mishustin \cite{furusawa18a}.
The vapor volume can be calculated as $V'=V(1-\sum_{AZ} n_{AZ} A/n_{sAZ})$.
The summation covers all nuclides listed in the tabulated data, theoretical mass data or partition functions, as discussed later and shown in Fig.~\ref{fig_nuclide}.
The quantities $M_{AZ}$, $g_{AZ}$, $n_{AZ}$, and $n_{sAZ}$ are the mass, internal degree of freedom, number density, and saturation density of the nucleus with the mass number $A$ and the proton number $Z$.
The translational energy of the nuclei is calculated as for an ideal Boltzmann gas, with excluded volume effects represented by  $\kappa=1-n_B/n_0$.

\subsection{Nuclear matter calculation for dripped nucleons and bulk energies of nuclei}
The bulk properties at the nuclear saturation density  $n_0$ for symmetric matter are
characterized by the following parameters: the binding energy at saturation $\omega_0$, the incompressibility $K_0$,
 symmetry energy at saturation $S_0$, and  symmetry energy slope parameter $L$.
I use the three parameter sets \cite{oyamatsu03,oyamatsu07} listed in Tab~\ref{tab_bulk}.
Parameters B and E share the same value of $K_0$, while parameters D and  E are based on the same value of the saturation slope parameter  $y=-K_0 S_0/(3 n_0 L)$, which gives  the saturation point of asymmetric nuclear matter at about $x=0.5$ as $n_0 +  (1-2x)^2 S_0/y$.
The parameters other than $K_0$ and $y$ are optimized to reproduce the masses and radii of stable nuclei by using Thomas--Fermi calculations \cite{oyamatsu03,oyamatsu07}.
The free energy per baryon for uniform nuclear matter, $\omega(n_B,T,x)$, is
 based on the  finite-temperature calculation \cite{furusawa18a} of  Oyamatsu and Iida \cite{oyamatsu03,oyamatsu07}
and is given by:
\begin{eqnarray}%
 \label{eq:para1}
 \omega(n_B,T,x) &= &  \omega_{int}+\omega_{kin} , \\
\omega_{int}&=&4 x(1-x) v_s(n_B) /n_B + (1-2 x)^2 v_n(n_B) /n_B ,  \nonumber \\
\omega_{kin}&=& \frac{T}{2 \pi^2 n_B}  \left(\frac{2m_pT}{\hbar^2} \right)^{3/2}   F_{3/2} (\eta_p) +   \frac{T}{2 \pi^2 n_B}  \left( \frac{2m_nT}{\hbar^2} \right)^{3/2}  F_{3/2} (\eta_n)  , \nonumber \\
v_s(n_B) & = &a_1 n_B^2 + \frac{a_2 n_B^3}{1+a_3 n_B}   , \\
v_n(n_B) & = &b_1 n_B^2 + \frac{b_2 n_B^3}{1+b_3 n_B},
\end{eqnarray}
where
$\omega_{int}$ and $\omega_{kin}$ are interaction and kinetic energies,  $v_s$ and $v_n$ are the energy densities for symmetric matter and pure-neutron matter, respectively, $b_3=$1.58632~fm$^3$ and
$\eta_{p/n}=\left[ \mu^0_{p/n} - \partial (\omega_{int} n^0_{p/n})/ \partial n^0_{p/n} \right]/T
$.
The number densities of protons, $x n_B$, and neutrons, $(1-x) n_B$, correspond to the chemical potentials, $\mu^0_{p/n}$,
and $F_{k}(\eta)$ is defined by $F_{k}(\eta)=\int_0^{\infty} u^{k} \left[ 1+\exp{(u-\eta)} \right]^{-1} du $.
The parameters $a_1$, $a_2$, $a_3$, $b_1$, and $b_2$ for the potential energies
 are fitted to reproduce the bulk parameters in Table~\ref{tab_bulk} \{see Eqs.~(10--14) in Ref.~\cite{furusawa18a}\}.

Figure~\ref{fig_blk} presents the free-energy densities for the three parameter sets.
Differences appear at $n_B\approx 10^{-2}$~fm$^{-3}$.
For symmetric matter  with $x=0.5$, the bulk free energies for parameter sets B and E are almost identical to each other as they have the same value of  $K_0$.
In contrast,  the parameter set D has the soft property of a more gradual
 growth of $\omega$ because of the smaller value of  $K_0$.
For asymmetric matter with $x=0.1$,
the smaller values of $L$ make the densities higher, at which $\omega$ increases steeply.

\subsection{Nuclear mass data and model for internal free energy}
The nuclides taken into account in the statistical EOSs are presented in Fig~\ref{fig_nuclide}.
I utilize experimental mass data from AME12  \cite{audi12}
 for those nuclei with data available and  theoretical nuclear mass data from FRDM \cite{moller16}, HFB24 \cite{goriely13}, KTUY \cite{koura05}, or WS4 \cite{wang14}  for the other nuclei.
The middle character---F, H, K, or W---in the names of models 2-4 denotes the specific dataset.
These data also optimize the surface energy parameters for model 2, as discussed later.
Model 1 includes the nuclides  with FRDM  mass data available, although I do not use the values of nuclear masses.

In nuclear statistical equilibrium,
the number densities of all nuclei are given by
\begin{eqnarray}
n_{AZ}  =  \kappa \sum_i g_{AZi} \left(\frac{M_{AZ} T}{2\pi \hbar ^2}  \right)^{3/2} {\rm exp} \left(\frac{\mu_{AZ}-M_{AZ}}{T} \right)  ,
\label{eq:num0}
\end{eqnarray}
where  the index $i$ runs over the ground state and all excited states, $g_{AZi}$
is the degree of freedom of state $i$,
$\mu_{AZ}= (A-Z) \mu_n + Z \mu_p$, and
 $\mu_n$ and $\mu_p$ are the chemical potentials of neutrons and protons, respectively.
Unfortunately, the experimental data of  the excitation energies or level densities are insufficient at present.
Instead of including  all excited states explicitly, we often introduce
 internal degrees of freedom, $g_{AZ}(T)$, and/or internal free energy, $F_{AZ}(T)$, as a function of $T$:
\begin{eqnarray}
n_{AZ}  & \simeq &  \kappa g_{AZ}(T) \left(\frac{M_{AZ}T}{2\pi \hbar ^2}  \right)^{3/2} {\rm exp} \left(\frac{\mu_{AZ}-M_{AZ}}{T} \right)  , \label{nd34} \\
          & \simeq & \kappa g_{AZ}^s \left(\frac{M_{AZ}T}{2\pi \hbar ^2}  \right)^{3/2} {\rm exp} \left(\frac{\mu_{AZ}-F_{AZ}(T)}{T} \right).  \label{nd12}
\end{eqnarray}
The former equation can be obtained from the relations $\mu_{AZ}=\partial f/ \partial n_{AZ}$ and Eq.~(\ref{total}).
Models 3 and  4  employ
 $g_{AZ}(T)$ and use cold nuclear masses $M_{AZ}$ that include finite-density effects only,
as in  Eq.~(\ref{nd34}) and explained later.
In the models 1 and 2, I substitute
$F_{AZ}(T)$
and simplified factors $g_{AZ}^s$
for $M_{AZ}$ and $g_{AZ}(T)$, as in  Eq.~(\ref{nd12}).
For models 1B, 1D and 1E,
I employ a statistical model of the SMSM EOS \cite{botvina04, botvina10,buyukcizmeci14}
 with the nuclear bulk parameters and the formulation of Coulomb energies modified.
Models 2FB, 2FD, 2FE, 2HE, 2KE,  and 2WE, are based on
the liquid-drop model (LDM) \cite{furusawa18a} with the shell washout effect introduced in the FYSS EOS \cite{furusawa11,furusawa13a,furusawa17a,furusawa17d}.

The internal free energy of model 1 with proton numbers $Z>5$ and neutron numbers  $N>5$
 is composed of nucleon rest masses and bulk, Coulomb, and  surface  energies as follows:
\begin{eqnarray}
F_{AZ} &=& Z m_p + (A-Z) m_n +  F_{AZ}^{B} +  F_{AZ}^{C} + F_{AZ}^{Sf},  \label{eq:mass1} \\
F_{AZ}^{B}&=&  A \left( -\omega_0 - \frac{T^2}{\epsilon} \right) +S_0 \frac{(A-2Z)^2}{A} , \\
F_{AZ}^C&=&\displaystyle{\frac{3}{5}\left(\frac{3}{4 \pi}\right)^{-1/3}  e^2 n_{sAZ}^2 \left(Z/A - n'_p/n_{sAZ} \right)^2 \left({A /n_{sAZ}}\right)^{5/3} D(u_{AZ})} , \label{eq:cl}  \\
F_{AZ}^{Sf}&=& 18 \left( \frac{T^2_{cAZ}-T^2}{T^2_{cAZ}+T^2}\right)^{5/4} A^{2/3}.
\end{eqnarray}
In Eq.~(\ref{eq:cl}),
 $D(u_{AZ})=1-\frac{3}{2}u_{AZ}^{1/3}+\frac{1}{2}u_{AZ}$, where the filling factor $u_{AZ}=(n_e-n'_p)/(Z n_{sAZ} /A -n'_p) $,
 $n_e(=Y_p n_B)$ is the electron number density, and   $e$ is the elementary charge.
In the model, I ignore the iso-spin dependencies of the critical temperature for nuclear vaporization and
 nuclear saturation densities as $T_{cAZ}=18$ MeV and $n_{sAZ}=n_0$.

In model 2,

I consider nuclear shell effects, $F_{AZ}^{Sh}$, and dependencies of $n_{sAZ}$ on $T$ and on $Z/A$ and employ other formulations for bulk and surface energies.
The internal free energy is given by
\begin{eqnarray}
F_{AZ} &=& Z m_p + (A-Z) m_n +  F_{AZ}^{B} +  F_{AZ}^{C} + F_{AZ}^{Sf} +   F_{AZ}^{Sh},  \label{eq:mass2} \\
F_{AZ}^{B}&=&  A \{\omega(n_{sAZ},T,Z/A) \},  \label{eq:bulk} \\
F_{AZ}^{Sf}&=&\left(\displaystyle{\frac{36 \pi A^2}{n_{sAZ}^2}} \right)^{1/3}  \sigma_{AZ}(T)
\left(1-\displaystyle{\frac{n'_p+n'_n}{n_{sAZ}}} \right)^2 ,
 \label{eq:sf}  \\
 F_{AZ}^{Sh}&=&F_{AZ0}^{Sh} \displaystyle{\frac{\tau_{AZ}}{{\rm sinh}\tau_{AZ}}} . \label{eq:sh}
\end{eqnarray}
The quantity $F_{AZ}^B$ is based on the same calculations for dripped nucleons as Eq.~(\ref{eq:para1}), and
$n_{sAZ}(T)$ is defined as the density at which the free energy, $\omega (n_B,T,Z/A)$, takes its local minimum value around $n_0$.
The Coulomb energy, $F_{AZ}^{C}$, is calculated by  Eq.~(\ref{eq:cl}).
For $F_{AZ}^{Sf}$ in Eq.~(\ref{eq:sf}), $\sigma_{AZ}(T) =\sigma_0 [(T^2_{cAZ}-T^2)/(T^2_{cAZ}+T^2)]^{5/4} (16 + C_s) /[(1-Z/A)^{-3} + (Z/A)^{-3}  +C_s ] $ and
 $T_{cAZ}$, is defined by using  the bulk pressure, $P_{bulk}=n_B^2 \partial \omega(n_B,T,x=Z/A)/\partial n_B$.
This is  the temperature at which
$(\partial P_{bulk}/\partial n_B)|_{x}=0$ and $(\partial^2 P_{bulk}/\partial n_B^2)|_{x}=0$  simultaneously and  is shown in Fig.~\ref{fig_ct} for each bulk parameter set.
The last factor in Eq.~(\ref{eq:sh}) expresses the washout effects, where $\tau_{AZ} =2 \pi^2 T/(41 A^{-1/3})$ \cite{furusawa13a,furusawa17a}.

The shell energies of the ground states are set equal to the mass data, $M_{AZ}^{data}$,  minus the
LDM mass-energy in a zero-density medium, i.e., $ F_{AZ0}^{Sh}=M_{AZ}^{data} - [ Z m_p + (A-Z) m_n +  F_{AZ}^{B} +   F_{AZ}^{C} + F_{AZ}^{Sf} ]_{n_e=0,n'_p=0,n'_n=0,T=0}$.
This assumption allows the internal free energy of Eq.~(\ref{eq:mass2}) to  reproduce  exactly the mass data at low temperatures and densities.
Here, I assume  the shell energy to be positive to avoid  negative entropy production.
Nuclear excitations usually increase the internal degrees of freedom,
which corresponds to a free-energy reduction associated with a temperature increase and to positive entropy production.
This assumption for the sign of shell
correction, however, is not typical and, in some references \cite{egidy05}, the shell energies of magic nuclei take negative values.
The surface tension $\sigma_0$ and isospin-dependence parameter $C_s$ are optimized to minimize the total deviation of the LDM mass-energy (internal free energy in zero-temperature limit) per baryon from the mass data or the sum of the shell energies per baryon, $\sum_{Z>5,n>5}F^{Sh}_{AZ0}/A$.
The resulting values of these quantities are listed in Table~\ref{tab_model}.

Figure~\ref{figshll} illustrates the temperature dependence of the shell energies for models 2FD, 2FE, 2KE, and 2WE. The combination of mass data and bulk parameters leads to a set of shell energies and surface tension parameters. At $T=3$ MeV, the nuclear shell effects are significantly reduced, especially for nuclei with large mass numbers.
Note that this formulation is still rough and may overestimate the shell damping  as discussed later.

In models 3 and 4,
 the temperature dependence  is encapsulated in internal degrees of freedom $g_{AZ}(T)$
as in  Eq.~(\ref{nd34}).
To evaluate the nuclear masses, I take into account only Coulomb screening effects: $\Delta E^C_{AZ}=F_{AZ}^C(n_e,n'_p)-F_{AZ}^C(0,0)$ from Eq.~(\ref{eq:cl}),
with $M_{AZ}=M_{AZ}^{data} +\Delta E^C_{AZ}$.
In models 3 and 4,  I set $n_{sAZ}=n_0$ for all nuclei.

\subsection{Internal degrees of freedom}
In model 3, I use a semi-empirical function for internal degrees of freedom
 that depends only on the temperature and  mass number \cite{fai82}.
It is utilized in the HS EOSs \cite{hempel10,steiner13}
 and can be given by
\begin{equation}
\label{eq:ex}
g_{A}(T)=g_{A}^0 +\frac{c_1}{A^{5/3}}\int_0^{16.2 A} dE e^{-E/T}\exp\left(\sqrt{2 a(A) E}\right) , \label{eq:fai}
\end{equation}
where $a(A)=(A/8)(1-c_2 A^{-1/3})$MeV$^{-1}$, $c_1=0.2$MeV$^{-1}$, $c_2=0.8$,
and $g_{A}^0=1$ for even nuclei and $g_{A}^0=3$ for odd nuclei, and
the upper bound of the integral is set to be a typical bulk energy, $16.2 A$, as in the HS EOS.
There are several options for the upper limit according to
each model of energy integral: infinity; nuclear binding energies; the smaller of
neutron and proton separation energies with in-medium modifications \cite{gulminelli15,fowler78}.
The choice of the upper bound, however,  does not significantly change the results.
Model 4FE employs Raucher's tabulated data for the internal partition functions $g_{AZ}(T)$ for each nucleus \cite{rauscher03}, which are used in the SRO EOS  \cite{schneider17}.
The nuclei listed in these data are limited, as shown in Fig~\ref{fig_nuclide};
therefore, I assume no excitations for the nuclei,
for which partition function data are unavailable, even if the mass data are available.

In  models 1 and 2,
I consider the excited-state contributions to the quantities
 of $F_{AZ}^{B}$,  $F_{AZ}^{Sf}$, and  $F_{AZ}^{Sh}$  in Eqs.~(\ref{eq:mass1}) and (\ref{eq:mass2}).
The values of $g_{AZ}^s$ in model 1 are set to be unity as in the original EOS.
In model 2, the ground-state spin factors, $g_{AZ}^0$, are taken from Ref. \cite{rauscher00} for those nuclei for which data are available and are set to $g_{AZ}^0=1$ for other even nuclei and $g_{AZ}^0=3$ for other odd nuclei.
The internal degrees of freedom
are assumed to be washed out according to $g_{AZ}^s= (g_{AZ}^0 -1) \tau_{AZ}/{\rm sinh}\tau_{AZ}  +1$, approaching unity with
the reduction in the shell energy because of increase in the excited states \cite{furusawa17a}.

For comparison with the internal degrees of freedom for models 3 and 4, I define effective internal degrees of freedom for  models 1 and 2 in refer to Hempel \cite{hempel10d}:
 \begin{eqnarray}
 g^*_{AZ}(T) =g_{AZ}^s {\rm exp} \left(\frac{F^0_{AZ}-F^*_{AZ}}{T} \right)  . \label{eqeff}
\end{eqnarray}
Here, $F^*_{AZ}$ is the internal free energy in the zero-density-medium limit,
 but at finite temperature, $F_{AZ}(T)=F_{AZ}(T,n_e=0,n'_p=0,n'_n=0)$.
The internal free energy in vacuum limit ($T=0$, $n_e=0$, $n'_p=0$, and, $n'_n=0$),  $F^0_{AZ}$, is identified with $M_{AZ}^{data}$ in model 2.
This formulation leads to the expression for the number density of nuclei in the zero-density-medium limit:
 \begin{eqnarray}
n_{AZ}(T) =  g^*_{AZ}(T) \left(\frac{F_{AZ}^{0} T}{2\pi \hbar ^2}  \right)^{3/2} {\rm exp} \left(\frac{\mu_{AZ}-F_{AZ}^{0}}{T} \right)  ,
\end{eqnarray}
which is similar to Eq.~(\ref{nd34}) for models 3 and 4.
The internal degrees of freedom  for these four approaches are presented in Fig.~\ref{figlvl}.
Model 1 generally agrees with Rasucher's data (model 4)  at $T=1$~MeV and
tends to exhibit fewer internal degrees of freedom than model 2 at any temperature.
For model 2, nuclei with small mass numbers are more likely to be excited,
as compared with the other models,
as their shell and surface energies per baryon are large relative to nuclei with large mass numbers.
The semi-empirical internal degrees of freedom for model 3 are clearly fewer  than those given by
the other three models.

For all models, light clusters with  $Z \leq5$ or  $N \leq5$ are calculated as ideal Boltzmann gases using the  excluded volume, $\kappa$,
  Coulomb screening effects,  $F_{AZ}=M_{AZ}=M_{AZ}^{data} +\Delta E^C_{AZ}$, and no excitation, $g_{AZ}=g_0$ as in Ref. \cite{furusawa18a}.

\section{nuclear statistical equilibrium in collapsing cores \label{sec:res}}
I next compare some models listed in Table~\ref{tab_model}
 along core-collapse trajectories for stars of 15M$_\odot$ and 25M$_\odot$ \cite{juodagalvis10}.
Figure~\ref{fig_pro} presents the temperature and  charge fraction as a function of density for  these supernova progenitors.
The progenitor with 25M$_{\odot}$ exhibits higher temperature and lower charge fractions for the same density.
Note that these quantities depend on the EOS and the results of core-deleptonization in actual core-collapse simulations, although I assume the same thermodynamic conditions for all EOSs in this comparison.
I focus mainly on the nuclear composition of heavy nuclei at $\rho_B \approx10^{11-12}$~g/cm$^3$,
because weak interactions involving heavy nuclei play fundamental roles in core-deleptonization immediately before neutrino sphere formation,
and the impacts of nucleons and light clusters are not dominant.

Figures~\ref{fig_dis25} and \ref{fig_dis15} display
the mass fractions of all heavy nuclei, $X_{AZ}=A n_{AZ}/n_B$,
 for models 1E, 2FE, 3FE, and 4FE under thermodynamical conditions
at $(\rho_B, T, Y_p)=(1.9 \times 10^{11}$ g/cm$^3$, 1.3 MeV, 0.36)  for the 25M$_\odot$-star and
$(2.0 \times 10^{12}$ g/cm$^3$, 1.8 MeV, 0.28) for the 15M$_\odot$-star.
Model 1E shows smoother mass distribution than the other models,
because of the lack of nuclear shell effects.
Other models generally have a monomodal or bimodal structure around neutron magic numbers ($N=$28, 50, and 82).
 Model 3FE assumes the semi-empirical function for internal degrees of freedom and cold nuclear mass and, as a result,
overestimates the mass fractions of these magic nuclei.
In model 2FE, the shell energies $F_{AZ}^{Sh}$ are washed out,
 and, consequently, non-magic nuclei are also abundant.
In model 4FE, Rauscher's data for the internal degrees of freedom take into account the shell effects; therefore, non-magic nuclei are more easily  excited and have larger values of $g_{AZ}(T)$ than the magic nuclei.
Therefore,  models 2FE and 4FE produce  wider and smoother distributions
 of  nuclei with large abundances $X_{AZ}$ in the  $(N,Z)$ plane than model 3FE.
Unfortunately, model 4FE lacks the partition function data for some important nuclei close to the abundance peaks, e.g.,
$(N,Z)=(85, 28)$ and (75, 35), as displayed in Fig.~\ref{fig_nuclide}.

The average mass numbers and proton numbers of heavy nuclei are
shown in Fig.~\ref{fig_maz} for almost all the models.
The impact on these quantities of the mass data and bulk properties is not large,
 while the choice of partition function alters them greatly.
The average mass numbers for models 1B, 1D, and 1E are smaller and grow more gradually  than the other models  owing to the lack of shell effects.
For models 2FB, 2FE, and 2KE, nuclei with smaller mass numbers
are more highly excited than those in the other models,
 as discussed above and shown in Fig.~\ref{figlvl}.
Therefore, they tend to produce smaller average mass numbers than do models 3 and 4.
As discussed above and displayed in Figs.~\ref{fig_dis25} and \ref{fig_dis15},
nuclear shell effects remain at high temperatures in model 3FE.
Hence, only that model leads to a step-wise growth
 in the average mass numbers and proton numbers,
owing to the neutron-magic numbers: at about  $\rho_B=10^{12}$~g/cm$^3$,
the dominant constituents of nuclear matter abruptly change nuclei with $N\approx50$ to those with  $N\approx82$.

The bulk parameter does not affect significantly the results among models 2.
This is due to  the internal free energies being almost equal to the mass data at low temperatures owing to the introduction of shell energies
 and to the parameters for nuclear bulk energies other than $K_0$  and $y$ and for  surface tensions also being optimized to reproduce the similar nuclear properties, such as nuclear masses and radii.
The differences arising from the bulk properties are more visible among models 1.
 The larger values of $S_0$ and $L$  reduce
the binding energies and the average mass numbers
in models 1B and 2FB compared with  1D, 1E, and  2FE.
The choice of theoretical mass data does not significantly influence the results either
in models 2, as nuclei for which the experimental mass data exist are dominant.
 Differences from the theoretical mass data appear only for neutron-rich nuclei at high temperatures, where the shell effects derived from the mass data no longer survive.
Models 3FE, 3HE, 3KE, and 3WE display larger differences  due to the mass data than models 2KE and 2FE,
because they do not take into account the shell damping.

Figure~\ref{fig_xh} shows the total mass fraction of heavy nuclei, $X_h=\sum_{Z>5,N>5}X_{AZ}$.
It is found that the mass fraction is generally large for a model that has many
internal degrees of freedom.
The difference is greater for the 25M$_\odot$ star in which the central temperature  is higher than that in the 15M$_\odot$ star.
Model 4FE always reproduces larger mass fractions than model 3FE.
At the beginning of core collapse, models 2FB, 2FE, and 2KE yield larger mass fractions than models 3FE or 4FE
 because  of the larger internal degrees of freedom,
 while the mass fraction for model 4FE becomes larger than those of models 2
 at $\rho_B\approx10^{12}{\rm g/cm}^3$,  because  nuclei with large mass numbers  are
 populated in model 4FE, as shown in Figs.~\ref{figlvl} and \ref{fig_dis15}.
The mass fractions for models 1B, 1D, and 1E deviate from those for the other models,  because shell effects are not considered.
At low densities in models 2 and 3, the choice of mass data hardly affects the mass fraction
 because of the fact that the internal free energies of the nuclei are almost identical to
 their experimental values of nuclear masses and  the nuclei are not excited at low temperatures.
At $\rho_B\approx10^{12}{\rm g/cm}^3$, slight differences arising from the mass data appear.
It is found that softer bulk properties increase the mass fractions of heavy nuclei,
in the order of models 1D, 1E, and  1B  (2FD, 2FE, and 2FB),
which corresponds to the ascending order of  $S_0$ or $L$.

The central densities and temperatures of  the collapsing cores increase adiabatically,
and the entropy is an essential quantity for determining the dynamics of core collapse.
I therefore compare the baryonic entropies per baryon in Fig.~\ref{fig_ent}.
These are given by $s_B=-\partial f/\partial T|_{n_B,Y_p}/n_B$ and are essentially determined by the kinetic terms of the nuclear components and the temperature derivatives of the internal degrees of freedom, $\partial g_{AZ}(T)/\partial T$, and  of  the internal free energy, $\partial F_{AZ}(T)/\partial T$.
Note from the figure that the entropies differ among the models based on
different finite-temperature modeling, even in  the initial stages of core collapse.
For the models with smaller $X_h$, such as for models 1B, 1D, and 1E,
the entropy is more likely to be high
because of  the increase in the population of dripped neutrons.
Models 2FB, 2FE, and 2FK at  at $\rho_B\approx10^{10}{\rm g/cm}^3$ and model 4FE at $\rho_B\approx10^{12}{\rm g/cm}^3$  tend to yield  high entropies
because of the large values of  $\partial g^*_{AZ}(T)/\partial T$ and  $\partial g_{AZ}(T)/\partial T$.
 These differences in entropy would be  influential for the fate of core-collapse supernovae  \cite{suwa16, furusawa17d,nagakura18b}.

\section{Summary and Discussion \label{sec:conc}}
I have compared the 14 statistical EOSs for the core collapse of a massive star by systematically changing the bulk properties of nuclear matter,
the masses of neutron-rich nuclei, and  the treatment of finite-temperature effects in the nuclear model.
Overall, the temperature dependencies of the internal degrees of freedom
 and of the internal free energy are paramount  to determining the entropy and nuclear composition during the core-collapse phase.
 Differences in these quantities among EOSs would affect thermodynamic conditions
 and ensemble-averaged weak-interaction rates of nuclear electron-captures and neutrino-nucleus coherent scattering \cite{furusawa17b};
 therefore, the dynamics of the core collapse would be sensitive to the treatment of the finite-temperature effects.
The parameters of bulk nuclear matter and the theoretical mass data  for neutron-rich nuclei do not greatly influence the average  mass  number or the total mass fraction of heavy nuclei.

Semi-empirical   expression for internal degrees of freedom that ignore the temperature dependence of shell effects may not be appropriate for discussions of core-collapse nuclei,
because they overestimate  the number densities of magic nuclei.
More precise calculations for internal partition functions are required.
The individual level densities of medium-mass, neutron-rich nuclei  however,  have not
been well studied.
Actually, some important nuclei  lack the data \cite{rauscher03}
used in the present work.
Phenomenological models with nuclear shell washout
 tend to increase the populations of nuclei with small mass numbers that have
large electron-capture rates; they would therefore reduce the charge fraction and entropy significantly in collapse simulations \cite{furusawa17b,nagakura18b}.
The washout formulation for the nuclear shell energies, however, is very simple and some parts, especially its mass-number dependence, should be improved.

Some studies based on EOSs with semi-empirical expression for internal degrees of freedom \cite{sullivan16,titus18}
indicate that  nuclei above  the double-magic nuclei, $(N,Z)=(50,28)$ in the nuclear chart are primary targets for studies of electron-capture during core-collapse.
In all the EOSs discussed in this works, such nuclei  as $(N,Z)=(50,30)$ are certainly populated, while non-magic nuclei are also abundant
 in the EOSs with  more-reliable finite-temperature modeling.
Non-magic nuclei, such as $(N,Z)=(40,25)$ and  $(60,35)$, may also contribute
to the deleptonization of the core  \cite{furusawa17b}.
I conclude that further investigations of the internal partition functions of nuclei with proton numbers 25--45 and the neutron numbers 40--85 at $T\approx$0.5--3.0 MeV are needed to remove one of the serious ambiguities in the input physics for core-collapse supernova studies.
Finite-temperature effects may also change the electron-capture rates themselves in addition to the nuclear abundances  \cite{langanke03,dzhioev10,furusawa17b,titus18}.

It is surprising that  the choice of   finite-temperature modeling  produces significant
 differences in  the entropy and nuclear composition, even at the beginning of core  collapse.
 To perform  supernova simulations consistent with stellar evolution calculations, the same internal partition functions should be utilized in both stages.
Even in the postprocessing nucleosynthesis calculations for the ejecta of  supernova explosions,
the thermodynamic conditions at low temperatures and densities should be shared with the dynamical simulations.
The setup of mass data and internal degrees of freedom also needs to  be better unified for  the stages immediately following  supernova explosions.

Note that this comparative study is incomplete to cover all nuclear uncertainties.
For instance, the bulk free energy of uniform nuclear matter depends not only
on bulk properties at nuclear saturations but also on the theoretical approach such as relativistic mean-field theory, variational method, or chiral effective theory \cite{holt17}.
In addition, there are many other works for the nuclear physics inputs, such as
a comprehensive study of nuclear level densities based on two backshifted Fermi gas models and a constant-temperature model  \cite{egidy05}.
 Very recently, Raduta and Gullemineli \cite{raduta18} have constructed EOS data for supernova simulations based on
 Gullemineli and Raduta \cite{gulminelli15} with the level-density data.
Further comparison of the EOSs  with a focus on the level densities would be interesting.

In this work, I have  systematically compared several ingredients as independent inputs for EOS models.
In real, the mass data, partition functions, and nuclear matter calculations
should  be related to each other.
For instance, Rauscher's partition function is calculated by using the nucleon separation energies in the old version of the FRDM mass data, and
the WS4 mass data is based on a nuclear model with individual bulk properties, such as
$S_0$=30.16 MeV.
Level-density data consistent with HFB mass data
are also provided, which are based on Hartree-Fock-Bogoliubov theory  \cite{goriely08}.
Ideally, comparisons of self-consistent EOSs based on a specific nuclear model
should be performed.
At present, however, first-principles calculations of heavy nuclei have not been completed, even for ground-state nuclei.
More realistic supernova simulations, and  predictions that can be compared with  neutrino and gravitational wave observations,  will require
step-by-step improvements in the calculations of nuclear matter, nuclear masses, internal degrees of freedom and weak-interaction rates as nuclear physics inputs.

\begin{acknowledgments}
 The author acknowledges M. Hempel, I. N. Mishustin, K. Sumiyoshi, S. Yamada, H. Suzuki, S. Typel, G. Martinez-Pinedo, H. Sotani, Y. Yamamoto, H. Nagakura, H. Togashi, Z. Niu, and C. Kato for fruitful discussion.
The author is grateful to the anonymous referee for his or her careful
reading of the manuscript and helpful comments.
This work was supported by JSPS KAKENHI (Grant No. JP17H06365)
and  HPCI Strategic Program of Japanese MEXT  (Project ID: hp170304, 180111).
A part of the numerical calculations was carried out on  PC cluster at Center
for Computational Astrophysics, National Astronomical Observatory of Japan and XC40 at YITP in Kyoto University.
\end{acknowledgments}