\section{Introduction}

	Let $n$ and $k$ be two positive integers with $n\ge k$.
	For an $n$-element set $X$, denote the set of all $k$-subsets of $X$ by $\binom{X}{k}$.
	Given a positive integer $t$, we say a family $\mf\subset\binom{X}{k}$ is \emph{$t$-intersecting} if $|A\cap B|\ge t$ for any $A,B\in\mf$.

	A $t$-intersecting family is called \emph{trivial} if every element of this family contains a fixed $t$-subset of $X$.
	The famous Erd\H{o}s-Ko-Rado theorem \cite{EKR,Fn,Wn} states that a $t$-intersecting family with maximum size is trivial when $n>(t+1)(k-t+1)$.
	The structure of non-trivial $t$-intersecting families with maximum size was completely determined in 1996 \cite{SHMT,SHM}. Recently, some other maximal non-trivial $t$-intersecting families were described in \cite{LSET2,LSET3,LSET1}.

	Let $V$ be an $n$-dimensional vector space over the finite field $\ff$, where $q$ is a prime power, and ${V\b k}_q$ denote the family of all $k$-dimensional subspaces of $V$.
	We usually replace ``$k$-dimensional subspace" with ``$k$-subspace" for short.
	Let $\mf\subset{V\b k}_q$ be a family with $\dim(F_1\cap F_2)\ge t$ for any $F_1,F_2\in\mf$, which we say is \emph{$t$-intersecting}. In particular, $\mf$ is called \emph{trivial} if there exists a $t$-subspace contained in each element of $\mf$.
	The Erd\H{o}s-Ko-Rado Theorem for vector space \cite{VEKR2,VEKR1,VEKR3} shows that a $t$-intersecting subfamily of ${V\b k}_q$ with maximum size is trivial when $\dim V>2k$.
	The structure of maximum sized non-trivial $t$-intersecting subfamily of ${V\b k}_q$ was determined via the parameter ``$t$-covering number" \cite{VHM1,VHM2}.

	Let $f$ be a non-degenerate sesquilinear form defined on $V$. An $m$-subspace $M$ of $V$ is called \emph{totally isotropic} if $f(x,y)=0$ holds for any $x,y\in M$. Denote the set of all totally isotropic subspaces of $V$ with respect to $f$ by $\mp$. Equipped with the inclusion relation, $\mp$ is a \emph{classical polar space}, denoted by the same symbol $\mp$. There are six kinds of classical polar spaces \cite{CPS}. In this paper, we consider the case that $\mp$ is a \emph{symplectic polar space}, where $n=2\nu$ and $f$ is an alternating bilinear form.

	For $0\le m\le\nu$, denote the set of all $m$-subspaces in $\mp$ by $\mp_m$. Let $\mf$ be a subfamily of $\mp_m$. If $\dim(F_1\cap F_2)\ge t$ for any $F_1, F_2\in\mf$, we say $\mf$ is $t$-\emph{intersecting}.

	If there exists a fixed $t$-subspace contained in each member of $\mf$, we say $\mf$ is \emph{trivial}.
	The maximum sized $t$-intersecting subfamily of $\mp_m$ was proven to be trivial if $\nu\ge m$ for $t=1$ \cite{SPEKR2,SPEKR1}, or $2\nu\ge3m+1$ with $(2\nu,q)\neq(3m+1,2)$ for $t\ge2$ \cite{SPEKRT}.

	For a subspace $A$ of $V$, let $\mm(m;A)$ denote the set of all members of $\mp_m$ contained in $A$. For $T\in\mp_t$ and $U\in\mp_m$ with $\dim(T\cap U)=m-1$, write
	$$\mh_1(T,U)=\{F\in\mp_m: T\subset F, \dim(F\cap U)\ge t\}\cup\mm(m;T+U).$$
	For $Z\in{V\b t+2}$, write
	$$\mh_2(Z)=\{F\in\mp_m: \dim(F\cap Z)\ge t+1\}.$$

	Our main result states the structure of non-trivial $t$-intersecting subfamilies of $\mp_m$ with maximum size under certain conditions.

	\begin{thm}\label{SSHMT}
		Let $\nu,m$ and $t$ be integers such that $\nu\ge m+1=t+2$, or $2\nu\ge\max\{3m+2,2m+2t+6\}$ with $m\ge t+2$.
		Suppose $\mf$ is a maximum sized non-trivial $t$-intersecting subfamily of $\mp_m$.
		Then one of the following holds.
		\begin{itemize}
			\item[\rm{(\ro1)}] If $1\le t\le\frac{m}{2}-1$, then $\mf=\mh_1(T,U)$ for some $T\in\mp_t$ and $U\in\mp_m$ with $T+U\in\mp_{m+1}$.
			\item[\rm{(\ro2)}] If $\frac{m}{2}-1<t\le m-1$, then $\mf=\mh_2(Z)$ for some $Z\in\mp_{t+2}$,
			or $\mf=\mh_1(T,U)$ for some $T\in\mp_1$ and $U\in\mp_3$ with $T+U\in\mp_{4}$ if $(m,t)=(3,1)$.
		\end{itemize}
	\end{thm}

	\section{Inequalities for sizes of two $t$-intersecting families}

	In this section, we always auusme that $m\ge t+2$. We prove bounds for the families $\mh_1(T,U)$ for $T\in\mp_t$ and $U\in\mp$ with $\dim(T\cap U)=t-1$. Then we compare the size of $\mh_1(T,U)$ with the size of $\mh_2(Z)$ under the condition that $T+U\in\mp_{m+1}$ and $Z\in\mp_{t+2}$.

	For positive integers $a$ and $b$ with $a\ge b$, define the \emph{Gaussian binomial coefficient} by
	$${a\b b}_q:=\prod_{0\le i<b}\dfrac{q^{a-i}-1}{q^{b-i}-1},$$
	and set ${a\b0}_q=1$.
	It is routine to check that the size of ${V\b k}_q$ is ${n\b k}_q$.
	From now on, we will omit the subscript $q$. We present some properties for the Gaussian binomial coefficient and the symplectic polar space, which will be used throughout our proof.

	For positive integers $m,i$ with $m>i$, we have
	\begin{equation}\label{FFS}
		q^{m-i}<\frac{q^m-1}{q^i-1}<q^{m-i+1},\quad q^{i(m-i)}\le{m\b i}<q^{i(m-i+1)}.
	\end{equation}
	For positives $m_1,m$ with $m_1\le m\le\nu$, let $N(\nu,m,m_1)$ be the number of members of $\mp_{m}$ containing a fixed member of $\mp_{m_1}$. From \cite[Theorem 3.38]{JS}, we derive that
	\begin{equation}\label{JS}
		N(\nu,m,m_1)=\prod\limits_{i=1}^{m-m_1}\dfrac{q^{2(\nu-m+i)}-1}{q^{i}-1}.
	\end{equation}
	For $1\le m\le \nu-1$, let $Q$ be an $m$-subspace of $V$ and $\alpha_1,\dots,\alpha_m$ a basis of $Q$. Note that the rank of matrix $(f(\alpha_i,\alpha_j))_{m\times m}$ does not depend on the choose of the basis. It is routine to check that its rank is $2s$ for some non-negative integer $s$. Then we say $Q$ has \emph{type $(m,s)$}. Note that $Q$ has type $(m,0)$ if and only if $Q$ is totally isotropic. By \cite[Theorem 3.27]{JS}, the number of members of $\mp_m$ contained in a fixed $(m+1,1)$-type subspace is $q+1$.

	For $\mf\subset\mp_m$ and a subspace $A$ of $V$, let $\mf_A$ denote the set of members of $\mf$ containing $A$.  Write
	$$f_0(\nu,m,t)={m-t+1\b1}N(\nu,m,t+1)-q{m-t+1\b2}N(\nu,m,t+2).$$

	\begin{lem}\label{PC}
		Let $\nu,m,t$ be positive integers with $m\ge t+2$ and $2\nu\ge3m-t+3$.
		Suppose $T\in\mp_t$ and $U\in\mp_m$ with $\dim(T\cap U)=t-1$.  Then the following hold.
		\begin{itemize}
			\item[\rm{(\ro1)}] If $T+U\in\mp_{m+1}$, then $|\mh_1(T,U)|>f_0(\nu,m,t)$.
			\item[\rm{(\ro2)}] If $T+U\not\in\mp_{m+1}$, then $|\mh_1(T,U)|<f_0(\nu,m,t)$.
		\end{itemize}
	\end{lem}
	\begin{proof}
		Observe that $\dim(T+U)=m+1$.

		(\ro1)	Write all $(t+1)$-subspaces of $T+U$ containing $T$ as $W_1,\dots,W_{{m-t+1\b1}}$.
		For each $i=1,\dots,{m-t+1\b1}$, let $\mw_i$ denote the set of all $m$-subspaces in $\mp$ containing $W_i$.
		For $F\in\mp_m$ with $T\subset F$, if $\dim(F\cap(T+U))=l\ge t+1$, it appears ${l-t\b1}$ times in $\mw_1,\dots,\mw_{{m-t+1\b 1}}$. For each $i=1,2,\dots,m-t$, write $\me_i=\{F\in\mp_m: T\subset F,\ \dim(F\cap(T+U))=t+i\}$.
		Then we have
		\begin{equation*}
			\begin{aligned}
				{m-t+1\b1}N(\nu,m,t+1)&=\sum_{i=1}^{m-t+1\b1}|\mw_i|=\sum_{i=1}^{m-t}{i\b 1}|\me_i|.
			\end{aligned}
		\end{equation*}
		Similarly, we have
		$$q{m-t+1\b2}N(\nu,m,t+2)=\sum_{i=2}^{m-t}q{i\b 2}|\me_i|.$$
		From ${i\b1}<q{i\b2}$ for $i\ge3$, we obtain
		\begin{equation}\label{PC2}
			\begin{aligned}
				f_0(\nu,m,t)\le\left|\left\{F\in\mp_m: T\subset F,\ t+1\le\dim(F\cap(T+U))\le t+2\right\}\right|.
			\end{aligned}
		\end{equation}
		Then (\ro1) follows from (\ref{PC2}) and the definition of $\mh_1(T,U)$.

		(\ro2) Since $T\in\mp_t$, $U\in\mp_m$ and $T+U\not\in\mp_{m+1}$, there exist $\alpha\in T\bs U$ and $\gamma\in U\bs T$ such that $f(\alpha,\gamma)\neq0$.
		Therefore,
		$$|\{S\in\mm(t+1;T+U): T\subset S\}|\le\left|\left\{S\in{T+U\b t+1}: T\subset S\right\}\bs\{\langle T,\gamma\rangle\}\right|={m-t+1\b1}-1.$$
		Notice that
		$$\{F\in\mp_m: T\subset F, \dim(F\cap(T+U))\ge t+1\}\subset\bigcup_{S\in\mm(t+1;T+U),\ T\subset S}\left(\mh_1(T,U)\right)_S$$
		and $T+U$ has type $(m+1,1)$. Then we have
		$$|\mh_1(T,U)|\le\left({m-t+1\b1}-1\right)N(\nu,m,t+1)+(q+1).$$
		By (\ref{FFS}), we get
		\begin{equation*}
			\begin{aligned}
				\dfrac{|\mh_1(T,U)|-f_0(\nu,m,t)}{N(\nu,m,t+2)}&\le-\dfrac{q^{2(\nu-t-1)}-1}{q^{m-t-1}-1}+q{m-t+1\b2}+\dfrac{q+1}{N(\nu,m,t+2)}\\
				&<-q^{2\nu-m-t-1}+q^{2(m-t)+1}+q+1\\
				&<0,
			\end{aligned}
		\end{equation*}
		as desired.
	\end{proof}

	Let $T\in\mp_t$, $U\in\mp_m$ with $T+U\in\mp_{m+1}$, and $Z\in\mp_{t+2}$. From \cite[Lemma 3.11]{JS}, we derive that $|\mh_1(T,U)|$ and $|\mh_2(Z)|$ do not depend on the choice of $T$, $U$ and $Z$. Write
	$$h_1(\nu,m,t)=|\mh_1(T,U)|,\quad h_2(\nu,m,t)=|\mh_2(Z)|.$$
	Notice that
	$$h_2(\nu,m,t)={t+2\b1}N(\nu,m,t+1)-q{t+1\b1}N(\nu,m,t+2).$$

	\begin{lem}\label{C2
		Let $\nu,m$ and $t$ be positive integers with $m\ge t+2$ and $2\nu\ge 3m+1$. Then the following hold.
		\begin{itemize}
			\item[\rm{(\ro1)}] If $1\le t\le\dfrac{m}{2}-1$, then $h_2(\nu,m,t)<h_1(\nu,m,t)$.
			\item[\rm({\ro2)}] If $\dfrac{m}{2}-1< t\le m-2$, then $h_2(\nu,m,t)\ge h_1(\nu,m,t)$ and equality holds if and only if $(m,t)=(3,1)$.
		\end{itemize}
	\end{lem}

	\begin{proof}

		(\ro1) By (\ref{FFS}), (\ref{JS}) and Lemma \ref{PC}, we have
		\begin{equation*}
			\begin{aligned}
				&\dfrac{h_1(\nu,m,t)-h_2(\nu,m,t)}{N(\nu,m,t+2)}\\
				>&\left({m-t+1\b1}-{t+2\b1}\right)\dfrac{q^{2(\nu-t-1)}-1}{q^{m-t-1}-1}-q\left({m-t+1\b2}-{t+1\b1}\right) \\
				>& q^{t+2}\cdot q^{m-2t-2}\cdot q^{2\nu-t-m-1}-q^{2(m-t)+1}\\

				>&0.
			\end{aligned}
		\end{equation*}
		Hence $h_1(\nu,m,t)>h_2(\nu,m,t)$.

		(\ro2) Let $U\in\mp_m$, $T\in\mp_t$ with $T+U\in\mp_{m+1}$. Suppose $\dfrac{m}{2}-1< t\le m-3$ and $\mf=\mh_1(T,U)$.
		Write
		$$\ma_j=\left\{F\in\mp_m: T\subset F,\ \dim(F\cap(T+U))=j\right\}$$
		for each $t+1\le j\le m$, and
		$$\ml=\left\{(I,F)\in\mp_{t+1}\times\mp_m: T\subset I\subset T+U,\ I\subset F\right\}.$$
		Let $\mb=\bigcup\limits_{j=t+1}^m\ma_j$.
		We have
		$$|\mf\bs\mb|=|\{F\in\mm(m;T+U): T\not\subset F\}|=q^{m-t+1}{t\b1}.$$
		By double counting $|\ml|$, we have
		\begin{equation*}\label{BT1}
			{m-t+1\b1}N(\nu,m,t+1)=|\ml|=\sum_{j=t+1}^m|\ma_j|\cdot{j-t\b1}=|\mb|+\sum_{j=t+2}^mq{j-t-1\b1}|\ma_j|.
		\end{equation*}
		Since $m\ge t+3$, by \cite[Theorem 2.10]{QY}, we have
		\begin{equation*}
			\begin{aligned}
				|\ma_{m-1}|={m-t+1\b2}\left((q^{\nu-m+1}+q^2){\nu-m-1\b1}+q^{2(\nu-m)}{2\b1}\right).
			\end{aligned}
		\end{equation*}
		Then by (\ref{FFS}), we get
		$$\dfrac{|\ma_{m-1}|}{|\mf\bs\mb|}>\dfrac{q^{2(m-t-1)}\cdot q^{2(\nu-m)-1}}{q^{m-t+1}\cdot q^t}\ge q^{3m-2t-2}\cdot q^{-(m+1)}>1.$$
		Observe that $\mf=(\mf\bs\mb)\cup\mb$. It is routine to check that
		\begin{equation*}
			\begin{aligned}
				h_1(\nu,m,t)=|\mf|=|\mb|+|\mf\bs\mb|<|\ml|-(q-1)|\ma_{t+2}|.
			\end{aligned}
		\end{equation*}
		If $2t>m-1$, since $2\nu\ge3m+1$, by (\ref{FFS}), (\ref{JS}) and Lemma \ref{PC}, we obtain
		\begin{equation*}
			\begin{aligned}
				\dfrac{h_2(\nu,m,t)-h_1(\nu,m,t)}{N(\nu,m,t+2)}&>\left({t+2\b1}-{m-t+1\b1}\right)\dfrac{q^{2(\nu-t-1)}-1}{q^{m-t-1}-1}-q{t+1\b1}\\
				&\ge q^{m-t+1}\cdot q^{2t-m}\cdot q^{2\nu-m-t-1}-q^{t+2}\\

				&>0.
			\end{aligned}
		\end{equation*}
		If $2t=m-1$, from (\ref{FFS}), (\ref{JS}), $2\nu\ge3m+2$, $t\ge2$ and \cite[Theorem 2.10]{QY}, we obtain
		\begin{equation*}
			\begin{aligned}
			&h_2(\nu,m,t)-h_1(\nu,m,t)\\>&(q-1)|\ma_{t+2}|-q{t+1\b1}N(\nu,m,t+2)\\
			>&(q-1){t+2\b 2}\cdot q^{\frac{t(t-1)}{2}+2(t-1)(\nu-2t-2)+2(t-1)}-q^t{t+1\b1}q^{2(t-1)(\nu-2t-2)+\frac{t(t-1)}{2}+2(t-1)}\\
			>&0.
			\end{aligned}
		\end{equation*}
		Hence $h_1(\nu,m,t)<h_2(\nu,m,t)$.

		Now suppose $m=t+2$. Let $M\in\mm(t+2;T+U)$ containing $T$.
		Write
		\begin{equation*}
			\begin{aligned}
				\mb_0&=\{F\in\mm(t+2;T+U): T\not\subset F\},\\
				\mb_1&=\{F\in\mp_{t+2}: T\not\subset F,\ \dim(F\cap M)=t+1\},\\
				\mb_2&=\{F\in\mp_{t+2}: T\subset F,\ \dim(F\cap M)\ge t+1\},\\
				\mb_3&=\{F\in\mp_{t+2}: F\cap M=T,\ \dim(F\cap(T+U))=t+1\}.
			\end{aligned}
		\end{equation*}
		Observe that
		$|\mb_0|={t+3\b1}-{3\b1}=q^3{t\b1}$.
		By \cite[Theorem 2.10]{QY}, we have
		$$
		|\mb_1|=\left((q^{\nu-t+1}+q^3){\nu-t-2\b1}+q^{2(\nu-t)-1}\right){t\b1},
		$$
		$$
		|\mb_3|=(q^{\nu-t+1}+q^4){\nu-t-3\b1}+q^{2(\nu-t-1)}{2\b1}.
		$$
		Notice that $\mh_1(T,U)=\mb_0\overset{.}{\cup}\mb_2\overset{.}{\cup}\mb_3$, $\mh_2(M)=\mb_1\overset{.}{\cup}\mb_2$, which imply that $h_2(\nu,m,t)-h_1(\nu,m,t)=|\mb_1|-|\mb_0|-|\mb_3|$.
		When $t=1$, it is routine to chech that $h_2(\nu,3,1)-h_1(\nu,3,1)=0$.
		If $t\ge2$, from $2\nu\ge3m+1=3t+7$ and (\ref{FFS}), we have
		\begin{equation*}
			\begin{aligned}
				&h_2(\nu,t+2,t)-h_1(\nu,t+2,t)\\
				\ge&\left((q^{\nu-t+1}+q^3)q^{\nu-t-3}+q^{2(\nu-t)-1}\right)q^{t-1}-\left((q^{\nu-t+1}+q^4)q^{\nu-t-3}+q^{2(\nu-t)}\right)-q^{t+3}\\

				\ge&q^{2\nu-2t-2}-q^{t+3}\\
				>&0,
			\end{aligned}
		\end{equation*}
		as desired.
	\end{proof}

	\section{Proof of Theorem \ref{SSHMT}}

	Suppose $\mf\subset\mp_m$ is a maximum sized non-trivial $t$-intersecting family. When $m=t+1$ with $\nu\ge m+1$, by \cite[Remarks (\ro2) in Section 9.3]{T1}, there exists $Z\in{V\b t+1}$ such that $\mf\subset{Z\b t+1}$. Observe that $Z\in\mp_{t+2}$ or $Z$ has type $(t+2,1)$. Then it is routine to check that $\mf=\mh_2(Z)$ where $Z\in\mp_{t+2}$. To finish the proof of Theorem \ref{SSHMT}, in view of Lemmas \ref{PC} and \ref{C2}, it is sufficient to show that $|\mf|<f_0(\nu,m,t)$ if $m\ge t+2$ with $2\nu\ge\max\{3m+2,2m+2t+6\}$ and $\mf$ is not one of the expected families in Theorem \ref{SSHMT}. We first prove the following lemma.

	\begin{lem}\label{FS}
		Let $A\in\mm(a,0;2\nu)$ with $a\le m-1$, and $\my\subset\mp_m$ be a $t$-intersecting family. Suppose that there exists $Y\in\my$ such that $\dim(A\cap Y)=r\le t-1$. Then there exists $R\in\mp_{a+t-r}$ with $A\subset R$ such that
		$$|\my_A|\le{m-r\b t-r}|\my_R|.$$
	\end{lem}
	\begin{proof}
		W.l.o.g., assume that $\my_A\neq\emptyset$. Write
		$$\mh=\{H\in\mm(a+t-r,0;A+Y): A\subset H\}.$$
		Since $\my$ is $t$-intersecting, for each $G\in\my_A$, we have $\dim(G\cap Y)\ge t$, which implies that $\dim(G\cap(A+Y))\ge a+t-r$. Consequently, $\mh\neq\emptyset$.
		Then there exists $R\in\mh$ such that $|\my_B|\le|\my_R|$ for any $B\in\mh$.
		Observe that
		$$|\mh|\le{m-r\b t-r},\quad\my_A\subset\bigcup_{B\in\mh}\my_B.$$
		Then the desired result holds.
	\end{proof}

	For $S\in\mp_s$, we say $S$ is a $t$-\emph{cover} of $\mf$ if $\dim(S\cap F)\ge t$ for any $F\in\mf$.
	Write
	$$\tau_t(\mf):=\min\{\dim S: S\ \text{is a}\ t\text{-cover of}\ \mf\}.$$
	Since $\mf$ is non-trivial and each member of $\mf$ is a $t$-cover of $\mf$, we have $t+1\le\tau_t(\mf)\le m$. Notice that each member of $\kn$ containing $S$ is also a member of $\mf$ because of the maximality of $\mf$. The proof in the following is divided into two cases.

	\subsection{The case $\tau_t(\mf)=t+1$}

	\noindent\textbf{Assumption 1.} Let $\nu,m,t$ be positive integers with $2\nu\ge\max\{3m+2,2m+2t+6\}$ and $m\ge t+2$. Suppose $\mf$ is a maximum sized non-trivial $t$-intersecting subfamily of $\mp_m$ with $\tau_t(\mf)=t+1$ and $\ms$ is the set of all $t$-cover of $\mf$ with dimension $t+1$. Set $X=\sum\limits_{S\in\ms}S$.

	\begin{lem}\label{XJ}
		Let $\nu,m,t,\mf,\ms$ and $X$ be as in Assumption 1. Then $\ms$ is a $t$-intersecting family with $t\le\tau_t(\ms)\le t+1$. Moreover, if $\tau_t(\ms)=t$, then $t+1\le\dim X\le m+1$ and $\dim(F\cap X)=\dim X-1$ for any $F\in\mf\bs\mf_T$, where $T\in\mp_t$ is contained in each member of $\ms$.
	\end{lem}
	\begin{proof}
		Suppose for contradiction that there exist $S_1,S_2\in\ms$ such that $\dim(S_1\cap S_2)<t$.
		By (\ref{JS}), the number of $A\in\mp_{t+2}$ containing $S_2$ is $N(\nu,t+2,t+1)={2(\nu-t-1)\b1}$. Observe that the number of members of $\mm(t+2;S_1+S_2)$ containing $S_2$ is at most ${t+1\b1}$. Since $2\nu>3(t+1)$, we have ${2(\nu-t-1)\b1}>{t+1\b1}$. Then
		\begin{equation*}
			\begin{aligned}
				\{A\in\mp_{t+2}: S_2\subset A\not\subset S_1+S_2\}=&\{A\in\mp_{t+2}: S_2\subset A\}\bs\{A\in\mp_{t+2}: A\subset S_1+S_2\}\neq\emptyset.
			\end{aligned}
		\end{equation*}
		For $S_2'\in\mp_{t+2}$ with $S_2\subset S_2'\not\subset S_1+S_2$, it is routine to check that $S_2'\cap S_1=S_2\cap S_1$. Repeating this discussion several times, we obtain $U_1, U_2\in\mp_m$ containing $S_1$ and $S_2$, respectively, satisfying that $U_1\cap U_2=S_1\cap S_2$. By the maximality of $\mf$, we have $U_1,U_2\in\mf$, a contradiction to the fact that $\mf$ is $t$-intersecting, as desired. From $\ms$ is $t$-intersecting, we get $t\le\tau_t(\ms)\le t+1$.

		Suppose $\tau_t(\ms)=t$ and $T$ is a member of $\mp_t$ contained in each member of $\ms$. 	For each $F\in\mf\bs\mf_T$ and $S\in\ms$, observe that $\dim(S\cap F)=t$. Then
		$$m+1\le\dim(T+F)\le\dim(S+F)=m+1.$$
		Note that $T+F\subset S+F$. We get $S+F=T+F$, which implies that $T+F=X+F$. From $\dim(T+F)=m+1$, we get $\dim(F\cap X)=\dim X-1$ and $\dim X\le m+1$. Together with $\dim X\ge \dim S=t+1$, the desired result holds.
	\end{proof}

	\begin{pr1}
		Let $\nu,m,t,\mf,\ms$ and $X$ be as in Assumption 1. If $\tau_t(\ms)=t$ and $\dim X=m+1$, then $\mf=\mh_1(T,U)$ for some $T\in\mp_t$ and $U\in\mp_m$ with $T+U\in\mp_{m+1}$.
	\end{pr1}

	\begin{proof}
		Suppose $T\in\mp_t$ is contained in each member of $\ms$.  For each $U\in\mf\bs\mf_T$, by Lemma \ref{XJ}, we have $\dim(U\cap X)=m$ and $X=T+U$, which imply that
		$$\mf\bs\mf_T\subset\left\{F\in\mm(m;T+U): T\not\subset F\right\}.$$
		For each $G\in\mf_T$, since $\dim(U\cap G)\ge t$ and $\dim(T\cap U)=t-1$, we have $\dim(G\cap(T+U))\ge t+1$, which implies that
		$$\mf_T\subset\left\{F\in\mp_m: T\subset F,\ \dim(F\cap(T+U))\ge t+1\right\}.$$
		Thus $\mf\subset\mh_1(T,U)$. Observe that $\mh_1(T,H)$ is non-trivially $t$-intersecting. By the maximality of $\mf$, we have $\mf=\mh_1(T,U)$.

		If $X\not\in\mp_{m+1}$, by Lemma \ref{PC}, we have $|\mf|<h_1(\nu,m,t)$. Note $\mf$ is maximum sized. Then we get  $T+U\in\mp_{m+1}$.
	\end{proof}

	\begin{pr1}\label{t+1,t,1}
		Let $\nu,m,t,\mf,\ms$ and $X$ be as in Assumption 1. If $\tau_t(\ms)=t$ and $t+1\le\dim X\le m$, then $|\mf|<f_0(\nu,m,t)$.
	\end{pr1}
	\begin{proof}
		We investigate $\mf$ in two cases.

		\noindent\textbf{Case 1. $\dim X=t+1$.}

		In this case, we have $|\ms|=1$. Let $S$ be the unique member of $\ms$. Observe that
		\begin{equation}\label{t+1,t,1,1}
			\mf=\mf_S\cup\left(\bigcup_{R\in\mm(t;S)}\mf_R\bs\mf_S\right).
		\end{equation}
		For each $R\in\mm(t;S)$, since $\tau_t(\mf)=t+1$, there exists $F\in\mf\bs\mf_S$ such that $\dim(F\cap R)\le t-1$.  Note that $\dim(F\cap S)=t$.
		We get $\dim(F\cap R)=t-1$.
		Let $H_1=F+R$.
		We have $\dim H_1=m+1$.
		It is routine to check that $S\subset H_1$ and $\dim(F'\cap H_1)\ge t+1$ for any $F'\in\mf_R\bs\mf_S$.
		Then
		\begin{equation}\label{t+1,t,1,2}
			\mf_R\bs\mf_S\subset\bigcup_{I\in\mm(t+1;H_1)\bs\{S\},\ R\subset I}\mf_I.
		\end{equation}

		For each $I\in\mm(t+1;H_1)\bs\{S\}$ with $R\subset I$,  there exists $F''\in\mf$ such that $\dim(F''\cap I)\le t-1$. Then from $R\subset I$ and $\dim(S\cap F'')\ge t$, we have $\dim(F''\cap I)=t-1$. By Lemma \ref{FS}, we have
		$$|\mf_I|\le{m-t+1\b1}N(\nu,m,t+2).$$
		Note that
		$$\left|\left\{I\in\mm(t+1;H_1): R\subset S, I\neq S\right\}\right|\le{m-t+1\b1}-1=q{m-t\b1},$$
		$$|\mf_S|\le N(\nu,m,t+1),\quad\left|\mm(t;S)\right|={t+1\b1}.$$
		It follows from (\ref{t+1,t,1,1}) and (\ref{t+1,t,1,2}) that
		$$|\mf|\le N(\nu,m,t+1)+q{m-t\b1}{t+1\b1}{m-t+1\b1}N(\nu,m,t+2).$$
		By (\ref{FFS}), (\ref{JS}) and $2\nu\ge 2m+2t+6$, we obtain
		\begin{equation*}
			\begin{aligned}
				&\dfrac{f_0(\nu,m,t)-|\mf|}{N(\nu,m,t+2)}\\
				\ge&q{m-t\b1}\cdot\dfrac{q^{2(\nu-t-1)}-1}{q^{m-t-1}-1}-q{m-t+1\b2}-q{m-t\b1}{t+1\b1}{m-t+1\b1}\\

				>&q^{2\nu-2t-1}-q^{2(m-t)+1}-q^{2m-t+3}\\

				>&0,
			\end{aligned}
		\end{equation*}
		as desired.

		\noindent{\bf Case 2. $t+2\le\dim X\le m$.}

		Write $k=\dim X$. Suppose that $T\in\mp_t$ is contained in each member of $\ms$. For $G\in\mf\bs\mf_T$, let $H_2=T+G$.
		For each $G'\in\mf_T$, if $G'\cap X=T$, since $\dim(G\cap G')\ge t$ and $T\not\subset G$, we have $\dim(G'\cap H_2)\ge t+1$.
		On the other hand, if $T\subsetneq G'\cap X$, then there exists a member of $\mm(t+1;X)$ containing $T$.
		Write
		$$\ma:=\left\{A\in\mm(t+1;H_2): T\subset A\not\subset X\right\},\quad\mb:=\left\{B\in\mm(t+1;X): T\subset B\right\}.$$
		We have
		\begin{equation}\label{t+1,t,2,1}
			\mf_T=\left(\bigcup_{A\in\ma}\mf_A\right)\cup\left(\bigcup_{B\in\mb}\mf_B\right).
		\end{equation}
		For each $A\in\ma$, by assumption, $A$ is not a $t$-cover of $\mf$.
		Then there exists $G''\in\mf\bs\mf_T$ such that $\dim(G''\cap A)\le t-1$.
		Together with $\dim(G''\cap T)=t-1$,
		we have $\dim(G''\cap A)=t-1$.
		By Lemma \ref{FS}, we have
		$$|\mf_A|\le{m-t+1\b1}N(\nu,m,t+2).$$
		By Lemma \ref{XJ}, we get $\dim(T\cap G)=t-1$, which implies that $\dim H_2=m+1$ and $X=T+(G\cap X)\subset H_2$. Then
		$$|\ma|\le{m-t+1\b 1}-{k-t\b1}=q^{k-t}{m-k+1\b1},\quad|\mb|\le{k-t\b1}.$$
		Together with $|\mf_B|\le N(\nu,m,t+1)$ for each $B\in\mb$, from (\ref{t+1,t,2,1}) we obtain
		\begin{equation}\label{S}
			|\mf_T|\le{k-t\b1}N(\nu,m,t+1)+q^{k-t}{m-k+1\b1}{m-t+1\b1}N(\nu,m,t+2).
		\end{equation}
		Write $\mc=\{C\in\mm(k-1;X): T\not\subset C\}$. Suppose $X\not\in\mp_k$, then $X$ has type $(k,1)$. Note that $|\mc|\le q+1$.
		Then
		$$|\mf\bs\mf_T|\le(q+1)N(\nu,m,k-1).$$
		If $X\in\mp_k$, we have $|\mc|={k\b1}-{k-t\b1}=q^{k-t}{t\b1}$. Then
		$$|\mf\bs\mf_T|\le q^{k-t}{t\b1}(N(\nu,m,k-1)-N(\nu,m,k)).$$
		Since $2\nu\ge 3m+2$ and $m\ge k\ge t+2\ge3$, by (\ref{FFS}) and (\ref{JS}), we have
		\begin{equation*}
			\begin{aligned}
				&\dfrac{q^{k-t}{t\b1}(N(\nu,m,k-1)-N(\nu,m,k))}{N(\nu,m,k-1)}\\
				=&q^{k-t}{t\b1}\left(1-\dfrac{q^{m-k+1}-1}{q^{2(\nu-k+1)}-1}\right)\\
				>&q^{k-t}\cdot q^{t-1}\cdot\left(1-q^{m+k-2\nu-1}\right)\\

				\ge&q+1,
			\end{aligned}
		\end{equation*}
		which implies that
		\begin{equation}\label{C}|\mf\bs\mf_T|\le q^{k-t}{t\b1}(N(\nu,m,k-1)-N(\nu,m,k))\end{equation}
		whenever $X$ is a member of $\mp_k$ or not.

		Suppose $k\ge t+3$.
		For each $t+3\le l\le m$, write
		$$u(l)=q^{l-t}{t\b1}(N(\nu,m,l-1)-N(\nu,m,l)).$$
		For each $t+3\le l<m$, since $2\nu\ge3m+2>m+l$, by (\ref{FFS}), (\ref{JS}) and $u(l+1)\le q^{l-t+1}{t\b1}N(\nu,m,l)$, we have
		\begin{equation*}
			\begin{aligned}
				\dfrac{u(l)}{u(l+1)}&\ge\dfrac{1}{q}\cdot\left(\dfrac{q^{2(\nu-l+1)}-1}{q^{m-l+1}-1}-1\right)>q^{2\nu-m-l}-\dfrac{1}{q}>1.
			\end{aligned}
		\end{equation*}
		Hence
		\begin{equation*}
			\begin{aligned}
				&\dfrac{f_0(\nu,m,t)-|\mf|}{N(\nu,m,t+2)}\\
				\ge&q^{k-t}{m-k+1\b1}\cdot\dfrac{q^{2(\nu-t-1)}-1}{q^{m-t-1}-1}-q{m-t+1\b2}-q^{k-t}{m-k+1\b1}{m-t+1\b1}-u(t+3)\\
				>&q^{2\nu-2t-1}-q^{2(m-t)+1}-q^{2m-2t+2}-q^{t+3}\\

				\ge&0
			\end{aligned}
		\end{equation*}
		follows from (\ref{FFS}), (\ref{JS}), (\ref{S}), (\ref{C}) and $2\nu\ge\max\{3m+2,2m+2t+6\}$, as desired.

		Now suppose that $k=t+2$. For each $C\in\mc$, since $T\not\subset C$, $C$ is not a $t$-cover of $\mf$. Thus there exists $E\in\mf$ such that $\dim(E\cap C)\le t-1$. Notice that $\dim(E\cap X)\ge t$, which implies that $\dim(E\cap C)=t-1$.
		Then $|\mf_C|\le{m-t+1\b1}N(\nu,m,t+2)$ follows from Lemma \ref{FS}. Together with $|\mc|\le\max\left\{q+1,{t+2\b1}-{2\b1}\right\}=q^2{t\b1}$, we get
		\begin{equation}\label{D}
			|\mf\bs\mf_T|\le q^2{t\b1}{m-t+1\b1}N(\nu,m,t+2).
		\end{equation}
		Then from (\ref{FFS}), (\ref{JS}), (\ref{S}), (\ref{D}) and $2\nu\ge\max\{3m+2,2m+2t+6\}$, we obtain
		\begin{equation*}
			\begin{aligned}
				&\dfrac{f_0(\nu,m,t)-|\mf|}{N(\nu,m,t+2)}\\
				\ge& q^{k-t}{m-k+1\b1}\cdot\dfrac{q^{2(\nu-t-1)-1}}{q^{m-t-1}-1}-q{m-t+1\b2}-\left(q^{k-t}{m-k+1\b1}+q^2{t\b1}\right){m-t+1\b1}\\
				>&q^{2\nu-2t-1}-q^{2(m-t)+1}-q^{2m-2t+2}-q^{m+3}\\

				\ge&0,
			\end{aligned}
		\end{equation*}
		as desired.
	\end{proof}

	\begin{pr1}
		Let $\nu,m,t$ and $\mf,\ms$ be as in Assumption 1. If $\tau_t(\ms)=t+1$, then $\mf=\mh_2(Z)$ for some $Z\in\mp_{t+2}$.
	\end{pr1}

	\begin{proof}
		By $\tau_t(\ms)=t+1$ and Lemma \ref{XJ}, $\ms$ is a non-trivial $t$-intersecting subfamily of ${V\b t+1}$. Therefore, by \cite[Remarks (\ro2) in Section 9.3]{T1}, there exists $Z\in{V\b t+2}$ such that $\ms\subset{Z\b t+1}$.

		Pick different members $A,B$ of $\ms$. Let $E=A\cap B$. Since $\ms$ is non-trivial, there exists $C\in\ms$ such that $E\not\subset C$. For any $F\in\mf$, if $E\subset F$, together with $\dim(F\cap C)\ge t$, we get $\dim(F\cap Z)\ge t+1$. If $E\not\subset F$, then $F\cap A$ and $F\cap B$ are two different subspaces of $Z$ with dimension at least $t$, which implies that $\dim(F\cap Z)\ge t+1$. Therefore, $\mf\subset\mh_2(Z)$. By the maximality of $\mf$ and $\mh_2(Z)$ is $t$-intersecting, we have $\mf=\mh_2(Z)$.

		Observe that $A+B=Z$ and $\dim(A\cap B)=t$. If $Z\not\in\mp_{t+2}$, $Z$ has type $(t+2,1)$. Then $|\mf|=(q+1)N(\nu,m,t+1)$.
		To finish the proof, it is sufficient to show that
		$$(q+1)N(\nu,m,t+1)<h_2(\nu,m,t).$$
		Since $2\nu\ge3m+2$, by (\ref{FFS}) and (\ref{JS}), we have
		\begin{equation*}
			\begin{aligned}
				\dfrac{h_2(\nu,m,t)}{N(\nu,m,t+1)}&={t+2\b1}-q{t+1\b1}\cdot\dfrac{q^{m-t-1}-1}{q^{2(\nu-t-1)}-1}>q^{t+1}-q^{-3}>q+1,

			\end{aligned}
		\end{equation*}
		as desired.
	\end{proof}

	\subsection{The case $\tau_t(\mf)\ge t+2$}

	\begin{pr1}\label{t+2}
		Let $\nu,m,t$ be positive integers with $2\nu\ge\max\{3m+2,2m+2t+6\}$ and $m\ge t+2$. Suppose $\mf$ is a maximum sized non-trivial $t$-intersecting subfamily of $\mp_m$ with $\tau_t(\mf)\ge t+2$. Then $|\mf|<f_0(\nu,m,t)$.
	\end{pr1}
	\begin{proof}
		Assume that $\tau_t(\mf)=r$.
		For each $S\in\ms$, observe that
		$$\mf\subset\bigcup_{W\in\mm(t;S)}\mf_W.$$
		For each $W\in\mm(t;S)$, w.l.o.g., assume that $\mf_W\neq\emptyset$. It is routine to check that ${m-a\b t-a}\le{m-t+1\b1}^{t-a}$ for each $0\le a\le t-1$. Then by Lemma \ref{FS}, there exist some totally isotropic subspaces $W_1,\dots,W_u$ such that $W=W_0\subset W_1\subset\cdots\subset W_u$ with $\dim W_u\ge r$, $\dim W_{u-1}<r$ and $$|\mf_{W_{i-1}}|\le{m-t+1\b1}^{\dim W_i-\dim W_{i-1}}|\mf_{W_i}|$$
		for each $i=1,\dots,u$. Notice that $\dim W_u\le m$ from $$|\mf_{W_u}|\ge|\mf_W|/{k-t+1\b1}^{\dim W_u-t}>0.$$
		When $t+2\le a\le m-1$ and $2\nu\ge3m-t$, we have
		\begin{equation*}
			\dfrac{N(\nu,m,a)}{N(\nu,m,a+1)}=\dfrac{q^{2(\nu-a)-1}}{q^{m-a}-1}\ge q^{2\nu-m-a}\ge q^{m-t+1}\ge{m-t+1\b1}.
		\end{equation*}
		Consequently,
		$$|\mf_W|\le{m-t+1\b1}^{\dim W_u-t}N(\nu,m,\dim W_u)\le{m-t+1\b1}^{r-t}N(\nu,m,r).$$
		Together with $|\mm(t;S)|\le{r\b t}$, we get
		$$|\mf|\le {r\b t}{m-t+1\b1}^{r-t}N(\nu,m,r).$$
		For each integer $s$ with $t\le s\le m$, write
		$$g(s)={s\b t}{m-t+1\b1}^{s-t}N(\nu,m,s).$$
		For $t\le s\le m-1$, from (\ref{FFS}), (\ref{JS}) and $2\nu\ge3m+2$, we obtain
		$$\dfrac{g(s+1)}{g(s)}=\dfrac{(q^{s+1}-1)(q^{m-t+1}-1)(q^{m-s}-1)}{(q^{s-t+1}-1)(q-1)(q^{2(\nu-s)}-1)}<q^{2m+s+2-2\nu}\le q^{3m+1-2\nu}<1.$$
		Then
		$$|\mf|\le{t+2\b2}{m-t+1\b1}^2N(\nu,m,t+2).$$
		Since $2\nu\ge 2m+2t+6$, by (\ref{FFS}) and (\ref{JS}), we get
		\begin{equation*}
			\begin{aligned}
				&\dfrac{f_0(\nu,m,t)-|\mf|}{N(\nu,m,t+2)}\\
				\ge&{m-t+1\b1}\cdot\dfrac{q^{2(\nu-t-1)-1}}{q^{m-t-1}-1}-q{m-t+1\b2}-{t+2\b2}{m-t+1\b1}^2\\
				>&q^{2\nu-2t-1}-q^{2(m-t)+1}-q^{2(m+2)}\\

				\ge&0,
			\end{aligned}
		\end{equation*}
		as desired.
	\end{proof}

	\vskip0.1in
	\noindent\textsc{Acknowledgement.} This research is supported by NSFC (12071039).


\section{Introduction}

	Conventional multi-atom lasers operate in the regime of weak coupling between the atoms and the light modes, possess a threshold and emit classical light. Single-atom lasers with strong coupling have been gaining increased attention since 1990s, as their behavior was predicted to  be drastically different from usual lasers \cite{oneatomlasers, Ginzel, Pellizzari, Pellizzari1, Horak, Loffer, Meyer, Meyer1, phstat, comment, Kilin, Nori}.  Among their non-standard features are presence or absence \cite{comment, phstat} of the pumping threshold depending on the parameters \cite{natureblatt}; self-quenching, when the cavity population decreases when the pumping rate is increased \cite{oneatomlasers}; photon blockade prohibiting coherent pumping of more than a single photon in the strong-coupling regime \cite{blockade, natureblatt, blockade97, blockadequbits1, blockadequbits2}; bistability of the Wigner function \cite{wignersset, multistability}; sub-Poissonian statistics of the emitted radiation leading to amplitude squeezing and purely quantum phenomena such as photon antibunching or super-Poissonian statistics with high intensity fluctuations and phase squeezing \cite{oneatomlasers}. Owing to these peculiar properties, a single-atom laser may be used for a wide range of applications as a source of non-classical radiation \cite{NielsenChuang,squeezedlight, qmeasurements, qimaging,qimaging1,qcommunications}.

	Experimentally, the first single-atom maser has been implemented in 1985 on Rydberg atoms \cite{firstmaser}. Then, several other platforms for single-atom lasing have been proposed, including trapped ions \cite{natureblatt, iontrap}, trapped Cs atom \cite{firstCs,Cs}, quantum dots \cite{bistabledots, thresholddyn, Nomura}, superconducting single-electron transistors \cite{wignersset, Astafiev} and flux qubits \cite{Neilinger, fluxqubit, qubitampl}. The single-atom masers based on superconducting cQED devices are particularly interesting since they can be used in tandem with superconducting artificial atoms \cite{guide}.

	In this work, we present a blueprint for a maser in this architecture using a transmon-type artificial atom \cite{koch}. The transmon may be regarded as a $\Xi$-system with the eigenstates $\ket{g}, \ket{e} ,\ket{f}$ referred to the ground, first and second excited states, respectively. The lasing transition $\ket{e} \rightarrow \ket{g}$ is in resonance with the reservoir lasing cavity coupled to the transmon. A classical lasing scheme in quantum optics is based on a Lambda-configuration of a three-level atom \cite{oneatomlasers}. However, transmon is a cascade system and cannot be excited directly to the second excited state, and thus the pumping is done via two-photon excitation of the $\ket{g} \leftrightarrow \ket{f}$ transition. The pumping processes for population inversion may be described as $\ket{g} \xleftrightarrow{2\hbar\omega} \ket{f} \rightarrow \ket{e}$. However, due to the small non-linearity of the transmon the $ \ket{f} \rightarrow \ket{e} $ process is only $\sqrt{2}$ times faster than $ \ket{e} \rightarrow \ket{g} $. Therefore, to increase the effective pumping rate, we suggest to couple the transmon to another low-Q cavity resonant with the $\ket{f} \rightarrow \ket{e}$ transition, very similar to what was done before for a chain of transmons \cite{ma2019dissipatively}. This cavity plays a role of artificial bath \cite{cavitybath, artificialbath, bathengineering, bathengineering1}, which relaxes the system to $\ket{e}$, ensuring the metastability of $\ket{f}$. Then, we solve the master equation numerically for the full system and study its transient and steady-state behavior for various combinations of parameters.

 	The suggested model is different from those studied before at least in two ways. First, we use a two-photon process for coherent pumping. Second, we employ engineered $\ket{f} \rightarrow \ket{e}$ dissipation via an auxiliary resonator, which at the same time significantly modifies the system spectrum. Actually, the latter fact alleviates the effect of the photon blockade \cite{blockade} caused by the large size of the vacuum Rabi splittings in the strong coupling regime \cite{Pellizzari}, similar to the case of the cascade laser proposed in 1995 \cite{cascade}. This effect allows to pump more than 40 photons in the strong-coupling regime. A similar breakdown of the photon blockade was investigated before in a simpler Jaynes-Cummings system \cite{BlockadeBreakdown}.

	Additionally, our simulations predict that depending on its parameters, the system can be operated in two distinct regimes. In one regime we observe a well-defined lasing threshold in our architecture which we also associate with the unusual energy level structure and blockade breakdown. In the other, no threshold is present, and there is a phase transition between the regimes with a sharp increase of the steady-state cavity population and the intensity fluctuations. A similar transition between different lasing regimes was previously predicted for the strong coupling \cite{wignersset, phases}.

	\section{Modeling the system}\label{sec:modeling}

	The scheme of the system is shown conceptually in \autoref{scheme}~(a). A tunable transmon is coupled to two cavities: the reservoir (right) and the auxiliary (left) with the coupling strengths $g_{r}$ and $g_{a}$, respectively. The reservoir has high internal and external Q-factors and should accumulate a considerable number of photons being in resonance with the $\ket{g} \leftrightarrow \ket{e}$ transition of the transmon, at frequency $\omega_{ge}$. The second low-Q cavity is resonant with the $\ket{e} \leftrightarrow \ket{f}$ transition at frequency $\omega_{ef}$, and thus provides an engineered metastability of the $\ket{f}$-state. The cavities do not interact directly with each other. To describe this system we apply the Jaynes-Cummings model \cite{jaynes-cummings} with an obvious modification to include two cavities, and in the rotating wave approximation (RWA) the resulting Hamiltonian reads:
	\begin{align}
	\label{hamilt}
	  \hat{H} &=  \hat{H}_{t} + \hat H_{d} + \sum_{\lambda=r,a}( \hat{H}^{(\lambda)}_{c} +  \hat{H}^{(\lambda)}_{i}),\\
	  \hat{{H}}_{t} &= \hbar \Delta
	  b^{\dagger} b +\frac{\hbar \alpha}{2}
	  b^{\dagger} b( b^{\dagger} b-1),\\
	  \hat{H}_{d} &= \frac{\hbar \Omega}{2} ( b + b^\dag),\\
	  \hat{H}^{(\lambda)}_{c}& = \hbar\Delta_c^{(\lambda)} (a_{\lambda}^\dag a_\lambda + 1/2),\\
	  \hat{H}^{(\lambda)}_{i} &= \hbar g_{\lambda}  ( b a_{\lambda}^\dag + b^\dag  a_\lambda),
	\end{align}

	\begin{figure}[]
		\centering\includegraphics[width=0.5\textwidth]{scheme}
		\caption{\textbf{(a)} Schematic configuration of the correctly aligned energy levels without coupling between subsystems. The reservoir is resonant with the $\ket{g} \leftrightarrow \ket{e}$ transition, and the auxiliary cavity is resonant with the $\ket{e} \leftrightarrow \ket{f}$ transition. $\ket{g},\ket{e},\ket{f} $ are in blue (medium gray), orange (light gray) and red (dark gray), respectively. \textbf{(b)} Schematic of the circuit. To the left is the auxiliary and to the right is the reservoir cavity. $C_g^{ge}$ ($C_g^{ef}$) is a capacitor that couples transmon with reservoir (auxiliary) cavity, and $C_{\kappa}^{ge}$ ($C_{\kappa}^{ef}$)~-- capacitor that provides the desired external quality factor}\label{scheme}
	\end{figure}

	\noindent where $\hat{H}_{t},\ \hat H_{d},\ \hat{H}^{(\lambda)}_{c}$ and $\hat{H}^{(\lambda)}_{i}$ are transmon, drive, cavity and interaction Hamiltonians, respectively, and $a_{r}$, $ a_{a}$, $b$ are the standard bosonic annihilation operators for the reservoir, auxiliary resonator and the transmon, respectively. The transmon parameters are: $\Delta = \omega_{ge} - \omega_d$, the detuning of drive frequency $\omega_d$ from the $\ket{g} \leftrightarrow \ket{e}$ transition; $\alpha < 0$, the anharmonicity; $\Omega$, the microwave drive amplitude. Next, $\Delta_c^{(r, a)} = \omega_{r,a} - \omega_d$ are the detunings of the reservoir and auxiliary cavities, respectively. In the ideal configuration, when the reservoir is resonant with the $\ket{g} \leftrightarrow \ket{e}$ transition, the auxiliary cavity is resonant with the $\ket{e} \leftrightarrow \ket{f}$ transition, and the drive is at the two-photon frequency $\omega_d = \omega_{gf}/2 = \omega_{ge} + \alpha/2$, we can find $\Delta_c^{r} = \Delta = - \alpha/2$ and $\Delta_c^{a} = \alpha/2$.

	Since relaxation processes are essential for lasing, modeling the system requires solution of the full Lindbladian master equation: for the transmon, the reservoir, and the auxiliary cavity, we use the relaxing collapse operators $\sqrt{\gamma} \hat b,\ \sqrt{\kappa_{r}} \hat a_{r},\ \sqrt{\kappa_{a}}\hat a_{a}$, respectively.

	To demonstrate the feasibility of the proposed device, we evaluate also its physical parameters. The optimal values for $g_{r, a},\ \kappa_{r, a},\ \omega_d$ were found with a Nelder-Mead algorithm aimed at maximizing the emitted power from the reservoir calculated by the input-output theory \cite{gardiner1984master}:
	\begin{equation}\label{power}
		P = \kappa_{r} N_{ss} \cdot \hbar \omega_{ge},
	\end{equation}
	where $N_{ss}$ is the steady-state number of photons in the reservoir.
	\begin{figure*}[]
		\centering\includegraphics[width=1\textwidth]{sim}
		\caption{ \textbf{(a)} Population of the transmon levels under a two-photon drive, $\Omega/2\pi = 25$ MHz. Upper panel: $\gamma = 0$, $g_{a} = 0$, middle panel: $\gamma = 1$ $\mu\text{s}^{-1}$, $g_{a} = 0$, lower panel: $\gamma = 1$ $\mu\text{s}^{-1}$, $g_{a}/2\pi = 8$ MHz, $\kappa_{a} = 138$ $\mu\text{s}^{-1}$.
		\textbf{(b)} Low-energy subspace of the system spectrum. Dashed arrows show two-photon pumping, wavy arrows~-- relaxation processes. Blue (medium gray) levels correspond to $\ket{g}$, orange (light gray)~-- $\ket{e}$, red (dark gray)~-- $\ket{f}$.
		\textbf{(c)} A single pumping stage showing transition from $N-1$ to $N$ photons in the reservoir and the participating eigenstates.}\label{sim}
	\end{figure*}
	In our simulations, truncating the reservoir Fock space at 60 photons, we find the optimal parameters to be the following: $g_{r}/2\pi = 6.5$ MHz, $g_{a}/2\pi = 23.5$ MHz, $\kappa_{r} = 0.31$ $\mu\text{s}^{-1}$, $\kappa_{a} = 138$ $\mu\text{s}^{-1}$ and $\omega_d = \omega_{gf}/2$, for $\Omega/2\pi = 25$ MHz. The algorithm converges to the same values starting from different initial parameters, so there is probably a single optimum in the parameter space. Since $\kappa_{a}$ is much larger than the two-photon Rabi frequency (see below for details), the number of photons in auxiliary cavity stays close to zero on average, so we truncate its Fock space to one photon. For the optimization, we use fixed realistic values for $\gamma = 0.1$ $\mu\text{s}^{-1}$, $\omega_{ge}/2\pi = 6$ GHz, $\alpha/2\pi = -200$ MHz, $\omega_r/2\pi = 6$ GHz, and $\omega_a/2\pi = 5.8$ GHz. For the optimal parameters we found $N_{ss} = 49$, which corresponds to $-132$ dBm of emitted power.

	Fabricating the superconducting cavities and the transmon with megahertz accuracy in frequency and anharmonicity should be experimentally attainable, and then the ideal configuration can be reached by tuning the transmon into the resonance with the reservoir via the external magnetic flux. Schematically, the circuit model of the device is shown in \autoref{scheme}~(b).

	The reservoir population is tolerant to small changes of the parameters, which can occur during fabrication. $N_{ss}$ becomes lower, but still enough to measure, due to mismatch between the auxiliary resonator frequency and $\ket{e} \leftrightarrow \ket{f}$ transition up to 30 MHz, while $\ket{g} \leftrightarrow \ket{e}$ frequency of the transmon can be always tuned directly to resonance with the reservoir cavity. The changes of $\kappa_{a}$ do not affect $N_{ss}$ significantly: higher $\kappa_{a}$ does not change $N_{ss}$ at all, while up to 40\% lower $\kappa_{a}$ still allows for pumping enough photons to measure. Variation of $g_{r}$ and $g_{a}$ affects the population significantly but slowly (see below for details), so their deviation from ideal values also should not be a problem. The only parameter whose deviation significantly affects the maser is $\kappa_{r}$: the lower $\kappa_{r}$ allows for pumping much more photons, but higher $\kappa_{r}$ may completely prohibit lasing. However, higher $\kappa_{r}$ can be coped with by increasing $g_{r}$.

	In addition, we have verified that adding dephasing to the transmon with $\gamma_\phi$ in the range of one $\mu\text{s}^{-1}$ and even tens of $\mu\text{s}^{-1}$ does not alter the number of pumped photons significantly.

	\section{Pumping dynamics}\label{sec:splittings}

	To calculate numerically the dynamics of the system, we use the QuTiP Python package \cite{qutip}. Below we denote the factorized states of the tripartite system as $\ket{\mu, M, N}$ where $\mu$ is the transmon state ($g$, $e$, $f$), and $M,\ N$ are the number of photons in the auxiliary resonator and the reservoir, respectively. To facilitate the interpretation of the simulation results, we split the task into several steps.

	First of all, we simulate the transmon under a two-photon drive by turning off the coupling to the cavities and the dissipation. The simulated evolution is shown in the upper panel of \autoref{sim}~(a) (here and below we truncate the transmon subspace to four states). One can clearly see two-photon Rabi oscillations when the $\ket{f}$ energy level is fully populated by the two-photon process. There are also small oscillations of the $\ket{e}$-level population from the virtual processes.

	In the middle panel \autoref{sim}~(a), we add natural relaxation to the transmon. The rate of the $\ket{f} \rightarrow \ket{e}$ transition is higher $n_{12}/n_{01}$ times than of the $\ket{e} \rightarrow \ket{g}$ one, where $n_{ij}$ is matrix element of charge operator in the basis of transmon eigenstates \cite{koch}. In this case, some population inversion may be achieved without any additional effort. Nevertheless, when the coupling to the auxiliary cavity it turned on, the population inversion becomes significantly larger, as shown in the lower panel of \autoref{sim}~(a).

	Finally, we turn on the coupling to the reservoir cavity. The resulting level structure of the full tripartite system is depicted in \autoref{sim}(b), and \autoref{sim}(c) shows the core process of pumping an additional photon in the reservoir. One can see that there are two types of Rabi splittings in the spectrum. The first, caused by the interaction between the reservoir and the $\ket{g} \leftrightarrow \ket{e}$ transition of the transmon, leads to the formation of the
	\begin{equation}
		\ket{geN}^\pm = \frac{1}{\sqrt{2}}(\ket{g,0,N+1} \pm \ket{e,0,N})
	\end{equation}
	dressed states composed of the factorized states with different reservoir populations. To avoid ambiguity, we label them using the lower number $N$. The size of this splitting is calculated as
	\begin{equation}\label{splitting1}
		\Delta_{ge}(N) = 2 g_{r} \sqrt{N+1}.
	\end{equation}

	The second splitting is due to the coupling between three degenerate states $\ket{f,0,N}$, $\ket{e,1,N}$, and $\ket{g,1,N+1}$. When the degeneracy is lifted, triplet eigenstates appear:
	\begin{equation}
	\begin{gathered}
		\ket{gefN}^0 = \sin\theta \ket{g,1,N+1} - \cos\theta \ket{e,1,N},
	\end{gathered}
	\end{equation}
	\begin{equation}
	\begin{gathered}
		\ket{gefN}^\pm = \frac{1}{\sqrt{2}} (\cos\theta \ket{g,1,N+1} \pm \ket{e,1,N} + \\ + \sin\theta \ket{f,0,N}),
	\end{gathered}
	\end{equation}
	\noindent where
	\begin{equation}
		\tan\theta = \frac{g_{a}}{g_{r}\sqrt{N+1}}
	\end{equation}
	is the mixing angle. The energy difference between $\ket{gefN}^+$ and $\ket{gefN}^-$ is
	\begin{equation}\label{splitting2}
		\Delta_{gef}(N) = 2 \sqrt{(N+1) g_{r}^2 + g_{a}^2}.
	\end{equation}
	\noindent As before, we label the triplet states by the lowest $N$ among their factorized components.

	We note that the first splitting (\autoref{splitting1}) is ordinary for any coupled atom-cavity system and grows as $\sqrt{N}$ with photon population of the cavity \cite{sqrtN}. In contrast, the second one (\autoref{splitting2}) is the characteristic feature of our system due to the presence of the auxiliary resonator and a certain frequency configuration. As we show below, both these splittings play the crucial role in the dynamics and allow the device to function in principle for arbitrarily large reservoir populations.

	\section{Simplified analytical solution}\label{sec:approx}

	\begin{figure}[]
		\centering \includegraphics[width=0.5\textwidth]{analytical}
		\caption{\textbf{(a)} Estimating the effective relaxation rate $\kappa_a^\text{eff}$ due to the auxiliary cavity. \textbf{(b)} Estimating the effective pumping rate due to the two-photon drive $\Omega_{2ph}$ and the effective relaxation $\kappa_a^\text{eff}$. \textbf{(c)} Estimating the effective relaxation rate due to the reservoir with $N$ photons inside. \textbf{(d)} $\ket{g}$ and $\ket{e}$ levels of transmon with effective $\ket{e} \rightarrow \ket{g}$ dissipation, which act as a two-level laser.}\label{analytical}
	\end{figure}

	To calculate approximately the reservoir population in the steady state, we consider only the $\ket{g}$ and $\ket{e}$ states with the effective pumping rate $\Gamma$ and effective dissipation rate $\kappa_r^\text{eff}$ through the reservoir (\autoref{analytical}(d)). In this simplified view, the balance condition $\kappa_r^\text{eff} = \Gamma$ allows to determine the reservoir population when the pumping ceases since $\kappa_r^\text{eff}$ is proportional to $N$ and $\Gamma$ is constant. We prove this in several steps below.

	Firstly, we calculate the effective $\ket{f} \rightarrow \ket{e}$ relaxation rate due to the auxiliary cavity (\autoref{analytical}(a)). Consider two coupled two-level systems: the first two-level system (on the left) represents $\ket{e}$ and $\ket{f}$ levels of transmon, and the second one (on the right)~-- $\ket{0}$ and $\ket{1}$ levels of the auxiliary resonator. The relevant master equation for a Hamiltonian $H = g_{a} (b a_{a}^\dag + a_{a} b^\dag)$ with a collapse operator $\sqrt{\kappa_a} a_{a}$ assuming initial state $\ket{f}$ has separate analytical solutions for the underdamped and overdamped regimes.

	In the underdamped, or strong coupling regime when $4 g_{a} > \kappa_a$, the population of $\ket{f}$ is
	\begin{equation}
	\rho_{ff} = \frac{e^{-\frac{\kappa_a t}{2}}}{\cos^2\theta} \cos^2(g't+\theta),
	\end{equation}
	\noindent where $g' = \sqrt{16 g_{a}^2 - \kappa_a^2}/4$, $\tan \theta = \frac{\kappa_a}{g'}$. Using the decay rate of the vaccuum Rabi oscillations, one may assume that
	\begin{equation}\label{kappa_strong}
	\kappa_a^\text{eff} = \frac{\kappa_a}{2}.
	\end{equation}

	In the overdamped, or weak coupling regime where $4 g_{a} < \kappa_a$,
	\begin{equation}
	\rho_{ff} = \frac{e^{-\frac{\kappa_a t}{2}}}{\cosh^2\theta} \cosh^2(\gamma' t+\theta),
	\end{equation}
	\noindent where $\gamma' = \sqrt{\kappa_a^2 - 16 g_{a}^2}/4$, $\tan \theta = \frac{\kappa_a}{\gamma'}$.
	The decay rate is determined by the slowest exponent, so
	\begin{equation}\label{kappa_weak}
	\kappa_a^\text{eff} = \frac{\kappa_a - \sqrt{\kappa_a^2 - 16 g_{a}^2}}{2}.
	\end{equation}

	One can see that while $g_{a}$ is fixed, increasing $\kappa_a$ at first increases but then decreases $\kappa_a^\text{eff}$ upon transition from the underdamped to the overdamped regime ($\kappa_a^\text{eff} \rightarrow 4 g_a^2/\kappa_a$ there). Physically, this effect can be explained by the broadening of the $\ket{1}_{a}$ level, which prevents transfer of energy from $\ket{f}$.
	Thus, $\kappa_a^\text{eff}$ reaches its maximum when $\kappa_a = g_{a}$ which is consistent with our numerical simulations.

	Next, we calculate effective $\ket{g} \rightarrow \ket{e}$ pumping rate $\Gamma$. The simplified model of the system with two-photon pumping and $\ket{f} \rightarrow \ket{e}$ dissipation is shown on \autoref{analytical}(b). The two-photon process is replaced by the equivalent one-photon drive of strength $\Omega_{2ph}$ (see Appendix \ref{rate} for details). This step is only approximate as it does not take into account the full complexity of the energy spectrum and the possibility of the off-resonant two-photon drive.
	From \autoref{analytical}(b) one can see that the tree-level system with pumping and dissipation is completely equivalent to the two coupled two-level systems in \autoref{analytical}(a) with $g_{a}$ replaced by $\Omega_{2ph}/2$. The solution for the such a system is the same: in the strong coupling regime $2 \Omega_{2ph} > \kappa_a^\text{eff}$ the effective pumping rate is
	\begin{equation}
		\Gamma = \frac{\kappa_a^\text{eff}}{2},
	\end{equation}
	\noindent and for the weak coupling regime $2\Omega_{2ph} < \kappa_a^\text{eff}$
	\begin{equation}\label{gamma_eff}
	\Gamma = \frac{\kappa_a^\text{eff} - \sqrt{(\kappa_a^\text{eff})^2 - 4 \Omega_{2ph}^2}}{2}.
	\end{equation}

	In our case, the strong coupling is not achievable due to the weakness of the two-photon drive. Consequently, we take \autoref{gamma_eff} to calculate effective pumping rate. For the optimal parameters from \autoref{sec:modeling}, it leads to $\Gamma = 5.3$ MHz.

	The next step is the calculation of effective dissipation rate $\kappa_r^\text{eff}$ due to the reservoir populated with $N$ photons (\autoref{analytical}(c)). Again, the system is equivalent to \autoref{analytical}(a). Coupling strength $g_{r}$ should not be much lower then $g_{a}$, otherwise the Rabi splittings will disturb pumping (see next section for detail). Assuming this, for optimal $g_{a}$ and $\kappa_r<1$ $\mu \text{s}^{-1}$, the condition for the strong coupling regime $4\sqrt{N}g_{r} > N\kappa_r$ is met even for high reservoir populations ($N \sim 100$). Therefore,
	\begin{equation}
		\kappa_r^\text{eff} = \frac{N_{ss} \kappa_r}{2}.
	\end{equation}

	The energy balance condition is $\kappa_r^\text{eff} = \Gamma$ (\autoref{analytical}(d)). Consequently, the maximum number of pumped photons can be expressed as
	\begin{equation}
		N_{ss} = 2\frac{\Gamma}{\kappa_r}.
	\end{equation}

	For the optimal parameters, this equation gives $N_{ss} = 34$. The exact simulated value is $N_{ss} = 49$. The 30\% difference is probably caused by the differences in the real level structure due to the coupling to both cavities.

	\section{Overcoming photon blockade}\label{sec:blockade}

	The simplified model studied in the previous section does not take into account the effects of the photon blockade as it replaces the coherent pumping by an effective incoherent one. However, the auxiliary cavity actually allows to treat the system within the simplified model. To see this, one should first consider a simple device that does not include an auxiliary resonator, operates in the strong-coupling regime and is subject to photon blockade as follows. Due to the $\sqrt{N}$ growth of the $\Delta_{ge}(N)$, there is no single frequency that would satisfy the resonance condition of the two-photon pumping for an arbitrary population of the reservoir. The pumping at the $\ket{g,0,N} \leftrightarrow \ket{f,0,N}$ frequency becomes off-resonant when $\langle \hat N \rangle = \langle a_{r}^\dag a_{r} \rangle \geq 1$ \cite{Pellizzari} (see \autoref{sim}~(b)) as there would be no triplet Rabi splittings and the detuning would be $\delta = \Delta_{ge}(N)/2$. This effect is observed for all coherently pumped single-atom lasers: $\delta = g_{r} \sqrt{N} \rightarrow \infty$ for $N \rightarrow \infty$ \cite{blockade}.

	However, in the presence of the auxiliary resonator one can see from \autoref{sim}~(b,c) that the detuning of the two-photon drive is instead calculated as

	\begin{equation}\label{detuning}
	\delta = \left[\Delta_{gef}(N+1)-\Delta_{ge}(N)\right]/2
	\end{equation}

	One can see that $\delta \rightarrow 0$ for $N \rightarrow \infty$ (\autoref{sim}(c)), which means that it would be possible to use monochromatic pumping at $\omega_{gf}/2$ for sufficiently large values of $\langle \hat N \rangle$, even though it would be slightly off-resonant for low $\langle \hat N \rangle$. In other words, if it is possible to pump several photons in the reservoir off-resonantly, then the subsequent pumping will become resonant and similar to the simplified model. We demonstrate numerically that the latter is possible due to the large dissipative bandwidths of the states in the $\ket{gefN}^{\pm,0}$ splittings. The low Q-factor of the auxiliary resonator allows off-resonant transitions to them even if the detuning is of the order of several MHz. In \autoref{pumping}~(a, b), we compare the steady-state solution with the transient solutions at 1 an 2 $ \mu s $. The steady-state picture shows that pumping at a frequency close to $\ket{g,0,0} \leftrightarrow \ket{f,0,0}$ is indeed optimal in the long term, despite that at first the sideband $\ket{geN}^+ \leftrightarrow \ket{gefN}^+$ has a higher pumping rate, see Appendix \ref{rate} for details; the sidebands are subject to the photon blockade and eventually die off for large $\langle \hat N \rangle$.

	As a result, with optimal parameters, we can pump $\approx 50$ photons while maximizing output power (-132 dBm at the optimum). Maximizing the reservoir population allows pumping at least 90 photons in the steady state (above this point, the simulation becomes intractable due to the size of the necessary Hilbert space).

	\section{Different lasing regimes and phase transitions}\label{sec:regimes}

	While the simplified model is applicable for large reservoir populations, numerical solution is required to accurately characterize the behavior of the system for arbitrary parameter combinations. We find that the most important parameters are the coupling strengths between the transmon and the cavities, so in this section we study in detail how they affect the reservoir steady state properties.

	\begin{figure}[]
		\centering \includegraphics[width=0.5\textwidth]{pumping}
		\caption{\textbf{(a, b)} Average number of photons in the reservoir depending on the drive frequency $\omega$ for $g_{r}/2\pi = g_{a}/2\pi = 15$ MHz, $\kappa_{r} = 0.2$ $\mu\text{s}^{-1}$, $\kappa_{a} = 128$ $\mu\text{s}^{-1}$, $\Omega/2\pi = 20$ MHz (point A in (e, f)) in the steady state (a) and in 1 $\mu$s and 2 $\mu$s after the pumping was started (b). \textbf{(c, e)} Pumped photons and \textbf{(d, f)} logarithm of Fano factor vs. coupling strengths for $\kappa_{a} = 138$ $\mu\text{s}^{-1}$, $\kappa_{r} = 0.2$ $\mu\text{s}^{-1}$. \textbf{(c, d)} $\Omega/2\pi = 15$ MHz, \textbf{(e, f)} $\Omega/2\pi = 20$ MHz. The letters are different sets of the parameters used in the other subsections (O depicts the optimal parameters). Orange dashed line marks the transition to the threshold-bearing regime and blue dashed line the thresholdless regime.}\label{pumping}
	\end{figure}

	\autoref{pumping}(c-f) shows the reservoir population and the logarithm of the Fano factor $F = \frac{\langle N^2 \rangle - \langle N \rangle^2}{\langle N \rangle}$ in the steady state depending on $ g_{r, a} $ for two different values of $\Omega$. $F$ is a measure of intensity fluctuations of the laser field: $F<1$ indicates sub-Poissonian statistics and photon antibunching, $F=1$~-- coherent Poissonian field, and $F>1$~-- super-Poissonian statistics with high intensity fluctuations.

	For both driving strengths, one can see an area with a significant increase in the $N_{ss}$ and $F$ when the coupling strengths $ g_{r,a} $ become large enough. We find that this area marks the phase transition between two regimes of lasing: threshold-bearing and threshold-less (see \autoref{threshold-sec} for detail). The lasing threshold is said to exist if there is some pumping power value for which a rapid increase of intensity fluctuations and usually of the steady-state photon number in the cavity is observed \cite{comment}.

	\begin{figure*}[]
		\centering\includegraphics[width=1\textwidth]{threshold}
		\caption{ Reservoir population \textbf{(a)} and the Fano-factor \textbf{(b)} for the threshold-less and threshold-bearing regimes (points B-F in \autoref{pumping}, values for $g_{a,r}/2\pi$ in the legend are in MHz). For all panels $\kappa_{a} = 138$ $\mu\text{s}^{-1}$, $\kappa_{r} = 0.2$ $\mu\text{s}^{-1}$. \textbf{(c-d)} Reservoir Fock state probability for various $\Omega$: \textbf{(c)} point G, \textbf{(d)} point F. \textbf{(e-h)} Steady-state reservoir Wigner functions: \textbf{(e)} for (c), $\Omega/2\pi = 19.5$ MHz;  \textbf{(f)} for (c), $\Omega/2\pi = 21$ MHz;
	    \textbf{(g)} for (d), $\Omega/2\pi = 23$ MHz;
	    \textbf{(h)} for (d), $\Omega/2\pi = 26$ MHz.}\label{threshold}
	\end{figure*}

	For the smaller drive amplitude, the phase transition is not very pronounced; however, for the larger drive, it is sharp (orange dashed line). After the transition, the lasing is characterized by the large $ N_{ss} $ and $F<1$. Finally, the area of the threshold-bearing regime grows when the pumping rate is increased.

	Qualitatively, the location of the threshold-bearing regime in the parameter space can be explained as follows. Firstly, the large coupling constant $g_{r}$ allows pumping the reservoir quicker, which explains why the transition requires $g_{r}/2\pi$  to exceed 12 and 8 MHz, for the weak and strong drive, respectively. Secondly, when $g_{a}/2\pi < 8$ MHz, the coupling to the auxiliary resonator is too low to reach significant $\kappa_{a}^\text{eff}$ and to produce sufficient population inversion. Finally, increasing $g_{a}$ while  keeping $g_{r}$ fixed reduces $N_{ss}$ due to the increasing detuning discussed in \autoref{sec:blockade}.

	In the thresholdless regime in the bottom left corner of \autoref{pumping}(c, e) (dashed blue line)  it is still possible to pump a considerable number of photons (about 5-10) while $F \gtrsim 1$ in this area. This regime is standard for single-atom lasers in the intermediate and low-coupling regime. The coupling constants in the thresholdless regime are low, so the energy splittings $\Delta_{ge}(N)$ and $\Delta_{gef}(N)$ do not play any role due to the dissipative line broadening.

	\section{Threshold and bistability}\label{threshold-sec}

	Single-atom lasers previously discussed in literature were thresholdless in the strong coupling regime \cite{natureblatt, firstCs,comment}, but could have one or two thresholds in the regime of weak and intermediate coupling \cite{natureblatt, iontrap, coher}. Therefore, below we discuss why in our case it is possible to observe a threshold even in the strong coupling regime.

	\autoref{threshold}(a, b) shows the steady-state values for $N_{ss}$ and $F$ for five different combinations of coupling constants, corresponding to points B-F in \autoref{pumping} (e, f), depending on the drive amplitude $ \Omega $. For the first four sets of parameters (B-E), there is no well-defined threshold, even for $g_{r}/2\pi = 5$ MHz, $g_{a}/2\pi = 8$ MHz (E), which has a small maximum of $F$ at $\Omega/2\pi = 12$ MHz. The photon distribution in the regimes B-E can be super-Poissonian as well as sub-Poissonian depending on the drive power. The observed behavior is consistent with previous works \cite{phstat, coher}.

	For the combination B, $g_{r}/2\pi = 1$ MHz and $g_{a}/2\pi = 2$ MHz, we observe self-quenching probably caused by the dressing of the energy levels by the strong drive. This effect is still not clear to us, because self-quenching remains even if we manually tune pumping frequency into resonance with the correct transition to one of the dressed states, even though the maximum of pumped photons shifts to the larger $\Omega$.

	The most interesting feature in \autoref{threshold} (a),(b) is the presence of a well-defined threshold for the combination G: $g_{r}/2\pi = 8$ MHz and $g_{a}/2\pi = 15$ MHz and F: $g_{r}/2\pi = 8$ MHz and $g_{a}/2\pi = 11$ MHz. Our simulations show that such a threshold is present for any point above the orange dashed line at \autoref{pumping} (e, f). Above the threshold, for some parameters there is also no manifestations of the photon blockade~-- the average number of photons increases steadily with growing $\Omega$. As mentioned above, this is different from the previous results for single-cavity systems with coherent pumping in the strong-coupling regime \cite{iontrap,Cs}.

	\autoref{threshold}~(c, e, f) illustrate the cavity state in the vicinity of the threshold for the point G, and \autoref{threshold}~(d,g,h)~-- for the point F. One can see bistability which manifests itself as a double-ring structure of the Wigner function. At a sufficiently high driving power for the bistability to vanish, only the outer ring of Wigner function remains and the Fano factor is significantly lower than one.

	We explain the bistability in the studied system as follows. We consider the intermediate cavity states with $N$ approximately between 2 and 25 for \autoref{threshold} (c) and between 5 and 25 for \autoref{threshold} (d). These states can not be populated via the coherent pump due to the large detuning they have due to vacuum Rabi splittings (see Eq. \ref{detuning}), and the relaxation dominates their population dynamics. However, the upper states can be pumped resonantly, as discussed before, so the system can remain there indefinitely. Therefore, with increasing drive strength the higher levels gradually absorb more and more of the population from the lower ones as the population is leaking better and better through the ``blocked'' states. In the middle of the threshold transition, this results in the apparent bimodal distribution of photons and very high intensity fluctuations. Additionally, in Appendix \ref{wignerappendix}, we show that bistability would not be observed if the Rabi-splittings did not depend on the photon number.

	This mechanism resembles the one studied in \cite{bistablecoupling}, where the high-order multiphoton processes are becoming closer and closer to each other in frequency, while the single-photon processes are prone to photon blockade. Additionally, the similar mechanism was experimentally confirmed in \cite{BlockadeBreakdown}, where the photon blockade in the strong-coupled qubit-resonator system was broken by the weakly-coupled third level of the qubit, making resonant pumping of the cavity possible.

	\section{Conclusion}

   	In this work, we studied theoretically a single-atom maser based on two-photon coherent pumping of a transmon. We find that to attain high population inversion, one needs to artificially increase the relaxation rate of its $\ket{e} \rightarrow \ket{f}$ transition and show that this is feasible by coupling it to an auxiliary low-Q resonator.

    Our master equation simulations predict two distinct regimes of operation for the device with the presence or absence of the lasing threshold. The boundary between the regimes is determined by the coupling strengths between the transmon and the cavities.

    In the threshold-less regime, our system demonstrates non-classical behavior similar to the previously known for single-atom lasers, including both sub-Poissonian and super-Poissonian statistic, and self-quenching.

	\begin{figure*}[]
		\centering \includegraphics[width=1\textwidth]{wigner3}
		\caption{Visualized reservoir states for the system with the following parameters: $\kappa_{a} = 138$ $\mu\text{s}^{-1}$, $\kappa_{r} = 0.2$ $\mu\text{s}^{-1}$, $g_{r}/2\pi = 8$ MHz, $g_{a} = 11$ MHz (point F at \autoref{pumping}(e, f), for different $\Omega$.  \textbf{(a)} With ordinary interaction Hamiltonian, \textbf{(b)} with modified interaction Hamiltonian. \textbf{(c,d)} Number of photons in the reservoir (c) and Fano factor (d) for ordinary and modified interaction Hamiltonian. \textbf{(e)} Number of photons in the reservoir depending on the detuning of the drive frequency $\omega$ from a frequency of reservoir resonator $\omega_{ge}$. The system is modified such a way that allows only one-photon direct reservoir pumping: the number of the transmon states is truncated to 2, when $g_r/2\pi = 8$ MHz, $g_a = 0$ and the other parameters are the same as in \autoref{threshold}. }\label{wignerpeaks}
	\end{figure*}

	In the other regime, for which the coupling strengths are higher, we observe a well-defined lasing threshold, marked by an increase of intensity fluctuations. We associate its presence with the interplay between the photon blockade and the peculiarities of the energy spectrum caused by the auxiliary resonator, and the consequent bistability of Wigner function. The light statistics in the reservoir in this regime is highly sub-Poissonian, with $F<0.5$, and its population is of order of tens of photons.

	Finally, via numerical optimization, we find an optimal set of parameters for a device that allows pumping about 50 photons in the strong-coupling regime, which corresponds to $-132$ dBm of emitted power. The statistics in this regime is predicted to be sub-Poissonian with $F = 0.27$. We believe that due to the non-classical properties such maser may have a wide range of applications for microwave quantum optics on a chip \cite{qimaging,NielsenChuang,squeezedlight}.

	\section{Acknowledgments}

	The investigation was conducted with the support of Russian Science Foundation, Grant No. 16-12-00095. We also acknowledge support from the Ministry of Education and Science of the Russian Federation in the framework of the Increased Competitiveness Program of the National University of Science and Technology MISIS (Contract No. K2-2020-022).


\chapter{Introduction}

        \section{Introduction}
        The n-harmonic 2 dimesion has been studied by Doscharis in \cite{Doscharis}, and  neumann Problem 1-harmonic m dimesion have been solved by Otuf in \cite{Otuf}. We are intersted in n-harmonic equation m dimension which is in the form of
           \begin{equation}
           \triangle_{m}^{n} u(x)=0,
           \end{equation}
           where $  x=(x_1,x_2,........, x_m) . $
           This equations has high complexity since it is a partial differential equation of order 2n. However, separation of variables ia an intersting approch to solve this kind of problem. At begining,  we need to put our attention in driving new concept which decrease the complexity of solving this problem. Multinomial theorem of n-harmonic equation which is coming directly from the original multinomial theorem allows to wrtie the n-harmonic equation as summuation of ordinary differential equations. This lead to separate these equations and write them as sequnce of equations. Solving these equations will desgin an important gereralization. Moreover, we are intersted in solving unbounded boundary value problem that satisfying $ \triangle_m^{n} u(x)=0 $ and then application of its convolutions theorem.

           \section{ Multinomial Theorem of n-Harmonic Equation}

          We formalize multinomial theorem to desgin short approch to separate the variables of $ \triangle_m^{n} u(x) ,$  and
          write them as summutions of ordinary differential equations. Supose for   $
          y_i=  \frac{X_{i}^{(2)}}{X_{i}} ,$  we have $ y_i^{(n)} = \frac{X_{i}^{(2n)}}{  X_i } $ . We let
 $  u(x)= \prod_{i = 1}^{m-1} X_i (x_i)   X_m (x_m) ,$ then apply $ \triangle_m^{n}  $ to  $ u(x) $ and then divide by $ \prod_{i = 1}^{m-1} X_i(x_i)   X_m (x_m)  $  leads to the  multinomial theorem of n-harmonic equation which is as following

          \begin{equation}
           \left(   \sum_{r=1}^{m} y_i \right)^{(n)}  = \sum_{h_1 + h_2+......+h_m=n}  \binom {n} {h_1, h_2, ......h_m}  \prod_{i = 1}^{m} y_i^{h_i}.
          \end{equation}
          Where
          \begin{equation}
          \binom {n} {h_1, h_2, ...... , h_m}  =   \frac{n!}{h_1 ! h_2 ! .........h_m !},
          \end{equation}
           $ (n) $ is order of dervative and $ h_i \neq 0. $
          Furthermore, mutiplication of any two or more ordinary differential equations, yields  sum of their orders. For example
          $ X_{2}^{(2)} X_{2}^{(4)}= X_{2}^{(6)} $
      We list some examples
           For $ \triangle_2^{2} u(x)=0 , $ we have
          \begin{equation} ( y_1 + y_2 )^{(2)} =   y_1^{(2)}+ 2 y_1y_2  +  y_2^{(2)} = \frac{X_{1}^{(4)}}{X_{1}}
          + 2 \frac{X_{1}^{(2)}}{X_{1}} \frac{X_{2}^{(2)}}{X_{2}} + \frac{X_{2}^{(4)}}{X_{2}}   .
          \end{equation}
          For $ \triangle_3^{2} u(x)=0 , $ we have
     \begin{equation}
     ( y_1 + y_2  )^{(3)}= y_1^{(3)}+  3 y_1^{(2)} y_2  +  3 y_1 y_2^{(2)}  + y_2^{(3)} =\frac{X_{1}^{(6)}}{X_{1}} + 3 \frac{X_{1}^{(4)}}{X_{1}}  \frac{X_{2}^{(2)}}{X_{2}} + 3 \frac{X_{1}^{(2)}}{X_{1}} \frac{X_{2}^{(4)}}{X_{2}} + \frac{X_{1}^{(6)}}{X_{1}}  .
       \end{equation}

          For $ \triangle_3^{3} u(x)=0 , $ we have
       \begin{equation}
    ( y_1 + y_2  + y_3)^{(3)}=   \frac{X_{1}^{(6)}}{X_{1}} + \frac{X_{2}^{(6)}}{X_{2}} +  \frac{X_{3}^{(6)}}{X_{3}}
          +  3 \frac{X_{1}^{(4)}}{X_{1}}  \frac{X_{2}^{(2)}}{X_{2}} + 3 \frac{X_{1}^{(2)}}{X_{1}} \frac{X_{2}^{(4)}}{X_{2}}   \end{equation}
           \begin{equation}   +  3 \frac{X_{1}^{(4)}}{X_{1}}  \frac{X_{3}^{(2)}}{X_{3}} + 3 \frac{X_{1}^{(2)}}{X_{1}} \frac{X_{3}^{(4)}}{X_{3}} +  3 \frac{X_{3}^{(4)}}{X_{3}}  \frac{X_{2}^{(2)}}{X_{2}} + 3 \frac{X_{3}^{(2)}}{X_{3}} \frac{X_{2}^{(4)}}{X_{2}}  .
          \end{equation}

\section{Generalization Of Separation Of Variables}

           Let    $ x=(x_1,x_2,........x_m) $ be vector in $ \mathbb{R}^{m}. $ The n-harmonic equation is
             \begin{equation}
           \triangle_m^{n} u(x)=0.
           \end{equation}
                Let $ \lambda_{i_j} $ to be an constants in $ \mathbb{R},  $
             where $ i=1,2,3,.......m ,$ and $j=1,2,......,n.$
             We seek solution of the form
              \begin{equation}
               u(x)=\prod_{i = 1}^{m-1} X_i(x_i)   X_m (x_m) .    \label{eq-110}
              \end{equation}
              Apply method of sepration of varibles to get multinomial theorem of n-harmonic equation that is  \begin{equation}
           \left(   \sum_{r=1}^{m} y_i \right)^{(n)}  = \sum_{h_1 + h_2+......+h_m=n}  \binom {n} {h_1, h_2, ......h_m}  \prod_{i = 1}^{m} y_i^{h_i} .
          \end{equation}

          Where
          \begin{equation}
            y_i=  \frac{X_{i}^{(2)}}{X_{i}} ,
            \end{equation}

                 and
                 \begin{equation}
                  y_i^{(n)} = \frac{X_{i}^{(2n)}}{ {X_{i}}} .
                  \end{equation}

          Let \begin{equation}
           \lambda_{i_j}= \frac{X_{i}^{(2j)}}{X_{i}} ,
           \end{equation}
                       we get ordinary differential equations those are
               \begin{equation}
               X_{i}^{(2j)} - \lambda_{i_j} X_{i}=0,      \label{eq-111}
              \end{equation}
               and
              \begin{equation}
              \sum_{k=0}^{(n)} \binom {n} {k}  X_{m}^{(2n-2k)}  K_m^{(2k)}  = 0 .  \label{eq-112}
              \end{equation}
              Where   $ X_{m}^{(0)} = X ,      $
                $ (2j) $ is the order of dervative, and   $ k_m= \sum_{i =1}^{m-1} \lambda_{i_1}.  $
                Let solve \eqref{eq-111} if  $ j=1, $    we get
                 \begin{equation}
                 X_{i}^{(2)} - \lambda_{i_1} X_{i}=0.  \label{eq-113}
                \end{equation}
                 If $ \lambda_{i_1} < 0, $ we have
                 \begin{equation}
                 X_i (x_i)= a_i \cos \sqrt{\lambda_{i_1}} x_{i}  + b_i \sin \sqrt{\lambda_{i_1}} x_{i}.   \label{eq-1}
                  \end{equation}
                for  $ j > 1 . $  we will let $ \lambda_{i_j} = (\lambda_{i_1})^j   $ in   \eqref{eq-111}, then the solutions of
                this forms have two terms. First one is  \eqref{eq-1}, and second terms in which their constant forced to be
                zeros. The final result is \eqref{eq-1}.
                Similary if $ \lambda_{i_1} > 0,$    we  have
                 \begin{equation}
                 X_i (x_i) = a_i \cosh \sqrt{\lambda_{i_1}} x_{i}  + b_i \sinh \sqrt{\lambda_{i_1}} x_{i}  .
                 \label{eq-2}
                  \end{equation}
                  If $ \lambda_{i_1} = 0,$   we have
                   \begin{equation}
                   X_i(x_i) = a_i   + b_i x .
                      \label{eq-3}
                    \end{equation}

                The solution of \eqref{eq-112}, elementary differetial equation yields
                 \begin{equation}  \sum_{k=0}^{n} ( X_{m}^{2} + k_m )^{(n)}=0 .
                 \end{equation} Where $ X^{(0)}=X . $  This gives
                   \begin{equation}
                   \sum_{k=0}^{n} ( s^{2} + k_m )^{n}=0 .
                   \end{equation}
                  If $ k_m < 0 ,$  we have
                \begin{equation}
                 X_m(x_m)= \sum_{r=1}^{2n} x_m^{r-1}[ c_r \cosh \sqrt{- k_m} x_m  +  d_r \sinh \sqrt{ - k_m} x_m ].  \label{eq-4}
                \end{equation}
                 If $ k_m > 0 ,$  we have
                  \begin{equation}
                  X_m(x_m)= \sum_{r=1}^{2n} x_m^{r-1}[ c_r \cos  \sqrt{ k_m} x_m)  +  d_r
                 \sin  \sqrt{k_m}x_m ].  \label{eq-5}
                \end{equation}
                Finally if $ k_m = 0 ,$  we have
                \begin{equation}
                  X_m(x_m)= \sum_{r=1}^{2m} c_r x^{r-1} .   \label{eq-6}
                \end{equation}
                 The final solutions depends on choice of  $ \lambda_{i_1} $ for each $ i=1,2,.........m-1. $  We have
                   \begin{equation}
                    u(x)=\prod_{i = 1}^{m-1} X_i(x_i)   X_m (x_m) .    \label{eq-7}
                    \end{equation}
                Remark: We may have different signs for $ \lambda_{i_1} . $

                \section{Boundary Value Problem Unbounded Domain}
    We consider the Boundry value problem (BVP) replace $ x=(x_1, x_2, ........, x_{m-1},x_m) $ by  $ (x,x_m)=(x_1, x_2, ........, x_{m-1},x_m).  $

    \begin{equation}
    \triangle_m^{n} u(x,x_m) =0 , \quad   - \infty < x_i <  \infty;  \quad  L < x_m < \infty ,
    \end{equation}

 \begin{equation}
  u (x, 0)=f(x) ,  \  \  \mid u(x,x_m)\mid  \textless \  M  .
   \end{equation}
Write  $ u(x,x_m) =  \prod_{i = 1}^{m-1} X_i(x_i)   X_m (x_m). $
Apply  gererlazation of sepration of varibles to get

                \begin{equation}
               X_{i}^{(2j)} - \lambda_{i_j} X_{i} =0 ,
              \end{equation}
              and
             \begin{equation}
             \sum_{k=0}^{n} ( X_{m}^{2} + k_m )^{(n)}=0 .
                 \end{equation}
                 Where  $ X^{(0)}=X .  $
                 From the gereraliztion we need only $ j=1 $ i.e.
                 \begin{equation}
               X_{i}^{(2)} - \lambda_{i_1} X_{i}=0,
              \end{equation}
                      and by letting  $ \lambda_{i_1}= - w_{i}^2 .$
                     We get   \begin{equation}
               X_{i}^{(2)} + w_{i}^{2} X_{i}=0 .
              \end{equation}
              The solutions is
        \begin{equation}
                      X_i (x_i)= a_i \cos w_i x_{i}  + b_i \sin w_i x_{i} .
                    \end{equation}
                    The solutions for the mth variable is in the form

                    \begin{equation}
                    X_m(x_m)= \sum_{r=1}^{2n} x_m^{r-1}[ q_r e^{ - \sqrt{ - k_m} x_m}  +  f_r  e^{ \sqrt{ - k_m} x_m}  ].
                    \end{equation}

                    The final solutions
                  \begin{equation}
                   u(x,x_m) =  \prod_{i = 1}^{m-1} [a_i \cos w_i x_{i}  + b_i \sin w_i x_{i} ]   \sum_{r=1}^{2n} x_m^{r-1}[ q_r e^{ - \sqrt{ - k_m} x_m}  +  f_r  e^{ \sqrt{- k_m} x_m}  ] .
                   \end{equation}
                 Where $ k_m= - \sum_{i =1}^{m-1} w_i^{2} . $  Since  $ w_i  > 0  , $   and  $ x_m \rightarrow  \ \infty ,   $  we must have $ f_r = 0 .$  It follows that
              \begin{equation}
                u(x,x_m) =  \int_{(0,\infty)^{m-1} }  \sum_{r=1}^{2n}  \prod_{i = 1}^{m-1} [A_{i_r} x_m^{r-1} \cos w_i x_{i}
              +  x_m^{r-1}  B_{i_r} \sin w_i x_{i} ] [ e^{ - \sqrt{ - k_m} x_m}   ] \   dw_1 dw_2 \cdots dw_{m-1} .
           \end{equation}

                    Since    $ u (x, L)=f(x) ,  $  then
                           \begin{equation}
                            f(x)  =    \int_{(0,\infty)^{m-1} }\sum_{r=1}^{2n}   \prod_{i = 1}^{m-1} L^{r-1} e^{ - \sqrt{- k_m} L }    [ A_{i_r}   \cos w_i x_{i}  + A_{i_r}  \sin w_i x_{i} ]   \  d w_1 dw_2 \cdots  d w_{m-1} .
                            \end{equation}
                              Where $ w=(w_1, w_2, ........, w_{m-1}).  $
                    Thus, we have
                   \begin{equation}
                   A_{i_r} (w) =   \frac{ \sum_{r=1}^{2n} L^{1-r} e^{ \sqrt{ - k_m} L }}{  \pi^{m-1}} \int_{(-\infty,\infty)^{m-1} }  f(x)     \prod_{i = 1}^{m-1}  \cos w_i x_i \  d x_1 d x_2  \cdots   d x_{m-1} .
                 \end{equation}
                     \begin{equation}
                     B_{i_r} (w) = \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ \sqrt{- k_m} L } }{ \pi^{m-1}}  \int_{(-\infty,\infty)^{m-1} }  f(x)     \prod_{i = 1}^{m-1}  \sin w_i x_i \   d x_1 d x_2 \cdots      d x_{m-1} .
                     \end{equation}
                   Where      $ r=1,2,...........,2n .$

\newcommand{\mathbb{N}}{\mathbb{N}}

            \section{Application of Convolution Theorem}
             The case 1-harmonic 2 dimenstion has been done in \cite{Haberman} , we extended to  n-harmonic 2 dimenstion.

             Let us the result of previous problem when m=2,
              \begin{equation}
               u(x_1, x_2) =  \int_0^{\infty}  \sum_{r=1}^{2n}   [A_{1_r} x_2^{r-1}
              \cos w_1  x_{1}  +  x_2^{r-1}  B_{i_r} \sin w_1 x_{1} ] [ e^{ - w_1  x_2}   ] \ dw_1 .
               \end{equation}
            where
           \begin{equation}   A_{1_r} (w) =   \frac{ \sum_{r=1}^{2n} L^{1-r} e^{ w_1 L }}{  \pi } \int_{-\infty}^{\infty}  f(x_1)      \cos w_1 x_1  \ dx_1 .
            \end{equation}
               \begin{equation}
                B_{1_r} (w) = \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1  L } }{ \pi } \int_{-\infty}^{\infty}   f(x_1)     \sin   w_1 x_1 \  dx_1.
                \end{equation}
                After we plug the coefficients into the general formula, we get
                     \begin{equation}
                      u(x_1, x_2) =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1} L }{ \pi}  \int_{-\infty}^{\infty}  \int_{0}^{\infty}   \sum_{r=1}^{2n}   [  x_2^{r-1}
              \cos w_1  x_1 \cos w_1 z_1   f(z)  .
                 \end{equation}
              \begin{equation}   +   x_2^{r-1}   \sin w_1 x_{1} \sin w_1 z_1 f(z_1)   ] [ e^{ - w_1 x_2}   ] \ dw_1  dz_1
              \end{equation}
                 \begin{equation}  =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1 L } }{ \pi}  \int_{-\infty}^{\infty} \int_0^{\infty}   \sum_{r=1}^{2n}   [ x_2^{r-1} e^{ - w_1 x_2}
               \cos w_1 (z_1 -x_1)   f(z_1)  ]  \ dw_1  dz_1 .
                \end{equation}

               Let \begin{equation}
                g(x_1, x_2)=  \int_0^{\infty}  \sum_{r=1}^{2n}   \cos w_1 (x_1) x_2^{r-1} e^{ -  w_1 x_2}  \ dw_1
                \end{equation}
             \begin{equation}  =  \sum_{r=1}^{2n}    \int_0^{\infty}  \frac{e^{i w_1 x_1} +e^{-i w_1 x_1 }}{2}  e^{ - w_1 x_2} x_2^{r-1}   \ dw_1
             \end{equation}
           \begin{equation}    =  \sum_{r=1}^{2n}   x_2^{r-1} e^{ - w_1 x_2}  \frac{1}{2} \left( \frac{1}{x_2-ix_1} + \frac{1}{x_2+i x_1} \right)  .
            \end{equation}
                \begin{equation}
                        = \sum_{r=1}^{2n}   x_2^{r-1} e^{ - w_1 x_2}  \frac{x_2}{x_2^{2}+x_1^{2}} .
                \end{equation}
           Define the convolution theorem in the form of
               \begin{equation}
                u(x_1,x_2) =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1 L } }{ \pi }  \int_{-\infty}^{\infty}  f(z_1)  g(x_1-z_1 , x_2 ) dz_1.
                \end{equation}
           Let
   \begin{equation}
             f(z_1)= \left \{ \begin{array}{rcl} 1  & \mbox{if} & z \textgreater  \ 0   , \\  0   &       \mbox{if} &  z \textless \ 0   . \\ \end{array}  \right \}
    \end{equation}

           Now apply the theorem to get

           \begin{equation}
           u(x_1,x_2) =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1 L } }{ \pi}   \int_{-\infty}^{\infty}   \sum_{r=1}^{2n}   x_2^{r-1} e^{ - w_1 x_2}  \frac{x_2}{x_2^{2}+(x_1 - z_1)^{2}}   dz_1.
          \end{equation}

           \begin{equation}
            =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1 L } }{ \pi}    \sum_{r=1}^{2n}  x_2^{r-1} e^{ - w_1 x_m}  \int_{-\infty}^{\infty}   e^{ - \sqrt{ k_m} x_m}  \frac{x_m}{x_m^{2}+(x_i - z_i)^{2}}   dz_i.
              \end{equation}

                   \begin{equation}
                    =  \frac{\sum_{r=1}^{2n}  L^{1-r} e^{ w_1 L } }{ \pi }    \sum_{r=1}^{2n}
              x_2^{r-1} e^{ - w_1 x_2}  \left [  \frac{\pi }{ 2} + \tan^{-1} \left(
             \frac{x_1}{x_2} \right) \right ]  .
              \end{equation}

             \section{Parabolic Version Of n-Harmonic Equation}
             The parabolic version of n-harmonic equation is
           \begin{equation}
 \alpha \triangle_m^{n} u(x, t) = u_t (x, t).
   \end{equation}

 Where $ (x,t)=(x_1,x_2,........x_m, t) , $ $ t >0 ,$  $ \alpha \in \mathbb{R} , $ and  $ (x,t) \in \mathbb{R}^m \times (0, \infty ) . $
            The between  difference between the eillptic version of n harmonic equation and prarbolic one  is similar to difference between heat equation and  1-harmonic equation. The solution can presented as $  u(x)=\prod_{i = 1}^{m} X_i(x_i) T(t) .$   $X_i(x_i)  $ has same solutions as the one in n harmonic equation.   For $ T(t) $ we have
           $$ T^{'} (t)- \alpha k_{m+1}T(t)=0 .$$
           The solutions is $$ T(t)= A e^{ \alpha t k_{m+1}  } . $$

             \section{Hyperbolic Version Of n-Harmonic Equation}
             The hyperbolic version of n-harmonic equation is
       \begin{equation}
   \beta^2 \triangle_m^{n} u(x, t) = u_t (x, t).
 \end{equation}  Where $ (x,t)=(x_1,x_2,........x_m, t) , $ $ t >0 ,$  $ \alpha \in \mathbb{R} , $ and  $ (x,t) \in \mathbb{R}^m \times (0, \infty ) . $
             The  difference between the eillptic version of n harmonic equation and hyperbolic one is similar to difference between wave equation and  1-harmonic equation. The solution can be written  as $  u(x)=\prod_{i = 1}^{m} X_i(x_i) T(t) .$   $X_i(x_i)  $ has same solutions as the one in n harmonic equation.   For $ T(t) $ we have
           $$ T^{''} (t)- \beta^2 k_{m+1} T(t)=0. $$
           The solutions is  $$ T(t)= C e^{ \beta \sqrt {k_{m+1}} t }  +  D  e^{ - \beta t \sqrt {k_{m+1} t }}   . $$

        \section{Conclusion}

         Sepations of variables is powerful method that assist to reduce the diffcuilty of partial differetion equation problems and solve them. The n-harmonic m dimension has been generlaized. Moreover, unbounded boundary value problem has been solved and applications of its convolution. Hyperbolic version and parabolic version are dicussed.

\addcontentsline{toc}{section}{REFERENCES}

\singlespacing