\section{Introduction}\label{sec:intro}
The aim of this paper is to study the Dirichlet problem for Hessian quotient equation
\begin{equation}\label{eq:pro}
S_{k,l}(D^2u):=\frac{\sigma_k(\lambda(D^2u))}{\sigma_l(\lambda(D^2u))}=g(x)
\end{equation}
in the exterior domain $\mathbb{R}^n\setminus\overline{D}$, where $D$ is a bounded open set in $\mathbb{R}^n$ $(n\geq 3)$, $0\leq l<k\leq n$, $\lambda(D^2 u)=(\lambda_1,\cdots,\lambda_n)$ denotes the eigenvalue vector of the Hessian matrix $D^2u$,
\begin{equation*}
\sigma_j(\lambda(D^2u))=
\begin{cases}
\sum_{1\leq i_1<\cdots<i_j\leq n}\lambda_{i_1}\cdots\lambda_{i_j}, & j=1,\cdots, n,\\
1,& j=0,
\end{cases}
\end{equation*}
and $g\in C^0(\mathbb{R}^n\setminus D)$ is a positive function satisfying
\begin{equation}\label{eq:beta-g}
\inf_{\mathbb{R}^n\setminus D}g>0\quad\text{and}\quad\limsup_{|x|\to\infty}|x|^\beta|g(x)-1|<\infty
\end{equation}
for some constant $\beta>0$. We are mainly concerned with the existence and uniqueness of viscosity solutions to \eqref{eq:pro}, with prescribed boundary data and asymptotic behavior at infinity.

Note that equation \eqref{eq:pro} embraces several typical cases. When $l=0$, it is Poisson equation
\begin{equation}\label{eq:Poisson}
\sigma_1(\lambda(D^2u))=\Delta u=g
\end{equation}
if $k=1$ and fully nonlinear $k$-Hessian equation
\begin{equation}\label{eq:k-Hessian-g}
\sigma_k(\lambda(D^2u))=g
\end{equation}
if $2\leq k\leq n$, which in particular corresponds to the famous Monge--Amp\`ere equation
\begin{equation}\label{eq:M-A-g}
\det(D^2u)=g
\end{equation}
for $k=n$. It is well-known that, as a special class of nonlinear, second-order elliptic equations of the form $f(\lambda(D^2u))=g(x)$, the classical solvability of the interior Dirichlet problem for equation \eqref{eq:pro} has been extensively studied. We would like to mention the work of Caffarelli--Nirenberg--Spruck \cite{Caffarelli1984,Caffarelli1985} on equations \eqref{eq:k-Hessian-g}--\eqref{eq:M-A-g}, and that of Trudinger \cite{Trudinger1995} on general equation \eqref{eq:pro}; see also Ivochkina \cite{Ivochkina1983,Ivochkina1985}, Urbas \cite{Urbas1990}, Guan \cite{Guan1994}, and the references therein.

In exterior domains, the study of the Dirichlet problem for equation \eqref{eq:pro} also has received increasing attentions in recent years. It goes back to Meyers--Serrin \cite{Serrin1960}, who gave a sufficient condition for the source term $g$ to obtain the existence and uniqueness of solutions to \eqref{eq:Poisson} with an assigned limit at infinity. For nonlinear case of \eqref{eq:pro}, it started by Caffarelli--Li \cite{Caffarelli-Li-2003}, which is an extension of celebrated J\"{o}rgens--Calabi--Pogorelov theorem \cite{Calabi1958,Jorgens1954,Pogorelov1972} stating that any classical convex solution of
\begin{equation}\label{eq:M-A}
\det(D^2u)=1
\end{equation}
in $\mathbb{R}^n$ must be a quadratic polynomial (see also \cite{Caffarelli1995,Cheng-Yau-1986,Jost2001}). They showed that if $u$ is a convex viscosity solution of \eqref{eq:M-A} outside a bounded open convex subset of $\mathbb{R}^n$ ($n\geq3$), then there is a symmetric positive definite matrix $A\in\mathbb{R}^{n\times n}$ with $\det (A)=1$, a vector $b\in\mathbb{R}^n$ and a constant $c\in\mathbb{R}$ such that
\begin{equation}
\label{eq:C-Li}
\limsup_{|x|\to\infty}\left(|x|^{n-2}\left|u(x)-\left(\frac12x^TAx+b\cdot x+c\right)\right|\right)<\infty.
\end{equation}
Moreover, by Perron's method they established an existence and uniqueness theorem for exterior solutions of \eqref{eq:M-A} prescribed by \eqref{eq:C-Li}. Then Bao--Li--Zhang \cite{Bao-Li-Zhang-2015} extended that to general equation \eqref{eq:M-A-g}, assuming \eqref{eq:beta-g} with $\beta>2$. Recently, Li--Lu \cite{Li-Lu-2018} further obtained nonexistence results to the exterior problems considered in \cite{Bao-Li-Zhang-2015,Caffarelli-Li-2003}, giving a complete characterization for their solvability.

In spirit of \cite{Caffarelli-Li-2003}, when $g\equiv 1$ the exterior Dirichlet problem for Hessian equations \eqref{eq:k-Hessian-g} and Hessian quotient equations \eqref{eq:pro} was subsequently treated by Bao--Li--Li \cite{Bao-Li-Li-2014} and Li--Li \cite{Li-Li-2018} respectively, under proper prescribed asymptotic condition similar to \eqref{eq:C-Li}. Inspired by \cite{Bao-Li-Zhang-2015}, Cao--Bao \cite{Cao-Bao-2017} further extended the result in \cite{Bao-Li-Li-2014} to equation \eqref{eq:k-Hessian-g} with a general right hand side $g$ satisfying \eqref{eq:beta-g}. However, equation \eqref{eq:pro} with general $g\not\equiv 1$ is much more complicated to deal with than that since whenever $l>0$ elementary symmetric functions $\sigma_k$ and $\sigma_l$, of different order of homogeneities, are involved simultaneously. In particular, in this case we are confronted with the effect of irreducible fraction involved in \eqref{eq:pro}, as will be explained in detail in Remark \ref{rk:thm} (ii) below. More importantly, the comparison principle for viscosity solutions of equation \eqref{eq:pro} have not been found in a general setting of the function $g$, as far as we know, so that ones are unable to carry out the Perron process to construct the solution of the Dirichlet problem as in the Monge--Amp\`ere case and $k$-Hessian case.

Regarding more related studies on exterior domains, we refer the readers to \cite{Bao-Li-2013,Jiang-Li-Li-2020,Li-Bao-2014,Li2019} where other Hessian type equations associated with  \eqref{eq:pro}, including the special Lagrangian equations, are considered. Besides, the Liouville type results for equation \eqref{eq:pro} can be found in \cite{Bao2003,CY2010,Li-Li-Yuan-2019,Li-Ren-Wang-2016} and the references therein.

 In this paper, we shall prove the solvability of the exterior Dirichlet problem for equation \eqref{eq:pro}, under a prescribed quadratic condition of type \eqref{eq:C-Li}. Inspired by the discussions of Trudinger \cite{Trudinger1990} for the prescribed curvature equations, we investigate the comparison principles in a general setting of Hessian type fully nonlinear elliptic equations which may embrace \eqref{eq:pro} as a special case; see Appendix \ref{sec:comparison}. Then in the spirit of Crandall, Ishii, and Lions \cite{Ishii1992}, we are allowed to establish an existence and uniqueness theorem for viscosity solutions of \eqref{eq:pro} by applying an adapted Perron's method. This extends previously known results on Monge--Amp\`ere equation \eqref{eq:M-A-g} in \cite{Bao-Li-Zhang-2015} and on $k$-Hessian equation \eqref{eq:k-Hessian-g} in \cite{Cao-Bao-2017}, and also extends those obtained in \cite{Li-Li-2018} particularly for the case $g\equiv 1$ to the general $g$ fulfilling condition \eqref{eq:beta-g}. In order to state the result specifically, we introduce below some notations.

 Let $$\mathcal{A}_{k,l}=\left\{A\in S^+(n):S_{k,l}(A)=1\right\}$$ where $S^+(n)$ is the space of real $n\times n$ symmetric positive definite matrices. For $A\in S^+(n)$, we define
$$H_k(A)=\max_{1\leq i\leq n}\frac{\sigma_{k-1;i}(\lambda(A))\lambda_i(A)}{\sigma_k(\lambda(A))}$$
where $\sigma_{k-1;i}(\lambda):=\sigma_{k-1}(\lambda)|_{\lambda_i=0}$. Then $H_k(A)\in [\frac{k}{n},1)$ for $1\leq k\leq n-1$ and $H_n(A)\equiv 1$. Let $$\mathscr{A}_{k,l}=\left\{A\in\mathcal{A}_{k,l}: H_k(A)<\frac{k-l}{2}\right\}.$$ To avoid $\mathscr{A}_{k,l}$ to be empty, we require the indices $k$ and $l$ to satisfy
 \begin{equation}\label{eq:index-k-l}
\begin{cases}
0\leq l\leq n-3 &\text{if }k\geq l+2,\\
0\leq l<\frac{n}{2}-1 &\text{if }k=l+1.
\end{cases}
\end{equation}

Our main result is the following.
\begin{theorem}\label{thm:main}
Let $D$ be a smooth, bounded, strictly convex domain in $\mathbb{R}^n$ ($n\geq 3$) and let $\varphi\in C^2(\partial D)$. Assume that $g\in C^0(\mathbb{R}^n\setminus D)$ satisfies \eqref{eq:beta-g} with $\beta>2$. Then for any $A\in\mathscr{A}_{k,l}$ and $b\in\mathbb{R}^n$, there exists a constant $c_*$ depending only on $n,A,b,D,g$ and $\|\varphi\|_{C^2(\partial D)}$, such that for every $c>c_*$ there exists a unique viscosity solution $u\in C^0(\mathbb{R}^n\setminus D)$ to the problem
\begin{equation}\label{eq:asym-thm}
\begin{cases}
S_{k,l}(D^2u)=g \quad\text{in}\ \mathbb{R}^n\setminus\overline{D},\\
u=\varphi \quad\text{on}\ \partial D,\\
\limsup_{|x|\to\infty}\left(E(x)\Big{|}u(x)-(\frac12x^TAx+b\cdot x+c)\Big{|}\right)<\infty,
\end{cases}
\end{equation}
where
\begin{equation}\label{eq:intro-E}
E(x):=\left\{
\begin{array}{ll}
|x|^{\min\{\beta,\frac{k-l}{H_k(A)}\}-2} & \text{if }\beta\neq\frac{k-l}{H_k(A)},\\
|x|^{\frac{k-l}{H_k(A)}-2}(\ln |x|)^{-1} & \text{if }\beta=\frac{k-l}{H_k(A)}.
\end{array}
\right.
\end{equation}
\end{theorem}

\begin{remark}\label{rk:thm}
(i) Notice that when $l=0$ and $2\leq k\leq n$, $\mathscr{A}_{k,0}=\mathcal{A}_{k,0}$ by Proposition \ref{G-pro:A-k-l} and Theorem \ref{thm:main} was proved in \cite{Bao-Li-Zhang-2015} for equation \eqref{eq:M-A-g} and in \cite{Cao-Bao-2017} for equation \eqref{eq:k-Hessian-g}, respectively. We here extend to the general case \eqref{eq:pro} and are able to prove Theorem \ref{thm:main} in a systematic way for all of $(l,k)$ satisfying \eqref{eq:index-k-l}.

(ii) When $l\geq1$ and $g\equiv 1$, Theorem \ref{thm:main} was proved in \cite{Li-Li-2018} by rewriting \eqref{eq:pro} as $$\sigma_k(\lambda(D^2u))=\sigma_l(\lambda(D^2u)),$$ namely, not considering it in a sense of fraction. However, whenever $g\not\equiv 1$ one has to tackle the difficulty caused by the irreducible fraction in \eqref{eq:pro} when seeking subsolutions and supersolutions of \eqref{eq:pro} for carrying out the Perron process. Indeed, different from that adopted in \cite{Li-Li-2018}, here we turn to solve a new ODE to construct the desired subsolutions; see Remark \ref{G-rk:no-h-l} in Subsection \ref{G-sec:sub} for more details. Also, we can obtain the desired supersolutions in a parallel way, since unlike \cite{Li-Li-2018} we are unable to directly pick quadratic polynomials $\frac12x^TAx+b\cdot x+c$ as such ones unless assuming $\inf_{\mathbb{R}^n\setminus D}g=1$.
\end{remark}
\begin{remark}
It is not difficult to observe that Theorem \ref{thm:main} still holds with $\mathscr{A}_{k,l}$ adapted to a slightly larger set $
\mathscr{A}^*_{k,l}:=\{A\in\mathbb{R}^{n\times n}: A^*\in\mathscr{A}_{k,l}\}$
where $A^*=\frac{A+A^T}{2}$.
\end{remark}

The proof of Theorem \ref{thm:main} is based on Perron's method, formulated in Theorem \ref{thm:Perron-m} for the general fully nonlinear elliptic equation of form \eqref{app-eq:F-x-D2u} which includes \eqref{eq:pro} as a specific case. We remark that Theorem \ref{thm:Perron-m}  is an adaption of arguments as in \cite{Ishii1992,Ishii1989,Ishii-Lions-1990} once the uniqueness and comparison results are well established (see Theorem \ref{thm:comparison-b} and Corollary \ref{thm:comparison-ub}). With the help of that, the key ingredient of our proof lies in finding a family of appropriate subsolutions and supersolutions of \eqref{eq:pro} both with uniformly quadratic asymptotics at infinity. As a realization, we work with the so-called generalized symmetric functions (an extension of radial functions, specified in Definition \ref{G-def:G-Sym}), on which $k$-Hessian operators $S_{k,0}$ can be computed explicitly (see Lemma \ref{G-lem:k-He-G-Sym} below), thus helping us obtain the desired subsolutions and supersolutions by solving two corresponding second-order ODEs (\eqref{G-eq:ode-sub} and \eqref{G-eq:ode-super}). We would like to point out that the reason why we restrict $\beta>2$ and $A\in\mathscr{A}_{k,l}$ in Theorem \ref{thm:main} is to ensure the asymptotically quadratic property of subsolutions and supersolutions of \eqref{eq:pro} to be constructed in Propositions \ref{G-pro:sub} and \ref{G-pro:super}. One can see this from the computations performed in part \textbf{(c)} in Subsection \ref{G-sec:sub}.

The remainder of this paper is organized as follows. Section \ref{sec:G} is devoted to the construction of a family of subsolutions and supersolutions of \eqref{eq:pro} with asymptotically quadratic property. Then we apply Theorem \ref{thm:Perron-m} to prove Theorem \ref{thm:main} in Section \ref{sec:3}. In Appendix \ref{sec:comparison}, in a general setting of fully nonlinear elliptic equations related to the eigenvalues of the Hessian, we employ approximation arguments to prove comparison principles for viscosity solutions of the equations, both on bounded and unbounded domains. In Appendix \ref{sec:perron}, we state precisely Perron's method in Theorem \ref{thm:Perron-m} and present the proof, for the reader's convenience.

\section{Generalized symmetric subsolutions and supersolutions}\label{sec:G}
In this section, we first introduce the concept of generalized symmetric functions as well as their fine properties. Then by solving two second-order ODEs we obtain a family of generalized symmetric subsolutions and supersolutions of \eqref{eq:pro}, both of which are uniformly $k$-convex and asymptotically quadratic near infinity; see Propositions \ref{G-pro:sub} and \ref{G-pro:super}. This is essential to apply the Perron's arguments to prove Theorem \ref{thm:main} in the next section.

Throughout the section, we use the symbol $$\lambda(A):=(\lambda_1(A),\lambda_2(A),\cdots,\lambda_n(A))$$ with the order $\lambda_1(A)\leq\lambda_2(A)\leq\cdots\leq\lambda_n(A)$ to denote the eigenvalue vector of a real $n\times n$ symmetric matrix $A$.

\subsection{Preliminary}\label{G-sec:pre}
We present some definitions and related facts which will be used later to seek proper subsolutions and supersolutions of  \eqref{eq:pro}.

As in \cite{Bao-Li-Li-2014}, we define the generalized symmetric (abbreviated to G-Sym in the sequel) function in the following sense.
\begin{definition}\label{G-def:G-Sym}
For a symmetric matrix $A$, we call $u$ a G-Sym function with respect to $A$ if it is a function of $s=\frac12x^TAx$, $x\in\mathbb{R}^n$, that is $$u(x)=u(s):=u(\frac12x^TAx).$$
If $u$ is a subsolution (supersolution) of \eqref{eq:pro} and is also a G-Sym function, we say that $u$ is a G-Sym subsolution (supersolution) of \eqref{eq:pro}.
\end{definition}

There is an explicit formula \eqref{G-eq:k-H-G-Sym} below for $k$-Hessian operators acting on smooth G-Sym functions and its proof can be found in Bao--Li--Li \cite{Bao-Li-Li-2014}.
\begin{lemma}\label{G-lem:k-He-G-Sym}
Let $A=\mathrm{diag}(a_1,\cdots,a_n)$. If $w\in C^2(\mathbb{R}^n)$ is a G-Sym function with respect to $A$, then
\begin{equation}\label{G-eq:k-H-G-Sym}
\sigma_k(\lambda(D^2w))=\sigma_k(\lambda(A))(w')^k+w''(w')^{k-1}\sum_{i=1}^n\sigma_{k-1;i}(\lambda(A))(a_ix_i)^2,
\end{equation}
where $w':=\frac{dw}{ds}$ and $w'':=\frac{d^2w}{ds^2}$, $s=\frac12x^TAx$.
\end{lemma}

Taking advantage of the following properties of the $k$-th elementary symmetric function $\sigma_k$:
\begin{equation}\label{G-eq:sigma-k}
\sigma_k(\lambda)=\frac{1}{k}\sum_{i=1}^n\lambda_i\sigma_{k-1;i}(\lambda)=\sigma_{k;i}(\lambda)+\lambda_i\sigma_{k-1;i}(\lambda), \quad\forall\,
1\leq i\leq n,\forall\,\lambda\in\mathbb{R}^n,\\
\end{equation}
we can easily verify that
$$H_k(A)=\frac{\sigma_{k-1;n}(\lambda(A))\lambda_n(A)}{\sigma_k(\lambda(A))}$$ and
\begin{equation*
\frac{k}{n}\leq H_k(A)<1\text{ for }1\leq k\leq n-1; \ H_n(A)\equiv 1.
\end{equation*}

We close this preliminary by presenting the relation of $\mathcal{A}_{k,l}$ and $\mathscr{A}_{k,l}$ as follows, which can be directly verified.
\begin{proposition}\label{G-pro:A-k-l}
Suppose $0\leq l<k\leq n$ and $n\geq 3$.
\begin{enumerate}
\item[(i)] $c^*(k,l)I\in\mathcal{A}_{k,l}$, where
\begin{equation*
c^*(k,l):=\left(\frac{C_n^l}{C_n^k}\right)^{\frac{1}{k-l}},\quad C_n^j=\frac{n!}{(n-j)!j!},\  j=k,l.
\end{equation*}
\item[(ii)] If $k-l\geq 3$ or $k-l=2$ with $k<n$, then $\mathcal{A}_{k,l}=\mathscr{A}_{k,l}$.
\end{enumerate}
\end{proposition}

\subsection{Generalized symmetric subsolutions}\label{G-sec:sub}
We construct here the wanted G-Sym subsolutions of \eqref{eq:pro}.

Throughout this subsection, let $A =\text{diag}(a_1, a_2,\cdots, a_n)\in\mathcal{A}_{k,l}$, $a:=\lambda(A)$ and $w:=w(s)$ is a G-Sym function with respect to $A$, where $s=\frac12x^TAx=\frac12\sum_{i=1}^na_ix_i^2$, $x\in\mathbb{R}^n$. We denote
\begin{equation}\label{G-eq:D-s}
D(s):=\left\{x\in\mathbb{R}^n:\frac12x^TAx<s\right\}\quad\text{for }s>0.
\end{equation}

Since $g$ satisfies \eqref{eq:beta-g}, there exist $C_0$ and $s_0>1$ such that
\begin{equation}\label{G-eq:g-upper}
g(x)\leq \bar{g}(x)=\bar{g}(s):=1+C_0s^{-\frac{\beta}{2}}, \text{ when }s=\frac12 x^TAx\geq s_0.
\end{equation}
In order to make $w$ be a smooth subsolution of \eqref{eq:pro}, i.e. $S_{k,l}(D^2w)\geq g$, we consider the following ODE:
\begin{equation}\label{G-eq:ode-sub}
\begin{cases}
(w')^{k-l}+2H_k(A)w''(w')^{k-l-1}s=\bar{g}(s),& s>1,\\
w'(s)>0,\quad w''(s)\leq 0,&s\geq1.
\end{cases}
\end{equation}
 A uniformly $k$-convex (see \eqref{app-eq:k-convex}) solution of \eqref{G-eq:ode-sub} will be a G-Sym subsolution of equation \eqref{eq:pro} with respect to $A$ when $s\geq s_0$. Indeed, applying Lemma \ref{G-lem:k-He-G-Sym}, one has
\begin{align}
S_{k,l}(D^2w)&=\frac{\sigma_k(a)(w')^k+w''(w')^{k-1}\sum_{i=1}^n\sigma_{k-1;i}(a)(a_ix_i)^2}{\sigma_l(a)(w')^l+w''(w')^{l-1}\sum_{i=1}^n\sigma_{l-1;i}(a)(a_ix_i)^2}\notag\\
&\geq \frac{\sigma_k(a)(w')^k+w''(w')^{k-1}\sum_{i=1}^n\sigma_{k-1;i}(a)(a_ix_i)^2}{\sigma_l(a)(w')^l}\notag\\
&=(w')^{k-l}+w''(w')^{k-l-1}\sum_{i=1}^n\frac{\sigma_{k-1;i}(a)}{\sigma_l(a)}(a_ix_i)^2\notag\\
&\geq (w')^{k-l}+2H_k(A)w''(w')^{k-l-1}s=\bar{g}(s).\label{G-eq:S-D2w-sub}
\end{align}

For ODE \eqref{G-eq:ode-sub}, by the variation-of-constant formula we figure out it has a family of smooth solutions
\begin{equation}\label{G-eq:ode-sub-solu}
w_{c_1,c_2}(s):=c_2+\int_{s_0}^s\left(\eta^{-\mathcal{H}}\left(\int_1^\eta\mathcal{H}t^{\mathcal{H}-1}\bar{g}(t)\,dt+c_1\right)\right)^{\frac{1}{k-l}}d\eta
\end{equation}
where $c_i\in\mathbb{R}$ ($i=1,2$) are two constants and $$\mathcal{H}=\mathcal{H}(k,l,A):=\frac{k-l}{2H_k(A)}.$$
We next perform a series of careful calculations, which contains three parts, to see whether  $w_{c_1,c_2}$ given by \eqref{G-eq:ode-sub-solu} have the fine properties we expected.

\textbf{(a)} We claim $w_{c_1,c_2}'>0$ and $w_{c_1,c_2}''\leq 0$ on $[1,+\infty)$.

After a direct calculation, we have
$$w_{c_1,c_2}'(s)=\left(s^{-\mathcal{H}}\left(\int_1^s\mathcal{H}t^{\mathcal{H}-1}\bar{g}(t)\,dt+c_1\right)\right)^{\frac{1}{k-l}}>0$$
if $c_1>0$, and
 \begin{align*}
w_{c_1,c_2}''(s)&=-\frac{s^{-\mathcal{H}-1}}{2H_k(A)}(w')^{1-k+l}\left(\int_1^s\mathcal{H}t^{\mathcal{H}-1}\bar{g}(t)\,dt+c_1-s^{\mathcal{H}}\bar{g}(s)\right)\\
&=:-\frac{s^{-\mathcal{H}-1}}{2H_k(A)}(w')^{1-k+l}G(s),
\end{align*}
where
\begin{equation*}
G(s)=
\begin{cases}
c_1+\frac{C_0\beta}{2(\mathcal{H}-\frac{\beta}{2})}s^{\mathcal{H}-\frac{\beta}{2}}-1-\frac{C_0\mathcal{H}}{\mathcal{H}-\frac{\beta}{2}}&\text{if }\mathcal{H}\neq \frac{\beta}{2},\\
c_1+C_0\mathcal{H}\ln s-1-C_0&\text{if }\mathcal{H}=\frac{\beta}{2}.
\end{cases}
\end{equation*}
Hence, there is $\tilde{C}>0$, dependent on $\beta$ but independent of $s$, such that $w_{c_1,c_2}'(s)>0$ and $w_{c_1,c_2}''(s)<0$ for $s\geq1$ when $c_1>\tilde{C}$.

\textbf{(b)} Claim: $w_{c_1,c_2}(s)$ are uniformly $k$-convex on $[1,+\infty)$, provided $c_1>\tilde{C}$.

For any $1\leq m\leq k$ and $c_1>\tilde{C}$,
\begin{align*}
&\sigma_m(\lambda(D^2 w_{c_1,c_2}))\\&=\sigma_m(a)(w_{c_1,c_2}')^m+w_{c_1,c_2}'' (w_{c_1,c_2}')^{m-1}\sum_{i=1}^n\sigma_{m-1;i}(a)(a_ix_i)^2\\
&=\sigma_m(a)(w_{c_1,c_2}')^{m-1}\left(w_{c_1,c_2}'+w_{c_1,c_2}''\sum_{i=1}^n\frac{\sigma_{m-1;i}(a)}{\sigma_m(a)}(a_ix_i)^2\right)\\
&\geq \sigma_m(a)(w_{c_1,c_2}')^{m-1}\left(w_{c_1,c_2}'+2H_m(A)w_{c_1,c_2}''s\right).
\end{align*}
Since
\begin{align*}
\sigma_m(a)\sigma_{k-1;n}(a)&=\left(\sigma_{m;n}(a)+a_n\sigma_{m-1;n}(a)\right)\sigma_{k-1;n}(a)\\
&\geq \sigma_{m-1;n}(a)\sigma_{k;n}(a)+a_n\sigma_{m-1;n}(a)\sigma_{k-1;n}(a)\\
&=\sigma_{m-1;n}(a)\sigma_{k}(a)
\end{align*}
which follows from \eqref{G-eq:sigma-k} and the Newtonian inequality (see \cite{Hardy1952}):
$$\sigma_{k+1;n}(a)\sigma_{k-1;n}(a)\leq \sigma_{k;n}^2(a),\quad 1\leq k\leq n-1.$$
So that
$$H_m(A)=\frac{\sigma_{m-1;n}(a)a_n}{\sigma_m(a)}\leq \frac{\sigma_{k-1;n}(a)a_n}{\sigma_k(a)}=H_k(A)\quad \text{for }1\le m\le k.$$
Hence,
\begin{equation}\label{G-eq:k-convex-sub}
\sigma_m(\lambda(D^2w_{c_1,c_2}))\ge \sigma_m(a)(w_{c_1,c_2}')^{m-1}(w_{c_1,c_2}'+2H_k(A)w_{c_1,c_2}''s).\end{equation}
We claim $w_{c_1,c_2}'+2H_k(A)w_{c_1,c_2}''s>0$ for any $s\geq 1$. Indeed, from  \textbf{(a)} we obtain that
\begin{align*}
\frac{H_k(A)w_{c_1,c_2}''}{w_{c_1,c_2}'}&=-\frac{G(s)}{2(w_{c_1,c_2}')^{k-l}s^{\mathcal{H}+1}}\\
&=-\frac{1}{2s}\frac{\int_1^s\mathcal{H}t^{\mathcal{H}-1}\bar{g}(t)\,dt+c_1-s^{\mathcal{H}}\bar{g}(s)}{\int_1^s\mathcal{H}t^{\mathcal{H}-1}\bar{g}(t)\,dt+c_1}\\
&\geq -\frac{1}{2s}.
\end{align*}
 Consequently, \eqref{G-eq:k-convex-sub} shows $$\sigma_m(\lambda(D^2 w_{c_1,c_2}))>0, \quad 1\leq m\leq k.$$ Namely, $w_{c_1,c_2}$ are uniformly $k$-convex when $c_1>\tilde{C}$.

\textbf{(c)} We determine the asymptotic behavior of $w_{c_1,c_2}(s)$ as $s$ tends to infinity.

We start by showing that $w_{c_1,c_2}'$ is close to $1$ at infinity in both cases $\mathcal{H}\neq\frac{\beta}{2}$ and $\mathcal{H}=\frac{\beta}{2}$. When $\mathcal{H}\neq\frac{\beta}{2}$, we have
\begin{align}
w_{c_1,c_2}'(s)-1=&\left(s^{-\mathcal{H}}\left(t^{\mathcal{H}}\Big|^s_1+\frac{2C_0\mathcal{H}}{2\mathcal{H}-\beta}t^{\mathcal{H}-\frac{\beta}{2}}\Big|^s_1+c_1\right)\right)^{\frac{1}{k-l}}-1\notag\\
=&\,O\left(s^{-\min\{\frac{\beta}{2},\mathcal{H}\}}\right), \quad \text{as }s\rightarrow\infty.\label{G-eq:O-w'-sub-1}
\end{align}
When $\mathcal{H}=\frac{\beta}{2}$,
\begin{align}
w_{c_1,c_2}'(s)-1=&\left(s^{-\mathcal{H}}\left((1+C_0t^{\frac{\beta}{2}})t^{\mathcal{H}}\Big|^s_1+\frac{\beta}{2}\ln t\Big|^s_1+c_1\right)\right)^{\frac{1}{k-l}}-1\notag\\
=&O\left(s^{-\mathcal{H}}\ln s\right), \quad \text{as }s\rightarrow\infty.\label{G-eq:O-w'-sub-2}
\end{align}

We now rewrite $w_{c_1,c_2}(s)$ as:
\begin{align*}
w_{c_1,c_2}(s)&=c_2+\int_{s_0}^sw_{c_1,c_2}'(\eta)\,d\eta\\
&=c_2+s-s_0+\int_{s_0}^s\left(w_{c_1,c_2}'(\eta)-1\right)d\eta\\
&=s+\mu(c_1,c_2)-\int_{s}^\infty\left(w_{c_1,c_2}'(\eta)-1\right)d\eta
\end{align*}
where
\begin{equation}\label{G-eq:mu-c1}
\mu(c_1,c_2)=c_2-s_0+\int_{s_0}^\infty\left(w_{c_1,c_2}'(\eta)-1\right)d\eta.
\end{equation}
By \eqref{G-eq:O-w'-sub-1} and \eqref{G-eq:O-w'-sub-2}, $\mu(c_1,c_2)<\infty$ only if $\frac{\beta}{2}>1$ and $ \mathcal{H}>1$, and in this situation,
\begin{equation}\label{G-eq:asym-sub}
\begin{cases}
w_{c_1,c_2}(s)=s+\mu(c_1,c_2)+O(s^{1-\min\{\frac{\beta}{2},\mathcal{H}\}})&\text{if }\mathcal{H}\neq\frac{\beta}{2},\\
w_{c_1,c_2}(s)=s+\mu(c_1,c_2)+O(s^{1-\mathcal{H}}\ln s)&\text{if }\mathcal{H}=\frac{\beta}{2}.
\end{cases}
\end{equation}
This indicates when $\beta>2$ and $A\in\mathscr{A}_{k,l}$, $w_{c_1,c_2}$ given by \eqref{G-eq:ode-sub-solu} is asymptotically close to a quadratic polynomial at infinity.

Based on \textbf{(a)}--\textbf{(c)}, we conclude that
\begin{lemma}\label{G-lem:ode-sub-solu}
Assume $A\in\mathscr{A}_{k,l}$ and $\beta>2$. Then there exists $\tilde{C}>1$ dependent of $A$ and $\beta$, such that when $c_1>\tilde{C}$ the G-Sym function $w_{c_1,c_2}(s)$ given by \eqref{G-eq:ode-sub-solu} is a uniformly $k$-convex solution of problem \eqref{G-eq:ode-sub} and satisfies \eqref{G-eq:asym-sub}.
\end{lemma}

As a consequence of Lemma \ref{G-lem:ode-sub-solu}, we arrive at the following result.

\begin{proposition}\label{G-pro:sub}
Assume that $A\in\mathscr{A}_{k,l}$ is diagonal and $g$ satisfies \eqref{eq:beta-g} for some $\beta>2$. Let $s_0$ be as in \eqref{G-eq:g-upper} and $\tilde{C}$ be given in Lemma \ref{G-lem:ode-sub-solu}. Then when $c_1>\tilde{C}$ the G-Sym function $w_{c_1,c_2}(x)=w_{c_1,c_2}(s)$ given by \eqref{G-eq:ode-sub-solu} is a uniformly $k$-convex subsolution of equation \eqref{eq:pro}  in $\mathbb{R}^n\setminus D(s_0)$ and fulfills
\begin{equation*}
w_{c_1,c_2}(x)=\frac12x^TAx+\mu(c_1,c_2)+O\left(|x|^{2-\min\{\beta, \frac{k-l}{H_k(A)}\}}\right)\quad\text{as }|x|\to+\infty
\end{equation*}
if $\beta\neq\frac{k-l}{H_k(A)}$, or
 \begin{equation*}
w_{c_1,c_2}(x)=\frac12x^TAx+\mu(c_1,c_2)+O\left(|x|^{2-\frac{k-l}{H_k(A)}}\ln |x|\right)\quad\text{as }|x|\to+\infty
\end{equation*}
if $\beta=\frac{k-l}{H_k(A)}$, where $\mu(c_1,c_2)$ is as in \eqref{G-eq:mu-c1}.
\end{proposition}
\begin{proof}
It follows directly by combining \eqref{G-eq:g-upper}, \eqref{G-eq:S-D2w-sub} and Lemma \ref{G-lem:ode-sub-solu}.
\end{proof}

\begin{remark}\label{G-rk:no-h-l}
 We would like to make some comparisons of Proposition \ref{G-pro:sub} with related results available in literature. For $2\leq k\leq n$ and $l=0$, the G-Sym subsolutions given by Proposition \ref{G-pro:sub} also are previously constructed in \cite{Bao-Li-Zhang-2015,Cao-Bao-2017} for equations \eqref{eq:k-Hessian-g} and \eqref{eq:M-A-g}, while for $l>0$ and $g\equiv 1$ they are different from those obtained in \cite{Li-Li-2018} for equation \eqref{eq:pro} due to the ODE \eqref{G-eq:ode-sub} satisfied by subsolutions disagrees that derived in \cite{Li-Li-2018}. Actually, if $g\equiv1$, then \eqref{eq:pro} becomes $\sigma_{k}(\lambda(D^2u))=\sigma_l(\lambda(D^2u))$, from which via \eqref{G-eq:k-H-G-Sym} one can study the following ODE
\begin{equation}\label{G-eq:ode-sub-1}
(w')^{k}+2H_k(A)w''(w')^{k-1}s=(w')^{l-l}+2h_l(A)w''(w')^{l-1}s
\end{equation}
to seek subsolutions of \eqref{eq:pro}, where $h_l(A)\in[0,\frac{l}{n}]$ is defined by $$h_l(A)=\min_{1\leq i\leq n}\frac{\sigma_{l-1;i}(\lambda(A))\lambda_i(A)}{\sigma_l(\lambda(A))}.$$ Here introducing the quantities $H_k(A)$ and $h_l(A)$ at the same time is to strike a balance between different order of homogeneities and then \eqref{G-eq:ode-sub-1} could be analyzed explicitly as in \cite{Li-Li-2018}. However, whenever $g\not\equiv1$, the same idea is not helpful for removing the effect of irreducible fraction in \eqref{eq:pro}. In order to treat this, here we drop part of the denominator $\sigma_l(\lambda(D^2w))$ as handled in \eqref{G-eq:S-D2w-sub} (equivalent to putting $h_l(A)=0$ in \eqref{G-eq:ode-sub-1}). As a result, our admissible set of $A$, $\mathscr{A}_{k,l}$, may be a subset of the one in \cite{Li-Li-2018} which is given by
$$\tilde{\mathscr{A}}_{k,l}:=\left\{A\in\mathcal{A}_{k,l}: H_k(A)-h_l(A)<\frac{k-l}{2}\right\}.$$
Nevertheless, by Proposition \ref{G-pro:A-k-l} our method still works well in cases $k-l\geq2$ where it holds $$\mathscr{A}_{k,l}=\tilde{\mathscr{A}}_{k,l}=\mathcal{A}_{k,l}.$$
\end{remark}

\subsection{Generalized symmetric supersolutions}\label{G-sec:sup}
We will construct a family of G-Sym supersolutions of \eqref{eq:pro} coinciding at infinity with the G-Sym subsolutions we just obtained in Proposition \ref{G-pro:sub}, which is necessary in Perron's construction (Theorem \ref{thm:Perron-m}). Unlike the well-treated case $g\equiv 1$ in which one can directly pick $\frac12x^TAx+c$ to be desired supersolutions (see for example \cite{Bao-Li-Li-2014,Caffarelli-Li-2003,Jiang-Li-Li-2020,Li-Li-2018,Li-Bao-2014}), here we adapt the idea in Subsection \ref{G-sec:sub} to study the ODE:
\begin{equation}\label{G-eq:ode-super}
\begin{cases}
(w')^{k-l}+2H_k(A)w''(w')^{k-l-1}s=\underline{g}(s),& s>1,\\
w'(s)>0, \quad w''(s)\geq0,&s\geq1,
\end{cases}
\end{equation}
where $\underline{g}$ is an increasing smooth function of $s$ satisfying, without loss of generality,
\begin{equation}\label{G-eq:g-lower}
0<\underline{g}\leq g\text{ for}\ s\geq 1,\quad\text{and}\quad\underline{g}=1-C_0s^{-\frac{\beta}{2}}\ \text{when }s\geq s_0,
\end{equation}
where $C_0$ and $s_0$ are as in \eqref{G-eq:g-upper}. Replacing $\bar{g}$ by $\underline{g}$ in \eqref{G-eq:ode-sub-solu} and taking $c_1=0$, one can easily obtain the following uniformly $k$-convex solutions of \eqref{G-eq:ode-super}:
\begin{equation}\label{G-eq:ode-super-solu}
\overline{w}_{c_2}(s):=c_2+\int_{1}^s\left(\eta^{-\mathcal{H}}\left(\int_1^\eta\mathcal{H}t^{\mathcal{H}-1}\underline{g}(t)\,dt\right)\right)^{\frac{1}{k-l}}d\eta
\end{equation}
Moreover, thanks to \eqref{G-eq:g-lower}, arguing as in part \textbf{(c)} before we infer that $\overline{w}_{c_2}$ possess asymptotics \eqref{G-eq:asym-sub} where $\mu(c_1,c_2)$ need to be replaced by
\begin{equation}\label{G-eq:lmu-c1-c2}
\overline{\mu}(c_2)=c_2-1+\int_{1}^\infty\left(\overline{w}_{c_2}'(\eta)-1\right)d\eta.
\end{equation}

Similarly as showed in \eqref{G-eq:S-D2w-sub}, we can find for $s\geq 1$, $$S_{k,l}(D^2\overline{w}_{c_2})\leq (\overline{w}_{c_2}')^{k-l}+2H_k(A)\overline{w}_{c_2}''(\overline{w}_{c_2}')^{k-l-1}s=\underline{g}(s)\leq g.$$
Namely, $\overline{w}_{c_2}$ are supersolutions of \eqref{eq:pro}.

The following statement is a summary of above facts, which serves as a counterpart of Proposition \ref{G-pro:sub}.
\begin{proposition}\label{G-pro:super}
Assume that $A\in\mathscr{A}_{k,l}$ is diagonal and $g$ satisfies \eqref{eq:beta-g} for some $\beta>2$. Let $c_2\in\mathbb{R}$ and $s=\frac12x^TAx$. Then the G-Sym function $\overline{w}_{c_2}(x)=\overline{w}_{c_2}(s)$ given by \eqref{G-eq:ode-super-solu} is a uniformly $k$-convex supersolution of equation \eqref{eq:pro} in $\mathbb{R}^n\setminus D(1)$ and fulfills
\begin{equation*}
\overline{w}_{c_2}(x)=\frac12x^TAx+\overline{\mu}(c_2)+O\left(|x|^{2-\min\{\beta, \frac{k-l}{H_k(A)}\}}\right)\quad\text{as }|x|\to+\infty
\end{equation*}
if $\beta\neq\frac{k-l}{H_k(A)}$, or
 \begin{equation*}
\overline{w}_{c_2}(x)=\frac12x^TAx+\overline{\mu}(c_2)+O\left(|x|^{2-\frac{k-l}{H_k(A)}}\ln |x|\right)\quad\text{as } |x|\to+\infty
\end{equation*}
if $\beta=\frac{k-l}{H_k(A)}$, where $\overline{\mu}(c_2)$ is as in \eqref{G-eq:lmu-c1-c2}.
\end{proposition}

\begin{remark}
One may wonder whether there are smooth G-Sym solutions of \eqref{eq:pro}, provided that $g$ is of G-Sym. Generally, this cannot be expected for all $0\leq l<k\leq n$. When $g\equiv 1$, \cite[Proposition 1.1]{Li-Li-Zhao-2019} states that \eqref{eq:pro} admits a G-Sym solution of $C^2$ with respect to a diagonal $A\in\mathcal{A}_{k,l}$ if and only if $l=0$ and $k=n$, unless $A=c^*(k,l)I$ (see (i) of Proposition \ref{G-pro:A-k-l}). The same rigidity result was also exploited in \cite{Cao-Bao-2017} for \eqref{eq:k-Hessian-g} with $g=\bar{g}$, introduced in \eqref{G-eq:g-upper}. That is the reason why we here look for G-Sym subsolutions and G-Sym supersolutions of \eqref{eq:pro} separately.
\end{remark}

\section{Proof of Theorem \ref{thm:main}} \label{sec:3}
In view of Theorem \ref{thm:Perron-m}, to prove Theorem \ref{thm:main} it suffices to demonstrate the existence of a viscosity subsolution $\underline{u}$ of equation \eqref{eq:pro} attaining the prescribed boundary value and asymptotic behavior at infinity, as well as that of a viscosity supersolution $\bar{u}\geq\underline{u}$ but agreeing on $\underline{u}$ at infinity. Roughly speaking, such a subsolution can be obtained by splicing together the supremum of barriers over the boundary points of the domain and the G-Sym subsolution constructed in Proposition \ref{G-pro:sub}; the desired supersolution is prepared in Proposition \ref{G-pro:super}. The uniqueness of the solution is guaranteed by comparison principle, Corollary \ref{thm:comparison-ub}.

Set $\mathcal{G}=\sup_{\mathbb{R}^n\setminus D} g$. To process the boundary behavior of the solution, we need the following existence result of barrier functions.

\begin{lemma}\label{P-lem:w-xi}
Let $D$ be a bounded strictly convex domain of $\mathbb{R}^n$ ($n\geq3$) with $\partial D\in C^2$ and let $\varphi\in C^2(\partial D)$. For an invertible and symmetric matrix $A$, there exists some constant $C$, depending only on $n,\mathcal{G},\|\varphi\|_{C^2(\partial D)}$, the upper bound of $A$, the diameter and the convexity of $D$, and the $C^2$ norm of $\partial D$, such that for every $\xi\in\partial D$, there exists $\bar{x}(\xi)\in\mathbb{R}^n$ satisfying $$|\bar{x}(\xi)|\leq C\quad\text{and}\quad w_\xi<\varphi\quad\text{on}\quad \partial D\setminus\{\xi\},$$ where $$w_\xi(x)=\varphi(\xi)+\frac{\mathcal{G}^{\frac{1}{k-l}}}{2}\left[(x-\bar{x}(\xi))^TA(x-\bar{x}(\xi))-(\xi-\bar{x}(\xi))^TA(\xi-\bar{x}(\xi))\right]$$
for $x\in\mathbb{R}^n$.
\end{lemma}

\begin{proof}
As proved in \cite{Cao-Bao-2017} for $A\in\mathcal{A}_{k,0}$, it is a direct adaption of the arguments as in the proofs of \cite[Lemma 5.1]{Caffarelli-Li-2003} and \cite[Lemma 3.1]{Bao-Li-Li-2014} for the case $\mathcal{G}=1$. We thus omit it.
\end{proof}
We now present the proof of Theorem \ref{thm:main}.

\begin{proof}[Proof of Theorem \ref{thm:main}]
By an orthogonal transformation and by subtracting a linear function from $u$, we need only to prove for $A=\text{diag}(a_1,a_2,\cdots,a_n)$ and $b=0$; see \cite[Lemma 3.3]{Li-Li-2018} for a specific demonstration. Also, without loss of generality, we assume $D(1)\subset D\subset D(s_0)$, where $D(\cdot)$ is as in \eqref{G-eq:D-s} and $s_0$ is as in \eqref{G-eq:g-upper}.

\textbf{Step 1} We first construct a viscosity subsolution $\underline{u}$ of \eqref{eq:pro} with $\underline{u}=\varphi$ on $\partial D$ and the asymptotics as in \eqref{eq:asym-thm}.

Recalling $w_{c_1,c_2}(s)$ given in \eqref{G-eq:ode-sub-solu}, by Proposition \ref{G-pro:sub} it is a uniformly $k$-convex subsolution of \eqref{eq:pro} satisfying
\begin{equation}\label{P-eq:asym-E}
w_{c_1,c_2}(x)=\frac12x^TAx+\mu(c_1,c_2)+O\left(E^{-1}(x)\right),\quad |x|\to\infty,
\end{equation}
 provided $c_1>\tilde{C}$ and $x\in\mathbb{R}^n\setminus D(s_0)$. Here $E(x)$ is defined in \eqref{eq:intro-E}. We later will pick suitable $c_1$ to make \eqref{P-eq:asym-E} reach asymptotics \eqref{eq:asym-thm} as in \eqref{P-eq:mu-c1-c} below. Now regarding the boundary value, we are going to find another viscosity subsolution $\underline{w}$ attaining $\varphi$ on $\partial D$, such that, for some fixed $\bar{s}>s_0$,
\begin{equation}\label{P-eq:barrier}
\max_{\partial D(s_0)}w_{c_1,c_2}\leq\min_{\partial D(s_0)}\underline{w}\quad\text{and}\quad\min_{\partial D(\bar{s})}w_{c_1,c_2}\geq\max_{\partial D(\bar{s})}\underline{w}.
\end{equation}
For this aim, we set $$\underline{w}(x)=\max\{w_{\xi}(x)~|~\xi \in \partial D\},$$ where $w_\xi$ is introduced in Lemma \ref{P-lem:w-xi}. Clearly, $\underline{w}=\varphi$ on $\partial D$.  Since $$S_{k,l}(D^2w_\xi)=\mathcal{G}S_{k,l}(A)=\mathcal{G}\geq g\quad\text{in }\mathbb{R}^n\setminus\overline{D},$$
$w_\xi$ is a smooth convex subsolution of \eqref{eq:pro}. By Lemma \ref{lem:Perron-sub}, $\underline{w}$ is also a viscosity subsolution of \eqref{eq:pro}. Let
$$c_2=m_{s_0}:=\min\{w_{\xi}(x)~|~\xi\in\partial D, x\in \overline{D(s_0)}\setminus D\}.$$
Thus, by \eqref{G-eq:ode-sub-solu} $w_{c_1,m_{s_0}}$ satisfies the first condition in \eqref{P-eq:barrier}. To realize the second one, it suffices to choose a large $c_1$ (assume $c_1\geq\alpha>\tilde{C}$) since $w_{c_1,m_{s_0}}$ is monotonically increasing with respect to $c_1$.

We next fix $c_*$ such that given $c>c_*$ one can find $c_1(c)$ to fulfill
\begin{equation}\label{P-eq:mu-c1-c}
\mu(c_1(c),m_{s_0})=c.
\end{equation}
Notice from \eqref{G-eq:mu-c1} that $\mu(c_1,m_{s_0})$ is strictly increasing in $c_1$, and $$\lim_{c_1\to+\infty}\mu(c_1,m_{s_0})=+\infty.$$ Hence, if $c_*\geq\mu(\alpha,m_{s_0})$, then such $c_1(c)>\alpha$ exists. This means that, for $c>c_*$, $w_{c_1(c),m_{s_0}}$ satisfies \eqref{P-eq:barrier} and possesses asymptotic behavior:
\begin{equation}\label{P-eq:asym-sub}
w_{c_1(c),m_{s_0}}(x)=\frac12x^TAx+c+O\left(E^{-1}(x)\right),\quad |x|\to\infty.
\end{equation}

For $c>c_*$, we define
\begin{equation*}
\underline{u}(x)=
\begin{cases}
\underline{w}(x),\quad &x\in D(s_0)\backslash D,\\
\max\{w_{c_1(c),m_{s_0}}(x),\underline{w}(x)\},\quad &x\in D(\bar{s})\backslash D(s_0),\\
w_{c_1(c),m_{s_0}}(x),&x\in\mathbb{R}^n\backslash D(\bar{s}).
\end{cases}
\end{equation*}
Then from Definition \ref{app-def:visc} and Lemma \ref{lem:Perron-sub}, $\underline{u}$ is a viscosity subsolution of \eqref{eq:pro} satisfying \eqref{P-eq:asym-sub}, and $\underline{u}=\underline{w}=\varphi$ on $\partial D$.

\textbf{Step 2}
We construct a viscosity supersolution $\bar{u}$ of \eqref{eq:pro} to satisfy
\begin{equation}\label{P-eq:step-super}
\underline{u}\leq\bar{u}\ \text{ in }\mathbb{R}^n\setminus\overline{D}\quad\text{and}\quad\lim_{|x|\to\infty}(\bar{u}-\underline{u})(x)=0.
\end{equation}

  By Proposition \ref{G-pro:super}, $\overline{w}_{c_2}$ given in \eqref{G-eq:ode-super-solu} is a uniformly $k$-convex supersolution of \eqref{eq:pro} satisfying
\begin{equation}\label{P-eq:asym-super}
\overline{w}_{c_2}(x)=\frac{1}{2}x^TAx+\overline{\mu}(c_2)+O\left(E^{-1}(x)\right),\quad |x|\to\infty.
\end{equation}
Then for $c>c_*$ the second condition in \eqref{P-eq:step-super} holds for $\overline{w}_{c_2(c)}$ where $c_2(c)$ is determined by $$\overline{\mu}(c_2(c))=c,$$ which implies \eqref{P-eq:asym-super} agrees with \eqref{P-eq:asym-sub}. Actually, by definition of $\overline{\mu}$ (see \eqref{G-eq:lmu-c1-c2}), $$c_2(c)=c+1-\int_1^\infty\left[\left(s^{-\mathcal{H}}\int_1^s\mathcal{H}t^{\mathcal{H}-1}\underline{g}(t)\,dt\right)^{\frac{1}{k-l}}-1\right]\,ds\geq c\,;$$
here we have used the fact that $\underline{g}\leq 1$ for $s\geq 1$. Next, with the help of comparison principles proved in Appendix \ref{sec:comparison}, we show that $\overline{w}_{c_2(c)}$ also agrees the first condition in \eqref{P-eq:step-super} for proper $c$.

In \textbf{Step 1}, we fixed $c_*\geq\mu(\alpha,m_{s_0})$. We now further require that $c_*>M_{s_0}$, where
$$M_{s_0}:=\max\{w_{\xi}(x)~|~\xi\in\partial D, x\in \overline{D(s_0)}\setminus D\}.$$
Then for $c>c_*$,
$$\overline{w}_{c_2(c)}\geq c_2(c)\geq c>c_*>M_{s_0}\geq m_{s_0}\geq w_{c_1(c),m_{s_0}}\quad\text{on }\partial D(s_0),$$
and also $$\lim_{|x|\to\infty}(\overline{w}_{c_2(c)}-w_{c_1(c),m_{s_0}})(x)=0.$$
Applying Corollary \ref{thm:comparison-ub}, we thus deduce that
\begin{equation}\label{P-eq:compare-super-sub}
\overline{w}_{c_2(c)}\geq w_{c_1(c),m_{s_0}}\quad\text{in }\mathbb{R}^n\setminus\overline{D(s_0)}.
\end{equation}
On the other hand,
$$\overline{w}_{c_2(c)}\geq c_2(c)\geq c>c_*>M_{s_0}\geq \underline{w}\quad\text{on }\partial D,$$
and by \eqref{P-eq:barrier}, $$\overline{w}_{c_2(c)}\geq w_{c_1(c),m_{s_0}}\geq\underline{w}\quad\text{on }\partial D(\bar{s}).$$
Hence, applying Theorem \ref{thm:comparison-b} yields
\begin{equation}\label{P-eq:compare-super-barrier}
\overline{w}_{c_2(c)}\geq\underline{w}\quad\text{in }D(\bar{s})\setminus\overline{D}.
\end{equation}
By virtue of \eqref{P-eq:compare-super-sub} and \eqref{P-eq:compare-super-barrier}, we get $$\overline{w}_{c_2(c)}\geq \underline{u}\quad\text{in }\mathbb{R}^n\setminus\overline{D}$$ for $c>c_*$.

Based on the above arguments, we let $\bar{u}=\overline{w}_{c_2(c)}$ for $c>c_*$, which is a desired viscosity supersolution of \eqref{eq:pro}.

\textbf{Step 3} We show the existence and uniqueness of viscosity solutions to problem \eqref{eq:asym-thm}.

With $\underline{u}$ and $\bar{u}$, we define
\begin{align*}
u(x):=&\sup\{v(x)| v\in\mathrm{USC}(\mathbb{R}^n\setminus\overline{D}), S_{k,l}(D^2v)\geq g\text{ in }\mathbb{R}^n\setminus\overline{D}\text{ in the}\\
&\text{viscosity sense, with }\underline{u}\leq v\leq\bar{u}\text{ in }\mathbb{R}^n\setminus\overline{D}\text{ and }v=\varphi\text{ on }\partial D\}.
\end{align*}
Since $\underline{u}$ and $\bar{u}$ both satisfy \eqref{P-eq:asym-sub}, one has
$$\limsup_{|x|\to\infty}\left(E(x)\Big{|}u(x)-(\frac12x^TAx+c)\Big{|}\right)<\infty.$$
From Theorem \ref{thm:Perron-m}, we thus conclude that $u\in C^0(\mathbb{R}^n\setminus D)$ is a viscosity solution of \eqref{eq:asym-thm}.

Finally, the uniqueness of viscosity solutions to problem \eqref{eq:asym-thm} is forced by Corollary \ref{thm:comparison-ub}. This completes the proof.
\end{proof}

\begin{remark}
We remark that from the above demonstration the lowerbound $c_*$ of $c$ in Theorem \ref{thm:main} may not be removed. Actually, in \cite[Theorem 1.3]{Li-Lu-2018}, Li and Lu proved that when $l=0$ and $k=n$, there is a sharp constant $c_*$ such that problem \eqref{eq:asym-thm} admits a viscosity solution if and only if $c\geq c_*$. It would be interesting to see whether
such a sharp charaterization result holds for equation \eqref{eq:pro} with general $0\leq l<k\leq n$ and will be left to study in our future work.
\end{remark}


\section{Introduction}

Recent efforts in band structure calculations have shown that it might be possible to obtain many intermetallic compounds with topologically non-trivial surface states \cite{Zhang2019,Vergniory2019,Tang2019}. Among these systems, a few become superconducting at low temperatures. It has been shown that superconductivity induced by proximity into a topologically non-trival electronic surface state can lead to triplet correlations and the appearance of Majorana modes at the surface \cite{PhysRevB.84.195442,doi:10.1146/annurev-conmatphys-030212-184337,Alicea_2012,PhysRevLett.100.096407}.

Here we analyze the compound Au$_2$Pb. Au$_2$Pb is a binary intermetallic system which crystallizes at room temperature in a cubic Laves phase \cite{Schoop2015,PhysRevB.93.045118}. At low temperatures, Au$_2$Pb presents an orthorhombic phase which is superconducting with $T_c=1.1$ K \cite{HAMILTON1965665}. Band structure calculations of the room temperature cubic phase predict a bulk 3D Dirac cone due to symmetry allowed band crossings \cite{Schoop2015,PhysRevB.93.045118}. The presence of the bulk Dirac cone is confirmed by angle resolved photoemission studies (ARPES) \cite{Wu2018}. The lowered crystal symmetry of the low temperature phase lifts the band degeneracy and a gap opens at the bulk Dirac cone \cite{Schoop2015,Wu2018}. As we show in detail below and represent schematically in Fig.\,\ref{Schematics}(a), the parity of the bands is inverted at the gap, with the valence band having Au d-electron character and the conduction band Pb-sp$^3$ character.

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{FigScheme_1.pdf}
\caption{(a) Schematic representation of the band inversion occuring in Au$_2$Pb upon cooling into the orthorhombic low temperature phase \cite{Schoop2015,Wu2018}. The color code represents the relation between Pb-sp$^3$ and Au-d character from blue to yellow. At the surface, it is possible that the inverted gap closes (dashed black lines). (b) Schematic representation of our findings. The bulk of Au$_2$Pb is shown in blue and the surface in grey. We find a bulk superconducting gap with $\Delta=0.17$ meV and a surface with a rugosity (white to black spots on the top, following scale bar on the left) which is of the order of a unit cell size and consists of a random height variation. We also find that the surface presents a superconducting density of states with a finite value at zero energy. In (c) we show schematically the superconducting density of states of a system with a surface Dirac cone dispersion inside the superconducting gap (blue) and the localized states generated by scattering (dark grey) at a certain energy $E_{Dist}$ in the vicinity of the Fermi level $E_F$. In a system with a large coherence length and a random distribution of localized states at many different $E_{Dist}$, these smear the low energy part of the density of states into a gapless regime\cite{doi:10.1080/00018732.2014.927109}.
}
\label{Schematics}
\end{figure}

Experiments aiming to characterize surface states in the normal phase using optical conductivity have concluded that the contribution of the Dirac cone to the optical properties is masked by bulk bands crossing the Fermi energy \cite{Schoop2015,Wu2018,Kemmler_2018}. In the superconducting phase, low temperature specific heat and thermal conductivity measurements have suggested that the bulk superconducting state of Au$_2$Pb is fully gapped, agreeing with s-wave superconductivity \cite{Schoop2015, Yu2016}. Here we analyze the superconducting properties using Scanning Tunnelling Microscope (STM) experiments on the surface of Au$_2$Pb. We measure the superconducting density of states as a function of temperature and magnetic field. We find a superconducting gap value corresponding to a $T_c=1.1$ K within BCS theory and rough surfaces (Fig.\,\ref{Schematics}(b)). The zero bias tunneling conductance is larger than about a third of the normal state density of states. We analyze the possible origin for the finite zero bias tunneling conductance in terms of the influence of defects on surface superconducting states, creating in-gap states schematically represented Fig.\,\ref{Schematics}(c).

\section{Experiments and methods}

\begin{figure}[htbp]
\includegraphics[width =0.8\columnwidth]{figStructureCp.pdf}
\caption{(a) Powder x-ray data at room temperature of Au$_2$Pb (purple lines). We show the expected Bragg peaks by vertical lines on the bottom. A picture of the crystals is shown in the upper right inset. We also show a picture of atomic arrangements (Au is shown in gold color and Pb in grey) of the room temperature cubic phase. (b) Resistivity vs temperature of Au$_2$Pb is shown as colored points. The specific heat vs temperature is shown in the inset (with a zoom for the transition at 50 K as an additional inset). We use different colors for the different structural phases. The low temperature crystalline structure (points in blue) is shown in the upper left inset.}
\label{FigSample}
\end{figure}

We synthesized the sample with a Pb flux growth method \cite{doi:10.1080/13642819208215073,Canfield_2019,Schoop2015,Wu2018}. Lumps of gold (99.99\% Sempsa) and lead (99.999\% Goodfellow) with a molar ratio Au:Pb 42:58 were put inside frit disk alumina crucibles \cite{doi:10.1080/14786435.2015.1122248} and sealed in an evacuated silica ampoule. The ampoule was then heated in 12 hours to 1150$^{\circ}$C, fast-cooled to 375$^{\circ}$C in 12 hours, held at this temperature for 6 hours and then slow-cooled to 275$^{\circ}$C in 22 hours. At this point the excess flux was decanted using a centrifuge. The obtained crystals typically had sizes of a few millimeters and showed clean triangular dark-silvered facets. We show powder x-ray data in Fig.\,\ref{FigSample}(a) and a picture of the samples in the upper left inset. In Fig.\,\ref{FigSample}(b) we show the resistivity as a function of temperature, and highlight with colors the different crystalline structures\cite{PhysRevB.93.045118,Schoop2015}. In the lower right inset of Fig.\,\ref{FigSample}(b) we show the specific heat vs temperature. These data mostly coincide with previous results in high quality single crystals of Au$_2$Pb\cite{Schoop2015,Wu2018}. Furthermore, in the insets of Fig.\,\ref{FigSample}(a,b) we show the room and low temperature crystalline structures found in Au$_2$Pb. At room temperature (violet points in Fig.\,\ref{FigSample}(b)), the crystal structure is the cubic Laves phase $Fd\overline{3}m$ (no. 227) \cite{PhysRevB.93.045118}. It is based on the pyrochlore lattice, with Au atoms having a tetrahedric arrangement distributed inside a face centered cubic lattice of Pb atoms. The crystalline structures of the phase at intermediate temperatures (yellow and red points in Fig.\,\ref{FigSample}(b)) are unknown. The structure of the low temperature phase (blue points in Fig.\,\ref{FigSample}(b)) is derived from the cubic room temperature phase by small atomic distortions, giving orthorhombic $Pbcn$, No. 60 \cite{Schoop2015}.

We performed STM measurements using a home-made microscope described in Ref.\,\onlinecite{doi:10.1063/1.3567008} and the software described in Ref.\,\cite{doi:10.1063/5.0064511}. The energy resolution of the tunneling conductance measurement setup is of 10 $\mu$eV as shown measurements in Al with a similar set-up ($T_c=1.1$ K) \cite{doi:10.1063/5.0059394}. For image treatment, we use Refs.\,\cite{doi:10.1063/5.0064511,doi:10.1063/1.2432410}. We cut single crystalline samples in rectangular bars that are a few mm long and a fraction of a mm wide and mount them on the STM sample holder. The position of the sample holder can be modified to approach the long axis of the bar shaped sample to a blade in-situ at 4 K in such a way that the blade cleaves the sample and exposes a fresh surface. We exposed a surface perpendicular to a main crystalline axis of the cubic room temperature crystalline structure. Given the different structural transitions undergone upon cooling we can expect that we maintain a crystalline face at low temperatures. We observed large and flat areas over fields of view of a few hundreds of nm. We did not observe atomic resolution, but instead surfaces with elongated structures and a surface corrugation of order of the unit cell size (between 1 nm and 2 nm, see Fig.\,\ref{FigHisto}(a) and the Appendix).

First-principles calculations with spin-orbit coupling (SOC) were performed using the Quantum ESPRESSO package \cite{Giannozzi_2009}. In our calculations we employed the generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) \cite{PhysRevLett.77.3865} exchange-correlation functional. We used fully-relativistic norm-conserving Pseudo DoJo pseudopotentials \cite{VANSETTEN201839} and the electronic wave functions were expanded with well-converged kinetic energy cutoffs of 45 Ry and 180 Ry for the wavefunctions and charge density, respectively. The Brillouin zone was sampled with a Monkhorst-Pack \cite{PhysRevB.13.5188} ($8\times 8 \times 8$) $k$-point mesh. Dispersion interactions to account for van der Waals interactions between the layers were considered by applying semi-empirical Grimme DFT-D3 corrections \cite{https://doi.org/10.1002/jcc.20495}. Structures were fully optimized using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm \cite{doi:10.1063/1.462844} until the forces on each atom were smaller than $2\cdot 10^4$ Ry/au and the energy difference between two consecutive relaxation steps less than $10^4$ Ry. To analyze the surface bandstructure, we constructed a tight-binding model based on maximally localized Wannier functions \cite{PhysRevB.56.12847,PhysRevB.65.035109,RevModPhys.84.1419} with Au $sd$ and Pb $sp$ orbitals in the range $E_F \pm 1$ eV. We verified that we find the same results as in the first principles calculations within this energy range. To find the structure of surface states, we built surface spectral functions for a semi-infinite Au$_2$Pb(001) surface via the surface Green’s function method \cite{PhysRevB.23.4988,PhysRevB.23.4997,Sancho_1984,Sancho_1985} as implemented in the WannierTools package \cite{WU2018405}.

\section{Results}

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{figZeroFieldT.pdf}
\caption{(a) Tunneling conductance as a function of the bias voltage for different temperatures (colored lines, bottom is 100 mK, in light blue, and top is 1 K, in red, with 100 mK increases). The tunneling conductance is normalized to the value at bias voltages between 0.5 mV and 1 mV. (b) Local density of states (LDOS) as a function of energy for the same temperatures and using the same color scheme. Black lines in (a) provide the LDOS shown in (b) convoluted with temperature. (c) Temperature dependence of the superconducting gap, also shown by the arrows in (b) The BCS temperature dependence of the superconducting gap is shown by a black line.
}
\label{FigZeroField}
\end{figure}

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{FigHisto.pdf}
\caption{(a) Scanning Tunneling Microscopy image obtained with a constant current of 1 nA and at a bias voltage of 10 mV. The obtained height changes are shown as a grey scale, following the bar on the left. The height profile along the yellow line is shown in the bottom panel. (b) Tunneling conductance map at zero bias in the same field of view as in (a). Conductance changes are shown as a blue-red variation following the bar on the left. The conductance profile along the yellow line is shown in the bottom panel. White scale bars in (a,b) are 60 nm long. (c) Tunneling spectroscopy as a function of the bias voltage along the yellow lines in (a,b).
}
\label{FigHisto}
\end{figure}

\begin{figure*}[htbp]
\includegraphics[width = \textwidth]{figTeo1.pdf}
\caption{(a) The electronic bandstructure of the low temperature orthorhombic phase of Au$_2$Pb, obtained from density functional calculations is shown by the black lines. In the upper right inset we show as lines the partial densities of states of the orbitals marked in the legend as a function of the energy. We see that, at the Fermi level, the largest contributions are due to Pb (p), Au (s) and Au (d) orbitals. In the lower left inset we show a zoom over the $\Gamma$-$X$ direction, with the gap that opens in the low temperature phase and the associated band inversion \cite{Schoop2015,Wu2018}. The blue to yellow color scale shows the relation between Pb-sp$_3$ and Au d-orbital character, following the scale bar at the right. (b) Fermi surface of the orthorhombic low temperature phase of Au$_2$Pb is shown as colored surfaces. The first Brillouin zone, with the high symmetry directions, is shown in red. The Fermi surface is made by three bands, the 15$^{th}$ (in blue, (c)), the 17th$^{th}$ (in yellow, (d)) and the 19th$^{th}$ band (in red, (e)).}
\label{FigTeo1}
\end{figure*}

In Fig.\,\ref{FigZeroField} we show Scanning Tunneling Spectroscopy measurements of the superconducting gap of Au$_2$Pb as a function of temperature. We see that there is no fully opened superconducting gap, but instead a large zero bias conductance. The quasiparticle peaks are also supressed with respect to the expected high peaks for a usual s-wave BCS superconductor. We deconvolute the local density of states ($N(E)$) at each temperature, using the relation $G(V)=\int N(E)\frac{df(E-eV)}{dV} dE$ (see Refs.\,\cite{PhysRevLett.96.027003,PhysRevLett.101.166407,PhysRevB.97.134501,doi:10.1063/5.0064511} and the Appendix) and show the result in Fig.\,\ref{FigZeroField}(b). From the density of states, we obtain the superconducting gap $\Delta$ and its temperature dependence by following the quasiparticle peak positions in $N(E)$. We find a superconducting gap of $\Delta(T=0K)=0.17$ meV. The BCS expression $\Delta(T=0$ K$)=1.76k_BT_c$ provides $\Delta(T=0$ K$)\approx 0.17$meV using $T_c=1.1$ K, in close agreement with the value we find. Fig.\,\ref{FigZeroField}(c) shows the temperature dependence of the superconducting gap $\Delta(T)$ up to $T_c=1.1$ K. $\Delta(T)$ roughly follows the temperature dependence of a fully opened s-wave superconducting gap expected within BCS theory.

The variations of the zero bias conductance as a function of the position are shown in Fig.\,\ref{FigHisto}(b). We observe variations that occur randomly over the whole surface. These are small, below about 10\% of the normal state tunneling conductance. The magnetic field dependence is shown in the Appendix and is compatible with superconductivity having a low Ginzburg-Landau parameter $\kappa$ and a coherence length of about $\xi\approx$ 60 nm, as shown previously\cite{Schoop2015,Wu2018}. Thus, we find a gap value which is compatible with specific heat measurements \cite{Schoop2015,Xing2016,Wu2018}, but there is a sizeable finite density of states at the Fermi level which slightly varies as a function of the position.

We show the calculated electronic band structure of the low temperature orthorhombic phase in Fig.\,\ref{FigTeo1} (a). The overall shape coincides with previous calculations \cite{Schoop2015}. The Dirac cone of the room temperature cubic structure is opened by a gap along the $\Gamma-X$ line. The Fermi surface presents three bands crossings, the 15$^{th}$, the 17$^{th}$ and the 19$^{th}$ bands. All are spin degenerate and produce three dimensional pockets with an intricate structure (Fig.\,\ref{FigTeo1}(b)). The top and bottom of the bands crossing the Fermi level are just a few hundreds of meV from the Fermi level (Fig.\,\ref{FigTeo1}(a)), confirming that Au$_2$Pb is a low density semimetallic system. The 15$^{th}$ band (Fig.\,\ref{FigTeo1}(c)) consists essentially of two hole pockets located at the top side and the edges of the Brillouin zone. It has a mixed Pb-p and Au-s,d character. The 17$^{th}$ band (Fig.\,\ref{FigTeo1}(d)) consists of three electron pockets. One closed spherical pocket located well inside the Brillouin zone. Another one goes from the top to the bottom of the Brillouin zone, and the third one is located at the sides of the Brillouin zone. The 19$^{th}$ band (Fig.\,\ref{FigTeo1}(e)) produces two electron pockets that are nested inside the 17$^{th}$ band pockets located at the sides of the Brillouin zone. Actually, these two pockets are nearly degenerate, as shown by the bandstructure between the $T$ and $Y$ high symmetry points and both have a predominant Au-s,d character.

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{figTeo2.pdf}
\caption{Surface bandstructure of the low temperature orthorhombic phase of Au$_2$Pb, with the electronic density shown in blue to red. In (a) we show the result for the (010) surface as a function of energy. In blue are gapped regions without bulk or surface states. The green rectangle shows the region corresponding to the green rectangle in Fig.\,\ref{FigTeo1}(a), highlighted by Refs.\,\cite{Schoop2015,Wu2018}. In (b) we show the Fermi level cut of the surface bandstructure. The Brillouin zone high symmetry directions are shown as red points connected by dashed lines. The green bar shows approximately the region covered by the green rectangle in (a).}
\label{FigTeo2}
\end{figure}

\section{Discussion}

Our tunneling conductance measurements (Fig.\,\ref{FigZeroField}) show a single quasiparticle peak, showing that the gap magnitude is similar over the whole Fermi surface. However, we do not observe the zero density of states expected at a Fermi level for a usual isotropic s-wave superconductor. Furthermore, the measured $LDOS(E)$ at low temperatures (Fig.\,\ref{FigZeroField}(b)) has a continous V-shaped variation for low energies, instead of a U or squarish shaped curve expected for a usual s-wave superconductor. These two findings differ from previous bulk specific heat measurements, which suggest instead a fully opened superconducting gap \cite{Schoop2015, Yu2016}.

One possible origin for modified superconducting properties at the surface is chemical degradation. This can be excluded in our case, because samples were cleaved at low temperatures in fully inert cryogenic vacuum. Thus, the STM measurements show that superconductivity at or close to the surface in Au$_2$Pb is different than in the bulk.

To discuss this issue, we should first note that the surface roughness we find here is of order of the size of a unit cell (c-axis lattice parameter is of 1.1 nm). Previous ARPES experiments have found electronic surface states in polished surfaces of the room temperature cubic phase, whose roughness is above a unit cell size. Thus, the actual surface roughness does not disturb significantly the formation of surface states  \cite{Wu2018}.

As a first scenario to discuss our results, we consider group theory applied to the low temperature orthorhombic phase. The Dirac cone of the high temperature tetragonal phase along the $\Gamma-X$ line opens in the low temperature orthorhomic phase, due to symmetry breaking allowed band mixing. The influence of this situation in superconductivity was recently analyzed in detail for the case of Au$_2$Pb in Ref.\,\cite{Cheon2021}. For the $D_{2h}$ orthorhombic point group (low temperature crystalline phase), only the superconducting states belonging to the trivial $A_g$ and the triplet $A_u$ representations have no nodes in the bulk and are thus compatible with the available bulk specific heat. Out of these, the $A_u$ state presents a linear Dirac band dispersion inside the bulk superconducting gap at the surface. This surface state is topologically protected by a zero-dimensional topological number that can be defined using the $y-z$ mirror symmetry \cite{PhysRevLett.115.187001,PhysRevB.94.014510,Cheon2021}. The linear Dirac band dispersion inside the superconducting gap generates a linearly dependent superconducting density of states, with massless particles and thus zero states at the Fermi level. Atomic scale impurities within such a two-dimensional Dirac spectrum are known to produce distinct impurity resonances at an energy $E_{dist}$ varying with the strength in the impurity and located within the linear portion of the spectrum (see Fig.\,\ref{Schematics}(c)) \cite{doi:10.1080/00018732.2014.927109}. These resonances have been observed inside the Dirac in-gap dispersion of cuprate superconductors and graphene and are inherent to a Dirac dispersion relation \cite{doi:10.1126/science.285.5424.88,PhysRevLett.104.096804,doi:10.1080/00018732.2014.927109}. Thus, one possible origin for a finite zero bias tunneling conductance is that the surface corrugation generates in-gap states at different energies $E_{dist}$ that overlap and provide a finite zero bias tunneling conductance with slight variations as a function of the position. This scenario would imply, however, a triplet bulk superconducting order parameter belonging to the $A_u$ representation. In absence of bulk spin-dependent measurements, such as Knight shift, there is no firm ground to assume triplet superconductivity in Au$_2$Pb.

As a second scenario, we analyze in more detail the surface band structure in the normal state. We have calculated the surface Green's function using DFT on a semi-infinite Au$_2$Pb sample. We find well defined surface states located in the gaps of the bulk bandstructure for all surface terminations. Furthermore, the spin degeneracy of these surface states is lifted by spin orbit coupling. We show results on the (010) surface in Fig.\,\ref{FigTeo2}(a,b). The $\mathbb Z_2$ topological invariant is one on this plane. We can clearly see that the gap along $\Gamma-X$ is closed at the surface by a Dirac cone. In addition, there are three surface states inside the gap, two of which cross the Fermi level at the center of the Dirac cone. These two surface states cross again the Fermi level close to $X$. As we can see in Fig.\,\ref{FigTeo2}(b), there are a number of additional surface states crossing the Fermi level along different directions of the Brillouin zone. Ginzburg showed that it was in principle possible to find superconductivity on a surface two-dimensional electron gas \cite{Ginzburg2D}. It has been shown that surface states of Ag and Au deposited on top of a superconductor show a superconducting density of states with a non-zero value at the Fermi level \cite{PhysRevB.94.220503,PhysRevLett.122.247002}. Our measurements provide a similar result, suggesting that the opening of the superconducting gap at the surface states is incomplete.

We observe in our images that the surface corrugation is well above atomic size. The surface consists of irregular changes in height. However, absence of atomically flat surfaces does not modify surface states, as mentioned above and shown in Ref.\,\cite{Wu2018}. Still, it induces disorder, which in turn smears out any sharp features in the density of states. The density of states we find in our experiments has a strong V-shape that can also be compatible with gapless s-wave superconductivity at surface states, the triplet order parameter discussed above as well as other possibilities that are not compatible with bulk measurements, such as nodal d-wave superconductivity with impurities\cite{doi:10.1080/00018732.2014.927109}. In any case, the density of states we observe is strongly influenced by the disordered surface topography.

\section{Conclusion}

In conclusion, we have determined the superconducting gap and density of states of Au$_2$Pb. We find a superconducting gap magnitude of $0.17$ meV, compatible with $T_c=1.1$ K and BCS relation between the critical temperature and the gap. The superconducting density of states measured at the surface has a V-shape and a large finite value at zero bias, suggesting that there is a finite zero energy density of states. As possible origins for the finite density of states we discuss bulk triplet superconductivity and the concomitant Dirac states at the surface. We also discuss an intricate structure of surface states with incomplete superconducting gap opening.

\section*{Acknowledgements}

This work was supported by the Spanish Research State Agency (PID2020-114071RB-I00, FIS2017-84330-R, CEX2018-000805-M, RYC-2014-15093, MAT2017-87134-C2-2-R), by the Comunidad de Madrid through program NANOFRONTMAG-CM (S2013/MIT-2850), by the European Research Council PNICTEYES grant agreement 679080 and by the Swedish Research Council (Vetenskapsr\r{a}det Grant No. 2018-03488). We acknowledge collaborations through EU program Cost CA16218 (Nanocohybri). J.J.B. acknowledges the Marie Curie Fellowship program (H2020-MSCA-IF2016-751047) and the Generalitat Valenciana (CDEIGENT/2019/022). Ames Laboratory is operated for the U.S. Department of Energy by Iowa State University under Contract No. DE-AC02-07CH11358. We acknowledge SEGAINVEX at UAM for design and construction of cryogenic equipment. We thank Prof. Ángel Rubio for access to computational resources. We thank Marta Sánchez Lomana and Raquel S\'anchez Barquilla for support during crystal growth. We also thank Rafael \'Alvarez Montoya, Sara Delgado and Jos\'e Mar\'ia Castilla for technical support. We acknowledge SIDI at UAM for support in sample characterization.

\newpage

\section*{Appendix}

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{figTeo3.pdf}
\caption{In (a) we show the calculated surface bandstructure of the (100) as a function of the energy. The density of states is shown by the red-blue color scale. Surface states appear in bright red color and join areas with projected bulk states. In (b) we show the Fermi energy cut in the Brillouin zone high symmetry points indicated by red dots and dashed lines. The same figure are shown for the (001) surface in (c,d).}
\label{FigTeo3}
\end{figure}

\section*{Calculation of the band structure at the surface}

We show the surface bandstructure at the (100) and (001) surfaces in Fig.\,\ref{FigTeo3}. The gapped region around $U$ remains free of surface states close to the Fermi level. But there are numerous surface states at the Fermi level along all other directions. There is for instance a close surface contour that emerges from the pocket of the 17$^{th}$ band on top and bottom of the Brillouin zone. There is also another state connecting the 15$^{th}$ and 17$^{th}$ bands together along the $\Gamma$-U line and the same occurs along the $\Gamma$-X line.

The $\mathbb Z_2$ topological index of the (100) and (010) surfaces is one, while it is zero for the (001) surface.

\section*{Details on the crystalline structure and the surface topography}

In Fig.\,\ref{FigSurface}(a) we present an optical picture of the surface seen after having cleaved the sample at low temperatures. We see triangular like features. We present the room temperature structure of Au$_2$Pb schematically in the left panel of Fig.\,\ref{FigSurface}(b). This is a cubic Laves phase Fd$\overline 3$m with a cell parameter of $a=$ 7.9141(2) \AA \cite{Chen2016}. The structure in the low temperature phase is the Pbcn structure and is derived from the cubic Laves phase by small atomic distortions \cite{Schoop2015}. The intermediate temperature structures (between 40 K and 97 K) are still unknown. Further structural analysis has also lead to the structures Pca2$_1$ and $I\overline{4}2d$ which are also derived from the Laves phase and occur under pressure \cite{Wu2019}. These are shown for completeness in Fig.\,\ref{FigSurface}(b) (bottom panels). Notice that the crystalline transitions induced on cooling are quite varied, as shown in the specific heat (Fig.\,\ref{FigSample}(b)). The transitions at 97 K and at 40 K consist of jumps and the transition at 50 K consists of a large peak, with additional hysteresis, indicating latent heat. The actual atomic arrangements on cooling might lead to changes in the surface crystalline plane. This might influence cleaving, making it difficult to establish a clear cleaving plane at low temperatures, and leading to the absence of atomic size features, as we observe here. The optical picture obtained after cleaving at room temperature (Fig.\,\ref{FigSurface}(a)) shows triangular features, but no clear squarish or linear features as expected in a cubic structure, suggesting that the cleavage was not made along a clear crystalline direction.

We have cleaved two samples and analyzed in each over ten different fields of view, always finding surfaces as discussed in the main text. We often use 10 mV as a bias and work with tunneling currents between 0.1 nA and about 3 nA, finding the same images and tunneling spectroscopy, irrespective of the actual tunneling conductance. We present further images obtained in different fields of view in  Fig.\,\ref{FigSurface}(c,d).

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{figSurface.pdf}
\caption{(a) Optical picture of the surface of the sample made at room temperature after cleaving the sample at 4 K (scale bar is 10 $\mu$m long). (b) Schematic atomic arrangement, with Au in gold and Pb in grey, in the room temperature structure (cubic laves phase, Fd$\overline{3}$m), the low temperature structure (Pbcn) and two further structures that are found under pressure\cite{Wu2019}. All structures derive from the room temperature cubic laves phase. (c,d) STM images obtained at 100 mK on different fields of view. Scale bars are 20 nm long. The corrugation shown in the STM images is of about 2-3 nm, as shown by the color bars on top of (c,d).
}
\label{FigSurface}
\end{figure}

\section*{Tunneling conductance and density of states in a superconductor}

At exactly zero temperature, the tunneling current vs bias voltage obtained with a STM is proportional to the convolution of the densities of states of tip and sample within an energy interval between zero energy and eV, where V is the applied bias voltage and e the elementary charge\cite{Wolf}. As the superconducting gap $\Delta$ is usually much smaller than the Fermi energy $E_F$, $\Delta \ll E_F$, the density of states of the tip is featureless in the energy range of interest in superconducting samples. Thus, the tunneling current is just the integral between zero energy and eV of the density of states of the sample. Derivation leads to the well known relation for the tunneling conductance $G$, $G(V)\propto N_{Superc}(eV)$ where $N_{Superc}(eV)$ is the density of states of the superconductor. In a single gap s-wave superconductor, the density of states is given by $N_{Superc}(eV)=N_{BCS}(E)=\left(\frac{E}{\sqrt{E^2-\Delta^2}}\right)$ for $\vert{E}\vert>\Delta$ and zero elsewhere. At finite temperatures, the Fermi function smears the features of the density of states in the conductance, which is given by $G(V)\propto \int_{-\infty}^{\infty}N_{Superc}(eV)\frac{\partial f(E-eV)}{\partial V}dE$. At the same time, the gap $\Delta$ changes with temperature. For a s-wave BCS superconductor, we can try to follow the main feature of $N_{BCS}(eV)$, the superconducting gap $\Delta(T)$ as a function of temperature. To this end, one might try to trace the maximum in $G(V)$ versus temperature. However, while $\frac{\partial f(E-eV)}{\partial V}$ is a symmetric function at any $eV$, $N_{BCS}(E)$ is definitely not symmetric for $eV\approx \Delta$, because $N_{BCS}(E)\equiv 0$ for $E<\Delta$ and presents a sharp peak for $E>\Delta$. Thus, the maximum in $G(E)$ does not coincide with the maximum in $N_{BCS}(E)$ for finite temperatures $T$. This is valid for features usually found in any $N_{Superc}(E)$ in superconductors that do not follow the simple s-wave single band expression $N_{BCS}(E)$. For example when there are multiple gap peaks in the density of states of multiband superconductors, or for the linearly dispersing $N_{Superc}(E)$ proposed in Fig.\,\ref{Schematics}(c), or for the density of states found in d-wave superconductors. To obtain the temperature dependence of features in $N_{Superc}(E)$, we have to de-convolute $N_{Superc}(eV)$ from $G(V)\propto \int_{-\infty}^{\infty}N_{Superc}(eV)\frac{\partial f(E-eV)}{\partial V}dE$. Using this method, the temperature dependence of the main gap values were found in several multiband superconductors \cite{PhysRevLett.101.166407,PhysRevB.97.134501} and the gap was followed as a function of temperature in a magnetic superconductor with many states inside the gap\cite{PhysRevLett.96.027003}. The method is however best illustrated in a simple s-wave superconductor as Pb. The results in that case are described in Ref.\,\cite{doi:10.1063/5.0064511}.

The BCS shape of $N_{BCS}(eV)$ is found until close to $T_c$ in Pb, i.e. $N_{Superc}(eV)=N_{BCS}(eV)$. In Pb there is a small but non-zero contribution to the density of states from inelastic scattering. One can use then $N_{Dynes}(E)=\Re \left(\frac{E-i\Gamma}{\sqrt{(E-i\Gamma)^2-\Delta^2}}\right)$, with $\Gamma$ being very small in Pb, $\Gamma\approx 0.01$ meV, i.e. less than 1\% of the superconducting gap value, $\Delta=1.3$ meV\cite{PhysRevLett.41.1509}. As shown in Ref.\,\cite{doi:10.1063/5.0064511}, the influence of $\Gamma$ is practically negligible in the experiment, so that $G(V)$ can be traced as a function of temperature by using $N_{BCS}(E)=\left(\frac{E}{\sqrt{E^2-\Delta^2}}\right)$ until very close to $T_c$. It is important to note that, even close to $T_c$ and in spite of the smearing induced by the Fermi function, it is still possible to disentangle details in $N_{exp}(E)$. For example, introducing a $N_{Superc}(E)$ with a finite density of states at $E=0$ at 6 K ($T_c=$ 7.2 K) or taking a different position for the quasiparticle peaks in $N_{Superc}(E)$ does not lead to the measured $G(V)$ in Ref.\,\cite{doi:10.1063/5.0064511}. This also applies to the results shown here in Au$_2$Pb. For example, we show in Fig.\,\ref{FigTunnelTemp}(a) the tunneling conductance $G(V)$ obtained at 0.5 K as black disks. The tunneling conductance obtained by using the same $N_{Superc}(E)$ as the one used at 0.1 K is shown as a blue line. Clearly, there is a temperature induced modification. The tunneling conductance obtained by a modified $N_{Superc}(E)$ (also shown in Fig.\,\ref{FigZeroField}(b) for the curve at 0.5 K) is shown as a red line.

Other than the density of states $N(E)$ proposed in Fig.\,\ref{Schematics}(c), we can also analyze further possibilities. We show in Fig.\,\ref{FigTunnelTemp}(b) the tunneling conductance obtained in the experiment at 0.1 K as black disks. The $N_{Superc}(E)$ used to obtain the black line through convolution with the Fermi function at 0.1 K is shown as a red line. We see a sub-linear dispersion to very low energies and a close to linear increase above about 0.1 mV. We also see that $N_{Dynes}(E)$ (blue curve in Fig.\,\ref{FigTunnelTemp}(b)) leads to a flat energy dependence at low energies which does not reproduce what we find in the experiment. The $N_{Superc}(E)$ curve we find is also compatible with a d-wave order parameter, which provides a density of states that is very similar to the one shown as a red line in Fig.\,\ref{FigTunnelTemp}(b), provided that there are randomly distributed pair breaking states. Nevertheless, macroscopic experiments show an opened superconducting gap, there is no evidence for d-wave superconductivity in the bulk.

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{Fig10N.pdf}
\caption{(a) The normalized tunneling conductance obtained in the experiment at 0.5 K is shown as black disks. The conductance obtained from the convolution of $N_{Superc}(E)$, shown in Fig.\,\ref{FigZeroField}(b) at 0.5 K, and the derivative of the Fermi function is shown as a red line. In blue we show the convolution of the density of states chosen for the lowest temperatures in Fig.\,\ref{FigZeroField}(b) with the Fermi function at a temperature of 0.5 K. We see that the density of states has a temperature dependence that is well captured by the deconvolution procedure leading to Fig.\,\ref{FigZeroField}(b). (b) The tunneling conductance obtained at 0.1 K is shown as black discs. The black line is the convolution of the red line (which is $N_{Superc}(E)$ for 0.1 K in Fig.\,\ref{FigZeroField}(b)) with the Fermi function at a temperature of 0.1 K. For comparison, we provide in blue the density of states following Dynes inelastic scattering expression $N_{Dynes}(E)$ that best follows $N_{Superc}(E)$, with $\Gamma=$0.05 meV and $\Delta=0.16$ meV.
}
\label{FigTunnelTemp}
\end{figure}

\section*{Tunneling conductance in a magnetic field}

To further study the superconducting properties we measure the tunneling spectroscopy in a magnetic field. We observe that superconductivity disappears most often between 20 mT and 30 mT (Fig.\,\ref{FigField}), depending on the field of view. We could find no indications for vortices or of separated normal and superconducting areas within the fields of view we investigated under magnetic fields. However, we observe a clear hysteresis when measuring the magnetic field dependence of the tunneling conductance at a fixed position (Fig.\,\ref{FigField}). Previous bulk measurements have found a first critical field of about 30 mT. There are reports about the presence of a second critical field in the range between 20 mT and 80 mT \cite{Schoop2015,Xing2016,Yu2016,Wu2018}. From these data we can estimate a coherence length of at least $\xi\approx 60$ nm.

There is clearly Meissner field expulsion until several tens of mT at low temperatures \cite{Xing2016}. Thus, we can expect strong demagnetizing effects, similar to those that can be found in elemental Pb or in Nb at small magnetic fields \cite{Prozorov2008,PhysRevLett.102.136408}. Our sample is nearly plate-like, with a demagnetizing factor close to one. This implies that the magnetic field enters the sample well below the first critical field and that we can expect an inhomogeneous field distribution and hysteretic behavior when increasing and decreasing the magnetic field due to trapped flux \cite{RevModPhys.91.015004,BRANDT349,PhysRevB.3.2231}. The large value of the coherence length suggests that the separation of the sample into normal and superconducting areas, characteristic of the intermediate state \cite{Prozorov2008,PhysRevLett.102.136408,garciacampos2020visualization}, occurs over length scales that are much larger than our fields of view. The hysteresis shows however that such well separated normal and superconducting areas are present in our sample. We also notice that small angle neutron scattering experiments did not lead to the observation of a vortex lattice \cite{Emma22}.

\begin{figure}[htbp]
\includegraphics[width = \columnwidth]{figSTMField.pdf}
\caption{(a) Magnetic field dependence of the superconducting tunneling conductance, with the values of the magnetic field indicated. (b) Zero bias conductance as a function of the magnetic field for increasing (green) and decreasing (dark red) magnetic fields. (c,d) Same as (a,b) in a different field of view. The inset in (d) shows a tunneling conductance map at zero bias voltage at 100 mK and at a (decreasing) magnetic field of 0.01 T. Horizontal scale bar (black line) is 50 nm long.
}
\label{FigField}
\end{figure}

Clearly, Au$_2$Pb is a type I or weakly type II superconductor. This adds Au$_2$Pb to the rapidly increasing list of non-elemental materials that are type I or weak type II superconductors. We can find reports indicating such a behavior in Al$_6$Re, YNiSi$_3$, LuNiSi$_3$, SnAs, RuB$_2$, KBi$_2$, PdTe$_2$, CaBi$_2$, $\beta$-IrSn$_4$, RPd$_2$Si$_2$ ($R=Y, La, Lu$), LaRh$_2$Si$_2$, TaSe$_2$, Ag$_5$Pb$_2$O$_6$, LaRhSi$_3$, ScGa$_3$, LuGa$_3$, YbSb$_2$, BeAu, LiBi, potassium graphite intercalation compounds or SrRu$_2$O$_4$ (the latter two along a given field direction) \cite{PhysRevB.99.144519,PhysRevB.99.224505,PhysRevB.100.184514,PhysRevB.97.054506,Sun_2016,Chen2018,PhysRevB.96.220506,PhysRevB.94.144506,C6CP02856J,Tran_2013,doi:10.1021/acs.chemmater.0c00179,doi:10.1143/JPSJ.50.3245,PhysRevB.34.4566,PhysRevB.45.4803,PhysRevB.72.180504,PhysRevB.83.064522,PhysRevB.85.174514,doi:10.1143/JPSJ.56.419,PhysRevB.99.134509,TakashiAkima1999,garciacampos2020visualization}.