\section{Introduction}

Let $D^{2}u$ be the Hessian matrix of a $C^{2}$ (i.e., a twice continuously
differentiable) function $u$ defined over $\mathbb{R}^{N}$ ($N\geq 3$) and
\lambda \left( D^{2}u\right) =\left( \lambda _{1},...,\lambda _{N}\right) $
the vector of eigenvalues of $D^{2}u$. For $k=1,2,...,N$ is defined the $k
-Hessian operator as follow
\begin{equation*}
S_{k}\left( \lambda \left( D^{2}u\right) \right) =\underset{1\leq
i_{1}<...<i_{k}\leq N}{\sum }\lambda _{i_{1}}\cdot ...\cdot \lambda _{i_{k}}
\end{equation*
i.e., it is the $k^{th}$ elementary symmetric polynomial of the Hessian
matrix of $u$. In other words, $S_{k}\left( \lambda \left( D^{2}u\right)
\right) $ it is the sum of all $k\times k$ principal minors of the Hessian
matrix $D^{2}u$ and so is a second order differential operator, which may
also be called the $k$-trace of $D^{2}u$. Especially, it is easily to see
that the $N$-Hessian is the Monge-Amp\'{e}re operator and that the $1
-Hessian is the well known classical Laplace operator. Hence, the $k
-Hessian operators form a discrete collection of partial differential
operators which includes both the Laplace and the Monge-Amp\'{e}re operator.

In this paper we study the existence of radial solutions for the following
Hessian equation
\begin{equation}
S_{k}^{1/k}\left( \lambda \left( D^{2}u\right) \right) =p\left( \left\vert
x\right\vert \right) h\left( u\right) \text{ in }\mathbb{R}^{N}\text{ ,}
\label{11}
\end{equation
and syste
\begin{equation}
\left\{
\begin{array}{c}
S_{k}^{1/k}\left( \lambda \left( D^{2}u\right) \right) =p\left( \left\vert
x\right\vert \right) f\left( u,v\right) \text{ in }\mathbb{R}^{N}\text{ ,}
\\
S_{k}^{1/k}\left( \lambda \left( D^{2}v\right) \right) =q\left( \left\vert
x\right\vert \right) g\left( u,v\right) \text{ in }\mathbb{R}^{N}\text{ ,
\end{array
\right.  \label{11s}
\end{equation
where $k\in \left\{ 1,2,...,N\right\} $, the continuous functions $p$, $q
\textit{\ }$:\left[ 0,\infty \right) \rightarrow \left( 0,\infty \right) $,
h:\left[ 0,\infty \right) \rightarrow \left[ 0,\infty \right) $ and $f$, $g
\left[ 0,\infty \right) \times \left[ 0,\infty \right) \rightarrow \left[
0,\infty \right) $ satisfy some of the conditions:

(P1)\quad $p$, $q$ is a spherically symmetric function (i.e.\textit{\ }
p\left( x\right) =p\left( \left\vert x\right\vert \right) $, $q\left(
x\right) =q\left( \left\vert x\right\vert \right) );$

(P2)\quad $r^{N+\frac{N}{k}-2}p^{k}\left( r\right) $\textit{\ }is
nondecreasing for large $r$;

(P3)\quad $r^{N+\frac{N}{k}-2}\left[ p^{k}\left( r\right) +q^{k}\left(
r\right) \right] $\textit{\ }is nondecreasing for large $r$\textit{;}

(C1)\quad $h$ is monotone non-decreasing, $h(0)=0$ and $h\left( s\right) >0$
for all $s>0$;

(C2)\quad $f$, $g$ are monotone non-decreasing in each variable,
f(0,0)=g\left( 0,0\right) =0$ and $f(s,t)>0,$ $g\left( s,t\right) >0$ for
all $s,t>0$;

(C3)\quad $\int_{1}^{\infty }\frac{1}{\sqrt[k+1]{\left( k+1\right) H(t)}
dt=\infty $ for\ $H(t)=\int_{0}^{t}h^{k}(z)dz;$

(C4)\quad $\int_{1}^{\infty }\frac{1}{\sqrt[k+1]{\left( k+1\right) F(t)}
dt=\infty $ for\ $F(t)=\int_{0}^{t}\left( f^{k}(z,z)+g^{k}(z,z)\right) dz.$

The properties of the $k$-Hessian operator was well discussed in a numerous
papers written as a first author by Ivochkina (see \cite{III}-\cite{IF} and
others). Moreover, this operator appear as an object of investigation by
many remarkable geometers. For example, Viaclovsky (see \cite{V}, \cite{VI})
observed that the $k-$Hessian operator is an important class of fully
nonlinear operators which is closely related to a geometric problem of the
type (\ref{11}), where we cite the work of Bao-Ji-Li \cite{BAOII} for a more
detailed discussion. Moreover, equation (\ref{11}) arises via the study of
the quasilinear parabolic problem ( see for example the introduction of
Moll-Petitta \cite{MP}). In the present work we will limit ourselves to the
development of mathematical theory for (\ref{11}) and (\ref{11s}). The main
difficulty in investigating problems, such as (\ref{11}) or (\ref{11s}), in
which appear the $k$-Hessian operator is related to the fact that their
properties change depending on the subset of $C^{2}$ from where the solution
is taken. Our main objective here is to find functions in $C^{2}$ that are
strictly $k$-convex and verifies the problems (\ref{11}), (\ref{11s}), where
by strictly $k$-convex function $u$ we mean that all eigenvalues $\lambda
_{1},...,\lambda _{N}$ of the symmetric matrix $D^{2}u$ are in the so called
G\aa rdding open cone $\Gamma _{k}$ which is defined b
\begin{equation*}
\Gamma _{k}\left( N\right) =\left\{ \lambda \in \mathbb{R}^{N}\left\vert
S_{1}\left( \lambda \right) >0,....,S_{k}\left( \lambda \right) >0\right.
\right\}
\end{equation*
In the next we adopt the notation from Bao-Li \cite{BAO} for the space of
all admissible functions
\begin{equation*}
\Phi ^{k}\left( \mathbb{R}^{N}\right) :=\left\{ u\in C^{2}\left( \mathbb{R
^{N}\right) \left\vert \lambda \in \Gamma _{k}\left( N\right) \text{ for all
}x\in \mathbb{R}^{N}\right. \right\} .
\end{equation*
In our direction, there are some recently papers resolving existence for
blow-up solutions of (\ref{11}) and (\ref{11s}). Here we wish to mention the
works of Bao-Ji-Li \cite{BAOII}, Jacobsen \cite{J}, Bao-Li \cite{BAO}, Lazer
and McKenna \cite[(the case $k=N$)]{LM}, Salani \cite{S} and Zhang-Zhou \cit
{ZZ} which will be useful in our proofs. It is interesting to note that in
our results the dimension of the space $\mathbb{R}^{N}$ affect the
properties of the solution of the equation and system which in the case of
the classical Laplace operator and the Monge-Amp\'{e}re operator this
condition doesn't appear in any works.

Motivated by the recent work of Zhang-Zhou \cite{ZZ} we are interested in
proving the following theorems:

\begin{theorem}
\label{th1}Let $k\in \left\{ 1,2,...,\left[ N/2\right] \right\} $ if $N$ is
odd or $k\in \left\{ 1,2,...,\left[ N/2\right] -1\right\} $ if $N$ is even.
Suppose that \textrm{(P1), (P2)}, \textrm{(C1)}, \textrm{(C3)} are
satisfied. If there exists a positive number $\varepsilon $ such that
\begin{equation}
\quad \text{ }\int_{0}^{\infty }t^{1+\varepsilon +\frac{2\left( k-1\right) }
k+1}}\left( p\left( t\right) \right) ^{\frac{2k}{k+1}}dt<\infty ,  \label{5}
\end{equation
then system \textrm{(\ref{11})} has a nonnegative nontrivial radial bounded
solution $u\in \Phi ^{k}\left( \mathbb{R}^{N}\right) $.
\end{theorem}

\begin{theorem}
\label{th2}\textit{If }$p$ satisfy \textit{\textrm{(P1)}} and\ $f$ \textit
satisfy \textrm{(C1)}, \textrm{(C3)}, then the problem \textrm{(\ref{11})}
has a nonnegative nontrivial entire radial solution }$u\in \Phi ^{k}\left(
\mathbb{R}^{N}\right) $\textit{. Suppose furthermore that \textrm{(P2)}}
holds.\textit{\ }If $p$ satisfies
\begin{equation}
\int_{0}^{\infty }\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac{s^{N-1}}
C_{N-1}^{k-1}}p^{k}\left( s\right) ds\right) ^{1/k}dt=\infty ,\text{ }
\label{12}
\end{equation
then any nonnegative nontrivial radial solution $u\in \Phi ^{k}\left(
\mathbb{R}^{N}\right) $ of \textrm{(\ref{11})} is large. Conversely, \textit
if }\textrm{(\ref{11})} has a nonnegative entire large radial solution $u\in
\Phi ^{k}\left( \mathbb{R}^{N}\right) $, then one or both of the following
\begin{equation}
\begin{array}{ll}
1. & \int_{0}^{\infty }t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1
}\left( p\left( t\right) \right) ^{\frac{2k}{k+1}}dr=\infty \text{ for every
}\varepsilon >0\,; \\
2. & k\in \left\{ \left[ N/2\right] +1,...,N\right\} \text{ if }N\text{ is
odd or }k\in \left\{ \left[ N/2\right] ,...,N\right\} \text{ if }N\text{ is
even,
\end{array}
\label{13}
\end{equation
hold.
\end{theorem}

Regarding existence of solution to (\ref{11s}), we have the following
results.

\begin{theorem}
\label{th3}Let $k\in \left\{ 1,2,...,\left[ N/2\right] \right\} $ if $N$ is
odd or $k\in \left\{ 1,2,...,\left[ N/2\right] -1\right\} $ if $N$ is even.
Suppose that \textrm{(P1), (P3)}, \textrm{(C2)}, \textrm{(C4)} are
satisfied. If there exists a positive number $\varepsilon $ such that
\begin{equation}
\text{ }\int_{0}^{\infty }t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1
}\left( p^{k}\left( t\right) +q^{k}\left( t\right) \right) ^{\frac{2}{k+1
}dt<\infty \text{ },  \label{5s}
\end{equation
then system \textrm{(\ref{11s})} has a nonnegative nontrivial radial bounded
solution $\left( u,v\right) \in \Phi ^{k}\left( \mathbb{R}^{N}\right) \times
\Phi ^{k}\left( \mathbb{R}^{N}\right) $.
\end{theorem}

\begin{theorem}
\label{th4}\textit{If }$p,$ $q$ satisfy \textit{\textrm{(P1)}} and\textit{\
\ $f,$ $g$ \textit{satisfy \textrm{(C2)}, \textrm{(C4)}, then the problem
\textrm{(\ref{11})} has a nonnegative nontrivial entire radial solution.
Suppose furthermore that \textrm{(P3)} holds. }If $p$ satisfies
\begin{equation}
\int_{0}^{\infty }\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac{s^{N-1}}
C_{N-1}^{k-1}}p^{k}\left( s\right) ds\right) ^{1/k}dt=\infty \text{ and
\int_{0}^{\infty }\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac{s^{N-1}}
C_{N-1}^{k-1}}q^{k}\left( s\right) ds\right) ^{1/k}dt=\infty ,  \label{12s}
\end{equation
then any nonnegative nontrivial solution $\left( u,v\right) \in \Phi
^{k}\left( \mathbb{R}^{N}\right) \times \Phi ^{k}\left( \mathbb{R
^{N}\right) $ of \textrm{(\ref{11s})} is large. Conversely, if \textrm{(\re
{11s})} has a nonnegative entire large radial solution $\left( u,v\right)
\in \Phi ^{k}\left( \mathbb{R}^{N}\right) \times \Phi ^{k}\left( \mathbb{R
^{N}\right) $, then one or both of the following
\begin{equation}
\begin{array}{ll}
1. & \int_{0}^{\infty }t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1
}\left( p^{k}\left( t\right) +q^{k}\left( t\right) \right) ^{\frac{2}{k+1
}dr=\infty \text{ for every }\varepsilon >0\,; \\
2. & k\in \left\{ \left[ N/2\right] +1,...,N\right\} \text{ if }N\text{ is
odd or }k\in \left\{ \left[ N/2\right] ,...,N\right\} \text{ if }N\text{ is
even,
\end{array}
\label{13s}
\end{equation
hold.
\end{theorem}

For the readers' convenience, we recall the radial form of the $k$-Hessian
operator.

\begin{remark}
(see, for example, \cite{BAO}, \cite{S}) If $u:\mathbb{R}^{N}\rightarrow
\mathbb{R}$ is radially symmetric then a calculation sho
\begin{equation*}
S_{k}\left( \lambda \left( D^{2}u\left( r\right) \right) \right)
=r^{1-N}C_{N-1}^{k-1}\left[ \frac{r^{N-k}}{k}\left( u^{^{\prime }}\left(
r\right) \right) ^{k}\right] ^{\prime }\text{ ,}
\end{equation*
where the prime denotes differentiation with respect to $r=\left\vert
x\right\vert $ and $C_{N-1}^{k-1}=(N-1)!/\left[ (k-1)!(N-k)!\right] $.
\end{remark}

\section{Proofs of the main results}

In this section we give the proofs of Theorems \ref{th1} - \ref{th4}. The
main references for proving Theorems \ref{th1} - \ref{th2} is the work of
Lair \cite{LA} and Delano\"{e} \cite{DEL} see also Afrouzi-Shokooh \cite{AF}.

\subparagraph{Proof of the Theorem \protect\ref{th1}}

Assume that (\ref{5}) holds. We prove the existence of $w\in \Phi ^{k}\left(
\mathbb{R}^{N}\right) $ to the problem

\begin{equation}
\begin{array}{l}
S_{k}^{1/k}\left( \lambda \left( D^{2}w\left( \left\vert x\right\vert
\right) \right) \right) =p\left( \left\vert x\right\vert \right) h\left(
w\left( \left\vert x\right\vert \right) \right) \text{ in }\mathbb{R}^{N
\text{.
\end{array}
\label{6}
\end{equation
Observe that we can rewrite (\ref{6}) as follows:
\begin{equation*}
\left[ \frac{r^{N-k}}{k}\left( w^{^{\prime }}(r)\right) ^{k}\right] ^{\prime
}=\frac{r^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( r\right) h^{k}\left( w\left(
r\right) \right) ,\text{ }r=\left\vert x\right\vert .
\end{equation*
Then radial solutions of (\ref{6}) are any solution $w$ of the integral
equatio
\begin{equation*}
w\left( r\right) =1+\int_{0}^{r}\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac
s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left( w\left( s\right)
\right) ds\right) ^{1/k}dt.\text{ }
\end{equation*
To establish a solution to this problem, we use successive approximation.
Define sequence $\left\{ w^{m}\right\} ^{m\geq 1}$ on $\left[ 0,\infty
\right) $ b
\begin{equation*}
\left\{
\begin{array}{l}
w^{0}=1,\text{ }r\geq 0, \\
w^{m}\left( r\right) =1+\int_{0}^{r}\left( \frac{k}{t^{N-k}}\int_{0}^{t
\frac{s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left( w^{m-1}\left(
s\right) \right) ds\right) ^{1/k}dt
\end{array
\right.
\end{equation*
We remark that, for all $r\geq 0$ and $m\in N$
\begin{equation*}
w^{m}\left( r\right) \geq 1\text{.}
\end{equation*
Moreover, proceeding by induction we conclude $\left\{ w^{m}\right\} ^{m\geq
1}$ are non-decreasing sequence on $\left[ 0,\infty \right) $. We note that
\left\{ w^{m}\right\} ^{m\geq 1}$ satisfie
\begin{equation*}
\left\{ \frac{r^{N-k}}{k}\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right]
^{k}\right\} ^{\prime }=\frac{r^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( r\right)
h^{k}\left( w^{m-1}\left( r\right) \right) .
\end{equation*
By the monotonicity of $\left\{ w^{m}\right\} ^{m\geq 1}$ we have the
inequalitie
\begin{equation}
\left\{ \frac{r^{N-k}}{k}\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right]
^{k}\right\} ^{\prime }=\frac{r^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( r\right)
h^{k}\left( w^{m-1}\left( r\right) \right) \leq \frac{r^{N-1}}{C_{N-1}^{k-1}
p^{k}\left( r\right) h^{k}\left( w^{m}\left( r\right) \right) .  \label{8}
\end{equation
Choose $R>0$ so that $r^{N+\frac{N}{k}-2}p^{k}\left( r\right) $ are
non-decreasing for $r\geq R$. We are now ready to show that $w^{m}\left(
R\right) $ and $\left( w^{m}\left( R\right) \right) ^{\prime }$, both of
which are nonnegative, are bounded above independent of $m$. To do this, let
\begin{equation*}
\phi ^{R}=\max \{p^{k}\left( r\right) :0\leq r\leq R\}\text{.}
\end{equation*
Using this and the fact that $\left( w^{m}\right) ^{\prime }\geq 0$, we note
that (\ref{8}) yield
\begin{eqnarray*}
r^{N-k}\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k-1}\left(
w^{^{m}}(r)\right) ^{\prime \prime } &\leq &\frac{N-k}{k}r^{N-k-1}\left[
\left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k}+r^{N-k}\left[ \left(
w^{^{m}}(r)\right) ^{\prime }\right] ^{k-1}\left( w^{^{m}}(r)\right)
^{\prime \prime } \\
&\leq &\phi _{1}^{R}\frac{r^{N-1}}{C_{N-1}^{k-1}}h^{k}\left( w^{m}\left(
r\right) \right) ,
\end{eqnarray*
and moreove
\begin{equation*}
r^{N-k}\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k-1}\left(
w^{^{m}}(r)\right) ^{\prime \prime }\leq \phi ^{R}\frac{r^{N-k+k-1}}
C_{N-1}^{k-1}}h^{k}\left( w^{m}\left( r\right) \right) \leq R^{k-1}\phi ^{R
\frac{r^{N-k}}{C_{N-1}^{k-1}}h^{k}\left( w^{m}\left( r\right) \right) ,
\end{equation*
from which we hav
\begin{equation*}
\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k-1}\left(
w^{^{m}}(r)\right) ^{\prime \prime }\leq R^{k-1}\phi ^{R}\frac{1}
C_{N-1}^{k-1}}h^{k}\left( w^{m}\left( r\right) \right) .
\end{equation*
Multiply this by $\left( w^{^{m}}(r)\right) ^{\prime }$ we obtai
\begin{equation}
\left\{ \left[ \left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k+1}\right\}
^{\prime }\leq \frac{\left( k+1\right) R^{k-1}\phi ^{R}}{C_{N-1}^{k-1}
h^{k}\left( w^{m}\left( r\right) \right) \left( w^{^{m}}(r)\right) ^{\prime
}.  \label{ec}
\end{equation
Integrate (\ref{ec}) from $0$ to $r$ to ge
\begin{equation}
\left[ \left( w^{^{m}}(r)\right) ^{\prime }\right] ^{k+1}\leq \frac{\left(
k+1\right) R^{k-1}\phi ^{R}}{C_{N-1}^{k-1}}\int_{0}^{r}h^{k}\left(
w^{m}\left( s\right) \right) \left( w^{^{m}}(s)\right) ^{\prime }ds=\frac
\left( k+1\right) R^{k-1}\phi ^{R}}{C_{N-1}^{k-1}
\int_{1}^{w^{^{m}}(r)}h^{k}\left( s\right) ds  \label{9}
\end{equation
for $0\leq r\leq R$, which yield
\begin{equation*}
\int_{1}^{w^{^{m}}(R)}\left[ \int_{1}^{t}h^{k}\left( s\right) ds\right]
^{-1/\left( k+1\right) }dt\leq \sqrt[k+1]{\frac{\left( k+1\right) \phi ^{R}}
C_{N-1}^{k-1}}}\cdot R^{\frac{2k}{k+1}}.
\end{equation*
It follows from the above relation and by the assumption (C2) that
w_{1}^{m}\left( R\right) $ is bounded above independent of $m$. Using this
fact in (\ref{9}) shows that the same is true of $\left( w^{m}\left(
R\right) \right) ^{\prime }$. Thus, the sequences $w^{m}\left( R\right) $
and $\left( w^{m}\left( R\right) \right) ^{\prime }$ are bounded above
independent of $m$.

Finally, we show that the non-decreasing sequences $w^{m}$ is bounded for
all $r\geq 0$ and all $m$. Multiplying the equation (\ref{8}) by $r^{N+\frac
N}{k}-2}\left( w^{^{m}}(r)\right) ^{\prime }$, we ge
\begin{equation}
\left\{ \left[ r^{\frac{N}{k}-1}\left( w^{m}\left( r\right) \right) ^{\prime
}\right] ^{k+1}\right\} ^{\prime }=\frac{k+1}{C_{N-1}^{k-1}}p^{k}\left(
r\right) h^{k}(w^{m}\left( r\right) )r^{N+\frac{N}{k}-2}\left( w^{m}\left(
r\right) \right) ^{\prime }.  \label{ma14}
\end{equation
Integrating from $R$ to $r$ give
\begin{equation*}
\left[ r^{\frac{N}{k}-1}\left( w^{m}\left( r\right) \right) ^{\prime }\right]
^{k+1}=\left[ R^{\frac{N}{k}-1}\left( w^{m}\left( R\right) \right) ^{\prime
\right] ^{k+1}+\frac{k+1}{C_{N-1}^{k-1}}\int_{R}^{r}p^{k}\left( s\right)
h^{k}(w^{m}\left( s\right) )s^{N+\frac{N}{k}-2}\left( w^{m}\left( s\right)
\right) ^{\prime }ds,
\end{equation*
for $r\geq R$. Noting that, by the monotonicity of $s^{N+\frac{N}{k
-2}p^{k}\left( s\right) $ for $r\geq s\geq R$, we ge
\begin{equation*}
\left[ r^{\frac{N}{k}-1}\left( w^{m}\left( r\right) \right) ^{\prime }\right]
^{k+1}\leq C+\frac{k+1}{C_{N-1}^{k-1}}r^{N+\frac{N}{k}-2}p^{k}\left(
r\right) H\left( w^{m}\left( r\right) \right)
\end{equation*
where $C=\left[ R^{\frac{N}{k}-1}\left( w^{m}\left( R\right) \right)
^{\prime }\right] ^{k+1}$, which yield
\begin{equation*}
r^{\frac{N}{k}-1}\left( w^{m}\left( r\right) \right) ^{\prime }\leq C^{\frac
1}{k+1}}+\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}r^{\frac{N}
k}-\frac{2}{k+1}}p^{\frac{k}{k+1}}\left( r\right) H^{\frac{1}{k+1}}\left(
w^{m}\left( r\right) \right)
\end{equation*
or, equivalentl
\begin{equation*}
\left( w^{m}\left( r\right) \right) ^{\prime }\leq C^{\frac{1}{k+1}}r^{1
\frac{N}{k}}+\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}r^{1
\frac{2}{k+1}}p^{\frac{k}{k+1}}\left( r\right) H^{\frac{1}{k+1}}\left(
w^{m}\left( r\right) \right)
\end{equation*
and henc
\begin{equation}
\frac{d}{dr}\int_{w^{m}\left( R\right) }^{w^{m}\left( r\right) }\left[
H\left( t\right) \right] ^{-1/\left( k+1\right) }dt\leq C^{\frac{1}{k+1
}r^{1-\frac{N}{k}}H^{-\frac{1}{k+1}}\left( w^{m}\left( r\right) \right)
+\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}\left(
r^{k-1}p^{k}\left( r\right) \right) ^{\frac{1}{k+1}}.  \label{in}
\end{equation
Inequality (\ref{in}) combined wit
\begin{eqnarray*}
\frac{1}{\sqrt{2}}\sqrt{2\cdot \left( s^{k-1}p^{k}\left( s\right) \right) ^
\frac{2}{k+1}}} &=&\frac{1}{\sqrt{2}}\sqrt{2\cdot s^{\frac{1+\varepsilon }{2
}\left( s^{k-1}p^{k}\left( s\right) \right) ^{\frac{2}{k+1}}s^{\frac
1-\varepsilon }{2}}} \\
&\leq &\frac{1}{\sqrt{2}}\left[ s^{1+\varepsilon }\left( s^{k-1}p^{k}\left(
s\right) \right) ^{\frac{2}{k+1}}+s^{-1-\varepsilon }\right]
\end{eqnarray*
give
\begin{equation*}
\begin{array}{ll}
\int_{w^{m}\left( R\right) }^{w^{m}\left( r\right) }\left[ H\left( t\right)
\right] ^{-1/\left( k+1\right) }dt & \leq C^{\frac{1}{k+1}}\int_{R}^{r}t^{1
\frac{N}{k}}H^{-\frac{1}{k+1}}\left( w^{m}\left( t\right) \right) dt \\
& +\frac{1}{\sqrt{2}}\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}
\left[ \int_{R}^{r}t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1}}\left(
p\left( t\right) \right) ^{\frac{2k}{k+1}}dt+\int_{R}^{r}t^{-1-\varepsilon
}dt\right]  \\
& \leq C^{\frac{1}{k+1}}H^{-\frac{1}{k+1}}\left( w^{m}\left( R\right)
\right) \int_{R}^{r}t^{1-\frac{N}{k}}dt \\
& +\frac{1}{\sqrt{2}}\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}
\left[ \int_{R}^{r}t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1}}\left(
p\left( t\right) \right) ^{\frac{2k}{k+1}}dt+\frac{1}{\varepsilon
R^{\varepsilon }}\right] \text{.
\end{array
\end{equation*
The above relation is needed in proving the bounded of the function $\left\{
w^{m}\right\} ^{m\geq 1}$ in the following. Indeed, since for each
\varepsilon >0$ the right side of this inequality is bounded independent of
m$ (note that $w^{m}\left( t\right) \geq 1$), so is the left side and hence,
in light of (C2), the sequence $\left\{ w^{m}\right\} ^{m\geq 1}$ is a
bounded sequence and so $\left\{ w^{m}\right\} ^{m\geq 1}$ are bounded
sequence. Thus $\left\{ w^{m}\right\} ^{m\geq 1}\rightarrow w$ as
m\rightarrow \infty $ and the limit functions $w$ are positive entire
bounded solutions of equation (\ref{6}).

\subparagraph{Proof of the Theorem \protect\ref{th2}.}

We know that for any $a_{1}>0$ a solution o
\begin{equation*}
v\left( r\right) =a_{1}+\int_{0}^{r}\left( \frac{k}{t^{N-k}}\int_{0}^{t
\frac{s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left( v\left(
s\right) \right) ds\right) ^{1/k}dt,
\end{equation*
exists, at least, small $r$. Since $v^{\prime }\geq 0$, the only way that
the solution can become singular at $R$ is for $v\left( r\right) \rightarrow
\infty $ as $r\longrightarrow \infty $. Thus, we can show that, for each $R>0
$, there exists $C_{R}>0$ so that $v\left( R\right) \leq C_{R}$, we have
existence. To this end, let $M_{R}=\max \left\{ p\left( r\right) \left\vert
0\leq r\leq R\right. \right\} $ and consider the equatio
\begin{equation*}
w\left( r\right) =a_{2}+M_{R}\int_{0}^{r}\left( \frac{k}{t^{N-k}}\int_{0}^{t
\frac{s^{N-1}}{C_{N-1}^{k-1}}h^{k}\left( v\left( s\right) \right) ds\right)
^{1/k}dt
\end{equation*
where $a_{2}>a_{1}$. We next observe that the solution to this equation
exists for all $r\geq 0$ and of course, it is a solution to
S_{k}^{1/k}\left( \lambda \left( D^{2}w\left( r\right) \right) \right)
=M_{R}h\left( w\right) $ on $\mathbb{R}^{N}$ which is treated in \cite
(Theorem 1.1, p. 177)]{BAO}. We now show that $v\left( r\right) \leq w\left(
r\right) $ for all $0\leq r\leq R$ and hence we conclude the proof of
existence. Clearly $v\left( 0\right) <w\left( 0\right) $ so that $v\left(
r\right) <w\left( r\right) $ for at least all $r$ near zero. Let
\begin{equation*}
r_{0}=\sup \left\{ r\left\vert v\left( s\right) <w\left( s\right) \text{ for
all }s\in \left[ 0,r\right] \right. \right\} \text{.}
\end{equation*
If $r_{0}=R$, then we are done. Thus assume that $r_{0}<R$. It follows from
assumption $a_{2}>a_{1}$ that
\begin{eqnarray*}
v\left( r_{0}\right)  &=&a_{1}+\int_{0}^{r_{0}}\left( \frac{k}{t^{N-k}
\int_{0}^{t}\frac{s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left(
v\left( s\right) \right) ds\right) ^{1/k}dt \\
&<&a_{2}+M_{R}\int_{0}^{r_{0}}\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac
s^{N-1}}{C_{N-1}^{k-1}}h^{k}\left( v\left( s\right) \right) ds\right)
^{1/k}dt=w\left( r_{0}\right) .
\end{eqnarray*
Thus there exists $\varepsilon >0$ so that $v\left( r\right) <w\left(
r\right) $ for all $\left[ 0,r+\varepsilon \right) $, contradicting the
definition of $r_{0}$. Thus we conclude that $v<w$ on $\left[ 0,R\right] $
for all $R>0$ and hence $v$ is a nontrivial entire solution of (\ref{11}).
Now let $u$ be any nonnegative nontrivial entire solution of (\ref{11}) and
suppose $p$ satisfies
\begin{equation*}
\int_{0}^{\infty }\left( \frac{k}{t^{N-k}}\int_{0}^{t}\frac{s^{N-1}}
C_{N-1}^{k-1}}p^{k}\left( s\right) ds\right) ^{1/k}dt=\infty .
\end{equation*
Since $u$ is nontrivial and non-negative, there exists $R>0$ so that
u\left( R\right) >0$. On the other hand since $u^{\prime }\geq 0$, we get
u\left( r\right) \geq u\left( R\right) $ for $r\geq R$ and thus fro
\begin{equation*}
u\left( r\right) =u\left( 0\right) +\int_{0}^{r}\left( \frac{k}{t^{N-k}
\int_{0}^{t}\frac{s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left(
u\left( s\right) \right) ds\right) ^{1/k}dt,\text{ }
\end{equation*
since $u$ will satisfy that equation for all $r\geq 0,$ we get
\begin{equation*}
\begin{array}{l}
u\left( r\right) =u\left( 0\right) +\int_{0}^{r}\left( \frac{k}{t^{N-k}
\int_{0}^{t}\frac{s^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( s\right) h^{k}\left(
u\left( s\right) \right) ds\right) ^{1/k}dt \\
\text{ \ \ \ \ \ \ }\geq u\left( R\right) +h\left( u\left( R\right) \right)
\int_{R}^{r}\left( \frac{k}{t^{N-k}}\int_{R}^{t}\frac{s^{N-1}}{C_{N-1}^{k-1}
p^{k}\left( s\right) ds\right) ^{1/k}dt\rightarrow \infty \text{ as
r\rightarrow \infty
\end{array
\end{equation*
Conversely, assume that\ $h$ satisfy (C1), (C3) and that $w$ is a
nonnegative entire large solution of (\ref{11}). Note also, that $w$
satisfie
\begin{equation*}
\left[ \frac{r^{N-k}}{k}\left( w^{^{\prime }}(r)\right) ^{k}\right] ^{\prime
}=\frac{r^{N-1}}{C_{N-1}^{k-1}}p^{k}\left( r\right) h^{k}\left( w\left(
r\right) \right) .
\end{equation*
Using the monotonicity of \textit{\ }$r^{N+\frac{N}{k}-2}p\left( r\right)
\textit{\ } we can apply similar arguments used in obtaining Theorem \re
{th1} to ge
\begin{equation*}
\left( w\left( r\right) \right) ^{\prime }\leq C^{\frac{1}{k+1}}r^{1-\frac{
}{k}}+\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}r^{1-\frac{2}
k+1}}p^{\frac{k}{k+1}}\left( r\right) H^{\frac{1}{k+1}}\left( w\left(
r\right) \right) ,
\end{equation*
which we may rewrite a
\begin{equation}
\begin{array}{ll}
\int_{w\left( R\right) }^{w\left( r\right) }\left[ H\left( t\right) \right]
^{-1/\left( k+1\right) }dt & \leq C^{\frac{1}{k+1}}\int_{R}^{r}t^{1-\frac{N}
k}}H^{-\frac{1}{k+1}}\left( w\left( t\right) \right) dt \\
& +\frac{1}{\sqrt{2}}\left[ \int_{R}^{r}t^{1+\varepsilon +\frac{2\left(
k-1\right) }{k+1}}\left( p\left( t\right) \right) ^{\frac{2k}{k+1
}dt+\int_{R}^{r}t^{-1-\varepsilon }dt\right]  \\
& \leq C^{\frac{1}{k+1}}H^{-\frac{1}{k+1}}\left( w\left( R\right) \right)
\int_{R}^{r}t^{1-\frac{N}{k}}dt \\
& +\frac{1}{\sqrt{2}}\left( \frac{k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}
\left[ \int_{R}^{r}t^{1+\varepsilon +\frac{2\left( k-1\right) }{k+1}}\left(
p\left( t\right) \right) ^{\frac{2k}{k+1}}dt+\frac{1}{\varepsilon
R^{\varepsilon }}\right]  \\
& \leq C_{R}\int_{R}^{r}t^{1-\frac{N}{k}}dt+\frac{1}{\sqrt{2}}\left( \frac
k+1}{C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}\int_{R}^{r}t^{1+\varepsilon
\frac{2\left( k-1\right) }{k+1}}\left( p\left( t\right) \right) ^{\frac{2k}
k+1}}d
\end{array}
\label{final}
\end{equation
where
\begin{equation*}
C_{R}=C^{\frac{1}{k+1}}H^{-\frac{1}{k+1}}\left( w\left( R\right) \right)
\frac{1}{\sqrt{2}}\frac{1}{\varepsilon R^{\varepsilon }}\left( \frac{k+1}
C_{N-1}^{k-1}}\right) ^{\frac{1}{k+1}}.
\end{equation*
By taking $r\rightarrow \infty $ in (\ref{final}) we obtain (\ref{13}) since
$w$ is large and $h$ satisfies (C3). These observations completes the proof
of the theorem.

\subparagraph{\textbf{Proof of the Theorem \protect\ref{th3} and \protect\re
{th4}.}}

In order, to obtain the conclusion, combine the proof of \textbf{Theorem \re
{th1} }and\textbf{\ \ref{th2}} with some technical results from \cite{COV}
and \cite{COVS}.

\qed


\section{Introduction}
High intensity lasers ($>$\Wcm{17}) have enabled particle acceleration to high energies (keV to GeV) through a variety of methods and mechanisms, all of which exploit the strong electromagnetic fields that exist during intense laser plasma interactions. Some high intensity laser sources are capable of accelerating electrons to GeV energies through direct laser-field and/or laser wakefield acceleration \cite{esaray2009}. In addition, ultra-intense laser-plasma interactions can create bright X-ray sources from high harmonic generation \cite{Dromey2006}, bremsstrahlung collisions of electrons \cite{murnane1991ultrafast}, betatron oscillations \cite{kneip2010}, Compton scattering \cite{phuoc2012, chen2013} and other techniques. These applications benefit from advances in the peak power that ultra-intense laser systems can provide.

While these experiments are very interesting from a scientific perspective, for laser technology to provide a compelling alternative to conventional accelerators, and for accelerator-based light sources, it is an important goal for the field to produce ultra-intense laser pulses with high \emph{average} power \cite{hooker2013}. This involves creating intense laser pulses at a high repetition rate. Most laser facilities operate with a repetition rate between 0.03-1 Hz \cite{hooker2013}, which is orders of magnitude less than conventional accelerators. Moreover, ultra-intense laser experiments typically operate in a single-shot mode with diagnostics designed for single-shot operation. Next generation laser facilities are soon coming online, with laser systems boasting 10 Hz repetition rate at peak powers exceeding 1 PW \cite{rus2015}. However, these systems would still need high repetition rate target insertion and alignment systems where inexpensive targets can be brought in and aligned within a few ms between laser shots, and experimental diagnostics to extract data from such interactions at the operating duty cycle.

Millijoule-class tabletop short-pulse lasers, with reduced size, complexity, cost, and operating expenses, can provide a useful template for the next generation laser experimental systems to follow. Such a laser experimental system is discussed in this paper, which can operate at a high repetition rate (up to $\sim$1~kHz).  One example of 0.5~kHz laser-plasma electron acceleration is found in \etal{He}\citet{he2013}, with acceleration to 100~keV by laser-wakefield acceleration and demonstrated applications in ultrafast electron diffraction \cite{he2016}. Another high-repetition-rate laser wakefield acceleration study by \etal{Guénot}\citet{guenot2017} used a near-single-cycle laser pulse with temporal contrast better than 10$^{10}$ focused into an underdense flowing gas jet target to accelerate electrons to $\sim$5~MeV. Experiments with these mJ-class laser systems can produce GeV/cm accelerating fields and energetic bursts of electrons with MeV energies \cite{milchberg2015}.
This paper will present direct evidence for some of the most energetic, highest-efficiency MeV electron acceleration that has been observed in the literature to date from laser systems of this kind.

Laser-plasma interaction experiments at and below relativistic threshold intensity (\Wcm{18} for 800~nm light) typically show exponential decays in electron spectra with characteristic energies well below 100~keV (e.g. \cite{mordovanakis2010}), but some studies using dense targets with a controlled pre-plasma have shown higher-energy electron acceleration in oblique-incidence laser-target configurations. \etal{Mordovanakis}\citet{mordovanakis2009} used a 3~mJ laser to accelerate 7~pC of charge to between 100~keV and 1.4~MeV with an energy spectrum peaked at 0.8~MeV, and \etal{Mao}\citet{mao2012} used a 240~mJ laser to accelerate, in at least a 14\degree{} full-angle cone, 100~pC of electrons per shot with similar energy. \etal{Bocoum}\citet{bocoum2016} explored the role of pre-plasma scale length in electron confinement within the target, accelerating electrons with up to 600~keV energy in an estimated 11~pC bunch after careful optimization of pre-plasma scale length ($L\sim0.1\lambda$). In prior studies, the ejection of energetic electrons depends on \textit{p}-polarized laser pulses and large angles of incidence ($\geq$45\degree{}).
 Electrons tend to be emitted near the laser specular and normal directions with laser-to-electron energy conversion efficiency below 0.05\%. The conversion efficiencies in laser-wakefield and dense gas jet experiments likewise tend to fall below this level. A recent study by \etal{Salehi}\citet{salehi2017} demonstrated significant numbers of forward-going $>1$~MeV electrons from 4~mJ to 10~mJ laser interactions with a dense gas jet target. According to their measurements, the laser-to-electron conversion efficiency in those experiments was roughly 0.015\% or less.

In the presently studied case of normal-incidence laser-target interaction, \textit{p}-polarization-dependent electron acceleration mechanisms are suppressed. However, three recently published experimental and computational papers \cite{morrison2015backward, orban2015backward,ngirmang2016} have demonstrated how a wide spray ($\pm$40\degree{} cone) of electrons can be accelerated to $>$120~keV with $\gtrsim$1.5\% laser-to-electron energy conversion efficiency in the laser reflection direction for few-mJ, \Wcm{18}, normal-incidence interactions. A standing wave and direct laser acceleration theoretical framework was developed in \etal{Orban}\citet{orban2015backward} to provide a physical interpretation of this normal-incidence, extremely-high-efficiency electron acceleration by a mJ-class relativistic laser. Simulations in that paper showed $\gtrsim$0.5\% efficiency for electrons accelerated to $>$MeV through the same mechanism, without experimental confirmation.

This paper contributes two scientific advancements to the study of laser-based electron acceleration. First, $>$MeV electrons are experimentally observed in the explicit $-\bm{k}_{\text{laser}}$ (backward) direction for a normal-incidence laser-plasma-interaction, despite a $81$~keV relativistic ponderomotive energy \cite{wilks1992}. Until the present work there has been no direct experimental confirmation of electrons with $>$120~keV energies in experiments of this geometry and directionality. Second, the measured electrons/MeV/s.r. above 1~MeV generated using this low-energy laser (3~mJ) exceeds or equals that of related prior studies employing lasers at oblique incidence of similar intensity. This paper provides the first clear evidence, through direct electron detection in a magnetic spectrometer, that the highly efficient interactions of \etal{Morrison}\citet{morrison2015backward} produce large numbers of $>$MeV electrons.

This work also contributes novel methods for high-acquisition-rate High Energy Density (HED) experiment and normal-incidence laser-target setup. Due to the possibility of damaging delicate optical systems, it is difficult to find studies with ultra-intense laser interactions at normal incidence in the literature, and, perhaps due to the difficulty of obtaining a direct line-of-sight between an electron spectrometer and the mirror opening, it is even harder to find studies with careful investigation of electron acceleration in this configuration.
According to the authors' best knowledge, a 0\degree{} (direct laser back-reflection direction) measurement of super-ponderomotive electron energy in a normal-incidence laser configuration has no precedent in the literature. In this paper, an experimental electron energy characterization is presented which confirms electron acceleration up to 3~MeV in the laser back-reflection direction.   These measurements are performed while taking advantage of the kHz repetition rate of a mJ-class ultra-intense laser system. In the field of HED ultra-intense laser-matter interactions, shot-to-shot variations can frustrate experiments on low-repetition-rate laser systems where, by necessity, researchers are often forced to draw conclusions from only a few shots. A mean spectrum and standard deviation range using fifty electron spectra (each an average of ten shots) is presented, which was adequate for capturing the range of shot-to-shot variations in this experiment.
This energy characterization is discussed with new 3D Particle-in-cell (PIC) simulations and in a theoretical framework of normal-incidence, laser reflection super-ponderomotive electron acceleration at relativistic threshold intensity.

\section{Experimental Setup}

\begin{figure}[htbp]
    \centering
    \includegraphics[width=8cm]{holey9}
    \caption{\label{fig:setup}To measure the electron spectrum in the direction anti-parallel to the laser axis (back-reflection direction), a 3~mm hole was drilled into the 25.4~mm diameter 30\degree{} gold-coated f/1.1 OAP. The central portion of specular-direction electrons, accelerated by normal-incidence laser interactions with a flowing water target, pass through the hole. Electron energies are characterized by a magnetic spectrometer located 127~mm behind the target. The laser is polarized in the plane of the diagram.}

\end{figure}

The target area and experimental setup for this experiment are sketched in Fig. \ref{fig:setup}. A 1~kHz Ti:sapphire laser system delivers 4.8~mJ (2.9~mJ on target) with a pulse duration of 40~fs full-width at half-maximum (FWHM). The laser is focused onto a high-velocity water stream (a column 30~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} in diameter) \cite{morrison2015backward}, which enables the study of laser-plasma interactions at 1~kHz. The laser is focused using an off-axis paraboloid (OAP). In order to make the desired normal-incidence $-\bm{k}_{\text{laser}}$ (backward) experimental measurement, a 3~mm hole was drilled into its center. As shown in Fig.~\ref{fig:setup}, this hole allowed a portion of electrons accelerated from the target in a 6.3\degree{} full-angle cone in the laser reflection direction to pass through to an electron spectrometer placed 127~mm behind the target.

Modifying the OAP in this way to obtain such a measurement is atypical in experimental high-intensity laser-plasma interaction physics because it severely distorts the laser focus, reducing on-target intensity. An adaptive optics system, similar to the one pioneered by \etal{Albert}\citet{albert2000}, was custom designed and configured to compensate for intensity degradation due to the hole in the OAP. The deformable mirror was inserted into the laser chain after the laser compressor. The in-vacuum focal spot area after the OAP was measured using a microscope objective imaging system and CCD camera. The imaged laser focal spot area and quality was analyzed in real-time to close a feedback loop (at 60~Hz) with deformable mirror actuators through a genetic algorithm (written by the authors in LabVIEW).
Despite the energy loss and decreased optical quality due to the hole in the OAP, the adaptive optics setup converged within five minutes to a 2.2~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} FWHM in-vacuum focal spot, with 2~mJ of the 2.9~mJ on-target energy contained in the Gaussian laser focus (see Fig. \ref{fig:spots}(a)). The pulse energy stability was $3\%$ peak to peak with $<1\%$ rms.
Accounting for the hole and adaptive optics corrections, the total on-target laser energy in this experiment was 2.9~mJ and laser intensity was \Wcm[1.5]{18}.

The laser's intrinsic pre-pulse, which extends over several nanoseconds, is experimentally manipulated using electro-optic gating.  Pockels cell timing and extinction ratio are adjusted at four points throughout the laser amplification chain to reach the condition of high-number/high-energy electron acceleration presented in this paper. These adjustments are made in response to electron spectra that are acquired, analyzed and displayed in real-time. As a consequence of $\sim$3~ns rise time and temporal jitter of the Pockels cells, combined with the unsaturated gain of the pre-pulse in the amplifiers, the photodiode-measured pre-pulse energy stability was 15\% rms. Importantly, the optimal conditions do \emph{not} correspond to the highest level of pre-pulse suppression that the system is capable of.
Nanosecond-timescale laser contrast was measured using a Thorlabs DET10A photodiode and a 2.5~GHz oscilloscope. Dynamic range in the measurement was enhanced by recording and recombining multiple oscilloscope measurements, each with different levels of neutral-density optical filtering on the photodiode \cite{morrison2015backward}.  Due to the oscillator pulse-train and ghost reflections within the multi-pass laser amplifiers,  40-fs-duration short-pulse replicas are expected at fixed delay times in the nanoseconds preceding the main pulse. Evidence for these short-pulse replicas can be inferred from the nanosecond signal, and deconvolved using the impulse response of the photodiode/oscilloscope system. A scanning third-order cross-correlator measurement of the picosecond-timescale laser contrast confirms this conclusion with greater temporal precision.
The laser amplified spontaneous emission (ASE) contrast for this experiment is measured to be $5~\cdot~10^{-6}$ at {$t=-15$~ps} with a 40-fs-duration short pulse replica at $t=-6.2$~ns containing a shot-to-shot averaged 70~\textmu J} %% defined to make microjoules easier to write. Invoke via \uJ{{} at \Wcm{16}. A plot of the average laser pre-pulse profile is shown in Fig. \ref{fig:ppulse}(a). As mentioned earlier, these experiments could have been performed with more aggressive suppression of the $t=-6.2$~ns short pulse replica and the ASE. Instead, a relatively low-contrast condition was adopted because this empirically maximizes the number of high-energy electrons accelerated in the laser reflection direction.

To characterize the effect of the pre-pulse on the condition of the target, a hydrodynamic simulation of pre-plasma evolution was performed using the FLASH code \cite{flash2000, flash2009, flash2012}. The simulation was performed in 2D R-Z cylindrical geometry.
Because of the 2D cylindrical geometry, the water column target is represented as a water droplet of radius 15~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}. The simulation extent is $R=\pm$28~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} and $Z=$-85~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} to $Z=+$35~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}, with $R=Z=0$ at the front center surface of the water droplet. Cell spatial resolution is 0.25~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}. The timestep varies during the simulation according to the Courant limit, which is determined by the fastest sound speed in the simulation and the cell size. The laser focus convergence with the experimental f/1.1 and focal depth is also set to match the experiment (2.2~um{} Gaussian focus located 15~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} beneath the target surface). Because highly intense femtosecond-timescale laser-plasma-interaction physics can not be accurately modeled in a hydrodynamic simulation, the 40~fs replica pulses at $t$=-6.2~ns and $t$=-15.3~ps are excluded from the pre-pulse profile. Results of this simulation show that even without these features, a pre-plasma with tens-of-microns extent develops (Fig.~\ref{fig:ppulse}(b)).

\begin{figure}[htbp]
    \centering
    \includegraphics[width=8cm]{spots3}
    \caption{\label{fig:spots}
    a) The focal spot measured after deformable mirror optimization, with a calculated \Wcm[1.5]{18} intensity. b) Portion of the electron beam that passes through the hole of the OAP, imaged with a Lanex phosphor camera \cite{orban2015backward}. The black rectangle denotes the region (1~mm wide x 1.7~mm tall) that enters the magnetic spectrometer's slit and is energy-resolved. In this image, a diagnostic 250~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} thick stainless steel mesh sits in front of the camera to check spectral uniformity; both camera and mesh are removed before electron spectra measurements.
    }

    \end{figure}

\begin{figure}[htbp]
    \centering
    \includegraphics[width=13cm]{ppulse6}
    \caption{\label{fig:ppulse}a) Measured characteristic pre-pulse intensity on the nanosecond and picosecond timescales. The spikes in intensity at $t$~=~-6.2~ns and $t$~=~-15.3~ps are 40-fs-duration main-pulse replicas. The main pulse arrives at t~=~0~ps. Light blue segments of the picosecond trace have been modified to remove third-order cross-correlation artifacts. b) Pre-plasma density profile from the long-pulse component only of (a) as modeled in FLASH. Ablated material extends tens of microns from the target, with a steeper profile nearer to solid density (1~g/cc). Original target interface is at Z=0~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}. The density profile can be described by multiple exponential decays, indicated with their characteristic scale lengths; long scale lengths exist far from target. The 800~nm plasma critical density is indicated assuming eight free electrons per H$_2$O molecule (a typical value derived from the PIC simulations, due to ionization by the main pulse). c) Reference interferometric image for (d) and (e), with the laser blocked. d) Femtosecond-gated 420~nm probe beam interferometry captures some information about the target condition 1~ps before main pulse ionization. Fringe shifts are evident up to several microns away from the target, demonstrating the presence of pre-ablated material. According to the fringe shift direction, ionization is incomplete at this stage, which may affect the sensitivity of the interferometry to lower-density material further from the target. e) Interferometry delayed to 5~ps after main pulse ionization. The direction of the fringe shift indicates that the material is in the ionized state. Fringe shifts at tens of microns away from the target indicate the presence of material at this distance. Due to the brief time interval between the two interferometric measurements this is interpreted as diffuse material that was present before the ultra-intense laser interaction, but undetectable in interferometry due to its low density and partially ionized state.}

\end{figure}

Laterally sheared optical interferometry is used in Fig. \ref{fig:ppulse}(d) and \ref{fig:ppulse}e to show target state near the time of interaction. To maintain order-femtoseconds synchronicity between pump and the probe pulses, seed pulses from the same ultrafast oscillator are used for generating pump and probe. Separate amplifier systems are used. The amplified probe beam is frequency-shifted and frequency doubled, resulting in 420~nm probe which is optically isolated from the on-target main pulse second harmonic light (390~nm) using a narrow-linewidth optical filter (see \etal{Feister}\citet{feister2014}). The probe beam passes through the interaction region, as shown in Fig. \ref{fig:setup}, then enters a microscope objective for imaging purposes. The beam path is then is split for separate shadowgraphic and interferometric analysis. The interferometric beam path undergoes lateral shearing such that the vacuum region of the beam image is overlapped with the interaction region to create interference fringes, which change in angle based on integrated index of refraction differences between the overlapped regions.

A comprehensive characterization of electron energies in the laser-plasma-interaction was desired, and stochastic shot-to-shot fluctuations are a typical feature of ultra-high intensity laser plasma interactions with near solid density targets. To capture stochastic shot-to-shot fluctuations in accelerated electron energy, a high acquisition rate electron spectrometer was designed, built, and calibrated for this experiment. The spectrometer consists of a 1~mm slit followed by two yoked magnets to deflect electrons according to their energy, and encased in a steel enclosure. The electrons are imaged using a Tb-doped $\mathrm{Gd_2 O_2 S}$ phosphor screen (Lanex regular) that is placed directly upon a bare linear CCD chip (MighTex CCD camera, Toshiba 1304DG linear CCD) and made light-tight with a 25~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} thick covering of aluminum. An absolute calibration of this detector, in electrons per CCD count (adjusted for CCD bias) at 1~MeV, 2.0~MeV, and 2.5~MeV, was obtained using the University of Notre Dame Radiation Laboratory Van de Graaff 3~MeV electron accelerator.
Considerations were taken of the continuous-source nature of the Van de Graaff electron beam and the discretely pulsed nature of the experimental electron beam. The Lanex regular phosphor response is linear with electron number at these fluences, with a 1/e light decay time of 0.6 ms \cite{Nikl2006}. For a continuous source, the total light emitted from the screen during a single CCD exposure of fixed, arbitrarily short duration will include some decaying phosphorescent emission of electrons incident prior to CCD exposure, but this light will integrate equally to the missed phosphorescent decay of electrons incident during CCD exposure. Although there are shot-to-shot variations in the laser experiment, the same argument holds true for the 10~ms experimental exposures that were used to gather data for Figs.~\ref{fig:espec}-\ref{fig:specscan}. Due to the kHz repetition rate these 10~ms exposures capture integrated light from ten consecutive electron pulses.
Because the exposure time is much longer than the decay time, a shot occurring prior to exposure is unlikely to affect the recorded exposure by more than $\sim$2\%, even with 100\% deviation from the mean. Moreover, this effect should neither bias the measured fluence up or down, as mentioned earlier.
An MCNP model of the calibration setup at the Van de Graaff accelerator was used to correct for electron scattering within the beam's aluminum exit window from vacuum, the at-air environment of the spectrometer, and the Van-de-Graaff-facility Faraday cup. The absolute experimental result at these discrete energies, after corrections from MCNP modeling of the effects just mentioned, were used to determine the electron spectrometer response characteristics applied to generate the experimental electron spectra of Figs.~\ref{fig:espec} and \ref{fig:specscan}. The electron spectrometer's electron energy positions (at the detector) were verified in two ways. First, Hall probe magnetic field measurements were combined with relativistic particle deflection calculations. Second, energy positions were verified between 1~MeV and 2.5~MeV using the Van de Graaff accelerator.

Before performing energy characterization of back-reflection direction electrons, the laser focus was aligned to the target as follows: the cylindrical water stream was vertically translated such that the laser was incident 1~mm below the nozzle output, vertically angled for normal incidence (as seen in shadowgraphy/interferometry, through air breakdown by the laser above 150~Torr), horizontally translated for normal incidence (as seen by optical back-reflection), and focally translated for strongest backward-going electron signal on the electron spectrometer. The in-vacuum laser pointing and target position stability was measured to be $\sim$1~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}. The electron spectrometer is located 127~mm behind the target (100~mm behind the OAP). Before measurements, a separate permanent magnet is swung into place between the OAP and electron slit (magnetically but not physically blocking the electrons' paths) to check that recorded signal drops to zero. This ensures the measured signal is from electrons and not due to stray X-rays or laser light. The spectrometer is triggered and acquires electron traces at $\sim$100~Hz.

Figure~\ref{fig:espec} shows experimental and 3D Particle-in-cell (PIC) simulation results for laser-reflection direction electron acceleration. These 3D PIC simulations of femtosecond timescale laser-matter interactions were performed using the Large-Scale Plasma (LSP) code  \cite{Welch2004}. The incident laser parameters were similar to that of the highest energy experimental results shown: intensity \Wcm[1.5]{18}, laser wavelength $\lambda = 800$~nm, Gaussian spot diameter 2.3~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}, and 30~fs FWHM pulse duration. The simulation space was 35~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} along the longitudinal ($k$-vector) direction, 40~\textmu m} %% defined to make microns easier to write. Invoke via \um{{}~$\times~$40~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} in the transverse directions, and the cell size for these simulations was $\lambda /8 = 100$~nm with 27 macro-particles per cell per species (electrons, Oxygen, and Hydrogen ions). The density profile was assumed to be exponential with a 1.5~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} scale length. Due to computational limitations, the long-scale-length low-density component of the FLASH simulations was excluded. The simulations were run for more than 100~fs after the peak of the pulse reached the critical density surface of the pre-plasma. Backward-going electron spectra were obtained by tracking electron macroparticles and binning the charges and energies of these particles exiting the simulation space.

\section{Results and Discussion}
\begin{figure*}[htbp]
    \centering
    \includegraphics[width=13cm]{specmain4}
    \caption{\label{fig:espec}
    Energy characterization of electrons accelerated in the laser reflection direction for normal-incidence laser interaction. The y-axis (charge density) for the experimental data is absolutely calibrated and can be directly compared to the 3D PIC simulation results shown. This calibration also includes an adjustment for detector sensitivity, number of shots per spectrum, and electron scattering losses in the 20~Torr chamber pressure. The black curve shows the mean of fifty consecutive experimental electron spectra in which the in-vacuum laser intensity was measured to be \Wcm[1.5]{18}. The standard deviation range of the measured spectra brackets the curve in gray. 3D PIC electron spectra are plotted for two intensities, \Wcm[1.5]{18} and \Wcm[5.4]{17}, which highlights the rapid turn-on with intensity of this acceleration mechanism. For both PIC simulations, the pre-plasma density profile is the same and only electrons emitted in the laser reflection direction within a 6.3\degree{} full-angle cone are included.
    }
    \end{figure*}

\begin{figure*}[htbp]
    \centering
    \includegraphics[width=13cm]{multishot3}
    \caption{\label{fig:multishot} A waterfall plot presents fifty consecutive experimental electron spectra from \Wcm[1.5]{18} laser interactions used to compute the mean electron spectrum of Fig. \ref{fig:espec}. These measurements are enabled by the union of a high-repetition-rate laser, liquid targets, electronic detectors, and an automated acquisition system. Each spectrum, represented by an individual stream in the plot, is the result of a ten-millisecond CCD exposure (averaging over ten laser-target shots). This plot in total represents 0.5 seconds of continuous data acquisition at a 100~Hz rate. Large variations in the electron spectra are partly due to fluctuations in the laser pre-pulse and the sensitivity of the ultra-intense laser-matter interactions with near solid density targets to the properties of the pre-plasma.
    }
    \end{figure*}

\begin{figure*}[htbp]
    \centering
    \includegraphics[width=13cm]{specscan5}
    \caption{\label{fig:specscan} Experimental electron energy characterization for three laser intensities shows extremely steep intensity scaling. Intensity is reduced through use of a variable waveplate and polarizer, meaning the main pulse and its pre-pulse are equally attenuated. Fifty to one hundred consecutive electron spectra are acquired at each laser intensity, and a mean electron spectrum (solid line) as well as the standard deviation range (1$\sigma$ range, filled area) are plotted. As in Fig. \ref{fig:espec}, the experimental spectra are adjusted for detector sensitivity, number of shots per spectrum, and electron scattering losses due to the 20~Torr chamber pressure.}
    \end{figure*}

The data shown in Fig.~\ref{fig:espec} represents the synthesis of high repetition rate targetry, a high repetition rate intense laser and high acquisition rate (100~Hz) spectrometery. There are several features of Fig.~\ref{fig:espec} that merit discussion. First, due to the low repetition rate of most ultra-intense laser systems \cite{hooker2013},
it is relatively rare in studies of ultra-intense lasers interacting with dense targets to see plots with both a mean and standard deviation for electron spectra.  Each shaded infinitesimal area slice between a particular electron energy E\textsubscript{0}, and E\textsubscript{0} + $\Delta$E in Fig.~\ref{fig:espec} is determined from 50 spectral acquisitions (500~laser-interactions with the target). This capability is an important advantage of kHz ultra-intense laser systems because, as shown in Fig.~\ref{fig:multishot}, one finds large variations in the electron spectrum in laser plasma interactions with near solid targets. A major contribution to electron spectrum fluctuations is the sensitivity of ultra-intense laser interactions with solid-density targets to the properties of pre-pulse-generated pre-plasma \cite{bocoum2016}.  Few-shot experiments with low repetition rate ultra-intense laser systems should be similarly sensitive to these details. Secondly, the two higher show electrons accelerated to well above the ponderomotive scaling energy (81~keV for \Wcm[1.5]{18}) in large numbers. Results like these, which are insightful for understanding the limits of plasma acceleration processes,  may go unnoticed by few-shot experiments.

Third, compared with electron acceleration studies at oblique laser incidence, the measurements of  electrons/MeV/s.r. above 1.3~MeV for \Wcm[1.5]{18} exceed all previous observations in the literature for a few-mJ laser system. While there is tension between this measurement and the electron spectrum for the 3D PIC simulation, with the experiment indicating more electrons than expected from simulations, both the simulations and experiment qualitatively agree that large numbers of MeV electrons are ejected from the target. Figure~\ref{fig:espec} provides the first direct evidence that significant numbers of electrons can be accelerated to super-ponderomotive energies in water jet experiments with $\sim$\Wcm{18} intensities at near normal incidence. Although the experimental conditions in \etal{Morrison}\citet{morrison2015backward} do differ in some ways from the experimental conditions examined here, Fig.~\ref{fig:espec} also confirms an indirect measurement of the ejected electron energies from \etal{Morrison}\citet{morrison2015backward}.

Figure~\ref{fig:specscan} demonstrates the intensity sensitivity of the electron acceleration in these experiments. Earlier, it was shown in Fig.~\ref{fig:espec} that a 3D PIC simulation at \Wcm[5.4]{17} resulted in orders of magnitude lower numbers of ejected electrons despite laser intensity lowered by only about a factor of two from that of the previous one at \Wcm[1.5]{18}. In Fig.~\ref{fig:specscan} it is likewise found that, experimentally, the efficiency of the electron acceleration process is highly sensitive to the peak intensity in this range.

An important aspect of the experiments presented here is that to avoid freezing the water jet target, the target chamber pressure could only be reduced to 20~Torr. A natural question is whether this affected the properties of the ultra-intense laser pulse prior to reaching the target. An empirical test was performed without a water jet target to examine this possibility. The laser was allowed to propagate through background pressures between 1~Torr and 100~Torr, and the spectrum and beam mode after going through focus is observed. No spectral shift was apparent below 40~Torr, and no beam mode degradation was observed below 80~Torr, indicating that self-focusing was insignificant at the 20~Torr operating pressure of the water jet experiments presented here.
\textit{Geant4} \cite{agostinelli2003} modeling was also performed to assess whether electrons (especially at lower energies) would be deflected out of the electron spectrometer's acceptance angle due to scattering in ambient water vapor in the vacuum chamber.
The \textit{Geant4} simulations indicate that the 20~Torr ambient pressure in the target chamber diverts $\sim$22\% of $\sim$1~MeV electrons and 9\% of 2~MeV electrons from entering the spectrometer. These corrections were applied to the electron spectra plotted in Figs. \ref{fig:espec} and \ref{fig:specscan}.

A theoretical basis for the super-ponderomotive energies of backward-going electron acceleration for this experiment is described in \etal{Orban}\citet{orban2015backward}. A standing wave is established between the forward-going laser and its backward reflection. As a result, the peak electric and magnetic fields significantly exceed the values present in the incident laser pulse due to constructive interference.
Unlike oblique incidence standing waves, as described for example by \etal{Mordovanakis}\citet{mordovanakis2009}, normal incidence standing waves have stable node and anti-node positions, which enhances electron acceleration by developing quasi-static charge separations in a process known as ponderomotive steepening \cite{estabrook1983}. The standing wave provides a pathway for electron acceleration both into and away from the target \cite{kemp2009,orban2015backward}, a process that can be influenced by quasi-static fields that arise from the charge separation caused by the standing wave. Importantly, electrons accelerated into the target will continue without another strong interaction with the laser field, while those accelerated away from the target (laser back-reflection direction) can be injected into the reflected light and undergo direct laser acceleration to higher energies. This direct laser acceleration of ejected electrons by reflected laser light provides a compelling explanation for the MeV electron energies detected experimentally in the current study. The appendix of \etal{Orban}\citet{orban2015backward} also provides an argument for why this mechanism should be much less efficient for laser intensities significantly below \Wcm{18}.

As has been discussed, the measured values for the number of electrons/MeV/s.r. are large compared to prior literature using similar laser energies and intensities. The measured values are also large compared with the 3D PIC simulation presented earlier, and it is argued here that this is due to the limitations of the pre-plasma extent adopted in the simulation. In particular, a multi-exponential-scale-length pre-plasma was not included in the LSP simulation. Instead, a single 1.5~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} scale length was assumed which is similar to the $\sim$1.3~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} exponential scale length that the FLASH simulation results in Fig.~\ref{fig:ppulse}(b) indicate is present near the critical density. Since many prior studies have emphasized the importance of the characteristic exponential scale length of the pre-plasma near the critical density this is an important aspect of the laser-matter interactions to include.

Although the FLASH results also indicate a $\gtrsim 5$~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} scale length extending from the target and indications from interferometry also point to the existence of an extended pre-plasma, 3D PIC simulations with a more extended pre-plasma were not presently possible due to computational expense. This affected the results of the PIC simulations in important ways. It is known from performing 2D PIC simulations using pre-plasmas that are similar to the FLASH density profile in Fig.~\ref{fig:ppulse}(b) and similar laser pulse parameters that the presence or absence of a low density, long-scale length component of the pre-plasma can have a strong effect on the number of MeV back-directed electrons. In a test by the authors, the addition of a diffuse (sub-sixth-critical-density), long-scale length pre-plasma with a scale length of 7.2~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} increased the conversion efficiency of back-directed $>120$~keV electrons from 0.67\% to as much as 3.4\%, an increase by a factor of five\cite{feisterDPP2016}. This dramatic increase in efficiency can be explained by the "capacitor" model developed by \etal{Link}\citet{link2011}, which considers in detail the process of electron ejection from short-pulse laser illuminated targets that are charging up as energetic electrons depart. As these electrons leave the electric potential on the surface of the conducting plasma target increases, which strongly affects the energy threshold for escape and the final energy of escaped electrons. \etal{Link}\citet{link2011} also considers the effect of target size, concluding "More electrons escape for internal electron distributions with more energy, for electron distributions with higher hot electron temperatures, and for larger targets." The simple reason for this being that larger conductors naturally have a lower electric potential, and due to the generically steep decrease in the numbers of energetic electrons with increasing energy in the target, a lower electric potential will allow many more electrons to escape. \etal{Link}\citet{link2011} discusses this in the context of multi-Joule laser systems, but the essential plasma physics is the same in describing mJ-class experiments. It should be emphasized that the size of the target depends on the spatial extent of the pre-plasma, and that 2D PIC simulations confirm the essential reasoning that increasing the target size in this way can lead to increases in the number of MeV scale electrons by a large factor. Our 3D PIC simulations therefore describe the ejection of electrons via a standing-wave-acceleration process, but the details of electrostatic turnaround and charging of the target are mis-represented in a way that reduces the energy and, importantly, number of MeV scale electrons in important ways. Consequently, the 3D PIC simulations confirm the existence of significant numbers of MeV scale electrons but the simulations underpredict the measured numbers of ejected electrons by a large factor.

It is interesting to consider the properties of the 6.3\degree{} full-angle divergence sub-beam (shown in Fig.~\ref{fig:spots}(b)) which is created by clipping through the OAP hole.
These estimates do not include the electrons in the sub-beam with energies below 1~MeV.
Integrating the fifty electron spectral curves (for \Wcm[1.5]{18}) from 1~MeV to 3~MeV and over the solid angle 9.6~msr of the OAP hole, one finds in the 500-shot-mean an acceleration of 6.9~pC/shot and 0.33\%~efficiency into $>$MeV electrons in the sub-beam. Since the laser-plasma interaction occurs continuously at 1~kHz, this represents a time-averaged 6.9~nA current of $\gtrsim$MeV electron bursts powered by a 3~W average-power laser. The liquid target recovers in tens of microseconds \cite{feister2014}, so a higher-repetition-rate laser could be expected to generate linearly higher currents. The highest-efficiency spectrum indicates a ten-shot-mean of 16.4~pC/shot and 0.83\% efficiency into the sub-beam. It is worth noting that the laser-specular electron spray is clipped by the OAP hole to create this sub-beam, and therefore the total accelerated $>$MeV charge/shot into the \emph{full} electron spray is almost certainly higher than these quoted numbers (consistent with the observations of \etal{Morrison}\citet{morrison2015backward}).

From an applications point-of-view, two notes on this work should be made. First, this setup creates a high peak current density of MeV electrons with a large \emph{average} current using a relatively inexpensive target setup in a rough vacuum and without stringent constraints on the laser pre-pulse contrast. This demonstration opens the door to a ultra-intense laser-based compact MeV gamma ray generator. In fact, using this setup, but without a hole in the OAP, one can easily generate highly directional rad/hr hard x-ray and $\gamma$-ray radiation from electron scattering with the OAP \cite{morrison2015backward}. Since the electrons are accelerated on fs-timescales and the fluence is high, this source could be used as a compact electron or X-ray backlighter to probe the temporal evolution of various systems of interest with a kHz cadence.
Second, the shot-to-shot fluctuations could be significantly reduced in future studies and in applications through modifying the experimental setup. Whereas in the current study the pre-plasma was created by the deliberate amplification and shaping of the natural laser ``noise'' of ASE and spurious short-pulse-replicas, one could instead introduce a separately amplified pulse in order to reduce shot-to-shot variation.
Another source of shot-to-shot variation is due to the convex curvature of the water column target. Because of this curvature, $\sim$1~\textmu m} %% defined to make microns easier to write. Invoke via \um{{} transverse jitter in the laser focal spot position on the target changes the shot-to-shot effective depth of laser focus. This source of variation could be suppressed by changing the geometry of the flowing water from cylinder-like to sheet-like, and by widening the on-target pre-pulse diameter by defocusing any injected laser pre-pulse beam to increase transverse pre-plasma uniformity.

\section{Conclusion}
A direct measurement of kHz-pulsed super-ponderomotive electron acceleration in the back-reflection direction from normal incidence, 3~mJ laser interactions with solid-density liquid water at \Wcm[1.5]{18} has been performed and presented. The high-repetition-rate and low-laser-energy of this MeV~electron experiment is promising for electron-beam applications, and the results are scientifically novel in two ways: 1) the experiments were performed at normal incidence and the super-ponderomotive electrons were back-directed, details that required substantial modification of the experimental setup in order to study, and 2) the number of electrons/s.r./mJ (and electrons/s.r./minute) above 1.3~MeV exceeds all other studies for this laser energy and intensity. These measurements also directly confirm, for the first time, expectations from 2D and 3D PIC simulations that MeV electrons are indeed present in these interactions \cite{orban2015backward,ngirmang2016}. Since the experiment occurs at normal incidence and the ejection direction of super-ponderomotive electrons is laser specular, the mechanism is distinct from that of oblique-incidence electron acceleration. The number of electrons/MeV/s.r. detected at these energies also reinforces earlier conclusions that these ultra-intense laser interactions occur with a remarkably high conversion efficiency from laser energy to energetic electron energy \cite{morrison2015backward}. Methodologically, this work provides an unprecedented template for high-average-power HED experiment combining liquid-density targets, high acquisition rate plasma characterization, and large-sample statistical data analysis.
Future work will involve replacing the ASE/short-pulse pre-pulse with a controlled second-beam pre-pulse, co-timing acquisition from a variety of high-acquisition-rate laser/plasma diagnostics for statistical interrogations, and optimization of the electron beam divergence using adaptive laser optics.

\section*{Funding}

Air Force Office of Scientific Research (AFOSR). The  authors  acknowledge significant support from the Department of Defense High  Performance  Computing  Modernization  Program (DOD  HPCMP)  Internship  Program.

\section*{Acknowledgments}
This research was sponsored by the Quantum and Non-Equilibrium Processes Division of the Air Force Office of Scientific Research, under the management of Dr. Enrique Parra, Program Manager. Support was also provided by the DOD HPCMP Internship Program. The Notre Dame Radiation Laboratory and J.A. L. are supported by the Division of Chemical Sciences, Geosciences and Biosciences, Basic Energy Sciences, Office of Science, United States Department of Energy through grant number DE-FC02-04ER15533. This contribution is NDRL-6000 from the Notre Dame Radiation Laboratory. ABW: 88ABW-2015-3607.

\newpage
\section*{References}
\bibliographystyle{unsrtnat}


\section{Introduction}
Carl Seyfert noted in 1943, the existence of a small fraction of spiral galaxies that had bright star-like nuclei and peculiar emission line spectra populated by broad and narrow high-ionisation emission lines. We now know that these Seyfert galaxies are a subset of the class of sources hosting active galactic nuclei (AGN). The extreme luminosities of AGN, which can exceed the total light from the host galaxy by a factor of ten to a thousand, are understood to be a consequence of the release of gravitational energy as matter gets accreted onto a supermassive black hole (mass $\sim10^6-10^9$\,M$_\odot$). The emission lines are produced in dense fast-moving gas clouds (velocities $\sim100-10,000$\,km\,s$^{-1}$) that orbit around the black hole-accretion disk system, inside what are known as the broad and narrow line regions (BLR and NLR).
A small fraction of AGN exhibit double peaks in their emission line spectra (e.g., Figure~1), indicative of a different kind of activity in their central regions.

Three popular scenarios have been put forward to explain double peaked emission line AGN (DPAGN): (i) their emission line regions are in rotating disks, (ii) their emission line clouds are being pushed out by radio outflows originating in the central engine, and (iii) there are binary black holes in the AGN carrying their respective emission line regions with them \citep[e.g.,][]{Gaskell83,Rosario10,Smith12,Gabanyi14}. However, binary black holes are expected to be rare in spiral galaxies that typically host Seyfert nuclei. The hierarchical galaxy formation model \citep[e.g.,][]{Steinmetz02} proposes that spiral galaxies merge to form elliptical galaxies, while they themselves only grow through minor mergers which may or may not result in binary black holes \citep{Aguerri01}. $<$1\% of DPAGN are likely to have binary black holes with kiloparsec-scale separations in complete samples of ``radio-quiet'' AGN \citep{Rosario11}. The search for binary black holes with parsec-scale separations requires the technique of very long baseline interferometry (VLBI). Here we briefly describe the VLBI observations of two double peaked emission line Seyfert galaxies residing in spiral hosts.

\section{Seyfert Galaxies with Double Peaked Emission Lines}
We identified six DPAGN in the KPNO Internal Spectroscopic Survey Red \citep[KISSR;][]{Wegner03} comprising overall of 72 Seyfert and Low-Ionization Nuclear Emission-Line Region (LINER) galaxies, through their SDSS\footnote{Sloan Digital Sky Survey} optical spectra. The fraction of DPAGN in the KISSR Seyfert+LINER sample is therefore $\approx$8\%. However only three of these six Seyfert galaxies were detected in the VLA FIRST\footnote{Faint Images of the Radio Sky at Twenty-Centimeters} and NVSS\footnote{The NRAO VLA Sky Survey} surveys, and were suitable targets for carrying out a phase-referenced VLBI study.
Results from the VLBI study of KISSR\,1494 and KISSR\,1219 are presented below.

\noindent
{\bf KISSR\,1494:}
VLBI observations of the Seyfert galaxy, KISSR\,1494 ($z=0.0574$), revealed a single radio core of size $\sim7.5\times5$~milliarcsecond (=$8\times6$ parsec) with a flux density of $\sim650\,\mu$Jy at 1.6~GHz \citep{Kharb15}. There was no core detection at 5 GHz. The core brightness temperature ($T_b$) of $\sim1.4\times10^7$\,K was too high to be from a nuclear starburst which typically have $T_b<10^5$\,K, while the $1.6-5$ GHz radio spectral index of $-1.5\pm0.5$ was too steep to be from thermal free-free emission. These parameters are similar to those derived in ``normal'' Seyfert galaxies (ones with single peaked emission lines) that have parsec-scale radio outflows \citep[e.g.,][]{Kharb10,Kharb14}. The steep radio spectrum was also inconsistent with the KISSR\,1494 core being the unresolved base of a relativistic jet, which is expected to have a flat or inverted spectrum. We therefore concluded in \citet{Kharb15} that the radio emission was non-thermal synchrotron in origin that was likely coming from the base of an extended outflowing coronal wind. The NLR clouds were probably being pushed out by this broad nuclear outflow, giving rise to the double peaks in the optical spectrum.

\noindent
{\bf KISSR\,1219:}
VLBI observations of the second Seyfert galaxy, KISSR\,1219 ($z=0.0375$), took place in February 2015. A clear 100 milliarcsecond ($\sim$70 parsec) one-sided core-jet structure is revealed in our preliminary 1.6~GHz image (Kharb et al 2016, in preparation). The core but not the jet is detected at 5 GHz, indicating that the jet has a steep radio spectrum. Interestingly, this parsec-scale jet is pointing directly towards some lobe-like radio emission about 50 kiloparsec away, as seen in its VLA FIRST image (see Figure 2). If the parsec-scale and the kiloparsec-scale emission are connected, KISSR\,1219 will be the fourth only Seyfert galaxy with a large-scale jet ejected from a spiral galaxy \citep[see][]{Bagchi14}. New Expanded VLA observations of KISSR\,1219 are currently being analysed to search for this connection. Like in KISSR\,1494, there is an indication that the double peaked emission lines in KISSR\,1219 are a result of an AGN outflow impacting the emission line clouds.

\section{Summary}
Double peaked emission lines could indicate the presence of binary black holes, rotating BLR/NLR disks, or AGN outflows impacting emission-line clouds. Parsec-scale outflows or binary black holes with parsec-scale separations can be searched through the technique of VLBI. In this paper, we briefly describe results from phase-referenced VLBI observations at 1.6 and 5 GHz of two DPAGN from the KISSR Seyfert+LINER sample, residing in spiral host galaxies.

A single slightly resolved radio core was detected at 1.6~GHz in the Seyfert galaxy KISSR\,1494. The high brightness temperature and the steep $1.6-5$ GHz spectral index were consistent with the core being the base of an extended outflowing synchrotron-emitting coronal wind, rather than the unresolved compact base of a radio jet \citep{Kharb15}. A clear $\sim$70 parsec radio jet is revealed, on the other hand, in the 1.6~GHz VLBI image of the Seyfert galaxy KISSR\,1219.
The VLBI observations therefore do not reveal any signatures of binary black holes in the two Seyfert galaxies, but rather those of an AGN-related outflow. We conclude that an AGN outflow$-$emission line cloud interaction is the most likely cause of the splitting of emission line peaks in these Seyfert galaxies.

\noindent
{\bf Acknowledgements}
We thank the referee for helpful suggestions. PK would like to thank the RETCO-II meeting organisers and Dr. Indranil Chattopadhyay for this invited talk. ZP was supported by the International Space Science Institute in Bern.

\bibliographystyle{apj}