This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

On Horizontal Immersions of Discs in Fat Distributions of Type (4,6)(4,6)

Aritra Bhowmick Statistics and Mathematics Unit, Indian Statistical Institute
203, B.T. Road, Kolkata 700108, India
[email protected]
Abstract.

In this article we discuss horizontal immersions of discs in certain corank-22 fat distributions on 66-dimensional manifolds. The underlying real distribution of a holomorphic contact distribution on a complex 33 manifold belongs to this class. The main result presented here says that the associated nonlinear PDE is locally invertible. Using this we prove the existence of germs of embedded horizontal discs.

Key words and phrases:
h-principle, fat distribution, holomorphic contact distribution, elliptic PDE
2010 Mathematics Subject Classification:
58A30, 58J05, 58A15, 35J60, 53C23

1. Introduction

In subriemannian geometry one studies distributions on smooth manifolds. Bracket-generating distributions, which lie at the opposite end of the integrable ones, are the primary focus. A distribution 𝒟TM\mathcal{D}\subset TM is bracket generating if successive Lie brackets of local vector fields in 𝒟\mathcal{D} around any point xMx\in M span the tangent space TxMT_{x}M.

1.1. Immersions In a Manifold with a Distribution

Given a distribution 𝒟\mathcal{D}, we can consider smooth curves γ:𝕀M\gamma:\mathbb{I}\to M which are everywhere tangent to it. We shall call them 𝒟\mathcal{D}-horizontal curves or simply horizontal curves. These curves play an important role in understanding the distribution. In fact, if 𝒟\mathcal{D} is bracket generating, any two points of the manifold MM can be joined by a smooth horizontal curve ([Cho39]) and the space of horizontal curves joining two points has the same homotopy type as the space of smooth curves joining them ([Ge93, Gro96]).

More generally, one may consider horizontal immersions of kk-manifolds in a manifold MM endowed with a distribution 𝒟\mathcal{D}. Such maps can be thought of as the solutions to the differential operator,

𝔇:C(Σ,M)\displaystyle\mathfrak{D}:C^{\infty}(\Sigma,M) Ω1(Σ,p)=Γhom(TΣ,p)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{p})=\Gamma\hom(T\Sigma,\mathbb{R}^{p})
f\displaystyle f (fλ1,,fλp)\displaystyle\mapsto\big{(}f^{*}\lambda^{1},\ldots,f^{*}\lambda^{p}\big{)}

where we assume that 𝒟\mathcal{D} is a corank-pp distribution on MM, given as the common kernel, 𝒟=i=1pkerλi\mathcal{D}=\cap_{i=1}^{p}\ker\lambda^{i}, for λiΩ1(M),i=1,,p\lambda^{i}\in\Omega^{1}(M),i=1,\ldots,p. Gromov defines ([Gro86, pg. 338]) a 𝒟\mathcal{D}-horizontal immersion f:ΣMf:\Sigma\to M to be regular if the algebraic system

dλi(,uXj)=σij,1ip, 1jdimΣd\lambda^{i}\big{(}\partial,\;u_{*}X_{j}\big{)}=\sigma_{ij},\quad 1\leq i\leq p,\;1\leq j\leq\dim\Sigma

is solvable for (local) vector fields Γf𝒟\partial\in\Gamma f^{*}\mathcal{D}, for any given set of arbitrary smooth functions {σij}\{\sigma_{ij}\}. Here, {Xj}\{X_{j}\} is some fixed local framing of TΣT\Sigma. It follows that the linearization of 𝔇\mathfrak{D} at some f:ΣMf:\Sigma\to M is invertible provided ff is regular. In fact, the linearization operator at a regular map ff has a 0th0^{\text{th}}-order inversion. Then, using a version of the Nash implicit function theorem, Gromov ([Gro86]) proceeds to obtain the hh-principle for regular 𝒟\mathcal{D}-horizontal immersions. This problem has been revisited with further details in [Gro96, Pan16].

Now, for any horizontal immersion ff, Imdfx\operatorname{Im}df_{x} is an isotropic subspace of 𝒟f(x)\mathcal{D}_{f(x)}. So, accounting for this isotropy, the above algebraic system is underdetermined whenever

rk𝒟dimΣ>dimΣ×corank𝒟.\operatorname{rk}\mathcal{D}-\dim\Sigma>\dim\Sigma\times\operatorname{corank}\mathcal{D}.

Gromov proves ([Gro96, pg. 256]) that a generic distribution 𝒟\mathcal{D} on MM admits horizontal germs of immersions through generic points of MM, under the above inequality. Then he proceeds to state the hh-principle ([Gro96, pg. 258]) for 𝒟\mathcal{D}-horizontal immersions of ΣM\Sigma\to M, with some additional regularity condition, whenever

rk𝒟(dimΣ+1)(corank𝒟+1)\operatorname{rk}\mathcal{D}\geq(\dim\Sigma+1)(\operatorname{corank}\mathcal{D}+1)

holds. He also conjectures ([Gro96, pg. 259]) that one may be able to improve the inequality to dimM(dimΣ+1)(corank𝒟+1)\dim M\geq(\dim\Sigma+1)(\operatorname{corank}\mathcal{D}+1).

1.2. Horizontal Immersions in a Manifold with a Fat Distribution

Fat distributions (see subsection 2.2) form an interesting class of bracket generating distributions; a distribution 𝒟\mathcal{D} is called fat (or strongly bracket generating) if any non-vanishing local section of 𝒟\mathcal{D} Lie bracket generates the tangent space in 11-step. Contact distributions are the prime examples of fat distributions on an odd dimensional manifold MM; they are locally defined by 1-forms α\alpha such that dαd\alpha is non-degenerate on kerα\ker\alpha. Horizontal immersions for such distributions are well understood; in fact, they satisfy the parametric hh-principle ([Duc84, EM02]).

In corank 11, fatness is a generic property. However, for higher corank, fat distribution germs are never generic (see 2.7). In corank 22, the most prominent examples of fat distributions are given by the holomorphic analogue of contact structures (see 2.9). These manifolds are modeled on the holomorphic 1-jet space J1(n,)J^{1}(\mathbb{C}^{n},\mathbb{C}) like their real counterparts, and 1-jet prolongation of any holomorphic map n\mathbb{C}^{n}\to\mathbb{C} is a holomorphic Legendrian embedding. So, there are plenty of holomorphic horizontal submanifolds in any holomorphic contact manifold. In [FL18] the authors have shown that holomorphic Legendrian embeddings of an open Riemann surface Σ\Sigma into the standard holomorphic contact manifold (2n+1,dziyidxi)\big{(}\mathbb{C}^{2n+1},dz-\sum_{i}y_{i}dx_{i}\big{)} satisfy the parametric Oka principle. In particular, they prove that the space of Legendrian holomorphic embeddings Σ2n+1\Sigma\hookrightarrow\mathbb{C}^{2n+1} has the same homotopy type as the space of continuous maps Σ𝕊4n1\Sigma\to\mathbb{S}^{4n-1}. The authors further observe that such a global hh-principle type result may not be true for a general holomorphic contact manifold.

Now, recall that for a given (real) contact structure kerα\ker\alpha, we have the Reeb vector field RR, which is defined by

ιRdα|kerα=0,α(R)=1.\iota_{R}d\alpha|_{\ker\alpha}=0,\quad\alpha(R)=1.

Similarly, given a holomorphic contact structure Ξ\Xi on a complex manifold MM, the underlying real distribution 𝒟\mathcal{D} can be locally written as 𝒟=kerλ1kerλ2\mathcal{D}=\ker\lambda^{1}\cap\ker\lambda^{2} and furthermore, we can identify (local) vector fields Z1,Z2Z_{1},Z_{2}, such that TM=𝒟SpanZ1,Z2TM=\mathcal{D}\oplus\textrm{Span}\langle Z_{1},Z_{2}\rangle and Z1,Z2Z_{1},Z_{2} display Reeb-like properties. This motivates the definition of local Reeb like directions (Defn 2.10). Holomorphic contact distributions are the best known examples in the class of fat corank 22 distributions considered in this article, which admit local Reeb directions.

Gromov’s general theorem for regular horizontal immersions ΣM\Sigma\to M is not applicable to such a distribution 𝒟\mathcal{D} on a 66-dimensional manifold MM, whenever dimΣ=2\dim\Sigma=2, as there cannot exist any ‘regular’ map due to dimension constraints. And yet the results of [FL18] indicates that there is still a possibility of obtaining an hh-principle.

1.3. Main Results

In this article, we consider 𝒟\mathcal{D} to be a corank 22 fat distribution on a manifold of dimension 66, which admits (local) Reeb directions, and then we study horizontal immersions of the closed unit disc Σ=𝔻2\Sigma=\mathbb{D}^{2} into MM. The main theorem of this article may be stated as follows.

Theorem 1.1.

There exists an open subset 𝔘C(𝔻2,M)\mathfrak{U}\subset C^{\infty}(\mathbb{D}^{2},M) such that for every f:𝔻2Mf:\mathbb{D}^{2}\to M in 𝔘\mathfrak{U}, the linearization 𝔏f\mathfrak{L}_{f} of the operator 𝔇\mathfrak{D} at ff admits a tame right inverse.

The open subset 𝔘\mathfrak{U} consists of solutions to certain open first order differential relation in J1(Σ,M)J^{1}(\Sigma,M). We observe that the linearization 𝔏f\mathfrak{L}_{f} of 𝔇\mathfrak{D} at f𝔘f\in\mathfrak{U}, factors as a composition of two linear operators, one of which is elliptic. In particular, the inverse of 𝔏f\mathfrak{L}_{f} in Theorem 1.1 is not a differential operator unlike in the case when ff is a ‘regular’ map in the sense of Gromov.

Applying the Nash-Hamilton Implicit Function Theorem, due to R. Hamilton ([Ham82]), for smooth tame operators between tame Fréchet spaces we obtain the following result.

Theorem 1.2.

𝔇\mathfrak{D} is locally invertible on 𝔘\mathfrak{U}.

Next, we obtain jets of infinitesimal solutions to 𝔇\mathfrak{D}. An application of the local invertibility then gives us the following local hh-principle.

Theorem 1.3.

Horizontal maps 𝔻2M\mathbb{D}^{2}\to M satisfy the local h-principle.

We also prove the existence of horizontal germs.

Theorem 1.4.

There exist germs of 𝒟\mathcal{D}-horizontal submanifolds of dimension 22.

The article is arranged as follows. In section 2 we first discuss the preliminaries of fat distributions. In section 3 we prove our main result : the local invertibility of the nonlinear differential operator 𝔇\mathfrak{D}. Then in section 4 we derive the local h-principle for horizontal maps and prove the existence of germs of such maps. In order to make the article self-contained, we briefly outline in section 5 the background and the statement of the Implicit Function Theorem for differential operators between Frechet spaces following [Ham82], as this has been crucially used in section 3.

2. Basic Notions

In this section we recall the notion of fat distribution and then focus on distributions which are of corank 22.

2.1. Bracket-Generating Distributions

A distribution 𝒟\mathcal{D} on a manifold MM is a subbundle of the tangent bundle TMTM. A distribution 𝒟\mathcal{D} can equivalently be identified with its sheaf of sections Γ𝒟\Gamma\mathcal{D}. Hence by X𝒟X\in\mathcal{D} we will mean that XX is a local section of 𝒟\mathcal{D}. We denote by [𝒟,𝒟][\mathcal{D},\mathcal{D}] the sheaf of all vector fields on MM which are obtained by taking Lie brackets of two sections in 𝒟\mathcal{D}. The sheaf [𝒟,𝒟][\mathcal{D},\mathcal{D}] need not be associated to a distribution, since it may fail to have constant rank. We define recursively, for i1i\geq 1,

𝒟i+1=𝒟i+[𝒟,𝒟i],𝒟1=𝒟.\mathcal{D}^{i+1}=\mathcal{D}^{i}+[\mathcal{D},\mathcal{D}^{i}],\quad\mathcal{D}^{1}=\mathcal{D}.
Definition 2.1.

𝒟\mathcal{D} is said to be bracket generating if, at each point xMx\in M, we have, 𝒟r|x=TxM\mathcal{D}^{r}|_{x}=T_{x}M, for some positive integer rr, possibly depending on xx.

Example 2.2.

A corank 1 distribution ξTM\xi\subset TM on MM is defined locally as the kernel of a 1-form α\alpha on MM. The distribution ξ\xi is said to be a contact distribution if α(dα)n0\alpha\wedge(d\alpha)^{n}\neq 0. The property does not depend on the choice of defining 1-form α\alpha. By the Darboux theorem ([Gei08]), α\alpha has a local normal form α=dzi=1nyidxi\alpha=dz-\sum_{i=1}^{n}y_{i}dx_{i}. It is then easy to check that,

ξ=loc.Spanyi,xi+yiz.\xi\underset{loc.}{=}\textrm{Span}\langle\partial_{y_{i}},\partial_{x_{i}}+y_{i}\partial_{z}\rangle.

Since [yi,xi+yiz]=z[\partial_{y_{i}},\partial_{x_{i}}+y_{i}\partial_{z}]=\partial_{z}, we see that TM=ξ+[ξ,ξ]=ξ2TM=\xi+[\xi,\xi]=\xi^{2}. Hence, contact distributions are bracket generating.

In this article we will be interested in distributions which are 11-step bracket generating, i.e, TM=𝒟2TM=\mathcal{D}^{2}. If rk𝒟=n\operatorname{rk}\mathcal{D}=n and corank𝒟=p\operatorname{corank}\mathcal{D}=p then we say that 𝒟\mathcal{D} is of type (n,n+p)(n,n+p). As observed, every contact distribution is of type (2n,2n+1)(2n,2n+1) for some n1n\geq 1.

2.2. Fat Distributions

We can associate to any distribution 𝒟\mathcal{D}, its curvature form Ω\Omega which is a TM/𝒟TM/\mathcal{D}-valued 2-form on 𝒟\mathcal{D}. The curvature form plays a very important role in subriemannian geometry.

Definition 2.3.

Given a distribution 𝒟\mathcal{D}, we have the quotient map λ:TMTM/𝒟\lambda:TM\to TM/\mathcal{D}. The curvature form of the distribution is a map Ω:Λ2𝒟TM/𝒟\Omega:\Lambda^{2}\mathcal{D}\to TM/\mathcal{D} defined as follows:

Ω(X,Y):=[X,Y]mod𝒟=λ([X,Y]),for local sections X,Y𝒟.\Omega(X,Y):=-[X,Y]\mod\mathcal{D}=-\lambda([X,Y]),\quad\text{for local sections $X,Y\in\mathcal{D}$.}

Observe that Ω\Omega is C(M)C^{\infty}(M)-linear. 𝒟\mathcal{D} is 11-step bracket-generating precisely when Ω\Omega is surjective. Now, there is a special class of bracket generating distributions, called fat. These can be defined in several equivalent ways as described below.

Strong Bracket Generation :

A distribution 𝒟\mathcal{D} on a manifold MM is called fat (or strongly bracket generating) at xMx\in M if, for every nonzero vector v𝒟xv\in\mathcal{D}_{x}, we have that

TxM=𝒟x+[V,𝒟]x,T_{x}M=\mathcal{D}_{x}+[V,\mathcal{D}]_{x},

where VV is some (local) section of 𝒟\mathcal{D} with Vx=vV_{x}=v and [V,𝒟]x[V,\mathcal{D}]_{x} is the subspace of TxMT_{x}M defined as follows:

[V,𝒟]x={[V,X]x|X𝒟}[V,\mathcal{D}]_{x}=\big{\{}[V,X]_{x}\big{|}X\in\mathcal{D}\big{\}}

The distribution is fat if it is fat at every point xMx\in M.

Nondegeneracy of Curvature Form :

Suppose we identify (TM/𝒟)\big{(}TM/\mathcal{D}\big{)}^{*} with the annihilator bundle 𝒟TM\mathcal{D}^{\perp}\subset T^{*}M and define the dual curvature form ω:𝒟Λ2𝒟\omega:\mathcal{D}^{\perp}\to\Lambda^{2}\mathcal{D}^{*} by ω(α)=dα|𝒟\omega(\alpha)=d\alpha|_{\mathcal{D}}. We say that 𝒟\mathcal{D} is fat at xMx\in M if and only if ω(α)\omega(\alpha) is a nondegenerate 22-form on 𝒟\mathcal{D}, for each 0α𝒟x0\neq\alpha\in\mathcal{D}^{\perp}_{x}. The distribution 𝒟\mathcal{D} is fat if it is fat at every xMx\in M

Remark 2.4.

Clearly, any fat distribution is 11-step bracket generating. Another important consequence of fatness is that dα|𝒟d\alpha|_{\mathcal{D}} is nondegenerate for any non-vanishing, (local) 1-form α\alpha annihilating 𝒟\mathcal{D}. Indeed, for a corank pp distribution, the strongly bracket generating property implies (and is implied by) that for any 0v𝒟x0\neq v\in\mathcal{D}_{x}, the map

Φv:𝒟x\displaystyle\Phi_{v}:\mathcal{D}_{x} TxM/𝒟xp\displaystyle\to T_{x}M/{\mathcal{D}_{x}}\cong\mathbb{R}^{p}
u\displaystyle u (ιvω1(u),,ιvωp(u))\displaystyle\mapsto\big{(}\iota_{v}\omega_{1}(u),\ldots,\iota_{v}\omega_{p}(u)\big{)}

is surjective where 𝒟=loc.i=1pkerλi\mathcal{D}\underset{loc.}{=}\cap_{i=1}^{p}\ker\lambda^{i} and ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}}. This, in particular, implies that ωi\omega_{i} are non-degenerate for each i=1,,pi=1,\ldots,p.

Example 2.5.

A contact distribution is fat. In fact, contact distributions are precisely the fat distributions in the corank 1 case.

Since fatness is an open condition, fat distribution germs form an open set in the space of all germs of distributions of a fixed rank and corank. In general, fatness imposes strong numerical constraints on the rank and corank of the distribution.

Theorem 2.6 ([Ray68, Mon02]).

Suppose 𝒟\mathcal{D} is a rank kk-distribution on MM with dimM=n\dim M=n. If 𝒟\mathcal{D} is fat then the following numeric constraints hold.

  • kk is divisible by 22 and if k<n1k<n-1, then kk is divisible by 44

  • k(nk)+1k\geq(n-k)+1

  • The sphere Sk1S^{k-1} admits nkn-k linearly independent vector fields

Conversely, given any pair (k,n)(k,n) satisfying the above, there is a germ of fat distribution of type (k,n)(k,n).

Remark 2.7.

Though contact distributions happen to be generic in the corank 11 situation, fatness in higher corank is not a generic property. For corank 3\geq 3, this follows easily from the Determinacy theorem ([Mon02, pg. 65]). As for corank 22, we can define the signature of a (germ of) distribution of type (4,6)(4,6) ([Mon02, pg. 92]). We have two disjoint open classes of distributions germs of type (4,6)(4,6) : the elliptic and the hyperbolic type ([Zan15]). The elliptic type corresponds to fat distribution germs, whereas a typical example of a hyperbolic type is given by the product of two contact structures of type (2,3)(2,3) as considered in [D’A94, BH05].

It follows that the fat distributions are of the type (4n,4n+p)(4n,4n+p) when p>1p>1. We will now focus on p=2p=2.

2.3. Corank 22 Fat Distributions

Suppose 𝒟\mathcal{D} is a corank 22 fat distribution on MM. For simplicity, let us first assume that 𝒟\mathcal{D} is cotrivial. Hence there exist 11-forms λ1,λ2Ω1(M)\lambda^{1},\lambda^{2}\in\Omega^{1}(M) such that the quotient map λ:TMTM/𝒟M×2\lambda:TM\to TM/\mathcal{D}\cong M\times\mathbb{R}^{2} is given as, λ=(λ1,λ2)\lambda=(\lambda^{1},\lambda^{2}) and

𝒟=kerλ1kerλ2.\mathcal{D}=\ker\lambda^{1}\cap\ker\lambda^{2}.

Moreover, the curvature 22-form Ω:Λ2𝒟TM/𝒟M×2\Omega:\Lambda^{2}\mathcal{D}\to TM/\mathcal{D}\cong M\times\mathbb{R}^{2} is given by

Ω(X,Y)=(ω1(X,Y),ω2(X,Y)),\Omega(X,Y)=\big{(}\omega_{1}(X,Y),\omega_{2}(X,Y)\big{)},

where ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}}. Since 𝒟\mathcal{D} is fat, we have from 2.4 that ω1,ω2\omega_{1},\omega_{2} are nondegenerate. Therefore, we can define an automorphism A:𝒟𝒟A:\mathcal{D}\to\mathcal{D} by the following rule :

ω1(u,Av)=ω2(u,v),u,v𝒟.\omega_{1}(u,Av)=\omega_{2}(u,v),\forall u,v\in\mathcal{D}.

Explicitly, we have,

A=Iω11Iω2,A=-I_{\omega_{1}}^{-1}\circ I_{\omega_{2}},

where Iωi:𝒟𝒟I_{\omega_{i}}:\mathcal{D}\to\mathcal{D}^{*} is the induced isomorphism Iωi(v)=ιvωiI_{\omega_{i}}(v)=\iota_{v}\omega_{i}.

For any subspace V𝒟xV\subset\mathcal{D}_{x}, denote the symplectic complement of VV with respect to ωi\omega_{i} by ViV^{\perp_{i}}, i=1,2i=1,2.

Vi={w𝒟x|dλi(v,w)=0,vV}.V^{\perp_{i}}=\big{\{}w\in\mathcal{D}_{x}|d\lambda^{i}(v,w)=0,\forall v\in V\big{\}}.

It is easy to deduce that

V2=(AV)1,V1=A(V2).V^{\perp_{2}}=\big{(}AV\big{)}^{\perp_{1}},\quad V^{\perp_{1}}=A\big{(}V^{\perp_{2}}\big{)}.

For a general corank2\operatorname{corank}2 fat distribution 𝒟\mathcal{D}, not necessarily cotrivializable, the automorphism AA can only be defined locally, since it depends on the choice of annihilating forms for 𝒟\mathcal{D}. We observe a criteria for fatness, in the corank 22 situation.

Proposition 2.8.

Suppose 𝒟\mathcal{D} is a corank 22 distribution on MM defined locally by a pair of 1-forms λ1,λ2\lambda^{1},\lambda^{2}. Then 𝒟\mathcal{D} is fat if and only if the following conditions are satisfied:

  • ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}} is nondegenerate for i=1,2i=1,2.

  • The (local) automorphism A:𝒟𝒟A:\mathcal{D}\to\mathcal{D} relating ω1,ω2\omega_{1},\omega_{2} has no real eigenvalue.

Proof.

First suppose that 𝒟\mathcal{D} is strongly bracket generating at xx. This means that for any 0v𝒟x0\neq v\in\mathcal{D}_{x}, the map

Φv:𝒟x\displaystyle\Phi_{v}:\mathcal{D}_{x} TxM/𝒟x2\displaystyle\to T_{x}M/{\mathcal{D}_{x}}\cong\mathbb{R}^{2}
u\displaystyle u (ιvω1(u),ιvω2(u))\displaystyle\mapsto\big{(}\iota_{v}\omega_{1}(u),\iota_{v}\omega_{2}(u)\big{)}

is surjective. This, in particular, implies that ωi\omega_{i} are non-degenerate. On the other hand, if ω1\omega_{1} and ω2\omega_{2} are non-degenerate then

Φv\Phi_{v} is onto \Leftrightarrow corank(v1v2)=2\operatorname{corank}(v^{\perp_{1}}\cap v^{\perp_{2}})=2

Since v2=(Av)1v^{\perp_{2}}=(Av)^{\perp_{1}}, this is equivalent to AA having no real eigenvalue. This completes the proof.∎

Example 2.9.

A holomorphic 11-form Θ\Theta on a complex manifold MM with dimM=2n+1\dim_{\mathbb{C}}M=2n+1, is called holomorphic contact if it satisfies ΘdΘn0\Theta\wedge d\Theta^{n}\neq 0. By the holomorphic Darboux theorem ([AFL17]), we have holomorphic coordinates (z,x1,,xn,y1,,yn)(z,x_{1},\ldots,x_{n},y_{1},\ldots,y_{n}) on MM such that holomorphic contact form is given as, Θ=loc.dzj=1nyjdxj\Theta\underset{loc.}{=}dz-\sum_{j=1}^{n}y_{j}dx_{j}. If we identify, 2n+1\mathbb{C}^{2n+1} with 4n+2\mathbb{R}^{4n+2} and write z=z1+ιz2,xj=xj1+ιxj2,yj=yj1+ιyj2z=z_{1}+\iota z_{2},\,x_{j}=x_{j1}+\iota x_{j2},\,y_{j}=y_{j1}+\iota y_{j2}, then Θ\Theta can be expressed as Θ=λ1+ιλ2\Theta=\lambda^{1}+\iota\lambda^{2}, where

λ1=dz1j=1n(yj1dxj1yj2dxj2),λ2=dz2j=1n(yj2dxj1+yj1dxj2).\lambda^{1}=dz_{1}-\sum_{j=1}^{n}\big{(}y_{j1}dx_{j1}-y_{j2}dx_{j2}\big{)},\qquad\lambda^{2}=dz_{2}-\sum_{j=1}^{n}\big{(}y_{j2}dx_{j1}+y_{j1}dx_{j2}\big{)}.

This gives us a corank 22 distribution 𝒟=kerλ1kerλ2TM\mathcal{D}=\ker\lambda^{1}\cap\ker\lambda^{2}\subset TM, which is canonically isomorphic to the holomorphic contact subbundle kerΘT(1,0)M\ker\Theta\subset T^{(1,0)}M. We can explicitly write down a local frame 𝒟=SpanXj1,Xj2,Yj1,Yj2\mathcal{D}=\textrm{Span}\langle X_{j1},X_{j2},Y_{j1},Y_{j2}\rangle, where

Xj1=xj1+yj1z1+yj2z2,Xj2=xj2yj2z1+yj1z2,Yj1=yj1,Yj2=yj2.X_{j1}=\partial_{x_{j1}}+y_{j1}\partial_{z_{1}}+y_{j2}\partial_{z_{2}},\quad X_{j2}=\partial_{x_{j2}}-y_{j2}\partial_{z_{1}}+y_{j1}\partial_{z_{2}},\quad Y_{j1}=\partial_{y_{j1}},\quad Y_{j2}=\partial_{y_{j2}}.

Then the connecting automorphism A:𝒟𝒟A:\mathcal{D}\to\mathcal{D} defined by dλ1(u,Av)=dλ2(u,v)d\lambda^{1}(u,Av)=d\lambda^{2}(u,v) for each u,v𝒟u,v\in\mathcal{D}, satisfies the following equations :

AXj1=Xj2,AXj2=Xj1,AYj1=Yj2,AYj2=Yj1.AX_{j1}=-X_{j2},\quad AX_{j2}=X_{j1},\qquad AY_{j1}=-Y_{j2},\quad AY_{j2}=Y_{j1}.

In particular, we have A2=IdA^{2}=-Id, that is, AA induces a complex structure on 𝒟\mathcal{D}. Since AA has no real eigenvalue, by 2.8 the distribution 𝒟\mathcal{D} is fat.

Furthermore, the tangent bundle of MM splits as the direct sum TM=𝒟SpanZ1,Z2TM=\mathcal{D}\oplus\textrm{Span}\langle Z_{1},Z_{2}\rangle, where Zi=ziZ_{i}=\partial_{z_{i}} are two vector fields satisfying the relations below :

[Z1,Z2]=0,λi(Zj)=δij,ιZidλj|𝒟=0,i,j=1,2.[Z_{1},Z_{2}]=0,\quad\lambda^{i}(Z_{j})=\delta_{ij},\quad\iota_{Z_{i}}d\lambda^{j}|_{\mathcal{D}}=0,\qquad i,j=1,2.

Motivated by this, we consider the following.

Definition 2.10.

A corank 22 distribution 𝒟\mathcal{D} on MM is said to admit (local) Reeb directions Z1,Z2Z_{1},Z_{2}, if 𝒟=loc.kerλ1kerλ2\mathcal{D}\underset{loc.}{=}\ker\lambda^{1}\cap\ker\lambda^{2} and TM=𝒟SpanZ1,Z2TM=\mathcal{D}\oplus\textrm{Span}\langle Z_{1},Z_{2}\rangle such that,

  1.   (a)

    λ1(Z1)=1,λ1(Z2)=0\lambda^{1}(Z_{1})=1,\;\lambda^{1}(Z_{2})=0

  2.   (b)

    λ2(Z1)=0,λ2(Z2)=1\lambda^{2}(Z_{1})=0,\;\lambda^{2}(Z_{2})=1

  3.   (c)

    ιZidλj|𝒟=0\iota_{Z_{i}}d\lambda^{j}|_{\mathcal{D}}=0 for i,j=1,2i,j=1,2

  4.   (d)

    [Z1,Z2]=0[Z_{1},Z_{2}]=0

As observed in 2.9, the real distribution associated to the holomorphic contact structure, admits (local) Reeb directions.

Given any corank 22 fat distribution 𝒟\mathcal{D} on a manifold MM of dimension 4n+24n+2, we may find ([Ge92]) a coordinate system (x1,,x4n,z1,z2)(x_{1},\ldots,x_{4n},z_{1},z_{2}) and 11-forms

λi=dzij,kΓjkixjdxk+Ri,i=1,2,\lambda^{i}=dz_{i}-\sum_{j,k}\Gamma^{i}_{jk}x_{j}dx_{k}+R_{i},\quad i=1,2,

such that 𝒟=loc.kerλ1kerλ2\mathcal{D}\underset{loc.}{=}\ker\lambda^{1}\cap\ker\lambda^{2}. Here Ri=j=12fijdzj+j=14ngijdxjR_{i}=\sum_{j=1}^{2}f_{ij}dz_{j}+\sum_{j=1}^{4n}g_{ij}dx_{j} is a 11-form such that, fij,gijO(|x|2+|z|2)f_{ij},g_{ij}\in O(|x|^{2}+|z|^{2}), and {Γjki}\{\Gamma^{i}_{jk}\} constitute the structure constants of some nilpotent Lie algebra, known as the nilpotentization ([Mon02, Tan70]), associated to the distribution 𝒟\mathcal{D}. In particular Γjki=Γkji\Gamma^{i}_{jk}=-\Gamma^{i}_{kj}. Observe that, if we take fij=0f_{ij}=0 and gijg_{ij} to be functions of xkx_{k}’s only, then any such tuple of forms (λ1,λ2)(\lambda^{1},\lambda^{2}) above gives a corank 22 distribution, which admits local Reeb directions (z1,z2)(\partial_{z_{1}},\partial_{z_{2}}).

From the classification results of [CFS05], we see that the only possible Lie algebra that can arise as the nilpotentization of a corank 22 fat distribution on a 66 dimensional manifold is the complex Heisenberg Lie algebra.

Question 2.11.

Is every (germ of) corank 22 fat distribution on 6\mathbb{R}^{6}, which admits local Reeb directions, diffeomorphic to the germ of the distribution associated to a holomorphic contact structure as in 2.9?

For a general corank 2 fat distribution, the answer is clearly no. From a result of Montgomery ([Mon93]), it follows that a generic distribution germ of type (4,6)(4,6) cannot admit a local frame generating a finite dimensional Lie algebra. This differs from a holomorphic contact distribution, which does admit a local frame generating the complex Heisenberg Lie algebra (see 2.9). Since the set of germs of fat distributions of type (4,6)(4,6) is open, there are plenty of fat distributions, non-diffeomorphic to the contact holomorphic one. But it is not clear whether any of these fat distributions admits (local) Reeb directions.

Note that if the answer to the above question is in the affirmative, we can characterize germs of horizontal immersions, given by the 11-jet prolongation of holomorphic maps \mathbb{C}\to\mathbb{C}. But we suspect that the answer is negative as the 11-forms obtained in [Ge92] indicates the presence of function moduli. Let us now study the question of horizontal immersions.

3. Horizontal immersions in a corank 2 fat distribution

Suppose 𝒟\mathcal{D} is a corank 22 fat distribution on MM defined by a pair of 11-forms λ1,λ2\lambda^{1},\lambda^{2}. Hence ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}} are non-degenerate and the connecting homomorphism A:𝒟𝒟A:\mathcal{D}\to\mathcal{D} defined by

ω2(u,v)=ω1(u,Axv),u,v𝒟x,xM\omega_{2}(u,v)=\omega_{1}(u,A_{x}v),\quad\forall u,v\in\mathcal{D}_{x},\;x\in M

has no real eigenvalue. We further assume that the distribution 𝒟\mathcal{D} admits local Reeb directions.

Now for a fixed manifold Σ\Sigma, consider the partial differential operator,

𝔇:C(Σ,M)\displaystyle\mathfrak{D}:C^{\infty}(\Sigma,M) Ω1(Σ,2)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{2})
f\displaystyle f (fλ1,fλ2)\displaystyle\mapsto\big{(}f^{*}\lambda^{1},f^{*}\lambda^{2}\big{)}

The CC^{\infty}-solutions of 𝔇(f)=0\mathfrak{D}(f)=0 are precisely the 𝒟\mathcal{D}-horizontal maps since the derivative of ff maps TΣT\Sigma into 𝒟\mathcal{D}. Furthermore, horizontality implies that fdλ1=0=fdλ2f^{*}d\lambda^{1}=0=f^{*}d\lambda^{2}; hence dfx:TxΣ𝒟f(x)df_{x}:T_{x}\Sigma\to\mathcal{D}_{f(x)} is an isotropic map with respect to both the forms ω1\omega_{1} and ω2\omega_{2} on 𝒟\mathcal{D} for every xΣx\in\Sigma. Now, linearizing 𝔇\mathfrak{D} at an fC(Σ,M)f\in C^{\infty}(\Sigma,M) we get the linear differential operator 𝔏f\mathfrak{L}_{f} as follows :

𝔏f:ΓfTM\displaystyle\mathfrak{L}_{f}:\Gamma f^{*}TM Ω1(Σ,2)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{2})
\displaystyle\partial (d(λi)+fιdλi)i=1,2\displaystyle\mapsto\Big{(}d\big{(}\lambda^{i}\circ\partial\big{)}+f^{*}\iota_{\partial}d\lambda^{i}\Big{)}_{i=1,2}

Restricting 𝔏f\mathfrak{L}_{f} to Γf𝒟\Gamma f^{*}\mathcal{D} we have,

f:Γf𝒟\displaystyle\mathcal{L}_{f}:\Gamma f^{*}\mathcal{D} Ω1(Σ,2)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{2})
\displaystyle\partial (fιdλi)i=1,2\displaystyle\mapsto\big{(}f^{*}\iota_{\partial}d\lambda^{i}\big{)}_{i=1,2}

Observe that f\mathcal{L}_{f} is a C(Σ)C^{\infty}(\Sigma)-linear map and hence is induced by a bundle map f𝒟TΣ2f^{*}\mathcal{D}\to T^{*}\Sigma\otimes\mathbb{R}^{2}.

An horizontal immersion ff is said to be regular if this bundle map f\mathcal{L}_{f} is surjective; this is referred to as Ω\Omega-regularity in [Gro96]. One then gets that the operator 𝔇\mathfrak{D} is infinitesimally invertible over the set of regular maps and an appeal to Gromov’s general theorems ([Gro86]). But for such a map ff to exist, i.e, for the existence of a regular, common isotropic subspace V𝒟V\subset\mathcal{D}, we must have the inequality ([Gro96])

rk𝒟dimΣ2dimΣ,\operatorname{rk}\mathcal{D}-\dim\Sigma\geq 2\dim\Sigma,

as common isotropic subspaces are necessarily in the kernel of the map f\mathcal{L}_{f}. We will be focusing on 𝒟\mathcal{D}-horizontal immersions of discs 𝔻2\mathbb{D}^{2} in 66-dimensional manifold, where rk𝒟=4\operatorname{rk}\mathcal{D}=4. Clearly, rk𝒟=46=3.2=3dimΣ\operatorname{rk}\mathcal{D}=4\not\geq 6=3.2=3\dim\Sigma. Hence there is no possibility of an regular horizontal map 𝔻2M\mathbb{D}^{2}\to M to exist and so Gromov’s method does not apply directly.

3.1. Inversion of 𝔏f\mathfrak{L}_{f} at 𝒟\mathcal{D}-horizontal Immersions

We now denote Σ=𝔻2\Sigma=\mathbb{D}^{2} and M=6M=\mathbb{R}^{6}. Suppose 𝒟=kerλ1kerλ2\mathcal{D}=\ker\lambda^{1}\cap\ker\lambda^{2} is a given corank 22 fat distribution, which admits local Reeb directions. Since Σ\Sigma is a compact manifold with boundary, we have (see 5.2),

Observation 3.1.

The spaces Γ(fTM)\Gamma(f^{*}TM) and Ω1(Σ,2)\Omega^{1}(\Sigma,\mathbb{R}^{2}) are tame Fréchet spaces.

As before we have the linearization map, 𝔏f:Γf𝒟Ω1(Σ,2)\mathfrak{L}_{f}:\Gamma f^{*}\mathcal{D}\to\Omega^{1}(\Sigma,\mathbb{R}^{2}). Since 𝔏f\mathfrak{L}_{f} is a linear partial differential operator of order 11, we have (see 5.4),

Observation 3.2.

𝔏f\mathfrak{L}_{f} is a tame linear map of order 11.

This sets the problem into the framework of differential operators between Fréchet spaces for studying the existence of local inversion. We refer to the appendix (section 5) for relevant details. We first prove the following result.

Proposition 3.3.

If ff is a smooth horizontal immersion, then 𝔏f\mathfrak{L}_{f} admits a tame inverse 𝔐f\mathfrak{M}_{f}.

Note that we are assuming the existence of 𝒟\mathcal{D}-horizontal immersions in the above proposition. In fact, in the next section, we shall prove the inversion for an open set of maps (Theorem 3.9). Let us first prove the following.

Lemma 3.4.

If V𝒟xV\subset\mathcal{D}_{x} is common isotropic with respect to ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}} and dimV=2\dim V=2, then V=AVV=AV.

Proof.

Since VV is common isotropic,

VV1V2=(V+AV)1dim(V+AV)dimV=2V\subset V^{\perp_{1}}\cap V^{\perp_{2}}=(V+AV)^{\perp_{1}}\;\Rightarrow\;\dim(V+AV)^{\perp}\geq\dim V=2

and so, dim(V+AV)dim𝒟x2=2\dim(V+AV)\leq\dim\mathcal{D}_{x}-2=2. On the other hand, dim(V+AV)dimV=2\dim(V+AV)\geq\dim V=2. Hence, dim(V+AV)=2=dimV\dim(V+AV)=2=\dim V, which is only possible if V=AVV=AV. ∎

Proposition 3.5.

If ff is a smooth 𝒟\mathcal{D}-horizontal immersion, given any (P,Q)Ω1(Σ,2)(P,Q)\in\Omega^{1}(\Sigma,\mathbb{R}^{2}), the equation 𝔏f()=(P,Q)\mathfrak{L}_{f}(\partial)=(P,Q) admits a unique solution =𝔐f(P,Q)\partial=\mathfrak{M}_{f}(P,Q), subject to a boundary condition. The process of obtaining the solution depends on a choice of complex structure JJ on 𝒟\mathcal{D}.

Proof.

First, choose an almost complex structure JJ on 𝒟\mathcal{D}, compatible with ω1=dλ1|𝒟\omega_{1}=d\lambda^{1}|_{\mathcal{D}}, i.e, the assignment (u,v)ω1(u,Jv)(u,v)\mapsto\omega_{1}(u,Jv) is a nondegenerate symmetric form. Such a JJ always exists ([CdS01, pg. 86]). Clearly, JAJ\neq A.

Since ff is 𝒟\mathcal{D}-horizontal we have,

fλi=0fdλi=0.f^{*}\lambda^{i}=0\Rightarrow f^{*}d\lambda^{i}=0.

Thus, for σΣ\sigma\in\Sigma, Imdfσ\operatorname{Im}df_{\sigma} is common isotropic with respect to both ωi=dλi|𝒟\omega_{i}=d\lambda^{i}|_{\mathcal{D}}. In particular, Imdfσ\operatorname{Im}df_{\sigma} is JJ-totally real, since JJ is ω1\omega_{1}-compatible. Also since ff is an immersion, dimImdfσ=2\dim\operatorname{Im}df_{\sigma}=2. Then by 3.4, we have that

A(Imdfσ)=Imdfσ,for σΣ.A(\operatorname{Im}df_{\sigma})=\operatorname{Im}df_{\sigma},\quad\text{for $\sigma\in\Sigma$.}

Let us denote, X=f(x),Y=f(y)X=f_{*}(\partial_{x}),Y=f_{*}(\partial_{y}), where x,y\partial_{x},\partial_{y} are the coordinate vector fields on Σ=𝔻2\Sigma=\mathbb{D}^{2}. We thus have

SpanAX,AY=SpanX,Y.\textrm{Span}\langle AX,AY\rangle=\textrm{Span}\langle X,Y\rangle.

Hence, AA restricts to an automorphism on SpanX,Y\textrm{Span}\langle X,Y\rangle :

A0=A|SpanX,Y.A_{0}=A|_{\textrm{Span}\langle X,Y\rangle}.

Let us write,

(*) AX=pX+qY,AY=rX+sYAX=pX+qY,\quad AY=rX+sY

for some functions p,q,r,sC(Σ)p,q,r,s\in C^{\infty}(\Sigma). Then we have that A0=(pqrs)A_{0}=\begin{pmatrix}p&q\\ r&s\end{pmatrix} with respect to the basis (X,Y)(X,Y). Since AA has no real eigenvalue, A0A_{0} also has no real eigenvalue. This means that the characteristic polynomial

λ2(p+s)λ+(psqr)\lambda^{2}-(p+s)\lambda+(ps-qr)

of A0A_{0} has negative discriminant, i.e.,

(p+s)24(psqr)=(ps)2+4qr<0.(p+s)^{2}-4(ps-qr)=(p-s)^{2}+4qr<0.

Now let us consider the equation

𝔏f()=(P,Q),\mathfrak{L}_{f}(\partial)=(P,Q),

where P,QΩ1(Σ)P,Q\in\Omega^{1}(\Sigma). We write

=0+aZ1+bZ2,\partial=\partial_{0}+aZ_{1}+bZ_{2},

where 0f𝒟\partial_{0}\in f^{*}\mathcal{D} and Z1,Z2Z_{1},Z_{2} are the Reeb directions associated to (λ1,λ2)(\lambda^{1},\lambda^{2}), pulled back along ff. We then have,

𝔏f()=(da+fι0dλ1,db+fι0dλ2)\mathfrak{L}_{f}(\partial)=\Big{(}da+f^{*}\iota_{\partial_{0}}d\lambda^{1},db+f^{*}\iota_{\partial_{0}}d\lambda^{2}\Big{)}

Also let us write

P=P1dx+P2dy,Q=Q1dx+Q2dy.P=P_{1}dx+P_{2}dy,\quad Q=Q_{1}dx+Q_{2}dy.

Evaluating both sides on x,y\partial_{x},\partial_{y} and using properties (a), (b), (c) of Defn 2.10, we have the system,

(1) {xa+dλ1(0,X)=P1ya+dλ1(0,Y)=P2\displaystyle\left\{\quad\begin{aligned} \partial_{x}a+d\lambda^{1}(\partial_{0},X)&=P_{1}\\ \partial_{y}a+d\lambda^{1}(\partial_{0},Y)&=P_{2}\end{aligned}\right.
(2) {xb+dλ2(0,X)=Q1yb+dλ2(0,Y)=Q2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b+d\lambda^{2}(\partial_{0},X)&=Q_{1}\\ \partial_{y}b+d\lambda^{2}(\partial_{0},Y)&=Q_{2}\end{aligned}\right.

Now using (*3.1) we have,

dλ2(0,X)=dλ1(0,AX)=pdλ1(0,X)+qdλ1(0,Y)d\lambda^{2}(\partial_{0},X)=d\lambda^{1}(\partial_{0},AX)=p\,d\lambda^{1}(\partial_{0},X)+q\,d\lambda^{1}(\partial_{0},Y)
dλ2(0,Y)=dλ1(0,AY)=rdλ1(0,X)+sdλ1(0,Y)d\lambda^{2}(\partial_{0},Y)=d\lambda^{1}(\partial_{0},AY)=r\,d\lambda^{1}(\partial_{0},X)+s\,d\lambda^{1}(\partial_{0},Y)

This transforms (2) into the following system of PDEs :

(22^{\prime}) {xb+pdλ1(0,X)+qdλ1(0,Y)=Q1yb+rdλ1(0,X)+sdλ1(0,Y)=Q2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b+p\,d\lambda^{1}(\partial_{0},X)+q\,d\lambda^{1}(\partial_{0},Y)&=Q_{1}\\ \partial_{y}b+r\,d\lambda^{1}(\partial_{0},X)+s\,d\lambda^{1}(\partial_{0},Y)&=Q_{2}\end{aligned}\right.

Using (1) we eliminate 0\partial_{0} from (22^{\prime}) and get

(2′′2^{\prime\prime}) {xbpxaqya=Q1pP1qP2ybrxasya=Q2rP1sP2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b-p\partial_{x}a-q\partial_{y}a&=Q_{1}-pP_{1}-qP_{2}\\ \partial_{y}b-r\partial_{x}a-s\partial_{y}a&=Q_{2}-rP_{1}-sP_{2}\end{aligned}\right.

Since (ps)2+4qr<0(p-s)^{2}+4qr<0, the system of PDEs given by (2′′2^{\prime\prime}) is elliptic. Hence, the Dirichlet problem (2′′2^{\prime\prime}) with the boundary condition

(3) a|Σ=a0,b|Σ=b0,\displaystyle a|_{\partial\Sigma}=a_{0},\quad b|_{\partial\Sigma}=b_{0},

will have a unique solution

(a,b)=Mf(P,Q,a0,b0).(a,b)=M_{f}(P,Q,a_{0},b_{0}).

Now consider an auxiliary system of equations :

(4) {dλ1(0,JX)=0dλ1(0,JY)=0\displaystyle\left\{\quad\begin{aligned} d\lambda^{1}(\partial_{0},JX)&=0\\ d\lambda^{1}(\partial_{0},JY)&=0\end{aligned}\right.

Then using the solution (a,b)=Mf(P,Q,a0,b0)(a,b)=M_{f}(P,Q,a_{0},b_{0}), we get from (1), (4), the system

(5) {dλ1(0,X)=P1xadλ1(0,Y)=P2yadλ1(0,JX)=0dλ1(0,JY)=0\displaystyle\left\{\quad\begin{aligned} d\lambda^{1}(\partial_{0},X)&=P_{1}-\partial_{x}a\\ d\lambda^{1}(\partial_{0},Y)&=P_{2}-\partial_{y}a\\ d\lambda^{1}(\partial_{0},JX)&=0\\ d\lambda^{1}(\partial_{0},JY)&=0\end{aligned}\right.

Since Imdfσ\operatorname{Im}df_{\sigma} is JJ-totally real, (X,Y,JX,JY)(X,Y,JX,JY) is a local framing of 𝒟\mathcal{D}, and since dλ1|𝒟d\lambda^{1}|_{\mathcal{D}} is nondegenerate, (5) can be uniquely solved for 0\partial_{0}. Thus, 𝔏f()=(P,Q)\mathfrak{L}_{f}(\partial)=(P,Q) has a unique solution

=𝔐f(P,Q,a0,b0)\partial=\mathfrak{M}_{f}(P,Q,a_{0},b_{0})

subject to satisfying the auxiliary system (4) and the boundary condition (3). ∎

Remark 3.6.

It can be easily seen from 2.9 that for our model case (M,𝒟)(M,\mathcal{D}) of holomorphic contact structure, we have A=J|𝒟A=-J|_{\mathcal{D}}, where JJ is the (integrable) almost complex structure on MM. 3.4 can then be interpreted as follows : common isotropic 22-subspaces of 𝒟\mathcal{D} are complex subspaces. In particular, the left hand side of (2′′2^{\prime\prime}) can then be compared to the usual Cauchy-Riemann equations for the tuple of functions (a,b)(a,b) on Σ\Sigma.

We can now prove 3.3

Proof of 3.3.

From 3.5 we have that 𝔏f\mathfrak{L}_{f} admits unique solution 𝔐f\mathfrak{M}_{f}, whenever ff is a 𝒟\mathcal{D}-horizontal immersion. As in 3.5, MfM_{f} is obtained as a solution to a Dirichlet problem and hence it is tame (see 5.4 (2)). Then 𝔐f\mathfrak{M}_{f} is obtained from MfM_{f} by solving a linear system, which is again tame. Hence the inverse 𝔐f\mathfrak{M}_{f} is tame, as composition of tame maps is tame. ∎

Remark 3.7.

In fact, the operator 𝔐f\mathfrak{M}_{f} above is tame of degree 11. Indeed, the proof of tameness for elliptic boundary value problems ([Ham82, pg. 161]) suggests that MfM_{f} is tame of degree 0. Next, to get 𝔐f\mathfrak{M}_{f} from MfM_{f}, the linear system (5) involves taking first order differentials and hence it is tame of degree 11. Thus, 𝔐f\mathfrak{M}_{f} is tame of degree 11.

3.2. Local Inversion of 𝔇\mathfrak{D}

From 3.3 we see that the linearization 𝔏f\mathfrak{L}_{f} admits a right inverse 𝔐f\mathfrak{M}_{f}, provided ff is a 𝒟\mathcal{D}-horizontal immersion. But in order to apply the Implicit Function Theorem due to Hamilton (Theorem 5.9), we need to show that there is an open set of maps 𝔘C(Σ,M)\mathfrak{U}\subset C^{\infty}(\Sigma,M) such that the family {𝔏f|f𝔘}\{\mathfrak{L}_{f}\;|\;f\in\mathfrak{U}\} admits a smooth tame inverse. We now identify this set 𝔘\mathfrak{U}.

We first restrict ourselves to a collection 𝔘0\mathfrak{U}_{0} of maps f:ΣMf:\Sigma\to M satisfying the following conditions :

  • ff is an immersion, and

  • Imdf\operatorname{Im}df is transverse to SpanZ1,Z2\textrm{Span}\langle Z_{1},Z_{2}\rangle.

This collection 𝔘0C(Σ,M)\mathfrak{U}_{0}\subset C^{\infty}(\Sigma,M) is clearly open, since it is defined by open conditions. Now we have a canonical projection

π𝒟:TM=𝒟SpanZ1,Z2𝒟.\pi_{\mathcal{D}}:TM=\mathcal{D}\oplus\textrm{Span}\langle Z_{1},Z_{2}\rangle\to\mathcal{D}.

For any f𝔘0f\in\mathfrak{U}_{0} we see that the image π𝒟(Imdf)\pi_{\mathcal{D}}(\operatorname{Im}df) has dimension 22 at each point of Σ\Sigma. Let us choose an almost complex structure J:𝒟𝒟J:\mathcal{D}\to\mathcal{D}, compatible with dλ1|𝒟d\lambda^{1}|_{\mathcal{D}}, as in 3.5. Then the set

{(X,Y)Fr2𝒟|V=SpanX,Y is J-totally real}\Big{\{}(X,Y)\in\mathrm{Fr}_{2}\mathcal{D}\;\Big{|}\;\text{$V=\textrm{Span}\langle X,Y\rangle$ is $J$-totally real}\Big{\}}

is open in the 22-frame bundle Fr2𝒟\mathrm{Fr}_{2}\mathcal{D}, since the totally real condition VJV=0V\cap JV=0 is open. For any such tuple (X,Y)(X,Y) we have the framing (X,Y,JX,JY)(X,Y,JX,JY) of 𝒟\mathcal{D} and we can write

A=(A11A12A21A22)A=\begin{pmatrix}A_{11}&A_{12}\\ A_{21}&A_{22}\end{pmatrix}

with respect to this basis. Let 𝒪xFr2𝒟x\mathcal{O}_{x}\subset\mathrm{Fr}_{2}\mathcal{D}_{x} be the set of those (X,Y)Fr2𝒟x(X,Y)\in\mathrm{Fr}_{2}\mathcal{D}_{x} such that,

  • V=SpanX,YV=\textrm{Span}\langle X,Y\rangle is JJ-totally real, and

  • The matrix A11A_{11} as above is negative definite.

Since both are open conditions, we see that 𝒪x\mathcal{O}_{x} is open in Fr2𝒟x\mathrm{Fr}_{2}\mathcal{D}_{x}.

We now define,

Definition 3.8.

A map f:ΣMf:\Sigma\to M is said to be admissible if it satisfies the following.

  • f𝔘0f\in\mathfrak{U}_{0}, i.e, ff is an immersion with ImdfSpanZ1,Z2\operatorname{Im}df\pitchfork\textrm{Span}\langle Z_{1},Z_{2}\rangle.

  • Imdfσ=Spanfx,fyπ𝒟1(𝒪f(σ))\operatorname{Im}df_{\sigma}=\textrm{Span}\langle f_{*}\partial_{x},f_{*}\partial_{y}\rangle\in\pi_{\mathcal{D}}^{-1}\big{(}\mathcal{O}_{f(\sigma)}\big{)} for each σΣ\sigma\in\Sigma.

Denote by 𝔘C(Σ,M)\mathfrak{U}\subset C^{\infty}(\Sigma,M) the set of admissible maps.

In fact we have defined an open relation 𝒜J1(Σ,M)\mathcal{A}\subset J^{1}(\Sigma,M) such that 𝔘\mathfrak{U} is exactly the smooth holonomic solutions of 𝒜\mathcal{A}. Since 𝒜\mathcal{A} is an open relation, we have that 𝔘\mathfrak{U} is open in C(Σ,M)C^{\infty}(\Sigma,M). It is apparent that any 𝒟\mathcal{D}-horizontal immersion is admissible. We now prove the following.

Theorem 3.9.

The linearization 𝔏f\mathfrak{L}_{f} admits a smooth tame inverse 𝔐f\mathfrak{M}_{f} for every f𝔘f\in\mathfrak{U}.

Proof.

Suppose f𝔘f\in\mathfrak{U}. We have Imdf=Spanfx,fy\operatorname{Im}df=\textrm{Span}\langle f_{*}\partial_{x},f_{*}\partial_{y}\rangle. Let us write

fx=X+a1Z1+a2Z2,fy=Y+b1Z1+b2Z2f_{*}\partial_{x}=X+a_{1}Z_{1}+a_{2}Z_{2},\quad f_{*}\partial_{y}=Y+b_{1}Z_{1}+b_{2}Z_{2}

where X=π𝒟(fx),Y=π𝒟(fy)X=\pi_{\mathcal{D}}(f_{*}\partial_{x}),Y=\pi_{\mathcal{D}}(f_{*}\partial_{y}). By assumption (X,Y)Fr2𝒟(X,Y)\in Fr_{2}\mathcal{D} so that, (X,Y,JX,JY)(X,Y,JX,JY) is a frame of 𝒟\mathcal{D}. Hence we can write

(6) {AX=pX+qY+pJX+qJYAY=rX+sY+rJX+sJY\displaystyle\left\{\quad\begin{aligned} AX&=pX+qY+p^{\prime}JX+q^{\prime}JY\\ AY&=rX+sY+r^{\prime}JX+s^{\prime}JY\end{aligned}\right.

The matrix of AA has the form

(prqsprqs)\begin{pmatrix}p&r&*&*\\ q&s&*&*\\ p^{\prime}&r^{\prime}&*&*\\ q^{\prime}&s^{\prime}&*&*\end{pmatrix}

and by the hypothesis on 𝔘\mathfrak{U}, A11=(pqrs)A_{11}=\begin{pmatrix}p&q\\ r&s\end{pmatrix} is negative definite, which is equivalent to

(ps)2+4qr<0.(p-s)^{2}+4qr<0.

Now, we wish to solve 𝔏f()=(P,Q)\mathfrak{L}_{f}(\partial)=(P,Q), as we did in 3.5, where

𝔏f:ΓfTM\displaystyle\mathfrak{L}_{f}:\Gamma f^{*}TM Ω1(Σ,2)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{2})
\displaystyle\partial (d(λi)+fιdλi)i=1,2\displaystyle\mapsto\Big{(}d\big{(}\lambda^{i}\circ\partial\big{)}+f^{*}\iota_{\partial}d\lambda^{i}\Big{)}_{i=1,2}

Let =0+aZ1+bZ2\partial=\partial_{0}+aZ_{1}+bZ_{2}, where 0f𝒟\partial_{0}\in f^{*}\mathcal{D}. Since [Z1,Z2]=0[Z_{1},Z_{2}]=0 (by (d) of Defn 2.10), we have,

dλ1(Z1,Z2)=Z1(λ1(Z2))Z2(λ1(Z1))λ1([Z1,Z2])=Z1(0)Z2(1)λ1(0)=0.d\lambda^{1}(Z_{1},Z_{2})=Z_{1}(\lambda^{1}(Z_{2}))-Z_{2}(\lambda^{1}(Z_{1}))-\lambda^{1}([Z_{1},Z_{2}])=Z_{1}(0)-Z_{2}(1)-\lambda^{1}(0)=0.

and similarly, dλ2(Z1,Z2)=0d\lambda^{2}(Z_{1},Z_{2})=0. Hence,

dλ1(,fx)=dλ1(0+aZ1+bZ2,X+a1Z1+a2Z2)=dλ1(0,X)d\lambda^{1}(\partial,f_{*}\partial_{x})=d\lambda^{1}(\partial_{0}+aZ_{1}+bZ_{2},X+a_{1}Z_{1}+a_{2}Z_{2})=d\lambda^{1}(\partial_{0},X)

and similarly the remaining ones. Thus, we get a system as before :

(7) {xa+dλ1(0,X)=P1ya+dλ1(0,Y)=P2\displaystyle\left\{\quad\begin{aligned} \partial_{x}a+d\lambda^{1}(\partial_{0},X)&=P_{1}\\ \partial_{y}a+d\lambda^{1}(\partial_{0},Y)&=P_{2}\\ \end{aligned}\right.
(8) {xb+dλ1(0,AX)=Q1yb+dλ1(0,AY)=Q2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b+d\lambda^{1}(\partial_{0},AX)&=Q_{1}\\ \partial_{y}b+d\lambda^{1}(\partial_{0},AY)&=Q_{2}\end{aligned}\right.

We add the linear equations

(9) dλ1(0,JX)=0=dλ1(0,JY)\displaystyle d\lambda^{1}(\partial_{0},JX)=0=d\lambda^{1}(\partial_{0},JY)

to (7),(8). Then using (6) and (9), the system (8) becomes

(88^{\prime}) {xb+pdλ1(0,X)+qdλ1(0,Y)=Q1yb+rdλ1(0,X)+sdλ1(0,Y)=Q2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b+p\,d\lambda^{1}(\partial_{0},X)+q\,d\lambda^{1}(\partial_{0},Y)&=Q_{1}\\ \partial_{y}b+r\,d\lambda^{1}(\partial_{0},X)+s\,d\lambda^{1}(\partial_{0},Y)&=Q_{2}\end{aligned}\right.

Using (7) we can eliminate 0\partial_{0} in (88^{\prime}) and get

(8′′8^{\prime\prime}) {xbpxaqya=Q1pP1qP2ybrxasya=Q2rP1sP2\displaystyle\left\{\quad\begin{aligned} \partial_{x}b-p\partial_{x}a-q\partial_{y}a&=Q_{1}-pP_{1}-qP_{2}\\ \partial_{y}b-r\partial_{x}a-s\partial_{y}a&=Q_{2}-rP_{1}-sP_{2}\end{aligned}\right.

Since (ps)2+4qr<0(p-s)^{2}+4qr<0, we have that (8′′8^{\prime\prime}) is elliptic. Hence given any arbitrary boundary condition a|Σ=a0,b|Σ=b0a|_{\partial\Sigma}=a_{0},\;b|_{\partial\Sigma}=b_{0}, we have the unique solution

(a,b)=Mf(P,Q,a0,b0).(a,b)=M_{f}(P,Q,a_{0},b_{0}).

Then as done in 3.5, we obtain unique solution

=𝔐f(P,Q,a0,b0)\partial=\mathfrak{M}_{f}(P,Q,a_{0},b_{0})

to the system given by (7), (8) and (9). Thus whenever f𝔘f\in\mathfrak{U}, we have a solution 𝔐f\mathfrak{M}_{f} for the linearized equation 𝔏f=(P,Q)\mathfrak{L}_{f}=(P,Q). As argued in the proof of 3.3, both 𝔏f\mathfrak{L}_{f} and 𝔐f\mathfrak{M}_{f} are tame operators. ∎

Since 𝔏f\mathfrak{L}_{f} is surjective for every f𝔘f\in\mathfrak{U} and the family of right inverses 𝔐:𝔘×Ω1(Σ,2)C(Σ,M)\mathfrak{M}:\mathfrak{U}\times\Omega^{1}(\Sigma,\mathbb{R}^{2})\to C^{\infty}(\Sigma,M) is a smooth tame map we obtain the following by Theorem 5.9.

Theorem 3.10.

The operator 𝔇\mathfrak{D} restricted to 𝔘\mathfrak{U} is locally right invertible. Given any f0𝔘f_{0}\in\mathfrak{U}, there exists an open neighborhood UU of f0f_{0} and a smooth tame map 𝔇f01:𝔇(U)U\mathfrak{D}_{f_{0}}^{-1}:\mathfrak{D}(U)\to U such that 𝔇𝔇f01=Id\mathfrak{D}\circ\mathfrak{D}_{f_{0}}^{-1}=\textrm{Id}.

The proof of the Implicit Function Theorem, in fact, implies that there exists a positive integer r0r_{0} such that the following holds true.

Theorem 3.11.

Let f0𝔘f_{0}\in\mathfrak{U} and g0=𝔇(f0)g_{0}=\mathfrak{D}(f_{0}). Let ϵ>0\epsilon>0 be any positive number. Then there exists a δ>0\delta>0 and an integer r0r_{0}, such that for αr0\alpha\geq r_{0} and for every gΩ1(Σ,2)g\in\Omega^{1}(\Sigma,\mathbb{R}^{2}) with |g|α<δ|g|_{\alpha}<\delta there is an f=𝔇f01(g0+g)𝔘f=\mathfrak{D}_{f_{0}}^{-1}(g_{0}+g)\in\mathfrak{U} satisfying the following conditions :

𝔇(f)=g0+g and |ff0|α+2<ϵ.\mathfrak{D}(f)=g_{0}+g\quad\text{ and }\quad|f-f_{0}|_{\alpha+2}<\epsilon.

4. Existence of Horizontal Germs and the Local hh-Principle

Since we are only interested in germs, without loss of generality, we assume that M=6M=\mathbb{R}^{6} and Σ=2\Sigma=\mathbb{R}^{2}. Suppose, we have a corank 22 fat distribution 𝒟\mathcal{D} on MM, which admits Reeb directions (Defn 2.10). Consider the relation 𝒜J1(Σ,M)\mathcal{A}\subset J^{1}(\Sigma,M), as in section 3, such that the set of admissible maps (Defn 3.8) 𝔘\mathfrak{U} are exactly the smooth holonomic sections of 𝒜\mathcal{A}, i.e, we have 𝔘=Sol\mathfrak{U}=\operatorname{Sol}\mathcal{R}. We have shown that the operator

𝔇:𝔘C(Σ,M)\displaystyle\mathfrak{D}:\mathfrak{U}\subset C^{\infty}(\Sigma,M) Ω1(Σ,2)\displaystyle\to\Omega^{1}(\Sigma,\mathbb{R}^{2})
f\displaystyle f (fλ1,fλ2)\displaystyle\mapsto\big{(}f^{*}\lambda^{1},f^{*}\lambda^{2}\big{)}

is locally invertible over 𝔘\mathfrak{U}.

Now, following Gromov ([Gro86]), we can get the (parametric) local h-principle. One crucial thing to observe is that the inversion of 𝔇\mathfrak{D} as we have obtained, does not conform to the notion of locality as considered by Gromov ([Gro86, pg. 117-118]). But we observe that the proof of the local hh-principle goes through, without the locality property of 𝒟1\mathcal{D}^{-1}. For the sake of completeness, we reproduce the proof following Gromov.

Definition 4.1.

([Gro86, pg. 118]) A germ f:ΣMf:\Sigma\to M at σΣ\sigma\in\Sigma is called an infinitesimal solution of order α\alpha of 𝔇(f)=0\mathfrak{D}(f)=0 if

j𝔇(f)α(σ)=0j^{\alpha}_{\mathfrak{D}(f)}(\sigma)=0

i.e, the section 𝔇(f)\mathfrak{D}(f) has zero αth\alpha^{\text{th}}-jet at the point σ\sigma.

Observe that, since 𝔇\mathfrak{D} has order 11, the property that ff is an infinitesimal solution of order α\alpha, only depends on the jet jfα+1(σ)j^{\alpha+1}_{f}(\sigma). Consider the relation α=α(𝔇,0,𝒜)Jr+1(Σ,M)\mathcal{R}_{\alpha}=\mathcal{R}_{\alpha}(\mathfrak{D},0,\mathcal{A})\subset J^{r+1}(\Sigma,M) consisting of jets jfα+1(σ)j^{\alpha+1}_{f}(\sigma) represented by CC^{\infty} germs f:ΣMf:\Sigma\to M at σ\sigma, so that

j𝔇(f)α(σ)=0andjfα+1(σ)𝒜.j^{\alpha}_{\mathfrak{D}(f)}(\sigma)=0\quad\text{and}\quad j^{\alpha+1}_{f}(\sigma)\in\mathcal{A}.

Then one sees that Cα+1C^{\alpha+1} holonomic sections of α\mathcal{R}_{\alpha} are exactly the admissible Cα+1C^{\alpha+1}-solutions of 𝔇=0\mathfrak{D}=0. In particular, for α0\alpha\geq 0, the CC^{\infty}-solutions α\mathcal{R}_{\alpha} are all same, namely the CC^{\infty}-solutions of 𝔇=0\mathfrak{D}=0 which are admissible, i.e, we have

Solα\operatorname{Sol}\mathcal{R}_{\alpha} is the set of 𝒟\mathcal{D}-horizontal immersions, for any α0\alpha\geq 0.

We then prove the following.

Theorem 4.2.

If α\alpha is sufficiently large, then for any jet jfα(σ)αj^{\alpha}_{f}(\sigma)\in\mathcal{R}_{\alpha}, represented by some f:𝒪p(σ)Mf:\operatorname{\mathcal{O}p}(\sigma)\to M, we have a homotopy ft:𝒪p(σ)Mf_{t}:\operatorname{\mathcal{O}p}(\sigma)\to M, such that f0=ff_{0}=f on some 𝒪p(σ)\operatorname{\mathcal{O}p}(\sigma) and f1f_{1} is a 𝒟\mathcal{D}-horizontal admissible solution, i.e, 𝔇(f1)=0\mathfrak{D}(f_{1})=0. Furthermore, the jet jftα+1(σ)j^{\alpha+1}_{f_{t}}(\sigma) belongs to α\mathcal{R}_{\alpha}, for all t[0,1]t\in[0,1].

Proof.

Suppose ff is defined on an open ball VΣV\subset\Sigma about σ\sigma. Since f𝔘f\in\mathfrak{U} and 𝒜\mathcal{A} is open, we can get a neighborhood V0V_{0} of σ\sigma, such that σV0V\sigma\in V_{0}\subset V and f|V0f|_{V_{0}} is a solution of 𝒜\mathcal{A}. In other words, f|V0f|_{V_{0}} is admissible. Denote, g0=𝔇(f|V0)g_{0}=\mathfrak{D}(f|_{V_{0}}).

Since jfα+1(σ)αj^{\alpha+1}_{f}(\sigma)\in\mathcal{R}_{\alpha}, we have jg0α(σ)=j𝔇(f)α(σ)=0j^{\alpha}_{g_{0}}(\sigma)=j^{\alpha}_{\mathfrak{D}(f)}(\sigma)=0. Hence, for any given ϵ>0\epsilon>0, there exists a neighborhood WV0W\subset V_{0} of σ\sigma such that |g0|α<ϵ|g_{0}|_{\alpha}<\epsilon on WW. We can get some gϵg_{\epsilon} on V0V_{0} so that,

  • gϵ=g0g_{\epsilon}=-g_{0} on some neighborhood WV0W\subset V_{0} of σ\sigma, and

  • gϵg_{\epsilon} is ϵ\epsilon-small in CαC^{\alpha}-norm, i.e, |gϵ|α<ϵ|g_{\epsilon}|_{\alpha}<\epsilon on V0V_{0}.

Now let us apply Theorem 3.11 for the domain V0V_{0}. Since y0:=f|V0y_{0}:=f|_{V_{0}} is admissible, we have that 𝔇y0\mathfrak{D}_{y_{0}} admits a local inverse. In particular, there exists some ϵ,δ>0\epsilon,\delta>0 such that for any |g|α<ϵ|g|_{\alpha}<\epsilon we have unique yy such that 𝔇(y)=𝔇(y0)+g\mathfrak{D}(y)=\mathfrak{D}(y_{0})+g and |yy0|α+1<δ|y-y_{0}|_{\alpha+1}<\delta. Here we require that α\alpha to be sufficiently large. Now, in particular, for this ϵ=ϵ(y0,α)\epsilon=\epsilon(y_{0},\alpha), we can get WW and gϵg_{\epsilon} as above. Since,

|tgϵ|α<tϵ<ϵ,for t[0,1],|tg_{\epsilon}|_{\alpha}<t\epsilon<\epsilon,\quad\text{for $t\in[0,1]$,}

we have unique solutions

ft=𝔇y01(tgϵ),f_{t}=\mathfrak{D}_{y_{0}}^{-1}(tg_{\epsilon}),

over V0V_{0}, satisfying |fty0|α+1<ϵ|f_{t}-y_{0}|_{\alpha+1}<\epsilon for t[0,1]t\in[0,1]. Now,

𝔇(ft)=𝔇(y0)+tgϵ=𝔇(f|V0)+tgϵ=g0+tgϵ.\mathfrak{D}(f_{t})=\mathfrak{D}(y_{0})+tg_{\epsilon}=\mathfrak{D}(f|_{V_{0}})+tg_{\epsilon}=g_{0}+tg_{\epsilon}.

In particular, we have 𝔇(f0)=g0\mathfrak{D}(f_{0})=g_{0} and hence f0=u|V0f_{0}=u|_{V_{0}} from uniqueness. On the other hand, over WW,

𝔇(f1)=g0+gϵ=g0g0=0.\mathfrak{D}(f_{1})=g_{0}+g_{\epsilon}=g_{0}-g_{0}=0.

Thus f1f_{1} is a solution 𝔇(f1)=0\mathfrak{D}(f_{1})=0, over WW. Furthermore, ftf_{t} is admissible and

j𝔇(ft)α(σ)=jg0+tgϵα(σ)=0, as gϵ=g0 on some 𝒪p(σ).j^{\alpha}_{\mathfrak{D}(f_{t})}(\sigma)=j^{\alpha}_{g_{0}+tg_{\epsilon}}(\sigma)=0,\text{ as $g_{\epsilon}=-g_{0}$ on some $\operatorname{\mathcal{O}p}(\sigma)$.}

Thus, jftα+1(σ)αj^{\alpha+1}_{f_{t}}(\sigma)\in\mathcal{R}_{\alpha} for all t[0,1]t\in[0,1]. ∎

We now have a (parametric) local hh-principle for α\mathcal{R}_{\alpha} ([Gro86, pg. 119])

Corollary 4.3.

For α\alpha large enough, the jet map jα+1:SolαΓαj^{\alpha+1}:\operatorname{Sol}\mathcal{R}_{\alpha}\to\Gamma\mathcal{R}_{\alpha} is a local weak homotopy equivalence.

In order to prove the existence of a horizontal germ, i.e, a local solution of 𝔇=0\mathfrak{D}=0, we need to show that α\mathcal{R}_{\alpha}\neq\emptyset at some σ\sigma. One issue with Theorem 4.2 is that we do not specify the higher jet order α\alpha that is crucial in order to get a local solution. We now show that, in fact, we can get a lift to any arbitrary higher jet from the first jet relation of isotropic horizontal maps. Recall that given any map ff satisfying fλi=0f^{*}\lambda^{i}=0 we have, taking derivatives, that fdλi=0f^{*}d\lambda^{i}=0. That is, Imdf\operatorname{Im}df is dλid\lambda^{i}-isotropic. Now from 3.3, we have that every solution is automatically admissible. On the other hand, we have the relation 0J1(Σ,M)\mathcal{R}\subset\mathcal{R}_{0}\subset J^{1}(\Sigma,M) consisting of jets (x,y,F:TxΣTyM)(x,y,F:T_{x}\Sigma\to T_{y}M) such that, Fdαs=0F^{*}d\alpha^{s}=0 for s=1,2s=1,2. That is, sections of \mathcal{R} are bundle maps F:TΣTMF:T\Sigma\to TM, which is a formal isotropic 𝒟\mathcal{D}-horizontal immersion. Observe that,

Sol=Solα,for any α0.\operatorname{Sol}\mathcal{R}=\operatorname{Sol}\mathcal{R}_{\alpha},\;\text{for any $\alpha\geq 0$.}

We have the following result.

Lemma 4.4.

For any α1\alpha\geq 1, the jet projection map p=p1α+1:Jα+1(Σ,M)J1(Σ,M)p=p^{\alpha+1}_{1}:J^{\alpha+1}(\Sigma,M)\to J^{1}(\Sigma,M) maps α|(x,y)\mathcal{R}_{\alpha}|_{(x,y)} surjectively onto |(x,y)\mathcal{R}|_{(x,y)}, for any (x,y)Σ×M(x,y)\in\Sigma\times M. Furthermore, the fiber of pp over any jet in |(x,y)\mathcal{R}|_{(x,y)} is contractible and consequently, the induced map ΓαΓ\Gamma\mathcal{R}_{\alpha}\to\Gamma\mathcal{R} is a local weak homotopy equivalence.

We defer the details of the proof to subsection 4.1. Let us first get the local hh-principle.

Proof of Theorem 1.3.

From 4.3, we have that for α\alpha sufficiently large, the jet map jα+1:SolΓαj^{\alpha+1}:\operatorname{Sol}\mathcal{R}\to\Gamma\mathcal{R}_{\alpha} is a local weak weak homotopy equivalence. On the other hand, by 4.4, the jet projection p1α+1:ΓαΓp^{\alpha+1}_{1}:\Gamma\mathcal{R}_{\alpha}\to\Gamma\mathcal{R} is a weak homotopy equivalence. Hence the composition

j1=p1α+1jα+1:SolΓj^{1}=p^{\alpha+1}_{1}\circ j^{\alpha+1}:\operatorname{Sol}\mathcal{R}\to\Gamma\mathcal{R}

is a local weak homotopy equivalence. In other words, 𝒟\mathcal{D}-horizontal immersions satisfy the (parametric) local hh-principle. ∎

Next, we get the existence of germs of 𝒟\mathcal{D}-horizontal 22-submanifolds.

Proof of Theorem 1.4.

Suppose 𝒟=kerλ1kerλ2\mathcal{D}=\ker\lambda^{1}\cap\ker\lambda^{2} for some local 11-forms λi\lambda^{i} around some yMy\in M. Pick some arbitrary 0v𝒟x0\neq v\in\mathcal{D}_{x} and set u=Avu=Av, where AA is the (local) automorphism. Then observe that,

dλ1(u,v)=dλ1(Av,v)=dλ2(v,v)=0and,dλ2(u,v)=dλ1(u,Av)=dλ1(u,u)=0.d\lambda^{1}(u,v)=d\lambda^{1}(Av,v)=d\lambda^{2}(v,v)=0\quad\text{and,}\quad d\lambda^{2}(u,v)=d\lambda^{1}(u,Av)=d\lambda^{1}(u,u)=0.

In other words, Spanu,v𝒟x\textrm{Span}\langle u,v\rangle\subset\mathcal{D}_{x} is Ω\Omega-isotropic. Now, consider the jet σ=(0,y,F:T0𝔻2TyM)J1(𝔻2,M)\sigma=(0,y,F:T_{0}\mathbb{D}^{2}\to T_{y}M)\in J^{1}(\mathbb{D}^{2},M), given by,

F(x)=u,F(y)=v.F(\partial_{x})=u,\;F(\partial_{y})=v.

We clearly have σ|(0,y)\sigma\in\mathcal{R}|_{(0,y)} by construction. But then an application of the the local hh-principle gives us that there exists a 𝒟\mathcal{D}-horizontal immersion f:𝒪p(0)Mf:\operatorname{\mathcal{O}p}(0)\to M. Since ff is an immersion, it is a local diffeomorphism and thus we have a (germ of a) 𝒟\mathcal{D}-horizontal submanifold of dimension 22. ∎

4.1. Proof of 4.4

In this section, we discuss the proof of 4.4. Instead of proving it only for fat distributions of rank 44 on 66-dimensional manifold, we consider an arbitrary fat distribution of corank pp on manifolds of dimension NN. Since we are only considering jets of maps, let us consider Σ=2\Sigma=\mathbb{R}^{2} and M=NM=\mathbb{R}^{N}, with fixed coordinates {x1,x2}\{x^{1},x^{2}\} on Σ\Sigma and {y1,,yN}\{y^{1},\ldots,y^{N}\} on MM. Suppose 𝒟TM\mathcal{D}\subset TM is a corank pp fat distribution, given as the common kernel of 11-forms λ1,,λp\lambda^{1},\ldots,\lambda^{p} where we have λs=λisdyi\lambda^{s}=\lambda^{s}_{i}dy^{i}. For any f:ΣMf:\Sigma\to M we have the operator

𝔇:f(fλs)s=1,,p\mathfrak{D}:f\mapsto\big{(}f^{*}\lambda^{s}\big{)}_{s=1,\ldots,p}

We need to understand the relation αJα+1(Σ,M)\mathcal{R}_{\alpha}\subset J^{\alpha+1}(\Sigma,M), which consists of jets jfα+1(x)j^{\alpha+1}_{f}(x), where f:𝒪p(x)Mf:\operatorname{\mathcal{O}p}(x)\to M satisfies j𝔇(f)α(x)=0j^{\alpha}_{\mathfrak{D}(f)}(x)=0. We also have the first jet relation J1(Σ,M)\mathcal{R}\subset J^{1}(\Sigma,M) consisting of jets (x,y,F:TxΣTyM)(x,y,F:T_{x}\Sigma\to T_{y}M), where FF is an injective map taking TσΣT_{\sigma}\Sigma to a common isotropic subspace of 𝒟y\mathcal{D}_{y}. We prove the following stronger version of 4.4.

Lemma 4.5.

The jet projection map p=p1α+1:Jα+1(Σ,M)J1(Σ,M)p=p^{\alpha+1}_{1}:J^{\alpha+1}(\Sigma,M)\to J^{1}(\Sigma,M) maps α|(x,y)\mathcal{R}_{\alpha}|_{(x,y)} surjectively onto |(x,y)\mathcal{R}|_{(x,y)}. Furthermore, the fiber of pp over any fixed jet is contractible and consequently p:ΓαΓp:\Gamma\mathcal{R}_{\alpha}\to\Gamma\mathcal{R} is a weak homotopy equivalence.

First we need to understand the equation j𝔇(f)α(x)j^{\alpha}_{\mathfrak{D}(f)}(x) in terms of jets. We write down, fλs=ηasdxaf^{*}\lambda^{s}=\eta^{s}_{a}dx^{a}, where we have the functions

ηas=fλs(a)=(λisf)afi,1sp,1a2.\eta^{s}_{a}=f^{*}\lambda^{s}(\partial_{a})=(\lambda^{s}_{i}\circ f)\partial_{a}f^{i},\qquad 1\leq s\leq p,1\leq a\leq 2.

Since 𝔇(f)=(fλ1,,fλp)\mathfrak{D}(f)=\big{(}f^{*}\lambda^{1},\ldots,f^{*}\lambda^{p}\big{)}, we have that the jet j𝔇(f)α(x)j^{\alpha}_{\mathfrak{D}(f)}(x) is nothing but jηasα(x)j^{\alpha}_{\eta^{s}_{a}}(x). We need to find out its higher derivatives. Let us fix our conventions first.

Convention for Multi-Indices :

By a multi-index of order rr on the coordinates {x1,x2}\{x^{1},x^{2}\} we will consider an ordered tuple I=(i1ir)I=(i_{1}\leq\ldots\leq i_{r}) where ij{1,2}i_{j}\in\{1,2\}. We denote,

I=i1ir=xi1xir.\partial_{I}=\partial_{i_{1}}\ldots\partial_{i_{r}}=\partial_{x^{i_{1}}}\ldots\partial_{x^{i_{r}}}.

A typical multi-index of order α\alpha looks like (1,,1a,2,,2b)\big{(}\underbrace{1,\ldots,1}_{a},\underbrace{2,\ldots,2}_{b}\big{)}, for some integers a,b0a,b\geq 0 satisfying α=a+b\alpha=a+b. We denote |I||I| as the order of a multi-index. A multi-index of order 11 will be written without the parentheses. Let us denote by 2r\mathbb{N}_{2}^{r} the set of all multi-indices of order rr over the coordinates {x1,x2}\{x^{1},x^{2}\}.

For any subsequence III^{\prime}\subset I, we will denote III-I^{\prime} as the multi-index obtained by taking the complimentary sequence. In particular for a given I=(i1,,ir)I=(i_{1},\ldots,i_{r}), we have

Ij:=Iij=(i1,,ij1,ij+1,,ir).I_{j}:=I-i_{j}=\big{(}i_{1},\ldots,i_{j-1},i_{j+1},\ldots,i_{r}\big{)}.

Also, for given two multi-indices II and JJ, by I+JI+J we will mean the multi-index obtained by first concatenating the two lists and then reordering in increasing order. In particular, observe that for any iIi\in I we have, (Ii)+i=I(I-i)+i=I.

Now, for any multi-index I=(i1,,iα)I=(i_{1},\ldots,i_{\alpha}) of order α\alpha, we have the higher order partial derivative,

Iηas\displaystyle\partial_{I}\eta^{s}_{a} =(λisf)I+afi\displaystyle=(\lambda^{s}_{i}\circ f)\partial_{I+a}f^{i}
+(νλisf)Ifνafi+l=1α(νλisf)ilfνIil+afi\displaystyle\qquad+(\partial_{\nu}\lambda^{s}_{i}\circ f)\partial_{I}f^{\nu}\partial_{a}f^{i}+\sum_{l=1}^{\alpha}(\partial_{\nu}\lambda^{s}_{i}\circ f)\partial_{i_{l}}f^{\nu}\partial_{I-i_{l}+a}f^{i}
+terms involving partial derivatives of f of order <α\displaystyle\qquad+\text{terms involving partial derivatives of $f$ of order $<\alpha$}

Then, the condition j𝔇(f)α(x)=0j^{\alpha}_{\mathfrak{D}(f)}(x)=0 is equivalent to the set of equations

{Iηas(x)=0| 1sp, 1a2,|I|α}.\big{\{}\partial_{I}\eta^{s}_{a}(x)=0\;\big{|}\;1\leq s\leq p,\;1\leq a\leq 2,\;|I|\leq\alpha\}.

We now translate this system in the language of jets.

Consider an α+1\alpha+1-jet

σ=(x,y,Pi:SymiTxΣTyM,i=1,,α+1)J(x,y)α+1(Σ,M).\sigma=\big{(}x,\;y,\;P_{i}:\operatorname{Sym}^{i}T_{x}\Sigma\to T_{y}M,\;i=1,\ldots,\alpha+1\big{)}\in J^{\alpha+1}_{(x,y)}(\Sigma,M).

Note that the ithi^{\text{th}} symmetric tensor PiP_{i} is completely determined by the values Pik(I)P^{k}_{i}(\partial_{I}), where PikP^{k}_{i} are the components of PiP_{i} and {I|I2i}\{\partial_{I}\;|\;I\in\mathbb{N}_{2}^{i}\} forms a basis of the symmetric space SymiTxΣ\operatorname{Sym}^{i}T_{x}\Sigma. Now, suppose the jet σ\sigma is represented as jfα+1(x)j^{\alpha+1}_{f}(x) for some f:𝒪p(x)Mf:\operatorname{\mathcal{O}p}(x)\to M. Then for some I=(i1,,ir)I=(i_{1},\ldots,i_{r}) where rαr\leq\alpha and some 1a2, 1sp1\leq a\leq 2,\;1\leq s\leq p, the equation Iηas(x)=0\partial_{I}\eta^{s}_{a}(x)=0 gives us the following.

Iηas\boxed{\partial_{I}\eta^{s}_{a}} λks(y)Pr+1k(I+a)+νλks(y)P1k(a)Prν(I)+l=1rνλks(y)P1ν(il)Prk(Iil+a)+terms involving P<r=0\displaystyle\begin{aligned} \lambda^{s}_{k}(y)P_{r+1}^{k}(\partial_{I+a})+\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{a})P_{r}^{\nu}(\partial_{I})+\sum_{l=1}^{r}\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{\nu}(\partial_{i_{l}})P_{r}^{k}(\partial_{I-i_{l}+a})+\text{terms involving $P_{<r}$}=0\end{aligned}

Note that we have labeled these equations by (Iηas)(\partial_{I}\eta^{s}_{a}) as well. In particular, we have that the relation α\mathcal{R}_{\alpha} consists of jets σJα+1(Σ,M)\sigma\in J^{\alpha+1}(\Sigma,M) which satisfy the jet equations (Iηas)(\partial_{I}\eta^{s}_{a}) for each 1sp,1a21\leq s\leq p,1\leq a\leq 2 and for each multi-index |I|α|I|\leq\alpha.

Now, a jet (x,y,P1:TxΣTyM)J(x,y)1(Σ,M)(x,y,P_{1}:T_{x}\Sigma\to T_{y}M)\in J^{1}_{(x,y)}(\Sigma,M) belongs to 0|(x,y)\mathcal{R}_{0}|_{(x,y)} precisely when the equations

λks(y)P1k(a)=0,1sp,1a2\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{a})=0,\quad 1\leq s\leq p,1\leq a\leq 2

hold. This is nothing but the condition that P1(TxΣ)𝒟yP_{1}(T_{x}\Sigma)\subset\mathcal{D}_{y}. We also have the relation 0J1(Σ,M)\mathcal{R}\subset\mathcal{R}_{0}\subset J^{1}(\Sigma,M) consisting of jets σ=(x,y,P1)\sigma=(x,y,P_{1}) which further satisfy the isotropy condition. That is for each v,wTxΣv,w\in T_{x}\Sigma we have, dλs(P1(v),P1(w))=01spd\lambda^{s}(P_{1}(v),P_{1}(w))=0\quad 1\leq s\leq p. This isotropy condition can be expressed more explicitly as,

(νλμs(y)μλνs(y))P1ν(a)P1μ(b)=0,for each 1a,b2 and each 1sp.\Big{(}\partial_{\nu}\lambda^{s}_{\mu}(y)-\partial_{\mu}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{\nu}(\partial_{a})P_{1}^{\mu}(\partial_{b})=0,\quad\text{for each $1\leq a,b\leq 2$ and each $1\leq s\leq p$.}

For αβ\alpha\geq\beta, we have the jet projection map p=pβ+1α+1:Jα+1(Σ,M)Jβ+1(Σ,M)p=p^{\alpha+1}_{\beta+1}:J^{\alpha+1}(\Sigma,M)\to J^{\beta+1}(\Sigma,M) and clearly, pp maps α\mathcal{R}_{\alpha} into β\mathcal{R}_{\beta}, since every defining equation for β\mathcal{R}_{\beta} is also used to define α\mathcal{R}_{\alpha}. We now proceed with the proof of 4.5, which is essentially done via induction on α1\alpha\geq 1.

Induction Statement 𝒫(α)\mathcal{P}(\alpha) :

For given α1\alpha\geq 1, the map p=p1α+1p=p^{\alpha+1}_{1} maps α\mathcal{R}_{\alpha} into \mathcal{R}. The system of equations involved in defining the relation α\mathcal{R}_{\alpha} can be replaced by another system such that the following holds.

  • The solution space remains unchanged.

  • The new system is a triangular one. That is, we can solve P1,P2,P_{1},P_{2},\ldots etc successively.

  • The highest order symmetric tensor Pα+1P_{\alpha+1} in the system appears as follows :

    λα+1J\boxed{\lambda_{\alpha+1}^{J}} (λ11(y)λN1(y)λ1p(y)λNp(y))(Pα+11(J)Pα+1N(J))=p×1-vector involving Pα terms\begin{pmatrix}\lambda^{1}_{1}(y)&\ldots&\lambda^{1}_{N}(y)\\ \vdots&&\vdots\\ \lambda^{p}_{1}(y)&\ldots&\lambda^{p}_{N}(y)\end{pmatrix}\begin{pmatrix}P_{\alpha+1}^{1}(\partial_{J})\\ \vdots\\ P_{\alpha+1}^{N}(\partial_{J})\end{pmatrix}=\text{a $p\times 1$-vector involving $P_{\leq\alpha}$ terms}

    for each multi-index |J|=α+1|J|=\alpha+1. Note that we label these systems by (λαJ)(\lambda_{\alpha}^{J}), since the coefficient matrix is the matrix of the p\mathbb{R}^{p}-valued 11-form λ=(λ1,,λp)\lambda=(\lambda^{1},\ldots,\lambda^{p}).

  • The system is consistent and admits solutions. Furthermore, the solution space is contractible.

The induction relies on two facts : We will see that 𝒫(1)\mathcal{P}(1) holds because by assumption the 11-forms λ1,,λp\lambda^{1},\ldots,\lambda^{p} are independent at each point and this independence enables us to solve certain affine system. Whereas, for any α1\alpha\geq 1, to get 𝒫(α+1)\mathcal{P}(\alpha+1) from 𝒫(α)\mathcal{P}(\alpha), we need to utilize the fact the distribution 𝒟\mathcal{D} is in fact fat. Let us now proceed with the details.

Induction Base Step 𝒫(1)\mathcal{P}(1) :

We focus on the relation 1J2(Σ,M)\mathcal{R}_{1}\subset J^{2}(\Sigma,M). Consider some jet,

σ~=(x,y,P1:TxΣTyM,P2:Sym2TxΣTyM)|(x,y)\tilde{\sigma}=(x,y,P_{1}:T_{x}\Sigma\to T_{y}M,P_{2}:\operatorname{Sym}^{2}T_{x}\Sigma\to T_{y}M)\in\mathcal{R}|_{(x,y)}

For 1a<b21\leq a<b\leq 2, i.e, for a=1,b=2a=1,b=2, we have the equations defining 1\mathcal{R}_{1} as follows :

aηbs\boxed{\partial_{a}\eta^{s}_{b}} λks(y)P2k(ab)+νλks(y)P1k(b)P1ν(a)=0\displaystyle\lambda^{s}_{k}(y)P_{2}^{k}(\partial_{ab})+\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{b})P_{1}^{\nu}(\partial_{a})=0
bηas\boxed{\partial_{b}\eta^{s}_{a}} λks(y)P2k(ab)+νλks(y)P1k(a)P1ν(b)=0\displaystyle\lambda^{s}_{k}(y)P_{2}^{k}(\partial_{ab})+\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{a})P_{1}^{\nu}(\partial_{b})=0

Since σ~1|(x,y)\tilde{\sigma}\in\mathcal{R}_{1}|_{(x,y)} must satisfy both (aηbs)(\partial_{a}\eta_{b}^{s}) and (bηas)(\partial_{b}\eta_{a}^{s}) above, we perform (aηbs)(bηas)(\partial_{a}\eta^{s}_{b})-(\partial_{b}\eta^{s}_{a}) and get

0=(νλks(y)kλνs(y))P1ν(a)P1k(b).\displaystyle 0=\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{\nu}(\partial_{a})P_{1}^{k}(\partial_{b}).

But observe that the isotropy condition on P1P_{1} is precisely

0=dλs|y(P1(a),P1(b))=(νλks(y)kλνs(y))P1ν(a)P1k(b).0=d\lambda^{s}|_{y}\big{(}P_{1}(\partial_{a}),P_{1}(\partial_{b})\big{)}=\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{\nu}(\partial_{a})P_{1}^{k}(\partial_{b}).

Hence we see that p(σ~)=(x,y,P1)|(x,y)p(\tilde{\sigma})=(x,y,P_{1})\in\mathcal{R}|_{(x,y)}, as P1P_{1} satisfies the isotropy condition. Thus, pp maps 1\mathcal{R}_{1} into \mathcal{R}.

Next, assume that σ=(x,y,P1:TxΣTyM)|(x,y)\sigma=(x,y,P_{1}:T_{x}\Sigma\to T_{y}M)\in\mathcal{R}|_{(x,y)} is given. We need to look for a jet σ~=(x,y,P1,P2)1|(x,y)\tilde{\sigma}=(x,y,P_{1},P_{2})\in\mathcal{R}_{1}|_{(x,y)}. That is, we need to find out P2P_{2} satisfying the equations {aηbs|1sp,1ab2}\{\partial_{a}\eta^{s}_{b}\big{|}1\leq s\leq p,1\leq a\leq b\leq 2\}. Now for a<ba<b, we see that the terms P2k(ab)P_{2}^{k}(\partial_{ab}) appear linearly in both the systems {aηbs|1sp}\{\partial_{a}\eta^{s}_{b}|1\leq s\leq p\} and {bηas|1sp}\{\partial_{b}\eta^{s}_{a}|1\leq s\leq p\}, with identical coefficients. Thus we run into a question of consistency. But we have seen that the expression (aηbs)(bηas)(\partial_{a}\eta^{s}_{b})-(\partial_{b}\eta^{s}_{a}) is identical to the isotropy condition for P1P_{1}. Since the jet σ\sigma\in\mathcal{R}, we know that P1P_{1} satisfies the isotropy condition. Thus we have that the equations (aηbs)(bηas)(\partial_{a}\eta^{s}_{b})-(\partial_{b}\eta^{s}_{a}) is zero and so for each tuple aba\leq b, we can keep the equation labeled by aηbs\partial_{a}\eta^{s}_{b} and remove the equation bηbs\partial_{b}\eta_{b}^{s}. We are left with the system

λ2ab\boxed{\lambda_{2}^{ab}} (λ11(y)λN1(y)λ1p(y)λNp(y))(P21(ab)P2N(ab))=p×1 vector involving P1\begin{pmatrix}\lambda^{1}_{1}(y)&\ldots&\lambda^{1}_{N}(y)\\ \vdots&&\vdots\\ \lambda^{p}_{1}(y)&\ldots&\lambda^{p}_{N}(y)\end{pmatrix}\begin{pmatrix}P_{2}^{1}(\partial_{ab})\\ \vdots\\ P_{2}^{N}(\partial_{ab})\end{pmatrix}=\text{a $p\times 1$ vector involving $P_{1}$}

Clearly the solution space remains unchanged.

Lastly, to show that P2P_{2} can be solved, note that the above affine system has full rank coefficient matrix, since the rows are nothing but the 11-forms λs\lambda^{s}, which are given to be independent. Thus, the system admits a solution. We have proved that p:1|(x,y)|(x,y)p:\mathcal{R}_{1}|_{(x,y)}\to\mathcal{R}|_{(x,y)} is surjective with affine fibers.

Induction Hypothesis 𝒫(α)\mathcal{P}(\alpha) :

Assume that for some α1\alpha\geq 1, the statement 𝒫(α)\mathcal{P}(\alpha) holds.

Induction Step 𝒫(α)𝒫(α+1)\mathcal{P}(\alpha)\Rightarrow\mathcal{P}(\alpha+1) :

First observe that the jet map pp maps α+1\mathcal{R}_{\alpha+1} into \mathcal{R}; since we have already proved this for 1\mathcal{R}_{1} and the equations involved in 1\mathcal{R}_{1} are also present in α+1\mathcal{R}_{\alpha+1}. We prove that pp is surjective, with contractible fiber.

Fix a jet σ=(x,y,P1)|(x,y)\sigma=(x,y,P_{1})\in\mathcal{R}|_{(x,y)}. Since the equation system defining α\mathcal{R}_{\alpha} is included in α+1\mathcal{R}_{\alpha+1}, using the induction hypothesis 𝒫(α)\mathcal{P}(\alpha), we replace this (sub)system with the triangular system, keeping the equations involving Pα+2P_{\alpha+2} untouched. Next, we solve for the tensors PrP_{r} with rαr\leq\alpha from this transformed system. Note that, we could also solve for the tensors Pα+1P_{\alpha+1} using the induction hypothesis, but we defer this for later. As we will see that while solving for the tensor Pα+2P_{\alpha+2} we will run into some consistency problem, which will introduce new sets of equations for Pα+1P_{\alpha+1}. The following ladder like diagram gives a schematic representation of this step :

Jα+2(Σ,M){J^{\alpha+2}(\Sigma,M)}Jα+1(Σ,M){J^{\alpha+1}(\Sigma,M)}Jα(Σ,M){J^{\alpha}(\Sigma,M)}{\cdots}J2(Σ,M){J^{2}(\Sigma,M)}J1(Σ,M){J^{1}(\Sigma,M)}α+1{\mathcal{R}_{\alpha+1}}α{\mathcal{R}_{\alpha}}α1{\mathcal{R}_{\alpha-1}}{\cdots}1{\mathcal{R}_{1}}0{\mathcal{R}_{0}}𝒮α{\mathcal{S}_{\alpha}}Sα1{S_{\alpha-1}}S0={S_{0}=\mathcal{R}}pα+1α+2\scriptstyle{p^{\alpha+2}_{\alpha+1}}pαα+1\scriptstyle{p^{\alpha+1}_{\alpha}}p12\scriptstyle{p^{2}_{1}}pα+1α+2\scriptstyle{p^{\alpha+2}_{\alpha+1}}\subset\subset\subset\subset\subsetadditionalaffine eqnslift usingindependence of formslift usingfatness of 𝒟\scriptstyle\mathcal{D}isotropy(quadratic)lift to α\scriptstyle\mathcal{R}_{\alpha} using induction hypothesisand then project to 𝒮α1\scriptstyle\mathcal{S}_{\alpha-1}

where, we have denoted Sβ=pβ+1β+2(β+1)S_{\beta}=p^{\beta+2}_{\beta+1}(\mathcal{R}_{\beta+1}) as the image. So, using the induction hypothesis, we first get a lift of the jet σ\sigma to 𝒮α1=pαα+1(α)\mathcal{S}_{\alpha-1}=p^{\alpha+1}_{\alpha}(\mathcal{R}_{\alpha}). We now need to identify 𝒮α=pα+1α+2(α+1)\mathcal{S}_{\alpha}=p^{\alpha+2}_{\alpha+1}(\mathcal{R}_{\alpha+1}), which is defined by the new set of equations coming from the consistency.

Let us first fix some more notations for the multi-indices. Recall that 2r\mathbb{N}_{2}^{r} is the set of all multi-indices II on {x1,x2}\{x^{1},x^{2}\}, with order |I|=r|I|=r. Now for some J2r+1J\in\mathbb{N}_{2}^{r+1} we denote

(J)={(I,a)|I2r, 1a2,J=I+a}\mathfrak{I}(J)=\Big{\{}(I,a)\;\Big{|}\;I\in\mathbb{N}_{2}^{r},\;1\leq a\leq 2,J=I+a\Big{\}}

Then observe that for any (I,a)(J)(I,a)\in\mathfrak{I}(J) and some 1sp1\leq s\leq p, the equation labeled by Iηas\partial_{I}\eta^{s}_{a} involves the terms Pr+1k(J)P_{r+1}^{k}(\partial_{J}). Explicitly, if we consider J=(a0,a1,,ar)J=(a_{0},a_{1},\ldots,a_{r}) with 1a0ar21\leq a_{0}\leq\ldots\leq a_{r}\leq 2, then we get

(J)=i=0r{(Ii,ai)|Ii:=Iai=(a0,,ai^,,ar)}\mathfrak{I}(J)=\bigcup_{i=0}^{r}\Big{\{}(I_{i},a_{i})\;\Big{|}\;I_{i}:=I-a_{i}=\big{(}a_{0},\ldots,\hat{a_{i}},\ldots,a_{r}\big{)}\Big{\}}

Since we may have repetitions in the aia_{i}’s, we see that |(J)||J|=r+1|\mathfrak{I}(J)|\leq|J|=r+1. In fact, since we have only two indices {x1,x2}\{x^{1},x^{2}\}, we can see that |(J)|2|\mathfrak{I}(J)|\leq 2.

We make the following observation. Fix some J=(a0,,aα+1)2α+2J=(a_{0},\ldots,a_{\alpha+1})\in\mathbb{N}_{2}^{\alpha+2} with 1a0aα+121\leq a_{0}\leq\ldots\leq a_{\alpha+1}\leq 2. Consider (I0,a0)(I_{0},a_{0}) and (Ij,aj)(I_{j},a_{j}) in (J)\mathfrak{I}(J) for some 1jr+11\leq j\leq r+1 fixed. Then for each 1sp1\leq s\leq p, we have the equations I0ηa0s\partial_{I_{0}}\eta^{s}_{a_{0}} and Ijηajs\partial_{I_{j}}\eta^{s}_{a_{j}} as follows :

I0ηa0s\boxed{\partial_{I_{0}}\eta^{s}_{a_{0}}} λks(y)Pα+2k(J)+νλks(y)P1k(a0)Pα+1ν(I0)+i=1α+1νλks(y)P1ν(ai)Pα+1k(Jai)+terms involving Pα=0\displaystyle\begin{aligned} &\lambda^{s}_{k}(y)P_{\alpha+2}^{k}(\partial_{J})+\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{a_{0}})P_{\alpha+1}^{\nu}(\partial_{I_{0}})+\sum_{i=1}^{\alpha+1}\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{\nu}(\partial_{a_{i}})P_{\alpha+1}^{k}(\partial_{J-a_{i}})\\ &\hskip 142.26378pt+\text{terms involving $P_{\leq\alpha}$}=0\end{aligned}
Ijηajs\boxed{\partial_{I_{j}}\eta^{s}_{a_{j}}} λks(y)Pα+2k(J)+νλks(y)P1k(aj)Pα+1ν(Ij)+i=0ijα+1νλks(y)P1ν(ai)Pα+1k(Jai)+terms involving Pα=0\displaystyle\begin{aligned} &\lambda^{s}_{k}(y)P_{\alpha+2}^{k}(\partial_{J})+\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{k}(\partial_{a_{j}})P_{\alpha+1}^{\nu}(\partial_{I_{j}})+\sum_{\begin{subarray}{c}i=0\\ i\neq j\end{subarray}}^{\alpha+1}\partial_{\nu}\lambda^{s}_{k}(y)P_{1}^{\nu}(\partial_{a_{i}})P_{\alpha+1}^{k}(\partial_{J-a_{i}})\\ &\hskip 142.26378pt+\text{terms involving $P_{\leq\alpha}$}=0\end{aligned}

If a0=aja_{0}=a_{j}, then I0=IjI_{0}=I_{j} and hence the equations I0ηa0s\partial_{I_{0}}\eta^{s}_{a_{0}} and Ijηajs\partial_{I_{j}}\eta^{s}_{a_{j}} becomes identical. So without a loss of generality, we assume that a0<aja_{0}<a_{j}. If no such jj exists, then we must have J=(a0,,a0α+2-times)J=(\underbrace{a_{0},\ldots,a_{0}}_{\text{$\alpha+2$-times}}) which gives a single equation I0ηa0s\partial_{I_{0}}\eta^{s}_{a_{0}} to consider. Note that the only way a0<aja_{0}<a_{j} can hold is for some jj is a0=1,aj=2a_{0}=1,a_{j}=2.

Now, observe that the two systems {I0ηa0s,1sp}\{\partial_{I_{0}}\eta^{s}_{a_{0}},1\leq s\leq p\} and {Ijηajs,1sp}\{\partial_{I_{j}}\eta^{s}_{a_{j}},1\leq s\leq p\} both look like

λα+2J\boxed{\lambda_{\alpha+2}^{J}} (λ11(y)λN1(y)λ1p(y)λNp(y))(Pα+21(J)Pα+2N(J))=p×1 vector involving Pα+1\begin{pmatrix}\lambda^{1}_{1}(y)&\ldots&\lambda^{1}_{N}(y)\\ \vdots&&\vdots\\ \lambda^{p}_{1}(y)&\ldots&\lambda^{p}_{N}(y)\end{pmatrix}\begin{pmatrix}P_{\alpha+2}^{1}(\partial_{J})\\ \vdots\\ P_{\alpha+2}^{N}(\partial_{J})\end{pmatrix}=\text{a $p\times 1$ vector involving $P_{\leq\alpha+1}$}

And thus we run into the question of consistency from these two systems. To address this issue, we first perform (I0ηa0s)(Ijηajs)(\partial_{I_{0}}\eta^{s}_{a_{0}})-(\partial_{I_{j}}\eta^{s}_{a_{j}}) and get

(I0ηa0s)(Ijηajs)\boxed{(\partial_{I_{0}}\eta^{s}_{a_{0}})-(\partial_{I_{j}}\eta^{s}_{a_{j}})} (νλks(y)kλνs(y))(P1k(a0)Pα+1ν(I0)P1k(aj)Pα+1ν(Ij))\displaystyle\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}\Big{(}P_{1}^{k}(\partial_{a_{0}})P_{\alpha+1}^{\nu}(\partial_{I_{0}})-P_{1}^{k}(\partial_{a_{j}})P_{\alpha+1}^{\nu}(\partial_{I_{j}})\Big{)}
+terms involving Pα=0\displaystyle\hskip 113.81102pt+\text{terms involving $P_{\leq\alpha}$}=0

Note that the difference does not involve any Pα+2P_{\alpha+2} terms at all. We keep the equations labeled by (I0ηa0s\partial_{I_{0}}\eta^{s}_{a_{0}}) and for each j=1,,α+1j=1,\ldots,\alpha+1, we replace the equation (Ijηajs\partial_{I_{j}}\eta^{s}_{a_{j}}) by the equation ((I0ηa0s)(Ijηajs)(\partial_{I_{0}}\eta^{s}_{a_{0}})-(\partial_{I_{j}}\eta^{s}_{a_{j}})); provided a0<aja_{0}<a_{j}. Clearly this does not change the solution space, but introduces new set of equations involving Pα+1P_{\leq\alpha+1}. The system is still affine. These added set of equations, together with the original system involving Pα+1P_{\alpha+1}, now define 𝒮α=pα+1α+2(α+1)\mathcal{S}_{\alpha}=p^{\alpha+2}_{\alpha+1}(\mathcal{R}_{\alpha+1}) (see diagram (4.1)). Note that, at this point we are left with exactly one system that involves Pα+2k(J)P_{\alpha+2}^{k}(\partial_{J}), which looks like (λα+2J\lambda^{J}_{\alpha+2}) as above.

Let us now fix the dictionary order \prec on 2α+1\mathbb{N}_{2}^{\alpha+1}, induced by the obvious ordering of the coordinate indices. Then we have IjI0I_{j}\prec I_{0}, since the first position they differ must be larger than a0a_{0}. We rewrite the equation ((I0ηa0s)(Ijηajs)(\partial_{I_{0}}\eta^{s}_{a_{0}})-(\partial_{I_{j}}\eta^{s}_{a_{j}})) as

(I0ηa0s)(Ijηajs)\boxed{(\partial_{I_{0}}\eta^{s}_{a_{0}})-(\partial_{I_{j}}\eta^{s}_{a_{j}})} (νλks(y)kλνs(y))P1k(a0)Pα+1ν(I0)\displaystyle\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{k}(\partial_{a_{0}})P_{\alpha+1}^{\nu}(\partial_{I_{0}}) =(νλks(y)kλνs(y))P1k(aj)Pα+1ν(Ij)\displaystyle=\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{k}(\partial_{a_{j}})P_{\alpha+1}^{\nu}(\partial_{I_{j}})
+terms with Pα\displaystyle\qquad+\text{terms with $P_{\leq\alpha}$}

and add these to the set of equations that is used to solve for Pα+1k(I0)P^{k}_{\alpha+1}(I_{0}), which, from from the induction hypothesis, is given as the system (λα+1I0\lambda^{I_{0}}_{\alpha+1}).

Observe that we are adding at most one new system of pp-many equations for each tensor Pα+1k(I)P_{\alpha+1}^{k}(\partial_{I}), since we have |(I)|2|\mathfrak{I}(I)|\leq 2. In fact, it is clear that for any fixed |I|=α+1|I|=\alpha+1, the system that we are adding looks like

(νλks(y)kλνs(y))P1k(1)Pα+1ν(I)=known terms with Pα and Pα+1(I) with II,for s=1,,p.\displaystyle\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{k}(\partial_{1})P_{\alpha+1}^{\nu}(\partial_{I})=\text{known terms with $P_{\leq\alpha}$ and $P_{\alpha+1}(\partial_{I^{\prime}})$ with $I^{\prime}\prec I$,}\quad\text{for $s=1,\ldots,p$.}

Furthermore, we are not adding any equation when I=I^:=(1,,1α+1)I=\hat{I}:=(\underbrace{1,\ldots,1}_{\alpha+1}), which is the least element in the ordering \prec.

We claim that we are able to solve the Pα+1P_{\alpha+1}-terms in a triangular fashion, ordered by \prec. To see this, first observe that the only equation system that involves Pα+1(I^)P_{\alpha+1}(\partial_{\hat{I}}) looks like (λα+1I^\lambda^{\hat{I}}_{\alpha+1}). But this affine system admits solutions, since the coefficient matrix has full rank. Now, inductively assume that for some multi-index II, with |I|=α+1|I|=\alpha+1 and I^I\hat{I}\prec I, we have solved the tensor Pα+1(I)P_{\alpha+1}(\partial_{I^{\prime}}) for any III^{\prime}\prec I. For this II, we have added the following set of equations,

(νλks(y)kλνs(y))P1k(a)Pα+1ν(I)=known terms with Pα and Pα+1(I) with II,\displaystyle\Big{(}\partial_{\nu}\lambda^{s}_{k}(y)-\partial_{k}\lambda^{s}_{\nu}(y)\Big{)}P_{1}^{k}(\partial_{a})P_{\alpha+1}^{\nu}(\partial_{I})=\text{known terms with $P_{\leq\alpha}$ and $P_{\alpha+1}(\partial_{I^{\prime}})$ with $I^{\prime}\prec I$,}

for each 1sp1\leq s\leq p, to the system (λα+1I\lambda^{I}_{\alpha+1}). That is, we have an affine system in Pα+1(I)P_{\alpha+1}(\partial_{I}) given as follows :

dλα+1I\boxed{d\lambda_{\alpha+1}^{I}} (λ11(y)λN1(y)λ1p(y)λNp(y)(νλ11(y)1λν1(y))P1ν(1)(νλN1(y)Nλν1(y))P1ν(1)(νλ1p(y)1λνp(y))P1ν(1)(νλNp(y)Nλνp(y))P1ν(1))(Pα+11(I)Pα+1N(I))=2p×1-vector of known terms involving Pα, and Pα+1(I) with II\begin{pmatrix}\lambda^{1}_{1}(y)&\ldots&\lambda^{1}_{N}(y)\\ \vdots&&\vdots\\ \lambda^{p}_{1}(y)&\ldots&\lambda^{p}_{N}(y)\\ \\ \big{(}\partial_{\nu}\lambda^{1}_{1}(y)-\partial_{1}\lambda^{1}_{\nu}(y)\big{)}P_{1}^{\nu}(\partial_{1})&\ldots&\big{(}\partial_{\nu}\lambda^{1}_{N}(y)-\partial_{N}\lambda^{1}_{\nu}(y)\big{)}P_{1}^{\nu}(\partial_{1})\\ \vdots&&\vdots\\ \big{(}\partial_{\nu}\lambda^{p}_{1}(y)-\partial_{1}\lambda^{p}_{\nu}(y)\big{)}P_{1}^{\nu}(\partial_{1})&\ldots&\big{(}\partial_{\nu}\lambda^{p}_{N}(y)-\partial_{N}\lambda^{p}_{\nu}(y)\big{)}P_{1}^{\nu}(\partial_{1})\end{pmatrix}\begin{pmatrix}P_{\alpha+1}^{1}(\partial_{I})\\ \vdots\\ P_{\alpha+1}^{N}(\partial_{I})\end{pmatrix}=\parbox{71.13188pt}{\centering\small A $2p\times 1$-vector of known terms involving $P_{\leq\alpha}$, and $P_{\alpha+1}(\partial_{I^{\prime}})$ with $I^{\prime}\prec I$\@add@centering}

Note that the rows of the (2p×N)(2p\times N)-sized coefficient matrix above are, respectively, the 11-forms on TyMT_{y}M,

λ1|y,,λp|y,ιP(1)dλ1|y,,ιP(1)dλp|y,\lambda^{1}\big{|}_{y},\ldots,\lambda^{p}\big{|}_{y},\quad-\iota_{P(\partial_{1})}d\lambda^{1}\big{|}_{y},\ldots,-\iota_{P(\partial_{1})}d\lambda^{p}\big{|}_{y},

written with respect to the basis {dy1,,dyN}\{dy^{1},\ldots,dy^{N}\}. Hence the matrix is full rank precisely when we have the wedge

(λ1λp)(ιP1(1)dλ1)(ιP1(1)dλp)\big{(}\lambda^{1}\wedge\ldots\wedge\lambda^{p}\big{)}\wedge\big{(}\iota_{P_{1}(\partial_{1})}d\lambda^{1}\big{)}\wedge\ldots\wedge\big{(}\iota_{P_{1}(\partial_{1})}d\lambda^{p}\big{)}

is non-zero. But since 𝒟\mathcal{D} is taken to be fat, for any vector 0v𝒟y0\neq v\in\mathcal{D}_{y}, and so in particular for v=P1(1)v=P_{1}(\partial_{1}), we have that

λ1λpιvdλ1ιvdλp0,\lambda^{1}\wedge\ldots\wedge\lambda^{p}\wedge\iota_{v}d\lambda^{1}\wedge\ldots\wedge\iota_{v}d\lambda^{p}\neq 0,

which follows from 2.4. Thus the coefficient matrix in the system (dλα+2Id\lambda^{I}_{\alpha+2}) above indeed has full rank. We can then inductively solve the tensor Pα+1P_{\alpha+1} completely. At this point, we have a lift of the jet σ\sigma to 𝒮α\mathcal{S}_{\alpha}.

Lastly, for each multi-index JJ of order α+2\alpha+2, we can easily solve Pα+2(J)P_{\alpha+2}(\partial_{J}) from the affine systems (λα+2J\lambda^{J}_{\alpha+2}), which has full rank coefficient matrix. Thus we have obtained a jet σ~α+1\tilde{\sigma}\in\mathcal{R}_{\alpha+1} so that p(σ~)=σp(\tilde{\sigma})=\sigma. Clearly the solution space p1(σ)p^{-1}(\sigma) is contractible, since at each stage we have solved affine system of equations, in a triangular fashion. Thus we have proved that 𝒫(α+1)\mathcal{P}(\alpha+1) is true.

This concludes the induction. Furthermore, it is clear from the algorithmic approach above that we can get lift of arbitrary sections of \mathcal{R} to α\mathcal{R}_{\alpha}, along the map pp, using local triviality arguments. Then, from the sheaf theoretic argument presented in [Gro86, pg. 76-78], we have that p:ΓαΓp:\Gamma\mathcal{R}_{\alpha}\to\Gamma\mathcal{R} is a weak homotopy equivalence. This completes the proof of 4.5.

5. Appendix : Hamilton’s Implicit Function Theorem

Nash’s Implicit Function Theorem [Nas56] in the context of CC^{\infty}-isometric immersions has been generalized by several authors. Here we recall Hamilton’s formalism of an infinite dimensional implicit function theorem that works for smooth differential operators between Fréchet spaces. This theorem is used crucially in order to get the local h-principle (Theorem 3.9) of horizontal maps into corank 22 fat distributions which admit Reeb directions. To begin with, we discuss the basic notions of tame spaces and tame operators from the exposition by Hamilton ([Ham82]).

Definition 5.1.

[Ham82, pg. 67] A Fréchet space is a complete, Hausdorff, metrizable, locally convex topological vector space.

In particular the topology of a Fréchet space FF is given by a countable collection of semi-norms {||n}\{|\cdot|_{n}\}, such that a sequence fjff_{j}\to f if and only if |fjf|n0|f_{j}-f|_{n}\to 0 for all nn, as jj\to\infty. A choice of this collection of norms is called a grading on the space and we say (F,{||n})(F,\{|\cdot|_{n}\}) is a graded Fréchet space.

Example 5.2.

Many naturally occurring spaces are in fact Fréchet spaces.

  1. (1)

    Every Banach space (X,||X)(X,|\cdot|_{X}) is a Fréchet space. It may also be graded if we set ||n=||X|\cdot|_{n}=|\cdot|_{X} for all nn ([Ham82, pg. 68]).

  2. (2)

    Given a compact manifold XX, possibly with boundary, the function space C(X)C^{\infty}(X) is a graded Fréchet space. More generally, given any vector bundle EXE\to X, the space of sections Γ(E)\Gamma(E) is also a graded Fréchet space. The CkC^{k}-norms on the sections give a possible grading ([Ham82, pg. 68]).

  3. (3)

    Given a Banach space (X,||X)(X,|\cdot|_{X}), denote by Σ(X)\Sigma(X) the space of exponentially decreasing sequences of XX, which consists of sequences {xk}\{x_{k}\} of elements of XX, such that,

    |{xk}|n=k=0enk|xk|X<,n0.|\{x_{k}\}|_{n}=\sum_{k=0}^{\infty}e^{nk}|x_{k}|_{X}<\infty,\quad\forall n\geq 0.

    Then Σ(X)\Sigma(X) is a graded Fréchet space with the norms defined above ([Ham82, pg. 134]).

Definition 5.3.

[Ham82, pg. 135] A linear map L:FGL:F\to G between Fréchet spaces F,GF,G is said to satisfy tame estimates of degree rr and base bb if there exists a constant c=c(n)c=c(n) such that,

|Lf|nC|f|n+r,nb,fF.|Lf|_{n}\leq C|f|_{n+r},\quad\forall n\geq b,\quad\forall f\in F.

LL is said to be tame if it satisfies the tame estimates for some nn and rr.

Example 5.4.

We have that a large class of operators are in fact tame.

  1. (1)

    A linear partial differential operator L:C(X)C(X)L:C^{\infty}(X)\to C^{\infty}(X) of order rr satisfies the tame estimate |Lu|n|u|n+r|Lu|_{n}\leq|u|_{n+r} for all n0n\geq 0 and hence LL is tame of degree rr ([Ham82, pg. 135]).

  2. (2)

    Inverses of elliptic, parabolic, hyperbolic and sub-elliptic operators are tame maps ([Ham82, pg. 67]). In particular, the solution of an elliptic boundary value problem is tame ([Ham82, pg. 161]).

  3. (3)

    Composition of two tame maps is again tame ([Ham82, pg. 136]).

Definition 5.5.

[Ham82, pg. 136] Given graded Fréchet spaces F,GF,G, we say FF is a tame direct summand of GG if there are tame linear maps L:FGL:F\to G and M:GFM:G\to F such that the composition ML:FFML:F\to F is the identity.

Definition 5.6.

[Ham82, pg. 136] A Fréchet space FF is said to be tame if FF is a tame direct summand of Σ(X)\Sigma(X), for some Banach space XX.

Example 5.7.

Given a compact manifold XX, possibly with boundary, and a vector bundle EXE\to X, the section space Γ(E)\Gamma(E) is a tame Fréchet space ([Ham82, pg. 139]).

Definition 5.8.

[Ham82, pg. 143] A map P:UFGP:U\subset F\to G between Fréchet spaces FF and GG, defined over some open set UFU\subset F, is said to be a smooth tame map if PP is smooth and all the derivatives DkPD^{k}P are tame linear maps.

We now state the inverse function theorem.

Theorem 5.9.

[Ham82, pg. 171] Consider tame Fréchet spaces F,GF,G and a tame smooth map P:UGP:U\to G, where UFU\subset F is open. Suppose for the derivative DP(f)DP(f) at fUf\in U, the equation DP(f)h=kDP(f)h=k admits unique solution h=VP(f)kh=VP(f)k for each kGk\in G. Furthermore, assume that VP:U×GFVP:U\times G\to F is a smooth tame map. Then PP is locally invertible and each local inverse P1P^{-1} is smooth tame.

Remark 5.10.

Unlike the inverse function theorem for Banach spaces, one needs to have that the derivative DPDP is invertible on an open set UFU\subset F.

Acknowledgment

The author would like to thank Mahuya Datta, Adi Adimurthi and Partha Sarathi Chakraborty for fruitful and enlightening discussions. The author would also like to thank the anonymous referee for many valuable comments and suggestions.

References

  • [AFL17] Antonio Alarcón, Franc Forstnerič, and Francisco J. López. Holomorphic Legendrian curves. Compositio Mathematica, 153(9):1945–1986, 2017. doi:10.1112/S0010437X1700731X.
  • [BH05] Gianluca Bande and Amine Hadjar. Contact pairs. The Tohoku Mathematical Journal. Second Series, 57(2):247–260, 2005.
  • [CdS01] Ana Cannas da Silva. Lectures on symplectic geometry, volume 1764 of Lecture Notes in Mathematics. Springer-Verlag, Berlin, 2001. doi:10.1007/978-3-540-45330-7.
  • [CFS05] Sergio Console, Anna Fino, and Evangelia Samiou. The moduli space of six-dimensional two-step nilpotent Lie algebras. Annals of Global Analysis and Geometry, 27(1):17–32, 2005. doi:10.1007/s10455-005-2569-2.
  • [Cho39] Wei-Liang Chow. Über Systeme von linearen partiellen Differentialgleichungen erster Ordnung. Mathematische Annalen, 117:98–105, 1939. doi:10.1007/BF01450011.
  • [D’A94] Giuseppina D’Ambra. Induced subbundles and Nash’s implicit function theorem. Differential Geometry and its Applications, 4(1):91–105, 1994. doi:10.1016/0926-2245(94)00003-4.
  • [Duc84] Tom Duchamp. The classification of Legendre immersions. 1984.
  • [EM02] Y. Eliashberg and N. Mishachev. Introduction to the hh-principle, volume 48 of Graduate Studies in Mathematics. American Mathematical Society, Providence, RI, 2002. doi:10.1090/gsm/048.
  • [FL18] Franc Forstnerič and Finnur Lárusson. The Oka principle for holomorphic Legendrian curves in 2n+1\mathbb{C}^{2n+1}. Mathematische Zeitschrift, 288(1-2):643–663, 2018. doi:10.1007/s00209-017-1904-1.
  • [Ge92] Zhong Ge. Betti numbers, characteristic classes and sub-Riemannian geometry. Illinois Journal of Mathematics, 36(3):372–403, 1992. doi:10.1215/ijm/1255987416.
  • [Ge93] Zhong Ge. Horizontal path spaces and Carnot-Carathéodory metrics. Pacific Journal of Mathematics, 161(2):255–286, 1993. doi:10.2140/PJM.1993.161.255.
  • [Gei08] Hansjörg Geiges. An introduction to contact topology, volume 109 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 2008. doi:10.1017/CBO9780511611438.
  • [Gro86] Mikhael Gromov. Partial differential relations, volume 9 of Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. Springer-Verlag, Berlin, 1986. doi:10.1007/978-3-662-02267-2.
  • [Gro96] Mikhael Gromov. Carnot-Carathéodory spaces seen from within. In Sub-Riemannian geometry, volume 144 of Progr. Math., pages 79–323. Birkhäuser, Basel, 1996. URL: https://www.ihes.fr/~gromov/wp-content/uploads/2018/08/carnot_caratheodory.pdf.
  • [Ham82] Richard S. Hamilton. The inverse function theorem of Nash and Moser. American Mathematical Society. Bulletin. New Series, 7(1):65–222, 1982. doi:10.1090/S0273-0979-1982-15004-2.
  • [Mon93] Richard Montgomery. Generic distributions and Lie algebras of vector fields. Journal of Differential Equations, 103(2):387–393, 1993. doi:10.1006/jdeq.1993.1056.
  • [Mon02] Richard Montgomery. A tour of subriemannian geometries, their geodesics and applications, volume 91 of Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 2002. doi:10.1090/surv/091.
  • [Nas56] John Nash. The imbedding problem for Riemannian manifolds. Annals of Mathematics. Second Series, 63:20–63, 1956. doi:10.2307/1969989.
  • [Pan16] Pierre Pansu. Submanifolds and differential forms on carnot manifolds, after m. gromov and m. rumin. 2016. arXiv:1604.06333.
  • [Ray68] C. B. Rayner. The exponential map for the Lagrange problem on differentiable manifolds. Philosophical Transactions of the Royal Society of London. Series A. Mathematical and Physical Sciences, 262:299–344, 1967/68. doi:10.1098/rsta.1967.0052.
  • [Tan70] Noboru Tanaka. On differential systems, graded Lie algebras and pseudogroups. Journal of Mathematics of Kyoto University, 10:1–82, 1970. doi:10.1215/kjm/1250523814.
  • [Zan15] Chiara De Zanet. Generic one-step bracket-generating distributions of rank four. Archivum Mathematicum, (5):257–264, 2015. doi:10.5817/am2015-5-257.