\section{Introduction}

Some preliminary work regarding Finsler metrics was carried out in \cite{Chavez2021} where the authors considered a lake in the form of the unit disks $\mathbb{B}^2,$ with a concentric and symmetric wind current given by the vector field $W(x_1,\,x_2)=(-x_1,\,-x_2)$. The distance function (time of displacement) in this framework, is given by
\begin{align}\label{111}
	d_F(P,Q)=\ln\left(\frac{\sqrt{\langle P, Q-P \rangle^2+(1-\Vert P\Vert^2)\Vert Q-P\Vert^2}-\langle P, Q-P \rangle}{\sqrt{\langle P, Q-P\rangle^2+(1-\Vert P\Vert^2)\Vert Q-P\Vert^2}-\langle Q, Q-P \rangle}
	\right).
\end{align}
In \cite{Chavez2021} were also obtained equations for the circle, and formulas for the distance from a point to a line and from a line to a point.\par

Some questions naturally arise in this model: Imagine an isle and a coast in the form of a straight line. Suppose a boat parting from the isle to a point in the beach (or the contrary), but it has to stop in some middle point along the way to fill its tank with fuel.\par
Now imagine two boats wanting to meeting in some point at the ocean, one coming from the beach (of the coast in straight line ) and the other one from the isle, both traveling with the same speed and arriving at the same time, then they should meet in a mid point between them.\par
Finally it can be thought of two boats parting from the same point, one arriving to an isle, the other to the beach at the same time.\par
The physics problems presented above motivate the study of an equivalent to parabolas in the Euclidean plane, that we will call Funk's Parabolas (see Definition \ref{def_parabola}). Note that in each problem we must consider the two directions because, due to the wind, there is no symmetry going from one point to another.\par
Once we fix the external force as been concentric and symmetric, the problems aforementioned involve situations that can be solved with the help of Funk's Parabolas.

\section{Preliminary}

In this section some definitions and results, required for the development of our work, are introduced. We adopted the definitions given in \cite{Chavez2021}, which can also be found in \cite{Manfredo1}, such as inner product, norm, regular curve, arc length (which will be called usual or Euclidean) and vector field.

\begin{defi}\label{FunkMetric}{\em
	Let $ x=(x_1,x_2)\in\Omega=\mathbb{B}^2$ and $y=(y_1,y_2)\in\mathbb{R}^2.$
	The function
	\begin{align*}
		F=&\sqrt{a_{11}(x)y_1^2+ 2a_{12}(x)y_1y_2 + a_{22}(x)y^2_2}+b_1(x)y_1+b_2(x)y_2,
	\end{align*}
	where $b_1= \dfrac{x_1}{1-x_1^2-x_2^2}$, $b_2= \dfrac{x_2}{1-x_1^2-x_2^2}$
	and
	\begin{displaymath}
		[a_{ij}]=\frac{1}{(1-x_1^2-x_2^2)^2}\left(\begin{array}{cc}
			{1-x_2^2}  & {x_1x_2} \\
			{x_1x_2} & {1-x_1^2}
		\end{array}\right),
	\end{displaymath}
	it is called \textit{Funk metric on the unit disks} $ \mathbb{B}^2=\{x\in\mathbb{R}^2;\;\Vert x\Vert<1\}$.}
\end{defi}

This metric models a boat sailing with a unitary speed in $\mathbb{B}^2$, where a wind current given by $W_x=(-x_1,-x_2)$ is present with speed lower that $1$. (see Section 3 in \cite{Chavez2021})

The Funk metric on $\mathbb{B}^2, $ is a special case of Randers metrics. Some properties where studied in \cite{Shen2001} for more general cases.

\begin{defi}[Funk's type Arc Length]{\em
	Let $c:[a,b]\rightarrow \Omega\subset\mathbb{R}^2$ a piecewise regular curve. The \textit{arc length (Funk's type) of $c$} is defined by
	\[ \mathscr{L}_F(c):=\int_a^bF(c(t),c'(t))dt.\]
	For any points $p, q \in \Omega \subset \mathbb{R}^2,$ we define the \textit{distance} from $p$ to $q$ induced by $F$, as
	\[ d_F(p,q):=\operatorname{inf}_c\mathscr{L}_F(c),\]
	where the infimum is taken over the set of all piecewise regular curves $c$ such that $c(a)=p$ and $ c(b)=q $.}
\end{defi}

It is known that the \textit{shortest paths} in this metric are straight lines (see Example 9.2.1 in \cite{Shen2001}), i.e., the optimal path is always a straight line. This property makes this metric manageable in relation to distance calculation. More details about geodesics and its relation with shortest paths can be found in Section 3.2 in \cite{ChernRF2005}, and  Section 2.3 in \cite{Cheng2012}.\par

In Remark 5.1 in \cite{Chavez2021} was proved that $d_F$ given by \eqref{111} is non-reversible ($ d_F(P,Q)\neq d_F(Q,P) $) and it is non invariant by translations, but it is invariant by rotations.

If we consider \begin{align}\label{eq:def_r}
r=&\frac{\sqrt{k}-\langle P,Q-P\rangle}{\sqrt{k}-\langle Q,Q-P\rangle},
\end{align} where
\[ k = \langle P, Q-P \rangle^2+(1-\Vert P\Vert^2)\Vert Q-P\Vert^2,\]
the equation \eqref{111} can be more manageable using the following version of Theorem 5.1 in \cite{Chavez2021}.
\begin{theo}\label{maintheorem}{\em
	Let $P,\;Q$ be points in $\mathbb{B}^2$ and $r\geq1$ a real number, then:
	\begin{align}\label{13}
		\left\Vert \dfrac{P}{r}-Q\right\Vert&=\dfrac{r-1}{r}
	\end{align}
	it is equivalent to
	$d_F (P, Q) = \ln r.$
}
\end{theo}

It is worth mentioning that the previous theorem can be easily generalized for points $P$ and $Q$ in $\mathbb{R}^n$ at the unit ball (see the proof of the theorem in \cite{Chavez2021}).
In order to facilitate calculations we define the function $ \rho:[0,1)\times [0,1)\to \mathbb{R} $ as,
\begin{align}\label{rxy}
\rho(\epsilon, \eta)=\frac{1+\operatorname{sgn}(\epsilon-\eta)\epsilon}{1+\operatorname{sgn}(\epsilon-\eta)\eta}
\end{align}
where $\operatorname{sgn}$ is the function sign:
\[
\operatorname{sgn}(x) = \left\{
\begin{array}{rll}
1, & \hbox{if} & x>0\\
0, & \hbox{if} & x=0\\
-1, & \hbox{if} & x<0.
\end{array}
\right.\]
The following formulas can be found in Section 5.2 e 5.3 in \cite{Chavez2021}, and play a fundamental role in the present work:

\begin{propo}[\cite{Chavez2021}]\label{dist_ponto_reta_const}{\em
Let $P=(a,b)$ be a point and $s:\;x_2=c$ a constant line. The Funk distance $d_F(P,s)$ from the point $P$ to the line $s$ and Funk distance $d_F(s,P)$ from the line $s$ to the point $P,$ are, respectively
\begin{align}
	d_{F}(P,s)=&d_F\left(P,\left(\frac{a}{\rho(b,c)},c\right)\right)=\ln\left(\rho(b,c)\right),\label{dist_pt_reta_const}\\
	d_{F}(s,P)=&d_F\left(\left(a\cdot \rho(c,b),c\right),P\right)=\ln\left(\rho(c,b)\right),\label{dist_reta_const_pt}
\end{align}

where $\rho(\epsilon,\eta)$ is given by \eqref{rxy}.
}\end{propo}

\section{Funk's Parabolas in $\mathbb{B}^2$}

As the Funk distance is non-reversible, we obtain for types of parabolas.

\begin{defi} \label{def_parabola}{\rm

Let $s \subset \mathbb{B}^2$ be a fixed line, called \textit{the directrix}, and $\mathcal{F}$ a point in $\mathbb{B}^2 \subset \mathbb{R}^2$, called \textit{the focus}. We define a \textit{Funk's Parabola}, with directrix $s$ and focus $\mathcal{F}$, as the set of points $P \in \mathbb{B}^2$ that satisfies one of the following properties:
\begin{multicols}{2}
\begin{enumerate}
	\item[(P1)] $d_F(\mathcal{F},P)=d_F(s,P);$
	\item[(P2)] $d_F(P,\mathcal{F})=d_F(P,s);$
	\item[(P3)] $d_F(\mathcal{F},P)=d_F(P,s);$
	\item[(P4)] $d_F(s,P)=d_F(P,\mathcal{F}).$
\end{enumerate}
\end{multicols}
}
\end{defi}
The set of point satisfying (P$k$) is called {\bf Funk's parabola of type {\boldmath$k$}} with $k=1,2,3,4$.

\begin{remark} Note that in the previous definition we are also considering the degenerate case, i.e., when the focus $\mathcal{F}$          belongs to the directrix $s$.
\end{remark}
Due to the invariance of the Funk distance by rotations, we will consider first {\bf the case where the directrix is parallel to the  {\boldmath$x$} axis}.\par
Throughout this work we will use the following notations
\[ s: y=y_0,\quad \mathcal{F}=(f_0,\,g_0)\;\text{ and }\;P=(x,\,y). \]
\subsubsection{Funk's Parabola type 1 - Particular Case}
By Theorem \ref{maintheorem} and \eqref{dist_reta_const_pt} we have that $P=(x,\,y) $ satisfies $P1$ if, and only if, \begin{align}\label{eq:tipe1norm}
\displaystyle\left\Vert \frac{\mathcal{F}}{\rho(y_0,y)}-P\right\Vert^2 = \left\Vert \frac{Q}{\rho(y_0,y)}-P \right\Vert^2,
\end{align}
where $ Q=\left(x.\rho(y_0,y),y_0\right)$ and $\rho$ is given by \eqref{rxy}.

First, suppose the degenerate case, when the focus $\mathcal{F}=(f_0,g_0)$ is in the line $s=y_0$, i.e., when $y_0=g_0$. Therefore the equation \eqref{eq:tipe1norm} is equivalent to \[x-f_0\cdot (\rho(y_0,y))^{-1}=0.\]
This implies that $P=(x,\,y)$, is either in the line segment $ (1-y_0)x-f_0(1-y)=0, $ for $ y\geq y_0, $ or in the line segment $(1+y_0)x-f_0(1+y)=0, $ for $ y<y_0$. If $f_0=0$, then $ P $ is in the line $x=0$. If $f_0\neq 0$ and $ x $ goes to zero, thus $ y $ approximates to $1$ or $-1$ independently of $f_0$ value.\par
In what follows we characterize non degenerate Funk's Parabola of type 1.
\begin{propo}
Let   $s: y=y_0 \subset \mathbb{B}^2$ be the directrix, $\mathcal{F}=(f_0,g_0) \in \mathbb{B}^2$ be the focus with $ g_0\neq y_0$ and $ \mathcal{E} \subset \mathbb{R}^2$ an Euclidean ellipse described by
\begin{align}\label{p1}
x^2 +B{x}\overline{y} +C{\overline{y}^2} +E{\overline{y}} = 0,
\end{align}
where
\begin{align}
\overline{y}=& 1+\operatorname{sgn}(y_0-g_0)y\\
B=&-\frac{2f_0}{\sigma_1},\label{defB}\\
C=&\frac{f_0^2 + (y_0-g_0)^2+2\vert y_0-g_0\vert}{\sigma_1^2},\label{defC}\\
E=&-\frac{2\vert y_0-g_0\vert}{\sigma_1},\label{defE}
\end{align}
and
\begin{align}\label{defsigma}
\sigma_1= 1+\operatorname{sgn}(y_0-g_0)y_0.
\end{align}
Hence the locus of points of the Funk's Parabola of type 1 is given by
\[ \mathcal{E} - \{(0,-\operatorname{sgn}(y_0-g_0))\}. \]
\end{propo}

\begin{proof}
The line $s$ divides $\mathbb{B}^2$ in two domains. Suppose $\mathcal{F}\notin s$. We will show that the Funk's Parabola of type 1 and the focus $\mathcal{F}$ are in the same subspace generated by $s$. In fact, suppose that $\mathcal{F}$ and $P$ are in different subspaces generated by $s$.  Let $Q$ and $Q_0$ be in $s$, such that $d_F(s,P)=d_F(Q,P) $ and $ Q_0$ is in the line segment $ \overline{\mathcal{F}P}, $ then
\[d(s,P)=d(Q,P)\leq d(Q_0,P)<d(\mathcal{F},Q_0)+ d(Q_0,P) = d(\mathcal{F},P).\]
Suppose that there is a point $P=(x,\,y)$ in the Funk's parabola of type 1 in the line $s$, i.e., there is a $y$ such that $y=y_0$, hence by (P1) we have that $ \mathcal{F}=P$.
Therefore, the sign of $ y_0-y $ is equal to the sign of $ y_0-g_0 $, thus we have that $\rho(y_0,y)$ can be written as
\[\rho(y_0,y)=\frac{1+\operatorname{sgn}(y_0-g_0)y_0}{1+\operatorname{sgn}(y_0-g_0)y}.\]
Using the properties of absolute value and inner product, equation \eqref{eq:tipe1norm} can be written as
\begin{align}\label{p0}
2\rho(y_0,y)\langle P, \mathcal{F}-Q\rangle + \Vert Q\Vert^{2} - \Vert \mathcal{F}\Vert^{2} = 0.
\end{align}
Due to the fact that $ \operatorname{sgn}(y_0-g_0)\cdot(y_0-g_0)=\vert y_0-g_0\vert $, and replacing $ \mathcal{F}, P$ and $ Q$, we have that equation \eqref{p0} is equivalent to:
\begin{align*}
x^2 -2f_0\frac{x}{\rho(y_0,y)} +\left[f_0^2 + (y_0-g_0)^2+2\vert y_0-g_0\vert\right]\frac{1}{\rho^2(y_0,y)} - 2\vert y_0-g_0\vert\frac{1}{\rho(y_0,y)} = 0,
\end{align*}
which is just the equation \eqref{p1}.

Equation \eqref{p1} is a second-degree equation in $x$ and $1+\operatorname{sgn}(y_0-g_0)y$, thus in $x$ and $y$, and whose discriminant $I$ is given by
\[I=B^2-4C=-\frac{4}{\sigma_1^2}[(y_0-g_0)^2+2\vert y_0-g_0\vert]<0 \]
what show us that \eqref{p1} is elliptical in Euclidean geometry (i.e., the empty set, an ellipse, a point or a circle in $\mathbb{R}^2 $).\par
Note that the point $(0,-\operatorname{sgn}(y_0-g_0))$ satisfies the equation \eqref{p1}. Then, the locus of \eqref{p1} can not be a closed curve. We will see that in fact it is an Euclidean ellipse without the point $(0,-\operatorname{sgn}(y_0-g_0))$.\par
Techniques for determine the conic from a given general equation of second degree can be found in \cite{EdwardsPenney1990,Leithold1976}.\\
{\bf Case 1:} {\boldmath$f_0=0$.}
When $f_0=0$, then the change of variables
\begin{align}
x=&\overline{x},\label{xf0}\\
y=&\left(y-\frac{E}{2C}-1\right)\operatorname{sgn}(y_0-g_0)\label{yf0}
\end{align}
where $\frac{E}{2C} = \frac{\sigma}{2+\vert y_0-g_0\vert}$, transforms the equation \eqref{p1} in the canonical form:
\begin{align}\label{canonicaf0}
\frac{\overline{x}^2}{\frac{E^2}{4C}} + \frac{\overline{y}^2}{\frac{E^2}{4C^2}}=1.
\end{align}
From \eqref{xf0},  \eqref{yf0} and \eqref{canonicaf0}, we have that the Funk's parabola of type 1 is a section of an Euclidean ellipse with center $ \mathcal{C}=\left(0, \frac{g_0-\operatorname{sgn}(y_0-g_0)}{2+\vert y_0-g_0\vert}\right) $ and vertices:
$ (\pm \frac{E}{2\sqrt{C}},-(1+\frac{E}{2C})\operatorname{sgn}(y_0-g_0)), (0, -(1+\frac{E}{2C}\pm\frac{E}{2C})).$\\
{\bf Case 2:} {\boldmath$f_0\neq 0$}. In \eqref{p1}, consider the change of variable
\begin{align}
x=&\alpha \overline{x} - \beta \overline{y}-\left(\alpha\frac{D}{2A} + \beta\frac{E}{2C}\right),\label{eq:xf0neq0}\\
1+\operatorname{sgn}(y_0-g_0)y=&\beta \overline{x} + \alpha\overline{y} - \left(\beta\frac{D}{2A}+\alpha\frac{E}{2C}\right),
\end{align} where $ \alpha=\frac{1}{\sqrt{2}}\sqrt{1-\operatorname{sgn}(f)\frac{(1-C)}{\sqrt{B^2+(1-C)^2}}}, $ $ \beta=\frac{1}{\sqrt{2}}\sqrt{1+\operatorname{sgn}(f)\frac{(1-C)}{\sqrt{B^2+(1-C)^2}}}$ we have that equation \eqref{p1} is reduced to the canonical form:
\begin{align}
\frac{\overline{x}^2}{\frac{1}{4}\left(\frac{\overline{D}^2}{\overline{A}^2} + \frac{\overline{E}^2}{\overline{A}\overline{C}}\right)} + \frac{\overline{y}^2}{\frac{1}{4}\left(\frac{\overline{D}^2}{\overline{A}\overline{C}} + \frac{\overline{E}^2}{\overline{C}^2}\right)}=1,
\end{align}
where
\begin{align}
\overline{A}=&\alpha^2 + B\alpha\beta + C\beta^2,\\
\overline{C}=&\beta^2-B\alpha\beta + C\alpha^2,\\
\overline{D}=&E\beta,\\
\overline{E}=&E\alpha.\label{eq:overE}
\end{align}
Note that $ \overline{A}=\left(\alpha + \frac{\beta}{2}B\right)^2  + \frac{(y_0-g_0)^2+2\vert y_0-g_0\vert}{\sigma_1^2}>0.$ Analogously, we can show that $ \overline{C}>0$.
From \eqref{eq:xf0neq0}-\eqref{eq:overE}, we have that the Funk's parabola of type 1 is a section of an Euclidean ellipse with center in $ \mathcal{C}=-\left(\alpha\frac{D}{2A} + \beta\frac{E}{2C}, \operatorname{sgn}(y_0-g_0)\left(1+ \beta\frac{D}{2A}+\alpha\frac{E}{2C}\right)\right) $ and vertices
$ \mathcal{C}\pm \frac{1}{2}\sqrt{\frac{\overline{D}^2}{\overline{A}^2}+\frac{\overline{E}^2}{\overline{A}\overline{C}}}\left(\alpha,\beta\operatorname{sgn}(y_0-g_0)\right),  $
$ \mathcal{C}\pm \frac{1}{2}\sqrt{\frac{\overline{D}^2}{\overline{A}\overline{C}}+\frac{\overline{E}^2}{\overline{C}^2}}\left(-\beta,\alpha\operatorname{sgn}(y_0-g_0)\right). $
\end{proof}

\subsubsection{Funk's Parabola type 2 - Particular Case}
By Theorem \ref{maintheorem} and \eqref{dist_pt_reta_const} we have that the point $P=(x,\,y)$ satisfies (P2) if, and only if, \begin{align}\label{eq:tipe2norm}
\displaystyle\left\Vert \frac{\mathcal{F}}{\rho^{-1}(y,y_0)}-P\right\Vert^2 = \left\Vert \frac{Q}{\rho^{-1}(y,y_0)}-P \right\Vert^2,
\end{align}
where $ Q=\left(x.\rho^{-1}(y,y_0),y_0\right)$ and $\rho $ is given by \eqref{rxy}.\par
The development of equations \eqref{eq:tipe1norm} and \eqref{eq:tipe2norm} are analogous. Hence we have, first, the degenerate case of the Funk's parabola of type 2 $(y_0=g_0)$, is characterized by the line
\[x-f_0\rho(y,y_0)=0,\] which is just a reflection of a Funk's parabola of type 1 (degenerate) in relation to the line $s$.\par
Second, the non-degenerate Funk's parabolas of type 2 are characterized by the following equation.
\begin{propo}\label{prop1}
Let $s: y=y_0 \subset \mathbb{B}^2$ be a directrix, $\mathcal{F}=(f_0,g_0) \in \mathbb{B}^2$ the focus with $g_0\neq y_0$ and $ \mathcal{H} \subset \mathbb{R}^2$ an Euclidean hyperbola characterized by the equation.
\begin{align}\label{p2}
x^2 +B{x}{\overline{y}} +C{\overline{y}^2} +E{\overline{y}} = 0,
\end{align}
where
\begin{align}
\overline{y}=&1-\operatorname{sgn}(y_0-g_0)y\\
B=&-\frac{2f_0}{\sigma_2},\label{defB2}\\
C=&\frac{f_0^2 + (y_0-g_0)^2-2\vert y_0-g_0\vert}{\sigma_2^2},\label{defC2}\\
E=&\frac{2\vert y_0-g_0\vert}{\sigma_2},\label{defE2}
\end{align}
and
\begin{align}\label{defsigma2}
\sigma_2= 1-\operatorname{sgn}(y_0-g_0)y_0.
\end{align}
Then the locus of the Funk's parabola of type 2 is the convex curve \[ \mathcal{H} \cap \mathbb{B}^2. \]
\end{propo}
\begin{proof}
Note that the discriminant $I=B^2-4C$ is given by
\[I=4\vert y_0-g_0\vert(2-\vert y_0-g_0\vert) >0.\]
what shows us that the second-degree equation \eqref{p1} is hyperbolic in the Euclidean geometry (i.e., two concurrent lines, or a hyperbola in $ \mathbb{R}^2 $). As $ E\neq 0$, the two lines option is disregarded.\par
Analogous to what was done in Proposition \ref{prop1}, it is possible to show that the focus $\mathcal{F}$ and the points $P$ of the Funk's parabola of type 2 are in an open subspace generated by the directrix $s:y=y_0$. Therefore, only one of the curves of the Euclidean hyperbola intercepts the domain $ \mathbb{B}^2$. It remains to show that the set $ \mathcal{H} \cap \mathbb{B}^2$ its not empty, in fact, let $ Q*\in s $ be such that $d(F,s)=d(F,Q*)$, and $f:[0,1]\to \mathbb{R}$ the function defined by
\[f(t)=d(P_t,\mathcal{F})-d(P_t,Q*),\]
where $ P_t=\mathcal{F}+ t(Q*-\mathcal{F})\in \mathbb{B}^2$. The function $f$ is continuous and $f(0)=-d(P,Q*)<0, $ $ f(1)=d(Q*,\mathcal{F})>0$. By the Intermediate Value Theorem , we have that there exist $t_0\in (0,1)$ such that $f(t_0)=0.$
\end{proof}
\subsubsection{Examples - Funk's Parabolas type 1 and 2}
\begin{exa}
\textbf{Funk's parabolas of type 1 and 2 passing through the origin}
Let be the parabolas of type 1 and 2 passing through $(0,0)$, whose equations are given by \eqref{p1} and \eqref{p2}.\par
Substituting  $x=y=0$ in \eqref{p1} or \eqref{p2}, we obtain $C = -E$, in both cases, and we have that $C = -E$ if, and only if $\Vert \mathcal{F}\Vert^2 = y_0^2$. This equations says that, the Funk's parabolas of type 1 and 2 pass through the origin if, and only if, the Euclidean distance of $\mathcal{F}$ and $s$ to the origin are the same.
\begin{figure}[H]
\centering
\caption{Funk's parabola of type 1 with $\mathcal{F} = (0.3, -0.4)$ and $y_0 = 0.5$.}
\includegraphics[width=0.5\textwidth]{Parabola_tipo_1_Ndegenerado_ex01.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's parabola of type 1 with $\mathcal{F} = (0.3, 0.4)$ and $y_0 = - 0.5$.}
\includegraphics[width=0.5\textwidth]{Parabola_tipo_1_Ndegenerado_ex02.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's parabola of type 2 with $\mathcal{F} = (0.3, - 0.4)$ and $y_0 = 0.5$.}
\includegraphics[width=0.5\textwidth]{Parabola_tipo_2_Ndegenerado_ex03.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's parabola of type 2 with $\mathcal{F} = (0.3, 0.4)$ and $y_0 = - 0.5$.}
\includegraphics[width=0.5\textwidth]{Parabola_tipo_2_Ndegenerado_ex04.png}

\end{figure}
\end{exa}
\begin{exa}
\textbf{Funk's parabolas of type 1 and 2 with focus in the origin.}

\end{exa}

\subsubsection{Funk's Parabola type 3 - Particular Case}

By Theorem \ref{maintheorem} and \eqref{dist_pt_reta_const} we have that the point $P$ satisfies (P3) if, and only if,
\begin{align}\label{3p0}
\left\Vert\frac{\mathcal{F}}{\rho(y,y_0)}-P\right\Vert^{2}=\left\Vert\frac{P}{\rho(y,y_0)}-Q\right\Vert^{2},
\end{align}
where $ Q=(\frac{x}{\rho(y,y_0)},y_0)\in s $ and $ \rho $ is given by \eqref{rxy}.
Using the properties of absolute value and inner product, we obtain the following equation
\eqref{3p0}, in the coordinates $\mathcal{F}=(f_0,g_0),\, P=(x,\,y)$ and $ s: y=y_0, $ is given by:
\begin{align}\label{3p1}
\overline{y}^4 + x^2\overline{y}^2-2\overline{y}^3+A\overline{y}^2 + Bx\overline{y} + C\overline{y}+D=0,
\end{align}
where  \begin{align*}
\overline{y}=&1-\operatorname{sgn}(y_0-y)\cdot y,\\
A=&2\sigma_3g_0\operatorname{sgn}(y_0-y),&
B=&-2\sigma_3 f_0,\\
C=&2\sigma_3(1-\operatorname{sgn}(y_0-y)\cdot g_0)&
D=&\sigma_3^2(f_0^2+g_0^2-1),
\end{align*}
and $ \sigma_3=1-\operatorname{sgn}(y_0-y)\cdot y_0. $

\subsubsection{Funk's Parabola type 4 - Particular Case}
By Theorem \ref{maintheorem} and \eqref{dist_pt_reta_const} we have that the point $P$ satisfies (P4) if, and only if,
\begin{align}\label{4p0}
\left\Vert\frac{\mathcal{F}}{\rho^{-1}(y_0,y)}-P\right\Vert^{2}=\left\Vert\frac{P}{\rho^{-1}(y_0,y)}-Q\right\Vert^{2},
\end{align}
where $ Q=(\frac{x}{\rho^{-1}(y_0,y)},y_0)\in s $ and $ \rho $ is given by \eqref{rxy}.
Using the properties of absolute value and inner product, we obtain the equation \eqref{4p0}, in the coordinates $ \mathcal{F}=(f_0,g_0),\, P=(x,\,y) $ and $ s: y=y_0, $ is given by:
\begin{align}\label{4p1}
\overline{y}^4 + x^2\overline{y}^2-2\overline{y}^3+A\overline{y}^2 + Bx\overline{y} + C\overline{y}+D=0,
\end{align}
where \begin{align*}
\overline{y}=&1+\operatorname{sgn}(y_0-y)\cdot y,\\
A=&-2\sigma_4g_0\operatorname{sgn}(y_0-y),&
B=&-2\sigma_4 f_0,\\
C=&2\sigma_4(1+\operatorname{sgn}(y_0-y)\cdot g_0),&
D=&\sigma_4^2(f_0^2+g_0^2-1),
\end{align*}
and $ \sigma_4=1+\operatorname{sgn}(y_0-y)\cdot y_0. $

\begin{remark}\label{obs34equiv}
Note that the equation for parabolas of type 3 when $y < y_0$ (resp. $y>y_0$) is equal to the equation for parabolas of type 4 when $ y > y_0$ (resp. $ y<y_0 $).
\end{remark}

\subsubsection{Further Analysis of Funk's Parabolas type 3 and 4}

Due to Remark \ref{obs34equiv}, the graphs of parabolas of type 3 and 4 are equal, then we will study only the equation  \eqref{3p1}.
As can be noticed, the equation \eqref{3p1} is a quartic equation. Now, we will show that the right side of \eqref{3p1} is irreducible.
In order to do this, lets suppose that the equation \eqref{3p1} can be written as
\begin{align*}
\left(\overline{y}^2+Gx\overline{y} + Hx^2 + Ix+ J\overline{y} + K\right)\times\left(\overline{y}^2 + Lx\overline{y} + Mx^2+Nx + P\overline{y} + Q\right) = 0.
\end{align*}
Multiplying and grouping terms with $x$ and $y$ above we have
\begin{align*}
0=&\overline{y}^4+(L+G)\overline{y}^3x + (M+H+GL)\overline{y}^2x^2+(MG+HL)\overline{y}x^3 + HMx^4+\\
&+(N+I+GP+JL)\overline{y}^2x + (P+J)\overline{y}^3 + (NG + JM+HP+IL)\overline{y}x^2 + \\
&+(HN+IM)x^3 + (Q+K+JP)\overline{y}^2 + (GQ+IP+LK+NJ)x\overline{y} + \\
&+(QH+NI+MK)x^2 + (JQ+PK)\overline{y} + (KN+IQ)x + KQ.
\end{align*}
By comparing the equations above with \eqref{3p1} we obtain 14 equations, of which we are only interested in
\begin{align}
HM=0,\quad L+G=0 \quad M+H=1,\quad MG+HL=0,\quad HN+IM=0\label{eq001}\\
NG+JM+HP+IL=0,\quad QH+NI+MK=0 \quad JQ+PK=C.\label{eq002}
\end{align}
From  \eqref{eq001} and \eqref{eq002}, we have  $C=2\sigma_3(1-\operatorname{sgn}(y_0-y)\cdot g_0)=0, $ which is a contradiction.\par
On the other side, suppose that the equation \eqref{3p1} can be written as
\begin{align*}
\left(\overline{y}^3+G\overline{y}^2x+H\overline{y}x^2 + Ix^3+J\overline{y}^2+K\overline{y}x + Lx^2+My+Nx+P\right)(\overline{y} + Qx+R) = 0,
\end{align*}
Doing the multiplications on the left side, and comparing with \eqref{3p1}, we obtain 14 equations:
\begin{align}
IQ=0, \quad Q+G=0, \quad HQ+I=0, \quad GR+JQ+K=0 , \quad QG+H=1\label{eq003}\\R+J=-2,
\quad HR+KQ+L=0\quad LR+NQ=0, \label{eq004}\\IR+LQ=0, \quad JR+M=A,\quad                 	KR+MQ+N=B,\quad  \label{eq005}\\
\quad MR+P=C, \quad NR+PQ=0, \quad PR=D
\end{align}
Three cases can be considered: when $I = Q = 0,$ or $I=0,\;Q\neq0$ or $I\neq0,\;Q=0$.\medskip\\
\noindent\textbf{Case 1: $I=Q=0$.}\par
From \eqref{eq003} we have $ G=0, K=0, H=1, $ then, from \eqref{eq004} we have $ R+L=RL=0, $ from the last equation in \eqref{eq005} we conclude that $ D=0, $ which is absurd.\smallskip\\
\noindent \textbf{Case 2: $I=0$ and $Q\neq0$.}\par
From \eqref{eq003} we have $ H=0 $ and $ Q+G=0, QG=1, $ then there no exist $ Q $ and $ G $ real numbers, satisfying these conditions.\smallskip\\
\textbf{Caso 3: $I\neq0$ e $Q=0$.}\par
From the third equation in \eqref{eq003} we have $I=0,$ which is absurd.\medskip\\
Therefore, we conclude that it is impossible to write the right side of equation \eqref{3p1} as a product of lower degree polynomials.
\subsubsection{Examples - Funk's Parabolas type 3 and 4}
\begin{exa}\label{par3ex1}\hfill
\begin{figure}[H]
\centering
\caption{Funk's Parabolas type 3 with $\mathcal{F} = (0, 0)$ and $y_0 = 0.3$.}
\includegraphics[width=0.5\textwidth]{parabola_tipo_3_ex1.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's Parabolas type 3 with $\mathcal{F} = (0.6, 0)$ and $y_0 =  0.3$.}
\includegraphics[width=0.5\textwidth]{parabola_tipo_3_ex2.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's Parabolas type 3 with $\mathcal{F} = (0.7, 0.3)$ and $y_0 =  0.3$.}
\includegraphics[width=0.5\textwidth]{parabola_tipo_3_ex3.png}

\end{figure}
\begin{figure}[H]
\centering
\caption{Funk's Parabolas type 3 with $\mathcal{F} = (0, -0.4)$ and $y_0 =  -0.8$.}
\includegraphics[width=0.5\textwidth]{parabola_tipo_3_ex4.png}

\end{figure}
\end{exa}

\subsection*{Acknowledgements} We want to thank to Jonny Ardila for its help in the writing of this work.


\section{Introduction}

Over the years, Si-photonics has emerged as a promising CMOS-compatible material platform for fabricating low-cost, scalable integrated components for on-chip applications \cite{soref2006past,thomson2016roadmap}. Owing to the high-index contrast, high Kerr coefficient, transparency over a broad wavelength range (telecom to mid-IR), and enhanced device performance, Silicon-on-insulator (SOI) technologies in Si-photonics has been widely accepted for the future generation of CMOS integrated circuits (ICs). Researchers have  harnessed Si-photonics and employed in a plethora of applications, including high-speed data processing, sensing \cite{rogers2021universal,dorfner2008silicon}, nonlinear and quantum photonics \cite{harris2016large,leuthold2010nonlinear}, WDM systems \cite{fang2010wdm,moscoso20188}, and all-optical signal processing \cite{koos2007nonlinear}. However, despite the tremendous achievement of Si photonics, the small electronic band gap of Si (1.12 eV) \cite{leuthold2010nonlinear} imposes a fundamental limitation owing to the large two-photon (TPA) and free-carrier absorption (FCA) losses below 2.2 \(\mu\)m wavelength \cite{yin2007soliton}. These nonlinear losses are detrimental to the optical performance even at low power and prevent the widespread adoption of the Si platform in the telecommunication band \cite{yin2007impact,tsang2002optical,haldar2019free}. On the other hand, Si exhibits high thermo-optic coefficient (TOC) \((dn/dt = 1.86\) $\times$ \(10^{-4} K^{-1})\) which causes the SOI devices strongly thermally sensitive \cite{dong2010low}. The optical properties of Si-based devices can be significantly affected either by local heating causes due to the launching of high-power sources or changes in the environmental temperature. Therefore, Si based photonic devices are facing limitations in the application corresponding to high-temperature variation. Different approaches have been reported to reduce the thermal effect of Si-based devices, such as the incorporation of external metal heater for optimal heating \cite{manipatruni2008wide}, deposition of a cladding layer of Si with a material with negative TOC to compensate the positive TOC of Si \cite{han2007temperature,qiu2016athermal}, and so on. The limitation of the Si platform leads to exploring a new flexible CMOS-compatible platform carrying a refractive index lying between Si and \(SiO_{2}\), which can be realized for a multitude of photonics applications. In the quest for proper substitution of the Si platform, researchers are harnessing different material platforms to overcome the limitations associated with Si.

Over the years, stoichiometric silicon nitride $(Si_{3}N_{4})$ material has gained significant attention as a promising CMOS platform for building photonic integrated circuits (PICs) \cite{moss2013new} . Apart from its CMOS compatibility, it offers multiple advantages, such as a very large band gap with low propagation loss, high-refractive index contrast with silicon oxide $(SiO_{2})$; absence of TPA in the telecom wavelength range; and it offers relatively low TOC \((dn/dt = 2.45 \times 10^{-5} K^{-1})\) compared to both crystalline and amorphous Si \cite{arbabi2013measurements}. These properties of $(Si_{3}N_{4})$ led to employ it in a multitude of applications ranging from nonlinear photonics \cite{pfeiffer2017octave}, high Q resonator \cite{luke2015broadband}, and all-optical signal processing \cite{sharma2020review}. However, despite many conveniences, $Si_{3}N_{4}$ suffers from potential drawbacks, mainly due to the low refractive index contrast between core and cladding layers in the SOI platform, which in turn increases the device footprint and also low Kerr nonlinearity.

However, these bottlenecks have been effectively addressed by tuning the optical properties of SiN through proper control over the N/Si ratio of the material \cite{tan2018nonlinear}. Si-rich SiN (SRSN) material has attracted considerable interest due to its fabrication flexibility and the ability to tailor the intrinsic properties between those of $Si_{3}N_{4}$ and Si \cite{sohn2019optical}. By adjusting the ratio of the precursor gases (\(SiH_{4}/N_{2}\)) in the plasma-enhanced chemical vapor deposition (PECVD) process, the percentage of Si content in the material can be effectively tuned which provides effective route to engineer the optical properties of SRSN material \cite{tan2021nonlinear}. Proper N/Si ratio exhibits optical band gap variation from 2.7 to 5 eV and refractive index (n) varies over a wide range from 1.91 to 3.1 at 1550 nm wavelength \cite{nejadriahi2020thermo}. On the other hand, the third-order nonlinear coefficient increases effectively with the increase in Si content in the material \cite{kruckel2015linear}, which has already been demonstrated in a variety of applications such as parametric gain, broadband supercontinuum generation \cite{liu2016octave}, nonlinear signal processing \cite{dizaji2017silicon, lacava2017si}, and wavelength conversion \cite{yang2018octave,lacava2019intermodal}. The TOC as well as linear and nonlinear refractive index of SRSN material also increases linearly with the increase in the fractional composition of silicon over a range from that of silicon nitride to a-Si \cite{pruiti2020thermo,nejadriahi2020thermo}. Hence, careful control over the Si composition in the material is required to realize a thermally stable and high refractive index contrast integrated CMOS SRN platform.

\begin{figure}[t]
\centering
{\includegraphics[width=1\linewidth,height=4.5 cm]{refractive_index_f.eps}}
\caption{Refractive index as measured by ellipsometry for the SRN film.}
\label{fig:refractive_index_label}
\end{figure}

\begin{figure}[t]
\centering
   \begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1\linewidth ,height=4 cm]{2d_diagram}}
   \caption{}
   \label{fig:2d_diagram_label}
\end{subfigure}
\begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=0.9\linewidth ,height=3.5 cm]{modal_profile}}
   \caption{}
   \label{fig:mode_profile_label}
\end{subfigure}
\begin{subfigure}[]{0.4\textwidth}
\centering
 {\includegraphics[width=1\linewidth ,height=4 cm]{dispersion_SRN.eps}}
   \caption{}
   \label{fig:dispersion_label}
\end{subfigure}

\caption{(a) Schematic of the transverse view of the SRN waveguide (b) Modal profile of the fundamental TE mode of the proposed waveguide (c) Dispersion plot of the SRN waveguide structure. }
\end{figure}

In this manuscript, we demonstrate a broadband dispersion engineered racetrack ring resonator structure fabricated with the SRN material platform to realize on-chip applications in the telecommunication band. We designed and fabricated grating couplers at the input and output terminals to couple the light into the device and perform optical characterizations. We measure the TOC of PECVD deposited SRSN that exhibits a refractive index of 2.25 at 1550 nm. The TOC of the SRSN waveguide has been measured by determining the temperature-dependent resonance shift of the racetrack ring resonator structure in the transmission spectrum. Furthermore, we investigate power-dependent transmission spectra to characterize the thermal stability and nonlinear loss in the communication wavelength range.
\section{Design and fabrication of SRSN Waveguide}

\begin{figure}[t]
\centering
{\includegraphics[width=1\linewidth,height=4.5 cm]{coupling_efficiency_f.eps}}
\caption{Coupling efficiency of the grating coupler simulated through Lumerical FDTD software.}
\label{fig:coupling_efficiency}
\end{figure}

\begin{figure}[t]
\centering
\begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1.1\linewidth ,height=3.5 cm]{racetrack_ring_2d}}
   \caption{}
   \label{fig:racetrack_ring_label}
\end{subfigure}
   \begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1.1\linewidth ,height=3.5 cm]{racetrack_ring_3d_fff.eps}}
   \caption{}
   \label{fig:racetrack_ring_3d_label}
\end{subfigure}
\begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.6 cm]{racetrack_ring_microscope}}
   \caption{}
   \label{fig:racetrack_ring_microscope_label}
\end{subfigure}
\caption{(a) Schematic diagram of the racetrack ring resonator structure. (b) Perspective view of the racetrack ring resonator and the bus waveguide. (c) Optical microscopy image of the racetrack ring resonator along with the bus waveguide and the grating coupler.}
\end{figure}

For fabrication, we have employed the PECVD deposition technique to deposit SRSN on the top of the Si wafer with a 3 \(\mu\)m buried oxide layer. The deposition was done at a low temperature (\(350^{0}\) C) and comprises precursor gases like \(NH_{3}, SiH_{4}\) , and \(N_{2}\) with different flow rates. To characterize the real and imaginary parts of the refractive index of SRSN film, we deposit 100 nm of SRSN on the top of a clean Si wafer. After that, we carried out ellipsometry measurements with a broad range of wavelengths. Fig. \ref{fig:refractive_index_label} shows the dispersion of the real and imaginary parts of the refractive index of the deposited SRSN layer. It is observed that the value of the real part of the refractive index is significantly higher than that of silicon nitride material. We have employed this refractive index profile to design our SRSN-based racetrack ring resonator structure.

The schematic of the waveguide structure with optimized waveguide parameters is shown in Fig. \ref{fig:2d_diagram_label}. The waveguide parameters have been optimized using full-vectorial mode solver COMSOL software to achieve a flat dispersion profile in the telecom wavelength range.  The optimized height (H) and width (W) for the SRSN waveguide are 450 nm and 1000 nm, respectively. Fig. \ref{fig:mode_profile_label} shows the modal confinement of the fundamental TE mode at 1550 nm for the optimized waveguide parameters. Whereas the dispersion profile is shown in Fig. \ref{fig:dispersion_label}, which exhibits flat dispersion in the telecom band and can be useful for complex nonlinear applications.

\begin{figure*}[htp]
\centering
{\includegraphics[width=0.8\linewidth,height=8 cm]{coupling_efficiency_fff_f.eps}}
\caption{Schematic of experimental setup for characterization, SLD: Superluminescent Diode source, PC: polarization controller, DUT: Device under test, OSA: Optical spectrum analyzer.}
\label{fig:exp_setup}
 \end{figure*}

The optimized waveguide parameters have been employed for fabricating the racetrack ring resonator structure. The schematic of the racetrack ring resonator structure is shown in Fig. \ref{fig:racetrack_ring_label}. The racetrack resonator is designed for operating in the telecommunication wavelength range. The gap and coupling length have been chosen as 200 nm and 550 \(\mu\)m. To achieve maximum coupling, the grating coupler has been designed with an SRSN width of 285 nm, a shallow-etched depth of 250 nm, and a grating period of 870 nm. The coupling efficiency of the optimized design is shown in Fig. \ref{fig:coupling_efficiency}, which is simulated through Lumerical FDTD software. The 3-dB bandwidth of the grating coupler is calculated as 50 nm. To fabricate the device, we deposit a 3 \(\mu\)m \(SiO_{2}\) buried oxide layer on the top of a clean Si wafer by the PECVD method. After that, a 450 nm thick SRSN layer is deposited at low temperature using the PECVD method. After depositing the layers, the racetrack ring and the grating coupler are patterned employing electron-beam lithography (Raith-Eline) with a negative MaN 2401 resist, followed by reactive ion etching (RIE), with fluorine chemistry. Fig. \ref{fig:racetrack_ring_microscope_label} shows the microscopic view of the fabricated device.
\begin{figure}[htp]
\centering
{\includegraphics[width=1\linewidth,height=6 cm]{transmission_spectrum_f.eps}}
\caption{The measured transmission spectrum of the racetrack ring resonator at the through port for the TE mode.}
\label{fig:transmission}
\end{figure}

Device characterization is performed using a TEC-enabled broadband superluminescent compact laser diode (Thorlabs 1550S-A2), which is connected to the input probe with a single-mode fiber passing through a polarization controller and three-axis translation stage. The fiber incident angle is set at \(10 ^{0}\). The output signal was collected with an optical spectrum analyzer (Yokogawa
AQ6370D) ranging from 600 to 1700 nm. The schematic of the experimental setup is shown in Fig. \ref{fig:exp_setup}. The photograph of the actual coupling stage is shown on the inset of the diagram. The normalized TE transmission is shown in Fig. \ref{fig:transmission}. The measured spectrum shows periodic resonance dips with uniform spacing between the adjacent dips. The quality factor (Q) measured at 1553.797 nm is \(22.5 \times 10^{3}\), whereas the FWHM linewidth and free spectral range (FSR) are measured as 0.069 nm and 0.326 nm, respectively.

\section{Thermo-optic characterization}

To characterize the thermal stability of the device, we have performed TOC measurements of the SRSN platform by calculating the wavelength shift \(\Delta \lambda\) of the resonant peak with the change in the environmental temperature. The wavelength shift can be calculated as
\begin{equation}
\Delta\lambda = \bigg(\dfrac{\partial n_{eff}}{\partial T_{R}} + n_{eff}. \alpha_{sub}\bigg) \dfrac{\lambda}{n_{g}} T_{R}
 \label{eq1}
\end{equation}
where \(\lambda\) is the resonance wavelength at room temperature, \(\alpha_{sub}\) is the thermal expansion coefficient of the substrate, and \(\Delta T_{R}\) is the change in environmental temperature. \(n_{g}\) is the group index which can be expressed as \(n_{g} = \dfrac{\lambda^{2}}{L_{R}.FSR}\), \(L_{R}\) and FSR are the resonator length and free-spectral range, respectively. The effective TOC defines the change in the effective index as a function of temperature, which can be expressed as
\begin{equation}
\kappa_{eff} = \dfrac{\partial n_{eff}}{\partial T_{R}} = \bigg(\dfrac{n_{g}}{\lambda} \dfrac{\partial\lambda}{\partial T_{R}} - n_{eff}.\alpha_{sub} \bigg)
 \label{eq2}
\end{equation}
where \(\dfrac{\partial\lambda}{\partial T_{R}}\) defines the wavelength shift of the resonance dip as a function of temperature.
However, as the effective TOC involves the modal overlapping between the core of the SRSN waveguide and cladding oxide, eq. \ref{eq2} can be rewritten considering the overlapping factor as
\begin{equation}
\kappa_{eff} =  \Gamma_{SiN_{x}} \bigg(\dfrac{\partial n}{\partial T} \bigg)_{SiN_{x}} + \Gamma_{SiO_{2}} \bigg(\dfrac{\partial n}{\partial T} \bigg)_{SiO_{2}}
 \label{eq3}
\end{equation}
where \(\Gamma_{SiN_{x}}\) and \(\Gamma_{SiO_{2}}\) are the modal overlapping factor with SRSN waveguide core and cladding oxide, respectively.
\begin{figure}[htp]
\centering
   \begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.5 cm]{TOC_shift.eps}}
   \caption{}
   \label{fig:TOC_shift_label}
\end{subfigure}
\begin{subfigure}[]{0.4\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.1 cm]{TOC_shift_temp.eps}}
   \caption{}
   \label{fig:TOC_shift_temp_label}
\end{subfigure}
\caption{(a) Spectral response of the resonance dip at different temperature . (b) Plot of the resonance shift \(\Delta\lambda\) as a function of substrate temperature.}
\end{figure}

\begin{figure}[htp]
\centering
   \begin{subfigure}[]{0.48\textwidth}
 {\includegraphics[width=1\linewidth ,height=4 cm]{TPA_528nm.eps}}
   \caption{}
   \label{fig:TPA_528_label}
\end{subfigure}
\begin{subfigure}[]{0.48\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.7 cm]{temp_shift_528nm_f.eps}}
   \caption{}
   \label{fig:temp_shift_528nm_label}
\end{subfigure}
\caption{(a) Power-dependent transmission spectra for diode pump laser at 528 nm (b) corresponding spectral shift of resonance wavelength with the pump power variation.}
\end{figure}

\begin{figure}[htp]
\centering
   \begin{subfigure}[]{0.45\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.8 cm]{TPA_780nm.eps}}
   \caption{}
   \label{fig:TPA_780_label}
\end{subfigure}
\begin{subfigure}[]{0.45\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.5 cm]{temp_shift_780nm_f.eps}}
   \caption{}
   \label{fig:temp_shift_780nm_label}
\end{subfigure}
\caption{(a) Power-dependent transmission spectra for diode pump laser at 780 nm (b) corresponding spectral shift of resonance wavelength with the pump power variation.}
\end{figure}

\begin{figure}[htp]
\centering
   \begin{subfigure}[]{0.45\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.8 cm]{TPA_1310nm.eps}}
   \caption{}
   \label{fig:TPA_1310_label}
\end{subfigure}
\begin{subfigure}[]{0.45\textwidth}
 {\includegraphics[width=1\linewidth ,height=3.8 cm]{TPA_1550nm.eps}}
   \caption{}
   \label{fig:TPA_1550_label}
\end{subfigure}
\caption{Power-dependent transmission spectra for diode pump laser at (a) 1310 nm (b) 1550 nm wavelength.}
\end{figure}

To extract the TOC of the SRSN waveguide, we performed a temperature-dependent measurement of the racetrack ring resonator. The temperature-dependent transmission spectrum for the resonator structure is shown in Fig. \ref{fig:TOC_shift_label}, where the substrate temperature was varied across the range between \(40^{0}\)C to \(50^{0}\)C. The resonance dip at \(\sim\) 1554.32 nm exhibits red-shifting with a value of \(\sim\) 200 pm when the substrate temperature is increased from \(40^{0}\)C to \(50^{0}\)C. The \({\partial\lambda}/{\partial T_{R}}\) value has been calculated by performing linear fitting of the resonator wavelength shift against the substrate temperature as shown in Fig. \ref{fig:TOC_shift_temp_label} and the value obtained as 0.021 nm/\(^{0}C\). Taking \(L_{R}\) and FSR values as 2984 \(\mu\)m and 0.326 nm, respectively, the group index \(n_{g}\) value is calculated as 2.4829.  The calculated values of the overlapping factor for \(\Gamma_{SiN_{x}}\), \(\Gamma_{SiO_{2}}\) , effective mode indices \(n_{eff}\) are given by 0.846 and 0.096, and 1.854, respectively, which were obtained through the finite difference method. Furthermore, for the calculation, we consider \(\alpha_{sub}\) = \(2.6 \times 10^{-6}\). Putting all these values in eq. \ref{eq2}, we calculate the \(\kappa_{eff}\) as \(2.873\) $\times$ \(10^{-3}\) \(^{0} C^{-1}\). Considering \((\partial n/ \partial T)_{SiO_{2}}\) = 1 $\times$ \(10^{-5}\)  \(^{o} C^{-1}\), TOC of the SRSN waveguide is determined from eq. \ref{eq3} to be \(3.2825\) $\times$ \(10^{-5}\)\(^{o} C^{-1}\). It can be seen that the measured value of TOC is well below that of Si.

\section{TPA characterization}

Additionally, in order to characterize the thermal stability and nonlinear losses in the telecommunication band, we measure the power-dependent transmission spectrum. Optical characterization was performed using a tunable laser source coupled with a diode laser source at 528 nm, 780 nm, 1310 nm, and 1550 nm wavelengths. Finally, the output signal was collected by OSA. The transmission spectrum is measured as a function of input optical power. The power-dependent transmission spectrum for the diode laser pumped at 528 nm is shown in Fig. \ref{fig:TPA_528_label}, where a blue-shifted spectral response with an increase in input pump power is observed. The corresponding spectral shift with the incident pump power is shown in Fig. \ref{fig:temp_shift_528nm_label}. The fitting line shows a spectral shift of 3 pm/mW. It is important to note that the direction of the wavelength shifting depends on the dominant nonlinear process. The incident light gets absorbed via TPA, which in turn generates free carriers. The generated free carriers cause blue-shifting of the resonance peaks owing to free carrier dispersion (FCD). On the other hand, the generated free carriers trigger the free carrier absorption (FCA) loss. The combined effect of TPA, FCA, and the linear surface absorption leads to heating of the resonator, where the thermo-optic effect causes red-shifting of the resonance peak \cite{luo2012power,gao2017cavity}. As Si has a high positive TOC, thermal effects dominate and exhibit red-shifting in Si-based devices \cite{jeyaselvan2021mitigation, xiang2020effects}. In contrast, in our case, we observe a blue-shifting of the resonance peak. This is because the SRSN has a much lower TOC as measured in the previous section compared to Si. Hence, the FCD mechanism dominates over the thermo-optic effect, resulting in blue-shifting of the spectral peaks. Fig. \ref{fig:TPA_780_label} exhibits power-dependent transmission spectra for the pump wavelength at 780 nm. A small spectral blue shift is observed at that wavelength. The variation of the resonance shift with the input pump power is shown in Fig. \ref{fig:temp_shift_780nm_label}. The fitting line of the resonance shift exhibits a spectral shift of 0.21 pm/ mW at the wavelength of 780 nm wavelength, which is much lower than that observed at 580 nm. This is because the TPA effect gets reduced at 780 nm wavelength compared to 528 nm. The power-dependent transmission spectrum for a diode laser pumped at 1310 nm is shown in Fig. \ref{fig:TPA_1310_label}. The input power of the diode laser varied from 1.5 mW to 20 mW. It is observed that the SRSN-based device exhibits efficient thermal stability and suffers from no nonlinear losses as the resonance shape remains almost unchanged. On the other hand, the power dependent-transmission spectrum for diode laser at 1550 nm is shown in Fig. \ref{fig:TPA_1550_label}, which also exhibits remarkable thermal stability and the absence of nonlinear losses. Therefore, our results indicate that the SRSN device offers extensive thermal stability and the nonexistence of nonlinear losses in the telecommunication band.

\section{Conclusion}

In conclusion, we have demonstrated an SRSN-based broadband integrated racetrack ring resonator for on-chip applications in the telecommunication wavelength range. SRSN material exhibits an enhanced refractive index of 2.25 at 1550 nm. The TOC of the SRSN platform is characterized by measuring the temperature-dependent transmission spectra for the resonator structure, and the TOC is estimated to be \(3.2825\) $\times$ \(10^{-5}\)\(^{o}\) \(C^{-1}\). Furthermore, we have measured the power-dependent transmission spectra in the communication band. Transmission spectra almost remain unchanged for high input power, which reveals significant thermal stability, and the absence of nonlinear absorption. Our experimental findings pave the way toward the implementation of on-chip devices for linear and nonlinear applications in the telecommunication band.
\section*{Acknowledgments}

DST/Nano Mission and Ministry of Electronics and Information Technology.


\section{Introduction}
\label{sec:introduction}

Rapid development in \ac{CLT} has been achieved by pre-training transformer-based models in large multilingual corpora \cite{conneau_unsupervised_2020, xue-etal-2021-mt5}, where these models have state-of-the-art results in multilingual NLU benchmarks \cite{ruder-etal-2021-xtreme}.
Moreover, adapter-based fine-tuning \cite{houlsby2019parameterefficient, pfeiffer-etal-2020-mad} has been proposed to minimize the misalignment of multilingual knowledge (alignment) when \ac{CLT} is applied, especially in a zero-shot fashion, where the target language is unseen during training.
\ac{CLT} is severely understudied in legal \acs{NLP} applications except for \citet{chalkidis-etal-2021-multieurlex} who experimented with several methods for \ac{CLT} on MultiEURLEX, a newly introduced multilingual legal topic classification dataset, including EU laws.

To the best of our knowledge, \ac{CLT} has not been applied to the \ac{LJP} task \cite{aletras_predicting_2016,xiao_cail2018_2018, chalkidis-etal-2019-neural,malik-etal-2021-ildc}, where the goal is to predict the verdict (court decision)  given  the  facts  of  a  legal  case.
In this setting, positive impact of cross-lingual transfer is not as conceptually straight-forward as in other general applications (NLU), since there are known complications for sharing legal definitions and interpreting law across languages \cite{gotti-2014-translate,mcauliffe-2014-translate,robertson2016multilingual,ramos-2021-translate}.

\begin{figure}[t]
    \centering
    \resizebox{0.97\columnwidth}{!}{
    \includegraphics{progress.png}
    }
    \vspace{-2mm}
    \caption{
    Incremental performance improvement through several development steps.

    }.

    \label{fig:progress}
    \vspace{-10mm}
\end{figure}

Following the work of \citet{niklaus-etal-2021-swiss}, we experiment with their newly released trilingual \ac{SJP} dataset, containing cases from the \ac{FSCS}, written in three official Swiss languages (German, French, Italian).
The dataset covers four legal areas (public, penal, civil, and social law) and lower courts located in eight regions of Switzerland (Zurich, Ticino, etc.), which poses interesting new challenges on model robustness / fairness and the effect of cross-domain and cross-regional knowledge sharing. In their experiments, \citet{niklaus-etal-2021-swiss} find that the performance in cases written in Italian is much lower compared to the rest, while also performance varies a lot across regions and legal areas.

\subsection*{Main Research Questions}
We pose and examine four main research questions:

\noindent\textbf{RQ1}: \emph{Is cross-lingual transfer beneficial across all or some of the languages?}

\noindent\textbf{RQ2}: \emph{Do models benefit or not from cross-regional and cross-domain transfer?}

\noindent\textbf{RQ3}: \emph{Can we leverage data from another jurisdiction to improve performance?}

\noindent\textbf{RQ4}: \emph{How does representational bias (wrt. language, origin region, legal area) affect model's performance?}

\subsection*{Contributions}
The contributions of this paper are fourfold:
\vspace{-2mm}

\begin{itemize}[leftmargin=8pt, itemsep=0em]

    \item We explore, for the first time, the application of cross-lingual transfer learning in the challenging \ac{LJP} task in several settings (Section~\ref{sec:cross_lingual_exps}).

    We find that a pre-trained language model fine-tuned multilingually, outperforms its monolingual counterparts, especially when we use adapter-based fine-tuning and augment the training data with machine-translated versions of the original documents ($3\times$ larger training corpus) with larger gains in a low-resource setting (Italian).

    \item  We perform cross-domain and cross-regional analyses (Section~\ref{sec:cross_domain_exps}) exploring the effects of cross-domain and cross-regional transfer, i.e., train a model across domains, i.e., legal areas (e.g., civil, penal law), or regions (e.g., Zurich, Ticino). We find that in both settings (legal areas, regions), models trained across all groups perform overall better and more robustly; while always improving performance in the worst-case (region or legal area) scenario.

    \item We also report improved results when we apply cross-jurisdiction transfer (Section~\ref{sec:cross_jurisdiction_exps}) , where we further augment our dataset with Indian legal cases originally written in English.
    \item We release the augmented dataset (incl. 100K machine-translated documents) and our code for replicability and future experimentation.\footnote{\url{https://huggingface.co/datasets/swiss_judgment_prediction}}

\end{itemize}

The cumulative performance improvement amounts to 7\% overall and 16+\% in the low-resource Italian subset, compared to the best reported scores in \citet{niklaus-etal-2021-swiss}, while using cross-lingual and cross-jurisdiction transfer we improve for 2.3\% overall and 4.6\% for Italian over our strongest baseline (NativeBERTs).

\section{Dataset and Task description}
\subsection{Swiss Legal Judgment Prediction Dataset}
\label{sec:swiss_ljp}

We investigate the \ac{LJP} task on the Swiss-Judgment-Prediction (SJP) dataset \cite{niklaus-etal-2021-swiss}. The dataset contains 85K cases from the Federal Supreme Court of Switzerland (FSCS) from the years 2000 to 2020 written in German, French, and Italian. The court hears appeals focusing on small parts of the previous (lower court) decision, where they consider possible wrong reasoning by the lower court. The dataset provides labels for a simplified binary (\emph{approval}, \emph{dismissal}) classification task. Given the facts of the case, the goal is to predict if the plaintiff's request is valid or partially valid (i.e., the court \emph{approved} the complaint).

Since the dataset contains rich metadata, such as legal areas and origin regions, we can conduct experiments on the robustness of the models (see Section \ref{sec:cross_domain_exps}). The dataset is not equally distributed; in fact, there is a notable representation disparity where Italian have far fewer documents (4K), compared to German (50K) and French (31K). Representation disparity is also vibrant with respect to legal areas and regions. We refer readers to the work of \citeauthor{niklaus-etal-2021-swiss} for detailed dataset statistics.

\subsection{Indian Legal Judgment Prediction Dataset}
\label{sec:indian_ljp}

The \ac{ILDC} dataset \cite{malik-etal-2021-ildc} comprises 30K cases from the Indian Supreme Court in English. The court hears appeals that usually include multiple petitions and rules a decision (\emph{accepted} vs. \emph{rejected}) per petition. Similarly to \citet{niklaus-etal-2021-swiss}, \citeauthor{malik-etal-2021-ildc} released a simplified version of the dataset
with binarized labels. In effect, the two datasets (SJP, ILDC) target the very same task (partial or full approval of plaintiff's claims), nonetheless in two different jurisdictions (Swiss Federation and India). Our main goal, when we use ILDC as a complement of SJP, is to assess the possibility of cross-jurisdiction transfer from Indian to Swiss cases (see Section~\ref{sec:cross_jurisdiction_exps}), an experimental scenario that has not been explored so far in the literature.

\subsection{NMT-based Data Augmentation}
\label{sec:easy_nmt}

In some of our experiments, we perform data augmentation using machine-translated versions of the original documents, i.e., translate a document originally written in a single language to the other two (e.g., from German to French and Italian). We performed the translations using the EasyNMT\footnote{\url{https://github.com/UKPLab/EasyNMT}} framework utilizing the \emph{many-to-many} \ac{NMT}  model
of \citet{fan2020}.\footnote{The \emph{one-to-one} OPUS-MT \cite{tiedemann-thottingal-2020-opus} models did not have any model available from French to Italian (fr2it) at the time of the experiments.} A preliminary manual check of some translated samples showed sufficient translation quality to proceed forward. We release the machine-translated additional dataset for future consideration on cross-lingual experiments or quality assessment.

To the best of our knowledge, machine translation for data augmentation has not been studied in legal \ac{NLP}  applications, while it is generally a straight-forward, though under-studied idea. As we show in the experiments (see Section \ref{sec:cross_lingual_exps}), the translations are effective, leading to an average improvement of 1.6\% macro-F1 for standard fine-tuning and 0.8\% for adapter-based one (see Table \ref{tab:cross_lingual}). For the low-resource Italian subset, the improvement even amounts to 3.2\% and 1.6\%, respectively.

\begin{table*}[t]
\centering
    \resizebox{\textwidth}{!}{
    \begin{tabular}{lrc|ccc|cc}
    \toprule
    \bf Model & \bf \#D & \bf \#M   &             \bf German $\uparrow$ &  \bf French $\uparrow$ &   \bf Italian $\uparrow$ &  \bf All $\uparrow$ & (Diff. $\downarrow$)\\
    \midrule
    \multicolumn{8}{l}{A1. \emph{Monolingual: Fine-tune  on the \bf{tgt training set}} (\textbf{src} $=$ \textbf{tgt}) --- Baselines }  \\
    \midrule

    Prior SotA (\small \citeauthor{niklaus-etal-2021-swiss}) & 3-35K & N & 68.5 {\small ± 1.6} & 70.2 {\small ± 1.1} & 57.1 {\small ± 0.4} & 65.2 {\small ± 0.8} & ( 13.1 )\\
    \midrule
    NativeBERTs  & 3-35K & N & \underline{69.6} {\small ± 0.4} &  \underline{72.0} {\small ± 0.5} &   \underline{68.2} {\small ± 1.3} &  \underline{69.9} {\small ± 1.6} & ( 3.8 )\\
    XLM-R        & 3-35K & N &  68.2 {\small ± 0.3} &  69.9 {\small ± 1.6} &  65.9 {\small ± 1.2} &  68.0 {\small ± 2.0} & ( 4.0 )\\
    \midrule
    \multicolumn{8}{l}{A2. \emph{Monolingual: Fine-tune on the \bf{tgt training set incl. machine-translations}} (\textbf{src} $=$ \textbf{tgt}) } \\
    \midrule
     NativeBERTs   & 60K & N &  \underline{70.0} {\small ± 0.7} &  \underline{71.0} {\small ± 1.3} &   \underline{71.9} {\small ± 2.5} &  \underline{71.0} {\small ± 0.8} & ( 0.9 ) \\
    XLM-R       & 60K & N &  68.8 {\small ± 1.4} &  70.7 {\small ± 2.1} &   71.9 {\small ± 2.6} &  70.4 {\small ± 1.3} & ( 1.1 ) \\
    \midrule
    \multicolumn{8}{l}{B1. \emph{Cross-lingual: Fine-tune on \bf{all training sets}} (\textbf{src} $\subset$ \textbf{tgt})} \\
    \midrule
    XLM-R   & 60K & 1 &  68.9 {\small ± 0.3} &  71.1 {\small ± 0.3} & 68.9 {\small ± 1.4} &  69.7 {\small ± 1.0} & ( 2.2 ) \\
    XLM-R {\small+ Adapters} & 60K & 1 &  \underline{69.9} {\small ± 0.6} &  \underline{71.8} {\small ± 0.7} &  \underline{70.7} {\small ± 1.8} &  \underline{70.8} {\small ± 0.8} & ( 0.9 ) \\
    \midrule
    \multicolumn{8}{l}{B2. \emph{Cross-lingual: Fine-tune on \bf{all training sets incl. machine-translations}} (\textbf{src} $\subset$ \textbf{tgt})} \\
    \midrule
     XLM-R   & 180K & 1 &  70.2 {\small ± 0.5} &  71.5 {\small ± 1.1} &  72.1 {\small ± 1.2} &  71.3 {\small ± 0.7} & ( 1.9 ) \\
    XLM-R {\small+ Adapters} & 180K & 1 &  \underline{\textbf{70.3}} {\small ± 0.9} &  \underline{\textbf{72.1}} {\small ± 0.8} &  \underline{\textbf{72.3}} {\small ± 2.1} &  \underline{\textbf{71.6}} {\small ± 0.8} & ( 2.0 ) \\
    \midrule
    \multicolumn{8}{l}{C. \emph{Zero-shot Cross-lingual: Fine-tune on \bf{all training sets excl. tgt language}} (\textbf{src} $\neq$ \textbf{tgt})} \\
    \midrule    XLM-R   & 25-57K & 1 &  58.4 {\small ± 1.2} &  58.7 {\small ± 0.8} &  \underline{68.1} {\small ± 0.2} &  61.7 {\small ± 4.5} & ( 9.7 ) \\
    XLM-R {\small+ Adapters} & 25-57K & 1 &  \underline{62.5} {\small ± 0.6} &  \underline{58.8} {\small ± 1.5} &  67.5 {\small ± 2.2} &  \underline{62.8} {\small ± 3.7} & ( 8.7 ) \\
    \bottomrule
    \end{tabular}
     }
    \caption{Test results for all training set-ups (monolingual w/ or w/o translations, multilingual w/ or w/o translations, and zero-shot) w.r.t source (src) and target (tgt) language. Best overall results are in \textbf{bold}, and best per setting (group) are \underline{underlined}.
    \#D is the number of training documents used. \#M is the number of models trained/used.
    The mean and standard deviation are computed across random seeds and across languages for the last column. Diff. shows the difference between the best and the worst performing language. \textbf{\emph{The adapter-based multilingually fine-tuned XLM-R model including machine-translated versions  ($3\times$ larger corpus) has the best overall results}}.}
    \label{tab:cross_lingual}
    \vspace{-4mm}
\end{table*}

\section{Experiments}
\label{sec:experiments}

\subsection{Hierarchical BERT}
\label{sec:hier_bert}
Since the examined dataset (\ac{SJP}) contains many documents with more than 512 tokens (90\% of the documents are up to 2048), we use Hierarchical BERT models \cite{chalkidis-etal-2019-neural,niklaus-etal-2021-swiss, dai-etal-2022-revisit} to encode up to 2048 tokens per document (4$\times$512 blocks).

We split the text into consecutive blocks of 512 tokens and feed the first 4 blocks to a shared standard BERT encoder. Then, we aggregate the block-wise \texttt{CLS} tokens by passing them through another 2-layer transformer encoder, followed by max-pooling and a final classification layer.

We re-use and expand the implementation released by \citet{niklaus-etal-2021-swiss},\footnote{\url{https://github.com/JoelNiklaus/SwissJudgementPrediction}} which is based on the Hugging Face library \cite{wolf_transformers_2020}. Notably, we first improve the masking of the blocks. Specifically, when the document has less than the maximum number (4) of blocks, we pad with extra sequences of \texttt{PAD} tokens, without the use of special tokens (\texttt{CLS}, \texttt{SEP}), as was previously performed. This minor technical improvement seems to affect the model's performance at large (group A1 Prior SotA vs. NativeBERTs –– Table~\ref{tab:cross_lingual}).

We experiment with monolingually pre-trained BERT models (aka NativeBERTs) and the multilingually pre-trained XLM-R
of \citet{conneau_unsupervised_2020}.
Specifically, for monolingual experiments (Native BERTs), we use German-BERT \cite{chan_deepset_2019} for German,
CamemBERT \cite{martin_camembert_2020} for French,
and UmBERTo \cite{parisi_umberto_2020} for Italian,
similar to \citet{niklaus-etal-2021-swiss}.

In our multilingual experiments, we also assess the effectiveness of adapter-based fine-tuning \cite{houlsby2019parameterefficient, pfeiffer-etal-2020-mad}, in comparison to standard full fine-tuning. In this setting, adapter layers are placed after all feed-forward layers of XLM-R and are trained together with the parameters of the layer-normalization layers. The rest of the model parameters remain untouched.

\subsection{Experimental Set Up}
\label{sec:experimental_setup}

We follow \citet{niklaus-etal-2021-swiss} and report macro-averaged F1 score to account for the high class-imbalance in the dataset (approx.\ 20/80 approval/dismissal ratio). We repeat each experiment with 3 different random seeds and report the average score and standard deviation across runs (seeds). We perform grid-search for the learning rate and report test results, selecting the hyper-parameters with the best development scores.\footnote{Additional details on model configuration, training, and hyper-parameter tuning can be found in Appendix~\ref{sec:hyperparam_tuning}.}

\subsection{Cross-lingual Transfer}
\label{sec:cross_lingual_exps}

We first examine \emph{cross-lingual transfer}, where the goal is to share (transfer) knowledge across languages, and we compare models in three main settings:
(a) \emph{Monolingual} (see Section \ref{sec:monolingual}): fine-tuned per language, using either the documents originally written in the language, or an augmented training set including the machine-translated versions of all other documents (originally written in another language),
(b) \emph{Cross-lingual} (see Section \ref{sec:crosslingual}): fine-tuned across languages with or without the additional translated versions, and
(c) \emph{Zero-shot cross-lingual} (see Section \ref{sec:zeroshot_crosslingual}): fine-tuned across a subset of the languages excluding the target language at a time.
We present the results in Table~\ref{tab:cross_lingual}.

\subsubsection{Mono-Lingual Training}
\label{sec:monolingual}
We observe that the baseline of \emph{monolingually} pre-trained and fine-tuned models (NativeBERTs) have the best results compared to the \emph{multilingually} pre-trained but \emph{monolingually} fine-tuned XLM-R (group A1 -- Table~\ref{tab:cross_lingual}). Representational bias across languages (Section~\ref{sec:swiss_ljp}) seems to be a key part of performance disparity, considering the performance of the least represented language (Italian) compared to the rest (3K vs. 21-35K training documents). However, this is not generally applicable, i.e., French have better performance compared to German, despite having approx.\ 30\% less training documents.

Translating the full training set provides a $3\times$ larger training set (approx.\ 180K in total) that ``equally'' represents all three languages.\footnote{Representational equality with respect to number of training documents per language, but possibly not considering text quality, since we use \ac{NMT}  to achieve that goal.} Augmenting the original training sets with translated versions of the documents (group A2 -- Table~\ref{tab:cross_lingual}), originally written in another language, improves performance in almost all (5/6) cases (languages per model). Interestingly, the performance improvement in Italian, which has the least documents (less than 1/10 compared to German), is the largest across languages with 3.7\% for NativeBERT (68.2 to 71.9) and 6\% for XLM-R (65.9 to 71.9) making Italian the best performing language after augmentation. Data augmentation seems more beneficial for XLM-R, which does not equally represent the three examined languages.\footnote{Refer to \citet{conneau_unsupervised_2020} for resources per language  used to pre-train XLM-R (50\% less tokens for Italian).}

\begin{table*}[t]
    \centering
    \resizebox{\textwidth}{!}{
    \begin{tabular}{lrl|cccccccc|c}
\toprule
\bf Origin Region & \bf \#D & \bf \#L &  \bf ZH &  \bf ES &   \bf CS &  \bf NWS  &   \bf EM &    \bf RL &   \bf TI &     \bf FED &  \bf All \\
\midrule
\multicolumn{12}{c}{Region-specific fine-tuning with MT data augmentation} \\
\midrule
Zürich (ZH)                     & 26.4K  & de        & \underline{65.5} &  65.6 &  63.7 &      68.2 &  62.0 &  57.9 &  63.2 &  54.8 &  62.6 \\
Eastern Switzerland (ES)        & 17.1K  & de        & 62.9 &  \underline{66.9} &  62.8 &      65.2 &  62.2 &  60.2 &  57.8 &  55.1 &  61.6 \\
Central Switzerland (CS)        & 14.4K  & de        & 62.5 &  65.5 &  \underline{63.2} &      65.1 &  60.7 &  57.8 &  60.5 &  55.9 &  61.4 \\
Northwestern Switzerland (NWS)  & 17.1K  & de        & 66.0 &  68.6 &  65.2 &      \underline{67.9} &  61.6 &  57.0 &  57.1 &  55.5 &  62.4 \\
Espace Mittelland (EM)          & 24.9K  & de,fr     & 64.1 &  66.6 &  63.3 &      66.7 &  \underline{64.0} &  66.8 &  63.2 &  58.4 &  64.1 \\
Région Lémanique (RL)           & 40.2K & fr,de     & 61.0 &  64.7 &  60.2 &      63.7 &  63.4 &  \underline{69.8} &  67.6 &  54.3 &  63.1 \\
Ticino (TI)                     & 6.9K  & it        & 55.0 &  56.3 &  53.2 &      54.5 &  56.0 &  54.7 &  \underline{66.0} &  53.1 &  56.1 \\
Federation (FED)                & 3.9K  & de,fr,it  & 57.5 &  59.6 &  56.8 &      58.9 &  55.0 &  56.5 &  53.5 &  \underline{54.9} &  56.6 \\
\midrule
\multicolumn{12}{c}{Cross-regional fine-tuning w/o MT data augmentation} \\
\midrule
XLM-R                           & 60K & de,fr,it & 68.5 &  71.3  &  67.7  &  71.2  &  69.0 &  71.4  &  67.4  &  64.6 &  68.9  \\
XLM-R {\small+ Adapters}        & 60K & de,fr,it & \textbf{69.2} &  \textbf{73.9}  &  67.9 &      72.6  &  69.0  &  \textbf{72.1} &  70.1  & 64.2  &  69.9  \\
\midrule
\multicolumn{12}{c}{Cross-regional fine-tuning with MT data augmentation} \\
\midrule
NativeBERTs                     & 180K & de,fr,it &  69.0 &  72.1 &   68.6 &      72.0 &   69.9 &   71.9 &  68.8 &  64.8 &  69.6 \\
XLM-R                           & 180K & de,fr,it & \bf 69.2 &  72.9 &  68.3 &       \bf 73.3 &    69.9 &  71.7 &   70.4 &   \bf 65.0 &   70.1 \\
XLM-R {\small+ Adapters} & 180K & de,fr,it & \textbf{69.2}  &  73.3  &  \textbf{69.9}  &      73.0 &  \textbf{70.3} &  \textbf{72.1}  &  \textbf{70.9} &  63.8  &  \textbf{70.3} \\
\bottomrule
\end{tabular}
}
    \vspace{-2mm}
    \caption{Test results for models trained per region or across all regions. Best overall results are in \textbf{bold}, and in-domain are \underline{underlined}.
    \#D is the total number of training examples. \#L are the languages covered.
    \textbf{\emph{Cross-regional transfer is beneficial for all regions and has the best overall results. The shared multilingual model trained across all languages and regions slightly outperforms the baseline (NativeBERTs).}}
    }
    \label{tab:cross_domain_origin_regions}
    \vspace{-4mm}
\end{table*}

\subsubsection{Cross-Lingual Training}
\label{sec:crosslingual}
We now turn to the \emph{cross-lingual transfer} setting, where we train XLM-R across all languages in parallel. We observe that cross-lingual transfer (group B1 -- Table~\ref{tab:cross_lingual}) improves performance (+4.5\% p.p.) across languages compared to the same model (XLM-R) fine-tuned in a monolingual setting (group A1 -- Table~\ref{tab:cross_lingual}). This finding suggests that cross-lingual transfer (and the inherited benefit of using larger multilingual corpora) has a significant impact, despite the legal complication of sharing legal definitions across languages.
Augmenting the original training sets with the documents translated across all languages, further improves performance (group B2 -- Table~\ref{tab:cross_lingual}).

\subsubsection{Zero-Shot Cross-Lingual Training}
\label{sec:zeroshot_crosslingual}
We also present results in a \emph{zero-shot cross-lingual} setting (group C -- Table~\ref{tab:cross_lingual}), where XLM-R is trained in two languages and evaluated in the third one (unseen in fine-tuning). We observe that German has the worst performance (approx. 10\% drop), which can be justified as German is a \emph{Germanic} language, while both French and Italian are  \emph{Romance} and share a larger part of the vocabulary.

Contrarily, in case of Italian, the low-resource language in our experiments, the model strongly benefits from zero-shot cross-lingual transfer, leading to 2.2\% p.p. improvement, compared to the monolingually trained XLM-R. In other words, training XLM-R with much more (approx\ $20\times$) out-of-language (57K in German and French) data is better compared to training on the limited (3K) in-language (Italian) documents (68.1 vs. 65.9).

\subsubsection{Fine-tuning with Adapters}
\label{sec:adapters}
Across all cross-lingual settings (groups B-C -- Table~\ref{tab:cross_lingual}), the use of Adapters improves substantially the overall performance. The multilingual adapter-based XLM-R in group B1 (Table~\ref{tab:cross_lingual}) has comparable performance to the NativeBERTs models of group A2, where the training dataset has been artificially augmented with machine translations. In a similar setting (group B2 -- Table~\ref{tab:cross_lingual}),  the multilingual adapter-based XLM-R in group B2 has the best overall results, combining the benefits of both cross-lingual transfer and data augmentation.

With respect to \emph{cross-lingual performance parity}, the adapter-based XLM-R model has also the highest performance parity (least diff. in the last column of Table~\ref{tab:cross_lingual}), while augmenting the dataset with \ac{NMT}  translations leads to both the worst-case (language) performance and best performance for the least represented language (Italian).

In conclusion, cross-lingual transfer with an augmented dataset comprised of the original and machine-translated versions of all documents, has the best overall performance with a vibrant improvement (3\% compared to our strong baselines -- second part of Group A1 in Table~\ref{tab:cross_lingual}) in Italian, the least represented language.

\begin{table*}[ht]
\centering
    \resizebox{\textwidth}{!}{
    \begin{tabular}{lr|cccc|c}
    \toprule
    \bf Legal Area & \bf \#D &                 \bf Public Law &                   \bf Civil Law &                    \bf Penal Law &                  \bf Social Law &                      \bf All \\
    \midrule
    \multicolumn{7}{c}{Domain-specific fine-tuning with MT data augmentation} \\
    \midrule

    Public Law & 45.6K & \underline{56.4} {\small ± 2.2} &  52.2 {\small ± 2.0} &   59.7 {\small ± 4.9} &  60.1 {\small ± 5.8} &   57.1 {\small ± 3.2} \\
    Civil Law  & 34.5K & 44.4 {\small ± 7.9} &  \underline{64.2} {\small ± 0.6} &  45.5 {\small ± 13.1} &  43.6 {\small ± 5.2} &   49.4 {\small ± 8.6} \\
    Penal Law  & 35.4K & 40.8 {\small ± 10.1} &  55.8 {\small ± 2.9} &   \textbf{\underline{84.5}} {\small ± 1.3} &  61.1 {\small ± 7.5} &  60.6 {\small ± 15.7} \\
    Social Law & 29.1K  & 52.6 {\small ± 4.2} &  56.6 {\small ± 2.0} &   69.0 {\small ± 5.5} &  \underline{70.2} {\small ± 2.0} &   62.1 {\small ± 7.6} \\
    \midrule
    \multicolumn{7}{c}{Cross-domain fine-tuning w/o MT data augmentation} \\
    \midrule
    XLM-R & 60K & 57.4 {\small ± 2.0} &  66.1 {\small ± 3.1} &   81.4 {\small ± 1.4} &  70.8 {\small ± 2.0}  &   68.9 {\small ± 8.7} \\
     XLM-R {\small+ Adapters} & 60K &   58.4 {\small ± 2.5} &  66.1 {\small ± 2.4} &  83.1 {\small ± 1.2} &  71.1 {\small ± 1.4} &   69.7 {\small ± 9.0} \\
    \midrule
    \multicolumn{7}{c}{Cross-domain fine-tuning with MT data augmentation} \\
    \midrule
    NativeBERTs & 180K & 58.1 {\small ± 3.0} &  64.5 {\small ± 3.7} &   83.0 {\small ± 1.3} & 71.1 {\small ± 4.3} &   69.2 {\small ± 9.2} \\
     XLM-R         & 180K &   58.0 {\small ± 3.0} &  \textbf{67.2} {\small ± 1.6} &   84.4 {\small ± 0.2} &  70.2 {\small ± 1.3} &   \textbf{70.0} {\small ± 9.5} \\
     XLM-R {\small+ Adapters}         &  180K &   \textbf{58.6} {\small ± 2.7} &  66.8 {\small ± 2.8} &  83.1 {\small ± 1.3} &  \textbf{71.3} {\small ± 2.4} & 69.9 {\small ± 8.8} \\
    \bottomrule
    \end{tabular}
    }
    \caption{Test results for models (XLM-R with MT unless otherwise specified) \textbf{fine-tuned} per legal area (domain) or across all legal areas (domains). Best overall results are in \textbf{bold}, and in-domain are \underline{underlined}. The mean and standard deviations are computed across languages per legal area and across legal areas for the right-most column. \#D is the total number of training examples. \textbf{\emph{Cross-domain transfer is beneficial for 3 out of 4 legal areas and has the best overall results.}} The shared multilingual model trained across all languages and legal areas outperforms the baseline (monolingual BERT models).}
    \label{tab:cross_domain_legal_areas}
    \vspace{-4mm}
\end{table*}

\subsection{Cross-Domain/Regional Transfer Analysis}
\label{sec:cross_domain_exps}

Further on, we examine the benefits of transfer learning (knowledge sharing) in other dimensions. Hence, we analyze model performance with respect to origin regions and legal areas (domains of law).

\subsubsection{Origin Regions}
\label{sec:cross_domain_origin_regions}

In Table \ref{tab:cross_domain_origin_regions} we present the results for \emph{cross-regional} transfer.
In the top section of the table, we present results with region-specific multilingual (XLM-R) models evaluated across regions (in-region on the diagonal, zero-shot otherwise). We observe that the cross-regional models (two lower groups of Table~\ref{tab:cross_domain_origin_regions}) always outperform the region-specific models. Moreover, cross-lingual transfer is beneficial across cases, while adapter-based fine-tuning further improves results in 5 out of 8 cases (regions). Data augmentation is also beneficial in most cases.

In the top part of Table~\ref{tab:cross_domain_origin_regions}, in 60\% of the cases (regions: ZH, ES, CS, NWS, TI), a ``zero-shot'' model, i.e., trained in the cases of another region, slightly outperforms the in-region model.
In other words, in almost every case (target region), there is another \emph{monolingual} region-specific model that outperforms the in-region one.

We consider two main factors that may explain these results: (a) the region-wise \emph{representational bias} considering the number of cases per region, and (b) the cross-regional \emph{topical similarity} of the training and test subsets across different regions. To approximate the cross-regional topical similarity, we consider the distributional similarity (or dissimilarity) w.r.t.\ legal areas (Table~\ref{tab:legal_area_distribution_distances} in Appendix~\ref{sec:distances}). None of these factors can fully explain the results. Although in 3 out of 5 cases, the best performing (out-of-region) model has been trained on more data compared to the in-region one. There are also other confounding factors (e.g., language), i.e., models trained on the cases of either Espace Mittelland (EM) or Région Lémanique (RL), both bilingual with 8-10K cases, have the best results across all single-region models, hence a further exploration of the overall dynamics is needed.

\begin{table*}[ht]
\centering
    \resizebox{\textwidth}{!}{
    \begin{tabular}{llr|ccc|cc}
    \toprule
    \bf Model & \bf Training Dataset & \bf \#D   &   \bf German $\uparrow$   &  \bf French $\uparrow$ & \bf Italian $\uparrow$ &  \bf All & (Diff. $\downarrow$)\\
    \midrule
    \multicolumn{7}{c}{Cross-lingual fine-tuning w/ or w/o MT data augmentation} \\
    \midrule
     XLM-R   & Original & 60K &  68.9 {\small ± 0.3} &  71.1 {\small ± 0.3} & 68.9 {\small ± 1.4} &  69.7 {\small ± 1.0} & ( 2.2 )\\
    XLM-R {\small+ Adapters} & Original & 60K &  \underline{69.9} {\small ± 0.6} &  \underline{71.8} {\small ± 0.7} &  \underline{70.7} {\small ± 1.8} &  \underline{70.8} {\small ± 0.8}  & ( 0.9 ) \\
    \midrule
    XLM-R &  + MT Swiss     & 180K                & 70.2 {\small ± 0.5} &  71.5 {\small ± 1.1} & \underline{72.1} {\small ± 1.2} & 71.3 {\small ± 0.7}  & ( 1.9 ) \\
    XLM-R {\small+ Adapters} & + MT Swiss  & 180K  & \underline{70.3} {\small ± 0.8} &  \underline{72.1} {\small ± 0.8} & \underline{72.1} {\small ± 1.2} & \underline{71.5} {\small ± 0.9} & ( 1.8 ) \\
    \midrule
    \multicolumn{7}{c}{Cross-jurisdiction fine-tuning w/ MT data augmentation} \\
    \midrule
    XLM-R &  + MT \{Swiss, Indian\} & 276K               & 70.5 {\small ± 0.4} & 71.8 {\small ± 0.3} & \textbf{\underline{73.5}} {\small ± 1.4} & 72.0 {\small ± 0.9}  & ( 3.0 ) \\
    XLM-R {\small+ Adapters} & + MT \{Swiss, Indian\}   & 276K     & \textbf{\underline{71.0}} {\small ± 0.4}  & \textbf{\underline{73.0}} {\small ± 0.6} & 72.6 {\small ± 1.1} & \textbf{72.2} {\small ± 1.2}  & ( 2.0 )\\
    \midrule
    \multicolumn{7}{c}{Cross-jurisdiction zero-shot fine-tuning w/ MT data augmentation} \\
    \midrule
    XLM-R & MT Indian         & 96K             & 50.4 {\small ± 1.5} & 47.9 {\small ± 1.0} & 49.5 {\small ± 1.3} & 49.3 {\small ± 1.0}  & ( 2.5)  \\
    XLM-R {\small+ Adapters} & MT Indian    & 96K   & \underline{51.6} {\small ± 2.9} & \underline{49.7} {\small ± 1.4} & \underline{50.1} {\small ± 1.4} & \underline{50.5} {\small ± 1.0}  & ( 1.9 ) \\
    \bottomrule
    \end{tabular}
    }
    \caption{Test results for cross-jurisdiction transfer. We present results in four settings: \emph{standard} (Original)  \emph{augmented} (+ MT Swiss), \emph{further augmented incl. cross-jurisdiction} (+ MT Swiss + MT Indian) and \emph{zero-shot} (MT Indian). Best results are in \textbf{bold}. Diff. shows the difference between the best performing language and the worst performing language (max - min). \textbf{\emph{Further augmenting with translated Indian cases is overall beneficial.}}
    }
    \label{tab:cross_jur}
    \vspace{-4mm}
\end{table*}

\subsubsection{Legal Areas}
\label{sec:cross_domain_legal_areas}
In Table \ref{tab:cross_domain_legal_areas} we present the results for \emph{cross-domain} transfer between legal areas (domains of law).
The results on the diagonal (\underline{underlined}) are in-domain, i.e., fine-tuned and evaluated in the same legal area. We observe that for each domain, the models trained on in-domain data have the best results in the respective domain compared to the rest.

Interesting to note is that the best results (\textbf{bold}) are achieved in the cross-domain setting in 3 out of 4 legal areas. Such an outcome is not anticipated based on the current trends in law industry, where legal experts (judges, lawyers) over-specialize and excel in specific legal areas, e.g., criminal defense lawyers. Penal law poses the only exception where the domain-specific model is on par with the cross-domain model. Again, the results per area do not correlate with the volume of training data (\emph{cross-domain representational bias}), and suggest that other qualitative characteristics (e.g., the idiosyncrasies of criminal law) affect the task complexity.

Similarly to the cross-regional experiments, the shared multilingual model (XLM-R) trained across all languages and legal areas with an augmented dataset outperforms the NativeBERTs models trained in a similar setting, giving another indication that the performance gains from cross-lingual transfer and data augmentation via machine translation are robust across domains as well.

\subsection{Cross-Jurisdiction Transfer}
\label{sec:cross_jurisdiction_exps}

We, finally, ``ambitiously'' stretch the limits of transfer learning in \ac{LJP} and we apply \emph{cross-jurisdiction} transfer, i.e., use of cases from different legal systems, another form of cross-domain transfer. For this purpose, we further augment the SJP dataset of \ac{FSCS} cases, with cases from the \ac{SCI}, published by \citet{malik-etal-2021-ildc}.\footnote{Although the \ac{SCI} rules under the Indian jurisdiction (law), while the \ac{FSCS} under the Swiss one, we hypothesize that the fundamentals of law in two modern legal systems are quite common and thus transferring knowledge could potentially have a positive effect. We discuss this matter in Section~\ref{sec:cross_jur_motivation}.} We consider and translate all (approx.\ 30K) Indian cases ruled up to the last year (2014) of our training dataset, originally written in English, to all target languages (German, French, and Italian).\footnote{We do not use the original documents written in English, as English is not one of our target languages.}

In Table~\ref{tab:cross_jur}, we present the results for two cross-jurisdiction settings: \emph{zero-shot}  (Only MT Indian), where we train XLM-R on the machine-translated version of Indian cases, and \emph{further augmented} (Original + MT Swiss + MT Indian), where we further augment the (already augmented) training set of Swiss cases with the translated Indian ones. While zero-shot transfer clearly fails; interestingly, we observe improvement for all languages in the further augmented setting. This opens a fascinating new direction for \ac{LJP} research.

Similar to our results in Section~\ref{sec:cross_lingual_exps} with respect to cross-lingual performance parity, the standard adapter-based XLM-R model has also the highest performance parity (least diff. on Table~\ref{tab:cross_jur}), while the same model trained on the fully augmented dataset leads to the worst-case (language; German) performance and best performance for the least represented language (Italian).
\vspace{3mm}

\noindent The cumulative improvement from all applied enhancements adds up to 7\% macro-F1 compared to the XLM-R baseline and 16\% to the best method by \citet{niklaus-etal-2021-swiss} in the low-resource Italian subset, while using cross-lingual and cross-jurisdiction transfer we improve for 2.3\% overall and 4.6\% for Italian over our strongest baseline (NativeBERTs).

Since our experiments present several incremental improvements, we assess the stability of the performance improvements with statistical significance testing by comparing the most crucial settings in Appendix \ref{sec:statistical}.

\section{Related Work}
\label{sec:related_work}

\paragraph{Legal Judgment Prediction}
(\ac{LJP}) is the task, where given the facts of a legal case, a system has to predict the correct outcome (legal judgement).
 Many prior works experimented with some forms of \ac{LJP},
however, the precise formulation of the \ac{LJP} task is non-standard as the jurisdictions and legal frameworks vary.
\citet{aletras_predicting_2016,Medvedeva2018, chalkidis-etal-2019-neural} predict the plausible violation of \ac{ECHR} articles of the \ac{ECtHR}. \citet{xiao_cail2018_2018,xiao2021} study Chinese criminal cases where the goal is to predict the ruled duration of prison sentences and/or the relevant law articles.

Another setup is followed by \citet{sulea_predicting_2017,malik-etal-2021-ildc,niklaus-etal-2021-swiss}, which use cases from Supreme Courts (French, Indian, Swiss, respectively), hearing appeals from lower courts relevant to several fields of law (legal areas). Across tasks (datasets), the goal is to predict the binary verdict of the court (approval or dismissal of the examined appeal) given a textual description of the case.
None of these works have explored neither cross-lingual nor cross-jurisdiction transfer, while the effects of cross-domain and cross-regional transfer are also not studied.

\paragraph{Cross-Lingual Transfer}

(\ac{CLT}) is a flourishing topic with the application of pre-trained transformer-based models trained in a multilingual setting \cite{devlin_bert_2019, lample_conneau_2019,conneau_unsupervised_2020, xue-etal-2021-mt5} excelling in NLU benchmarks \cite{ruder-etal-2021-xtreme}. Adapter-based fine-tuning \cite{houlsby2019parameterefficient,pfeiffer-etal-2021-adapterfusion} has been proposed as an anti-measure to mitigate misalignment of multilingual knowledge when \ac{CLT} is applied, especially in a zero-shot fashion, where the target language is unseen during training (or even pre-training).

Meanwhile, \ac{CLT} is understudied in legal \ac{NLP}  applications. \citet{chalkidis-etal-2021-multieurlex} experiment with standard fine-tuning, while they also examined the use of adapters \cite{houlsby2019parameterefficient} for zero-shot \ac{CLT} on a legal topic classification dataset comprising \ac{EU} laws. They found adapters to achieve the best tradeoff between effectiveness and efficiency. Their work did not examine the use of methods incorporating translated versions of the original documents in any form, i.e., translate train documents or test ones. Recently, \citet{xenouleas-etal-2022} used an updated, unparalleled version of \citeauthor{chalkidis-etal-2021-multieurlex} dataset to study \ac{NMT} -augmented CLT methods.  Other multilingual legal \ac{NLP}  resources \cite{galassi-etal-2020-cross,drawzeski-etal-2021-corpus} have been recently released, although \ac{CLT} is not applied in any form.

\section{Motivation and Challenges for Cross-Jurisdiction Transfer}
\label{sec:cross_jur_motivation}

Legal systems vary from country to country. Although they develop in different ways, legal systems also have some similarities based on historically accepted justice ideals, i.e., the rule of law and human rights. Switzerland has a civil law legal system \cite{walther_2001}, i.e., statutes (legislation) is the primary source of law, at the crossroads between Germanic and French legal traditions.

Contrary, India has a hybrid legal system with a mixture of civil, common law, i.e., judicial decisions have precedential value, and customary, i.e., Islamic ethics, or religious law \cite{indian-law-reuters}. The legal and judicial system derives largely from the British common law system, coming as a consequence of the British colonial era (1858-1947) \cite{singh-etal-2019}.

Based on the aforementioned, cross-jurisdiction transfer is challenging since the data (judgments) abide to different law standards.
Although the Supreme Court of India (SCI) rules under the Indian jurisdiction (law), while the Federal Supreme Court of Switzerland (FSCS) under the Swiss one, we hypothesize that the fundamentals of law in two modern legal systems are quite common and thus transferring knowledge could potentially have a positive effect, and thus it is an experiment worth considering, while we acknowledge that from a legal perspective equating legal systems is deeply problematic, since the legislation, the case law, and legal practice are different.

Our empirical work and experimental results shows that cross-jurisdiction transfer in this specific setting (combination of Swiss and Indian decisions) has a positive impact in performance, but we cannot provide any profound hypothesis neither we are able to derive any conclusions on the importance of this finding on legal literature and practice. We leave these questions in the hands of those who can responsibly bear the burden, the legal scholars.

\section{Conclusions and Future Work}
\label{sec:conclusions_and_future_work}

\subsection{Answers to the Research Questions}
Following the experimental results (Section~\ref{sec:experiments}), we answer the original predefined research questions:
\vspace{1mm}

\noindent\textbf{RQ1}: \emph{Is cross-lingual transfer beneficial across all or some of the languages?}
In Section~\ref{sec:cross_lingual_exps}, we find that vanilla \ac{CLT} is beneficial in a low-resource setting (Italian), with comparable results in the rest of the languages. Moreover, \ac{CLT} leveraging \ac{NMT} -based data augmentation is beneficial across all languages. Overall, our experiments lead to a single multi-lingual cross-lingually ``fairer'' model.\vspace{1mm}

\noindent\textbf{RQ2}: \emph{Do models benefit or not from cross-regional and cross-domain transfer?}
In Section~\ref{sec:cross_domain_exps}, we find that models benefit from cross-regional transfer across all cases, since they are exposed to (trained in) many more documents (cases). We believe cross-regional diversity is not a significant aspect, compared to the importance of the increased data volume and language diversity. Cross-domain transfer is beneficial in three out of four cases (legal areas), with comparable results on penal (criminal) law, where the application of law seems to be more straight-forward / standardized (higher performing legal area). Cross-regional and cross-domain transfer lead to more robust models.\vspace{1mm}

\noindent\textbf{RQ3}: \emph{Can we leverage data from another jurisdiction to improve performance?}
In Section~\ref{sec:cross_jurisdiction_exps}, we find that cross-jurisdiction transfer in our specific setup, i.e.,  very similar \ac{LJP} tasks, is beneficial. Again, we believe that this is mostly a matter of additional unique data (cases), rather than a matter of jurisdictional similarity. Cross-jurisdiction transfer leads to a better performing model.\vspace{1mm}

\noindent\textbf{RQ4}: \emph{How does representational bias (wrt. language, origin region, legal area) affect model's performance?}
We observe that representational bias -- in non-extreme cases (e.g., w.r.t. language) -- does not always explain performance disparities across languages, regions, or domains, and other characteristics also need to be considered.

\subsection{Conclusions - Summary}
We examined the application of \acf{CLT} in \acf{LJP} for the very first time, finding a multilingually trained model to be superior when augmenting the dataset with \ac{NMT}.
Adapter-based fine-tuning leads to even better results.
We also examined the effects of cross-domain (legal areas) and cross-regional transfer, which is overall beneficial in both settings, leading to more robust models.
Cross-jurisdiction transfer by augmenting the training set with machine-translated Indian cases further improves performance.

\subsection{Future Work}
In future work, we would like to explore the use of a legal-oriented multilingual pre-trained model by either continued pre-training of XLM-R, or pre-training from scratch in multilingual legal corpora. Legal \ac{NLP}  literature \cite{chalkidis-etal-2022-lexglue, Zheng2021} suggests that domain-specific language models positively affect performance.

In another interesting direction, we will consider other data augmentation techniques \cite{feng-etal-2021-survey, ma2019nlpaug} that rely on textual alternations (e.g., paraphrasing, etc.). We would also like to further investigate cross-jurisdictional transfer, either exploiting data for similar \ac{LJP}  tasks, or via multi-task learning on multiple \ac{LJP}  datasets with dissimilar task specifications.

\vspace{-2mm}
\section{Ethics Statement}
\vspace{-2mm}
The scope of this work is to study \ac{LJP} to broaden the discussion and help practitioners to build assisting technology for legal professionals and laypersons. We believe that this is an important application field, where research should be conducted \cite{tsarapatsanis-aletras-2021-ethical} to improve legal services and democratize law, while also highlight (inform the audience on) the various multi-aspect shortcomings seeking a responsible and ethical (fair) deployment of legal-oriented technologies.

In this direction, we study how we could better exploit all the available resources (from various languages, domains, regions, or even different jurisdictions). This combination leads to models that improve overall performance -- more robust models --, while having improved performance in the worst-case scenarios across many important demographic or legal dimensions (low-resource language, worst performing legal area and region).

Nonetheless, irresponsible use (deployment) of such technology is a plausible risk, as in any other application (e.g., online content moderation) and domain (e.g., medical). We believe that similar technologies should only be deployed to assist human experts (e.g., legal scholars in research, or legal professionals in forecasting or assessing legal case complexity) with notices on their limitations.

The main examined dataset, \acf{SJP}, released by \citet{niklaus-etal-2021-swiss}, comprises publicly available cases from the \ac{FSCS}, where cases are pre-anonymized, i.e., names and other sensitive information are redacted. The same applies for the second one, \acf{ILDC} of \citet{malik-etal-2021-ildc}.

\section*{Acknowledgements}
This work has been supported by the Swiss National Research Program “Digital Transformation” (NRP-77)\footnote{\url{https://www.nfp77.ch/en/}} grant number 187477. This work is also partly funded by the Innovation Fund Denmark (IFD)\footnote{\url{https://innovationsfonden.dk/en}} under File No.\ 0175-00011A.  This research has been also co‐financed by the European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation, under the call RESEARCH – CREATE – INNOVATE (Τ2ΕΔΚ-03849).

We would like to thank Thomas Lüthi for his legal advice, Mara Häusler for great discussions regarding the evaluation process of the models, and Phillip Rust and Desmond Elliott for providing valuable feedback on the original draft of the manuscript.