\section{Introduction and Statement of the Results}\label{Introduction}

	This article has been motivated by the following open problem: how can we construct a harmonic map explicitly?

	Our aim is to construct harmonic maps between Riemann surfaces, when the curvature in the target is a negative constant, say $-1$.  For that, we develop some of the machinery introduced in \cite{FotDask}, which connects the harmonic map problem with elliptic versions of the sinh-Gordon and sine-Gordon equation.

	Let $w=w(x,y)$ be a solution of the sinh-Gordon equation \begin{equation}\label{sinh-Gordon eqn}
		\Delta w=2 \sinh(2w)
	\end{equation}
	and let  $u=u(z, \bar{z})$ be a solution of the Beltrami equation
	\begin{equation}\label{Beltrami eqn}
		\frac{\partial_{\bar{z}}u}{\partial_{z}u}=e^{-2w},
	\end{equation}
	where $z=x+iy$. Then, $u$ is a harmonic map  between Riemann surfaces, when the curvature of the target is $-1$ \cite[Theorem 1, Corollary 2]{FotDask}.   In this case, \textit{we say that $u=R+iS$ corresponds to $w$.} From now on, consider the target  surface fixed and realized as the upper half hyperbolic plane. Using a specific coordinate system  on the domain (see Section \ref{sec:Preliminaries}  for more details),  the harmonic map equation under study becomes
	\begin{equation}\label{eq:HarmHyper}
		\frac{\partial_z u \partial_z\bar{u}}{S^2}=1.
	\end{equation}

	In \cite{FotDask}, it is proved that the system
	\begin{align}
		\partial_x w-\partial_y \theta&=-2\sinh w \sin \theta \label{Back1intro}\\
		\partial_y w+\partial_x\theta&=-2\cosh w \cos \theta\label{Back2intro},
	\end{align}
	is a \textit{B{\"a}cklund transformation} that connects a solution $w=w(x,y)$ of the sinh-Gordon equation (\ref{sinh-Gordon eqn}) and a solution $\theta=\theta(x,y)$ of the sine-Gordon equation
	\begin{equation}	\label{sin-Gordon eqn}
		\Delta \theta=-2 \sin(2\theta).
	\end{equation}
	If $(w, \theta)$ satisfy (\ref{Back1intro})-(\ref{Back2intro}), then \textit{we say that the pair $(w,\theta)$ is in the {class $(BT)$}}. From now on we assume that $(w,\theta)\in (BT)$.

	Set
	\begin{align*}I_1&=I_1(x)=\int_0^x\cosh w(t,0)\sin\theta(t,0)dt, \\
		I_2&=I_2(x,y)=\int_0^y\sinh w(x,s)\cos\theta(x,s)ds\\
		I_3&=I_3(x)=\int_0^x e^{2 I_1 (t)}\cosh w(t,0)\cos\theta(t,0)dt,\\
		I_4&=I_4(x,y)=e^{2 I_1 (x)}\int_0^y e^{2 I_2 (x,s)}\sinh w(x,s)\sin\theta(x,s)ds.
	\end{align*}
	We first prove the following result.
	\begin{proposition}\label{Theorem2} Define the function $S$ by
		\begin{align*}
			S(x,y)=S(0,0)e^{2(I_1+I_2)},
		\end{align*}
		and the function $R$ by
		\begin{align*}
			R(x,y)=R(0,0)+2S(0,0)(I_3-I_4).
		\end{align*}
		Then,
		\[
		u(x,y)=R(x,y)+iS(x,y)
		\]
		is a harmonic map that corresponds to $w$.
	\end{proposition}
	Therefore, given a pair of functions $(w,\theta)$ that satisfies the B{\"a}cklund transformation, there is an implicit formula for a harmonic map to the hyperbolic plane that involves the integrals $I_1, I_2, I_3, I_4$.  We provide an example of a harmonic map, first studied in \cite[Section 7]{FotDask}, that we recover by the algorithm of Proposition \ref{Theorem2}.

	We next characterize the solutions $w$ of the sinh-Gordon equation that are of the form
	\begin{equation}\label{separable}
		w(x,y)=2\arctanh({F(x)}{G(y)}),
	\end{equation}
	for some  non-constant functions $F$ and $G$. These functions, in general, turn out to be elliptic (only in certain cases they reduce to elementary functions).  We also prove that the initial value problem $\Delta w=2\sinh(2w)$, $\partial_{y}w(x,0)=0$, admits a solution of the form   (\ref{separable}), see Proposition \ref{Kenmotsu} for a class of such solutions. In this case,  we determine $\theta$ by considering a subclass of $(BT)$: \textit{we say that $(w,\theta)\in (BT_0)$ when the pair is in $(BT)$ and  $\partial_{y}w(x,0)=0$,  $\theta(0,0)=\frac{\pi}{2}$.} For this selection of $(w, \theta)$, and $w$ a solution of the form (\ref{separable}), we apply Proposition \ref{Theorem2} in order to prove our first main result,  which provides an explicit formula for a harmonic map that corresponds to $w$. This is a special case of the more general problem of finding all solutions of the sinh-Gordon equation, and the corresponding harmonic maps. Thus, it is a step (although small) towards the solution of the harmonic map problem.

	Next, we state our first main result.

	\begin{theorem}\label{new class} Assume that {$w(x,y)=2\arctanh({F(x)}{G(y)})$,  $b(x)=\frac{F^{\prime}(x)}{2F(x)}$} and $(w,\theta)\in (BT_0)$.
		Define the function $S$ by
		\begin{equation*}
			S(x,y)=S(0,0) \frac{e^{2X(x) }(\sin{\theta(x,y)}+b(x))}{1+b(x)},
		\end{equation*}
		and the function $R$ by
		\begin{equation*}
			R(x,y)=R(0,0)+S(0,0)\frac{e^{2X(x)}\cos{\theta(x,y)}}{1+b(x)},
		\end{equation*}
		where $X=X(x)=\int_0^x \cosh w(t,0)dt$.
		Then  $u(x,y)=R(x,y)+iS(x,y)$ is a harmonic map that corresponds to $w$.
	\end{theorem}

	Using the machinery of Theorem \ref{new class}, we construct an entirely new harmonic map to the hyperbolic plane. At this point, it is worth mentioning that the initial condition $\theta(0,0)=\frac{\pi}{2}$ was chosen just to facilitate computations; any other initial condition for $\theta$ at the origin would provide a corresponding harmonic map.

	In \cite{FotDask}, the approach pursued by the authors on the construction of harmonic maps, is to first consider a one-soliton (i.e., $w=w(x)$) solution of the sinh-Gordon equation and then solve the Beltrami equation (\ref{Beltrami eqn}) to find a harmonic map that corresponds to $w$. However, if $w(x,y)=2\arctanh({F(x)}{G(y)})$, and both $F, G$ are elliptic functions, then this Beltrami equation turns out to be equivalent to a Riccati equation. Therefore, explicit construction of harmonic maps by this approach is rather implausible. Instead, in place of the Beltrami equation, we use the B{\"a}cklund transformation combined with Theorem \ref{new class}, to provide an explicit formula for a harmonic map in terms of the functions $w$ and $\theta$.

	Next, instead of considering solutions to the sinh-Gordon equation first, we investigate the case of solutions to the sine-Gordon equation as a starting point. In particular, we prove that if $\theta=\theta(x)$ is a solution to the sine-Gordon equation (\ref{sin-Gordon eqn}), then the function $w$ determined by the B{\"a}cklund transformation is of the form (\ref{separable}).  We then provide an explicit formula for the corresponding harmonic map.

	Our second main result is the following.

	{\begin{theorem}\label{thetaxthm}
			Let $\theta=\theta(x)$ be a solution to the sine-Gordon equation and let $w$ be the associated solution to the B{\"a}cklund transformation. Set $a = 2\cos\theta(0) + \theta'(0)$, $b = 2\cos\theta(0) - \theta'(0)$, and $k\in \mathbb{C}$  such that
			$\tanh \frac{w(0,0)}{2}=- {\sqrt{\frac{a}{b}}\tan(k)}$. Then, the harmonic map $u=R+iS$ that corresponds to $w$ is given by
			\begin{align*}
				R(x,y)&= R(0,0)+S(0,0)\cosh^2\frac{w(0,0)}{2} J\left(\frac{2\cos\theta(x) - \theta'(x)}{b}\right)\\
				&+ 2S(0,0)\frac{\sin\theta(x)(\cos(\sqrt{ab}y + 2k) - \cos(2k))}
				{2\cos\theta(0)\cos(2k) - \theta'(0)}, \\
				S(x,y)&=S(0,0)\frac{2\cos\theta(x)\cos(\sqrt{ab}y + 2k) - \theta'(x)}
				{2\cos\theta(0)\cos(2k) - \theta'(x)},
			\end{align*}
			where
			\[ J(t)=\int_{1}^{t}
			\frac{u^{2} + \tanh^2\frac{w(0,0)}{2}}{u^{2}}
			\frac{bu^{2} + a}
			{\sqrt{2(8-ab)u^{2} - a^{2} - b^{2}u^{4}}}du.
			\]
		\end{theorem}

		As an application of Theorem \ref{thetaxthm}, we construct an entirely new example of a harmonic map to the hyperbolic plane.

		There are  many examples of harmonic diffeomorphisms, see for example \cite{C-T, Li-Ta2,S-T-W, Wolf1, Wolf2, X-A}. For instance, a minimal surface in $\mathbb{H}^2\times \mathbb{R}$ projects to a harmonic map to $\mathbb{H}^2$. The study of harmonic diffeomorphisms between Riemannian surfaces is central to the theory of harmonic maps (see for example \cite{S-Y, Hitchin} and observe that a twisted harmonic map, is a harmonic map in a local sense and that each harmonic map induces a solution of Hitchin self-duality equations). The most interesting case is when the surfaces are of constant curvature (see for example \cite{A-R,Kal,Li-Ta2,M2,Minsky, P-H,S-T-W,  Wolf1, Wolf2, Wolf3} and the references therein). In \cite{Han, Minsky, Wolf1, Wolf2, Wolf3}, the geometry of harmonic maps between hyperbolic surfaces is studied, while in \cite{M} the author proves a conjecture of R. Schoen on  harmonic diffeomorphisms between hyperbolic spaces.   For further results about the harmonic map equations see also  \cite{Han, H-T-T-W, Minsky, S-Y, Wolf1, Wolf2, Wolf3}. There is a similar analysis of the Wang equation, which is applied in the study of affine spheres in \cite{DW, Loftin, OT21} and \cite{Li19} explains the relation between the harmonic map and Wang equations with CMC surfaces and hyperbolic affine spheres in $\mathbb{R}^3$. Note finally that \cite{Tam19} extends \cite{DW} and explains the relation between harmonic map equation and anti-de Sitter geometry, see also \cite{BS20}.

		Harmonic maps to surfaces of constant curvature are closely related to the elliptic sinh-Gordon equation. The sinh-Gordon  and sine-Gordon equations have many applications and they have both been the subject of extensive study, \cite{F-P,H}. Furthermore, the sinh-Gordon equation was  crucial to the breakthrough work \cite{Wente} on the Wente torus. There is a  close relation to the theory of constant mean curvature surfaces as well, \cite{Joa,K}.

		The outline of this paper is as follows. In Section \ref{sec:Preliminaries} we  recall some basic notations and review certain known results. In Section \ref{sec:Harmonic maps to the upper half-plane model} we prove Proposition \ref{Theorem2} and in Section \ref{sec:The sinh-Gordon equation} we construct new solutions of the sinh-Gordon equation that are of the form $w(x,y)=2\arctanh({F(x)}{G(y)})$. In  Section \ref{sec:A new family of harmonic maps},  for this class of $w$, we find the explicit formula of a harmonic map that corresponds to $(w,\theta)\in (BT_0)$. Finally, in Section \ref{app}, we discuss one-soliton solutions of the sine-Gordon equation and the construction of the corresponding harmonic map.

		\vspace{10pt}\noindent\textbf{Acknowledgments.}
		The second author is supported by the Hellenic Foundation for Research and Innovation, Project HFRI-FM17-1733.

		\section{Preliminaries}\label{sec:Preliminaries}

		In this section we discuss some necessary preliminaries.

		Let $u : M \rightarrow N$ be a map between
		Riemann surfaces $(M, g)$, $(N, h)$. The map $u$ is locally represented by
		$u = u(z, \bar{z}) = R + iS$, where $z=x+iy$. From now on, we use the standard notation:
		\[
		\partial_z=\frac{1}{2}(\partial_x-i\partial_y), \quad \partial_{\bar{z}}=\frac{1}{2}(\partial_x+i\partial_y).
		\]

		Recall that isothermal coordinates on a Riemannian manifold are local coordinates where the metric is conformal to the Euclidean metric. The existence of  isothermal coordinates on an arbitrary surface with a real analytic metric is a well-known fact, first proved by Gauss (see for instance \cite[Section 8, p.396]{J}).
		Consider an isothermal coordinate system $(x, y)$ on $M$ such that
		\[g =  e^{f(z,\bar{z})}|dz|^2,\]
		where $z = x + iy$, and an isothermal coordinate system $(R, S)$ on
		$N$ such that
		\begin{equation*}\label{targetconf}h =  e^{F(u,\bar{u})}|du|^2,
		\end{equation*}
		where $u = R + iS$.
		The Gauss curvature on the target is given by the formula
		\[
		K_N(u, \bar{u})= -\frac{1}{2}\Delta F(u, \bar{u}) e^{-F(u, \bar{u})}.
		\]

		In isothermal coordinates,  a map between two surfaces is harmonic if and only if it satisfies
		\begin{equation*}
			\partial_{z\bar{z}}u + \partial_{u}F(u,\bar{u}) \partial_{z}u\partial_{\bar{z}}u=0,
		\end{equation*}
		see \cite[Section 8, p.397]{J}. Note that this equation only depends on the conformal structure of $N$.

		The Hopf differential of $u$ is given by
		\begin{equation*}\label{eq:Hopf_theorem}
			\Lambda(z) dz^2=\left( e^{F(u,\bar{u})}\partial_{z} u \partial_{z} \bar{u}\right)dz^2.
		\end{equation*}
		The following proposition is well-known.
		\begin{proposition}[\cite{J}]\label{prop:Hopf}
			A necessary and sufficient condition for  $u$ to be a harmonic map, is that the Hopf differential of $u$ is holomorphic, i.e.
			\begin{equation}\label{eq:Hopf}
				e^{F(u,\bar{u})} \partial_z u    \partial_z\bar{u}= e^{-\mu (z)},
			\end{equation}
			where $\mu (z)$ is a holomorphic function.\end{proposition}

		A fundamental result in the theory of harmonic maps, which connects the harmonic map problem with the sinh-Gordon equation, is the following. \begin{proposition}[\cite{Minsky}, \cite{Wolf1}] \label{PropMisnkyWolf}
			Let $u:M\rightarrow N$ be a harmonic map. Then, it satisfies the Beltrami
			equation
			\begin{equation*}\label{Beltrami eqnOLD}
				\frac{\partial_{\bar{z}}u}{\partial_{z}u}=e^{-2w+i\phi},
			\end{equation*}
			and $\phi$ is a harmonic function, i.e. $\phi_{z\bar{z}}=0$. Furthermore, if $\psi$ is the conjugate harmonic function to $\phi$, then
			\begin{equation*}
				K_N=-\frac{2w_{z\bar{z}}}{\sinh 2w}e^{\psi},
			\end{equation*}
			where $K_N$ is the curvature of the target manifold $N$.
		\end{proposition}

		On the domain surface $M$, we can choose a specific coordinate system in order to considerably facilitate the calculations. This \textit{specific}  system is defined by the conformal transformation
		\begin{equation*}\label{eq:specific}
			Z= \int e^{-\mu(z)/2} \, dz.
		\end{equation*}
		In this specific system the above equations simplify by considering $\mu(z)=0$.
		In particular, the harmonic map condition (\ref{eq:Hopf}) becomes \begin{equation}\label{Hopfnew}
			e^{F(u, \bar{u})}\partial_{z} u\partial_{z}\bar{u}=1.
		\end{equation}

		It follows from Proposition \ref{PropMisnkyWolf} that the harmonic map $u$  determined by (\ref{Hopfnew}) satisfies the Beltrami equation (\ref{Beltrami eqn}), where in turn the function $ w(z, \bar{z})$  occurring in the Beltrami coefficient} satisfies the sinh-Gordon equation (\ref{sinh-Gordon eqn}).  Conversely, as we have already mentioned, if $w$ is a solution of the sinh-Gordon equation (\ref{sinh-Gordon eqn}) and  $u$ is a solution of the Beltrami equation (\ref{Beltrami eqn}), then, $u$ is a harmonic map,  \cite{FotDask}.
	As a result, there is a classification of harmonic diffeomorphisms via the classification of the solutions of the sinh-Gordon equation. For further details about the association of the harmonic map problem to the sinh-Gordon equation, we refer to \cite[Section 3]{FotDask}.
	From now on, we   use the aforementioned specific coordinate system.

	We assume that the target is the upper half hyperbolic plane
	\[
	\mathbb{H}^{2}=\{R+iS: \; R\in \mathbb{R},\; S>0\},
	\] equipped with the hyperbolic metric
	\[h=\frac{dR^2+dS^2}{S^2}\]
	of curvature $-1$.
	This is by no means restrictive: it is well known that if $u$ is a harmonic map then for every conformal function $f$, the map $U=f(u(z, \bar{z}))$ is also a harmonic map, and the metric on the target is still of constant curvature $-1$. As mentioned earlier, in the specific coordinate system, the map $u=R+iS$ is harmonic if
	(\ref{eq:HarmHyper}) holds true, or equivalently if the following system is satisfied:
	\begin{equation}\label{Cartesian1}
		\partial_{x}R \;\partial_{y}R + \partial_{x}S\;\partial_{y}S = 0
	\end{equation}
	\begin{equation}\label{Cartesian2}
		\frac{(\partial_{x}R)^{2} + (\partial_{x}S)^{2} - (\partial_{y}R)^{2} - (\partial_{y}S)^{2}}
		{S^{2}} = 4.
	\end{equation}

	Finally, as we have described, our approach has as its starting point the B{\"a}cklund transformation (\ref{Back1intro})-(\ref{Back2intro}). For computational reasons, let us point out that it is more convenient to verify that a pair $(w, \theta)$ is a solution if the following system holds:
	\begin{align*}
		\partial_x W\;(1+\Theta^2) - \partial_y\Theta\;(1-W^2)&=-4W\Theta,\\
		\partial_y W\;(1+\Theta^2) + \partial_x\Theta\;(1-W^2)&=-(1+W^2)(1-\Theta^2),
	\end{align*}
	where $W=\tanh\frac{w}{2}$ and $\Theta=\tan\frac{\theta}{2}$.

	\section{Proof of Proposition \ref{Theorem2}}\label{sec:Harmonic maps to the upper half-plane model}

	In this section we consider $(w,\theta)\in (BT)$.  In other words, $w$ and $\theta$ are related by the B{\"a}cklund transformation.  Our aim is to construct a harmonic map that corresponds to $w$. As usual, the target is assumed to be the upper half-plane equipped with the hyperbolic metric.

	Our first goal is to prove the following result, which plays a key role for the proof of Proposition \ref{Theorem2}.

	\begin{lemma}Assume that $(w,\theta)\in (BT)$. Define the functions $R$ and $S$ by
		\begin{align}
			\partial_xS&=2S\cosh w\sin \theta \label{Sx}\\
			\partial_yS&=2S\sinh w\cos \theta \label{Sy}\\
			\partial_xR&=2S\cosh w\cos \theta \label{Rx}\\
			\partial_yR&=-2S\sinh w\sin \theta \label{Ry}.
		\end{align}
		Then,
		\[
		u(x,y)=R(x,y)+iS(x,y)
		\]
		is a harmonic map to the hyperbolic plane that corresponds to $w$.
	\end{lemma}
	\begin{proof}First, note that the existence of the functions $R$ and $S$ is ensured by the B{\"a}cklund transformation of $w$ and $\theta$. Indeed, the conditions $\partial_{xy}^2 S=\partial_{yx}^2 S$, $\partial_{xy}^2 R=\partial_{yx}^2 R$ are equivalent to  the B{\"a}cklund transformation equations (\ref{Back1intro})-(\ref{Back2intro}).

		Consider \[
		u=R+iS,
		\]
		and observe that
		\[
		\partial_zu=\frac{1}{2}((\partial_xR+\partial_yS)+i(\partial_xS-\partial_yR)).
		\]
		Then (\ref{Sx})-(\ref{Ry}), imply that
		\begin{equation}\label{uz}
			\partial_zu=S(\cosh w+\sinh w)(\cos\theta+i\sin\theta)=Se{^w} e^{i\theta}.
		\end{equation}
		Similarly, we compute
		\begin{equation}\label{uzbar}
			\partial_{\bar{z}}u=Se^{-w} e^{i\theta}.\end{equation}
		By (\ref{uz}) and (\ref{uzbar}), it follows that
		\[
		\frac{\partial_{\bar{z}}u}{\partial_zu}=e^{-2w}.
		\]
		Furthermore, by (\ref{Hopfnew}), the conformal factor in the target metric is equal to \[e^{F(R,S)}=\frac{1}{\partial_zu\partial_z\bar{u}}=\frac{1}{S^2}.\]
		Thus, $u$ is a harmonic map to the hyperbolic upper half hyperbolic plane that corresponds to $w$.
	\end{proof}

	Given now $w$ and $\theta$, we   find the associated harmonic map $u=R+iS$. More precisely, we provide the implicit formulas for $R$ and $S$ that involve the integrals $I_1, I_2, I_3, I_4$.

	\begin{proof}[End of the proof of Proposition \ref{Theorem2}]
		{{Our task is now to} solve {the }system (\ref{Sx})-(\ref{Ry}). {We} begin with (\ref{Sx})} for $y=0$, which yields
		\[
		\frac{\partial_xS}{S}(x,0)=2\cosh w(x,0)\sin\theta(x,0).
		\]
		Integrating, we obtain
		\begin{equation*}\label{S(x,0)}
			S(x,0)=S(0,0)e^{2\int_0^x \cosh w(t,0)\sin\theta(t,0)dt}.
		\end{equation*}
		So, (\ref{Sy}) implies that
		\begin{equation*}
			S(x,y)=S(0,0)e^{2\int_0^x \cosh w(t,0)\sin\theta(t,0)dt+2\int_0^y \sinh w(x,s)\cos\theta(x,s)ds}.
		\end{equation*}
		Thus
		\begin{align*}
			S(x,y)=S(0,0)e^{2(I_1+I_2)}.
		\end{align*}

		In order now to compute $R(x,y)$, we consider (\ref{Rx}) for $y=0$, which yields
		\begin{equation*}\label{R(x,0)}
			R(x,0)=R(0,0)+2\int_0^xS(t,0)\cosh w(t,0)\cos\theta(t,0)dt.
		\end{equation*}
		So, (\ref{Ry}) implies that
		\begin{align*}\
			R(x,y)=R(0,0)&+2\int_0^xS(t,0)\cosh w(t,0)\cos\theta(t,0)dt\notag\\&-2\int_0^yS(x,s)\sinh w(x,s)\sin\theta(x,s)ds,
		\end{align*}
		thus
		\begin{align*}
			R(x,y)=R(0,0)+2S(0,0)(I_3-I_4).
		\end{align*}
	\end{proof}

	We now provide a concrete example of a harmonic map to show how the algorithm in Proposition \ref{Theorem2} can be implemented. In particular, we recover the harmonic map obtained in \cite[Section 7]{FotDask}.
	\begin{example}
		Suppose that $\tanh{\frac{w(x,y)}{2}}=\frac{2y}{\cosh{2x}}$ and consider the case $R(0,0)=0$, $S(0,0)=-\frac{1}{4}$ as in {\cite[p.22]{FotDask}}.  Using the B{\"a}cklund transformation we find that
		\[
		\tan\frac{\theta(x,y)}{2}=\coth{x}.
		\]

		Then, a lengthy calculation yields
		\begin{align*}
			I_1 &= \frac{1}{2} \log{\cosh{2x}},\\
			I_{2}&=\frac{1}{2}\log{\left( 1-\frac{4 y^2}{\cosh^2{2x}}\right)},\\
			I_3 &= -x,\\
			I_4 &=2y^2 \tanh{2x},
		\end{align*}
		{and} therefore,
		\begin{align*}
			R&= \frac{x}{2}+y^2 \tanh{2x} \\
			S&= \frac{y^2}{\cosh{2x}}-\frac{\cosh{2x}}{4}.
		\end{align*}
		One can verify by (\ref{Cartesian1})-(\ref{Cartesian2}) that $u=R+iS$ is a harmonic map to the hyperbolic plane that corresponds to $w$.
	\end{example}

	\section{The sinh-Gordon equation}\label{sec:The sinh-Gordon equation}

	In this section we provide a new family of solutions of the sinh-Gordon equation.

	The motivation for obtaining solutions of the sinh-Gordon equation is their correspondence to harmonic maps between surfaces, as described in Section \ref{sec:Preliminaries}.  In addition, it is exactly the form of the special solutions to the sinh-Gordon equation considered in this section, that will allow explicit expressions of harmonic maps to hyperbolic plane. We therefore postpone our main goal, the construction of harmonic maps, to the next section.

	In \cite[Section 5]{FotDask}, the authors construct harmonic maps when $w$ is a one-soliton solution to (\ref{sinh-Gordon eqn}). In some sense, they obtain solutions to this sinh-Gordon equation that are of the form
	\[
	w(x,y)=2\arctanh({F(x)}).
	\]
	In what follows,  we focus on solutions to the sinh-Gordon equation (\ref{sinh-Gordon eqn}) that are of the form (\ref{separable}), i.e. $w(x,y)=2\arctanh({F(x)}{G(y)})$.
	Without loss of generality, we assume that $F$ and $G$ are non-constant functions. Condition  (\ref{separable}) is equivalent to $b=\frac{\partial_{x}w(x,y)}{2\sinh{w(x,y)}}=b(x)$.
	Note that condition (\ref{separable}) is also equivalent to the differential equation
	\[
	\partial_{xy}^{2}w(x,y)=\partial_{x}w(x,y)\partial_{y}w(x,y)\coth{w(x,y)}.
	\]
	In \cite{Abr}, the authors were lead to this equation for geometric reasons related to CMC surfaces. It is an open problem to relate this equation with the geometry of harmonic maps.

	The basic idea is that under condition (\ref{separable}), we can solve this sinh-Gordon equation by a separation of variables argument. Then, the problem is  reduced to solving certain non-linear ODEs, which  in turn can be solved by the use of elliptic functions.

	We now prove the following result.
	\begin{proposition}\label{w thm}
		{If $w$ is} a solution of the sinh-Gordon equation (\ref{sinh-Gordon eqn})
		of the form
		\begin{equation*}
			w(x,y)=2\arctanh({F(x)}{G(y)}),
		\end{equation*}
		then the functions $F,G$ satisfy the differential equations
		\begin{align}
			(F^{\prime}(x))^2&=AF^4(x)+BF^2(x)+C \label{odef}\\
			(G^{\prime}(y))^2&=-CG^4(y)-(B-4)G^2(y)-A \label{odeg},
		\end{align}
		where $A,B,C$ are arbitrary constants.
	\end{proposition}

	\begin{proof}
		Set $G(y)=\frac{1}{H(y)}$.	 Writing $w=2\arctanh\frac{F(x)}{H(y)}$, we find from {(\ref{sinh-Gordon eqn})} that
		\begin{align}
			&\;\frac{F^{\prime\prime}(x)}{F(x)}\left({H^2(y)}
			-{F^2(x)}\right)+2F^{\prime}(x)^2
			-\frac{H^{\prime\prime}(y)H(y)-2H^{\prime}(y)^2}{H(y)^2}\times \notag \\
			&\times\left({H^2(y)}-{F^2(x)}\right)+2F^2(x)\frac{H^{\prime}(y)^2}{H^2(y)} \label{eq:sinhSOL} \\&=4\left(F^2(x)+H^2(y)\right).\notag
		\end{align}
		Differentiating with respect to $x$ and $y$, we obtain
		\[
		2H(y)H^{\prime}(y)\left( \frac{F^{\prime\prime}(x)}{F(x)} \right)^{\prime}+2F(x)F^{\prime}(x)\left( \frac{H^{\prime\prime}(y)}{H(y)} \right)^{\prime}=0,
		\]
		or
		\begin{equation*}\label{7}
			\frac{1}{F(x)F^{\prime}(x)}\left( \frac{F^{\prime\prime}(x)}{F(x)} \right)^{\prime}=-\frac{1}{H(y)H^{\prime}(y)}\left( \frac{H^{\prime\prime}(y)}{H(y)} \right)^{\prime}=4A=\text{constant}.
		\end{equation*}
		Therefore,
		\begin{align*}
			\left( \frac{F^{\prime\prime}(x)}{F(x)} \right)^{\prime}=4AF(x)F^{\prime}(x), \quad\left( \frac{H^{\prime\prime}(y)}{H(x)} \right)^{\prime}=4AH(y)H^{\prime}(y).
		\end{align*}
		Integrating twice and taking into account (\ref{eq:sinhSOL}), we obtain
		\begin{equation*}\label{8}
			F^{\prime}(x)^2=AF^4(x)+BF^2(x)+C,
		\end{equation*}
		and
		\begin{equation*}
			H^{\prime}(y)^2=-AH^4(y)-(B-4)H^2(y)-C.
		\end{equation*}
		Finally, we find
		\[
		G^{\prime}(y)^2=-CG^4(y)-(B-4)G^2(y)-A.
		\]
	\end{proof}

	An argument similar to the one used in \cite{K}, reveals that the following result, which is a special case of Proposition \ref{w thm}, holds true.
	\begin{proposition}\label{Kenmotsu}
		Let $w_0>0$, and $\alpha,\beta>0$, such that
		\[
		\alpha+\beta=\cosh {w_0}>1.
		\]
		Consider $f(x)$, $g(y)$ such that
		\begin{align*}
			(f')^2&=f^4-4\left(1+\alpha^2-\beta^2\right)f^2+4^2  \alpha^2
			\\
			(g')^2&=g^4-4\left(1+\beta^2-\alpha^2\right)g^2+4^2  \beta^2,
		\end{align*}
		with $f(0)=0$, $f'(0)=-4\alpha$, $g(0)=0$, $g'(0)=-4 \beta$.
		Then the function $w(x,y)$  given by
		\begin{equation*}
			\tanh\frac{w(x,y)}{2}=\tanh\frac{w_0}{2}e^{-\int_0^xf(t)dt}e^{-\int_0^yg(s)ds},
		\end{equation*}
		is such that  $\Delta w=2\sinh(2w)$ and $w(0,0)=w_0$.
	\end{proposition}

	Finally, note that in general, $F$ and $G$ in Proposition \ref{w thm}, and $f$ and $g$ in Proposition \ref{Kenmotsu} respectively, are given in terms of elliptic functions, but reduce to elementary functions for certain choices of the coefficients of the ODEs.  In particular, we observe that the only case where one of the functions $f$, $g$ is elementary is when $|\alpha-\beta|=1$. In this case actually both functions $f$, $g$ end up to be elementary.

	\section{A new family of harmonic maps}\label{sec:A new family of harmonic maps}\

	In this section we  assume that $(w,\theta)\in (BT_0)$  and $w$ is of the form (\ref{separable}). In other words,  the functions $w=2\arctanh(F(x)G(y))$ and $\theta$ are related by the B{\"a}cklund transformation and $\partial_{y}w(x,0)=0$, $\theta(0,0)=\frac{\pi}{2}$. We start with a given solution $w$ of the sinh-Gordon equation and we find a formula for $\theta$ using the B{\"a}cklund transformation. We finally find a formula for the harmonic map $u$ that corresponds to $w$.

	\subsection{B{\"a}cklund transformation} \label{sec:Backlund transform of the sinh-Gordon equation}

	Our first task is to solve the system (\ref{Back1intro})-(\ref{Back2intro}). Suppose that $w$ is a solution of the sinh-Gordon equation, of the form
	\begin{equation}\label{omega1}
		\tanh \frac{w}{2}=F(x)G(y).
	\end{equation}

	Define
	\begin{align}
		X&=X(x)=\int_0^x \cosh w(t,0)dt=\int_0^x \frac{1+F^2(t)G^{2}(0)}{1-F^2(t)G^{2}(0)}dt, \label{Xvar}\\
		Y&=Y(x,y)=\int_0^y \sinh w(x,s)ds=\int_0^y \frac{2F(x)G(s)}{1-F^2(x)G^2(s)}ds\label{Yvar}.
	\end{align}

	The new variable $Y$ appears in the following lemma, which serves as an intermediate step for the formula of the harmonic map in Theorem \ref{new class}, where the other variable $X$ appears. We claim that we can find explicit expressions of $X$ and $Y$ in terms of $x$ and $y$. Assuming for the moment that this claim is true, we complete the construction of $\theta$ as follows, in the special case when {$\partial_{y}w(x,0)=0,$ $\theta(0,0)=\frac{\pi}{2}$.}

	\begin{lemma}\label{BackLemma}
		If $	w(x,y)=2\arctanh({F(x)}{G(y)})$, $b(x)=\frac{F^{\prime}(x)}{2F(x)}$ and $(w,\theta)\in (BT_0)$, then
		\begin{equation*}
			\tan\frac{\theta(x,y)}{2}=\frac{1}{b(x)}\left( \sqrt{b^2(x)-1}\tan({J_1(x)+J_2(x,y)}) -1 \right)
		\end{equation*}
		where{
			\begin{align*}
				J_1=J_1(x)=\arctan \left( \frac{b(x)+1}{\sqrt{b^2(x)-1}}  \right),
		\end{align*}}
		and
		\begin{align*}
			J_2=J_2(x,y)=\sqrt{b^2(x)-1}\;Y(x,y).
		\end{align*}
	\end{lemma}

	\begin{proof}
		By assumption, we have
		\begin{equation}\label{wy0}
			\partial_y w(x,0)=0.
		\end{equation}
		Taking (\ref{Back2intro}) for $y=0$ and applying (\ref{wy0}), gives
		\begin{equation*}
			\partial_x\theta(x,0)=-2\cosh w(x,0)\cos\theta(x,0).
		\end{equation*}
		Given the initial condition $\theta(0,0)=\frac{\pi}{2}$, we can take
		\begin{equation*}\label{tanthetaX0}
			\theta(x,0)=\frac{\pi}{2}.
		\end{equation*}
		Next, observe that since $w(x,y)=2\arctanh({F(x)}{G(y)})$, we have
		\[
		\partial_x w(x,y)=\frac{F^{\prime}(x)}{F(x)}\sinh w(x,y),
		\]
		hence (\ref{Back1intro}) implies
		\begin{equation}\label{thetay}
			\partial_y\theta(x,y)=2\sinh w(x,y)\left(\sin\theta(x,y)+\frac{F^{\prime}(x)}{2F(x)}\right).
		\end{equation}
		An integration yields
		\begin{equation*}
			\tan\frac{\theta(x,y)}{2}=\frac{1}{b(x)}\left( \sqrt{b^2(x)-1}\tan({J_1+J_2}) -1 \right),
		\end{equation*}
		where
		\begin{align*}
			J_1=J_1(x)&=\arctan \left( \frac{1+b(x)}{\sqrt{b^2(x)-1}}  \right),
		\end{align*}
		and
		\begin{align*}
			J_2&=J_2(x,y)=\sqrt{b^2(x)-1}\;Y(x,y).
		\end{align*}
	\end{proof}
	Note that the only, but crucial use of the condition (\ref{omega1}) is that $\frac{\partial_{x}w(x,y)}{2\sinh{w(x,y)}}=\frac{F'(x)}{2F(x)}=b(x)$ is a function of $x$, which facilitates the integration in (\ref{thetay}).

	It remains  to find explicit formulas of $X$ and $Y$ in terms of $x$ and $y$. Recall that by (\ref{Xvar}) and (\ref{Yvar}), these expressions are given by
	\[
	X=\int_{0}^{x}\frac{1+F^{2}(t)G^2(0)}{1-F^{2}(t)G^2(0)}dt,
	\]
	and
	\[
	Y=2F(x)\int_{0}^{y} \frac{G(s)}{1-F^2(x) G^2(s)}ds,
	\]
	where $F$ and $G$ are determined by (\ref{odef}) and (\ref{odeg}), respectively. Therefore, $F$ and $G$ are elliptic functions in general, making the integrals in the computation of $X$ and $Y$ nontrivial. In order to compute these integrals,  we make use of (\ref{odef}) and (\ref{odeg}) respectively.

	Write
	\begin{align*}
		X&=\int_{0}^{x}\frac{1+F^{2}(t)G^2(0)}{1-F^{2}(t)G^2(0)}dt\\
		&=\int_{0}^{x}\frac{\left(1+F^{2}(t)G^2(0)\right) F^{\prime}(t)}{\left(1-F^{2}(t)G^2(0)\right)\sqrt{A F^{4}(t)+B F^2(t)+C}}dt.
	\end{align*}
	Hence,
	\begin{equation*}
		X=\int_{F(0)}^{F(x)}\frac{1+v^{2}G^2(0)}{\left(1-v^{2}G^2(0)\right)\sqrt{A v^{4}+B v^2+C}}dv.
	\end{equation*}
	This integral can be computed explicitly and it involves first and third kind elliptic functions.

	Similarly, we write
	\[
	Y=F(x)\int_{0}^{y} \frac{2G(s) G^{\prime}(s)}{\left(1-F^2(x) G^2(s)\right)\sqrt{-CG^4(s)-(B-4)G^2(s)-A }}ds.
	\]
	Hence,
	\begin{equation*}
		Y=F(x)\int_{G^2(0)}^{G^2(y)}\frac{du}{\left(1-F^2(x) u\right)\sqrt{-Cu^2-(B-4)u-A }}.
	\end{equation*}
	This integral can be computed explicitly and it is an elementary function of $G^2(y)$.

	\subsection{The corresponding harmonic maps.}
	Given Lemma \ref{BackLemma}  we  now complete the proof of our main result.

	\begin{proof}[End of the proof of Theorem \ref{new class}]
		Note that from Lemma \ref{BackLemma} we know the formula of $\theta$. Thus, it is left to express the harmonic map in terms of $\theta$ and $b=\frac{\partial_{x}w(x,y)}{2\sinh{w(x,y)}}=b(x)$. The proof relies on the formulas given by Proposition \ref{Theorem2}. The trick of the proof is to use (\ref{thetay}) in order to substitute\begin{equation*}
			\sinh{w(x,s)}=\frac{\partial_{s}\theta(x,s)}{2 (\sin{\theta(x,s)}+b(x))},
		\end{equation*}
		in the integrals $I_2$ and $I_4$. For $I_1$ and $I_3$, we use the initial condition $\theta(x,0)=\frac{\pi}{2}$. A direct calculation now yields
		\begin{align*}I_1&=X, \\
			I_2&=\frac{1}{2}\log{\frac{\sin{\theta(x,y)}+b}{1+b}},\\
			I_3&=0,\\
			I_4&=-\frac{1}{2}\frac{e^{2X}\cos{\theta(x,y)}}{1+b}.
		\end{align*}
		Then the result follows by Prosopition \ref{Theorem2}  and a straightforward computation of $R$ and $S$.
	\end{proof}

	We elaborate on this result by an example motivated by Proposition \ref{Kenmotsu}. We consider the case $\alpha=1, \beta=2$, $\tanh{\frac{w(0,0)}{2}=\frac{\sqrt{2}}{2}}$ in Proposition \ref{Kenmotsu}, which implies that both $f, g$ are elementary functions. Next, using Theorem \ref{new class} we provide an entirely new harmonic map to the hyperbolic plane.

	\begin{example}
		Consider $S(0,0)=1$, $R(0,0)=0$ and $\tanh{\frac{w(x,y)}{2}}=\frac{\sqrt{2}}{2}\frac{\cosh{(2\sqrt{2}y)}}{\cos{(2x)}}$.
		We find
		\[
		X=x+\arctanh{(\tan{(2x)})},\quad  b=b(x)=\tan{(2x)},
		\]
		\[
		Y= {\frac{1}{\sqrt{1-\tan^2(2x)}} }  \arctanh{\frac{\sinh{(2\sqrt{2}y)}}{\sqrt{1-2\sin^2{(2x)}}}},
		\]
		and
		\[
		\tan{\frac{\theta(x,y)}{2}}=\frac{\sinh{(2\sqrt{2}y)}+\cos{(2x)}-\sin{(2x)}}{\cos{(2x)}-\sin{(2x)}-\sinh{(2\sqrt{2}y)}}.
		\]
		After a lengthy computation we obtain
		\[{
			S(x,y)=e^{2x}\frac{1+2 \cos{(4x)}-\cosh{(4\sqrt{2}y)}}{1+\cosh{(4\sqrt{2}y})-2\sin{(4x})}}
		\]
		and
		\[{
			R(x,y)=-4e^{2x}\frac{\cos{(2x)}\sinh{(2\sqrt{2}y)}}{1+\cosh{(4\sqrt{2}y})-2\sin{(4x})}.}
		\]
		Therefore, $u=R+iS$ is a harmonic map that corresponds to $w$.
	\end{example}

	\section{One-soliton solutions of the Sine-Gordon equation}\label{app}
	In the previous section, we started with a special solution $w$ to the sinh-Gordon equation and then determined $\theta$ by the B{\"a}cklund transformation, in order to construct harmonic maps by the algorithm proposed in Proposition \ref{Theorem2}. We now follow the reverse direction, and start with a special solution $\theta$ to the sine-Gordon equation.

	More precisely, we are motivated by an example in \cite[p.18]{FotDask}, where a special solution $\theta=\theta(x)=\arcsin \tanh (2x)$ of the sine-Gordon equation is considered. Notice that in the aforementioned example, $\theta(x)$ is an elementary function. However, we prove in this section that in general $\theta(x)$ is an elliptic function. In addition, we prove that if $\theta=\theta(x)$ is a solution to the sine-Gordon equation (\ref{sin-Gordon eqn}), then the function $w$ determined by the B{\"a}cklund transformation is in the family of solutions that we have studied in Section \ref{sec:The sinh-Gordon equation},  i.e. $w=\arctanh (F(x)G(y))$. Next, we provide an  explicit formula for the corresponding harmonic map.

	For clarity reasons, we first present the definitions and properties of the elliptic Jacobi functions. The formulation used in this paper is taken from \cite{MilThom64}.

	The elliptic integral of the first kind $F(\phi \vert n)$  and the Jacobi elliptic function $\sn(v\vert n)$  are defined by the formula
	\begin{equation*}\label{eq:Jacobi_ampli}
		F(\phi\vert n) = v=
		\int_0^x \, \dfrac{dt}{   \sqrt{ \left( 1-t^2\right)\left( 1 -n t^2\right)}}=
		\sn^{-1}(x\vert n),
	\end{equation*}
	where
	\[
	\sn(v \vert n)=x=\sin \phi.
	\]
	{Define}
	\[\cn(v \vert n)=\cos \phi ,
	\quad \dn(v \vert n)= \sqrt{1-n  \sin^2\, \phi}.
	\]

	We first prove the following result, which provides the characterization for $w$.

	\begin{proposition}
		If $\theta = \theta(x)$ is a solution of the sine-Gordon equation
		\begin{equation*}
			\Delta\theta = - 2\sin(2\theta),
		\end{equation*}
		the associated solution $w$ of the sinh-Gordon equation is
		\begin{equation*}
			\tanh \frac{w(x,y)}{2} = F(x)G(y) =
			\frac{\sqrt{ab}\tan(\frac{\sqrt{ab}}{2}y + k)}{\theta'(x) - 2\cos\theta(x)},
		\end{equation*}
		where $a = 2\cos\theta(0) + \theta'(0)$, $b = 2\cos\theta(0) - \theta'(0)$, and $k\in \mathbb{C}$ is such that
		$\tanh \frac{w(0,0)}{2}=- {\sqrt{\frac{a}{b}}\tan(k)}$.
	\end{proposition}

	\noindent	\textbf{Remark.}	The above formula holds true also for the case $a b <0.$

	\begin{proof}
		Our first task is to determine the function $\theta$. Consider $\theta= \theta(x)$ to be a solution of the sine-Gordon equation
		$$
		\Delta\theta = -2\sin(2\theta).
		$$
		Then,
		$$
		\theta'' = -2\sin(2\theta)
		$$
		and therefore
		$$
		\theta'(x)^{2} = - 4\sin^{2}\theta(x) + c^{2},
		$$
		where $c^2 = \theta'(0)^{2} + 4\sin^2\theta(0)$.
		A simple calculation reveals that
		$$
		\int_{0}^{\theta(x)}\frac{d\psi}{\sqrt{1 - \frac{4}{c^{2}}\sin^2(\psi)}} = c_{1} + c x,
		$$
		which implies
		\begin{align}\label{sn-theta}
			\sin\theta(x) = \sn(c x + c_{1}| \frac{4}{c^{2}}), \; \text{ where } c_{1} = \int_{0}^{\theta (0)}\frac{d\psi}{\sqrt{1 - \frac{4}{c^{2}}\sin^2(\psi)}}.
		\end{align}
		Moreover, we have
		\begin{align}\label{cndn-theta}
			\cos\theta(x) = \cn(c x + c_{1}| \frac{4}{c^{2}}),\quad \theta'(x)=c\; \dn(c x + c_{1}| \frac{4}{c^{2}}),
		\end{align}
		from which we derive the following differential equation
		\begin{equation}\label{thetaODE}
			\theta'^2-4\cos^2 \theta=c^2-4.
		\end{equation}
		Note that this implies
		\begin{equation*}
			ab=4-c^2.
		\end{equation*}

		We next prove that the function $w$ is of the form $w=2\arctanh (F(x)G(y))$. For that, we turn to the B{\"a}cklund transformation, which rewrites as
		\begin{align}
			\partial_x w&=-2\sinh w \sin \theta, \label{Back1Sec5} \\
			\partial_y w+\theta'(x)&=-2\cosh w \cos \theta \label{Back2Sec5}.
		\end{align}
		Integrating \eqref{Back1Sec5}, we get
		$$
		\int_{w(0,y)}^{w(x,y)} \frac{d\phi}{\sinh(\phi)} + 2\int_{0}^{x}\sin\theta(s)ds = 0.
		$$
		It follows that
		$$
		\tanh(\frac{w(x,y)}{2}) = \tanh(\frac{w(0,y)}{2})e^{-2I(x)},
		$$
		where $I(x) = \int_{0}^{x}\sin\theta(s)ds$.

		We now apply Theorem \ref{w thm} to find an explicit formula for $G(y)$.
		Set
		\begin{equation*} F(x) = \exp(-2I(x)) \; \text{and} \;G(y) = \tanh(\frac{w(0,y)}{2}).
		\end{equation*}
		By Theorem \ref{w thm}, the function $F$ satisfies
		\begin{equation}\label{FODE}
			\frac{(F'(x))^{2}}{F^2(x)} = AF^2(x) + B+ \frac{C}{F^{2}(x)},
		\end{equation}
		for some constants $A, B, C$. We next determine these constants.

		Observe first that $I(x)$ can be computed explicitly:
		\begin{align}
			I(x)&=\int_0^x\sin \theta(s)ds=\int_0^x \sn(cs+c_1|\frac{4}{c^2})ds \label{I(x)def} \\
			&=\frac{1}{2}\log\frac{c\;\dn(c_1|\frac{4}{c^2})+2\cn(c_1|\frac{4}{c^2})}{c\;\dn(cs+c_1|\frac{4}{c^2})+2\cn(cs+c_1|\frac{4}{c^2})}. \notag
		\end{align}
		This implies that
		\begin{align}
			F(x)&=e^{-2I(x)}=\frac{c\;\dn(cs+c_1|\frac{4}{c^2})+2\cn(cs+c_1|\frac{4}{c^2})}{a} \label{Felliptic1},\\
			\frac{1}{F(x)}&=e^{2I(x)}=-\frac{c\;\dn(cs+c_1|\frac{4}{c^2})-2\cn(cs+c_1|\frac{4}{c^2})}{b}\label{Felliptic2}.
		\end{align}
		Differentiating $F(x)=e^{-2I(x)}$ and using \eqref{sn-theta}, we obtain $F'(x)/F(x)=-2\sn(cs+c_1|\frac{4}{c^2})$. Plugging this quotient into the left hand side of \eqref{FODE} and using the elliptic expressions (\ref{Felliptic1}) and (\ref{Felliptic2}) for the right hand side, we obtain
		$A = - \frac{a^{2}}{4}$, $B = \frac{c^{2} + 4}{2}=\frac{8-ab}{2}$ and
		$C = - \frac{b^{2}}{4}$.

		Finally, by Theorem \ref{w thm}, the function $G(y)$ satisfies the equation
		\begin{equation}\label{GODE}
			(G'(y))^{2} = \frac{b^{2}}{4}G^{4}(y) + \frac{ab}{2}G^{2}(y) +
			\frac{a^{2}}{4}= \frac{1}{4}(bG^{2}(y) + a )^{2}.
		\end{equation}
		Choose
		\begin{equation}\label{Gwithk}
			G(y) = -\sqrt{\frac{a}{b}}
			\tan(\frac{\sqrt{ab}}{2}y +k),
		\end{equation}
		where $\tanh \frac{w(0,0)}{2}=- {\sqrt{\frac{a}{b}}\tan(k)}$. Then, one can verify that $w(x,y)$ satisfies the second equation (\ref{Back2Sec5}) of the B{\"a}cklund transformation, so we have found the  associated $w(x,y)$ to $\theta(x)$.
	\end{proof}

	Having calculated $\theta(x)$ and $w(x,y)$, we are now ready to compute the corresponding harmonic map by calculating the integrals $I_{1}$, $I_{2}$, $I_{3}$ and $I_{4}$ in Proposition \ref{Theorem2}.

	\begin{proof}[Proof of Theorem \ref{thetaxthm}]
		First, observe that (\ref{cndn-theta}) with (\ref{Felliptic1}) and (\ref{Felliptic2}) imply that
		\begin{align}
			F(x)&=\frac{\theta'(x)+2\cos\theta(x)}{a}, \label{Ftheta}\\
			\frac{1}{F(x)}&=-\frac{\theta'(x)-2\cos\theta(x)}{b}.\label{Finvtheta}
		\end{align}

		Let us start from $I_1$. Since $I'(x)=\left( \int_{0}^{x}\sin\theta(s)ds \right)'=\sin \theta(x)$ and $\tanh\frac{w(x,y)}{2}=e^{-2I(x)}G(y)$, applying an elementary trigonometric equality we obtain
		\begin{align*}
			I_{1}(x) &= \int_{0}^{x}\cosh w(t,0)\sin\theta(t)dt =
			\int_{0}^{x} \frac{\exp(4I(t)) + G^{2}(0)}{\exp(4I(t)) - G^{2}(0)}dI(t)\\
			&=\frac{1}{2} \int_{1}^{\exp(2I(x))} \left(\frac{2u}{u^{2} - G^{2}(0)} - \frac{1}{u}\right)du,
		\end{align*}
		or
		\begin{equation}\label{exp(2I1)}
			\exp(2I_{1}(x)) = \frac{1 - G^{2}(0)F^{2}(x)}{F(x)(1 - G^{2}(0))}=\frac{\frac{1}{F(x)} - G^{2}(0)F(x)}{1 - G^{2}(0)}.
		\end{equation}
		Using (\ref{Finvtheta}), (\ref{Ftheta}) and (\ref{Gwithk}), we conclude
		\begin{equation}\label{exp2I1theta}
			\exp(2I_{1}(x)) =
			\frac{2\cos\theta(x)\cos(2k) - \theta'(x)}{2\cos\theta(0)\cos(2k) - \theta'(0)}.
		\end{equation}

		Next, we compute $I_2$. Once again, using that $\tanh\frac{w(x,y)}{2}=F(x)G(y)$ as well as an elementary trigonometric equality, we obtain
		\begin{align*}
			I_{2}(x,y) &=\cos\theta(x) \int_{0}^{y}\sinh(w(x,s))ds\\
			&=2\cos\theta(x)F(x) \int_{0}^{y} \frac{
				G(s)}{1 - F^{2}(x)G^{2}(s)}ds.
		\end{align*}
		An application of (\ref{GODE}) and a change of variables, yield
		\begin{align*}
			I_{2}(x,y) &= 4\cos\theta(x)F(x) \int_{0}^{y} \frac{
				G(s)}{1 - F^{2}(x)G^{2}(s)} \; \frac{G'(s)}{bG^{2}(s) + a}ds\\
			&=\frac{ 2\cos\theta(x)F(x)}{aF^{2}(x) + b}
			\log(\frac{1 - F^{2}(x)G^{2}(0)}{1 - F^{2}(x)G^{2}(y)}
			\;\frac{bG^{2}(y) + a}{bG^{2}(0) + a}
			)\\
			&=\frac{1}{2}\log(\frac
			{1 - F^{2}(x)G^{2}(y)}
			{1 - F^{2}(x)G^{2}(0)}
			\frac{bG^{2}(0) + a}
			{bG^{2}(y) + a}
			),
		\end{align*}
		where in the last step we used \eqref{Ftheta}, \eqref{thetaODE} and the fact that $ab=4-c^2$. A manipulation similar to the one used to derive \eqref{exp2I1theta} therefore gives
		\begin{equation}\label{expI2}
			\exp(2I_{2}(x,y)) = \frac
			{2\cos\theta(x)\cos(\sqrt{ab}y + 2k) - \theta'(x)}
			{2\cos\theta(x)\cos(2k) - \theta'(x)}.
		\end{equation}
		Then, by Proposition \ref{Theorem2}, we have
		\begin{align*}
			S(x,y) &= S(0,0)\exp(2I_{1}(x))\exp(2I_{2}(x,y))\\
			&=S(0,0)\frac{2\cos\theta(x)\cos(\sqrt{ab}y + 2k) - \theta'(x)}
			{2\cos\theta(0)\cos(2k) - \theta'(x)}.
		\end{align*}

		We proceed to  $I_3$. Due to (\ref{I(x)def}), it can be written as
		\begin{align*}
			I_{3}(x) &= \int_{0}^{x} \exp(2I_{1}(t))\cosh(w(t,0))\cos\theta(t)dt\\
			&=\int_{0}^{x} \exp(2I_{1}(t))\cosh(w(t,0))\cot\theta(t)dI(t).
		\end{align*}

		On the one hand, since $F(x) = e^{-2I(x)}$ satisfies (\ref{FODE}) with $A = - \frac{a^{2}}{4}$, $B = \frac{8-ab}{2}$ and
		$C = - \frac{b^{2}}{4}$,
		we get
		$$
		\sin\theta(x) = I'(x) = -\frac{1}{2}\frac{F'(x)}{F(x)}=\frac
		{\sqrt
			{2(8-ab)e^{4I(x)} - a^{2} - b^{2}e^{8I(x)}}
		}
		{4e^{2I(x)}},
		$$
		therefore
		$$
		\cos\theta(x) = \sqrt{1 - I'(x)^{2}} = \frac{be^{4I(x)} + a}{4e^{2I(x)}}.
		$$
		On the other hand, replacing $e^{2I_1(x)}$ by (\ref{exp(2I1)}), we can write \begin{align*}
			I_{3}(x) &= \int_{0}^{x} \frac{\exp(4I(t)) + G^{2}(0)}{(1 - G^{2}(0))\exp(2I(t))}\frac{be^{4I(x)} + a}
			{\sqrt
				{2(8-ab)e^{4I(x)} - a^{2} - b^{2}e^{8I(x)}}
			}dI(t)\\
			&=\frac{1}{2(1 - G^{2}(0))}
			\int_{1}^{\frac{2\cos\theta(x) - \theta'(x)}{b}}
			\frac{u^{2} + G^{2}(0)}{u^{2}}
			\frac{bu^{2} + a}
			{\sqrt{2(8-ab)u^{2} - a^{2} - b^{2}u^{4}}}
			du,
		\end{align*}
		where we used the change of variables $u=e^{2I(x)}=1/F(x)$ and (\ref{Finvtheta}).
		Thus, taking into account (\ref{Gwithk}), we find
		\[ I_3(x)=\frac{1}{2}\cosh^2\frac{w(0,0)}{2}\;J\left(\frac{2\cos\theta(x) - \theta'(x)}{b}\right),
		\]
		where $J(t)=\int_{1}^{t} \frac{u^{2} + G^{2}(0)}{u^{2}}
		\frac{bu^{2} + a}
		{\sqrt{2(8-ab)u^{2} - a^{2} - b^{2}u^{4}}}
		du.$ Note that $J$ can be computed by the use of elliptic integrals.

		Finally, for the last integral $I_4$ we have
		\begin{align*}
			I_{4}(x,y) &= e^{2I_{1}(x)}\sin\theta(x)\int_{0}^{y}\exp(2I_{2}(x,s))\sinh(w(x,s))ds\\
			&=\frac{1}{2}e^{2I_{1}(x)}\tan\theta(x)\int_{0}^{y}
			e^{(2\cos\theta(x)\int_{0}^{s}\sinh(w(x,\psi))d\psi)}2\cos\theta(x)\sinh(w(x,s))ds\\
			&= \frac{1}{2}e^{2I_{1}(x)}\tan\theta(x)(e^{2I_{2}(x,y)} - e^{2I_{2}(x,0)}) \\
			&=  \frac{1}{2}e^{2I_{1}(x)}\tan\theta(x)(e^{2I_{2}(x,y)} - 1)\\
			&=\frac{\sin\theta(x)(\cos(\sqrt{ab}y + 2k) - \cos(2k))}
			{2\cos\theta(0)\cos(2k) - \theta'(0)},
		\end{align*}
		where in the last step we used (\ref{exp2I1theta}) and (\ref{expI2}). The calculation of $R(x,y) = R(0,0) + 2S(0,0)(I_{3}(x) - I_{4}(x,y))$ is now complete.
	\end{proof}
	Next, using Theorem \ref{thetaxthm}, we provide an entirely new harmonic map to the hyperbolic plane.

	\begin{example}
		Consider the case when $\theta(0) = 0$, $\theta'(0) = 1$, $G(0) = 0=\tanh\frac{w(0,0)}{2}$, $S(0,0) = 1$ and $R(0,0) = 0$. Then, we compute
		\[\sin\theta(x) = \sn(x|4) = \frac{1}{2} \sn (2x|\frac{1}{4}),\]
		\[\tanh(\frac{w(x,y)}{2}) = \frac{\sqrt{3}\tan(\frac{\sqrt{3}}{2}y)}{\theta'(x) - 2\cos\theta(x)},\]
		and
		$$
		I(x) = \frac{1}{2}\log(2\cos\theta(x) - \theta'(x)).
		$$
		Therefore,
		\begin{align*}
			I_{1}(x) &= I(x),\\
			I_{2}(x, y) &= \frac{1}{2}\log(\frac{2\cos\theta(x)\cos(\sqrt{3}y) - \theta'(x)}
			{2\cos\theta(x) - \theta'(x)}),\\
			I_{3}(x) &= \frac{1}{2}\int_{1}^{2\cos\theta(x) - \theta'(x)} \frac{u^{2} + 3}{\sqrt{10u^{2} - u^{4} - 9}}du,\\
			I_{4}(x,y) &= \sin\theta(x)(\cos(\sqrt{3}y) - 1).
		\end{align*}
		Then, the corresponding harmonic map then is $u(x,y) = R(x,y) + iS(x,y)$, where
		$$
		R(x,y) = J(2\cos\theta(x)-\theta'(x))
		- 2\sin\theta(x)(\cos(\sqrt{3}y) - 1)
		$$
		and
		$$
		S(x,y) = 2\cos\theta(x)\cos(\sqrt{3}y) - \theta'(x).
		$$
	\end{example}


\section{Introduction}\label{sec: introduction}
	\subsection{Stopping games}
	Dynkin \cite{Dynkin1967} considered a game with two players who observe a stochastic process $(R_t,X_t)_{t=0}^\infty$. Player $1$ (resp. 2) picks a stopping time $\tau$ (resp. $\nu$) such that $P(X_{\tau}\leq0,X_{\nu}>0)=1$. The purpose of Player 1 (resp. 2) is to maximize (resp. minimize) the expected value of $R_{\min\{\tau,\nu\}}$. Dynkin proved that once $\sup_t R_t$ is integrable, the game has a value. Kifer \cite{Kiffer1971} and Neveu \cite{Neveu1975} considered a similar game in which two players observe a three-dimensional stochastic process $(X_t,Y_t,Z_t)_{t=0}^\infty$. Each player chooses a stopping time and the purpose of Player 1 (resp. 2) is to maximize (resp. minimize) the expected value of
	\begin{equation}
	X_{\tau}\textbf{1}_{\{\tau<\nu\}}+Y_{\nu}\textbf{1}_{\{\tau>\nu\}}+Z_{\tau}\textbf{1}_{\{\tau=\nu\}}\,.
	\end{equation}
	In particular, it is possible that the players choose stopping times $\tau$ and $\nu$ such that $P(\tau=\nu)>0$. In two-player nonzero-sum stopping games, the processes $(X_t)$, $(Y_t)$ and $(Z_t)$ are two dimensional; the first (resp. second) coordinate is interpreted as the payoff of Player 1 (resp. Player 2). In general, the existence of an equilibrium point in these games is ensured when the joint distribution of the payoff processes satisfies some conditions. For example, see Morimoto \cite{Morimoto1986}, Mamer \cite{Mamer1987}, Ohtsubo \cite{Ohtsubo1987} and Hamad\'ene \textit{et al} \cite{Hamadene2010}. Also, some relevant literature reviews are given by  Sakaguchi \cite{Sakaguchi1995}, Novak \textit{et al.} \cite{Novak1999} and Ja\'skiewicz \textit{et al.} \cite{Jaskiewicz2018}.

	Another strand of literature studies stopping games with random stopping times. See, e.g., Yasuda \cite{Yasuda1985}, Rosenberg \textit{et al.} \cite{Rosenberg 2001}, Shmaya \textit{et al.} \cite{Shmaya2004} for the discrete-time case and Laraki \textit{et al.} \cite{Laraki2013} for the continuous time case. Notably, there is an extensive literature about continuous time versions of stopping games. For example, consider the classical works of Bismut \cite{Bismut1977} and Lepeltier \textit{et al.} \cite{Lepeltier 1984}. More recent works in this direction are those of: Attard \cite{Attard2017,Attard2018}, De Angelis \textit{et al.} \cite{De Angelis2018a,De Angelis2018b}, Hamad\'ene \textit{et al.} \cite{Hamadene2010}. In addition, a continuous (resp. discrete) version of stopping games with more than two players was discussed by Hamad\'ene \textit{et al.} \cite{Hamadene2014b} (resp. Hamad\'ene \textit{et al.} \cite{Hamadene2014a}).

	Stopping games are motivated by various applications. Traditionally, a task that has stimulated mathematicians to analyse models with unique implications on optimal-stopping theory is the celebrated secretary's problem (for surveys, see e.g., Ferguson \cite{Ferguson1989} and Freeman \cite{Freeman1983}). Some  versions of this problem with more than one decision maker were analyzed by Bruss \textit{et al.} \cite{Bruss1998}, Chen \textit{et al.} \cite{Chen1997}, Enns \textit{et al.} \cite{Enns1985}, Fushimi \cite{Fushimi1981},  and Sakaguchi \cite{Sakaguchi1980}. Some other applications of stopping games are:  Shrinking markets (see, e.g., Ghemawat \textit{et al.} \cite{Ghemawat}), duels (see, e.g., Lang \cite{Lang1975}) and financial markets with a recallable option (see, e.g., Hamad\'ene \textit{et al.} \cite{Hamadene2006,Hamadene2010} and Kifer \cite{Kifer2000}).

	\subsection{Asymmetric information}
	Recently, there is a growing attention towards stopping games with asymmetric information. Lempa \textit{et al.} \cite{Lempa2013} considered a stopping game with a finite-horizon which is an exponentially distributed random variable, where only one player is exposed to the value of this random variable.  Esmaeeli \textit{et al.}    \cite{Esmaeeli2019} and Gr\"un \cite{Grun2013} presented two models in which the asymmetry of information is modelled as in the
	classical work of Aumann and Maschler \cite{Aumann1995}. Namely, there is a finite set of states of nature in which the game can take place. The true state is a random variable whose distribution is a common knowledge. The asymmetry in information is due to the assumption that the true state is exposed only to  Player 1. An extension of this framework was discussed by Renault \cite{Renault2006}. De Angelis \textit{et al.} \cite{De Angelis2021} considered a game with payoffs which are diffusion processes, where only one player knows the exact value of the drift and the diffusion coefficient. Mamer \cite{Mamer1987} and Ohtsubo \cite{Ohtsubo1991} considered a class of games with a monotone payoff structure where  players have information flows that are modelled as two possibly different filtrations. Gapaev \textit{et al.} \cite{Gapeev2021} considered a game with asymmetric information in the context of a financial market. In particular Gapaev \textit{et al.} \cite{Gapeev2021} refer to a solution concept which is called Stackelberg equilibrium (for an analysis of this solution concept in other stopping games see, e.g., the works of Ektr\"om \textit{et al.} \cite{Ekstrom2008} and Peshkir \cite{Peshkir2009}). Gensbittel \cite{Gensbittel2019} considered a two-player zero-sum stopping game with payoff processes which are determined by a Markov chain. The evolution of this chain is observed only by one of the Players.  Ashkenazi-Golan \textit{et al.} \cite{Ashkenazi-Golan2020} analyzed a two-player two-state Markov game such that only one of the players knows when the state changes.

	One major challenge in the analysis of stopping games with asymmetric information is due to the fact that they cannot be analysed by the classical backward-induction approach. To see why, notice that when an informed player decides not to stop, he reveals information to the uninformed players. Indeed, he reveals that his payoff by stopping in early stages is not too high. In particular, the decision of the uninformed player at time $t$ may not be determined by the play after time $t$, but also on the play before time $t$ (see also Skarupski \textit{et al.} \cite{Skarupski2018}). Evidently, as it was shown by Laraki \cite{Laraki2000}, a zero-sum stopping games with asymmetric information may not have a value even under the standard assumptions. However, there are also other cases in which an equilibrium result may be derived. For example, Laraki \textit{et al.} (see Section 2.8.4.3 of \cite{Laraki2015}) mentioned that for finite-horizon zero-sum stopping games with incomplete information on one side, the value exists in discrete-time using Sion's minmax theorem.

	\subsection{The current work}
	We study finite-horizon two-player nonzero-sum stopping games in discrete-time, where $Z\equiv X$: If the players stop simultaneously, the payoff is the same as if Player 1 stops alone. Similar assumptions appear in the works of Hamad\'ene \textit{et al.} \cite{Hamadene2010} and Ramsey \cite{Ramsey2007}. Formally,  consider an integer $0<T<\infty$  with some two discrete-time processes $(X_t^1,Y_t^1)_{t=0}^T$ and $(X_t^2,Y_t^2)_{t=0}^T$. We assume that Player 1 observes all four processes while Player 2 observes only $(X_t^2,Y_t^2)_{t=0}^T$. In addition, Player 1 (resp. 2) picks a stopping time $\tau$ (resp. $\nu)$ with respect to his information flow and his goal is to maximize the expectation of
	\begin{equation}
	X^1_\tau\textbf{1}_{\left\{\tau\leq\nu\right\}}+Y^1_\nu\textbf{1}_{\left\{\tau>\nu\right\}}\ \ \left(\text{resp. }X^2_\nu\textbf{1}_{\left\{\nu<\tau\right\}}+Y^2_\tau\textbf{1}_{\left\{\tau\leq\nu\right\}}\right).
	\end{equation}
	One main difficulty in the analysis of this game is the fact that Player 1 may pick a stopping time $\tau$ which is not a stopping time with respect to the information flow of Player 2. Consequently, it is possible that Player 2 should decide whether to stop or keep playing without knowing whether Player 1 has already stopped. It is reasonable to expect that if Player 2 knew that Player 1 has not stopped yet, he could learn something new which might help him to gain greater utility. With some extent of surprise, we will show  that this intuition is wrong. Namely, Player 2 may always assume that Player 1 has not stopped yet and by doing so he will get the maximal expected utility that could be obtained under the assumption that he is informed about the stopping epoch of Player 1.  Then, this result is applied in order to prove that once $X_t^1\geq Y_t^1$ for every $0\leq t\leq T$, a recursive construction in the style of Hamad\'ene \textit{et al.} \cite{Hamadene2010} yields a pure-strategy Nash equilibrium. To see a crucial difference between the current work and Hamad\'ene \textit{et al.} \cite{Hamadene2010}, recall that Hamad\'ene \textit{et al.} \cite{Hamadene2010} studied stopping games with symmetric information. Thus, their construction yields stopping times which are measurable with respect to (w.r.t.) the information of Player 1 but not with respect to the information of Player 2. Hence, a proper modification of the construction is needed.

	The rest of this paper is organized as follows: Section \ref{sec: description} contains a detailed presentation of the model assumptions with the precise statements of the main results. In addition, this section contains a sketch of the proof of the equilibrium result. This sketch is based on some benchmarks which are given with some intuitive explanations. For clarity, the technical proofs are supplied in Section \ref{sec: proof}.

	\section{Game description and the main results}\label{sec: description}
	Let $T\in(0,\infty)$ be a positive integer and consider a probability space $(\Omega,\mathcal{F},P)$ with two filtrations $\mathbb{F}\equiv\left(\mathcal{F}_t\right)_{t=0}^T$ and $\mathbb{G}=\left(\mathcal{G}_t\right)_{t=0}^T$ such that $\mathcal{G}_t\subseteq\mathcal{F}_t$ for every $0\leq t\leq T$. Denote the set of all $\mathbb{F}$ (resp. $\mathbb{G}$)-stopping times by $\mathcal{T}(\mathbb{F})$ $\left(\text{resp. }\mathcal{T}(\mathbb{G})\right)$. Let $\mathcal{D}$ be the space of all $\mathbb{F}$-adapted $\mathbb{R}$-valued  processes $\left(S_t\right)_{t=0}^T$ such that

	\begin{equation}\label{eq: integrability}
	E\left(\max_{0\leq t\leq T}|S_t|\right)<\infty\,.
	\end{equation}
	Let $X^1,X^2,Y^1,Y^2\in\mathcal{D}$ such that $X^2$ and $Y^2$ are adapted to $\mathbb{G}$. Consider a two-player stopping game $\Gamma$ in which:
	\begin{itemize}
		\item Player 1 picks $\tau\in\mathcal{T}(\mathbb{F})$ and Player 2 picks $\nu\in\mathcal{T}(\mathbb{G})$.

		\item For every $\tau\in\mathcal{T}(\mathbb{F})$ and $\nu\in\mathcal{T}(\mathbb{G})$ the expected utility of Player 1 is given by
		\begin{equation}\label{eq: J_1}
		J_1(\tau,\nu)\equiv E\left[X^1_\tau\textbf{1}_{\left\{\tau\leq\nu\right\}}+Y^1_\nu\textbf{1}_{\left\{\tau>\nu\right\}}\right]\,,
		\end{equation}
		and the expected utility of Player 2 is given by
		\begin{equation}
		J_2(\tau,\nu)\equiv E\left[X^2_\nu\textbf{1}_{\left\{\nu<\tau\right\}}+Y^2_\tau\textbf{1}_{\left\{\tau\leq\nu\right\}}\right]\,.
		\end{equation}
	\end{itemize}

	\begin{remark}
		\normalfont Note that the assumption that $X^2$ and $Y^2$ are adapted to $\mathbb{G}$ practically means that Player 2 observes the evolution of his payoff processes over time.
	\end{remark}

	\begin{remark}
		\normalfont By the definitions of $J_1$ and $J_2$, Player 1 has a priority in stopping: The payoff if the two players stop simultaneously is the same as Player 1 stops alone.
	\end{remark}

	\begin{definition}
		A pair $(\tau^*,\nu^*)\in\mathcal{T}(\mathbb{F})\times\mathcal{T}(\mathbb{G})$ constitutes a pure-strategy Nash equilibrium in $\Gamma$ if and only if (iff)
		\begin{equation}
		J_1(\tau,\nu^*)\leq J_1(\tau^*,\nu^*), \ \ \forall\tau\in\mathcal{T}(\mathbb{F})\,,
		\end{equation}
		and
		\begin{equation}
		J_2(\tau^*,\nu)\leq J_2(\tau^*,\nu^*), \ \ \forall\nu\in\mathcal{T}(\mathbb{G})\,.
		\end{equation}
	\end{definition}
	In the current work, we are going to prove that if
	\begin{equation}\label{eq: sufficient condition}
	X_t^1\geq Y_t^1\ \ , \ \ \forall0\leq t\leq T\ , \ P\text{-a.s.},
	\end{equation}
	then $\Gamma$ admits a pure-strategy Nash equilibrium which follows from a construction in the style of Hamad\'ene \textit{et al.} \cite{Hamadene2010}.

	In the rest of this section, there is an elaboration which is necessary in order to make this statement more precise. In addition, we will review some important benchmarks like the analysis of the best response of Player 2 and the recursive construction which yields the equilibrium. Especially, the results which appear in the following Subsection 2.1 and Subsection 2.2 do not require the condition which appears in \eqref{eq: sufficient condition}.
	\subsection{The best response of Player 2}\label{subsec: best-response}
	For every $0\leq t\leq T$, let $\mathcal{T}_t(\mathbb{F})$ (resp. $\mathcal{T}_t(\mathbb{G})$)  be the set of all $\mathbb{F}$ (resp. $\mathbb{G}$)-stopping times $\tau$ such that $P(t\leq\tau\leq T)=1$. Furthermore, for every $\tau\in\mathcal{F}$ define the filtration $\mathbb{G}^\tau\equiv\left(\mathcal{G}_t^\tau\right)_{t=0}^T$ as follows:
	\begin{equation}
	\mathcal{G}_t^\tau\equiv\mathcal{G}_t\vee\sigma\left\{\textbf{1}_{\left\{\tau\leq s\right\}};0\leq s\leq t\right\}, \ \ \forall 0\leq t\leq T\,.
	\end{equation}
	Note that once Player 1 picks $\tau\in\mathcal{T}(\mathbb{F})$,  the only difference between Player 2 and an observer with a filtration $\mathbb{G}^\tau$ is the knowledge of whether Player 1 has already stopped. Thus, the following theorem states that even if Player 2 knows whether Player 1 has already stopped, his expected utility will remain the same. This is because at any moment $t$ such that $P(\tau>t)>0$, Player 2 (whose information flow is described by the filtration $\mathbb{G}$) may conduct his decision-making under the assumption that Player 1 has not stopped yet.
	\begin{theorem}\label{theorem: conditional expectation}
		Let $\tau\in\mathcal{T}(\mathbb{F})$ and $t\in[0,T-1]$.  Then, there exists $\hat{\nu}_t\in\mathcal{T}_t(\mathbb{G})$ such that the next equality holds $P$-a.s.

		\begin{align}\label{eq: ess-sup}
		&E\left[X^2_{\hat{\nu}_t}\textbf{1}_{\left\{\hat{\nu}_t<\tau\right\}}+Y^2_{\tau}\textbf{1}_{\left\{\tau\leq\hat{\nu}_t\right\}}\big|\mathcal{G}_t^\tau\right]=\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G}^\tau)}E\left[X^2_\nu\textbf{1}_{\left\{\nu<\tau\right\}}+Y^2_{\tau}\textbf{1}_{\left\{\tau\leq\nu\right\}}\big|\mathcal{G}_t^\tau\right]\nonumber\\&=Y^2_\tau\textbf{1}_{\{\tau\leq t\}}+\textbf{1}_{\{\tau>t\}}\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G})}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\,.
		\end{align}
	\end{theorem}
	The proof of Theorem \ref{theorem: conditional expectation} is supplied in Subsection \ref{subsec: theorem 1}.
	\subsection{Recursive construction}\label{subsec: recursive construction}
	The recursive construction is as follows: Let $\tau_1=\nu_1=T$ and for every $n\geq1$ define:
	\begin{align}
	&W^{1,n}_t\equiv\\&\text{ess}\sup_{\tau\in\mathcal{T}_t(\mathbb{F})}E\left[X^1_\tau\textbf{1}_{\left\{\tau<\nu_n\right\}}+\left(X^1_T\textbf{1}_{\{\nu_n=T\}}+Y^1_{\nu_n}\textbf{1}_{\{\nu_n<T\}}\right)\textbf{1}_{\left\{\tau\geq\nu_n\right\}}\big|\mathcal{F}_t\right] ,\ \ \forall t\in[0,T]\,,\nonumber\\& \widetilde{\tau}_{n+1}\equiv\inf\left\{t\geq0;W^{1,n}_t=X_t^1\right\}\wedge \nu_n\,,\\&\tau_{n+1}\equiv\begin{dcases}
	\widetilde{\tau}_{n+1}, & \widetilde{\tau}_{n+1}<\nu_n, \\
	\tau_n, & \text{o.w.}
	\end{dcases}
	\end{align}

	\begin{align}
	&W^{2,n}_t\equiv\text{ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G}^{\tau_n})}E\left[X^2_\nu\textbf{1}_{\left\{\nu<\tau_n\right\}}+Y^2_{\tau_n}\textbf{1}_{\left\{\tau_n\leq\nu\right\}}\big|\mathcal{G}_t^{\tau_n}\right]\ \ ,\ \ \forall t\in[0,T]\,,\\& \widetilde{\nu}_{n+1}\equiv\inf\left\{t\geq0;W^{2,n}_t=X^2_t\right\}\wedge \tau_n\,,\\&\nu_{n+1}\equiv\begin{dcases}
	\widetilde{\nu}_{n+1}, & \widetilde{\nu}_{n+1}<\tau_n, \\
	\nu_n, & \text{o.w.}
	\end{dcases}
	\end{align}
	By construction, both of the sequences $(\tau_n)_{n=1}^\infty$ and $(\nu_n)_{n=1}^\infty$ are nonincreasing and hence they have limits $\tau^*$ and $\nu^*$ as $n\to\infty$. When $\mathbb{F}$ and $\mathbb{G}$ are the same, Hamad\'ene \textit{et al.} \cite{Hamadene2010} showed (in an analogue continuous-time model) that $(\tau^*,\nu^*)$ is a pure-strategy Nash equilibrium. In our model, the following Proposition \ref{lemma: Hamadene} is an analogue statement. Its proof is provided in Subsection \ref{subsec: proposition 1}.

	\begin{proposition}\label{lemma: Hamadene}
		There exist limits
		$\tau^*\equiv\lim_{n\to\infty}\tau_n$ and $\nu^*\equiv\lim_{n\to\infty}\nu_n$ such that:
		\begin{enumerate}
			\item $\tau^*\in\mathcal{T}(\mathbb{F})$ and $\nu^*\in\mathcal{T}(\mathbb{F})$.

			\item \begin{equation}
			J_1(\tau,\nu^*)\leq J_1(\tau^*,\nu^*), \ \ \forall\tau\in\mathcal{T}(\mathbb{F})\,.
			\end{equation}

			\item \begin{equation}
			J_2(\tau^*,\nu)\leq J_2(\tau^*,\nu^*), \ \ \forall\nu\in\mathcal{T}(\mathbb{G})\,.
			\end{equation}
		\end{enumerate}
	\end{proposition}
	Importantly, $\nu^*$ is not necessarily a $\mathbb{G}$-stopping time, and hence an equilibrium result will not follow from the above-mentioned construction unless some modifications are applied. Specifically, for each $n\geq1$ and $0\leq t\leq T$ define
	\begin{align}
	&V^n_t\equiv\begin{dcases}
	\text{ess}\sup_{\theta\in\mathcal{T}_t(\mathbb{G})}E_{P|\left\{\tau_n>t\right\}}\left[X^2_\theta\textbf{1}_{\left\{\theta<\tau_n\right\}}+Y^2_{\tau_n}\textbf{1}_{\left\{\tau_n\leq\theta\right\}}|\mathcal{G}_t\right], & P\left(\tau_n>t\right)>0, \\
	\infty, & \text{o.w.},
	\end{dcases}
	\\& \theta_{n+1}\equiv \inf\left\{t\geq0;V^n_t=X^2_t\right\}\wedge T\,.
	\end{align}
	In particular, $\theta^*\equiv\inf_{n\geq2}\theta_n$ is a $\mathbb{G}$-stopping time as an infimum of a sequence of $\mathbb{G}$-stopping times. Moreover, for each $n\geq1$, Theorem \ref{theorem: conditional expectation} implies that
	\begin{align}\label{eq: tilde nu}
	\widetilde{\nu}_{n+1}&=\inf\left\{t\geq0;W^{2,n}_t=X^2_t\textbf{1}_{\left\{t<\tau_n\right\}}+Y^2_{\tau_n}\textbf{1}_{\left\{\tau_n\leq t\right\}}\right\}\\&=\inf\left\{t\geq0;Y^2_{\tau_n}\textbf{1}_{\{\tau_n\leq t\}}+V_t^n\textbf{1}_{\left\{t<\tau_n\right\}}=X^2_t\textbf{1}_{\left\{t<\tau_n\right\}}+Y^2_{\tau_n}\textbf{1}_{\left\{\tau_n\leq t\right\}}\right\}\nonumber\\&=\inf\left\{t\geq0; V_t^n=X_t^2\right\}\wedge\tau_n\nonumber\\&=\theta_{n+1}\wedge\tau_n\ \ , \ \ P\text{-a.s.}\nonumber
	\end{align}
	and hence we shall deduce that
	\begin{equation}\label{eq: nu*}
	\nu^*=\inf_{n\geq1}\left\{\theta_{n+1};\theta_{n+1}<\tau_n\right\}\wedge T, \ \ P\text{-a.s.}
	\end{equation}
	This identity yields the next Corollary \ref{cor: equivalence} which is going to be useful for the proof of the equilibrium result in the next subsection.
	\begin{corollary}\label{cor: equivalence}
		\begin{equation}\label{eq: iff}
		\nu^*<\tau^*\Leftrightarrow \theta^*<\tau^*\ \ , \ \ P\text{-a.s.}
		\end{equation}
	\end{corollary}
	The proof for this corollary is given in Subsection \ref{subsec: corollary 1}.
	\subsection{Equilibrium result}\label{subsec: equilibrium}
	\begin{theorem}\label{thm: equilibrium}
		if the condition \eqref{eq: sufficient condition} is satisfied, then $(\tau^*,\theta^*)$ is a pure-strategy Nash equilibrium.
	\end{theorem}
	The following Lemma \ref{lemma: auxilary} states that it is enough to prove Theorem \ref{thm: equilibrium} under the assumption that $X_t^1>Y_t^1$, $P$-a.s. for every $0\leq t\leq T$. The proof for this Lemma appears in Subsection \ref{subsec: lemma 1}.

	\begin{lemma}\label{lemma: auxilary}
		Assume that the condition \eqref{eq: sufficient condition} is satisfied and let $\epsilon>0$. In addition, define a two-player stopping game $\Gamma_\epsilon$ which is identical to $\Gamma$ except that the expected utility of Player 1 is given by

		\begin{equation}
		J_{1,\epsilon}(\tau,\nu)\equiv E\left[\left(X^1_\tau+\epsilon\right)\textbf{1}_{\left\{\tau\leq\nu\right\}}+Y^1_\nu\textbf{1}_{\left\{\tau>\nu\right\}}\right], \ \ \forall\tau\in\mathcal{T}(\mathbb{F})\ , \ \nu\in\mathcal{T}(\mathbb{G})\,.
		\end{equation}
		If a pair $(\tau_0,\nu_0)\in\mathcal{T}(\mathbb{F})\times\mathcal{T}(\mathbb{G})$ is a pure-strategy Nash equilibrium in $\Gamma_\epsilon$, then it is also a pure-strategy Nash equilibrium in $\Gamma$.
	\end{lemma}
	\subsubsection*{Proof of Theorem \ref{thm: equilibrium}}
	Assume by contradiction that $P(\tau^*>\nu^*)>0$. Since $X_t^1>Y_t^1$ for every $t\in[0,T]$, then $J_1(\tau^*,\nu^*)<J_1(\nu^*\wedge\tau^*,\nu^*)$ and hence the second part of Proposition \ref{lemma: Hamadene} implies that $\nu^*\wedge\tau^*\notin\mathcal{T}(\mathbb{F})$. This is a contradiction with the first part of Proposition \ref{lemma: Hamadene} and hence $\tau^*\leq\nu^*$, $P$-a.s. As a result, Corollary \ref{cor: equivalence} and \eqref{eq: nu*} imply that $\tau^*\leq\theta^*\leq\nu^*$, $P$-a.s. and hence $J_2(\tau^*,\nu^*)=J_2(\tau^*,\theta^*)$. Thus, the third part of Proposition \ref{lemma: Hamadene} yields that  $\theta^*$ is a best response of Player 2 to $\tau^*$.

	Assume by contradiction that there exists $\tau_0\in\mathcal{T}(\mathbb{F})$ such that $J_1(\tau_0,\theta^*)>J_1(\tau^*,\theta^*)$. Denote $\tau'\equiv\tau_0\wedge\theta^*$, and note that $\tau'\in\mathcal{T}(\mathbb{F})$. Since $X_t^1> Y_t^1$ for every $0\leq t\leq T$, deduce that
	\begin{equation}
	J_1(\tau',\theta^*)\geq J_1(\tau_0,\theta^*)>J_1(\tau^*,\theta^*)\,.
	\end{equation}
	Now, $\tau'\leq\theta^*\leq\nu^*$, $P$-a.s., and hence it does not matter whether Player 2 picks either $\theta^*$ or $\nu^*$, when picking $\tau'$, Player 1 is the one who is responsible for stopping, i.e., $J_1(\tau',\theta^*)=J_1(\tau',\nu^*)$. On the other hand, we have already proved that $\tau^*\leq\theta^*\leq\nu^*$, $P$-a.s. and hence $J_1(\tau^*,\theta^*)=J_1(\tau^*,\nu^*)$. This leads to a contradiction  because $\tau^*$ is the best response of Player 1 to $\nu^*$.  $\blacksquare$

\section{Proofs of Theorem \ref{theorem: conditional expectation} and some auxiliary results }\label{sec: proof}

\subsection{Proof of Theorem \ref{theorem: conditional expectation}}\label{subsec: theorem 1}
Let $\tau\in\mathcal{T}(\mathbb{F})$ and $0\leq t\leq T$.
The following lemma includes a characterization of the structure of a general event in $\mathcal{G}_t^\tau$. Intuitively speaking, it states that an event belongs to $\mathcal{G}_t^\tau$ if and only if conditioning on the available information about $\tau$ yields an event which belongs to $\mathcal{G}_t$.
\begin{lemma}\label{lemma: G_t^tau}
	Let $\tau\in\mathcal{T}(\mathbb{F})$ and $0\leq t\leq T$. $A\in\mathcal{G}_t^\tau$ iff there exist $A^0\in\mathcal{G}_t$ and $\left\{A_s;0\leq s\leq t\right\}\subseteq\mathcal{G}_t$ for which

	\begin{equation}\label{eq: event decomposition}
	A=\left(\{\tau>t\}\cap A^0\right)\cup\left[\cup_{s=0}^t\left(\{\tau=s\}\cap A_s\right)\right]\,.
	\end{equation}
\end{lemma}
\begin{proof}
	If $A$ has the structure of \eqref{eq: event decomposition}, then immediately deduce that $A\in\mathcal{G}_t^\tau$ because $\{\tau>t\}\in\mathcal{G}_t^{\tau}$ and $\{\tau=s\}\in\mathcal{G}_t^\tau$ for every $s\in[0,t]$.  In order to prove the other direction,  observe that every
	\begin{equation}
	A\in \mathcal{P}\equiv\left\{B\cap\{\tau\leq s\};B\in\mathcal{G}_t,0\leq s\leq t\right\}\cup \mathcal{G}_t
	\end{equation}
	has the structure of \eqref{eq: event decomposition}. Specifically, if $A=B\cap\{\tau\leq s\}$ for some $B\in\mathcal{G}_t$ and $s\in[0,t]$, then $A_u= B$ for every $u\in[0,s]$ and $A^0=A_{s+1}=\ldots=A_t=\emptyset$. Similarly, if $A\in\mathcal{G}_t$, then it is possible to set $A^0=A_1=\ldots=A_t=A$. In addition, it may be verified that $\mathcal{P}$ is a $\pi$-system. Now, let $\mathds{D}$ be the collection of all $A\in\mathcal{F}$ which have the structure of \eqref{eq: event decomposition} and observe that:
	\begin{enumerate}
		\item An insertion of $A^0=A_0=A_1=\ldots=A_t=\Omega$ implies that $\Omega\in \mathds{D}$.

		\item Assume that $A_1,A_2\in\mathds{D}$ are such that
		\begin{equation}
		A_i=\left(\left\{\tau>t\right\}\cap A_i^0\right)\cup\left[\cup_{s=0}^t\left(\left\{\tau=s\right\}\cap A_{s,i}\right)\right], \ \ i=1,2
		\end{equation}
		and $A_1\subseteq A_2$. As a result, deduce that \begin{equation}
		A_2\setminus A_1=\left(\left\{\tau>t\right\}\cap \widetilde{A}^0\right)\cup\left[\cup_{s=0}^t\left(\left\{\tau=s\right\}\cap \widetilde{A}_{s}\right)\right]
		\end{equation}
		where $\widetilde{A}^0\equiv A_2^0\setminus A_1^0\in\mathcal{G}_t$ and $\widetilde{A}_s\equiv A_{s,2}\setminus A_{s,1}\in\mathcal{G}_t$ for every $s\in[0,t]$. Therefore, deduce that $A_2\setminus A_1\in\mathds{D}$.

		\item Assume that $A_1\subseteq A_2\subseteq\ldots$ is a sequence of events which belong to $\mathds{D}$ such that
		\begin{equation}
		A_i=\left(\left\{\tau>t\right\}\cap A_i^0\right)\cup\left[\cup_{s=0}^t\left(\left\{\tau=s\right\}\cap A_{s,i}\right)\right]\ \ , \ \ i\geq1\,.
		\end{equation}
		Then,
		\begin{equation}
		\cup_{i\geq1}A_i=
		\left(\left\{\tau>t\right\}\cap \widetilde{A}^0\right)\cup\left[\cup_{s=0}^t\left(\left\{\tau=s\right\}\cap \widetilde{A}_{s}\right)\right]
		\end{equation}
		such that $\widetilde{A}^0\equiv\cup_{i\geq1}A^0_i\in\mathcal{G}_t$ and $\widetilde{A}_s\equiv\cup_{i\geq1}A_{s,i}\in\mathcal{G}_t$ for every $s\in[0,t]$. Therefore, deduce that $\cup_{i\geq1}A_i\in\mathds{D}$.
	\end{enumerate}
	This means that $\mathds{D}$ is a Dynkin system such that $\mathcal{P}\subseteq\mathds{D}$ and hence Dynkin's $\pi-\lambda$ theorem yields that $\mathcal{G}_t^\tau=\sigma\left\{\mathcal{P}\right\}\subseteq \mathds{D}$ from which the result follows.
\end{proof}
\newline\newline
Now we are ready to provide the proof of Theorem \ref{theorem: conditional expectation}.
\subsubsection*{Proof of Theorem \ref{theorem: conditional expectation}}
For every $t\leq u\leq T-1$, let $\nu^*_u$ be the minimal-optimal stopping time which solves
\begin{equation}
\text{ess}\sup_{\nu\in\mathcal{T}_u(\mathbb{G}^\tau)}E\left[X^2_\nu\textbf{1}_{\left\{\nu<\tau\right\}}+Y^2_{\tau}\textbf{1}_{\left\{\tau\leq\nu\right\}}\big|\mathcal{G}_u^\tau\right]\,.
\end{equation}
Thus, Lemma \ref{lemma: G_t^tau} implies that for every $t\leq u\leq T-1$, there are $A_u^0\in\mathcal{G}_u$ and $\left(A_{s,u}\right)_{s=0}^u\subseteq\mathcal{G}_u$ such that
\begin{equation}\label{eq: structure}
\{\nu^*_u=u\}=\left(A_u^0\cap\{\tau>u\}\right)\cup\left[\cup_{s=0}^u\left(\{\tau=s\}\cap A_{s,u}\right)\right]\,.
\end{equation}
As a result, deduce that
\begin{equation}\label{eq: u=t}
	\{\tau>t\}\cap\{\nu_t^*=t\}=\{\tau>t\}\cap A^0_t\,.
\end{equation}
Now, consider the case when $u>t$ and observe that for every $0\leq l\leq s$ such that $t\leq s\leq u-1$:
\begin{align}
	&\{\tau>u\}\subseteq\{\tau\neq l\}\,,\label{eq:l inclusion}\\&\{\tau>u\}\cap\{\tau\leq s\}=\emptyset\,.\label{eq:s emptyset}
\end{align}
Therefore, if we use the notation $\overline{A}\equiv\Omega\setminus A$ for every $A\subseteq\Omega$, then an insertion of \eqref{eq: structure} and applying standard algebra of sets imply that:
\begin{align}\label{eq: master decomposition}
&\{\tau>u\}\cap\{\nu_t^*=u\}=\{\tau>u\}\cap\{\nu^*_u=u\}\cap\left(\cap_{s=t}^{u-1}\{\nu_s^*>s\}\right)\\&=\{\tau>u\}\cap\{\nu^*_u=u\}\cap\left(\cap_{s=t}^{u-1}\overline{\{\nu_s^*=s\}}\right)\nonumber\\&=\{\tau>u\}\cap A_u^0\cap\left[\cap_{s=t}^{u-1}\overline{\left[\left(A_s^0\cap\{\tau>s\}\right)\cup\left[\cup_{l=0}^s\left(\{\tau=l\}\cap A_{l,s}\right)\right]\right]}\right]\nonumber\\&=\{\tau>u\}\cap A_u^0\cap\left[\cap_{s=t}^{u-1}\left[\left(\overline{A_s^0}\cup\{\tau\leq s\}\right)\cap\left[\cap_{l=0}^s\left(\{\tau\neq l\}\cup \overline{A_{l,s}}\right)\right]\right]\right]\nonumber\\&=\{\tau>u\}\cap A_u^0\nonumber\\&\cap\left[\cap_{s=t}^{u-1}\left[\left[\overline{A_s^0}\cap\left[\cap_{l=0}^s\left(\{\tau\neq l\}\cup \overline{A_{l,s}}\right)\right]\right]\cup\left[\{\tau\leq s\}\cap\left[\cap_{l=0}^s\left(\{\tau\neq l\}\cup \overline{A_{l,s}}\right)\right]\right]\right]\right]\,.\nonumber
\end{align}
In particular, \eqref{eq:s emptyset} implies that for every $t\leq s\leq u-1$,
\begin{equation}\label{eq: help1}
	\{\tau>u\}\cap\left[\{\tau\leq s\}\cap\left[\cap_{l=0}^s\left(\{\tau\neq l\}\cup \overline{A_{l,s}}\right)\right]\right]\subseteq\{\tau>u\}\cap\{\tau\leq s\}=\emptyset\,.
\end{equation}
In addition, \eqref{eq:l inclusion} implies that for every $0\leq l\leq s$ such that $t\leq s\leq u-1$,
\begin{equation}
\{\tau>u\}\cap\left(\{\tau\neq l\}\cup\overline{A_{l,s}}\right)=\{\tau>u\}\,,
\end{equation}
and hence
\begin{equation}\label{eq: help2}
	\{\tau>u\}\cap\left[\overline{A_s^0}\cap\left[\cap_{l=0}^s\left(\{\tau\neq l\}\cup \overline{A_{l,s}}\right)\right]\right]=\{\tau>u\}\cap \overline{A_s^0}\,.
\end{equation}
Combining \eqref{eq: master decomposition},\eqref{eq: help1} and \eqref{eq: help2} all together yields that
\begin{equation}\label{eq: u>t}
	\{\tau>u\}\cap\{v_t^*=u\}=\{\tau>u\}\cap A^0_u\cap\left(\cap_{s=t}^{u-1}\overline{A^0_s}\right)\,.
\end{equation}
Therefore, \eqref{eq: u=t} and \eqref{eq: u>t} imply that for every $t\leq u\leq T-1$, we get
\begin{equation}
	\{\tau>u\}\cap\{v_t^*=u\}=\{\tau>u\}\cap K_u
\end{equation}
where
\begin{equation}
K_u\equiv A_u^0\setminus\cup_{s=t}^{u-1}A_s^0\,.
\end{equation}
In particular, observe that
\begin{equation}
K_u\in\mathcal{G}_u\ \ , \ \ \forall t\leq u\leq T-1
\end{equation}
and
\begin{equation}
K_{u_1}\cap K_{u_2}=\emptyset\ \ , \ \ \forall t\leq u_1<u_2\leq T-1\,.
\end{equation}
Therefore, it is possible to define a random variable $\hat{\nu}_t:\Omega\rightarrow[t,T-1]$ such that
\begin{equation}
\{\hat{\nu}_t=u\}=K_u\ \ , \ \ \forall t\leq u\leq T-1
\end{equation}
and $\{\hat{\nu}_t=T\}=\overline{\{t\leq\hat{\nu}_t\leq T-1\}}$.
As a result, by construction we get that  $\hat{\nu}_t\in\mathcal{T}_t(\mathbb{G})$ and
\begin{equation}\label{eq: event equality}
\{\tau>u\}\cap\{\nu^*_t=u\}=\{\tau>u\}\cap\{\hat{\nu}_t=u\}\ \ , \ \ \forall t\leq u\leq T-1\,.
\end{equation}
In addition, the identity
\begin{align}
X^2_\nu\textbf{1}_{\left\{\nu<\tau\right\}}+Y^2_{\tau}\textbf{1}_{\left\{\tau\leq\nu\right\}}=\sum_{u=t}^{T-1}X^2_u\textbf{1}_{\{\tau>\nu=u\}}+Y^2_\tau\left(1-\sum_{u=t}^{T-1}\textbf{1}_{\{\tau>\nu=u\}}\right)
\end{align}
holds for every $\nu\in\mathcal{T}_t(\mathbb{G}^\tau)$ (especially, it holds for both $\nu^*_t$ and $\hat{\nu}_t$) and hence \eqref{eq: event equality} implies the first part of \eqref{eq: ess-sup}.

For the proof of the second part of \eqref{eq: ess-sup}, consider some arbitrary $\nu\in\mathcal{T}_t\left(\mathbb{G}^\tau\right)$ and assume that $P(\tau>t)>0$.	Let $A\in\mathcal{G}_t$ and observe that:
\begin{align}\label{eq: conditional expectation11}
& E\left\{\textbf{1}_A\textbf{1}_{\{\tau>t\}}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\right\}\\&=E\left\{\textbf{1}_{\{\tau>t\}}E_{P|\{\tau>t\}}\left[\left(X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\right)\textbf{1}_A\big|\mathcal{G}_t\right]\right\}\nonumber\\&=P\left(\tau>t\right)E_{P|\{\tau>t\}}\left\{E_{P|\{\tau>t\}}\left[\left(X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\right)\textbf{1}_A\big|\mathcal{G}_t\right]\right\}\nonumber\\&=P\left(\tau>t\right)E_{P|\{\tau>t\}}\left\{\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\right]\textbf{1}_A\right\}\nonumber\\&=E\left\{\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\right]\textbf{1}_A\textbf{1}_{\{\tau>t\}}\right\}\ \ , \ \ P\text{-a.s.}\nonumber
\end{align}
In addition,  $\nu\geq t$, $P$-a.s. and hence $P$-a.s., $\tau\leq t$ implies that $\tau\leq\nu$. This and \eqref{eq: conditional expectation11} yield that
\begin{align}\label{eq: conditional expectation}
&E\left\{\textbf{1}_A\left[Y^2_\tau\textbf{1}_{\{\tau\leq t\}}+\textbf{1}_{\{\tau>t\}}E_{P|\{\tau>t\}}\left(X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right)\right]\right\}\\&=E\left\{\textbf{1}_A\left[Y^2_\tau\textbf{1}_{\{\tau\leq t\}}\textbf{1}_{\{\tau\leq\nu\}}+\textbf{1}_{\{\tau>t\}}E_{P|\{\tau>t\}}\left(X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right)\right]\right\}\nonumber\\&=E\left\{\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\right]\textbf{1}_A\right\}\,.\nonumber
\end{align}
Note that in the last equality is valid because of the fact that $\nu\geq t$, $P$-a.s. and hence $\textbf{1}_{\{\nu<\tau\}}\textbf{1}_{\{\tau>t\}}=\textbf{1}_{\{\nu<\tau\}}$, $P$-a.s. Furthermore, notice that \eqref{eq: conditional expectation} holds for every
\begin{equation}
A\in \mathcal{P}\equiv\left\{B\cap\{\tau\leq s\};B\in\mathcal{G}_t,0\leq s\leq t\right\}\cup \mathcal{G}_t\,.
\end{equation}
In addition, let $\mathds{D}$ be the set of all $A\in\mathcal{F}$ for which \eqref{eq: conditional expectation} holds. Then, observe that:

\begin{enumerate}
	\item $\Omega\in\mathcal{G}_t$ and hence \eqref{eq: conditional expectation} implies that $\Omega\in \mathds{D}$.

	\item Assume that $A\in \mathds{D}$. Since $\textbf{1}_{\Omega\setminus A}=\textbf{1}_{\Omega}-\textbf{1}_A$ and $\Omega,A\in\mathds{D}$, then the linearity of conditional expectation yields that $\Omega\setminus A\in \mathds{D}$.

	\item Assume that $\left(A_n\right)_{n=1}^\infty\subset \mathds{D}$ such that $A_i\cap A_j=\emptyset$ for every $i\neq j$. Then, $\textbf{1}_{\cup_{n=1}^\infty A_n}=\sum_{n=1}^\infty\textbf{1}_{A_n}$ and hence the linearity of conditional expectation yields that $\cup_{n=1}^\infty A_n\in \mathds{D}$.
\end{enumerate}
Thus, $\mathds{D}$ is a Dynkin system. Since $\mathcal{P}$ is a $\pi$-system which is a subset of $\mathds{D}$, Dynkin's $\pi-\lambda$ theorem implies that $\mathcal{G}^{\tau}_t=\sigma\{\mathcal{P}\}\subseteq \mathds{D}$. As a result,  the definition of conditional expectation  implies that
\begin{align}\label{eq: conditional expectations}
&E\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t^\tau\right]=Y^2_\tau\textbf{1}_{\{\tau\leq t\}}\\&+\textbf{1}_{\{\tau>t\}}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\ \ , \ \ P\text{-a.s.}\nonumber
\end{align}
Now, recall that it has already been proved  that $\hat{\nu}_t$ solves the essential supremum problem of the LHS of \eqref{eq: conditional expectations} over the domain $\mathcal{T}_t(\mathbb{G}^\tau)$. Thus, $P$-a.s. $\hat{\nu}_t$ solves the essential supremum problem of the RHS of \eqref{eq: conditional expectations} on the same domain, i.e.,
\begin{align}
&E_{P|\{\tau>t\}}\left[X^2_{\hat{\nu}_t}\textbf{1}_{\{\hat{\nu}_t<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\hat{\nu}_t\}}\big|\mathcal{G}_t\right]\\&=\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G}^\tau)}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\ \ , \ \ P\text{-a.s.}\nonumber
\end{align}
Consequently, since $\hat{\nu}_t\in\mathcal{T}_t(\mathbb{G})\subseteq\mathcal{T}_t(\mathbb{G}^\tau)$, deduce that
\begin{align}\label{eq: redundant information}
&E_{P|\{\tau>t\}}\left[X^2_{\hat{\nu}_t}\textbf{1}_{\{\hat{\nu}_t<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\hat{\nu}_t\}}\big|\mathcal{G}_t\right]\\&=\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G})}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\ \ , \ \ P\text{-a.s.}\nonumber
\end{align}
This, with the first equality in \eqref{eq: ess-sup}, implies that the following equation holds $P$-a.s.
\begin{align}
&\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G}^\tau)}E\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t^\tau\right]=Y^2_\tau\textbf{1}_{\{\tau\leq t\}}\\&+\textbf{1}_{\{\tau>t\}}\text{\normalfont ess}\sup_{\nu\in\mathcal{T}_t(\mathbb{G})}E_{P|\{\tau>t\}}\left[X^2_\nu\textbf{1}_{\{\nu<\tau\}}+Y^2_\tau\textbf{1}_{\{\tau\leq\nu\}}\big|\mathcal{G}_t\right]\nonumber
\end{align}
and the result follows. $\blacksquare$

\subsection{Proof of Proposition \ref{lemma: Hamadene}}\label{subsec: proposition 1}
Proposition \ref{lemma: Hamadene} follows from the same arguments which were made by Hamad\'ene \textit{et al.} \cite{Hamadene2010} with certain modifications. Now, we elaborate on these modifications via several lemmata.
	\begin{lemma}
		For each $n\geq1
		$, $\tau_n\in\mathcal{T}(\mathbb{F})$ and $\nu_n\in\mathcal{T}(\mathbb{G}^\tau)$. In addition, both sequences $\left(\tau_n\right)_{n\geq1}$ and $\left(\nu_n\right)_{n\geq1}$ are $P$-a.s. nonincreasing.
	\end{lemma}

	\begin{proof}
		For $n=1$, the claim is immediate. Assume that it holds for some $n\geq1$ and notice that this assumption yields that $\nu_n\in\mathcal{T}(\mathbb{F})$. Now,
		define

		\begin{equation}
		U^{1,n}_t\equiv X_t^1\textbf{1}_{\{t\leq\nu_n\}}+\left(X^1_T\textbf{1}_{\{\nu_n=T\}}+Y^1_{\nu_n}\textbf{1}_{\{\nu_n<T\}}\right)\textbf{1}_{\left\{t\geq\nu_n\right\}} , \ \ \forall t\in[0,T]
		\end{equation}
		and $(U^{1,n}_t)_{t=0}^T$ satisfies the pre-conditions of Theorem 1.2 in \cite{Peshkir2006} w.r.t. the filtration $\mathbb{F}$. Thus,  deduce that $\widetilde{\tau}_{n+1}$ is the minimal-optimal  stopping time of the problem
		\begin{equation}
		\text{ess}\sup_{\tau\in\mathcal{T}(\mathbb{F})}E\left(U_\tau^{1,n}|\mathcal{F}_0\right)\,,
		\end{equation}
		and hence $\widetilde{\tau}_{n+1}\in\mathcal{T}(\mathbb{F})$. Then, showing that $\tau_{n+1}\in\mathcal{T}(\mathbb{F})$  follows by the same arguments which appear in the proof of Lemma 3.1 in \cite{Hamadene2010}. Similarly, define

		\begin{equation}
		U^{2,n}_t\equiv X_t^2\textbf{1}_{\{t<\tau_n\}}+Y_t^2\textbf{1}_{\{t\leq\tau_n\}} , \ \ \forall t\in[0,T]\,,
		\end{equation}
		and notice that the current model assumptions  yield that $(U^{2,n}_t)_{t=0}^T$  satisfies the pre-conditions of Theorem 1.2 in \cite{Peshkir2006} w.r.t. $\mathbb{G}^{\tau_n}$. Then, showing that $\widetilde{\nu}_{n+1}\in\mathcal{T}(\mathbb{G}^{\tau_n})$ and hence $\nu_{n+1}\in\mathcal{T}(\mathbb{G}^{\tau_n})$ follows by the same arguments which appear in the proof of Lemma 3.1 in Hamad\'ene \textit{et al.}  \cite{Hamadene2010}. Finally, the monotonicity of the sequences $\left(\tau_n\right)_{n\geq1}$ and $\left(\nu_n\right)_{n\geq1}$ is justified by construction and it is explained in the proof of Lemma 3.1 in Hamad\'ene \textit{et al.}  \cite{Hamadene2010}.
	\end{proof}

	\begin{lemma}
		For every $n\geq2$, one has:
		\begin{enumerate}
			\item \begin{equation}
			J_1(\tau,\nu_n)\leq J_1(\tau_n,\nu_n), \ \ \forall\tau\in\mathcal{T}(\mathbb{F})\,.
			\end{equation}

			\item \begin{equation}
			J_2(\tau_n,\nu)\leq J_2(\tau_n,\nu_n), \ \ \forall\nu\in\mathcal{T}(\mathbb{G})\,.
			\end{equation}
		\end{enumerate}
	\end{lemma}

	\begin{proof}
		The proof is the same as the  proof of Lemma 3.3 in \cite{Hamadene2010}.  The only thing that should be noticed (for the proof of the second part) is that $\nu\in\mathcal{T}(\mathbb{G})$ implies that $\nu\wedge\tau_n\in\mathcal{T}(\mathbb{G}^{\tau_n})$.
	\end{proof}
	\newline\newline
	Observe that the discrete-time modelling with the condition \eqref{eq: integrability} in the definition of $\mathcal{D}$ implies that a standard convergence theorem argument might be used in order to justify the following convergences:
	\begin{enumerate}
		\item $J_1(\tau_n,\nu)\rightarrow J_1(\tau^*,\nu)$ as $n\to\infty$  for every $\nu\in\mathcal{T}(\mathbb{G})$.

		\item $J_2(\tau,\nu_n)\rightarrow J_2(\tau,\nu^*)$ as $n\to\infty$ for every $\tau\in\mathcal{T}(\mathbb{F})$.

		\item For each $i=1,2$, $J_i(\tau_n,\nu_n)\rightarrow J_i(\tau^*,\nu^*)$ as $n\to\infty$.
	\end{enumerate}
	Finally, combining these results with the previous lemmata in this subsection yields the proof of Proposition \ref{lemma: Hamadene}. $\blacksquare$

	\subsection{Proof of Corollary \ref{cor: equivalence}}\label{subsec: corollary 1}
	Observe that the first direction is an immediate consequence of \eqref{eq: nu*}. In order to prove the other direction, assume that $\theta^*<\tau^*$. Thus, there exists $n'\geq1$ (which is a random variable) such that $\theta_{n'+1}<\tau_n$ for every $n\geq1$. In addition, according to \eqref{eq: tilde nu}, $\nu_{n'+1}=\theta_{n'+1}$, $P$-a.s. Recall that $(\tau_n)_{n\geq1}$ and $(\nu_n)_{n\geq1}$  are nonincreasing sequences. As a result, by construction, $\tau_{n'}=\tau_{n'+1}=\ldots$ and hence
	\begin{equation}
	\tau^*=\lim_{n\to\infty}\tau_n>\theta_{n'+1}\geq\nu^*, \ \ P\text{-a.s.}
	\end{equation}
	which makes \eqref{eq: iff} follows. $\blacksquare$

	\subsection{Proof of Lemma \ref{lemma: auxilary}}\label{subsec: lemma 1}
	The best response of Player 2 is the same in both games and hence it is sufficient to show that $\tau_0$ is a best response of Player 1 to $\nu_0$ in $\Gamma$. To this end, observe that \eqref{eq: sufficient condition} yields that $X^1_t+\epsilon>Y_t^1$ for every $0\leq t\leq T$. Therefore, since Player 1 has a priority in stopping and $\tau_0$ is a best response of Player 1 to $\nu_0$ in $\Gamma_\epsilon$, deduce that $\tau_0\leq\nu_0$, $P$-a.s. (otherwise $\tau_0\wedge\nu_0$ is a better response). Now, assume by contradiction that there exists $\tau'\in\mathcal{T}(\mathbb{F})$ such that $J_1(\tau_0,\nu_0)<J_1(\tau',\nu_0)$. Then, \eqref{eq: sufficient condition} with the priority in stopping of Player 1 implies that $\tau'\wedge\nu_0$ is also a better response of Player 1, i.e.,
	\begin{align}
	J_{1,\epsilon}(\tau'\wedge\nu_0,\nu_0)-\epsilon&=J_1(\tau'\wedge\nu_0,\nu_0)\\&\geq J_1(\tau',\nu_0)>J_1(\tau_0,\nu_0)=J_{1,\epsilon}(\tau_0,\nu_0)-\epsilon\nonumber
	\end{align}
	from which a contradiction follows. $\blacksquare$
	\begin{equation*}
	\end{equation*}
\textbf{Acknowledgement:}
The author would like to express his deep gratitude to Yan Dolinsky and Eilon Solan for insightful discussions about the problem as well as for their valuable comments on early versions of the current work. In addition, special thanks are devoted to Reviewer 1 for detecting a mathematical gap in the original manuscript and to Reviewer 2 whose suggestions lead to some major changes in the structure.