\subsection{Equivalence of first order PDE}\label{sec:1stpde}

Consider a general first order partial differential equation for a scalar function $u(x, y)$ of two variables,
\[
u_y=f(x, y, u, u_x).
\]
We shall study the equivalence of such equations under the pseudo-group, $\mathcal{G}$, of \emph{point transformations}.

\begin{Rem}
We shall not go into too much detail regarding the actual, nitty-gritty, computations in this example. A forthcoming paper is dedicated to the computational advantages of the involutive moving frame.
\end{Rem}

Writing $u_x=p$ and $u_y=q$, the elements of the pseudo-group are local diffeomorphisms
\[
(x, y, u, p, q)\overset{\varphi}{\mapsto} (X, Y, U, P, Q).
\]
Since the above is the first prolongation of $(x, y, u)\mapsto(X, Y, U)$, we have $X_p=X_q=Y_p=Y_q=U_p=U_q=0$ and
\begin{equation}\label{eq:pqmat}
\begin{bmatrix} P \\ Q \end{bmatrix}=\begin{bmatrix} X_x+pX_u & X_y+qX_y \\ Y_x+pY_u & Y_y+qY_u \end{bmatrix}^{-T}\begin{bmatrix} U_x+pU_u \\ U_y + qU_u\end{bmatrix}.
\end{equation}
These are the basic first order equations for $\mathcal{G}$ from which all other defining equations are found, by prolongation and projection. Note that differentiating both sides of (\ref{eq:pqmat}) w.r.t. $q$ will give determining equations $P_q=\cdots$ and $Q_q=\cdots$ where the right hand sides are first order, and we easily determine that this Lie pseudo-group is quasi horizontal of horizontal order 1.

Since the entries of the matrix above appear so frequently throughout we denote them by
\[
\begin{bmatrix} g_{11} & g_{12} \\ g_{21} & g_{22} \end{bmatrix}:=\begin{bmatrix} X_x+pX_u & X_y+qX_y \\ Y_x+pY_u & Y_y+qY_u \end{bmatrix}.
\]
We normalize $(X, Y, U, P, Q)$ to zero so that
\[
U_x=-pU_u,\quad U_y=-qU_u
\]
and the recurrence formula gives $\mu^u_X=-\o^p$ and $\mu^u_Y=-\sum \widehat{Q}_i\omega^i$. The recurrence formula for the first order lifted invariants $\widehat{Q}_i$ gives
\[
\aligned
d\widehat{Q}_X&=\sum \widehat{Q}_{Xi}\omega^i +\o^p \widehat{Q}_U-\widehat{Q}_P \mu^p_X+\widehat{Q}_P \widehat{Q}_X \mu^y_X+\widehat{Q}_X \mu^u_U-\widehat{Q}_X \mu^x_X-\widehat{Q}_X\mu^y_Y-\widehat{Q}_Y \mu^y_X+\mu^p_Y,\\
d\widehat{Q}_Y&=\sum \widehat{Q}_{Yi}\omega^i+\o^q \widehat{Q}_U-\widehat{Q}_P \mu^p_Y+\widehat{Q}_P \widehat{Q}_Y \mu^y_X+\widehat{Q}_Y\mu^u_U-\widehat{Q}_X \mu^x_Y-2 \widehat{Q}_Y \mu^y_Y+\mu^q_Y,\\
d\widehat{Q}_U&=\sum \widehat{Q}_{Ui}\omega^i-\widehat{Q}_P \mu^p_U+\widehat{Q}_P \widehat{Q}_U\mu^y_X-\widehat{Q}_X \mu^x_U-\widehat{Q}_U \mu^y_Y-\widehat{Q}_Y \mu^y_U+\mu^q_U,\\
d\widehat{Q}_P&=\sum \widehat{Q}_{Pi}\omega^i+ \mu^x_X+\widehat{Q}_P^2\mu^y_X-\widehat{Q}_P \mu^y_Y-\mu^x_Y,
\endaligned
\]
where, in the sums above, $i$ runs through the lifted horizontal base coordinates $X, Y, U, P$. We normalize $\widehat{Q}_i=0$ for $i=X, Y, U, P$ to obtain
\[
\aligned
\mu^p_Y&=-\sum \widehat{Q}_{Xi}\omega^i,\\
\mu^q_Y&=-\sum \widehat{Q}_{Yi}\omega^i,\\
\mu^q_U&=-\sum \widehat{Q}_{Ui}\omega^i,\\
\mu^x_Y&=-\sum \widehat{Q}_{Pi}\omega^i.
\endaligned
\]
Also, the recurrence formula $0=d\widehat{Q}=\sum \widehat{Q}_i\omega^i+\mu^u_Y$ now becomes $\mu^u_Y=0$.

The next step in the equivalence procedure is to compute the structure equations $d\omega^i$ on our equivariant moving frame $\widetilde{\mathcal{B}}_{1,1}$. Doing so we encounter as structure coefficients the lifted second order invariants $\widehat{Q}_{ij}$. We manage to normalize all of them except for $\widehat{Q}_{PP}$ to zero. At this point, there are $r^{(2)}=13$ second order group parameters we have \emph{not} solved for. These 13 parameters along with $X_x, X_u, Y_x, Y_y, Y_u, U_u, P_u, P_x$ and the second order jets $j^2q\restrict{(x,y,u,p)}$ parametrize the equivariant moving frame $\widetilde{\mathcal{B}}_{2,2}$.

Meanwhile,
\[
\widehat{Q}_{PP}=\frac{g_ {11}^2 q_{pp}}{U_u(g_{22}-g_{21} q_p)}
\]
and so this equivalence problem branches at this point, depending on whether $q_{pp}=0$ or not.

\subsubsection{Branch 1}
When $q_{pp}=0$ the original PDE has the form
\[
u_y=g(x, y, u)u_x+h(x, y,u).
\]
On this branch $\widehat{Q}_{PP}=0$ and there are no more normalizable structure coefficients in $d\omega^i$. We must then check for involutivity of the $\omega^i$. We have the structure equations
\[
\aligned
d\o^x&= -\o^x\wedge\mu^x_X-\o^u\wedge\mu^x_U,\\
d\o^u&= -\o^x\wedge\mu^y_X-\o^u\wedge\mu^y_Y-\o^u\wedge\mu^y_U,\\
d\o^u&= \o^x\wedge\o^p-\o^u\wedge\mu^u_U,\\
d\o^p&= -\o^x\wedge\mu^p_X-\o^u\wedge\mu^p_U-\o^p\wedge\mu^u_U+\o^p\wedge\mu^x_X.
\endaligned
\]
To perform Cartan's test for involutivity we compute the reduced Cartan characters $s^{(1)}_1, s^{(1)}_2, s^{(1)}_3, s^{(1)}_4$ for the above system. The first Cartan character $s^{(1)}_1$ is the maximum rank of the matrix
\[
\left(
\begin{array}{cccccccc}
 -a & 0 & 0 & -c & 0 & 0 & 0 & 0 \\
 0 & -a & -b & 0 & -c & 0 & 0 & 0 \\
 0 & 0 & 0 & 0 & 0 & -c & 0 & 0 \\
 d & 0 & 0 & 0 & 0 & -d & -a & -c \\
\end{array}
\right)
\]
where $a, b, c, d$ are arbitrary real numbers. This maximum rank is achieved for example when $a=c=1$ and $b=d=0$ and is $s^{(1)}_1=4$. Some simple computations reveal that $s^{(1)}_2=3$, $s^{(1)}_3=1$ and $s^{(1)}_4=0$. Now, Cartan's test for involution asks whether
\[
13=r^{(2)}=\sum_i i\cdot s^{(1)}_i=1\cdot 4+2\cdot 3+3\cdot1+4\cdot0=13.
\]
Since this is true, the frame $\{\o^x,\o^u,\o^u,\o^p\}$ is involutive and \emph{all} first order PDE of the form
\[
u_y=g(x, y, u)u_x+h(x, y,u)
\]
are equivalent under a point transformation and the general equivalence map depends on one arbitrary function of three variables, as per Cartan's theory of involution. A canonical form of a PDE in this class is the trivial $u_y=0$.

\subsubsection{Branch 2}
When $q_{pp}\neq0$ we normalize $U_u$ to obtain
\[
U_u= \frac{g_ {11}^2 q_{pp}}{g_{22}-g_{21} q_p}
\]
and the recurrence formula gives
\[
\mu^u_U=\sum \widehat{Q}_{PPi}\omega^i+2 \mu^x_X-\mu^y_Y.
\]
Going back to the structure equations $d\omega^i$ on $\widetilde{\mathcal{B}}_{1,2}$ and studying the structure coefficients, we find that
\[
d\widehat{Q}_{PPP}=\sum \widehat{Q}_{PPPi}\omega^i-\mu^x_X \widehat{Q}_{PPP}+\mu^y_Y \widehat{Q}_{PPP}+3 \mu^x_U+3 \mu^y_X
\]
and that we can normalize $\widehat{Q}_{PPP}=0$ and solve for $X_u$. The explicit formula for this normalization is
\[
X_u= \frac{g_{11} g_{21} q_{pp}}{g_{21} q_p-g_{22}}-\frac{g_{11} \
q_{ppp}}{3 q_{pp}}.
\]
At the next stage, computing $d\omega^i$ on $\widetilde{\mathcal{B}}_{1,3}$ we find the structure function $\widehat{Q}_{PPPP}$ and the recurrence formula gives
\[
d\widehat{Q}_{PPPP}=\sum \widehat{Q}_{PPPPi}\omega^i+\frac{2\widehat{Q}_{PPPP}}{3} (\mu^x_X -\mu^y_Y)-2 \mu^y_U.
\]
Therefor $Y_u$ can be normalized from $\widehat{Q}_{PPPP}=0$. The formula we obtain is
\[
Y_u= \frac{-9 g_ {21}^2 q_ {pp}^4+6 g_{21} q_{ppp} q_ {pp}^2 (g_{21} \
q_p-g_{22})-3 q_{pppp} q_{pp} (g_{22}-g_{21} q_p)^2+4 q_ {ppp}^2 \
(g_{22}-g_{21} q_p)^2}{18 q_ {pp}^3 (g_{22}-g_{21} q_p)}.
\]
At the next stage, computing $d\omega^i$ on $\widetilde{\mathcal{B}}_{1,4}$, we find $\widehat{Q}_{5P}$ and the recurrence formula tells us that
\[
d\widehat{Q}_{5P}=\sum \widehat{Q}_{5P,i}\omega^i-3\widehat{Q}_{5P}(\mu^x_X-\mu^y_Y)
\]
and some further computations reveal that
\[
\widehat{Q}_{PPPPP}=-\frac{(40 q_{ppp}^3-45 q_{pp} q_{pppp} q_{ppp}+9 q_{pp}^2 q_{ppppp}) (g_{22}-g_{21} q_p)^3}{54 g_{11}^3 q_{pp}^6}.
\]
This branch of the equivalence problem therefore branches again at this juncture depending on the vanishing or non-vanishing of the expression
\[
40 q_{ppp}^3-45 q_{pp} q_{pppp} q_{ppp}+9 q_{pp}^2 q_{ppppp}.
\]

In the case this expression is identically zero, for example for the PDE
\[
u_y=u_x^2,
\]
then, after we finish normalizing all possible structure coefficients from $d\omega^i$, we find that we have normalized all second order group parameters except for the $r^{(2)}=4$ parameters $P_{uu}, P_{ux}, X_{ux}$ and $P_{xx}$. We again check for involution. This time we have the structure equations
\[
\aligned
d\o^x&=-\o^x\wedge\mu^x_X-\o^u\wedge\o^p+\o^u\wedge\mu^y_X,\\
d\o^y&=-\o^x\wedge\mu^y_X-\o^u\wedge\mu^y_Y,\\
d\o^u&=\o^x\wedge\o^p-2 \o^u\wedge\mu^x_X+\o^u\wedge\mu^y_Y,\\
d\o^p&=-\o^x\wedge\mu^p_X-\o^u\wedge\mu^p_U-\o^p\wedge\mu^x_X+\o^p\wedge\mu^y_Y.
\endaligned
\]
This leads to reduced Cartan characters $4, 1, 0, 0$ and Cartans test gives
\[
r^{(2)}=4<1\cdot4+2\cdot1=6
\]
and so the coframe of $\omega^i$'s is not involutive. Prolonging to the next order we compute the structure equations $d\omega^i$ \emph{and} the exterior derivative of the not-yet normalized first order Maurer-Cartan forms $\mu^p_U,\mu^p_X,\mu^x_X,\mu^y_X,\mu^y_Y$. Normalizing all the lifted invariants appearing in this set of equations (notice however that the equation for $d\omega^i$ are unchanged) we normalize many third order group parameters and no new second order group parameters. After all is said and done we find that the ones we did not normalize are the $r^{(3)}=4$ third order group parameters
\[
P_{uuu}, P_{uux}, P_{uxx}, P_{xxx},
\]
while the structure equations are
\[
\aligned
d\mu^p_U&=\o^x\wedge\mu^p_{UX}+\o^u\wedge\mu^p_{UU}+\o^p\wedge\mu^x_{UX}+\mu^p_U\wedge\mu^x_X-\mu^p_X\wedge\mu^y_X,\\
d\mu^p_X&=\o^x\wedge\mu^p_{XX}+\o^u\wedge\mu^p_{UX}+\o^q\wedge\mu^x_{UX}+\o^p\wedge\mu^p_U+\mu^p_X\wedge\mu^y_Y,\\
d\mu^x_X&=\o^u\wedge\mu^x_{UX}+\o^x\wedge\mu^p_U-\o^u\wedge\mu^p_X,\\
d\mu^y_X&=-\o^x\wedge\mu^x_{UX}+\o^u\wedge\mu^p_U-\mu^x_X\wedge\mu^y_X-\mu^y_X\wedge\mu^y_Y,\\
d\mu^y_Y&=\o^x\wedge\mu^p_U-2\o^u\wedge\mu^p_X-\o^p\wedge\mu^y_X.
\endaligned
\]
This time the reduced Cartan characters are $4,0,0,0,0$ and Cartan's test for involution is satisfied: $r^{(3)}=4=1\cdot4 +0\cdot2+0\cdot3+0\cdot4+0\cdot5$. We say that the equivariant moving frame $\widetilde{\mathcal{B}}_{2,5}$ we have computed is \emph{involutive}. Notice that \emph{all} PDEs
\[
u_y=f(x,y,u,p)
\]
where $q_{pp}\neq0$ and
\[
40 q_{ppp}^3-45 q_{pp} q_{pppp} q_{ppp}+9 q_{pp}^2 q_{ppppp}=0
\]
are point-equivalent to $u_y=u_x^2$.

We stop here, but mention that for PDE such that neither $\widehat{Q}_{PP}$ nor $\widehat{Q}_{PPPP}$ vanish, we continue by normalizing $\widehat{Q}_{PPPP}=1$ and solving for $g_{11}$, recompute structure equations and so on.

\section*{Appendix: Eventually free Lie pseudo-groups}

As previously mentioned, the foundational work, \cite{OP-2008, OP-2009}, on the equivariant moving frame for pseudo-groups focused on pseudo-group actions that were \emph{eventually free}, meaning that at some order, the groupoid action, see (\ref{eq:hatw}), of $\H_q$ on some subset $\S_q\subset J^q(\mathcal{E})$ was free. Equivalently, each group parameter of order $\leq q$ can be normalized in the lifted invariants $\widehat{U}^\alpha_J$, $|J|\leq q$, restricted to $\S_q$. The key result in \cite{OP-2009}, dubbed \emph{persistence of freeness}, says that if the action is free at order $q$ in $\S_q$, it is also free at order $q+1$ on the pre-image $(\pi^{q+1}_q)^{-1}(\S_q)$ where $\pi^{q+1}_q:J^{q+1}(\mathcal{E})\to J^q(\mathcal{E})$ is the canonical projection. In this section we prove that eventually free actions are quasi-horizontal and that persistence of freeness is a trivial consequence thereof. We used persistence of freeness briefly leading up to the Lie-Tresse theorem above.

As a quick warm up for the general case, consider a Lie pseudo-group acting on graphs of maps $x\mapsto u(x)$, $x\in\mathbb{R}$, in $\mathbb{R}^2$. The total derivative on $J^\infty(\mathbb{R}^2)$ is
\[
D_x=\pder{x}+\sum u_{k+1}\pder{u_k},
\]
where $u_k$ is the $k\th$ $x$-derivative of u. The first prolongation of a map $\varphi(z)=\varphi(x,u)=(X,U)=Z$ is
\[
\widehat{U}_X=\frac{D_xU}{D_xX}=\frac{U_x+u_xU_u}{X_x+u_xX_u},
\]
but notice that the nominator and denominator are the $\displaystyle\frac{\partial}{\partial y}$ derivatives of $U$ and $X$ in flat coordinates for the graph of $x\mapsto u(x)$ since, on the graph of $u$,
\[
\displaystyle\frac{\partial}{\partial y}=\displaystyle\frac{\partial}{\partial x}+u_x\displaystyle\frac{\partial}{\partial u},
\]
(cf. (\ref{eq:dy})).

It is easily seen that in general, on the graph of $x\mapsto u(x)$, we have that
\begin{equation}\label{eq:dy}
\frac{d^k}{dy^k}Z=D_x^kZ.
\end{equation}
Defining
\[
D_{X}=\frac{1}{D_xX}D_x,
\]
the higher order prolongations of $\varphi$ are
\[
\widehat{U}_k=D_{X}U,
\]
which are functions of the $y$-derivatives of $\varphi$ only. In arbitrary many variables, $n+m$, for $\mathcal{E}=\mathcal{X}\times\mathcal{U}$, the formulas for $\widehat{U}_J$, restricted to a graph of $x\mapsto u(x)$, become rather complicated but they too are functions of the expressions
\[
\frac{\partial^{|K|}}{\partial y^K}Z^a=D^KZ^a,\quad K\in\mathbb{Z}_{\geq0}^n,
\]
where
\[
\frac{\partial}{\partial y^i}=\frac{\partial}{\partial x^i}+u\av_i\frac{\partial}{\partial u\av}.
\]

Now assume the action of $\H$ on $\mathcal{E}$ is free at order $q$ in $\mathcal{V}_q\subset J^q(\mathcal{E})$ and $j^qu\restrict{x} \in\mathcal{V}_q$. Let $x\mapsto u(x)$ be a representing map with $q$-jet $j^qu\restrict{x}$ at $x$ and such that the entire $q$ jet $j^qu$ lies in $\mathcal{V}_q$. Freeness means that it is possible to solve for every group parameter of order $\leq q$ in a set of equations (restricted to $x\mapsto u(x)$) of the form
\begin{equation}\label{eq:ucon}
\widehat{U}^{\alpha_i}_{J^i}=c^{\alpha_i}_{J^i}=\text{constant},
\end{equation}
for $1\leq i\leq \text{dim}\H_q$. We shall prove that $\H$ is necessarily quasi-horizontal at some order $\leq q$. Writing the determining equations in flat coordinates $(y,w)$ w.r.t $x\mapsto u(x)$, we write $\pmb{Y}^{(q)}$ for those pseudo-group jet coordinates $Z^a_K$, $|K|\leq q$, that are pure $y$-derivatives and $\pmb{W}^{(q)}$ for at most $q\th$ order pseudo-group jet coordinates that are not. We consider two sets of equations; the determining equations for $\H$ of order $\leq q$ and (\ref{eq:ucon}), in flat coordinates for $u$. These can be written
\begin{align*}
F(z,Z, \pmb{Y}^{(q)}, \pmb{W}^{(q)})&=0,\\
G(z, Z, \pmb{Y}^{(q)})&=0,
\end{align*}
since (\ref{eq:ucon}) only involve $\pmb{Y}^{(q)}$. The fact that this system can be uniquely solved for \emph{all} $\pmb{Y}^{(q)}$ and $\pmb{W}^{(q)}$ directly implies that
\[
F(z, Z, \pmb{Y}^{(q)}, \pmb{W}^{(q)})=0
\]
can be solved for $\pmb{W}^{(q)}$. In particular, \emph{each} pseudo-group jet involving a $w$-derivative can be taken as principal. Thus, eventually free Lie pseudo-group actions are quasi-horizontal and they have horizontal order equal to the order of the pseudo-group. In addition, each pseudo-group jet involving a $u^\alpha$-derivative can be taken as principal.

\begin{theorem}\label{thm:efqh}
Let $\H$ be an eventually free Lie pseudo-group, of order $t^*$, acting on the bundle $\mathcal{E}\to\mathcal{X}$. Then $\H$ is quasi-horizontal of horizontal order $t^*$ and all pseudo-group parameters can be taken to be pure $\mathcal{X}$ derivatives.
\end{theorem}

\begin{Ex}\label{ex:stab}
Consider the Lie group action
\[
(x, u)\mapsto (\lambda x + a, \lambda^{k+1}u + P_k(x)),\quad \lambda>0,
\]
where $P_k(x)$ is an arbitrary polynomial of order $k$. The determining equations for this pseudo group are
\[
X_u=0,~~X_{xx}=0,~~ U_u=X_x^{k+1},~~ \frac{\partial^{k+1}}{\partial x^{k+1}}U=0.
\]
Notice that, since Lie group actions are eventually free almost everywhere, in accordance with the above theorem, \emph{each} pseudo-group jet involving a $u$-derivative can be taken as principal. The order of the pseudo-group is $k+1$.
\end{Ex}

We can now present a rather trivial proof of \emph{persistence of freeness} of Lie pseudo-groups. Let $j^qu\restrict{x}\in\mathcal{V}_q\subset J^q(\mathcal{E})$ lie in a set $\mathcal{V}_q$ on which $\H_q$ acts freely, as above. In flat coordinates w.r.t. a representing map $x\mapsto u(x)$, each pseudo-group jet involving a $w$ derivative can be taken as principal. Let us denote the pseudo-group jet coordinates in $\pmb{Y}^{(q)}$ that have order equal to $q$ by $\pmb{Y}^q$. Then in the determining equations $F^{(q)}(z, Z^{(q)})=0$, ignoring those equations with principal derivatives having a $w$-derivative and those of order less than $q$, we have some subset of $q\th$ order determining equations $\bar{F}^{(q)}(z, \pmb{Y}^{q})=0$ that only involve pure $y$-derivatives of order $q$. Then the $y$-prolongations
\[
\left(\partial_{y_1}\bar{F}^{(q)}\right)(z, \pmb{Y}^{q+1})=0,\ldots, \left(\partial_{y_n}\bar{F}^{(q)}\right)(z, \pmb{Y}^{q+1})=0
\]
give a system $\bar{F}^{(q+1)}(z, \pmb{Y}^{q+1})=0$ of determining equations of order $q+1$ that involve only pseudo-group jets of pure $y$-derivatives.

By freeness, on $\mathcal{V}_q$, there is a collection $\widehat{U}^{\alpha_i}_{J^i}$, $1\leq i\leq (\text{dim}\H_q - \text{dim}\H_{q-1})$, of lifted invariants with $|J^i|=q$, from which all $q\th$ order parametric derivatives can be normalized. That is, the collection $\widehat{U}^{\alpha_i}_{J^i}$, $1\leq i\leq (\text{dim}\H_q - \text{dim}\H_{q-1})$, is full rank in the $q\th$ order parametric derivatives (all of which are pure $y$-derivatives since $\H$ is eventually free). Hence, the set of functions on $\widetilde{\H}^q$,
\begin{equation}\label{eq:fuh}
\left\{\bar{F}^{(q)}(z,  \pmb{Y}^{q}), \widehat{U}^{\alpha_i}_{J^i}\right\}\restrict{1\leq i\leq (\text{dim}\H_q - \text{dim}\H_{q-1})}
\end{equation}
is full rank in the pseudo-group jets $\pmb{Y}^{q}$ when evaluated at jets in $\mathcal{V}_q$. Now, for any lifted invariant $\widehat{U}^\alpha_J$, the two collections (once evaluated at $j^qu\restrict{x}$)
\begin{equation}\label{eq:2col}
\aligned
\{D_{X^i}\widehat{U}^\alpha_J\}_{1\leq i\leq n},\quad \{\partial_{y_i}\widehat{U}^\alpha_J\}_{1\leq i\leq n}
\endaligned
\end{equation}
are invertible linear combinations of each other with coefficients depending on first order jets. Therefore the two collections (\ref{eq:2col}) are equal in rank in the top order pseudo-group jets. Going back to (\ref{eq:fuh}), this means that prolonging all functions in (\ref{eq:fuh}) w.r.t $\partial_{y_i}$ for $1\leq i\leq n$ will give a collection having the same rank in the $\pmb{Y}^{q+1}$ as the system
\begin{equation}\label{eq:fuh2}
\left\{\bar{F}^{(q+1)}(z,  \pmb{Y}^{q+1}), \widehat{U}^{\alpha_i}_{J^i, j}\right\}\restrict{1\leq i\leq (\text{dim}\H_q - \text{dim}\H_{q-1}), 1\leq j\leq n}.
\end{equation}
But since (\ref{eq:fuh}) was full rank in $\pmb{Y}^q$, its prolongation w.r.t. the $y$ variables will be full rank in $\pmb{Y}^{q+1}$. Since (\ref{eq:fuh2}) has the same rank in $\pmb{Y}^{q+1}$ we have proved persistence of freeness.

\section{Appendix: Eventually free Lie pseudo-groups}

As previously mentioned, the foundational work, \cite{OP-2008, OP-2009}, on the equivariant moving frame for pseudo-groups focused on pseudo-group actions that were \emph{eventually free}, meaning that at some order, the groupoid action, see (\ref{eq:hatw}), of $\H_q$ on some subset $\mathcal{V}_q\subset J^q(\mathcal{E})$ was free. Equivalently, each group parameter of order $\leq q$ can be normalized in the lifted invariants $\widehat{U}^\alpha_J$, $|J|\leq q$, restricted to $\mathcal{V}_q$. The key result in \cite{OP-2009}, dubbed \emph{persistence of freeness}, says that if the action is free at order $q$ in $\mathcal{V}_q$, it is also free at order $q+1$ on the pre-image $(\pi^{q+1}_q)^{-1}(\mathcal{V}_q)$ where $\pi^{q+1}_q:J^{q+1}(\mathcal{E})\to J^q(\mathcal{E})$ is the canonical projection. In this section we prove that eventually free actions are quasi-horizontal and that persistence of freeness is a trivial consequence thereof.

As a quick warm up for the general case, consider a Lie pseudo-group acting on graphs of maps $x\mapsto u(x)$, $x\in\mathbb{R}$, in $\mathbb{R}^2$. The total derivative on $J^\infty(\mathbb{R}^2\to\mathbb{R})$ is
\[
D_x=\pder{x}+\sum u_{k+1}\pder{u_k},
\]
where $u_k$ is the $k\th$ $x$-derivative of u. The first prolongation of a map $\varphi(z)=\varphi(x,u)=(X,U)=Z$ is
\[
\widehat{U}_X=\frac{D_xU}{D_xX}=\frac{U_x+u_xU_u}{X_x+u_xX_u},
\]
but notice that the nominator and denominator are the $\displaystyle\frac{\partial}{\partial y}$ derivatives of $U$ and $X$ in flat coordinates for the graph of $x\mapsto u(x)$ since, on the graph of $u$,
\[
\displaystyle\frac{\partial}{\partial y}=\displaystyle\frac{\partial}{\partial x}+u_x\displaystyle\frac{\partial}{\partial u}.
\]

It is easily seen that in general, on the graph of $x\mapsto u(x)$, we have that
\begin{equation}\label{eq:dy}
\frac{d^k}{dy^k}\varphi=D_x^k\varphi.
\end{equation}
Defining
\[
D_{\mathrm{X}}=\frac{1}{D_xX}D_x,
\]
the higher order prolongations of $\varphi$ are
\[
\widehat{U}_k=D_{\mathrm{X}}^kU,
\]
which are functions of the $y$-derivatives of $\varphi$ only. In arbitrary many variables, $n+m$, for $\mathcal{E}=\mathcal{X}\times\mathcal{U}$, and a local diffeomorphism $\varphi(x,u)=(X,U)$, the formulas for $\widehat{U}_J$, restricted to a graph of $x\mapsto u(x)$, become rather complicated but they too are functions of the expressions
\[
\frac{\partial^{|K|}}{\partial y^K}\varphi^a=D^K\varphi^a,\quad K\in\mathbb{Z}_{\geq0}^n,
\]
where
\[
\frac{\partial}{\partial y^i}=\frac{\partial}{\partial x^i}+u\av_i\frac{\partial}{\partial u\av}.
\]

Now assume the action of $\H$ on $\mathcal{E}$ is free at order $q$ in $\mathcal{V}_q\subset J^q(\mathcal{E})$ and $j^qu\restrict{x} \in\mathcal{V}_q$. Let $x\mapsto u(x)$ be a representing map with $q$-jet $j^qu\restrict{x}$ at $x$ and such that the entire $q$ jet $j^qu$ lies in $\mathcal{V}_q$. Freeness means that it is possible to solve for every group parameter in a set of equations (restricted to $x\mapsto u(x)$) of the form
\begin{equation}\label{eq:ucon}
\widehat{U}^{\alpha_r}_{J^r}=c^{\alpha_r}_{J^r}=\text{constant},
\end{equation}
for $1\leq r\leq \delta(q)$, where $\delta(q)$ is equal to the number of parametric derivatives in $\H_q$. Write the determining equations in flat $(y,w)$ coordinates w.r.t. $x\mapsto u(x)$. We shall prove that $\H$ is necessarily quasi-horizontal. We write $\pmb{Y}^{(q)}$ for those pseudo-group jet coordinates $Z^a_K$, $|K|\leq q$, that are pure $y$-derivatives and $\pmb{W}^{(q)}$ for those $Z^a_K$, $|K|\leq q$ that are nor. We consider two sets of equations; the determining equations for $\H$ of order $\leq q$ and (\ref{eq:ucon}), in flat coordinates for $u$. These can be written
\begin{align*}
F(z,Z, \pmb{Y}^{(q)}, \pmb{W}^{(q)})&=0,\\
G(z, Z, \pmb{Y}^{(q)})&=0,
\end{align*}
since (\ref{eq:ucon}) only involve $\pmb{Y}^{(q)}$. The fact that this system can be uniquely solved for \emph{all} $\pmb{Y}^{(q)}$ and $\pmb{W}^{(q)}$ directly implies that
\[
F(z, Z, \pmb{Y}^{(q)}, \pmb{W}^{(q)})=0
\]
can be solved for $\pmb{W}^{(q)}$. In particular, \emph{each} pseudo-group jet involving a $w$-derivative can be taken as principal. Thus, eventually free Lie pseudo-group actions are quasi-horizontal and they have horizontal order equal to the order of the pseudo-group. In addition, each pseudo-group jet involving a $u^\alpha$-derivative can be taken as principal. The free pseudo-groups have therefore are rather special structure.

\begin{theorem}\label{thm:efqh}
Let $\H$ be a $t^*$ order, eventually free Lie pseudo-group acting on the bundle $\mathcal{E}\to\mathcal{X}$. Then $\H$ is quasi-horizontal and all pseudo-group parameters can be taken to be pure $\mathcal{X}$ derivatives.
\end{theorem}

\begin{Ex}\label{ex:stab}
Consider the Lie group action
\[
(x, u)\mapsto (\lambda x + a, \lambda^{k+1}u + P_k(x)),\quad \lambda>0,
\]
where $P_k(x)$ is an arbitrary polynomial of order $k$. The determining equations for this pseudo group are
\[
X_u=0,~~X_{xx}=0,~~ U_u=X_x^{k+1},~~ \frac{\partial^{k+1}}{\partial x^{k+1}}U=0.
\]
Notice that, in accordance with the above theorem, \emph{each} pseudo-group jet involving a $u$-derivative can be taken as principal. The order of the pseudo-group is $k+1$.
\end{Ex}

We can now present a rather trivial proof of \emph{persistence of freeness} of Lie pseudo-groups. Let $j^qu\restrict{x}\in\mathcal{V}_q\subset J^q(\mathcal{E})$ lie in a set $\mathcal{V}_q$ on which $\H_q$ acts freely, as above. In flat coordinates w.r.t. a representing map $x\mapsto u(x)$, each pseudo-group jet involving a $w$ derivative can be taken as principal. Let us denote the pseudo-group jets coordinates in $\pmb{Y}^{(q)}$ that have order equal to $q$ by $\pmb{Y}^q$. Then in the determining equations $F^{(q)}(z, Z^{(q)})=0$, ignoring those equations with principal derivatives having a $w$-derivatives and those of order less than $q$, we have some subset of $q\th$ order determining equations $\bar{F}^{(q)}(z, \pmb{Y}^{q})=0$ that only involve pure $y$-derivatives of order $q$. Then the $y$-prolongations
\[
\left(\partial_{y_1}\bar{F}^{(q)}\right)(z, \pmb{Y}^{q+1})=0,\ldots, \left(\partial_{y_n}\bar{F}^{(q)}\right)(z, \pmb{Y}^{q+1})=0
\]
give a system $\bar{F}^{(q+1)}(z, \pmb{Y}^{q+1})=0$ of determining equations of order $q+1$ that involve only pseudo-group jets of pure $y$-derivatives.

By freeness, on $\mathcal{V}_q$, there is a collection $\widehat{U}^{\alpha_s}_{J^s}$, $1\leq s\leq \delta^*(q)$, of lifted invariants with $|J^s|=q$, from which all $q\th$ order parametric derivatives can be normalized, where $\delta^*(q)$ is the number of pseudo-group parameters of order exactly $q$. That is, the collection $\widehat{U}^{\alpha_s}_{J^s}$, $1\leq s\leq \delta^*(q)$, is full rank in the $q\th$ order parametric derivatives (all of which are pure $y$-derivatives since $\H$ is eventually free). Hence, the set of functions on $\widetilde{\H}^q$,
\begin{equation}\label{eq:fuh}
\left\{\bar{F}^{(q)}(z,  \pmb{Y}^{q}), \widehat{U}^{\alpha_s}_{J^s}\right\}\restrict{1\leq i\leq \delta^*(q)}
\end{equation}
is full rank in the pseudo-group jets $\pmb{Y}^{q}$ when evaluated at jets in $\mathcal{V}_q$. Now, for any lifted invariant $\widehat{U}^\alpha_J$, the two collections (once evaluated at $j^qu\restrict{x}$)
\begin{equation}\label{eq:2col}
\aligned
\{D_{X^i}\widehat{U}^\alpha_J\}_{1\leq i\leq n},\quad \{\partial_{y_i}\widehat{U}^\alpha_J\}_{1\leq i\leq n}
\endaligned
\end{equation}
are invertible linear combinations of each other with coefficients depending on first order jets. Therefore the two collections (\ref{eq:2col}) are equal in rank in the top order pseudo-group jets. Going back to (\ref{eq:fuh}), this means that prolonging all functions in (\ref{eq:fuh}) w.r.t $\partial_{y_i}$ for $1\leq i\leq n$ will give a collection having the same rank in the $\pmb{Y}^{q+1}$ as the system
\begin{equation}\label{eq:fuh2}
\left\{\bar{F}^{(q+1)}(z,  \pmb{Y}^{q+1}), \widehat{U}^{\alpha_s}_{J^s, j}\right\}\restrict{1\leq i\leq \delta^*(q), 1\leq j\leq n}.
\end{equation}
But since (\ref{eq:fuh}) was full rank in $\pmb{Y}^q$, its prolongation w.r.t. the $y$ variables will be full rank in $\pmb{Y}^{q+1}$. Since (\ref{eq:fuh2}) has the same rank in $\pmb{Y}^{q+1}$ we have proved persistence of freeness.

\section{Introduction}
This paper continues the project, begun in \cite{IMF}, of relating Cartan's equivalence method to the equivariant moving frame for pseudo-groups. The major step was accomplished in \cite[Theorem 3.12]{IMF} where Cartan's classical solution to the congruence problem for submanifolds of finite dimensional Lie groups was extended to the groupoids associated to Lie pseudo-groups. Cartan's solution in the finite dimensional case underlies his classical method of moving frames as well as the modern equivariant moving frame, \cite{FO-1999} and its infinite dimensional analog harmonized the finite and infinite dimensional methods of equivalence. However, the key theorem \cite[Theorem 3.12]{IMF} was most directly applicable to so-called horizontal actions which left out some important equivalence problems, particularly those for partial differential equations. The current paper completes this picture by clarifying precisely which equivalence problems can be solved using some version of \cite[Theorem 3.12]{IMF} and how the theorem applies.

There are two major consequences of this, both of which indicate that our framework is fruitful for equivalence problems for Lie pseudo-groups. First, we obtain a new, constructive and intuitive proof of the Lie-Tresse theorem on the finiteness of generators of the differential invariant algebra of Lie pseudo-groups. Our proof is much simpler than previous modern proofs, \cite{KL-2006, KL-2016, OP-2009}, and, interestingly, makes critical use of the finite dimensional theory. Second, there are enormous computational advantages to the combination of the equivariant moving frame and Cartan's geometry. We study, in Section \ref{sec:1stpde}, the general equivalence problem of first order PDE
\[
u_y=f(x, y, u, u_x)
\]
under point transformations and characterize those equations that are point-equivalent to $u_y=0$ and $u_y=u_x^2$. The branching found here are new. A forthcoming paper is dedicated to the computational advantages of our framework.

The paper's setup is as follows. In Section \ref{sec:LG} we give a quick review of Lie pseudo-groups and equivariant moving frames and refine the definition of a partial moving frame from \cite{IMF}. We explain how the key results from \cite{IMF} carry over to the class of ``vertically integrable'' pseudo groups, leading to the equivalence method in Theorem \ref{thm:total}. In the following section we nail down the precise condition on pseudo-groups for which Theorem \ref{thm:total} is most natural and we demonstrate this application by studying first order PDE under point-transformations. In the final section we give a new, constructive and intuitive proof of the Lie-Tresse theorem. This version of the theorem also gives general upper bounds on the minimal number of generating invariants and we consider a few examples of this bound. In the Appendix we consider the fundamental concept from \cite{OP-2009} of \emph{persistence of freeness} of Lie pseudo-groups. This property of pseudo-groups was proven in \cite{OP-2009} by rather complex means, but we show that it becomes essentially trivial in our framework.

\subsection{Upper bounds on the size of the generating set}

Theorem \ref{thm:LieTresse} provides an upper bound on the number of generating invariants. We have that
\begin{equation}\label{genset}
\{I^a_{A; l}\}_{1\leq l\leq n,~ 1\leq a\leq n+m,~ |A|<r}
\end{equation}
is a set of generators, where
\[
\mu^a_A\mapsto I^a_{A; l}\o^l.
\]
Since those $\mu^a_A$ corresponding to principal derivatives in the determining equations for $\H$ are constant-coefficient linear combinations of parametric ones, we can focus only on the latter. Therefore, to count the number of $I$'s we simply have to count the number of parametric derivatives of order $<r$ in $\H$ and multiply by $n$. The former number is equal to the dimension of the groupoid $\H_{r-1}$. So a general upper bound for the \emph{minimal number of generating invariants} is
\[
n\cdot \left(\text{dim}\H_{r-1}\right).
\]
In the case that $\H$ is a finite dimensional Lie group its determining equations are maximally over-determined and so the dimension of the Lie group is equal to $\text{dim}\H_{r-1}$. We therefore have, for a Lie group $G$ acting freely on $\S_q$ the upper bound
\[
n\cdot \left(\text{dim}G\right).
\]
\begin{Rem}
This upper bound, in the Lie group case, has also be deduced using the equivariant moving frame, \cite{FO-1999}. There, the members of the generating set (\ref{genset}) were called \emph{Maurer-Cartan invariants}.
\end{Rem}

Further assuming that $\H$ acts transitively on $\mathcal{E}$ so that normalizing the $X$'s gives $\mu^i=-\o^i$ and so there are $n$ Maurer-Cartan forms that do not provide any invariants we obtain the upper bound
\begin{equation}\label{upperb2}
n\cdot \left( \text{dim}\H_{r-1}-n\right)
\end{equation}
on the number of generating invariants. Similarly, if $\H$ acts transitively on $J^1(\mathcal{E})$, normalizing all the first order (and zero order) lifted invariants $\widehat{U}^\alpha_i\equiv e^\alpha_i$ to constants $e^\alpha_i$, the recurrence formula gives
\[
0=d\widehat{U}^\alpha=\widehat{U}^\alpha_i\omega^i+\text{Maurer-Cartan forms}=e^\alpha_i\omega^i+\text{Maurer-Cartan forms}
\]
By transitivity we can solve these equations for $m$ Maurer-Cartan forms which then become \emph{constant coefficient} combinations of $\o$'s. We can hence subtract a further $n\cdot m$ from our upper bound (\ref{upperb2}). More generally, using the same argument, if $\H_{r}$ acts transitively on $J^{q^*}(\mathcal{E})$ we have the upper bound
\[
n\cdot \left(\text{dim}\H_{r-1}-\text{dim}J^{q^*-1}(\mathcal{E})\right).
\]
Now, $\text{dim}J^{q^*-1}=n+m\binom{q^*-1+n}{n}$ and so the above bound is equal to
\begin{equation}\label{upperb3}
n\cdot \left(\text{dim}\H_{r-1}-n-m\binom{q^*-1+n}{n}\right).
\end{equation}
For a free Lie-group action, we replace $ \text{dim}\H_{r-1}$ with the dimension of the Lie group.

\begin{Ex}
Consider the Lie group action
\[
(x, u)\mapsto (X, U)= (\lambda x + a, \lambda u + b),\quad \lambda>0.
\]

The defining equations are
\[
X_u=X_{xx}=U_x=U_{uu}=0,\quad U_u=X_x
\]
and this pseudo-group is quasi horizontal of order $2$. By Theorem \ref{thm:LieTresse}, a second order moving frame will provide us with a generating set of invariants. Normalizing $X=U=0$ the recurrence formula gives
\[
\mu^x=-\o^x=-X_xdx,\quad \mu^u=-\widehat{U}_X\o^x
\]
where $\widehat{U}_X=u_x$ is an invariant. Moving up to the next order, we have
\[
du_x=\widehat{U}_{XX}\o^x=\frac{u_{xx}}{X_x}\o^x
\]
and so (assuming $u_{xx}>0$) we normalize $\widehat{U}_{XX}=1$ to obtain
\[
X_x=u_{xx}
\]
and
\[
0=d\widehat{U}_{XX}=\widehat{U}_{XXX}\o^x-\mu^x_X.
\]
Solving for $\mu^x_X$ we find that
\[
\mu^x_X=\frac{u_{xxx}}{u_{xx}^2}\o^x.
\]
By Theorem \ref{thm:LieTresse}, $u_x$ and $\displaystyle \frac{u_{xxx}}{u_{xx}^2}$ generate the invariant differential algebra. Notice that $u_x$ alone is not a generator, since
\[
du_x=\o^x
\]
and so applying the invariant derivative to $u_x$ generates no new invariants. Hence, we have a minimal set of generators of size 2.

Now, the upper bound (\ref{upperb2}) is
\[
1\cdot\left(3-1\right)=2.
\]
Therefore, the upper bound (\ref{upperb2}) is actually optimal in complete generality.
\end{Ex}

Although, as the last example showed, the upper bound (\ref{upperb2}) is optimal it is far from robust; in most cases it will overshoot the actual size of a minimal generating set by quite a lot.

\begin{Ex}
Consider the Lie-group action
\[
(x, u)\mapsto (x+a, u+p(x))
\]
where $x, a\in\mathbb{R}^n$, $u\in\mathbb{R}^m$, and $p(x)$ is an arbitrary polynomial of total degree at most $d$. This Lie group has dimension
\[
n+m\binom{n+d}{n}
\]
and acts transitively of $J^d(\mathcal{E})$. The upper bound (\ref{upperb3}) is
\begin{equation}\label{dmax}
n\cdot\left(n+m\binom{n+d}{n}-n-m\binom{n+d-1}{n}\right)=nm\binom{n+d-1}{n-1}.
\end{equation}
It is easily seen that a \emph{minimal} set of generating invariants for this Lie group action is $u^\alpha_J$ where $J\in\mathbb{Z}_{\geq0}^n$ with $|J|=d+1$. This minimal generating set has
\[
m\binom{n+d}{d+1}=m\binom{n+d}{n-1}
\]
elements. Hence, the upperbound (\ref{dmax}) is just the simplistic upper bound on the number of monomials in $n$ variables of degree $d+1$ obtained by taking the number of monomials of degree $d$ and multiplying by $n$.
\end{Ex}

\end{section}

\section{Lie pseudo-groups and moving frames}\label{sec:LG}

This section gives a rapid introduction to the formalism and constructions of the groupoid approach to Lie pseudo-groups, \cite{OP-2005}, and the equivariant moving frame, \cite{OP-2008}. The setup mirrors that of the first paper \cite{IMF}.

\subsection{Pseudo-groups}
Let $\mathcal{E}\to\mathcal{X}$ be a smooth fiber bundle with fibers $\mathcal{U}$ and local coordinates $z=(x,u)$, $x\in\mathbb{R}^n$, $u\in\mathbb{R}^m$. Our results are mostly local and it will essentially suffice to consider $\mathcal{E}=\mathcal{X}\times\mathcal{U}\to\mathcal{X}$ where $\mathcal{X}\subset\mathbb{R}^n$ and $\mathcal{U}\subset\mathbb{R}^m$ are open sets. Regarding the local coordinates, we shall call the $x$-coordinates \emph{horizontal source coordinates} and the $u$-coordinates \emph{vertical source coordinates}. The \emph{diffeomorphism} pseudo-group on $\mathcal{E}$, $\mathcal{D}(\mathcal{E})$ (or $\mathcal{D}$ when $\mathcal{E}$ is implied), is the collection of all (locally defined) invertible smooth maps, $\varphi$, between open sets in $\mathcal{E}$. We abuse notations slightly and write $\varphi:\mathcal{E}\to\mathcal{E}$ for these maps, even though they are only locally defined. This convention will hold throughout this paper and will apply to all locally defined objects.

The pseudo-group $q$-jets $j^q\varphi\restrict{z}$, $\varphi\in\mathcal{D}$, form a groupoid $\mathcal{D}_q$ fibered over $\mathcal{E}$ (for all $q\leq\infty$) with source and target maps (let $j^q\varphi\restrict{z}=(z, Z, \ldots, Z^a_A, \ldots)\in\mathcal{D}_q$)
\[
\sigma(j^q\varphi\restrict{z})=z,\quad \tau(j^q\varphi\restrict{z})=Z
\]
and groupoid operation $j^q\psi\restrict{Z}\cdot j^q\varphi\restrict{z}=j^q(\psi\raise 1pt \hbox{$\,\scriptstyle\circ\,$}\varphi)\restrict{z}$ where the composition is the standard composition operation on jets and the source of $j^q\psi\restrict{Z}$ is equal to the target of $j^q\varphi\restrict{z}$. Each $\psi\in\mathcal{D}$ acts on the groupoid elements $\mathcal{D}_q$ that have source within the domain of definition of $\psi$ by
\begin{equation}\label{eq:act}
R_\psi\cdot j^q\varphi\restrict{z}:=j^q\varphi\restrict{z}\cdot j^q\psi^{-1}\restrict{\psi(z)}.
\end{equation}
Note that $R_\psi$ agrees with the action $z\to\psi(z)$ on $\mathcal{E}$. For $\varphi\in\mathcal{D}$ we denote its jets by
\[
j^q\varphi\restrict{z}=(z, Z, \ldots, Z^a_A, \ldots)_{|A|\leq q},
\]
where $1\leq a\leq n+m$ and $A\in\mathbb{N}_0^{n+m}$ is a multi-index. We define the vector of jets $Z^{(q)}$ by $(z, Z^{(q)}):=(z, Z, \ldots, Z^a_A, \ldots)_{|A|\leq q}$.

\begin{Rem}
In the following, in lieu of explicit mention, the lower case Latin letters, $i,j,k,l$ will denote indices between $1$ and $n$, while lower case Greek letters, $\alpha, \beta$, etc. will denote indices between $1$ and $m$. The upper case Latin letters $J,K,L$ denote multi-indices in $\mathbb{N}_0^n$. The letters $a, b, c$ are reserved for superscripts in the jets $Z$ and run from $1$ to $m+n$ and their upper case counterparts $A,B,C$ denote multi-indices in $\mathbb{N}_0^{n+m}$.
\end{Rem}

The contact 1-forms on $\mathcal{D}_\infty$ (those 1-forms that are annihilated by the (pull-back of the) maps $z\mapsto j^\infty\varphi\restrict{z}$) have a basis that is invariant under the action (\ref{eq:act}). These forms are called the Maurer-Cartan forms of $\mathcal{D}$. They are denoted $\mu^a_A$ and the finite collections $\{\mu^a_A\}_{|A|<q}$ are a basis of the contact 1-forms on $\mathcal{D}_q$. Their structure equations are
\begin{equation}\label{eq:Dstr}
d\mu^a_A=\sum_{1\leq b\leq n+m}\o^b\wedge\mu^a_{A, b}+\sum_{\substack{B+C=A \\ |C|\geq1}}\sum_{1\leq b\leq n+m}\mu^a_{B, b}\wedge\mu^b_C,
\end{equation}
where $\o^b:=\sum_cZ^b_cdz^c$ are invariant forms on the base $\mathcal{E}$. See \cite{OP-2005}

Now let $\H$ be a sub Lie pseudo-group of $\mathcal{D}(\mathcal{E})$ with sub groupoids $\H_p\subset\mathcal{D}_p$. This means that in local coordinates the elements of $\H$ are local solutions to a system of PDEs (called the \emph{determining, or defining, equations}),
\begin{equation}\label{eq:def}
F^{(t)}(z, Z^{(t)})=0
\end{equation}
where each expression in $F^{(t)}$ is assumed to be genuinely $t\th$ order and full rank in the $t\th$ order jets. We assume this system is \emph{formally integrable} and real analytic and therefore locally solvable. Also, we assume that (\ref{eq:def}) is a complete system for $\H$, in the sense that
\[
\H_p=\left\{F^{(t)}=0~|~ t\leq p\right\}.
\]

\begin{Rem}\label{rem:order}
The \emph{order} of the Lie pseudo-group $\H$ is the smallest number $t^*$ such that all determining equations are differential consequences of the system $\left\{F^{(t)}=0~|~ t\leq t^*\right\}$. The important property of $t^*$ is that any local solution to $\left\{F^{(t)}=0~|~ t\leq t^*\right\}$ will define a pseudo-group element, as that solution will solve all the determining equations.
\end{Rem}

 Linearizing (\ref{eq:def}) at the \emph{identity section} $\one^t\restrict{z}$ (the $t$-jet of the identity solution $z\mapsto z$) gives the linear system
\[
L^{(t)}(z, Z^{(t)})=\sum_{\substack{Z^a_A\\ 0\leq|A|\leq t}}\left(\frac{\partial F^{(t)}}{\partial Z^a_A}\restrictbig{\one^t\restrict{z}}\right)Z^a_A.
\]
and restricting (or pulling back) the Maurer-Cartan forms $\mu^a_A$, to $\H_t$ will introduce linear dependencies among them precicely prescribed by these linear equations such that on $\H_p$ we have (see \cite{OP-2005})
\[
\sum_{\substack{Z^a_A\\ 0\leq|A|\leq t}}\left(\frac{\partial F^{(t)}}{\partial Z^a_A}\restrictbig{\one^t\restrict{z}}\right)_{z\leftrightarrow Z}\mu^a_A=0,
\]
where, in the coefficients of $\mu^a_A$, we have replaced source coordinates $z$ by target coordinates $Z$.

By prolongation, $\H$ acts on the (submanifold) jet-bundles $J^q(\mathcal{E})$, whose elements we denote by $j^q u\restrict{x}=(x, u,\ldots, u^\alpha_J,\ldots)$, $1\leq\alpha\leq m$ and $J\in\mathbb{N}_0^n$. Recall the total derivative on $J^\infty(\mathcal{E})$,
\[
D_i=\pder{x_i}+\sum_{\beta, K}u^\beta_{K,i}\pder{u^\beta_K},
\]
s.t. $u^\alpha_J=D^Ju^\alpha$. We denote the target coordinates of the prolonged action with capitalized, hat-wearing, letters:
\begin{equation}\label{eq:hatw}
j^\infty\varphi\restrict{z}\cdot j^\infty u\restrict{x}=(X, U, \ldots, \widehat{U}^\alpha_J, \ldots),
\end{equation}
where $j^\infty\varphi\restrict{z}$ and  $j^\infty u\restrict{x}$ have the same source coordinate in $\mathcal{E}$, under $\H_\infty\to\mathcal{E}$ and $J^\infty(\mathcal{E})\to\mathcal{E}$, respectively.

The pull-back of $\H_q\to\mathcal{E}$ by the canonical projection $J^q(\mathcal{E})\to\mathcal{E}$, along with the prolonged action, gives the doubly fibered space
\[
\begin{tikzpicture}
  \node (A0) at (0,0) {$\widetilde{\H}_q$};
  \node (A1) at (240:2cm) {$J^q(\mathcal{E})$};
  \node (A2) at (300:2cm) {$J^q(\mathcal{E})$,};
  \draw[->,font=\scriptsize]
  (A0) edge node[left] {$\tilde{\sigma}$} (A1)
  (A0) edge node[right] {$\tilde{\tau}$} (A2);

\end{tikzpicture}
\]
where, in local coordinates,
\[
\tilde{\sigma}(j^q\varphi\restrict{z}, j^q u\restrict{x})=j^qu \restrict{x},\quad \tilde{\tau}(j^q\varphi\restrict{z}, j^q u\restrict{x})=j^q\varphi\restrict{z}\cdot j^q u\restrict{x}.
\]
The pseudo-group $\H$ acts on the bundles $\widetilde{\H}_q$ by \emph{right-regularization}, by
\begin{equation}\label{eq:rr}
\psi\cdot (j^q\varphi\restrict{z}, j^q u\restrict{x}):=(j^q\varphi\restrict{z}\cdot j^q\psi^{-1}\restrict{\varphi(z)}, j^q\psi\restrict{z}\cdot j^q u\restrict{x})
\end{equation}
where the domain of definition of $\psi$ includes $z$. Importantly, the target map $\tilde{\tau}$ is invariant under this action and so the pull-back of any differential form of $J^q(\mathcal{E})$ under $\tilde{\tau}$ is invariant on $\widetilde{\H}_q$. In particular, the \emph{lifted invariants} are the pull-backs of the coordinate functions on $J^\infty(\mathcal{E})$, $\tilde{\tau}(u^\alpha_J)=\widehat{U}^\alpha_J$. The exterior derivative of lifted invariants is given by the fundamental \emph{recurrence formula},
\begin{equation}\label{eq:rec}
d\widehat{U}^\alpha_J=\sum_i^n\widehat{U}^\alpha_{J, i}\omega^i+\tilde{\tau}^*\left(\pder[\widehat{U}\av_J]{Z^a_A}\restrictbig{\one}\right)\mu^a_A+\text{cont.}
\end{equation}
where ``cont.'' are contact forms on the bundle $J^{\infty}(\mathcal{E})$ and automatically vanish when we restrict to a section of $\mathcal{E}$. Notice that, in the above formula, the first sum is only taken over the first $n$ base forms $\o^i$, $1\leq i\leq n$.

It often happens that, when constructing a moving frame (see below), we find a lifted invariant $\widehat{U}^\alpha_J$ that depends on the (submanifold) jets $J^q(\mathcal{E})$ but on the pseudo-group parameters in $\H_p$, for $q>p$. To accommodate this feature, we define the following spaces. Let $q>p$ and consider the pull-back of $\H_p\to\mathcal{E}$ by $J^q(\mathcal{E})\to\mathcal{E}$. Call this space $\widetilde{\H}_{p,q}$. The canonical projection $\H_q\to\H_p$ induces a projection $\widetilde{\H}_q\to\widetilde{\H}_{p,q}$. The local coordinates of $\widetilde{\H}_{p,q}$ are
\[
(j^p\varphi\restrict{z}, j^q u\restrict{x})
\]
and the the pseudo-group action $\H$ on $\widetilde{\H}_q$ reduces to an action of $\H$ on $\widetilde{\H}_{p,q}$,
\begin{equation}\label{eq:rrpq}
\psi\cdot (j^p\varphi\restrict{z}, j^q u\restrict{x}):=(j^p\varphi\restrict{z}\cdot j^p\psi^{-1}\restrict{\varphi(z)}, j^q\psi\restrict{z}\cdot j^q u\restrict{x}).
\end{equation}

\subsection{Equivariant moving frames}

A \emph{partial moving frame} $\widetilde{\mathcal{B}}_{p,q}$ is any smooth $\H$-invariant subspace of $\widetilde{\H}_{p,q}$ under the right-regularized action (\ref{eq:rrpq}). We refer to the $p$ as the \emph{group order} of the partial moving frame and $q$ as its \emph{jet order}. The domain of definition of $\widetilde{\mathcal{B}}_{p,q}$ is the projection
\[
\S_q=\tilde{\sigma}\left(\widetilde{\mathcal{B}}_{p,q}\right)\subset J^q(\mathcal{E})
\]
where $\tilde{\sigma}$ is the source map in $\widetilde{\H}_q$. Note that $\S_q$ is an invariant set.

In practice, a partial moving frame is constructed as the level set of some of the lifted invariants $\widehat{U}^\alpha_J$, viewed as functions on some $\widetilde{\H}_{p,q}\to\mathbb{R}$. The recurrence formula (\ref{eq:rec}) then tells us what linear dependencies among the Maurer-Cartan forms are introduced by restricting them to $\widetilde{\mathcal{B}}_{p,q}$.

Let $\widetilde{\mathcal{B}}_{p,q}\to\S_q$ be a partial moving frame with domain of definition $\S_q$. Ultimately we are interested in the local congruence of (local) sections of $\mathcal{E}\to\mathcal{X}$ under $\H$. For a section $S$ of $\mathcal{E}$ denote the jets of $S$ by $j^qS$ and assume $im(j^qS)\subset\S_q$. Pulling the bundle $\widetilde{\mathcal{B}}_{p,q}\xrightarrow{\tilde{\sigma}}\S_q$ back by $j^qS$ will give a bundle\footnote{We continue to use the convention that functions and sections need only be defined on subsets of their domain space, and we extend this convention to bundles.} $\mathcal{B}^{j^qS}_p\to \mathcal{X}$. This is a subbundle of the pull-back bundle $S^*(\H_p\to\mathcal{E})$ over $\mathcal{X}$.

\begin{Rem}\label{rem:emb}
Note that since $S:\mathcal{X}\to\mathcal{E}$ is one-to-one, $S^*(\H_p\to\mathcal{E})$ and $\mathcal{B}^{j^qS}_p\to\mathcal{X}$ are subbundles of $\H_p\to\mathcal{E}$ under an embedding we denote $\bar{S}$.
\end{Rem}

The solution to the congruence problem for \emph{horizontal actions} in \cite{IMF} rested on results on the equivalence of sections of the bundle $S^*(\H_p\to\mathcal{E})$. However, the horizontal case was made simple by the fact that $S^*(\H_p\to\mathcal{E})$ was isomorphic to the groupoid $\mathcal{G}_p\to\mathcal{X}$ of a certain isomorphic Lie pseudo-group $\mathcal{G}\simeq\H$. (In the language of \cite{S-2000} $\H$ was a one-to-one prolongation of $\mathcal{G}$.) In the following we shall write $\H^S_p\to\mathcal{X}$ instead of $S^*(\H_p\to\mathcal{E})$.

Unfortunately, in the non-horizontal case the bundle $\H^S_p\to\mathcal{X}$ is no longer isomorphic to the groupoid $\mathcal{G}_p\to\mathcal{X}$ of some known Lie pseudo-group as in \cite{IMF}. However, the proof of the key theorem \cite[Theorem 3.12]{IMF} carries over, word-for-word, to prove the following modification. A quick comment on notation: As in the \cite{IMF}, let $\H_p\irestrict{Z_0}\to\mathcal{E}$ denote the subbundle of $\H_p\to\mathcal{E}$ given by $\tau^{-1}(Z_0)$. This is the bundle we obtain when we normalize all zero order lifted invariants to $Z_0$ in the process of building a moving frame. Assuming we can normalize all zero order lifted invariants is not necessary but does simplify the discussion (see \cite{IMF}).

\begin{theorem}\label{GcoX}
Let $\H$ be a Lie pseudo-group, of order $t^*$, of transformations on the base manifold $\mathcal{E}$, and let $S_1$ and $S_2$ be sections of $\mathcal{E}$. Let $s_j$ be a section of $\H^{S_j}_p\irestrict{Z_0}\to \mathcal{X}$ for $j=1,2$, and $p\geq t^*$. A local transformation $f:\mathcal{X}\to \mathcal{X}$ preserves all the pulled-back Maurer-Cartan forms (recall Remark \ref{rem:emb} for the definition of $\bar{S}$),
\begin{equation}\label{eq:GcoX}
f^*s_2^*\bar{S}_2^*\mu^a_A=s_1^*\bar{S}_1^*\mu^a_A,\quad |A|<p,
\end{equation}
if and only if there exists a section $f_p$ of $\H^{S_1}_p\to \mathcal{X}$ such that
\[
f_p^*\bar{S}_1^*\mu^a_A=0,\quad |A|<p,
\]
and $R_{\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p}\cdot(\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} s_1)=\bar{S}_2\raise 1pt \hbox{$\,\scriptstyle\circ\,$} s_2$.
\end{theorem}

\begin{Rem}
Notice that the above theorem is only true when $p$ is at least as great as the order of the Lie pseudo-group $\H$. This is because the proof relies on local solvability in the groupoid $\H_p$, i.e. for each point solution to
\[
F^{(t)}=0, \quad t\leq p
\]
there must be a pseudo-group element passing through that point. See \cite{IMF} for details.
\end{Rem}

Notice that $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p$ annihilates the contact codistribution of $\H_p$ of which the Maurer-Cartan forms $\mu^a_A$ are a basis and so there is hope that $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p$ can be \emph{integrated} to a bona-fide pseudo-group element of $\H$. That is, $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p$ could be given as an initial condition for the system of determining equations of $\H$. We make a general definition of the pseudo-groups for which this is possible, but hold off on characterizing convenient subclasses of these further until Section \ref{sec:qhpg} below.

\begin{definition}
Let $\H$ be a Lie pseudo-group, of order $t^*$, of transformations on the bundle $\mathcal{E}\to\mathcal{X}$. For $r\geq t^*$, the pseudo-group $\H$ is called \emph{$r\th$ order vertically integrable} if, for each section $S$ of $\mathcal{E}$, and each $p\geq r$, every map $f_p:\mathcal{X}\to\H^S_p$ such that $\bar{S}\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p$ annihilates all Maurer-Cartan forms $\mu^a_A$, $|A|<p$ on $\H_p\to\mathcal{E}$ and each point $f_p(x)$ on the image of $f_p$, there exists a neighborhood $V$ of $S(x)\in\mathcal{E}$  and a pseudo-group element $\varphi\in\H$ defined on $V$ such that $j^p\varphi\raise 1pt \hbox{$\,\scriptstyle\circ\,$} S=\bar{S}\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p$ on $S^{-1}(V)$, and $r$ is the smallest such number.
\end{definition}

For vertically-integrable pseudo-groups we have the following consequence of the previous theorem.

\begin{theorem}\label{GcoX2}
Let $\H$ be a $r\th$ order vertically integrable Lie pseudo-group of transformations on $\mathcal{E}\to\mathcal{X}$, and let $S_1$ and $S_2$ be sections of $\mathcal{E}$. Let $s_j$ be a section of $\H^{S_j}_p\irestrict{Z_0}\to \mathcal{X}$ for $j=1,2$ and $p\geq r$. A local transformation $f:\mathcal{X}\to \mathcal{X}$ preserves all the pulled-back Maurer-Cartan forms,
\begin{equation}\label{eq:GcoX}
f^*s_2^*\bar{S}_2^*\mu^a_A=s_1^*\bar{S}_1\mu^a_A,\quad |A|<p,
\end{equation}
if and only if there is a $\varphi\in\H$ such that $\varphi\raise 1pt \hbox{$\,\scriptstyle\circ\,$} S_1=S_2$.
\end{theorem}

Due to Theorem \ref{GcoX2} the involutive moving frame \cite{IMF} completely carries over to these pseudo-groups, including the following fundamental theorem on congruence in $\mathcal{E}$ under $\H$, which has, verbatim, the same proof as in the horizontal case, \cite[Theorem 4.15]{IMF}.

\begin{theorem}\label{thm:total}
Let $S_1$ and $S_2$ be sections of $\mathcal{E}$ and let $\H$ be an $r\th$ order vertically integrable Lie pseudo-group of transformations on $\mathcal{E}$. Let $\mathcal{B}_{p,q}\to\S_q$ be a partial moving frame with $p\geq r$, on which $U=\tilde{\tau}^*(u)=\text{constant}$. Then the following are equivalent.
\begin{enumerate}
\item $S_1$ and $S_2$ are locally congruent under some $\varphi\in\H$.
\item For each section $s_1$ of $\mathcal{B}^{j^qS_1}_p\to \mathcal{X}$ there is a section $s_2$ of $\mathcal{B}^{j^qS_2}_p\to \mathcal{X}$  and a map $f$ such that
\[
f^*s_2^*\bar{S}_2^*\mu^a_A=s_1^*\bar{S}_1\mu^a_A,\quad |A|<p.
\]
\item There exists a pair of sections, $s_j$ of $\mathcal{B}^{j^qS_j}_{p}\to \mathcal{X}$, $j=1,2$, and a map $f$ such that
\[
f^*s_2^*\bar{S}_2^*\mu^a_A=s_1^*\bar{S}_1\mu^a_A,\quad |A|<p.
\]
\end{enumerate}
\end{theorem}

This theorem presents us with two problems. First we must solve the coframe equivalence problem in part 3 of Theorem \ref{thm:total} to obtain a function $f:\mathcal{X}\to\mathcal{X}$ and its corresponding $f_p$ of Theorem \ref{GcoX}. Then we must integrate $f$ and $f_p$ to a pseudo-group element $\varphi\in\H$. In this paper we shall focus on those vertically integrable pseudo-groups where the section $f_p$ can be \emph{uniquely} integrated to a pseudo-group element $\varphi\in\H$. We shall call these pseudo-groups \emph{quasi-horizontal} and we now turn to characterizing them in terms of their determining equations. Note that for these quasi-horizontal pseudo-groups it is enough to find the $f$ of part 3 of Theorem \ref{thm:total} in order to completely solve the equivalence problem. Hence, the involutive moving frame method from \cite{IMF} will completely carry over to the quasi-horizontal case.

\subsection{Quasi-horizontal Lie pseudo groups}\label{sec:qhpg}

We now turn to the characterization of Lie pseudo-groups for which the sections $f_p$ of Theorem \ref{GcoX} can be uniquely integrated to pseudo-group elements.

\begin{definition}\label{def:quasi}
Let $\H$ be a Lie pseudo-group, of order $t^*$, of transformations on $\mathcal{E}$ with determining equations $F^{(t)}(t, Z^{(t)})=0$. Given local coordinates $z=(x,u)$ on $\mathcal{E}$, the elements of $\H$ are local diffeomorphisms $(x,u)=z\mapsto Z=(X,U)$, on $\mathcal{E}$. Then $\H$ is called \emph{quasi-horizontal} if there exists an $r\geq t^*$ such that the following holds.

In the determining equations of $\H$, each derivative of order $r$, $Z^a_A$, $|A|=r$, which has at least one $u^\alpha$ derivative can be taken as principal. In other words, for each $A\in\mathbb{Z}_{\geq0}^{r-1}$, $1\leq a\leq n+m$ and $1\leq\alpha\leq m$, one of the determining equations $F^{(r)}=0$ for $\H$ reads
\[
\frac{\partial^{r}Z^a}{\partialu^\alpha\partial z^A}=F^{a}_{\alpha, A}(z, Z^{(r)}).
\]

Notice that if the two conditions hold for some $r$, they also hold for all $r'\geq r$. We say $\H$ is \emph{quasi-horizontal of order} $r$ if $r$ is the smallest number for which the above holds.
\end{definition}

\begin{Rem}\label{rem:quasi}
Importantly, if $\H$ is quasi-horizontal in variables $(x, u)$ it is also quasi-horizontal (with the same horizontal order) after any fiber preserving change of variables in $\mathcal{E}$, $(x, u)\mapsto (y(x), w(x, u))$. The above definition is therefore consistent on a general bundle $\mathcal{E}\to\mathcal{X}$.
\end{Rem}

\begin{Rem}
Horizontal pseudo-groups, as defined in \cite{IMF}, are trivially quasi-horizontal.
\end{Rem}

Now let $f:\mathcal{X}\to\mathcal{X}$ and $f_p:\mathcal{X}\to\H^{S_1}_p$ satisfy the conditions of Theorems \ref{GcoX} and \ref{thm:total} where $\H$ is quasi horizontal (of order $r\leq p$) and assume the images of $S_j$ are the graphs of two functions $u_j:\mathcal{X}\to\mathcal{U}$, $j=1,2$. We can straighten $S_1$ out in the $x$-direction by a (fiber-preserving) change of variables,
\[
(y,w)\mapsto (y, u_1(y)+w).
\]
We call the coordinates $(y, w)$ \emph{flat} for $S_1$, and by Remark \ref{rem:quasi} $\H$ is also quasi-horizontal in $(y, w)$. Notice that on $S_1$ we have
\begin{equation}\label{eq:dy}
\pder{y}=\pder{x}+u_x\pder{u}.
\end{equation}
In the flat coordinates $S_1$ is the zero section $y\to (y, 0)$ and we shall treat $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p(y, 0)$ as initial conditions for a sequence of normal systems of partial differential equations for the vector
\[
\pmb{Z}:=Z^{(p-1)}=(Z^1, \ldots, Z^a_{A},\ldots)\restrict{|A|\leq p-1},
\]
where $A=(A_1, \ldots, A_{n+m})\in\mathbb{N}_0^{n+m}$ now denote derivatives w.r.t. $(y, w)$. The $p\th$ order determining equations for $\H$ involve some equations that only contain $y$-derivatives of $\pmb{Z}$, but these equations are satisfied by $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p(y, 0)$ by construction. The rest of the $p\th$ order equations involve at least one $w$-derivative. Since $\H$ is quasi-horizontal, using the standard trick of turning a PDE of arbitrary order into a first order system, we can write these latter equations as
\begin{equation}\label{eq:Zb}
\pder{w^1}\pmb{Z}=\pmb{F}_1(y, w, \pmb{Z},\pder{y_1}\pmb{Z},\ldots, \pder{y_n}\pmb{Z})
\end{equation}
for some function $\pmb{F}_1$.

Now, given the map $\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p(y, 0)$ consider the initial value problem (\ref{eq:Zb}) with initial value $\pmb{Z}(y, 0)=\bar{S}_1\raise 1pt \hbox{$\,\scriptstyle\circ\,$} f_p(y,0)$. This is obviously a normal system and uniquely integrable to an open set in $\mathbb{R}^{n+1}\times\{0\}^{m-1}$ by the Cauchy-Kovalevskaya theorem. We can then use the newly integrated $\pmb{Z}(y, w^1, 0)$ as an initial condition for a system
\[
\pder{w^2}\pmb{Z}=\pmb{F}_2(y, w, \pmb{Z},\pder{y_1}\pmb{Z}, \pder{y_n}\pmb{Z})
\]
and so on until we have integrated $\pmb{Z}$ to the prolongation of a pseudo-group element $\varphi\in\H$. This proves the following.

\begin{theorem}\label{thm:quastotal}
Let $\H$ be quasi horizontal of order $r$ and let $f$ and $f_p$ satisfy the conditions of Theorems \ref{GcoX} and \ref{thm:total} with $p\geq r$. Then $f$ can be uniquely integrated to a pseudo-group element $\varphi\in\H$.
\end{theorem}

Theorem \ref{thm:quastotal} allows us to transfer the equivalence method for horizontal actions from \cite{IMF} to the quasi horizontal case as outlined in the comments following Theorem \ref{thm:total}. We now demonstrate the use of Theorem \ref{thm:total} by studying the equivalence problem of first order partial differential equations.