\section{Introduction}{\label{Intro}}

Newton-Cartan gravity is the geometric reformulation of the Newtonian gravity that is described by a degenerate spatial metric $h^{\mu\nu}$, a temporal vielbein $\tau_\mu$ and a $U(1)$ connection $m_\mu$\footnote{In this paper, the Greek indices $(\mu,\nu,\ldots)$ refer to the coordinate frame and labels all spacetime coordinates with $x = (t, \vec{x})$ while the Latin alphabet letters $(a,b,\ldots)$ refer to the spatial local Galilean frame.}. This geometric setup of Newtonian gravity is realized based on two principles: (i) The geodesic equation is equivalent to the classical equation of motion of a massive particle, and (ii) The only non-vanishing component of the Riemann tensor gives rise to the Poisson equation of Newtonian gravity \cite{Cartan1,Cartan2}. If the Riemann tensor satisfies the Trautman \cite{Trautman} and the Ehlers \cite{Ehlers} conditions, then we are led to the so-called Newtonian connection that is the only non-vanishing component of the connection is given by $\Gamma^a_{00} = \delta^{ab} \partial_b \phi$, satisfying (i). Consequently, the only non-vanishing component of the Riemann tensor becomes the desired Poisson's equation for Newtonian gravity
\begin{eqnarray}
R^a{}_{0a0} (\Gamma)= \nabla^2 \phi = 4 \pi G \rho \,,
\label{NewtonCartan}
\end{eqnarray}
where $G$ is Newton's constant and $\rho$ is the mass density. An alternative derivation of this geometric formulation can be achieved by gauging the underlying symmetry algebra that is the Bargmann algebra \cite{EricPanda}. The Bargmann algebra consists of the generators of spatial rotations $J_{ab}$, spatial translations $P_a$, Galilean boosts $G_a$, time translations $H$ and a central charge $M$, corresponding to particle mass. In the gauging procedure, one first associates each generator with a gauge field
\begin{eqnarray}
J_{ab} \rightarrow \omega_\mu{}^{ab}, \quad P_a \rightarrow e_\mu{}^a, \quad G_a \rightarrow \omega_\mu{}^a, \quad H \rightarrow \tau_\mu, \quad M \rightarrow m_\mu \,,
\end{eqnarray}
and imposes constraints on the curvatures of $\tau_\mu, e_\mu{}^a$ and $m_\mu$
which connect the gauge transformations to general coordinate transformations and solve $\omega_\mu{}^{ab}$ and $\omega_\mu{}^a$ in terms of $e_\mu{}^a, \tau_\mu$ and $m_\mu$, leaving the triplet $(e_\mu{}^a, \tau_\mu,m_\mu)$ as the fundamental fields of the gauge theory formulation of Newton-Cartan gravity. Finally, one imposes two further constraints on the curvatures of $\omega_\mu{}^{ab}$ and $\omega_\mu{}^a$ and recover the Poisson equation of Newtonian gravity and the geodesic equation for a massive particle from a gauge theory viewpoint.

In three space-time dimensions, the Einstein gravity can be established as a Chern-Simons gauge theory in which case the necessary curvature constraints arise as an equation of motion \cite{3dWitten}. This line of thinking encourages one to build up a similar three-dimensional framework for the Newton-Cartan gravity. Indeed in \cite{EBG1} this problem was addressed and it was shown that the three-dimensional Bargmann algebra does not admit a non-degenerate invariant bi-linear form to construct a Chern-Simons theory of gravity but an extension of the Bargmann algebra with a second central charge $S$ is necessary. The resulting algebra is called the extended Bargmann algebra and the corresponding Chern-Simons model is called the extended Bargmann gravity (EBG) \cite{EBG1}. The EBG is in many ways different than Newton-Cartan gravity. First of all, while the Newton-Cartan gravity dictates that only the purely time-like component of the Ricci tensor is non-zero (\ref{NewtonCartan}), EBG allows all components of the Ricci tensor to be non-zero in the presence of matter coupling \cite{Bergshoeff EBG}. Second, the extended Bargmann algebra further allows a non-relativistic conformal extension and a corresponding action principle, called the extended Schr\"odinger gravity \cite{Hartong EBG}. This extension is noteworthy as it corresponds to a conformal non-projectable Horava-Lifshitz gravity \cite{Hartong EBG} through a map between the Newton-Cartan geometry and Horava-Lifshitz gravity \cite{Horava}. Finally, the extended Bargmann algebra admits a supersymmetric extension and a corresponding on-shell Chern-Simons supergravity action \cite{Bergshoeff EBG}. The existence of a supersymmetric action is particularly important as it signals the possibility of an off-shell action for the extended Bargmann supergravity which can serve as a starting point to construct non-relativistic field theories on curved backgrounds by means of localization \cite{Local1,Local2}.

The on-shell $\mathcal{N}=2$ extended Bargmann supergravity has the following field content \cite{Bergshoeff EBG}
\begin{eqnarray}
\{e_\mu{}^a, \tau_\mu, m_\mu, s_\mu, \psi_\mu^+, \psi_\mu^-, \rho_\mu \} \,,
\end{eqnarray}
where the fermionic fields $\psi^\pm_\mu$ and $\rho_\mu$ are Majorana fermions that are associated with the fermionic generators $Q^\pm$ and $R$, respectively. Here, we used the $\mathcal{N}=2$ terminology as the on-shell EBG can be derived by a Lie algebra expansion of the three-dimensional $\mathcal{N}=2$ Poincar\'e supergravity \cite{LAE5}. The construction of an off-shell extended Bargmann supergravity requires the addition of auxiliary fields and the application of localization techniques requires the knowledge of off-shell matter multiplet actions of extended Bargmann superalgebra. In the relativistic context, this is most straightforwardly achieved by applying superconformal tensor calculus \cite{SCTC1,SCTC2,SCTC3,SCTC4}, where one first constructs superconformal models then gauge fix the redundant conformal symmetries to obtain a super-Poincar\'e invariant theory. In the case of non-relativistic gravity, a methodology for a conformal tensor calculus was established for the bosonic and supersymmetric conformal extension of the Galilei algebra, known as the Schr\"odinger algebra, in \cite{AfsharBergshoeff} and \cite{SuperSch}, respectively.

In this paper, our purpose is to take a pioneering step towards an extended Schr\"odinger tensor calculus by establishing the Schr\"odinger extension of the extended Bargmann algebra and the corresponding extended Schr\"odinger supergravity. Due to the map between the Newton-Cartan geometry and Horava-Lifshitz gravity, this result also corresponds to a superconformal non-projectable Horava-Lifshitz gravity. To achieve the extended Schr\"odinger superalgebra, we first remind the reader in Section \ref{Section2} about the fundamentals of the extended Bargmann algebra and its cosmological extension, known as the extended Newton-Hooke algebra. In Section \ref{Section3}, we construct the cosmological extension of the extended Bargmann superalgebra. The cosmological extension of the algebra allows us to construct the extended Newton-Hooke and the extended exotic $\mathcal{N}=(2,0)$ Bargmann supergravity. Here, $\mathcal{N}=(2,0)$ terminology refers to the fact that the non-relativistic algebra can be understood as the Lie algebra expansion of the three-dimensional $\mathcal{N}=(2,0)$ AdS algebra. In Section \ref{Section4} we start adding dilatation and non-relativistic special conformal generators to the extended Bargmann algebra. In doing so, we first add the dilatations $D$ and a central charge $Y$. Perhaps unexpectedly, we found that the superalgebra can be closed by adding two $R-$symmetry generators and three fermionic generators $Q^\pm$ and $R$. The presence of the central charge $Y$ fixes the so-called dynamical critical exponent $z$ to $z=2$ and we establish $z=2$ extended Lifshitz supergravity. In the final step, we add the non-relativistic special conformal generator $K$ and a second central charge $Z$. The supersymmetric extension of the extended Schr\"odinger algebra requires two additional fermionic generators, which we call $F^\pm$. With the addition of the new fermionic generators, we close the superalgebra and construct the Chern-Simons action for the extended Schr\"odinger gravity.

\section{Extended Newton-Hooke Algebra and Chern-Simons Actions} \label{Section2}

In this first section, we briefly review the extended Newton-Hooke gravity and extended exotic Bargmann gravity as well as the underlying symmetry algebra that is the extended Newton-Hooke algebra. The Bargmann algebra consists of the generators of spatial rotations $J$, spatial translations $P_a$, Galilean boosts $G_a$, time translations $H$ and a central charge $M$. The extended Bargmann algebra extends the Bargmann algebra with a second central charge $S$ \cite{EBG1,Bergshoeff EBG,Hartong EBG}. The non-vanishing commutation relations for the extended Bargmann algebra are given by
\begin{align}
\left[ H, G_a \right]  &= - \epsilon_{ab} P^b\,,          &  \left[ J, P_a \right] &= - \epsilon_{ab} P^b\,,          &  \left[ J, G_a \right] &=- \epsilon_{ab} G^b\,, \nonumber\\
\left[ G_{a}, P_b \right]  &= \epsilon_{ab} M\,,        &  \left[ G_a, G_{b} \right]  &= \epsilon_{ab} S \,,
\label{EBA}
\end{align}
Here, we use the conventions of \cite{Bergshoeff EBG}. For reader’s convenience, we present a map between different conventions in Appendix \ref{AppB}. The cosmological extension of the extended Bargmann algebra requires two extra commutators
\begin{eqnarray}
&&   \left[ H, P_a \right] = - \frac{1}{\ell^2} \epsilon_{ab} G^b \,,\quad \left[ P_a, P_a \right] = \frac{1}{\ell^2} \epsilon_{ab} S \,.
\label{CEBA}
\end{eqnarray}
Note that we set $\Lambda = \frac{1}{\ell^2}$ for future purposes, however the opposite sign is equally well to construct the cosmological bosonic models.  The construction of actions for this symmetry algebra is based on the Chern-Simons action
\begin{eqnarray}
S = \frac{k}{4\pi} \int \rm{STr} \left(A \wedge dA + \frac23 A \wedge A \wedge A \right) \,,
\label{CS}
\end{eqnarray}
where $k$ is the Chern-Simons coupling constant, the gauge field A represents a Lie-algebra-valued one form and ``STr" represents the supertrace using the non-degenerate invariant bilinear form on the Lie algebra, which turns into the standard trace for the purely bosonic models. For the extended Bargmann algebra, the gauge field $A_\mu$ is given by
\begin{eqnarray}
A_\mu = \tau_\mu H + e_\mu{}^a P_a + \omega J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S \,.
\end{eqnarray}
The transformation rules for these gauge fields can be found by the standard rule
\begin{eqnarray}
\delta A_\mu^A = \partial_\mu \epsilon^A + f_{BC}{}^A \epsilon^C A_\mu^B \,,
\label{StandardRule}
\end{eqnarray}
where $\epsilon^A$ is the relevant gauge parameter and $f_{BC}{}^A$ are the structure constants. The covariant curvatures are given by
\begin{eqnarray}
R_{\mu\nu}(H) &=& 2 \partial_{[\mu} \tau_{\nu]}   \,,\nonumber\\
R_{\mu\nu}{}^a (P) &=&  2 \partial_{[\mu} e_{\nu]}{}^a + 2 \epsilon^{ab}\,  \omega_{[\mu} e_{\nu]b} - 2 \epsilon^{ab}\,  \omega_{ [\mu b} \tau_{\nu]} \,,\nonumber\\
R_{\mu\nu} (J) &=& 2 \partial_{[\mu} \omega_{\nu]} \,,\nonumber\\
R_{\mu\nu}{}^a (G) &=& 2 \partial_{[\mu} \omega_{\nu]}{}^a + 2 \epsilon^{ab} \omega_{[\mu} \omega_{\nu]b} + \frac{2}{\ell^2} \epsilon^{ab} \tau_{[\mu} e_{\nu]b}  \,,\nonumber\\
R_{\mu\nu}(M) &=& 2 \partial_{[\mu} m_{\nu]} + 2 \epsilon_{ab} \omega_{[\mu}{}^a e_{\nu]}{}^b   \,,\nonumber\\
R_{\mu\nu}(S) &=& 2 \partial_{[\mu} s_{\nu]} +  \epsilon_{ab} \omega_{[\mu}{}^a \omega_{\nu]}{}^b + \frac{1}{\ell^2}  \epsilon_{ab} e_{[\mu}{}^a e_{\nu]}{}^b  \,.
\label{Curvatures1}
\end{eqnarray}
Note that for $\ell \rightarrow \infty$, these curvatures become the curvatures of the extended Bargmann algebra.

As mentioned, the construction of a Chern-Simons action requires non-degenerate invariant bilinear forms which can simply be found by constructing a Casimir operator $C$. In the case of extended Newton-Hooke algebra the following bilinears are of interest
\begin{itemize}
\item {$(P_a,G_b) = \delta_{ab}\,, \quad (H,S) =-1\,, \quad (J,M) = -1$}
\item{$(P_a,P_b) =  \frac{1}{\ell} \delta_{ab} \,, \quad (G_a,G_b) = \ell  \delta_{ab}\,, \quad (H,M) =- \frac{1}{\ell} \,,\quad  (J,S) = -\ell\,, $}
\end{itemize}
In the first case, the cosmological constant does not appear in the bilinears hence it can also be used to construct the extended Bargmann gravity without a cosmological constant. On the other hand, the second set only exists in the presence of a cosmological constant. A more detailed discussion on the invariant bilinear forms, including the degenerate ones in context of extended Bargmann algebra, can be found in \cite{Hartong EBG}.

\subsection{Extended Newton-Hooke Gravity}

For the construction of the extended Newton-Hooke gravity we utilize the first set of invariant bilinear forms
\begin{eqnarray}
(P_a, G_b) = \delta_{ab} \,, \quad (H,S) =  -1 \,, \quad (J,M) = -1\,,
\label{Bilinear1}
\end{eqnarray}
and the structure constants of the extended Newton-Hooke algebra (\ref{EBA}) and (\ref{CEBA}). In this case, the Chern-Simons action (\ref{CS}) read
\begin{eqnarray}
S &=&\frac{k}{4 \pi} \int d^3x\, \epsilon^{\mu\nu\rho}  \left( e_\mu{}^a R_{\nu\rho a} (G) - \tau_\mu R_{\nu\rho}(S) - m_\mu R_{\nu\rho}(J) + \frac{2}{\ell^2} \epsilon_{ab} \tau_\mu e_\nu{}^a e_\rho{}^b \right) \,, \quad
\label{CEBG}
\end{eqnarray}
where the curvatures are as defined in (\ref{Curvatures1}). For $\ell \rightarrow \infty$, this action coincides with the extended Bargmann gravity of \cite{EBG1,Bergshoeff EBG,Hartong EBG}. Interestingly,  the action for the extended Newton-Hooke gravity can also be obtained by a contraction procedure, mimicking the Inonu-Wigner contraction  \cite{Bergshoeff EBG}\footnote{Extended Bargmann algebra can also be obtained by the In\"on\"u-Wigner contraction of the Poincar\'e $\otimes$ $U(1)^2$ algebra, see \cite{MaxwellContraction,Color} }. In order to do so, one needs to consider the cosmological Einstein-Hilbert action plus a Chern-Simons action for two abelian gauge fields $Z_{1\mu}$ and $Z_{2\mu}$ with an off-diagonal coupling \cite{Bergshoeff EBG}
\begin{eqnarray}
S &=&   \frac{k \omega}{4 \pi }\int  d^3 x \  \epsilon^{\mu \nu \rho} \left( \eta_{AB}  E_\mu{}^A R_ {\nu\rho}{}^B(\Omega) + \frac{1}{3L^2}\   \epsilon_{ABC}\ E_\mu{}^A   E_\nu{}^B E_\rho{}^C + 2 Z_{1 \mu} \partial_\nu Z_{2 \rho}   \right)\,, \label{EHV}
\end{eqnarray}
where
\begin{eqnarray}
R_{\mu\nu}{}^A (\Omega) &=& 2 \partial_{ [ \mu } \Omega_{\nu]}{}^A  + \epsilon^{ABC} \Omega_{\mu B} \Omega_{\nu C} \,.
\end{eqnarray}
In terms of the non-relativistic fields and parameters, the relativistic fields and parameters are  \cite{Bergshoeff EBG}
\begin{eqnarray}
&& E_\mu{}^a = e_\mu{}^a  \,, \quad E_\mu{}^0  = \omega  \tau_\mu + \frac{1}{2\omega} m_\mu  \,, \quad   \Omega_\mu{}^0= \omega_\mu + \frac{1}{2 \omega ^2 } s_\mu \,,\quad  L = \omega \ell \,, \nonumber \\
&&  \Omega_\mu{}^a  = \frac{1}{\omega}\omega_\mu{}^a\,, \quad Z_{1 \mu} = \omega \tau_\mu - \frac{1}{2 \omega } m_\mu \,, \quad Z_{2\mu} = \omega_\mu - \frac{1}{2 \omega^2}s_\mu \,.
\label{contraction1}
\end{eqnarray}
The off-diagonal coupling of the vector fields are essential for removing the $\omega^2$ divergencies. In the $\omega \rightarrow \infty$ limit, the action (\ref{EHV}) is precisely recover the extended Newton-Hooke gravity action (\ref{CEBG}).

\subsection{Extended Exotic Bargmann Gravity}

In three dimensions there is an alternative action that gives rise to the Einstein's equations in the presence of a cosmological constant, known as the exotic Einstein-Hilbert action \cite{3dWitten,TownsendExotic}
\begin{eqnarray}
S =  \frac{kL \omega  }{2 \pi }\int  d^3 x  \  \epsilon^{\mu \nu \rho} \left(\Omega_\mu{}^A \partial_ \nu \Omega_{\rho A}  + \frac{1}{3} \epsilon_{ABC} \Omega_\mu{}^A \Omega_\nu{}^B \Omega_\rho{}^C  +   \frac{1}{2 L^2} E_\mu{}^A  R_{\nu\rho A}(P)  \right) \,.
\label{ExoticEinstein}
\end{eqnarray}
This indicates that there might be an extended ``exotic" Bargmann gravity following the contraction procedure introduced in (\ref{contraction1}). As the exotic theory has diagonal coupling amongst the gauge fields of the Poincar\'e algebra, the vector fields that are neccessary to cancel out the $\omega^2$ divergences should also appear with a diagonal coupling. Thus, we consider the exotic Einstein-Hilbert action plus a diagonal Chern-Simon action for two abelian gauge fields
\begin{eqnarray}
S &=& \frac{kL\omega}{2 \pi }\int  d^3 x  \  \epsilon^{\mu \nu \rho} \Big(  \Omega_\mu{}^A \partial_ \nu \Omega_{\rho A}  + \frac{1}{3} \epsilon_{ABC} \Omega_\mu{}^A \Omega_\nu{}^B \Omega_\rho{}^C  +   \frac{1}{2 L^2} E_\mu{}^A  R_{\nu\rho A}(P) \nonumber\\
&&  \quad   + \frac{1}{L^2} Z_{1\mu} \partial_\nu Z_{1\rho}  + Z_{2\mu} \partial_\nu Z_{2\rho} \Big) \,,
\label{ExoticZZ}
\end{eqnarray}
where
\begin{eqnarray}
R_{\mu\nu}{}^A (P) &=& 2 \partial_{[\mu} E_{\nu]}{}^A + 2 \epsilon^{ABC} \Omega_{[\mu B} E_{\nu] C} \,.
\end{eqnarray}
Using the expressions in (\ref{contraction1}) in the action (\ref{ExoticZZ}), we find that the extended exotic Bargmann gravity is given by
\begin{eqnarray}
S &=& \frac{k\ell}{4 \pi} \int d^3x\, \epsilon^{\mu\nu\rho}  \Big( \omega_\mu{}^a R_{\nu\rho a} (G) -  2 s_\mu R_{\nu\rho}(J)  + \frac{1}{\ell^2}  \left( e_\mu{}^a R_{\nu\rho a} (P) - 2 m_\mu R_{\nu\rho} (H) \right)  \Big)\,. \qquad \quad
\label{EEBG}
\end{eqnarray}
At this stage, it is worthwhile to discuss the ``exotic" terminology in the context of relativistic and non-relativistic gravity. In the relativistic context, the fundamental fields $E^a$ and $\Omega^{ab}$ are parity even, thus the dual spin connection, $\Omega^a$ is parity odd. As a result the Einstein-Hilbert action (\ref{EHV}) is parity even while the exotic Einstein-Hilbert action (\ref{ExoticEinstein}) is parity odd. Here, the exotic terminology emphasizes the fact that although the exotic Einstein-Hilbert action is of odd-parity, it give rise to parity even equations of motion. In the non-relativistic context, we have the following parity assignment for the gauge fields
\begin{eqnarray}
\text{Even:} \quad   \{ \tau_\mu \,, \quad e_\mu{}^a\,, \quad m_\mu \} \,, \qquad \qquad \text{Odd:}  \quad  \{ \omega_\mu \,, \quad \omega_\mu{}^a\,, \quad s_\mu  \} \nonumber \,.
\end{eqnarray}
In this case, the extended Newton-Hooke gravity is parity even while the extended exotic Bargmann gravity is parity odd. As both of these theories have parity even equations of motion, we use the ``exotic" terminology for the action (\ref{EEBG}). It is also possible to obtain the extended exotic Bargmann gravity by using the second set of bilinear forms \cite{Hartong EBG}
\begin{eqnarray}
(P_a, P_b) =  \frac{1}{\ell^2} \delta_{ab} \,, \quad (G_a, G_b) = \delta_{ab} \,, \quad (H,M) = - \frac{1}{\ell^2}\,,\quad (J,S) = -1 \,,
\end{eqnarray}
and tracing the Chern-Simons action. The resulting action precisely coincide with the extended exotic Bargmann gravity (\ref{EEBG}).

Although the field equations for the extended Newton-Hooke gravity (\ref{CEBG}) and the  extended exotic Bargmann gravity (\ref{ExoticEBG}) are the same, these two models differ dramatically when we consider the matter couplings. In particular, let us first take a closer look at the $s_\mu, \omega_\mu{}^a$ and $\omega_\mu$ equations for the  extended Newton-Hooke gravity
\begin{eqnarray}
R_{\mu\nu}(H) = 0 \,, \qquad  R_{\mu\nu}{}^a (P) = 0 \,, \qquad R_{\mu\nu}(M) = 0 \,.
\end{eqnarray}
The first equation implies that this model is defined on non-relativistic space-times with torsionless Newton-Cartan geometry while the last two equations give rise to composite expressions for $\omega_\mu$ and $\omega_\mu{}^a$. The torsionless condition is imposed employing $s_\mu$ equation, which suggests that one needs to find a matter content that transforms non-trivially under $S$-transformations to enable torsional non-relativistic background geometries. On the other hand, the same torsion-free constraint is imposed by the $m_\mu$ equation in the case of extended exotic Bargmann gravity. Accordingly one needs to find a  matter content that transforms non-trivially under $M$-transformations to include torsion. As opposed to the extended Bargmann gravity, this can be simply achieved, for instance, by coupling a complex scalar field to extended exotic Bargmann gravity. Consequently, extended exotic Bargmann gravity provides a more useful setup to study three-dimensional torsional non-relativistic geometries.

\section{Extended Newton-Hooke and Exotic Bargmann Supergravity} \label{Section3}

The $\mathcal{N}=(2,0)$ supersymmetric extension of the extended Newton-Hooke algebra, (\ref{EBA}) and (\ref{CEBA}), requires three supercharges $Q^\pm_\alpha$ and $R_\alpha$  $(\alpha = 1,2)$ that are all Majorana spinors. It furthermore requires two extra bosonic generators $U_1$ and $U_2$. Both $U_1$ and $U_2$ act non-trivially on the spinors $Q^\pm$ and $R$ in the presence of the cosmological constant and they become central when cosmological constant vanishes. The $\mathcal{N}=(2,0)$ extended Newton-Hooke superalgebra has the following non-vanising $\left[B,B\right]$ and $\left[B,F\right]$ commutators
\begin{align}
\left[ H, G_a \right]  &= - \epsilon_{ab} P^b\,,          &  \left[ J, P_a \right] &= - \epsilon_{ab} P^b\,,          &  \left[ J, G_a \right] &=- \epsilon_{ab} G^b\,, \nonumber\\
\left[ G_{a}, P_b \right]  &= \epsilon_{ab} M\,,        &  \left[ G_a, G_{b} \right]  &= \epsilon_{ab} S \,, &    \left[ P_a, P_{b} \right]  &= \frac{1}{\ell^2} \epsilon_{ab} S   \,, \nonumber\\
\left[ H, P_a \right]  &= - \frac{1}{\ell^2} \epsilon_{ab} G^b\,, & [J, Q^\pm ] & =  - \frac{1}{2} \gamma_0 Q^\pm  \,,  & [J, R ] &= - \frac{1}{2} \gamma_0 R \,,\nonumber\\
[S, Q^+ ] & = - \frac{1}{2} \gamma_0 R \,, & [G_a, Q^+] &= - \frac{1}{2} \gamma_a Q^- \,,& [G_a, Q^-] &= - \frac{1}{2} \gamma_a R \,,\nonumber\\
[P_a, Q^+] &= - \frac{1}{2\ell} \gamma_a Q^-\,, & [P_a, Q^-] &= - \frac{1}{2\ell} \gamma_a R \,, & [H, Q^\pm] &= - \frac{1}{2\ell} \gamma_0 Q^\pm \,,\nonumber\\
[H, R] &= - \frac{1}{2\ell} \gamma_0 R \,, & [M, Q^+ ] &= - \frac{1}{2\ell} \gamma_0 R \,, &[U_1, Q^\pm] &= \mp \frac{1}{\ell}  \gamma_0 Q^\pm \,,\nonumber\\
[U_1, R] &= - \frac{1}{\ell}  \gamma_0 R\,, & [U_2, Q^+ ] &= - \frac{1}{\ell}  \gamma_0 R \,.
\label{CosmoEBG1}
\end{align}
while the non-vanishing $\{F, F\}$ anti-commutators are given by
\begin{eqnarray}
\{ Q^+_\alpha ,  Q^+_\beta  \} &=& (\gamma_0 C^{-1})_{\alpha \beta} H  +   \frac{1}{\ell}  (\gamma_0 C^{-1})_{\alpha \beta} J  -   (\gamma_0 C^{-1})_{\alpha \beta} U_1  \,,\nonumber\\
\{ Q^+_\alpha ,  Q^-_\beta  \}  &=& -(\gamma_a C^{-1})_{\alpha \beta} P^{a} -  \frac{1}{\ell}   (\gamma_a C^{-1})_{\alpha \beta} G^{a}\,,\nonumber\\
\{ Q^+_\alpha ,  R_\beta  \} &=& (\gamma_0 C^{-1})_{\alpha \beta} M  +  \frac{1}{\ell} (\gamma_0 C^{-1})_{\alpha \beta} S  -   (\gamma_0 C^{-1})_{\alpha \beta} U_2\,,\nonumber \\
\{ Q^-_\alpha ,  Q^-_\beta  \} &=& (\gamma_0 C^{-1})_{\alpha \beta} M  +   \frac{1}{\ell}  (\gamma_0 C^{-1})_{\alpha \beta} S  +  (\gamma_0 C^{-1})_{\alpha \beta} U_2  \,.
\label{CosmoEBG2}
\end{eqnarray}
The appearence of the extra bosonic generators $U_{1,2}$ can be understood as  the Lie algebra expansion of the $R$-symmetry generator of the $\mathcal{N}=(2,0)$ AdS superalgebra, which we defer the details to the Appendix \ref{AppA}. Here, our splitting of the indices of gamma matrices into its temporal and the spatial components as well as our treatment to gamma-matrices in general follow from \cite{SezBerg}. As in the bosonic case, the extended Newton-Hooke algebra admit two distinct non-degenerate bilinear forms. In the first case, it is given by
\begin{eqnarray}
&& (P_a,G_b) = \delta_{ab}\,, \quad (H,S) =-1\,, \quad (J,M) = -1\,, \quad (U_1, U_2) = \frac{2}{\ell}\,,\nonumber\\
&&  (Q^+_\alpha,R_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta}\,,\quad (Q^-_\alpha,Q^-_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta} \,,
\label{MetricCEBG1}
\end{eqnarray}
Taking the gauge field to be
\begin{eqnarray}
A_\mu &=& \tau_\mu H + e_\mu{}^a P_a + \omega_\mu J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S +r_{1\mu} U_1 + r_{2\mu} U_2 \nonumber\\
&& + \bar\psi_\mu^+ Q^+ + \bar\psi_\mu^- Q^- + \bar\rho_\mu R \,,
\label{GaugeField2}
\end{eqnarray}
we can construct a Chern-Simons action by using the invariant bilinear form (\ref{MetricCEBG1}). The resulting theory is the extended Newton-Hooke supergravity. This action is given by
\begin{eqnarray}
S &=&\frac{k}{4 \pi} \int d^3x\, \varepsilon^{\mu\nu\rho}  \Big( e_\mu{}^a R_{\nu\rho a} (G) - \tau_\mu R_{\nu\rho} (S) - m_\mu R_{\nu\rho}(J) + \frac{2}{\ell^2}  \epsilon_{ab} \tau_\mu e_\nu{}^a e_\rho{}^b \nonumber \\
&&  + \frac{4}{\ell} r_{1\mu} \partial_{\nu} r_{2\rho} +  \bar\psi_\mu^+ R_{\nu\rho} (R) +   \bar\rho_\mu R_{\nu\rho} (Q^+) +  \bar\psi_\mu^- R_{\nu\rho} (Q^-) \Big) \,.
\label{CosmoEBG}
\end{eqnarray}
The bosonic part of this action, excluding the $R$-symmetry part, matches with the extended Newton-Hooke action given in \cite{Hartong EBG}. Here, the bosonic curvatures are as given in (\ref{Curvatures1}) while the fermionic curvatures are given as
\begin{eqnarray}
R_{\mu\nu} (Q^+) &=& 2 \partial_{[\mu} \psi_{\nu]}^+ +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^+   + \frac{1}{\ell}  \tau_{[\mu} \gamma_0 \psi_{\nu]}^+ +   \frac{2}{\ell} r_{1[\mu} \gamma_0 \psi_{\nu]}^+  \,,\nonumber\\
R_{\mu\nu} (Q^-) &=& 2 \partial_{[\mu} \psi_{\nu]}^- +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^-  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^+    + \frac{1}{\ell} \tau_{[\mu} \gamma_0 \psi_{\nu]}^-  + \frac{1}{\ell} e_{[\mu}{}^a \gamma_a \psi_{\nu]}^+ -  \frac{2}{\ell}  r_{1[\mu} \gamma_0 \psi_{\nu]}^-  \,,\nonumber\\
R_{\mu\nu} (R) &=& 2 \partial_{[\mu} \rho_{\nu]} +  \omega_{[\mu} \gamma_0 \rho_{\nu]}  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^- +  s_{[\mu} \gamma_0 \psi_{\nu]}^+  + \frac{1}{\ell} \tau_{[\mu} \gamma_0 \rho_{\nu]}  + \frac{1}{\ell} m_{[\mu} \gamma_0 \psi_{\nu]}^+   \nonumber\\
&&  +  \frac{2}{\ell} r_{1[\mu} \gamma_0 \rho_{\nu]}  +   \frac{1}{\ell} e_{[\mu}{}^a \gamma_a \psi_{\nu]}^-  +  \frac{2}{\ell} r_{2[\mu} \gamma_0 \psi_{\nu]}^+  \,.
\end{eqnarray}
In the presence of a cosmological constant, there is a second non-degerate invariant bilinear form that is given by
\begin{eqnarray}
&&(P_a,P_b) =  \frac{1}{\ell} \delta_{ab} \,, \quad (G_a,G_b) = \ell  \delta_{ab}\,, \quad (H,M) =- \frac{1}{\ell} \,,\quad  (J,S) = -\ell\,, \nonumber \\
&&(U_1,U_2) = \frac{2}{\ell}\,, \quad (Q^+_\alpha,R_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta}\,,\quad (Q^-_\alpha,Q^-_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta} \,,
\end{eqnarray}
In this case, we obtain the ``exotic" extended Bargmann supergravity
\begin{eqnarray}
S &=& \frac{k \ell}{4\pi} \int d^3x\, \epsilon^{\mu\nu\rho}  \Big(\omega_\mu{}^a R_{\nu\rho a} (G) -  2 s_\mu R_{\nu\rho}(J) + \frac{1}{\ell^2}  e_\mu{}^a R_{\nu\rho a} (P) -  \frac{2}{\ell^2} \tau_\mu R_{\nu\rho} (M)  \nonumber \\
&&  \quad + \frac{8}{\ell^2} r_{2\mu} \partial_\nu r_{1\rho}  +  \frac{2}{\ell}  \bar\psi_\mu^+ R_{\nu\rho} (R)  +  \frac{2}{\ell}  \bar\rho_\mu R_{\nu\rho} (Q^+) + \frac{2}{\ell}   \bar\psi_\mu^- R_{\nu\rho} (Q^-)  \Big) \,.
\label{ExoticEBG}
\end{eqnarray}
The bosonic part of this action, excluding the $R$-symmetry part, matches with the extended exotic Bargmann action given in \cite{Hartong EBG}. The extended Newton-Hooke (\ref{CosmoEBG}) and the exotic extended Bargmann gravity (\ref{ExoticEBG}) are invraiant under the following supersymmetry transformation rules
\begin{eqnarray}
\delta \tau_\mu &=&  - \bar{\epsilon}^+ \gamma_0 \psi^+_\mu \,,\nonumber\\
\delta e_\mu{}^a &=&    \bar{\epsilon}^+ \gamma^{a} \psi_\mu^- +  \bar{\epsilon}^- \gamma^{a} \psi_\mu^+ \,,\nonumber\\
\delta \omega_\mu{}^a &=&  \frac{1}{\ell}(\bar{\epsilon}^+ \gamma^{a} \psi_\mu^- +  \bar{\epsilon}^- \gamma^{a} \psi_\mu^+) \,,\nonumber\\
\delta \omega_\mu &=& - \frac{1}{\ell} \bar{\epsilon}^+ \gamma_0 \psi^+_\mu \,,\nonumber\\
\delta m_\mu &=& - \bar{\epsilon}^- \gamma_0 \psi_\mu^- -  \bar{\epsilon}^+ \gamma_0 \rho_\mu -  \bar{\eta} \gamma_0 \psi_\mu^+  \,,\nonumber\\
\delta s_\mu  &=& - \frac{1}{\ell} \left( \bar{\epsilon}^- \gamma_0 \psi_\mu^- + \bar{\epsilon}^+ \gamma_0 \rho_\mu + \bar{\eta} \gamma_0 \psi_\mu^+\right) \,,\nonumber \\
\delta r_{1\mu}  &=&   \bar{\epsilon}^+ \gamma_0 \psi^+_\mu \,, \nonumber \\
\delta r_{2\mu}  &=&  -  \bar{\epsilon}^- \gamma_0 \psi_\mu^- +  \bar{\epsilon}^ + \gamma_0 \rho_\mu +  \bar{\eta} \gamma_0 \psi_\mu^+ \,, \nonumber \\
\delta \psi_\mu^+ &=& \partial_\mu \epsilon^+ + \frac12 \omega_\mu \gamma_0 \epsilon^+  +  \frac1{2\ell} \tau_\mu \gamma_0 \epsilon^+ + \frac{1}{\ell} r_{1\mu} \gamma_0 \epsilon^+ \,,\nonumber\\
\delta \psi_\mu^- &=&  \partial_\mu \epsilon^- + \frac12 \omega_\mu \gamma_0 \epsilon^- + \frac12 \omega_\mu{}^a \gamma_a \epsilon^+ +  \frac{1}{2\ell} \tau_\mu \gamma_0 \epsilon^-  +  \frac{1}{2\ell}  e_\mu^a \gamma_a \epsilon^+ -  \frac{1}{\ell}  r_{1\mu} \gamma_0 \epsilon^- \,,\nonumber\\
\delta \rho_\mu &=&  \partial_\mu \eta + \frac12 \omega_\mu \gamma_0 \eta + \frac12 \omega_\mu{}^a \gamma_a \epsilon^- + \frac12 s_\mu \gamma_0 \epsilon^+ +   \frac{1}{2\ell}  \tau_\mu \gamma_0 \eta + \frac{1}{2\ell}  m_\mu \gamma_0 \epsilon^+  \nonumber\\
&& +  \frac{1}{2\ell}  e_\mu^a \gamma_a \epsilon^- +   \frac{1}{\ell}  r_{1\mu} \gamma_0 \eta  +  \frac{1}{\ell}  r_{2\mu} \gamma_0 \epsilon^+  \,,
\label{CosmoGaugeEBG}
\end{eqnarray}
where $\epsilon^\pm$ and $\eta$ are the the parameters of the local $Q^\pm$ and $R$ transformations, respectively. In the supersymmetric extension of the extended Newton-Hooke gravity (or extended exotic Bargmann gravity), the curvature constraints that are given by the $s_\mu, \omega_\mu{}^a, \omega_\mu$ (or $m_\mu, e_\mu{}^a, \tau_\mu$) are given by
\begin{eqnarray}
R_{\mu\nu}(H) = 0 \,, \qquad R_{\mu\nu}{}^a (P) = 0 \,, \qquad R_{\mu\nu}(M) = 0 \,,
\end{eqnarray}
where
\begin{eqnarray}
R_{\mu\nu}(H) &=& 2 \partial_{[\mu} \tau_{\nu]}  + \bar\psi_\mu \gamma_0 \psi_\nu   \,,\nonumber\\
R_{\mu\nu}{}^a (P) &=&  2 \partial_{[\mu} e_{\nu]}{}^a + 2 \epsilon^{ab}\,  \omega_{[\mu} e_{\nu]b} - 2 \epsilon^{ab}\,  \omega_{ [\mu b} \tau_{\nu]} - 2 \bar\psi_{[\mu}^+ \gamma^a \psi_{\nu]}^-  \,,\nonumber\\
R_{\mu\nu}(M) &=& 2 \partial_{[\mu} m_{\nu]} + 2 \epsilon_{ab} \omega_{[\mu}{}^a e_{\nu]}{}^b + 2 \bar\psi_{[\mu}^+ \gamma_0 \rho_{\nu]} + \bar\psi_\mu^- \gamma_0 \psi_\nu^-   \,.
\end{eqnarray}
In this case, first equation imply a fermionic temporal torsion while the last two equations gives rise to composite expressions for $\omega_\mu$ and $\omega_\mu{}^a$ that includes a fermionic part. The back substitution of these composite expressions into the action (\ref{CEBG}) (or (\ref{ExoticEBG})) lead to quartic fermionic terms in the action.

\section{Extended Lifshitz and Schr\"odinger Supergravity} \label{Section4}

In this section, we take our final step towards the completion of the realm of three-dimensional $\mathcal{N}=2$ extended Bargmann supergravity. The missing pieces are models of supergravity with more symmetries such as dilatations and non-relativistic special conformal symmetry.  In the first subsection, we show that it is possible to extend the extended Bargmann algebra with dilatations and a central charge, $Y$, giving rise to the extended Lifshitz superalgebra and consequently the extended Lifshitz supergravity. We then proceed with the full Schr\"odinger superalgebra and construct a Chern-Simons action as an extended Schr\"odinger supergravity.

\subsection{Extended Lifshitz Supergravity}

The inclusion of dilatations ($D$) in the extended Bargmann algebra (\ref{EBA}) requires an additional central charge ($Y$)  in order to have a non-degenerate metric. The non-vanishing commutation relations for this algebra are given by
\begin{align}
\left[ H, G_a \right]  &= - \epsilon_{ab} P^b\,,          &  \left[ J, P_a \right] &= - \epsilon_{ab} P^b\,,          &  \left[ J, G_a \right] &=- \epsilon_{ab} G^b\,, \nonumber\\
\left[ P_{a}, G_b \right]  &= \epsilon_{ab} M + \delta_{ab} Y \,,        &  \left[ G_a, G_{b} \right]  &=  \epsilon_{ab} S \,, & \left[ D, S \right]  &=2 S \,,  \nonumber\\
\left[ D, G_a \right]  &= G_{a} \,,          &  \left[ D, P_a \right] &= -  P_{a} \,,          &  \left[ D, H \right] &=- 2 H\,, \nonumber\\
\left[ H, S \right]  &= - 2 Y  \,,
\label{LifBoson}
\end{align}
We refer to this algebra as the extended extended Lifshitz algebra. Note here that the Lifshitz symmetry usually allow the scaling dimension of the temporal translations $H$ to be arbitrary, which is referred to as the dynamical critical exponent $z$. However, the inclusion of the central charge $Y$ fixes the scaling dimension of $H$ to $z=2$. This algebra allows the following non-degenerate invariant bi-linear form
\begin{eqnarray}
(P_{a}, G_b) = \delta_{ab} \,, \quad (S,H) = -1 \,, \quad (M, J) = -1\,, \quad (Y,D) = -1 \,. \label{metricbosonicLifshitz}
\end{eqnarray}
Taking the gauge field to be
\begin{eqnarray}
A_\mu &=& \tau_\mu H + e_\mu{}^a P_a + \omega_\mu J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S + b_\mu D + y_\mu Y \,,
\label{GFLifBos}
\end{eqnarray}
 we can use to construct a Chern-Simons action for the extended Lifshitz algebra
\begin{eqnarray}
S &=&\frac{k}{4 \pi} \int d^3x\, \varepsilon^{\mu\nu\rho}  \Big( \omega_\mu{}^a  R_{\nu\rho a} (P) - 2 s_\mu  R_{\nu\rho} (H) - 2 m_\mu  R_{\nu\rho}(J) \nonumber\\
&& \quad + e_\mu{}^a  R_{\nu\rho a} (G)   - 2 y_\mu  R_{\nu\rho} (D)  \Big) \,.
\label{BosonicLifshitzAction}
\end{eqnarray}
Here, the curvatures of the gauge fields are given by
\begin{eqnarray}
R_{\nu\rho}{}^a (P) &=& 2 \partial_{[\mu} e_{\nu]}{}^a + 2 \epsilon^{ab}\,  \omega_{[\mu} e_{\nu]b} - 2 \epsilon^{ab}\,  \omega_{ [\mu b} \tau_{\nu]}  - 2 b_{[\mu} e_{\nu]}^a \,,\nonumber\\
R_{\nu\rho}{}^a (G) &=&   2 \partial_{[\mu} \omega_{\nu]}{}^a + 2  \epsilon{}^{ab} \omega_{[\mu}{} \omega_{\nu] b} + 2 b_{[\mu} \omega_{\nu]}^a \,\nonumber\\
R_{\nu\rho} (H) &=&  2 \partial_{[\mu} \tau_{\nu]}  - 4 b_{[\mu} \tau_{\nu]}\,,\nonumber\\
R_{\nu\rho}(J) &=&  2 \partial_{[\mu} \omega_{\nu]} \,,\nonumber\\
R_{\nu\rho} (D) &=& 2 \partial_{[\mu} b_{\nu]}   \,. \label{LifshitzBosonicCurvatures}
\end{eqnarray}
As the non-relativistic special conformal symmetry is absent, the supersymmetric extension of the extended Lifshitz algebra does not require a special supersymmetry generator but the three fermionic generators $Q^\pm$ and $R$ that we introduced in the previous sections are sufficient. The inclusion of the generators $U_{1,2}$ is still necessary, but while $U_1$ acts non-trivially on the fermionic generators, $U_2$ is a central charge in the extended Lifshitz superalgebra. With the assistance of the computer algebra program \textit{Cadabra} \cite{Cadabra1,Cadabra2}, we found that the non-vanishing $\left[B,B\right]$ commutators of the extended Lifshitz superalgebra are given by (\ref{LifBoson}) while the non-vanishing $\left[B,F\right]$ commutators are
\begin{align}
[J, Q^\pm ] &= - \frac{1}{2} \gamma_0 Q^\pm  \,, & [J, R] &= - \frac{1}{2} \gamma_0 R  \,, &[G_a, Q^+  ] &= - \frac{1}{2} \gamma_a Q^-  \,,\nonumber\\
[G_a, Q^-  ] &= - \frac{1}{2} \gamma_a R \,, & [U_1, Q^\pm ] &= \pm  \gamma_0 Q^\pm  \,, &  [U_1, R ] &=    \gamma_0 R \,,\nonumber\\
[D, Q^+] &= - Q^+ \,, & [D, R] &=  R\,, & [S, Q^+ ] &= - \frac{1}{2} \gamma_0 R
\end{align}
Finally, the non-vanishing $\{F,F\}$ anti-commutators are given by
 \begin{eqnarray}
 \{ Q^+_\alpha ,  Q^+_\beta  \} &=& (\gamma_0 C^{-1})_{\alpha \beta} H \,,\nonumber\\
\{ Q^+_\alpha ,  Q^-_\beta  \}  &=& -(\gamma_a C^{-1})_{\alpha \beta} P^{a} \,,\nonumber\\
 \{ Q^-_\alpha ,  Q^-_\beta  \}  &=& (\gamma_0 C^{-1})_{\alpha \beta} M  - (\gamma_0 C^{-1})_{\alpha \beta} U_2\,, \nonumber \\
\{ Q^+_\alpha ,  R_\beta  \}  &=& (\gamma_0 C^{-1})_{\alpha \beta} M + 2 C_{\alpha \beta} Y + (\gamma_0 C^{-1})_{\alpha \beta} U_2 \,.
\label{FermionicLifshitz}
 \end{eqnarray}
The construction of a Chern-Simons action for the extended Lifshitz supergravity can be achieved by using the following non-degenerate invariant bi-linear form
 \begin{eqnarray}
 && (P_{a}, G_b) = \delta_{ab} \,, \quad (S,H) = -1 \,, \quad (M, J) = -1\,, \quad (Y,D) = -1 \nonumber\,,\\
 &&  (U_1, U_2) =  2\,,\quad  (Q^+_\alpha,R_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta}\,,\quad (Q^-_\alpha,Q^-_\beta)= 2 \left( C^{-1}\right)_{\alpha \beta} \,.
 \label{metricfermionicLifshitz}
 \end{eqnarray}
 Taking the gauge field to be
 \begin{eqnarray}
 A_\mu &=& \tau_\mu H + e_\mu{}^a P_a + \omega_\mu J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S +r_{1\mu} U_1 \nonumber\\
 && + r_{2\mu} U_2 + b_\mu D + y_\mu Y + \bar\psi_\mu^+ Q^+ + \bar\psi_\mu^- Q^- + \bar\rho_\mu R  \,,
 \label{GFSuperLif}
 \end{eqnarray}
we can construct a supersymmetric action
 \begin{eqnarray}
 S &=&\frac{k}{4 \pi} \int d^3x\, \varepsilon^{\mu\nu\rho}  \Big( \omega_\mu{}^a  R_{\nu\rho a} (P) + e_\mu{}^a  R_{\nu\rho a} (G)  - 2 s_\mu  R_{\nu\rho} (H) - 2 m_\mu  R_{\nu\rho}(J)   \nonumber\\
 && \quad  - 2 y_\mu  R_{\nu\rho} (D)   +   4 r_{2\mu} R_{\nu\rho}(U_1) +  4  \bar\rho_\mu R_{\nu\rho} (Q^+)  + 2 \bar\psi_\mu^- R_{\nu\rho} (Q^-)  \Big) \,.
 \label{FermionicLifshitzAction}
 \end{eqnarray}
 We refer to this action as the extended Lifshitz supergravity. The supercovariant curvatures that appears in the action are given by
 \begin{eqnarray}
 R_{\mu\nu}(H) &=& 2 \partial_{[\mu} \tau_{\nu]}  - 4 b_{[\mu} \tau_{\nu]} + \bar\psi_\mu^+ \gamma_0 \psi_\nu^+  \,,\nonumber\\
 R_{\mu\nu}{}^a (P) &=&  2 \partial_{[\mu} e_{\nu]}{}^a + 2 \epsilon^{ab}\,  \omega_{[\mu} e_{\nu]b} - 2 \epsilon^{ab}\,  \omega_{ [\mu b} \tau_{\nu]}  - 2 b_{[\mu} e_{\nu]}^a- 2 \bar\psi_{[\mu}^+ \gamma^a \psi_{\nu]}^-  \,,\nonumber\\
 R_{\mu\nu}{}^a (G) &=& 2 \partial_{[\mu} \omega_{\nu]}{}^a + 2  \epsilon{}^{ab} \omega_{[\mu}{} \omega_{\nu] b}  + 2 b_{[\mu} \omega_{\nu]}^a  \,,\nonumber\\
 R_{\mu\nu} (J) &=& 2 \partial_{[\mu} \omega_{\nu]}  \,,\nonumber\\
 R_{\mu\nu} (D) &=& 2 \partial_{[\mu} b_{\nu]}  \,, \nonumber \\
R_{\mu\nu}(U_1) &=& 2 \partial_{[\mu} r_{1\nu]}   \,,\nonumber\\
 R_{\mu\nu}(Q^+) &=& 2 \partial_{[\mu} \psi_{\nu]}^+ +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^+  - 2 b_{[\mu} \psi_{\nu]}^+  -  2 r_{1[\mu} \gamma_0 \psi_{\nu]}^+  \,,\nonumber\\
R_{\mu\nu}(Q^-) &=& 2 \partial_{[\mu} \psi_{\nu]}^- +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^-  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^+  +    2   r_{1[\mu} \gamma_0 \psi_{\nu]}^-  \,,\nonumber\\
R_{\mu\nu}(R) &=& 2 \partial_{[\mu} \rho_{\nu]} +  \omega_{[\mu} \gamma_0 \rho_{\nu]}  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^- +  s_{[\mu} \gamma_0 \psi_{\nu]}^+  -   2 r_{1[\mu} \gamma_0 \rho_{\nu]}  + 2 b_{[\mu}  \rho_{\nu]}  \,.
 \label{Scalesupercurvatures}
 \end{eqnarray}
 The extended Lifshitz supergravity is invariant under the following supersymmetry transformation rules
 \begin{eqnarray}
 \delta \tau_\mu &=&  - \bar{\epsilon}^+ \gamma_0 \psi^+_\mu \,,\nonumber\\
 \delta e_\mu{}^a &=&    \bar{\epsilon}^+ \gamma^{a} \psi_\mu^- +  \bar{\epsilon}^- \gamma^{a} \psi_\mu^+ \,,\nonumber\\
 \delta m_\mu &=& - \bar{\epsilon}^- \gamma_0 \psi_\mu^- -  \bar{\epsilon}^+ \gamma_0 \rho_\mu -  \bar{\eta} \gamma_0 \psi_\mu^+  \,,\nonumber\\
 \delta r_{2\mu}  &=&   \bar{\epsilon}^- \gamma_0 \psi_\mu^- -  \bar{\epsilon}^+ \gamma_0 \rho_\mu -  \bar{\eta} \gamma_0 \psi_\mu^+ \,,\nonumber\\
 \delta y_\mu  &=&  2 \bar{\eta} \psi_\mu^+ - 2 \bar{\epsilon}^+ \rho_\mu  \,,\nonumber \\
 \delta \psi_\mu^+ &=& \partial_\mu \epsilon^+ + \frac12 \omega_\mu \gamma_0 \epsilon^+  - b_\mu \epsilon^+ -  r_{1\mu} \gamma_0 \epsilon^+ \,,\nonumber\\
 \delta \psi_\mu^- &=&  \partial_\mu \epsilon^- + \frac12 \omega_\mu \gamma_0 \epsilon^- + \frac12 \omega_\mu{}^a \gamma_a \epsilon^+  +   r_{1\mu} \gamma_0 \epsilon^- \,,\nonumber\\
 \delta \rho_\mu &=&  \partial_\mu \eta + \frac12 \omega_\mu \gamma_0 \eta + \frac12 \omega_\mu{}^a \gamma_a \epsilon^- + \frac12 s_\mu \gamma_0 \epsilon^+  -  r_{1\mu} \gamma_0 \eta    + b_\mu \eta   \,.
 \label{GaugeLifshitzSusy}
 \end{eqnarray}
\subsection{Extended Schr\"odinger Supergravity}

As mentioned, our final purpose is to establish the supersymmetric extension of the Schr\"odinger algebra. The bosonic sector of the Schr\"odinger algebra extends the Lifshitz algebra (\ref{LifBoson}) with non-relativistic special conformal transformations $(K)$ as well as a second central charge $(Z)$. The non-zero commutation relations are given by \cite{Hartong EBG}
\begin{align}
\left[ H, G_a \right]  &= - \epsilon_{ab} P^b\,,          &  \left[ J, P_a \right] &= - \epsilon_{ab} P^b\,,          &  \left[ J, G_a \right] &=- \epsilon_{ab} G^b\,, \nonumber\\
\left[ P_{a}, G_b \right]  &= \epsilon_{ab} M + \delta_{ab} Y \,,        &  \left[ G_a, G_{b} \right]  &= \epsilon_{ab} S \,, & \left[ P_{a}, P_b \right]  &= \epsilon_{ab} Z \,,  \nonumber\\
\left[ D, G_a \right]  &= G_{a} \,,          &  \left[ D, P_a \right] &= -  P_{a} \,,          &  \left[ D, H \right] &=- 2 H\,, \nonumber\\
\left[ D, S \right]  &=2 S \,,        &  \left[ D,K \right]  &= 2K \,, & \left[ D,Z\right]  &= - 2 Z \,, \nonumber\\
\left[ K, P_{a} \right]  &= - \epsilon_{ab} G^{b} \,,          &  \left[ K, H\right] &= -  D \,,          &  \left[ K, Y \right] &= S \,, \nonumber\\
\left[ K, Z \right]  &=2 Y \,,        &  \left[ H, S \right]  &= - 2 Y \,, & \left[ H, Y \right]  &= - Z \,,
\label{SchBoson}
\end{align}
The corresponding  bi-linear form invariant under the extended Schr\"odinger algebra is given by \cite{Hartong EBG}
\begin{eqnarray}
&&(P_a,G_b) = \delta_{ab}\,, \quad (H,S) =-1 \,, \quad (J,M) =-1 \,, \quad  (D,Y) = -1 \,, \quad (K,Z) =  -1\,, \label{BosonicSchrodingerMetric}
\end{eqnarray}
Taking the gauge field to be
 \begin{eqnarray}
A_\mu &=& \tau_\mu H + e_\mu{}^a P_a + \omega_\mu J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S +  b_\mu D + f_\mu K + y_\mu Y +z_\mu Z  \,,
\label{GFBosSch}
\end{eqnarray}
we can construct the extended Schr\"odinger gravity action \cite{Hartong EBG}
\begin{eqnarray}
S &=&\frac{k}{4 \pi} \int d^3x\, \varepsilon^{\mu\nu\rho}  \Big( \omega_\mu{}^a R_{\nu\rho a} (P) - 2 s_\mu R_{\nu\rho} (H) - 2 m_\mu R_{\nu\rho}(J) + e_\mu{}^a R_{\nu\rho a} (G) \nonumber\\
&&  \quad  - 2 y_\mu R_{\nu\rho} (D) - 2 z_\mu R_{\nu\rho} (K)    \Big) \,.
\label{BosonicSchrodinger}
\end{eqnarray}
The supersymmetric extension of the extended Schr\"odinger algebra requires the inclusion of two new supercharges, $F^\pm$ that are both Majorana spinors. With the assistance of the computer algebra program \textit{Cadabra} \cite{Cadabra1,Cadabra2}, we found that the non-vanishing $\left[B,B\right]$ commutators of the extended Schr\"odinger superalgebra are given by (\ref{SchBoson}) while the non-vanishing $\left[B,F\right]$ commutators are
\begin{align}
[J, Q^\pm ] & = - \frac{1}{2} \gamma_0 Q^\pm \,, & [J, R ] &= - \frac{1}{2} \gamma_0 R  \,, & [J, F^\pm ] & = - \frac{1}{2} \gamma_0 F^\pm  \,,\nonumber\\
[S, Q^+ ] &= - \frac{1}{2} \gamma_0 R \,, & [G_a, Q^+  ] & = - \frac{1}{2} \gamma_a Q^- \,, & [G_a, Q^-  ] &= - \frac{1}{2} \gamma_a R \,,\nonumber\\
[P_a, F^+] &= -\frac{1}{2} \epsilon_{ab} \gamma^b Q^-\,, & [P_a,Q^-] &= -\frac{1}{2} \epsilon_{ab} \gamma^b F^- \,, & [U_1, Q^\pm ] &= \pm  \gamma_0 Q^\pm \,,\nonumber\\
[U_1, R ] &=    \gamma_0 R \,,  &[U_1, F^\pm] &=    \gamma_0 F^\pm \,, &[U_2, Q^+ ] &=  - \frac{3}{4}   F^- \,,\nonumber\\
[U_2, F^+ ] &=   \frac{3}{4}  R \,, & [H, F^+] &=  - Q^+ \,, & [H, R ] &=    F^- \,, \nonumber\\
[D, Q^+] &= - Q^+ \,, & [D, F^\pm] &=  \pm F^\pm \,, & [D, R] &=  R \,,\nonumber\\
[K,Q^+] &= F^+ \,,  & [K,F^-] &= - R \,, &[M,Q^+] &= \frac{1}{4} F^- \,,\nonumber\\
 [M,F^+] &=  - \frac{1}{4} R \,, &[Y,Q^+] &= \frac{1}{4} \gamma_0  F^- \,, & [Y,F^+] &=  \frac{1}{4} \gamma_0  R \,,\nonumber\\
 [Z,F^+] &= -  \frac{1}{2} \gamma_0   F^-  \,,
\end{align}
while the non-vanishing $\{F,F\}$ anti-commutators are given by
\begin{eqnarray}
&&\{ Q^+_\alpha ,  Q^+_\beta  \} = (\gamma_0 C^{-1})_{\alpha \beta} H \,,\nonumber\\
&&\{ Q^+_\alpha ,  Q^-_\beta  \}  = -(\gamma_a C^{-1})_{\alpha \beta} P^{a} ,,\nonumber\\
&&\{ Q^-_\alpha ,  Q^-_\beta  \}  = (\gamma_0 C^{-1})_{\alpha \beta} M  - (\gamma_0 C^{-1})_{\alpha \beta} U_2\,, \nonumber \\
&&\{ Q^+_\alpha ,  R_\beta  \}  = (\gamma_0 C^{-1})_{\alpha \beta} M + 2 C_{\alpha \beta} Y + (\gamma_0 C^{-1})_{\alpha \beta} U_2\,, \nonumber \\
&&\{ Q^+_\alpha ,  F^-_{\beta}  \}  = - 2 C_{\alpha \beta} Z \,, \nonumber \\
&&\{ F^+_\alpha ,  F^+_\beta  \} = (\gamma_0 C^{-1})_{\alpha \beta} K \,,\nonumber\\
&&\{ F^+_\alpha ,  Q^+_\beta  \} =  - \frac{1}{2}(\gamma_0 C^{-1})_{\alpha \beta} D + \frac{1}{2} C_{\alpha\beta} J  + \frac{3}{4} C_{\alpha \beta} U_1 \,,\nonumber\\
&&\{ F^+_\alpha ,  Q^-_\beta  \} =  \epsilon_{ab}(\gamma^a C^{-1})_{\alpha \beta} G^b\,, \nonumber \\
&&\{ F^+_\alpha ,  R_\beta  \} =  2 C_{\alpha \beta} S\,, \nonumber \\
&&\{ F^+_\alpha ,  F^-_{\beta}  \} =  (\gamma_0 C^{-1})_{\alpha \beta} M - 2  C_{\alpha \beta} Y + (\gamma_0 C^{-1})_{\alpha \beta} U_2 \,.
\label{3dSchrodingerExtendedBargmann}
\end{eqnarray}
The supersymmetric extension of the extended  Schr\"odinger algebra admits the following invariant bilinear forms
\begin{eqnarray}
&&(P_a,G_b) = \delta_{ab}\,, \quad (H,S) =- 1\,, \quad (J,M) =-1 \,, \quad  (D,Y) = -1 \,,  \quad (K,Z) =  -1\,, \nonumber \\
&&(R_1,R_2) = 2  \,, \quad (Q^+_\alpha,R_\beta)= 2 C_{\alpha \beta}\,,\quad (Q^-_\alpha,Q^-_\beta)= 2 C_{\alpha \beta}  \,,\quad (F^+_\alpha,F_\beta^-)= 2 C_{\alpha \beta}\,. \label{FermionicSchrodingerMetric}
\end{eqnarray}
Taking the gauge field to be
\begin{eqnarray}
A_\mu &=& \tau_\mu H + e_\mu{}^a P_a + \omega_\mu J + \omega_\mu{}^a G_a + m_\mu M + s_\mu S +  b_\mu D + f_\mu K + y_\mu Y +z_\mu Z \nonumber\\
&&+r_{1\mu} U_1 +r_{2\mu} U_2 + \bar\psi^+_\mu Q^+   + \bar\psi^-_\mu Q^-   + \bar\rho_\mu R +   \bar\phi^+_\mu F^+   + \bar\phi^-_\mu F^-  \,,
\label{GFSuperSch}
\end{eqnarray}
we obtain the following Chern-Simons action
\begin{eqnarray}
S &=&\frac{k}{4 \pi} \int d^3x\, \varepsilon^{\mu\nu\rho}  \Big( \omega_\mu{}^a  R_{\nu\rho a} (P) + e_\mu{}^a  R_{\nu\rho a} (G)  - 2 s_\mu  R_{\nu\rho} (H)    \nonumber\\
&& \quad - 2 m_\mu  R_{\nu\rho}(J) - 2 y_\mu  R_{\nu\rho} (D)  - 2 z_\mu R_{\nu\rho}(K) +   4 r_{2\mu} R_{\nu\rho}(U_1)  \nonumber\\
&& \quad +  4  \bar\rho_\mu R_{\nu\rho} (Q^+) + 2 \bar\psi_\mu^- R_{\nu\rho} (Q^-)  + 4 \bar\phi_\mu^- R_{\nu\rho} (F^+) \Big) \,.
\label{SuperSch}
\end{eqnarray}
We refer to this action as the extended Schr\"odinger supergravity. The Schr\"odinger supercovariant curvatures are given as
\begin{eqnarray}
R_{\mu\nu}(H) &=& 2 \partial_{[\mu} \tau_{\nu]}  - 4 b_{[\mu} \tau_{\nu]} + \bar\psi_\mu^+ \gamma_0 \psi_\nu^-   \,,\nonumber\\
R_{\mu\nu}{}^a (P) &=&  2 \partial_{[\mu} e_{\nu]}{}^a + 2 \epsilon^{ab}\,  \omega_{[\mu} e_{\nu]b} - 2 \epsilon^{ab}\,  \omega_{ [\mu b} \tau_{\nu]}  - 2 b_{[\mu} e_{\nu]}^a- 2 \bar\psi_{[\mu}^+ \gamma^a \psi_{\nu]}^-  \,,\nonumber\\
R_{\mu\nu}(M) &=& 2 \partial_{[\mu} m_{\nu]} + 2 \epsilon_{ab} \omega_{[\mu}{}^a e_{\nu]}{}^b + 2 \bar\psi_{[\mu}^+ \gamma_0 \rho_{\nu]} + \bar\psi_\mu^- \gamma_0 \psi_\nu^-  + 2 \bar\phi_{[\mu}^+ \gamma_0 \varphi_{\nu]} \,,\nonumber\\
R_{\mu\nu}{}^a (G) &=& 2 \partial_{[\mu} \omega_{\nu]}{}^a + 2  \epsilon{}^{ab} \omega_{[\mu}{} \omega_{\nu] b}  + 2  \epsilon{}^{ab} f_{[\mu}{} e_{\nu] b} + 2 b_{[\mu} \omega_{\nu]}^a  - 2\epsilon^{a b} \bar \phi_{[\mu}^+ \gamma_b \psi_{\nu]}^- \,,\nonumber\\
R_{\mu\nu} (J) &=& 2 \partial_{[\mu} \omega_{\nu]} +  \bar\phi_{[\mu}^+  \psi_{\nu]}^- \,,\nonumber\\
R_{\mu\nu} (S) &=& 2 \partial_{[\mu} s_{\nu]}  +  \epsilon{}^{ab} \omega_{[\mu a} \omega_{\nu]b }   +  4 b_{[\mu} s _{\nu]} + 2 f_{[\mu} y_{\nu]}+ 4 \bar\phi^+_{[\mu}  \rho_{\nu]} \,, \nonumber \\
R_{\mu\nu} (K) &=& 2 \partial_{[\mu} f_{\nu]}   +  4 b_{[\mu} f _{\nu]} +  \bar\phi^+_{\mu} \gamma_0  \phi^+_{\nu} \,, \nonumber \\
R_{\mu\nu} (D) &=& 2 \partial_{[\mu} b_{\nu]}   +  2 \tau_{[\mu} f _{\nu]}  -  \bar\phi_{[\mu}^+ \gamma_0 \psi_{\nu]}^+ \,, \nonumber \\
R_{\mu\nu} (Y) &=& 2 \partial_{[\mu} y_{\nu]}  -  2  \omega _{[\mu a} e_{\nu] }^a +  4 f_{[\mu}  z_{\nu]} - 4 \tau_{[\mu} s _{\nu]} +  4 \bar\psi_{[\mu}^+  \rho_{\nu]} - 4 \bar\phi_{[\mu}^+  \phi^-_{\nu]} \,, \nonumber \\
R_{\mu\nu} (Z) &=& 2 \partial_{[\mu} z_{\nu]}  +  \epsilon{}^{ab} e_{[\mu a} e_{\nu]b }  - 4  b_{[\mu} z_{\nu]} - 2 \tau_{[\mu} y _{\nu]} - 4 \bar\psi_{[\mu}^+  \varphi_{\nu]}\,, \nonumber \\
R_{\mu\nu}(U_1) &=& 2 \partial_{[\mu} r_{1\nu]}   +  \frac{3}{2}\bar\phi^+_{[\mu}  \psi_{\nu]}^+ \,,\nonumber\\
R_{\mu\nu}(U_2) &=& 2 \partial_{[\mu} r_{2\nu]}  +  2 \bar\psi_{[\mu}^+ \gamma_0 \rho_{\nu]} -  \bar\psi_\mu^- \gamma_0 \psi_\nu^-  +  2 \bar\phi_{[\mu}^+ \gamma_0 \varphi_{\nu]} \,,\nonumber \\
R_{\mu\nu}(Q^+) &=& 2 \partial_{[\mu} \psi_{\nu]}^+ +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^+ - 2 \tau_{[\mu} \phi^+_{\nu]}  - 2 b_{[\mu} \psi_{\nu]}^+  -  2 r_{1[\mu} \gamma_0 \psi_{\nu]}^+  \,,\nonumber\\
R_{\mu\nu}(Q^-) &=& 2 \partial_{[\mu} \psi_{\nu]}^- +  \omega_{[\mu} \gamma_0 \psi_{\nu]}^-  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^+  + \epsilon^{a b} e_{[\mu a} \gamma_b \phi_{\nu]}^+  +    2   r_{1[\mu} \gamma_0 \psi_{\nu]}^-  \,,\nonumber\\
R_{\mu\nu}(R) &=& 2 \partial_{[\mu} \rho_{\nu]} +  \omega_{[\mu} \gamma_0 \rho_{\nu]}  + \omega_{[\mu}{}^a \gamma_a \psi_{\nu]}^- +  s_{[\mu} \gamma_0 \psi_{\nu]}^+  -   2 r_{1[\mu} \gamma_0 \rho_{\nu]} - \frac{3}{2} r_{2[\mu}  \phi^+_{\nu]}  + 2 b_{[\mu}  \rho_{\nu]} \nonumber \\
&&  - 2 f_{[\mu}  \phi^-_{\nu]}  - \frac{1}{2} m_{[\mu} \phi^+_{\nu]} -   \frac{1}{2} y_{[\mu} \phi^+_{\nu]}  \nonumber \\
R_{\mu\nu}(F^+) &=& 2 \partial_{[\mu} \phi^+_{\nu]} +  \omega_{[\mu} \gamma_0 \phi^+_{\nu]}  -   2 r_{1[\mu} \gamma_0 \phi^+_{\nu]} + 2 b_{[\mu}  \phi^+_{\nu]}  +  2 f_{[\mu}  \psi_{\nu]}^+   \nonumber \\
R_{\mu\nu}(F^-) &=& 2 \partial_{[\mu} \phi^-_{\nu]} +  \omega_{[\mu} \gamma_0 \phi^-_{\nu]}   -   2 r_{1[\mu} \gamma_0 \phi^-_{\nu]}  - 2 b_{[\mu}  \phi^-_{\nu]} +  \frac{3}{2} r_{2[\mu}  \psi_{\nu]}^+  +  2 \tau_{[\mu}  \rho_{\nu]}   \nonumber \\
&&  + \frac{1}{2} m_{[\mu} \psi_{\nu]}^+ -   \frac{1}{2} y_{[\mu} \psi_{\nu]}^+  + z_{[\mu}  \gamma_ 0 \phi^+_{\nu]} \,. \nonumber  \label{Conformalsupercurvatures}
\end{eqnarray}
Finally, the extended Schr\"odinger supergravity action is invariant under the following transformation rules
\begin{eqnarray}
\delta \tau_\mu &=&  - \bar{\epsilon}^+ \gamma_0 \psi^+_\mu \,,\nonumber\\
\delta e_\mu{}^a &=&    \bar{\epsilon}^+ \gamma^{a} \psi_\mu^- +  \bar{\epsilon}^- \gamma^{a} \psi_\mu^+ \,,\nonumber\\
\delta \omega_\mu{}^a &=&  \epsilon^{ab} \bar{\epsilon}^- \gamma_{b} \phi^+_\mu +  \epsilon^{ab} \bar{\zeta}^+ \gamma_{b} \psi_\mu^- \,,\nonumber\\
\delta \omega_\mu &=& \frac{1}{2} \bar{\epsilon}^+  \phi^+_\mu -\frac{1}{2} \bar{\zeta}^+ \psi_\mu \,,\nonumber\\
\delta m_\mu &=& - \bar{\epsilon}^- \gamma_0 \psi_\mu^- -  \bar{\epsilon}^+ \gamma_0 \rho_\mu -  \bar{\eta} \gamma_0 \psi_\mu^+   -  \bar{\zeta}^+ \gamma_0 \phi^-_\mu -  \bar{\zeta}^- \gamma_0 \phi^+_\mu\,,\nonumber\\
\delta s_\mu  &=&  2 \bar{\eta} \phi^+_\mu - 2 \bar{\zeta}^+ \rho_\mu \,,\nonumber \\
\delta r_{1\mu}  &=&   \frac{3}{4}\bar{\epsilon}^+  \phi^+_\mu -  \frac{3}{4}\bar{\zeta}^+  \psi_\mu^+ \,, \nonumber \\
\delta r_{2\mu}  &=&   \bar{\epsilon}^- \gamma_0 \psi_\mu^- -  \bar{\epsilon}^+ \gamma_0 \rho_\mu -  \bar{\eta} \gamma_0 \psi_\mu^+   -  \bar{\zeta}^+ \gamma_0 \phi^-_\mu -  \bar{\zeta}^- \gamma_0 \phi^+_\mu\,,\nonumber\\
\delta b_\mu  &=&  \frac{1}{2} \bar{\epsilon}^+ \gamma_0 \phi^+_\mu  +  \frac{1}{2} \bar{\zeta}^+ \gamma_0 \psi_\mu^+\,,\nonumber \\
\delta f_\mu &=&  - \bar{\zeta}^+ \gamma_0 \phi^+_\mu \,,\nonumber\\
\delta y_\mu  &=&  2 \bar{\eta} \psi_\mu^+ - 2 \bar{\epsilon}^+ \rho_\mu -2 \bar{\zeta}^- \phi^+_\mu + 2\bar{\zeta}^+ \phi^-_\mu \,,\nonumber \\
\delta z_\mu  &=&  -2 \bar{\zeta}^- \psi_\mu^+ +  2 \bar{\epsilon}^+ \phi^-_\mu \,,\nonumber\\
\delta \psi_\mu^+ &=& \partial_\mu \epsilon^+ + \frac12 \omega_\mu \gamma_0 \epsilon^+ -  \tau_\mu {\zeta}^+  - b_\mu \epsilon^+ -  r_{1\mu} \gamma_0 \epsilon^+ \,,\nonumber\\
\delta \psi_\mu^- &=&  \partial_\mu \epsilon^- + \frac12 \omega_\mu \gamma_0 \epsilon^- + \frac12 \omega_\mu{}^a \gamma_a \epsilon^+  +  \frac{1}{2} \epsilon^{ab} e_{\mu a} \gamma_b {\zeta}^+  +   r_{1\mu} \gamma_0 \epsilon^- \,,\nonumber\\
\delta \rho_\mu &=&  \partial_\mu \eta + \frac12 \omega_\mu \gamma_0 \eta + \frac12 \omega_\mu{}^a \gamma_a \epsilon^- + \frac12 s_\mu \gamma_0 \epsilon^+  - \frac{1}{4}  m_\mu {\zeta}^+  \nonumber\\
&& -  r_{1\mu} \gamma_0 \eta  +  \frac{3}{4}  r_{2\mu} {\zeta}^+   + b_\mu \eta - f_\mu {\zeta}^- - \frac{1}{4} y_\mu \gamma_0 {\zeta}^+   \,, \nonumber\\
\delta \phi_\mu^+ &=& \partial_\mu \zeta^+ + \frac12 \omega_\mu \gamma_0 \zeta^+   +  b_\mu \zeta^+  + f_\mu \epsilon^+-  r_{1\mu} \gamma_0 \zeta^+ \,,\nonumber\\
\delta \phi^-_\mu &=&  \partial_\mu {\zeta}^- + \frac12 \omega_\mu \gamma_0 {\zeta}^- + \frac12 \epsilon_{ab} e_\mu{}^a \gamma^b \epsilon^-   + \frac{1}{4}  m_\mu \epsilon^+ +   r_{1\mu} \gamma_0 {\zeta}^- \nonumber\\
&& - \frac{3}{4}r_{2\mu} \epsilon^+  +  \frac{1}{2}  z_{\mu} \gamma_0 {\zeta}^+  -  b_\mu {\zeta}^- +  \tau_\mu \eta -  \frac{1}{4} y_\mu \gamma_0 \epsilon^+  \,,
\label{GaugeSchrodingerSusy}
\end{eqnarray}
where $\zeta^\pm$ are the the parameters of the local $F^\pm$ transformations, respectively. Due to the map between the Newton-Cartan geometry and Horava-Lifshitz gravity \cite{Horava}, our action (\ref{SuperSch}) correspond to a superconformal non-projectable Horava-Lifshitz gravity \cite{Hartong EBG}.

\section{Discussion}

In this paper, we establish the supersymmetric extension of the extended Newton-Hooke, Lifshitz and Schr\"odinger algebras and construct the corresponding Chern-Simons supergravity models. The extended Newton-Hooke superalgebra admits two distinct non-degenerate invariant bi-linear form that gives rise to two different supergravity models with the same equations of motion. These two models are particularly different in terms of the parity of the bosonic actions. In particular, we showed that there is an exotic non-relativistic model such that parity-even field equations arise from a parity-odd Lagrangian. We then showed that it is possible to improve the extended Bargmann superalgebra with dilatations (without including non-relativistic special conformal symmetry) which we called the extended Lifshitz superalgebra and also established the Chern-Simons extended Lifshitz supergravity action. In the final step, we include the non-relativistic special conformal symmetry and establish the extended Schr\"odinger superalgebra and the corresponding Chern-Simons extended Schr\"odinger supergravity action.

We consider our paper as a first step to construct an off-shell formulation for the extended Bargmann supergravity and its matter couplings. Therefore, it would be very useful to study supermultiplet representations of the extended Schr\"odinger superalgebra. In particular, based on the bosonic construction \cite{AfsharBergshoeff}, it is natural to expect that a multiplet with a complex scalar field $\Psi$ as the lowest element can be used to construct a super-Schr\"odinger invariant gravity model. Taking multiple number of such multiplets and gauge fixing the Schr\"odinger supergravity action would give rise to the supergravity coupling of such a multiplet to off-shell supergravity, which is an important step towards non-relativistic localization.

The $\mathcal{N}=2$ supergravity models that we established here can be extended to higher number of supercharges. In this case, the extended Bargmann superalgebra with $\mathcal{N} > 2$ can be obtained by Lie algebra expansion of the relativistic three-dimensional superalgebras. However, the Lifshitz and the Schr\"odinger superalgebras cannot be accessed by an Inonu-Wigner contraction or a Lie algebra expansion. In particular, given that the Schr\"odinger algebra is a particular instance $(\ell=1/2)$ of the so called $\ell$-conformal Galilei algebra \cite{Negro,Henkel}, in which case its $\mathcal{N} > 2$ supersymmetric extensions exist \cite{Galajinsky1,Galajinsky2}, it would be very nice to initiate a systematic study to establish extended Schr\"odinger superalgebras for $\mathcal{N}>2$ supersymmetry. There are also other three-dimensional non-relativistic algebras with an action principle such as the extended Newtonian supergravity \cite{Ozdemir}. Therefore, such a systematic study would also be helpful to uncover the Schr\"odinger extension of other three-dimensional superalgebras and hopefully towards a catalog of non-relativistic supergravity, at least in three dimensions. As a final remark, we would like to point out an interesting relation between the extended Schr\"odinger algebra and an extended Poincar\'e algebra. The extended Schr\"odinger algebra (and extended $\ell$-conformal Galilean algebras) has a hidden relativistic structure such that they can be written in a manifestly 3d  Lorentz-covariant form \cite{Sorokin}. It might of interest to see if the superextension of the Schr\"odinger superalgebra and the corresponding supergravity theory that we constructed in this paper can also be recast in a manifestly relativistic form.

\section*{Acknowledgements}

We thank Eric Bergshoeff for useful comments and suggestions. The authors would like to thank the referee for insightful comments which helped to improve this paper. The work of M.O. is supported in part by TUBITAK grant 118F091. N.O and U.Z. are supported in parts by Istanbul Technical University Research Fund under grant number TDK-2018-41133.