Flows in Flatland: A Romance of Few Dimensions
Abstract
This paper is about gradientlike vector fields and flows they generate on smooth compact surfaces with boundary. We use this particular 2dimensional setting to present and explain our general results about nonvanishing gradientlike vector fields on $n$dimensional manifolds with boundary. We take advantage of the relative simplicity of 2dimensional worlds to popularize our approach to the Morse theory on smooth manifolds with boundary. In this approach, the boundary effects take the central stage.
Keywords
1 Introduction
This paper is about the gradient flows on compact 2dimensional surfaces, thus the reference to [Abbott1992]^{1}^{1}footnotemark: ^{2}^{2}In [Abbott1992], the drama unfolds in the Euclidean plane, the Flatland, while here we will tell a story that takes place in the 2Dland, a compact Riemannian surface. in the title. In other words, we are interested in solutions of ordinary differential equations (ODE’s) on compact surfaces $X$ with boundary. The ODE’s are represented by smooth vector fields $v$ on $X$. We assume that $v$ admits a global potential $f\colon X\to\mathbb{R}$ such that $df(v)>0$ away from the points $x$ where $v(x)=0$.
The article is an informal introduction into the philosophy and some key results from [Katz2009]; [Katz2014a]; [Katz2014b]; [Katz2014c]; [Katz2015]; [Katz2016a]; [Katz2016b]; [Katz2016c], as they manifest themselves in the dimension two. Our strongest results deal with so called traversing vector fields on smooth compact manifolds $X$ with boundary. (By Lemma 5.1, these are exactly the vector fields that are gradients of smooth functions $f\colon X\to\mathbb{R}$ with no critical points.)
Definition 1.1.
We say that a nonvanishing vector field $v$ on a compact surface/smooth manifold $X$ is traversing if all its trajectories are closed segments or singletons; the ends of these segments, as well as the singletons, reside in the boundary $\partial X$.
Thus, for a traversing vector field $v$, no trajectory belongs to the interior of $X$: it reaches the boundary $\partial X$ in both “finite positive and negative times”.
The main goal of our efforts is to understand the connection between the geometry of the integral curves of the $v$flow, as they interact with the boundary $\partial X$ of a given manifold $X$, and the topology of $X$. We are especially interested in describing the restrictions, imposed on $X$ by assuming that it admits a traversing vector field $v$ whose trajectories $\gamma$ intersect $\partial X$ in such a way that the resulting combinatorial tangency patterns (see Definition 7.1 and Lemma 7.1) are prescribed a priory.
Along the way, we produce various invariants that distinguish between such traversing flows on $X$^{3}^{3}The flows and the vector fields here are considered up to the natural action of the diffeomorphism group $\mathsf{Diff}(X)$ on them.. The invariants reflect complexities of flowinduced natural stratifications (see Fig. 5) of $X$. Then we employ these flow invariants to manufacture a variety of new smooth topological invariants of $X$ itself. In [Katz2009]; [Katz2014a]; [Katz2014b]; [Katz2014c]; [Katz2015]; [Katz2016a]; [Katz2016b]; [Katz2016c] and [Alpert and Katz2015], we use a variety of standard algebraic and geometric topology tools to compute or to estimate from below the new invariants. However, for this exposition, the reader is expected to be familiar only with the homology of compact surfaces and with their fundamental groups. We also assume some familiarity with the basic notions of the Morse Theory.
Although the class of manifolds with boundary is wider than the class of closed manifolds, the vector fields of the kind that we investigate do not exist on closed manifolds (a crude analogy of this phenomenon is provided by the Dirichlet boundary problems, which make little sense on closed Riemannian surfaces). Therefore it worth stressing that most of our results are not generalizations of some facts, valid for closed manifolds. In fact, often we drill holes in a given closed manifold $Y$ to manufacture a manifold $X$ with boundary and use some flowgenerated invariants of $X$ to produce new invariants of the closed $Y$.
Let us now concentrate on the flows on 2dimensional surfaces, the main subject of this paper. The remarkable convergence of topological, geometrical, and analytical approaches to the study of closed surfaces is widely recognized by the practitioners for more than a century. We will exhibit a similar convergence of different investigative approaches to gradient flows on surfaces with boundary. Because the “rigid” topology of surfaces is understood so well, no new invariants of surfaces will emerge; however, we will get new expressions of classical invariants in terms of the 2dimensional flows (see Lemma 3.1, Theorem 3.1, Corollaries 4.1, 4.2).
In particular, we take advantage of the relative simplicity of 2dimensional flows to illustrate and popularize the main ideas of our recent research of, so called, traversally generic flows (see Definition 5.1 and Definition 3.2 from [Katz2014b] for the multidimensional case) on manifolds with boundary.
When the results are specific to dimension two, their validation will be presented in detail. The multidimensional arguments that resist significant simplifications in dimension two will be described and explained in general terms.
Generically (that is, for a dense and open set in the space of all traversing vector fields on a given surface), where the flow line is tangent to the boundary of a surface, the boundary can be either convex or concave with respect to the flow (see Fig. 1). In particular, in dimension two, for a generic $v$, no $v$trajectory is cubically tangent to the boundary.
In what follows, we are motivated by few natural questions.
Question 1.1.
Given a compact connected surface $X$ with boundary, what is the minimal number of concave/convex tangency points that a gradient or a traversing vector flow on $X$ may have? How these minima can be expressed in terms of the topological invariants of $X$?
Lemma 3.1, Corollaries 4.1, 4.2 provide answers to Question 1.1. Theorem 1.1 below describes the flavor of these results. It is implied by Corollary 4.2 and Lemma 5.1.
Theorem 1.1.
Let a surface $X$ be formed by removing $k$ open disks from a closed orientable surface $Y$, the sphere with $g$ handles.
Then, for an open and dense set in the space of all traversing vector fields $v$ on $X$, the number of concave $v$flow tangencies to $\partial X$ is greater than or equal to $4g4+2k$.
Moreover, this inequality turns into an equality if and only if the trajectories of $v$ are tangent to $\partial X$ in a concave manner only (as in Fig. 1, the right diagram).
We denote by $gc^{+}(X)$ the minimal number of concave tangencies to $\partial X$ that a generic traversing vector flow on $X$ may have, and view $gc^{+}(X)$ as a measure of complexity of $X$. In fact, a similar invariant makes sense for $n$dimensional manifolds $X$: it is based on counting the flow trajectories that are simply tangent (see Definition 7.1) to $\partial X$ exactly $n1$ times, and then taking the minimum of these counts over all such flows on $X$ ([Alpert and Katz2015]; [Katz2016b]).
The next question is a close relative of Question 1.1.
Question 1.2.
Which compact connected surfaces with boundary admit traversing vector flows that have only concave/convex tangencies?
The next theorem provides a complete answer.
Theorem 1.2.

•
Any connected compact surface $X$ with boundary, but the disk and the Möbius band, admits a traversing vector field $v$ with only concave tangency points to $\partial X$.

•
The disk and the Möbius band do not admit such a vector field.

•
Any connected compact surface with boundary, that admits a traversing vector field with convex tangency points only is a disk or an annulus^{4}^{4}In the case of annulus, this set is empty.. So, the annulus is the only connected compact surface that admits both concave and convex traversing vector fields.
Theorems 6.1 and 6.2 generalize Theorem 1.1. They deliver more nuanced answers to Questions 1.1 and 1.2.
Now we are in position to formulate a “more existential” metaquestion:
Question 1.3.
Given a traversing vector field $v$ on a connected compact surface/ manifold $X$, are there some $v$flow generated structures on the boundary $\partial X$ that will allow for (perhaps, partial) reconstruction of $X$ and the flow?
We call such desirable structures holographic since the information, recorded on 1dimensional [in general, $(n1)$dimensional] screen, makes it possible to restore the topology of $X$ and the 2dimensional [in general, $n$dimensional] flow on it.
Let us describe one holographic structure that dominates the landscape. Take a compact connected surface $X$ with boundary and a traversing vector field $v$ on it. For each point $x\in\partial X$, where $v$ points inward of $X$ or is tangent to $\partial X$, consider the the $v$trajectory $\gamma$ through $x$. Let $C_{v}(x)\in\partial X$ be the first along $\gamma$ point, where $\gamma$ exits $X$ or touches its boundary. We think of the pointevent $x$ as the cause for the pointevent $C_{v}(x)$. This interpretation will be justified in Sect. 8, Example 8.1.
The correspondence $x\leadsto C_{v}(x)$ gives rise to a welldefined map
$\displaystyle C_{v}:\partial^{+}X(v)\longrightarrow\partial^{}X(v)$ 
from a portion $\partial^{+}X(v)$ of $\partial X$, where $v$ points inward of $X$ or is tangent to $\partial X$, to the closure $\partial^{}X(v)$ of the complementary portion $\partial X\setminus\partial^{+}X(v)$. For a generic^{5}^{5}That is, for an open and dense set in the space of all traversing vector fields on $X$. $v$, both portions share a common boundary: $\partial(\partial^{+}X(v))=\partial(\partial^{}X(v))$.
We may view $C_{v}$ as defining a partial order among points of $\partial X$ (see Sect. 8 for details).
The map $C_{v}$ is discontinuous; in fact, we assign the most important role to the discontinuities of $C_{v}$. (See Fig. 7 for a pictorial example of $C_{v}$.)
To get a feel for the “physics inspired” spirit of Question 1.3, imagine that we do not have access to the interior of a mystery surface $X$, but have positioned a set of detectors along its boundary, detectors that register where each particle of a traversing $v$flow enters $X$ and where its first exit point on the boundary $\partial X$ is located. The causality map $C_{v}$ records all such data. So we asking whether the data are sufficient for a reconstruction of the pair $(X,v)$.
Theorem 8.1 is a generalization of the following theorem.
Theorem 1.3.
(The Causal Holography Theorem) For an open and dense set in the space of all traversing vector fields $v$ on a compact connected surface $X$ with boundary, the causality map $C_{v}\colon\partial^{+}X(v)\to\partial^{}X(v)$ determines the topology of $X$ and the $v$oriented unparametrized $v$flow curves on it, up to a diffeomorphism of $X$ which is the identity on $\partial X$.
So Theorem 1.3 gives a positive answer to Question 1.3. Of course, if a surface $X$ and a vector field $v$ on it are real analytic and the boundary of $\partial X$ is analytic as well, then one should expect some rigid dependence of the flow in $X$ on the boundary data. Somewhat surprisingly, Theorem 1.3 tells us that, in the smooth category, a similar topological rigidity is valid!
2 On Morse Theory on Surfaces with Boundary and Beyond
In this section, we present our approach to the Morse Theory of gradient flows, in which the boundary effects (like the convexity and the concavity) rule.
Morse Theory, the classical book of [Milnor1965], starts with the canonical picture of a Morse function $f\colon T^{2}\to\mathbb{R}$ on a 2dimensional torus $T^{2}$ (see Fig. 2). It is portrayed as the height function $f$ on the torus $T^{2}$ residing in the space $\mathbb{R}^{3}$. The height $f$ has four critical points: $a,b,c$, and $d$ so that
$\displaystyle f(a)>f(b)>f(c)>f(d).$ 
Recall that a point $z$ is called critical if the differential $df$ of $f$ vanishes at $z$. In the vicinity of each critical point $z$, $T^{2}$ admits a pair of local coordinate functions, say $x$ and $y$, so that locally the function $f$ acquires the form
$\displaystyle f(x,y)=f(0,0)\pm x^{2}\pm y^{2},$ 
where the signs may form four possible combinations.
We call a vector field $v$, tangent to $T^{2}$, $f$gradientlike if $df(v)>0$ everywhere outside of the set $Cr(f)$ of critical points.
If the torus is “slightly slanted” with respect to the vertical coordinate $f$ in $\mathbb{R}^{3}$, then the following picture emerges. The majority of downward trajectories of the $f$gradientlike flow $\{\Phi_{t}\}_{t\in\mathbb{R}}$ that emanate from $a$, asymptotically reach $d$. There are two trajectories that asymptotically link $a$ with $b$, and two trajectories that link $a$ with $c$. No (unbroken) trajectory asymptotically connects $b$ to $c$.
Perhaps, a more transparent depiction of the gradientlike flow $\{\Phi_{t}\}_{t\in\mathbb{R}}$ is given in Fig. 3, where the torus is shown in terms of its fundamental domain, the square. To form $T^{2}$, the opposite sides of the square are identified in pairs.
The Morse Theory is concerned with the sets of constant level $\{f^{1}(\alpha)\}_{\alpha\in\mathbb{R}}$ and with the below constant level sets $\{f^{1}((\infty,\alpha))\}_{\alpha\in\mathbb{R}}$. The main observation is that the topology of these sets is changing in an essential way only when the rising level $\alpha$ crosses the critical values
$\displaystyle Cr(f)=\{f(a),f(b),f(c),f(d)\}.$ 
In the rest of this section, we assume that all the critical values are distinct.
Each such “critical crossing” results in an elementary surgery on the set $\{f^{1}((\infty,\alpha))\}_{\alpha\in\mathbb{R}}$, where $\alpha$ is just below a critical value $\alpha_{\star}\in Cr(f)$. For a small $\epsilon>0$, an elementary surgery
$\displaystyle f^{1}((\infty,\alpha_{\star}\epsilon))\Rightarrow f^{1}(( \infty,\alpha_{\star}+\epsilon))$ 
attaches the handle $f^{1}((\alpha_{\star}\epsilon,\alpha_{\star}+\epsilon))$ to the set $f^{1}((\infty,\alpha\epsilon))$. Eventually, when $\alpha$ rises above $f(a)$, the entire topology of torus $T^{2}$ is captured by a sequence of these elementary surgeries.
From a different angle, the knowledge of how the critical points $a,b,c,d$ interact via the trajectories of the $\Phi_{t}$flow is also sufficient for reconstructing the surface $T^{2}$, as the left diagram in Fig. 3 suggests [see [Cohen1991] for multidimensional generalizations].
Note that, in the vicinity of each critical point, the gradient flow exhibits discontinuity: small changes in the initial position of a point $z$, residing in the vicinity of a critical point, result in significant differences in the position of $\Phi_{t}(z)$ for big positive/small negative values of $t$ (see Fig. 4). In fact, this discontinuity of the gradient flow, expressed in terms of the stable and unstable manifolds of critical points [see [Milnor1965]], captures the topology of the surface!
As a result of gradient flow discontinuity, the space of trajectories $\mathcal{T}(v)$ is pathological (nonseparable). The space $\mathcal{T}(v)$ is constructed by declaring equivalent any two points that reside on the same trajectory.
Originally, the Morse theory of smooth functions and their gradient flows on compact manifolds $X$ with boundary has been studied (under general boundary conditions) in the papers of [Morse1929] and of [Morse and Van Schaack1934]. See [Barannikov1994] for the case $X$ being a ball.
The Morse Theory on manifolds with boundary can be viewed as a very special instance of the Morse Theory on stratified spaces (the two strata $\partial X$ and $X$ form the stratification). The latter was developed by Goresky and MacPherson in [Goresky and MacPherson1983a]; [Goresky and MacPherson1983b]; [Goresky and MacPherson1989]. See also [Hardon2004], [Laudenbach2011] for constructions of the MorseSmale differential complex on manifolds with conners.
If the Morse function $f\colon X\to\mathbb{R}$ is assumed to be constant on $\partial X$, then its gradient flow interacts with the boundary in a constrained way. Then the relative topology of the pair $(X,\partial X)$ can be captured in the ways analogous to the previous description of the Morse Theory on torus.
In this paper, we offer a different philosophy for the Morse Theory on compact surfaces/manifolds $X$ with boundary. To formulate it, let us revisit our favorite closed surface, the torus. By deleting from $T^{2}$ small disks, centered on the points of the critical set $Cr(f)$, we manufacture a surface $X$ whose boundary is a disjoint union of four circles. Evidently, $f\colon X\to\mathbb{R}$ has no critical points at all. Still it has a nontrivial topology! Can this topology be reconstructed from some data, provided by the criticalpointfree $f$ and its nonvanishing gradientlike vector field $v$? An experienced reader would notice that the restriction $f_{\partial X}$ has critical points (maxima and minima), some of which interact along the boundary (with the help of a gradientlike vector field $v^{\partial}$, tangent to $\partial X$). However, it is quite clear that these interactions are not sufficient for a reconstruction of the topology of $X$! In fact, a reconstruction of the surface $X$ becomes possible if one introduces additional interactions between the points of $Cr(f_{\partial X})$ that occur “through the bulk $X$” and are defined with the help of both vector fields $v$ and $v^{\partial}$. This observation has been explored by a number of authors ([Hardon2004]; [Laudenbach2011]), but it is not the world view that we are promoting here.
To dramatize further the situation we are facing, let us place four small disks, centered on the critical points of $f\colon T^{2}\to\mathbb{R}$, into a single open disk $D^{2}$ and form $X=T^{2}{\setminus}D^{2}$ (see Fig. 3, the right diagram). Again, $f_{X}\colon X\to\mathbb{R}$ has no critical points, the gradientlike vector field is nonvanishing on $X$, but its topology of $X$ is nontrivial. This time, the boundary $\partial X$ of the punctured torus $X$ is just a single circle! Let us keep this challenge in mind.
So can one propose a “Morse Theory” that is not centered on critical points? The answer is affirmative. It relies on the following observation. Typically, in the vicinity of $\partial X$, the $v$trajectories are interacting with the boundary in a number of very particular and stable ways: they are either transversal to $\partial X$, or are tangent to it in a concave or convex fashion^{6}^{6}It is possible to have a vector field $v$ for which some trajectories will be cubically tangent to the boundary, but the majority of 2dimensional vector fields $v$ avoid such cubic tangencies [Katz2014b]. (see Fig. 1). So the boundary $X$ may be “wiggly” with respect to the flow. We claim that this geometry of the $v$flow in relation to the boundary $\partial X$ is the crucial ingredient for reconstructions of $X$ in terms of the flow (see Sect. 8, especially Theorem 8.1).
In the vicinity of a concave tangency point, the $v$flow is discontinuous in the same sense as the gradient flow is discontinuous in the vicinity of its critical point: in time, close initial points become distant. In this context, the divergence of initially close points occurs due to very different travel times available to them; unlike the infinite travel time for the gradient flows of the Morse theory on closed surfaces, in the case of the nonsingular gradient flows on surfaces with boundary, every point exits the surface in finite time. In particular, the surface is not flowinvariant! And again, these discontinuities of the flow reflect the topology of the surface. Let us clarify this point.
Figure 5 shows a gradient flow $v$ on a surface $X\subset\mathbb{R}^{2}$, the disk with 4 holes. The nonsingular function $f\colon X\to\mathbb{R}$ is the vertical coordinate in $\mathbb{R}^{2}$. Each $v$trajectory is either a closed segment, or a singleton; so the flow is traversing. By collapsing each trajectory to a point, we create a quotient space $\mathcal{T}(v)$ of trajectories. Since the flow trajectories are closed segments or singletons, this time, the trajectory space $\mathcal{T}(v)$ is “decent”, a finite graph with verticies of valency 1 or 3 only. The verticies of valency 3 correspond to the points on $\partial X$ where the boundary is concave with respect to the flow, and the univalent verticies to the points on $\partial X$ where the flow is convex.
The obvious map $\Gamma\colon X\to\mathcal{T}(v)$ is cellular. Moreover, because the fibers of $\Gamma$ are contractible, $\Gamma$ is a homotopy equivalence. In particular, the fundamental groups $\pi_{1}(X)$ and $\pi_{1}(\mathcal{T}(v))$ are isomorphic with the help of $\Gamma$. So the trajectory spaces of generic nonvanishing gradientlike vector fields $v$ on connected surfaces $X$ with boundary deliver 1dimensional homotopy theoretical models of $X$. This observation is valid for generic $n$dimensional nonvanishing gradientlike flows [Katz2014c], but its justification requires an involved argument.
3 Vector Fields and Morse Stratifications on Surfaces
To achieve some uniformity of notations, put $\partial_{0}^{+}X=_{\mathsf{def}}X$ and $\partial_{1}X=_{\mathsf{def}}\partial X$.
Following [Morse1929], for any vector field $v$ on a compact surface $X$, such that the vector field does not vanish at the points of $\partial X$, we consider the closed locus $\partial_{1}^{+}X(v)$, where the vector field is pointing inside $X$ or is tangent to $\partial X$, and the closed locus $\partial_{1}^{}X(v)$, where it points outside or is tangent to $\partial X$. The intersection
$\displaystyle\partial_{2}X(v)=_{\mathsf{def}}\partial_{1}^{+}X(v)\cap\partial_ {1}^{+}X(v)$ 
is the locus where $v$ is tangent to the boundary $\partial X$. Points $z\in\partial_{2}X(v)$ come in two flavors: by definition, $z\in\partial^{+}_{2}X(v)$ when $v(z)$ points inside of the locus $\partial_{1}^{+}X(v)$; otherwise $z\in\partial^{}_{2}X(v)$.
Definition 3.1.
We say that a vector field $v$ on a compact surface $X$ is boundary generic if:

•
the vector field $v_{\partial X}$ does not vanish,

•
$v_{\partial X}$, viewed as a section of the normal 1dimensional (quotient) bundle
$\displaystyle n_{1}=_{\mathsf{def}}T(X)_{\partial X}\big/T(\partial X),$ is transversal to its zero section (at the points of the locus $\partial_{2}X(v)$).
Thus, for a boundary generic $v$, the loci $\partial_{1}^{\pm}X(v)$ are finite unions of closed intervals and circles, residing in $\partial X$; and the loci $\partial_{2}^{\pm}X(v)$ are finite unions of points in $\partial X$ (Fig. 5).
We denote by $\mathcal{V}^{\dagger}(X)$ the space of all boundary generic vector fields on a compact surface $X$. It is easy to check that $\mathcal{V}^{\dagger}(X)$ is an open and dense (in the $C^{\infty}$topology) subset of the space $\mathcal{V}(X)$ of all vector fields on $X$ ([Katz2014a]).
Let $\chi(Z)$ denote the Euler number of a space $Z$. Recall that $\chi(Z)$ is the alternating sum of dimensions of the homology spaces $\{H_{i}(Z;\mathbb{R})\}_{i}$.
Since for a connected surface $X$ with boundary $H_{2}(X;\mathbb{R})=0$, we get
$\displaystyle\chi(X)=1\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R})).$ 
For a closed connected surface,
$\displaystyle\chi(X)=2\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R})).$ 
Given a vector field $v$ with isolated zeros, we can associate an integer $\mathsf{ind}_{x}(v)$ with each zero $x$ of $v$. This integer is the degree of the map which, crudely speaking, takes each point $z$ on a small circle $C_{x}$ with its center at $x$ to the unit vector $v(z)/\v(z)\$. Then we define $\mathsf{Ind}(v)$, the (global) index of $v$, as the sum $\sum_{\{x\in\text{ zeros of }v\}}\mathsf{ind}_{x}(v)$.
The Morse formula ([Morse1929]), in the center of our investigation, computes the index $\mathsf{Ind}(v)$ of a given boundary generic vector field $v$ on a surface $X$ as the alternating sum of the Euler numbers of the Morse strata $\{\partial_{j}^{+}X(v)\}_{0\leq j\leq 2}$:
$\displaystyle\mathsf{Ind}(v)=\chi(X)\chi(\partial_{1}^{+}X(v))+\chi(\partial_ {2}^{+}X(v)).$  (3.1) 
In the case of a connected surface $X$ with boundary, $\chi(X)=1\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R}))$, and this formula reduces to
$\displaystyle\mathsf{Ind}(v)$  $\displaystyle=$  $\displaystyle 1\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R}))\#\{\mathsf{arcs}\,\mathsf{in}\; \partial_{1}^{+}X(v)\}+\#\{\partial_{2}^{+}X(v)\}$  
$\displaystyle=$  $\displaystyle 1\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R}))+\frac{1}{2}\big(\#\{\partial_{2}^{+}X(v )\}\#\{\partial_{2}^{}X(v)\}\big).$ 
In particular, if $v$ is a nonvanishing field in $X$, then $\mathsf{Ind}(v)=0$, and we get
$\displaystyle\frac{1}{2}\big(\#\{\partial_{2}^{+}X(v)\}\#\{\partial_{2}^{}X( v)\}\big)=\dim_{\mathbb{R}}(H_{1}(X;\mathbb{R}))1,$  (3.2) 
where the RHS of equality (3.2) is the topological invariant $\chi(X)$ of $X$. In contrast, the cardinality $\#\{\partial_{2}^{+}X(v)\}$ depends on $v$.
Lemma 3.1.
Let a surface $X$ be formed by removing $k$ open disks from a closed surface $Y$, the sphere with $g$ handles. Then, for any boundary generic nonvanishing vector field $v$ on $X$,
$\displaystyle\#\{\partial_{2}^{+}X(v)\}\,\geq\,4g4+2k.$ 
Moreover, $\#\{\partial_{2}^{+}X(v)\}=4g4+2k$ only when $\#\{\partial_{2}^{}X(v)\}=0$.
Proof.
The Euler number is additive under gluing surfaces along their boundary components. Therefore, if $k$ disks are removed from $Y$, the sphere with $g$ handles, then $\chi(X)=22gk$. Thus the formula (3.2) implies
$\displaystyle\#\{\partial_{2}^{+}X(v)\}\geq 4g4+2k$ 
for any nonvanishing $v$. Moreover, $\#\{\partial_{2}^{+}X(v)\}=4g4+2k$ if and only if $\#\{\partial_{2}^{}X(v)\}=0$, the main feature of the boundary concave vector fields (see Definition 4.1). $\square$
In particular, for any nonvanishing boundary generic vector field $v$ on a torus with a single hole, $\#\{\partial_{2}^{+}X(v)\}\geq 2$ (cf. Fig. 3).
Recall that an immersion is a smooth map of manifolds, whose differential has the trivial kernel.
Consider a smooth map $\alpha\colon X\to\mathbb{R}^{2}$, which is an immersion in the vicinity of $\partial X$. Any such $\alpha$ gives rise to the Gauss map $G\colon\partial X\to S^{1}$, defined by the formula $G(x)=\alpha_{\ast}(\tau_{x})/\\alpha_{\ast}(\tau_{x})\$, where $\tau_{x}$ is the tangent vector to $\partial X$ at $x$. The direction of $\tau_{x}$ is consistent with the orientation of $\partial X$, induced by the preferred orientation of $X$.
Let $\hat{v}$ be a nonvanishing vector field on $\mathbb{R}^{2}$. Since the kernel of the differential of $D\alpha\colon TX\to T\mathbb{R}^{2}$ is trivial along $\partial X$, the vector field $\hat{v}$ defines a vector field $\tilde{v}=\alpha^{\ast}(\hat{v})$ on $X$ in the vicinity of $\partial X$. The pullback vector field $\tilde{v}$ extends to a vector field $v$ on $X$, possibly with zeros (see [Gottlieb1996] for engaging discussions of vector field transfers and the GaussBonnet Theorem).
Then the degree of the Gauss map is given by a classical Hopf formula ([Hopf1937])
$\displaystyle\deg(G)=\chi(X)\mathsf{Ind}(v).$ 
When $\alpha\colon X\to\mathbb{R}^{2}$ is an immersion everywhere, the pullback vector field $v=\alpha^{\ast}(\hat{v})\neq 0$ everywhere. Thus $\mathsf{Ind}(v)=0$, and, for a connected $X$ with $\partial X\neq\emptyset$, we get
$\displaystyle\deg(G)=\chi(X)=1\dim(H_{1}(X;\mathbb{R})).$ 
So, for immersions $\alpha$, we get a new interpretation of formula (3.2):
$\displaystyle\deg(G)=\chi(X)=\frac{1}{2}\Big(\#\{\partial_{2}^{}X(v)\}\#\{ \partial_{2}^{+}X(v)\}\Big).$  (3.3) 
This globaltolocal formula has another classical geometrical interpretation. Let $\mathsf{g}=\alpha^{\ast}(\mathsf{g}_{E})$ be the Riemannian metric on $X$, the pullback of the Euclidean metric on $\mathbb{R}^{2}$. Let $K_{\nu}$ denote the normal curvature of $\partial X$ with respect to $\mathsf{g}$. Then
$\displaystyle\deg(G)=\frac{1}{2\pi}\int_{\partial X}K_{\nu}\,d\mathsf{g},$ 
which leads to another pleasing globaltolocal connection:
$\displaystyle\frac{1}{\pi}\int_{\partial X}K_{\nu}\,d\mathsf{g}=\#\{\partial_{ 2}^{}X(v)\}\#\{\partial_{2}^{+}X(v)\}.$ 
In particular, for a connected orientable surface $X$ of genus $g$ with a single boundary component,
$\displaystyle\chi(X)=12g=\frac{1}{2}\big(\#\{\partial_{2}^{}X(v)\}\#\{ \partial_{2}^{+}X(v)\}\big).$  (3.4) 
So the number of $v$trajectories $\gamma$ in $X$ that are tangent to $\partial X$, but are not singletons (they correspond to points of $\partial_{2}^{+}X(v)$), as a function of genus $g$, is greater than or equal to $4g2$!
On the other hand, when $\partial X$ is connected, by the Whitney index formula ([Whitney1937]), the degree of the Gauss map $G\colon\partial X\to S^{1}$ can be also calculated as $\mu+N^{+}N^{}$, where $N^{\pm}$ denotes the number of positive/negative selfintersections of the curve $\alpha(\partial X)\subset\mathbb{R}^{2}$, and $\mu=\pm 1$. Here is a brief description of the rule by which these selfintersections acquire polarities. Let $p\in\alpha(\partial X)$ be a point where the coordinate function $y\colon\mathbb{R}^{2}\to\mathbb{R}$ attends its minimum on the curve $\alpha(\partial X)$. If the tangent vector $\tau_{p}$ at $p$, which defines the orientation of $\alpha(\partial X)$, is $\partial_{x}$, then we put $\mu=+1$; if $\tau_{p}=\partial_{x}$, then $\mu=1$. Starting at $p$ and moving in the direction of $\tau_{p}$, we visit each selfintersection $a$ twice and in a particular order. The first visitation defines a tangent vector $\tau_{1}(a)$, the second visitation defines a tangent vector $\tau_{2}(a)$. When the ordered pair $(\tau_{1}(a),\tau_{2}(a))$ defines the clockwise orientation of the $xy$plane, then we attach “$$” to $a$. Otherwise, the polarity of $a$ is “$+$”.
Therefore we get a somewhat mysterious connection between the selfintersections of $\partial X$ under immersions $\alpha\colon X\to\mathbb{R}^{2}$ and the tangency patterns of the flows in $X$ that are the $\alpha$pullbacks of nonvanishing flows in the plane.
Theorem 3.1.
Let $\hat{v}$ be a nonvanishing vector field in the plane $\mathbb{R}^{2}$. Let $X$ be a connected orientable surface of genus $g$ and with a connected boundary. Consider an immersion $\alpha\colon X\to\mathbb{R}^{2}$ such that the map $\alpha\colon\partial X\to\mathbb{R}^{2}$ is 2to1 at most and the loop $\alpha(\partial X)$ has transversal selfintersections only. Assume that the pullback $v=\alpha^{\ast}(\hat{v})$ is a boundary generic vector field on $X$. Then
$\displaystyle\frac{1}{2}\Big(\#\{\partial_{2}^{+}X(v)\}\#\{\partial_{2}^{}X( v)\}\Big)=N^{+}N^{}\pm 1=2g1,$ 
$\displaystyle\frac{1}{2}\Big(\#\{\partial_{2}^{+}X(v)\}\#\{\partial_{2}^{}X( v)\}\Big)+2\,\leq N^{+}+N^{},$ 
the latter inequality being sharp by an appropriate choice of $\alpha$.
Proof.
The first formula is the result of combining the Whitney formula for $\deg(G)$ with formulas (3.3), (3.4).
By a theorem of [Guth2009], for any immersion $\alpha\colon X\to\mathbb{R}^{2}$, the total number of selfintersections of the loop $\alpha(\partial X)$ admits an estimate
$\displaystyle N^{+}+N^{}\geq\,2g+2.$ 
Moreover, this lower bound is realized by an immersion $\alpha\colon X\to\mathbb{R}^{2}$! Therefore, by formula (3.4), the Guth inequality is transformed into
$\displaystyle N^{+}+N^{}\geq 2+\frac{1}{2}\big(\#\{\partial_{2}^{+}X(v)\}\# \{\partial_{2}^{}X(v)\}\big).$ 
Moreover, for some optimal immersion $\alpha$,
$\displaystyle N^{+}+N^{}=2+\frac{1}{2}\big(\#\{\partial_{2}^{+}X(v)\}\#\{ \partial_{2}^{}X(v)\}\big)=2\frac{1}{2\pi}\int_{\partial X}K_{\nu}\,d\mathsf {g}.$ 
When a surface $X$ is oriented and a vector field $v$ is boundary generic, then the points from $\partial_{2}^{+}X(v)$ come in two new flavors: “$\oplus,\ominus$”. By definition, a point $a\in\partial_{2}^{+}X(v)$ has the polarity “$\oplus$” if the orientation of $T_{a}X$ determined by the pair $(\nu_{a},v(a))$, where $\nu_{a}$ is the inner normal to $\partial X$, agrees with the preferred orientation of $X$. Otherwise, the polarity of $a$ is defined to be “$\ominus$”. (These new polarities will play a crucial role in Theorem 9.2.)
Thus, for each choice of orientation of $X$ (and hence of $\partial X$) we get a partition
$\displaystyle\partial_{2}^{+}X(v)=\partial_{2}^{+,\oplus}X(v)\coprod\partial_{ 2}^{+,\ominus}X(v).$ 
Switching the orientation of $X$ switches the second polarities in the partition.
4 Convexity, Concavity, and Complexity of Flows in 2D
Definition 4.1.
We say that a boundary generic vector field $v$ is boundary convex if $\partial_{2}^{+}X(v)=\emptyset$, and is boundary concave if $\partial_{2}^{}X(v)=\emptyset$ (see Fig. 1).
The existence of a boundary convex vector field puts severe restrictions on the topology of the surface.
Lemma 4.1.
If a compact connected surface $X$ with boundary admits a boundary convex traversing vector field $v$, then $X$ is either a disk $D^{2}$, or an annulus $A^{2}$.
Proof.
The convexity of the vector field $v$ implies that $X$ admits a $(v)$directed continuous retraction on the locus $\partial_{1}^{+}X(v)$. Since $X$ is connected, it follows that $\partial_{1}^{+}X(v)$ is connected as well. Thus, $\partial_{1}^{+}X(v)$ is either a circle, or a segment. In the first case, $X$ is diffeomorphic to an annulus $A^{2}$; in the second case, $X$ is diffeomorphic to a disk $D^{2}$. $\square$
The same phenomenon occurs in any dimension: if a compact connected smooth $(n+1)$dimensional manifold $X$ with a connected boundary admits a boundary convex gradientlike nonvanishing vector field $v$, then $H_{n}(X;\mathbb{Z})=0$ ([Katz2014a]). In other words, $H_{n}(X;\mathbb{Z})\neq 0$ is a topological obstruction to the existence of a boundary convex nonvanishing gradientlike vector field on $X$.
In contrast, the boundary concave nonvanishing gradientlike vector fields are plentiful. For example, consider a radial vector field $v$ on an annulus $A^{2}$. Delete from $A^{2}$ any number of pairwise disjoint convex disks and restrict $v$ to the resulting 2disk with holes. The convexity of the removed disks implies that any disk with holes admits a boundary concave traversing vector field $v$.
Many other surfaces admit such concave vector fields as well. Here is an important example: consider a Morse function $f\colon Y\to\mathbb{R}$ on a closed surface $Y$ and its gradientlike vector field $v$. Then removing small convex (in the local Morse coordinates) disks, centered on the critical points, from $Y$, produces a boundary concave nonvanishing gradient vector field on $X$. In particular, if $Y$ is a sphere with $g$ handles, then one can find a Morse function with $2g+2$ critical points (see Fig. 2). So the orientable surface $X$, obtained from $Y$ by removing $2g+2$ balls, admits a nonvanishing boundary concave gradientlike vector field $v$.
In fact, by Theorem 6.2, any connected orientable surface with boundary, but the disk, admits a boundary concave nonvanishing gradientlike vector field!
We view the integer $c^{+}(v)=_{\mathsf{def}}\#(\partial_{2}^{+}X(v))$ as a measure of complexity of the $v$flow, subject to the condition $\mathsf{Ind}(v)=0$ or, alternatively, subject to the condition $v\neq 0$ .
We define the complexity of a compact connected surface $X$ with boundary as the minimum
$\displaystyle c^{+}(X)=\min_{\{v\neq 0\}}\{c^{+}(v)\},$ 
where $v$ runs over all nonvanishing boundary generic vector fields on $X$.
By varying $v$ within different spaces of vector fields, one may consider a variety of such minima; nonvanishing vector fields and nonvanishing gradientlike (equivalently, traversing) vector fields are the two most important cases. So we reintroduce^{7}^{7}see the paragraph that follows Theorem 1.1. the gradient complexity
$\displaystyle gc^{+}(X)=_{\mathsf{def}}\min_{\{v\neq 0\;\text{of the gradient type}\}}\{c^{+}(v)\},$ 
where $v$ runs over all nonvanishing gradientlike and boundary generic (equivalently, traversing and boundary generic) vector fields on $X$.
Evidently $gc^{+}(X)\geq c^{+}(X)$. Let $M$ denote the Möbius band. In Sect. 6, we will show that $gc^{+}(M)=1$, while $c^{+}(M)=0$. So the two notions of complexity are different!
When $c^{+}(v)\leq N$, then by Lemma 3.1, $4g+2k4\leq N$, which confines $(g,k)\in\mathbb{Z}_{+}^{2}$ to a triangle whose area is a quadratic function in the variable $N$. Hence, in terms of the complexity $c^{+}(v)$, we can restate the Lemma 3.1 as follows.
Corollary 4.1.
Let $X$ be a connected compact surface with boundary. Let $v$ be a boundary generic vector field on $X$, subject to the condition $\mathsf{Ind}(v)=0$.
Then the complexity of $v$ satisfies the inequality
$\displaystyle c^{+}(v)\geq 2\cdot\dim_{\mathbb{R}}H_{1}(X;\mathbb{R})2=2\cdot\chi(X).$ 
When $\chi(X)\leq 0$, this inequality turns into the equality $c^{+}(v)=2\cdot\chi(X)$ if and only if $v$ is boundary concave.
As a result, for any natural $N$, there are finitely many (counted up to a homeomorphism) connected compact surfaces of bounded complexity $c^{+}(X)\leq N$. In fact, the number of such surfaces is a quadratic function in $N$.
Example 4.1.
In fact, for any nonvanishing boundary concave vector field $v$ on the torus with a single hole, $\#\{\partial_{2}^{+}X(v)\}=2$. The constant vector field $v$, being restricted to the complement to a convex disk in $T^{2}$, is boundary concave and has the property $\#\{\partial_{2}^{+}X(v)\}=2$. Thus, by Corollary 4.1, $c^{+}(X)=2$.
Lemma 3.1 leads immediately to the following corollary:
Corollary 4.2.
Let $X$ be a sphere with $g$ handles and $k$ holes, where $g,k\geq 1$. If $X$ admits a nonvanishing boundary concave vector field $v$, then $c^{+}(X)=4g4+2k$.
Given a compact surface $X$ with boundary, we form its double $DX=_{\mathsf{def}}X\cup_{\partial X}X$ by attaching two copies of $X$ along their boundaries. Note that $\chi(DX)=2\cdot\chi(X)$. Therefore, $\chi(X)<0$ if and only if $\chi(DX)<0$.
Recall that any closed orientable surface with a negative Euler number admits a metric of constant negative curvature $1$, the hyperbolic metric. So if $\chi(X)<0$, then $DX$ admits such hyperbolic metric.
Let $vol(DX)$ denote the hyperbolic volume of $DX$, and let $vol(\Delta^{2})$ denote the volume of an ideal hyperbolic triangle $\Delta^{2}$ in the hyperbolic plane $\mathbf{H}^{2}$. In fact, $vol(\Delta^{2})=\pi$.
In $2D$, a remarkable convergence of topology and geometry takes place. In the spirit of this convergence, since $\chi(DX)=vol(DX)/vol(\Delta^{2})$, Corollary 4.1 admits a more geometric reformulation:
Theorem 4.1.
Let $v$ be a nonvanishing boundary generic vector field on a compact connected and orientable surface $X$ with boundary. Assume that $\chi(X)<0$^{8}^{8}This excludes the disk and the annulus.. Then the complexity of the $v$flow satisfies the inequality:
$\displaystyle c^{+}(v)\;\geq\;vol(DX)/\pi.$ 
Moreover, $c^{+}(v)=vol(DX)/\pi$ if and only if $v$ is boundary concave.
Theorem 4.1 admits far reaching multidimensional generalizations [see [Alpert and Katz2015], [Katz2016b]]. They are valid for, so called, traversally generic vector fields [see Definitions 5.1 and [Katz2014b]] on arbitrary smooth compact $(n+1)$dimensional manifolds $X$ with boundary. Such vector fields $v$ naturally generate stratifications of trajectory spaces $\mathcal{T}(v)$, whose strata are labeled by the combinatorial patterns of tangency from the universal partially ordered set $\Omega^{\bullet}_{{}^{\prime}\langle n]}$ [see the end of Sect. 6 and [Katz2016a]]. In high dimensions, we use the simplicial seminorms $\\sim\_{\Delta}$ of Gromov [Gromov1982] on the homology $H_{\ast}(X;\mathbb{R})$ and $H_{\ast}(DX;\mathbb{R})$ (as a substitute of the hyperbolic volume) to provide lower bounds on the number of connected components of the $\Omega^{\bullet}_{{}^{\prime}\langle n]}$strata of any given dimension. These estimates are sensitive to the algebraic nature of the fundamental groups $\pi_{1}(X)$, $\pi_{1}(DX)$. In particular, they are trivial for so called amenable groups [see [Gromov1982] for the definition].
5 On Spaces of Vector Fields
Recall that each trajectory $\gamma$ of a traversing vector field $v$ on a compact connected surface $X$ must reach the boundary $\partial X$ both in positive and negative times.
We denote by $\mathcal{V}_{\mathsf{trav}}(X)$ the space (in the $C^{\infty}$topology) of all traversing vector fields on a given compact surface $X$.
We denote by $\mathcal{V}_{\mathsf{grad}}(X)$ the space (in the $C^{\infty}$topology) of all gradientlike vector fields on $X$, and by $\mathcal{V}_{\neq 0}(X)$ the space of all nonvanishing vector fields on $X$.
The next lemma says that $v$ is traversal if and only if it is nonvanishing and gradientlike [(see [Katz2014a], Corollary 4.1, for the proof].
Lemma 5.1.
For any compact connected surface $X$ with boundary,
$\displaystyle\mathcal{V}_{\mathsf{trav}}(X)=\mathcal{V}_{\mathsf{grad}}(X)\cap \mathcal{V}_{\neq 0}(X),$ 
an open set in the space $\mathcal{V}_{\mathsf{grad}}(X)$.
The surfaces $X$ and vector fields $v$ we consider are all smooth. We can add an external collar to $X$ to form a diffeomorphic surface $\hat{X}\supset X$ and to extend $v$ to a smooth vector field $\hat{v}$ on $\hat{X}$. Let $\hat{\gamma}$ be a $\hat{v}$trajectory (or rather its germ) through a point $x$ of $\partial X$. We can talk about order of tangency of two smooth curves, $\hat{\gamma}$ and $\partial X$, at $x\in\hat{\gamma}\cap\partial X$ in $\hat{X}$ (see Definition 7.1). We say that the tangency of $\hat{\gamma}$ to $\partial X$ is simple if its degree is 2. When the two curves are transversal at $x$ we say that the order of tangency is 1. In fact, these notions depends only on $(X,v)$ and not on the extension $(\hat{X},\hat{v})$.
Definition 5.1.
A traversing vector field $v$ on a compact surface $X$ is called traversally generic, if two properties are valid:

(1)
when a trajectory $\gamma$ is tangent the boundary $\partial X$, then the tangency is simple,

(2)
no $v$trajectory $\gamma$ contains more then one simple point of tangency to $\partial X$.^{9}^{9}In particular, a traversally generic $v$ is boundary generic.
We denote by $\mathcal{V}^{\ddagger}(X)$ the space of all traversally generic vector fields on a compact surface $X$. In fact, the notion of traversally generic vector field is available in any dimension [see [Katz2014b]].
As the name suggests, the traversally generic vector fields are typical among all traversing vector fields; furthermore, a perturbation of any traversally generic vector field is traversally generic. This is the content of the next theorem. Its validation requires an involved argument, which even in $2D$ resists a significant simplification [Katz2014b].
Theorem 5.1.
For any compact connected surface $X$ with boundary, the space $\mathcal{V}^{\ddagger}(X)$ of traversally generic vector fields is open and dense in the space $\mathcal{V}_{\mathsf{trav}}(X)=\mathcal{V}_{\mathsf{grad}}(X)\cap\mathcal{V}_{ \neq 0}(X)$.
6 GraphTheoretical Approach to the Concavity of Traversing Vector Fields in 2D
Recall that $c^{+}(X)\leq gc^{+}(X)$. So, we start with a natural question.
Question 6.1.
Are there compact connected surfaces $X$ with boundary for which $c^{+}(X)<gc^{+}(X)$?
On many occasions we took advantage of the fact that, for traversally generic vector fields $v$, the trajectory spaces $\mathcal{T}(v)$ are finite graph whose verticies have valency 1 and 3 only (see Fig. 5). Moreover, for a traversally generic boundary concave vector field $v$, all the verticies of $\mathcal{T}(v)$ have valency 3. Now we will take a closer look at the graphtheoretical models of the boundary concave and traversally generic vector fields in 2D.
Let $G$ be a finite connected trivalent graph with $a$ verticies. We denote by $\beta G$ its barycentric subdivision: each edge $e$ of $G$ is divided by a new vertex $v_{e}$, its center. We consider the finite set $\mathsf{Tri}(G)$ of all colorings of the edges of $\beta G$ with three colors so that, at each vertex of $G$, exactly three distinct colors are applied. Thus, $\#\mathsf{Tri}(G)=6^{a}$.
Theorem 6.1.
Let $G$ be a finite connected trivalent graph. Each coloring $\alpha\in\mathsf{Tri}(G)$ produces (in a canonical way) a compact connected surface $X(G,\alpha)$ with boundary. The surface $X(G,\alpha)$ admits a traversally generic concave vector field $v(G,\alpha)$. The cardinality of the locus $\partial_{2}^{+}X(G,\alpha)\big(v(G,\alpha)\big)$ is the number of verticies in $G$.
Moreover, every connected surface with boundary, which admits a traversally generic concave vector field, can be produced in this way.
Proof.
Let $\mathsf{A},\mathsf{B},\mathsf{C}$ denote the three distinct colors, and $\mathcal{P}=\{\mathsf{A},\mathsf{B},\mathsf{C}\}$ the entire pallet.
Consider a 2dimensional space $Z=G\times(0,4)$. It has singularities in the form of binders of the 3page open books (see Fig. 6). The binders correspond to the vertices of $G$.
First, employing a given coloring $\alpha$, we will construct a piecewise linear surface $\hat{X}(G,\alpha)\subset Z$. The vector field on $\hat{X}(G,\alpha)$ will be induced by the product structure in $Z$.
For each edge $e\subset G$ and its barycenter $v_{e}\in\beta G$, we place the interval $v_{e}\times[2,3]\subset Z$ over $v_{e}$. Let $\hat{e}$ be half of the interval $e\subset G$, bounded two verticies $v_{e}\in\beta G$ and $w\in G$. Over $\hat{e}$, we place a strip $E\subset Z$; its construction depends on the color attached to the interval $[v_{e},w]$ as follows:

•
if the color of $[v_{e},w]$ is $\mathsf{A}$, then we link the vertex $v_{e}\times 2$ with the vertex $w\times 1$ by a line in the rectangle $R=[v_{e},w]\times(0,4)$, and the vertex $v_{e}\times 3$ with the vertex $w\times 2$ by another line in $R$;

•
if the color of $[v_{e},w]$ is $\mathsf{B}$, then we link by a line in $R$ the vertex $v_{e}\times 2$ with the vertex $w\times 2$, and the vertex $v_{e}\times 3$ with the vertex $w\times 3$ by another line;

•
if the color of $[v_{e},w]$ is $\mathsf{C}$, then we link by a line in $R$ the vertex $v_{e}\times 2$ with the vertex $w\times 1$, and the vertex $v_{e}\times 3$ with the vertex $w\times 3$ by another line.
By definition, $E(e,w)$ is the strip in $[v_{e},w]\times(0,4)$, bounded by the two lines whose construction has been described above. Thanks to the monotonicity of the bijections $A\colon\{2,3\}\to\{1,2,3\}$, $B\colon\{2,3\}\to\{1,2,3\}$, and $C\colon\{2,3\}\to\{1,2,3\}$ that correspond to the colors $\mathsf{A},\mathsf{B},\mathsf{C}$, the lines that bound the strips $E(e,w)$ do not intersect. We denote by $\hat{X}(G,\alpha)$ the union of all such strips.
The local model of each binder implies that $\hat{X}(G,\alpha)$ is indeed piecewiselinear surface, embedded in the singular space $Z$. Inside $Z$, one can smoothen the sharp edges of the boundary $\partial\hat{X}(G,\alpha)$ in order to get a smooth surface $X(G,\alpha)$ (we challenge the reader to visualize this smoothing in the vicinity of point $w\times 2$ from Fig. 6). The restriction of the product structure in $Z$ to its subspace $X(G,\alpha)$ produces a smooth nonvanishing vector field $v(G,\alpha)$ on $X(G,\alpha)$. Its trajectories (the vertical lines in $Z$) will be simply tangent to $\partial X(G,\alpha)$ exactly at the points of the type $w\times 2$, where $w$ runs over the set of verticies of $G$. By Theorem 5.1, this vector field $v(G,\alpha)$ is of the gradient type.
Conversely, any traversally generic and concave vector field $v$ on a connected compact surface $X$ with boundary, produces a map $\Gamma\colon X\to\mathcal{T}(v)$, where the space of trajectories is a finite trivalent graph. Its verticies are in 1to1 correspondence with the points of the locus $\partial_{2}^{+}X(v)$.
As a point $v_{e}$ in the open edge $e$ of the graph $\mathcal{T}(v)$ approaches a vertex $w$, the intersection of the $v$trajectory $\gamma=\Gamma^{1}(v_{e})$ with the boundary $\partial X$ defines a bijection of the $v$ordered set $\gamma\cap\partial X$ of cardinality 2 to a $v$ordered set $\Gamma^{1}(w)\cap\partial X$ of cardinality 3, the orders being respected by the bijections. This determines one of three colors we attach to the halfedge $[v_{e},w]$. Therefore the geometry of the flow determines a tricoloring of the graph $\beta{\mathcal{T}}(v)$. $\square$
Corollary 6.1.
The cardinality of the set $\partial_{2}^{+}X(v)$ is even for any boundary convex and traversally generic vector field $v$.
Proof.
For any boundary convex and traversally generic vector field $v$, the graph $\mathcal{T}(v)$ is trivalent. Its verticies are in 1to1 correspondence with the points of $\partial_{2}^{+}X(v)$. Since for any trivalent graph, $3\times\,\{\text{the number of verices}\}=2\times\,\{\text{the number of edges}\}$, we conclude that the number of vericies in $\mathcal{T}(v)$ must be even. $\square$
The next theorem provides complete answers to Question 1.1 and Question 1.2.
Theorem 6.2.

•
For any compact connected surface $X$ with boundary, but the disk and the Möbius band, the two complexities are equal and are given by the formula:
$\displaystyle gc^{+}(X)=c^{+}(X)=2\beta_{1}(X)2,$ where $\beta_{1}(X)$ denotes the first Betti number of $X$. Moreover, any such surface $X$, but the disk and the Möbius band, admits a boundary concave traversing^{10}^{10}Even traversally generic vector field.

•
In contrast, the disk $D$ does not admit a boundary concave traversing vector field. In fact, $c^{+}(D)=gc^{+}(D)=0$.

•
The Möbius band $M$ does not admit a boundary concave traversing vector field either. However, $M$ does admit a concave and convex nonvanishing vector field. In fact, $c^{+}(M)=0$, and $gc^{+}(M)=1$.

•
The annulus is the only compact connected surface that admits both concave and convex traversing vector fields.
Proof.
Recall the topological classification of closed surfaces. Any connected orientable closed surface is either a sphere $S$, or a torus $T$, or a connected sum of several tori $T$. Similarly, any connected nonorientable closed surface is a projective plane $P$, or a connected sum of several projective planes $P$. Evidently, any compact connected orientable surface with boundary can be obtained from a closed connected surface by deleting one or several open disks. Therefore, any compact connected orientable surface with boundary is either a disk $D$, or an annulus $A$, or a punctured torus $Q=T{\setminus}D^{\circ}$, or a boundary connected sum of such building blocks. Similarly, any compact connected nonorientable surface with boundary is either a Möbius band $M$, or a boundary connected sum of a several Möbius bands, or a connected sum of these surfaces with one or several annuli $A$.
By Corollary 4.1, for any compact connected $X$ with boundary,
$\displaystyle gc^{+}(X)\geq c^{+}(X)\geq 2\beta_{1}(X)2,$ 
where $\beta_{1}(X)$ denotes the first Betti number of $X$. If $X$ admits a boundary concave nonvanishing vector field $v$, then by the same corollary, $c^{+}(X)=c^{+}(v)=2\beta_{1}(X)2$. In particular, if $X$ admits a boundary concave traversing field, then
$\displaystyle gc^{+}(X)=c^{+}(X)=2\beta_{1}(X)2.$ 
Note that the first Betti numbers are additive under the connected sum operation:
$\displaystyle\beta_{1}(X_{1}\#_{\partial}X_{2})=\beta_{1}(X_{1})+\beta_{1}(X_{ 2}).$ 
The additivity of the first Betti numbers implies the relation
$\displaystyle 2\beta_{1}(X_{1}\#_{\partial}X_{2})2=[2\beta_{1}(X_{1})2]+[2 \beta_{1}(X_{2})2]+2.$  (6.1) 
Here is the beginning of a list of compact connected surfaces with boundary, arranged according the increasing values of their first Betti numbers:
$\displaystyle\beta_{1}=0;$  $\displaystyle=1;$  $\displaystyle\quad\;=2;\qquad\qquad\cdots$  
$\displaystyle\;D;$  $\displaystyle A,\,M;$  $\displaystyle\;Q,\;A\#_{\partial}A,\;A\#_{\partial}M,\;M\#_{\partial}M;\quad\cdots$  (6.2) 
Using the additivity of the first Betti numbers under the connected sum operation, for any compact $X$ with boundary, we get:
$\displaystyle\beta_{1}(A\#_{\partial}X)=1+\beta_{1}(X),\;\;\beta_{1}(M\#_{ \partial}X)=1+\beta_{1}(X),\;\;\beta_{1}(Q\#_{\partial}X)=2+\beta_{1}(X).$ 
Now, let us examine the behavior of the two complexities, $c^{+}(\sim)$ and $gc^{+}(\sim)$, under the connected sum operation $\#_{\partial}$.
Given two boundary generic vector fields $v_{1}$ and $v_{2}$ on $X_{1}$ and $X_{2}$, respectively, there exists a boundary generic vector field $w$ on $X_{1}\#_{\partial}X_{2}$ such that the cardinalities of the concave and convex loci satisfy the equations ^{11}^{11}Note that Eq. (6.3) resembles Eq. (6.1).:
$\displaystyle\partial_{2}^{+}(X_{1}\#_{\partial}X_{2})(w)=\partial_{2}^{+}( X_{1})(v_{1})+\partial_{2}^{+}(X_{2})(v_{2})+2,$  (6.3)  
$\displaystyle\partial_{2}^{}(X_{1}\#_{\partial}X_{2})(w)=\partial_{2}^{}( X_{1})(v_{1})+\partial_{2}^{}(X_{2})(v_{2}).$ 
Indeed, we may attach a 1handle $H$ to the locus
$\displaystyle\partial_{1}^{}X_{1}(v_{1})\coprod\partial_{1}^{+}X_{2}(v_{2})$ 
so that $H$ has a neck with respect to the extension $w$. Such vector field $w$ contributes two new points to the locus $\partial_{2}^{+}(X_{1}\#_{\partial}X_{2})(w)$ and no new points to $\partial_{2}^{}(X_{1}\#_{\partial}X_{2})(w)$. Of course, this construction fails when the set $\partial_{1}^{}X_{1}(v_{1})$ or $\partial_{1}^{+}X_{2}(v_{2})$ is empty; however, one may change $v_{1}$ to $v_{1}$, or $v_{2}$ to $v_{2}$ to arrange for $\partial_{1}^{}X_{1}(v_{1})\neq\emptyset$ and/or $\partial_{1}^{+}X_{2}(v_{2})\neq\emptyset$. These flips of the vector fields do no affect the loci $\partial_{2}^{\pm}X_{1}$ and $\partial_{2}^{\pm}X_{2}$.
In fact, this construction of the vector field $w$ validates the inequalities
$\displaystyle c^{+}(X_{1}\#_{\partial}X_{2})\leq c^{+}(X_{1})+c^{+}(X_{2})+2, \quad{\!\!}gc^{+}(X_{1}\#_{\partial}X_{2})\leq gc^{+}(X_{1})+gc^{+}(X_{2})+2.$ 
Note that when both vector fields $v_{1}$ and $v_{2}$ are boundary concave, then $w$ is boundary concave as well. As a result, the set $\mathcal{C}$ of compact connected surfaces^{12}^{12}Considered up to diffeomorphisms that admit nonvanishing boundary concave vector fields is closed under the boundary connected sum operation $\#_{\partial}$ ! Similarly, the set $\mathcal{CT}$ of compact connected surfaces that admit boundary concave and traversing vector fields is closed under $\#_{\partial}$. We will see soon that $\mathcal{C}$ consists of all compact connected surfaces with boundary, with a single exception — the disk $D$. At the same time, $\mathcal{CT}$ consists of all compact connected surfaces with boundary, except for $D$ and $M$.
Let $v_{1}$ and $v_{2}$ be some nonvanishing boundary generic (alternatively, traversing and boundary generic) vector fields which deliver the two complexities $c^{+}(X_{1})$ and $c^{+}(X_{2})$ (alternatively, $gc^{+}(X_{1})$ and $gc^{+}(X_{2})$). The previous construction of $w$ (extending $v_{1}$ and $v_{2}$ across the handle $H$) has an implication: if $c^{+}(X_{1})=2\beta(X_{1})2$ and $c^{+}(X_{2})=2\beta(X_{2})2$, then
$\displaystyle c^{+}(X_{1}\#_{\partial}X_{2})\leq c^{+}(X_{1})+c^{+}(X_{2})+2=2 \beta(X_{1}\#_{\partial}X_{2})2.$ 
Similarly, if $gc^{+}(X_{1})=2\beta(X_{1})2$ and $gc^{+}(X_{2})=2\beta(X_{2})2$, then
$\displaystyle gc^{+}(X_{1}\#_{\partial}X_{2})\leq gc^{+}(X_{1})+gc^{+}(X_{2})+ 2=2\beta(X_{1}\#_{\partial}X_{2})2.$ 
Since the reverse inequality holds by Corollary 4.1, we get
$\displaystyle c^{+}(X_{1}\#_{\partial}X_{2})=2\beta(X_{1}\#_{\partial}X_{2})2,$  (6.4) 
provided that $c^{+}(X_{1})=2\beta(X_{1})2$ and $c^{+}(X_{2})=2\beta(X_{2})2$.
Similarly, by Corollary 4.1, we get
$\displaystyle gc^{+}(X_{1}\#_{\partial}X_{2})=2\beta(X_{1}\#_{\partial}X_{2})2$  (6.5) 
when $gc^{+}(X_{1})=2\beta(X_{1})2$ and $gc^{+}(X_{2})=2\beta(X_{2})2$. So we conclude that, if $v_{1}$ delivers the complexity/gradient complexity of $X_{1}$, and $v_{2}$ delivers the complexity/gradient complexity of $X_{2}$, then $w$ delivers the the complexity/gradient complexity of $X_{1}\#_{\partial}X_{2}$.
Let us now compute the complexities of the surfaces $D,A,Q$, the basic building blocks in the representation (as boundary connected sums) of orientable connected surfaces with boundary, and the complexities of the basic blocks $M,A$, which participate in the representation of nonorientable connected surfaces with boundary.
Note that $c^{+}(D)=0=gc^{+}(D)$, since $D$ admits a convex traversing flow. By formula (3.2), $D$ does not admit a boundary concave nonvanishing vector field.
Also $c^{+}(A)=0=gc^{+}(A)$, the latter equality being delivered by the radial gradientlike vector field.
We claim that $c^{+}(Q)=2=gc^{+}(Q)$. Indeed, since $\beta_{1}(Q)=2$, by Corollary 4.1, we get $c^{+}(Q)\geq 2$. On the other hand, there exists a trivalent graph $G_{Q}$ with an appropriate tricoloring and exactly two verticies such that, applying the construction from Theorem 6.1, we produce a traversally generic vector field $v(G_{Q},\alpha)$ on the surface $X(G_{Q},\alpha)=Q$ with the cardinality 2 locus $\partial_{2}^{+}X(G_{Q},\alpha)(v(G_{Q},\alpha))$. As a result, both complexities of $Q$ equal to 2.
The Möbius band $M$ is different. We notice that $M$ admits a nonvanishing vector field $v$ with a single closed trajectory—the core of the Möbius band—and transversal to the boundary $\partial M$. Thus, $c^{+}(M)=0$.
Now consider a trivalent graph $G_{M}$ (“$\mathbf{*\bullet\,\big<\big>}$”) with a single vertex “ $\bullet$” of valency 3 and a single vertex “$*$” of valency 1 (this $G_{M}$ is a circle “$\mathbf{\big<\big>}$” to which a radius “$\mathbf{*\bullet}$” is attached). The construction from Theorem 6.1 applies to produce a remarkable embedding of the Möbius band in the product $G_{M}\times[0,4]$. So we conclude that $M$ admits a traversally generic vector field $v_{M}$ (not a concave one!) with $\partial_{2}^{+}M(v_{M})$ being a singleton ($\partial_{2}^{}M(v_{M})$ is a singleton as well). As a result, $gc^{+}(M)\leq 1$. On the other hand, any traversally generic vector field $v$ on $M$ must produce the graph $\mathcal{T}(v)$ which is homotopy equivalent to a circle, the homotopy type of $M$. If $\partial_{2}^{+}M(v)=\emptyset$, this graph $\mathcal{T}(v)$ has no trivalent verticies, in which case, $\mathcal{T}(v)$ is homeomorphic to a circle. So $M\to\mathcal{T}(v)$ must be a fibration whose fibers (the $v$trajectories) are segments. Moreover, thanks to the vector field $v$, this fibration is orientable, a contradiction with the nonorientability of $M$. Therefore, we conclude that $gc^{+}(M)=1$, while $c^{+}(M)=0$.
Let $K=_{\mathsf{def}}M\#_{\partial}M$, $L=_{\mathsf{def}}A\#_{\partial}M$, and $N=_{\mathsf{def}}A\#_{\partial}A$. In fact, Theorem 6.1 implies that $Q$, $K$, $L$, and $N$ are the only connected surfaces of the gradient complexity 2; all of them admit boundary concave traversally generic vector fields. Indeed, just start with the tree “$\big>\bullet\bullet\big<$” with two trivalent verticies and consider the ways in which one can identify its four leaves in pairs. Then consider all admissible tricologings of the resulting graphs $G$. These cases will deliver the model tricolored graphs $G_{Q},G_{K},G_{L},G_{N}$.
Therefore $gc^{+}(X)=c^{+}(X)=2$ only for $X=Q,K,L,N$.
Now the rest of the surfaces (of complexities higher than 2) can be obtained by applying repeatedly the connected sum operations $X\Rightarrow X\#_{\partial}T$, $X\Rightarrow X\#_{\partial}M$, and $X\Rightarrow X\#_{\partial}A$ to the surfaces $X=Q,K,L,N$ of complexity/gradient complexity 2.
When $c^{+}(X)=2\beta_{1}(X)2$ or $gc^{+}(X)=2\beta_{1}(X)2$, these changes are described by formulas (6.4) and (6.5), respecively.
First, let us consider the operation $X\Rightarrow X\#_{\partial}M$. According to these formulas, the inequalities
$\displaystyle 2\beta_{1}(X\#_{\partial}M)2\leq c^{+}(X\#_{\partial}M)\leq gc^ {+}(X\#_{\partial}M)\leq gc^{+}(X)+2$ 
can be reduced to
$\displaystyle 2\beta_{1}(X)\leq c^{+}(X\#_{\partial}M)\leq gc^{+}(X\#_{ \partial}M)\leq 2\beta_{1}(X),$ 
provided that $gc^{+}(X)=2\beta_{1}(X)2$. For such $X$, we get again the property
$\displaystyle c^{+}(X\#_{\partial}M)=gc^{+}(X\#_{\partial}M)=2\beta_{1}(X)=2 \beta_{1}(X\#_{\partial}M)2.$ 
Similar considerations hold for the operations $X\Rightarrow X\#_{\partial}T$ and $X\Rightarrow X\#_{\partial}A$.
Therefore, by an induction in complexities, for all the surfaces $X$, but $D$ and $M$, $c^{+}(X)=gc^{+}(X)=2\beta_{1}(X)2$. Moreover, all such surfaces admit boundary concave and traversing vector fields. $\square$
7 Combinatorics of Tangency for Traversing flows in 2D
Pick an extension $\hat{X}$ of a given compact surface $X$ by adding an external collar to $X$. Let $\hat{v}$ be an extension of a given vector field $v$ into $\hat{X}$. Pick a smooth auxiliary function $z\colon\hat{X}\to\mathbb{R}$ such that:

•
0 is a regular value of z,

•
$z^{1}(0)=\partial X$,

•
$z^{1}((\infty,0])=X$.
(7.1)
Definition 7.1.
Let $\hat{\gamma}$ be a $\hat{v}$trajectory through a point $x\in\partial X$. We say that $\hat{\gamma}$ has the order of tangency (multiplicity) $k$ to $\partial X$ at $x$, if $\mathcal{L}_{\hat{v}}^{\{j\}}(z)=0$ for all $j<k$, and $\mathcal{L}_{\hat{v}}^{\{k\}}(z)\neq 0$ at $x$ ^{13}^{13}This is equivalent to saying that the $(k1)$st jet at $x$ of $z_{\gamma}$ vanishes, but the $k$th jet does not.. Here $\mathcal{L}_{\hat{v}}^{\{j\}}(z)$ denotes the $j^{th}$ iterated $\hat{v}$directional derivative of the function $z$.
Given a traversally generic (see Definition 5.1) vector field $v$ on a compact connected surface $X$, we will attach the combinatorial pattern $(11)$ to typical $v$trajectory $\gamma\subset X$ that correspond to points in the interior of an edge of the graph $\mathcal{T}(v)$, the pattern $(121)$ to the trajectories that correspond to the trivalent verticies of $\mathcal{T}(v)$, and the pattern $(2)$ to the univalent verticies (see Fig. 5). In fact, the numbers 1 and 2 in these patterns reflect the order of tangency of the curves $\hat{\gamma}$ and $\partial X$ at the points of $\gamma\cap\partial X$ (see Definition 7.1).
On a given compact surface $X$, for traversally generic vector fields $v$ no other patterns (say, like $(1221)$ or $(13)$) occur. In $2D$, this conclusion follows from Definition 5.1.
The lemma below is another way to state this fact. Its proof, relying on the Malgrange Preparation Theorem [Malgrange1964], can be found in [Katz2014b], Lemma 3.4.
Lemma 7.1.
Let $v$ be a traversally generic vector field on $X$. Extend $(X,v)$ to a pair $(\hat{X},\hat{v})$. In the vicinity of each $v$trajectory $\gamma$, there exist special local coordinates $(u,x)$ in $\hat{X}$ and a real polynomial $P(u,x)$ of degree 2 or 4 such that:

•
each $\hat{v}$trajectory is given by the equation $\{x=const\}$,

•
the boundary $\partial X$ is given by the polynomial equation $\{P(u,x)=0\}$,

•
$X$ is given by the polynomial inequality $\{P(u,x)\leq 0\}$.
The polynomial $P(u,x)$ takes three canonical forms:

(1)
$u(u1)$, which corresponds to the combinatorial pattern $\mathbf{(11)}$,

(2)
$u^{2}x,$ which corresponds to the combinatorial pattern $\mathbf{(2)}$,

(3)
$u\big((u1)^{2}+x\big)(u2),$ which corresponds to the pattern $\mathbf{(121)}$.
To summarize, at $\partial_{2}X(v)$ the order of tangency is 2; the trajectories through $\partial_{2}^{+}X(v)$ have the combinatorial tangency pattern $(121)$, and through $\partial_{2}^{}X(v)$ the combinatorial tangency pattern $(2)$. The rest of trajectories have the pattern $(11)$.
We denote by $\Omega^{\bullet}_{{}^{\prime}\langle 1]}$ the partially ordered set whose elements are $(11),(2),(121)$ and the order is defined by $(11)\succ(2)$ and $(11)\succ(121)$. This combinatorics does not look impressive. However, in higher dimensions, traversally generic vector fields on $(n+1)$manifolds with boundary generate a rich partially ordered finite list $\Omega^{\bullet}_{{}^{\prime}\langle n]}$ of combinatorial tangency patters. The poset $\Omega^{\bullet}_{{}^{\prime}\langle n]}$ is universal in each dimension $n+1$. It is discussed in [Katz2016a].
8 Holography of Traversing Flows on Surfaces
Let $v$ be a traversally generic (see Definition 5.1) vector field on a compact connected surface $X$ with boundary. In Sect. 1, we have introduced the causality map
$\displaystyle C_{v}:\partial_{1}^{+}X(v)\to\partial_{1}^{}X(v).$ 
It is a distant relative of the classical Poincaré Return Map [e.g., see [Teschi2012]]^{14}^{14}In the forthcoming paper we will see that, for the geodesic flows $v^{g}$ on the tangent spherical bundle $SM\to M$ over a smooth connected Riemannian $n$manifold $(M,g)$ with boundary, $C_{v^{g}}$ is the scattering map, a version of the billiard map..
Alternatively, one can think of $C_{v}$ as determining a partial order “$x\prec C_{v}(x)$” among the points $x$ of the boundary $\partial X$.
The word “causality” in the name of $C_{v}$ is motivated by the following pivotal special case.
Example 8.1.
Let $w=w(\theta,t)$ be a smooth timedependent vector field on the circle $S^{1}$ (equipped with the angular coordinate $\theta$). It gives rise to a vector field $v=(w,1)$ on the cylinder $S^{1}\times\mathbb{R}$. We think about the factor $S^{1}$ as space and about the factor $\mathbb{R}$ as time. We denote by $t$ the time coordinate. So we call $S^{1}\times\mathbb{R}$ the space of events. Note that $v$ is a nonvanishing gradientlike vector field with respect to the time function $T\colon S^{1}\times\mathbb{R}\to\mathbb{R}$.
Pick any smooth compact surface $X\subset S^{1}\times\mathbb{R}$ with a boundary $\partial X$. We call $X$ the event domain, and its boundary $\partial X$ the event horizon.
Since the nonvanishing vector field $v$ is traversing in $X$, the map $C_{v}$ is welldefined. Then the map $C_{v}\colon\partial_{1}^{+}X(v)\to\partial_{1}^{}X(v)$ indeed gives rise the causality relation on the event horizon: the correspondence $C_{v}$ reflects the evolution of an event $z=(\theta,t)$ into the event $C_{v}(z)$.
Let $\mathcal{C}(\partial_{2}^{+}X(v))$ denotes the union of $v$trajectories through the points of the concavity locus $\partial_{2}^{+}X(v)$. The causality map is discontinuous at the points of the intersection $\mathcal{C}(\partial_{2}^{+}X(v))\cap\partial_{1}^{+}X(v)$ (see Fig. 7). On the positive side, the discontinuities of the causality map $C_{v}$ are not too bad: the map has the “left” and “right” limits.
Given a pair $(X,v)$, the $v$trajectories, viewed as unparametrized $v$oriented curves, produce an oriented 1dimensional foliation $\mathcal{F}(v)$ on $X$.
Theorem 8.1.
(The Causal Holography Principle in 2D)
Let $(X_{1},v_{1})$ and $(X_{2},v_{2})$ be two compact connected surfaces with boundaries, carrying traversally generic vector fields $v_{1}$ and $v_{2}$, respectively. Assume that there is a diffeomorphism $\Phi^{\partial}\colon\partial X_{1}\to\partial X_{2}$ which conjugates the two causality maps:
$\displaystyle C_{v_{2}}\circ\Phi^{\partial}=\Phi^{\partial}\circ C_{v_{1}}.$ 
Then $\Phi^{\partial}$ extends to a diffeomorphism $\Phi\colon X_{1}\to X_{2}$ which maps the oriented foliation $\mathcal{F}(v_{1})$ to the oriented foliation $\mathcal{F}(v_{2})$.
Proof.
We will only sketch the argument. A fully developed proof of the multidimensional analogue of this theorem is contained in [Katz2014c].
First, we notice that since $C_{v_{1}}$ and $C_{v_{2}}$ are $\Phi^{\partial}$conjugate, the diffeomorphism $\Phi^{\partial}$ induces a welldefined continuous map $\Phi_{\mathcal{T}}\colon\mathcal{T}(v_{1})\to\mathcal{T}(v_{2})$ of the trajectory spaces. Moreover, $\Phi_{\mathcal{T}}$ preserves the stratifications of the two trajectory spaces/graphs by the combinatorial type of trajectories. That is, the trivalent verticies of $\mathcal{T}(v_{1})$ are mapped to the trivalent verticies of $\mathcal{T}(v_{2})$, the univalent verticies are mapped to univalent verticies, and the interior of the edges to the interior of the edges.
Then we pick a smooth function $f_{2}\colon X_{2}\to\mathbb{R}$ such that $df_{2}(v_{2})>0$. With the help of $\Phi^{\partial}$, we pullback $f_{2}_{\partial X_{2}}$ to get a smooth function $f_{1}^{\partial}\colon\partial X_{1}\to\mathbb{R}$ such that
$\displaystyle f_{1}^{\partial}(z)<f_{1}^{\partial}(C_{v_{1}}(z))$ 
for all $z\in\partial_{1}^{+}X_{1}(v_{1})$.
Then we argue that $f_{1}^{\partial}$ extends to a smooth function $f_{1}\colon X_{1}\to\mathbb{R}$ such that $df_{1}(v_{1})>0$.
We use $f_{1}$ to embed $X_{1}$ in the product $\mathcal{T}(v_{1})\times\mathbb{R}$ by the formula
$\displaystyle\alpha_{\{v_{1},\,f_{1}\}}(z)=(\gamma_{z},f_{1}(z)),$ 
where $\gamma_{z}$, the $v_{1}$trajectory through $z$, is viewed as the point $\Gamma_{1}(z)$ of the graph $\mathcal{T}(v_{1})$. Similarly, we use $f_{2}$ to embed the surface $X_{2}$ in the product $\mathcal{T}(v_{2})\times\mathbb{R}$ with the help of the map $\alpha_{\{v_{2},f_{2}\}}$.
Finally, we employ $\Phi_{\mathcal{T}}$, $f_{1}$ and $f_{2}$ to construct a map
$\displaystyle\hat{\Phi}:\mathcal{T}(v_{1})\times\mathbb{R}\to\mathcal{T}(v_{2})\times\mathbb{R}$ 
by the formula
$\displaystyle\hat{\Phi}(\gamma,t)=\big(\Phi_{\mathcal{T}}(\gamma),\,f_{1}(f_{2 }^{1}(t))\big),$ 
where $t$ belongs to the $f_{2}$image of the trajectory $\Gamma_{2}^{1}\big(\Phi_{\mathcal{T}}(\gamma)\big)$.
Crudely, the restriction of $\hat{\Phi}$ to $\alpha_{\{v_{1},\,f_{1}\}}(X_{1})\subset\mathcal{T}(v_{1})\times\mathbb{R}$ is the desired diffeomorphism $\Phi\colon X_{1}\to X_{2}$.
Note that, in general, the image $\Phi_{\ast}(v_{1})$ is not $v_{2}$; so the parametrizations of the trajectories are not respected by the diffeomorphism $\Phi$, but the oriented 1foliations $\mathcal{F}(v_{1})$ and $\mathcal{F}(v_{2})$ are. $\square$
Theorem 8.1 leads instantly to the following corollary, stated as Theorem 1.3 in the introduction:
Corollary 8.1.
Let $X$ be a compact connected surface with boundary, and $v$ a smooth traversally generic vector field on it.
Then the knowledge of the causality map $C_{v}\colon\partial_{1}^{+}X(v)\to\partial_{1}^{}X(v)$ is sufficient for a reconstruction of the pair $(X,\mathcal{F}(v))$, up to a diffeomorphism of $X$ that is the identity on $\partial X$.
The world “holography” is present in the name of Theorem 8.1 since the surface $X$ and the $2D$dynamical system, produced by the $v$flow, are recorded on two 1dimensional screens, $\partial_{1}^{+}X(v)$ and $\partial_{1}^{}X(v)$.
Theorem 8.1 and Corollary 8.1 are valid in any dimension ([Katz2014c]).
Example 8.2.
Let $v$ be a traversally generic vector field on a connected surface $X$ whose boundary $\partial X$ is a single loop. Then the boundary $\partial X$ is divided into $q$ disjoint arcs $a_{1},\dots,a_{q}$ that form $\partial_{1}^{+}X(v)$ and $q$ complementary arcs $b_{1},\ldots,b_{q}$ that form $\partial_{1}^{}X(v)$. The causality map
$\displaystyle C_{v}:\coprod_{i=1}^{q}a_{i}\;\to\;\coprod_{i=1}^{q}b_{i}$ 
can be represented by its graph $G(C_{v})\subset\prod_{i,j}a_{i}\times b_{j}$.
The map $C_{v}$ (the curve $G(C_{v})$) is discontinuous at exactly $c^{+}(v)$ points in $\coprod_{i=1}^{q}a_{i}$ that correspond to the points of the intersection $\mathcal{C}(\partial_{2}^{+}X(v))\cap\partial_{1}^{+}X(v)$. There the map $C_{v}$ has distinct left and right limits.
According to the Corollary 8.1, the curve $G(C_{v})\subset\prod_{i,j}a_{i}\times b_{j}$ determines $X$ and the unparametrized dynamic of the $v$flow, up to a diffeomorphism $\Phi\colon X\to X$ that is the identity on $\partial X$. Note that the number $q$ alone is not sufficient even to determine the genus of the surface $X$.
Revisiting Example 8.2, we get the following interpretation of Corollary 8.1:
Corollary 8.2.
For any smooth timedependent vector field $w$ on the circle $S^{1}$, the causality relation on the event horizon $\partial X$ is sufficient for a reconstruction of the event domain $X$ and the unparametrized $(w,1)$flow, up to a diffeomorphism of $X$ that is the identity on $\partial X$.
In order to interpret the causality map as the Poincaré return map [see [Teschi2012]] of a dynamic system, we need to leave the Flatland and to sail into 3dimensional world.
Let us recall briefly the notion of a geodesic flow. Consider a compact surface/manifold $X$ with a Riemannian metric $g$. We denote by $SX$ the space of unit tangent vectors on $X$. Given a point $x\in X$ and a unit tangent vector $v\in T_{x}X$, take the germ of the unique geodesic curve $\gamma\colon(\epsilon,\epsilon)\to X$ through $x$ such that $\gamma(0)=x$ and $\dot{\gamma}(0)=v$. Then the pair $(\gamma(t),\dot{\gamma}(t))$ is a path in $SX$; its tangent vector $v^{g}=(\dot{\gamma},\ddot{\gamma})$ belongs to the tangent space $T_{(x,v)}(SX)$. The vector field $v^{g}$ on $SX$ is called the geodesic vector field, and the flow on $SX$ it generates is called the geodesic flow.
The theory of billiards on connected Riemmanian surfaces $(X,g)$ with boundary benefits from applying a 3dimensional version of the Causal Holography Theorem 8.1 to the geodesic flow $v^{g}$ on the 3dimensional manifold $SX$. See [Katz2014c] for some of these applications. In addition to geodesic billiards, they include the classic inverse geodesic scattering problems.
Let us outline how these applications work. When $\partial X$ is connected (i.e., is a circle), then $\partial(SX)$ is a 2torus. If the metric $g$ on $X$ is such that the geodesic flow $v^{g}$ is gradientlike^{15}^{15}Such metrics form an open set in the space of all Riemannian metrics., then the causality (scattering) map $C_{v^{g}}$ acts from the annulus $\partial_{1}^{+}SX(v^{g})$ to the annulus $\partial_{1}^{}SX(v^{g})$. Both annuli can be identified via a diffeomorphism $\tau\colon\partial_{1}^{}SX(v^{g})\to\partial_{1}^{+}SX(v^{g})$ that reflects each tangent vector to $X$ at a point of $\partial X$ with respect to the boundary $\partial X$. The composition $B_{v^{g}}=_{\mathsf{def}}\tau\circ C_{v^{g}}$, called the billiard map, is a selfmap of the annulus $\partial_{1}^{+}SX(v^{g})$. The iterations of the billiard map $B_{v^{g}}$ give rise to an interesting dynamics on the annulus, a focal point of many investigations.
In this context, the main open and challenging for us question is: For which Riemannian metrics $g$ on $X$, the knowledge of the scattering map $C_{v^{g}}\colon\partial_{1}^{+}SX(v^{g})\to\partial_{1}^{}SX(v^{g})$ (equivalently, of the billiard map) is sufficient for a reconstruction of $g$, say, up to a scalar?
9 Convex Envelops and Characteristic Classes of Traversing Flows
Let $X$ be a compact surface/manifold with boundary and $v$ a traversing vector field (see Definition 1.1). We have seen that, as a function of a point $x\in X$, the $v$trajectory $\gamma_{x}\subset X$ through $x$ exhibits a discontinuous behavior in the vicinity of any concavity point $x_{0}\in\partial_{2}^{+}X$. In order to get around this fundamental difficulty, we “envelop” the pair $(X,v)$ into a pair $(\hat{X},\hat{v})$, an ambient compact surface/manifold $\hat{X}\supset X$ ($\dim(\hat{X})=\dim(X)$) and a traversing vector field $\hat{v}$, such that (see also Definition 9.3):

(1)
$\partial\hat{X}$ is convex with respect to the $\hat{v}$flow on $\hat{X}$,

(2)
$\hat{v}_{X}=v$,

(3)
the $\hat{v}$flow is traversally generic with respect to $\partial X$.
Not any pair $(X,v)$ admits such convex envelop (see Lemma 9.1). However, when available, the convex envelop $(\hat{X},\hat{v})$ simplifies the analysis of $(X,v)$ greatly.
In this context, our ultimate goal is to study the pseudoisotopies and bordisms of convex envelops, by organizing them into monoids and groups, respectively, and to compute these algebraic structures for an a priori prescribed set of combinatorial tangency patterns. For $n$dimensional flows, this goal will be achieved in the forthcoming paper ([Katz2016c]) and in a chapter from the book ([Katz2017]). Although in dimension two these structures are quite primitive, they are not completely trivial either.
Recall that, in the study of manifolds, the universal classifying spaces like Grassmanians play a pivotal role. In the category of convex envelops, the role of universal objects (of “the new Grassmanians”) is played by various spaces of smooth functions $f\colon\mathbb{R}\to\mathbb{R}$, whose zeros and their multiplicities are modeled after the real divisors of real polynomials. The topology of these functional spaces with constrained zero divisors is interesting on its own right. One kind of these functional spaces has been studied in depth by [Arnold1989] and [Vasiliev1994].
In [Katz2016c], we will compute the homology of these functional spaces in terms of quite rich and universal combinatorics. This is reminiscent to the role played by the Schubert calculus in depicting the characteristic classes of classical Grassmanians.
In somewhat similar spirit, we apply our understanding of the algebraic topology of the universal functional spaces with restricted zero divisors on the number line to the traversing flows on manifolds with boundary. In particular, in [Katz2016c], we will introduce a rich variety of characteristic classes of convex envelops. In turn, this will allow us to manufacture, in a systematic way, examples of traversing flows with certain a priori prescribed combinatorial tangency patters in relation to the boundary.
Recall that, for boundary generic 2dimensional traversing flows, no tangencies of orders $\geq 3$ occur. In light of what has been outlined above, we should anticipate a link between traversally generic flows on surfaces and the spaces of smooth functions $f\colon\mathbb{R}\to\mathbb{R}$ (or even polynomials) that have no zeros of multiplicities $\geq 3$.
Let $\mathcal{F}$ denote the space (in the $C^{\infty}$topology) of smooths functions $f\colon\mathbb{R}\to\mathbb{R}$ that are identically 1 outside of a compact interval (the interval depends on a particular function). Let $\mathcal{F}_{<k}$ be its subspace, formed by functions that have zeros only of the multiplicities less than $k$. In [Arnold1989], Arnold calls such functions “functions with moderate singularities”. One can prove that any function with moderate singularities can be represented as a polynomial times a smooth positive function [Katz2016c]. In particular, the zero set $f^{1}(0)$ is finite for each $f\in\mathcal{F}_{<k}$.
An important theorem of Vassiliev [see [Vasiliev1994], Corollary on page 81 and The First and Second Main Theorems on pp. 78–79] describes the weak homotopy types of the spaces $\mathcal{F}_{<k}$ for all $k\geq 4$ and their homology types for all $k\geq 3$. In particular, the homology of the space $\mathcal{F}_{\leq 2}=_{\mathsf{def}}\mathcal{F}_{<3}$ is isomorphic to the homology of $\Omega S^{2}$, the space of loops on a 2sphere! On the other hand, Arnold proved that the fundamental group $\pi_{1}(\mathcal{F}_{\leq 2})\approx\mathbb{Z}$ ([Arnold1989]). Thus, $H^{1}(\mathcal{F}_{\leq 2};\mathbb{Z})\approx\mathbb{Z}$. For the 2dimensional traversing flows, this fact allows to define characteristic classes $\theta_{\alpha}$ of convex envelops $\alpha\colon(X,v)\subset(\hat{X},\hat{v})$, which reside in $H^{1}(\hat{X};\mathbb{Z})$.
Let $\hat{v}$ be a boundary convex traversing vector field on an annulus $A$. With the help of $\hat{v}$, we can introduce a product structure $A\approx S^{1}\times[0,1]$ so that the fibers of the projection $\theta\colon A\to S^{1}$ are the $\hat{v}$trajectories.
Definition 9.1.
Consider a collection $L$ of several smooth immersed loops in the annulus $A$, which intersect and selfintersect transversally and do not have triple intersections.
We say that a boundary convex traversing vector field $\hat{v}$ on $A$ is generic relative to $L$, if no $\hat{v}$trajectory $\gamma$ contains more than one point of selfintersection from $L$ and no more than one point of simple tangency to $L$, but not both.
For a given $L$, by standard techniques of the singularity theory, we can find a perturbation of $\hat{v}$ within the space $\mathcal{V}_{\mathsf{trav}}(A)$ so that the perturbed vector field is generic with respect to $L$ and still convex with respect to $\partial A$ ([Katz2014b]).
Since an immersion is a smooth map of manifolds whose differential has the trivial kernel, the immersions allow for a transfer of a given vector field on the target manifold to a vector field on the source manifold. The transfer of a nonvanishing vector field is a nonvanishing vector field.
Note that any orientable surface $X$ admits an immersion $\alpha:X\to A$ (or even in the plane $\mathbb{R}^{2}$) (see Fig. 8). We will use this fact to pullback to $X$ traversing convex vector fields on the target space $A$.
Definition 9.2.
Consider an immersion $\alpha:X\to A$ of a given compact orientable surface $X$ into an annulus $A$, equipped with a traversal boundary convex (“radial”) vector field $\hat{v}$. We call such $\alpha$ generic relative to $\hat{v}$, if $\hat{v}$ is generic relative to the curves $\alpha(\partial X)$ in the sense of Definition 9.1.
Given a transversally generic vector field $v$ on a connected compact surface $X$, we call a map $\alpha\colon(X,v)\to(A,\hat{v})$ a convex quasienvelop of $(X,v)$ if there exists an immersion $\alpha\colon X\to A$ which is generic relative to the radial vector field $\hat{v}$ on $A$, and $v=\alpha^{\ast}(\hat{v})$, the pullback of $\hat{v}$.
Given a boundary generic relative to $\hat{v}$ immersion $\alpha\colon X\to A$, the $\alpha$pullback (transfer) of the vector field $\hat{v}$ defines a nonvanishing vector field $v$ on $X$. Since $\alpha$ is an immersion, evidently the pullback $v$ is traversing on $X$. Moreover, $v$ is traversally generic in the sense of Definition 5.1, since no $v$trajectory $\gamma$ has more than one point of simple tangency to $\partial X$.
Definition 9.3.
Let $\alpha\colon X\subset A$ be a regular embedding of a given compact surface $X$ into an annulus $A$, carrying a traversal boundary convex vector field $\hat{v}$. We denote by $v$ the pullback of $\hat{v}$ under $\alpha$. If $\alpha$ is traversally generic relative to $\hat{v}$, then we say that the pair $(A,\hat{v})$ is a convex envelop of $(X,v)$.
The existence of a convex envelop puts significant restrictions on the topology of $X$: such orientable surfaces $X$ do not have 1handles. In other words, they are disks with holes.
Lemma 9.1.
If a compact connected surface $X$ with boundary has a pair of loops whose transversal intersection is a singleton, then no traversal flow on $X$ admits a convex envelop. In other words, if a connected surface $X$ with boundary has a handle, then no traversal flow on $X$ can be convexly enveloped.
Proof.
By Lemma 1.2, the space $\hat{X}$ of a convex envelop is either a disk or an annulus, both surfaces residing in the plane. No two loops in the plane intersect transversally at a singleton. Thus, for surfaces with a handle, no convex envelops exist. $\square$
To incorporate surfaces with handles into our constructions, we have introduced the notion of a convex quasienvelop (Definition 9.2).
Now we are in position to explore a connection between immersions $\alpha\colon(X,v)\to(A,\hat{v})$ of a given surface $X$ in the annulus $A$, such that $v=\alpha^{\ast}(\hat{v})$ and $\hat{v}$ is generic with respect to $\alpha(\partial X)$ on one hand, and loops in the functional spaces $\mathcal{F}_{\leq 2}$ on the other.
Let $\alpha(\partial X)^{\times}$ denote the set of selfintersections of the curves forming the image $\alpha(\partial X)$. Let $\alpha(\partial X)^{\circ}$ denote the set $\alpha(\partial X){\setminus}\alpha(\partial X)^{\times}$.
We associate an auxiliary smooth function $z_{\alpha}\colon A\to\mathbb{R}$ with the curves $\alpha(\partial X)$, subject to the following properties:

•
$z_{\alpha}^{1}(0)=\alpha(\partial X),$

•
0 is the regular value of $z_{\alpha}$ at the points of $\alpha(\partial X)^{\circ},$

•
in the vicinity of each crossing point $a\in\alpha(\partial X)^{\times}$, consider local coordinates $(x_{1},x_{2})$ such that $\{x_{1}=0\}$ and $\{x_{2}=0\}$ define the two intersecting branches of $\alpha(\partial X)$; then locally $z_{\alpha}=c\cdot x_{1}x_{2}$, where the constant $c\neq 0$.

•
$z_{\alpha}=1$ in the vicinity of $\partial A$,
(9.1)
The sign of $z_{\alpha}$ changes to the opposite, as a path crosses an arc from $\alpha(\partial X)^{\circ}$ transversally, thus providing a “checker board” coloring of the domains in $A\setminus\alpha(\partial X)$.
We denote by $A^{\circ}$ the interior of the annulus $A$. Let $\phi\colon A^{\circ}\to\mathbb{R}$ be a smooth function so that $d\phi(\hat{v})>0$ in $A^{\circ}$ and $\phi(\hat{\gamma}\cap A^{\circ})=\mathbb{R}$ for all $\hat{v}$–trajectories $\hat{\gamma}$ in $A$. Then, with the help of $z_{\alpha}$ and $\phi$, we get a map $J_{z_{\alpha}}\colon\mathcal{T}(\hat{v})\to\mathcal{F}_{\leq 2}$. We define the map $J_{z_{\alpha}}$ by the formula
$\displaystyle J_{z_{\alpha}}(\hat{\gamma})=(z_{\alpha}_{\hat{\gamma}})\circ( \phi_{\hat{\gamma}})^{1},$  (9.2) 
where, abusing notations, $\hat{\gamma}$ stands for both a $\hat{v}$trajectory in $A$ and for the corresponding point in the trajectory space $\mathcal{T}(\hat{v})\approx S^{1}$.
For a fixed $\alpha$, it is easy to check that the homotopy class $[J_{z_{\alpha}}]$ of $J_{z_{\alpha}}$ does not depend on the choice of the auxiliary function $z_{\alpha}$, subject to the four properties in (9.1). Indeed the space of such functions $z_{\alpha}$ is convex and thus contractible.
We pick a generator $\kappa\in\pi_{1}(\mathcal{F}_{\leq 2})\approx\mathbb{Z}$ and define the integer $J^{\alpha}$ by the formula $J^{\alpha}\cdot\kappa=[J_{z_{\alpha}}]$. As a result, any immersion $\alpha\colon X\to A$, which is generic with respect to $\hat{v}$, produces a homotopy class $[J_{z_{\alpha}}]\in\pi_{1}(\mathcal{F}_{\leq 2})$ and an integer $J^{\alpha}$.
The isomorphism $\pi_{1}(\mathcal{F}_{\leq 2})\approx\mathbb{Z}$ follows from [Arnold1989] by a slight modification of the arguments there.
Let us start with the key observation in [Arnold1989]. Generic loops in $\beta\colon S^{1}\to\mathcal{F}_{\leq 2}$ have an interpretation in terms of finite collections $C$ of smooth closed curves in the annulus $A$ with no inflection points with respect to their tangent lines of the form $\{\theta=const\}$ in the $(u,\theta)$coordinates. We call such tangent lines $\theta$vertical. Furthermore, the generic homotopy between such loops $\beta$ corresponds to some cobordism relation between the corresponding plane curves, the cobordism also avoids the $\theta$vertical inflections.
First, let us spell out the genericity requirements on the collections $C$ of closed curves in the annulus $A$ (carrying the angular coordinate $\theta$ and the radial coordinate $u$):

(1)
$C\subset A$ is a finite collection of closed smooth immersed curves $\{C_{j}\}_{j}$,

(2)
the projections $\{\theta\colon C_{j}\to S^{1}\}_{j}$ have the Morse type singularities only^{16}^{16}This allows for local maxima/minima only and excludes the $\theta$vertical inflections.,

(3)
the selfintersections and mutual intersections of the curves $\{C_{j}\}_{j}$ are transversal and no triple intersections are permited,

(4)
at each double intersection, the two banches of $C$ are not parallel to the $u$coordinate,

(5)
the $\theta$images of the intersections and of the critical values of $\{\theta\colon C_{j}\to S^{1}\}_{j}$ are all distinct in $S^{1}$.
(9.3)
Definition 9.4.
Given two collections $C_{0}$ and $C_{1}$ of immersed closed curves as in (9.3), we say that they are cobordant with no $\theta$vertical inflections, if there is a smooth function $F:A\times[0,1]\to\mathbb{R}$ such that:

•
0 is a regular value of $F$,

•
the restriction of the projection $T:A\times[0,1]\to[0,1]$ to the zero set $W=_{\mathsf{def}}F^{1}(0)$ is a Morse function,

•
$C_{0}=W\cap(A\times\{0\})$ and $C_{1}=W\cap(A\times\{1\})$,

•
for each $t\in[0,1]$, the section $C_{t}=_{\mathsf{def}}W\cap(A\times\{t\})$ has no $\theta$horizontal inflections^{17}^{17}Note that the second bullet excludes the triple intersections of $C_{t}$..
It is possible to verify that the cobordism with no $\theta$vertical inflections is an equivalence relation among collections of curves as in (9.3). Indeed, if $C$ is cobordant to $C^{\prime}$ with the help of $F$, and $C^{\prime}$ to $C^{\prime\prime}$ with the help of $F^{\prime}$, then there exists a piecewise smooth function $F\cup F^{\prime}\colon A\times[0,2]\to\mathbb{R}$ whose restriction to $A\times[0,1]$ is $F$ and to $A\times[1,2]$ is $1+F^{\prime}$. Smoothing $F\cup F^{\prime}$ along $A\times\{1\}$ in the normal direction and scaling down the interval $[0,2]$ to $[0,1]$, produces the desired functioncobordism $F\ast F^{\prime}\colon A\times[0,1]\to\mathbb{R}$.
So we can talk about the set of bordisms $\mathbf{B}_{\mathsf{no\,\thetainflect.}}$, based on collections of closed curves in the annulus with no $\theta$vertical inflections. This set is a group: the operation $C,C^{\prime}\Rightarrow C\ast C^{\prime}$ is defined by the union $\tilde{C}\cup\tilde{C}^{\prime}\subset A$, where $\tilde{C}\subset S^{1}\times(0,0.5)$ and $\tilde{C}^{\prime}\subset S^{1}\times(0.5,1)$ are the images of $C$ and $C^{\prime}$, scaled down in the $u$direction by the factor 0.5 and placed in subannuli of $A=S^{1}\times[0,1]$. The role of $C$ is played by the mirror image of $C$ with respect to a vertical (equivalently, horizontal) line, a fiber of $\theta\colon A\to S^{1}$.
Note that the operation $\ast$ may affect the maximal cardinalities $d$ and $d^{\prime}$ of the fibers $\theta\colon C\to S^{1}$ and $\theta:C^{\prime}\to S^{1}$ in a somewhat unpredictable way. In any case, the fiber cardinality of $\theta\colon C\ast C^{\prime}\to S^{1}$ has the upper boundary $d+d^{\prime}$.
The previous constructions deliver the following slight modification of Theorem from [Arnold1989].
Theorem 9.1.
The fundamental group $\pi_{1}(\mathcal{F}_{\leq 2})$ is isomorphic to the bordism group $\mathbf{B}_{\mathsf{no\,\thetainflect.}}$, based on finite collections of immersed loops with no $\theta$vertical inflections in the annulus $A$ and subject to the constraints (9.3). The isomorphism is induced by the correspondence
$\displaystyle K:\{\beta:S^{1}\to\mathcal{F}_{\leq 2}\}\;\Rightarrow\bigcup_{ \theta\in[0,2\pi]}(\theta,\beta(\theta)^{1}(0))\subset S^{1}\times[0,1].$ 
This theorem is a foundation of a graphic calculus that converts homotopies of loops in the functional space $\mathcal{F}_{\leq 2}$ into cobordisms of closed loop patterns with no $\theta$vertical inflections in the annulus $A$. The reader may glance at Figs. 10–14 from [Katz2015] to get a feel how the calculus works. These figures explain why any loop in $\mathcal{F}_{\leq 2}$ is homotopic to an integral multiple of a generator $\kappa\in\pi_{1}(\mathcal{F}_{\leq 2})$, represented by a model loop pattern $K\subset A$ as in Fig. 9, diagram (a) or (b).
Let us orient the annulus $A=S^{1}\times[0,1]$ so that the the $\theta$coordinate, corresponding to $S^{1}$, is the first, and the $u$coordinate, corresponding to $[0,1]$, is the second.
We fix an orientation of $X$, thus picking orientations for each component of $\partial X$. Given an orientationpreserving immersion $\alpha:(X,v)\to(A,\hat{v})$ such that $\alpha(\partial X)$ has the properties as in (9.3), we assign the polarity “$\oplus$” to a point $a\in\partial_{2}^{+}X(v)$ if $\alpha_{\ast}(\nu_{a})$ points in the direction of $\theta$. Here where $\nu_{a}$ is the inner normal to $\partial X$ at $a$. Otherwise, the polarity of $a\in\partial_{2}^{+}X(v)$ is defined to be “$\ominus$” (see Fig. 9). Thus, using new polarities, we get a partition $\partial_{2}^{+}X(v)=\partial_{2}^{+,\oplus}X(v)\coprod\partial_{2}^{+,\ominus }X(v)$.
Finally, we are in position to state the main result of Sect. 9. Its proof can be found in [Katz2015]. It is an application of the Arnold surgery calculus of loops with no $\theta$vertical inflections.
Theorem 9.2.
Any orientationpreserving immersion $\alpha\colon(X,v)\to(A,\hat{v})$, giving rise to a convex quasienvelop, such that $\hat{v}$ is generic with respect to $\alpha(\partial X)$, produces a map $J_{z_{\alpha}}\colon S^{1}\to\mathcal{F}_{\leq 2}$ (see (9.2)). Its homotopy class $[J_{z_{\alpha}}]=J^{\alpha}\cdot\kappa$, where $\kappa$ denots a generator of $\pi_{1}(\mathcal{F}_{\leq 2})\approx\mathbb{Z}$.
The integer $J^{\alpha}$ can be computed by the formula:
$\displaystyle J^{\alpha}=\#\{\partial_{2}^{+,\oplus}X(v)\}\#\{\partial_{2}^{+ ,\ominus}X(v)\}$ 
and thus does not depend on $\alpha$ (as long as the transfer $\alpha^{\ast}(\hat{v})=v$).
Moreover, $J^{\alpha}\leq c_{2}^{+}(v)$, the complexity of the $v$flow.
Remark 9.1.
It is interesting to notice that the invariant $J^{\alpha}=\#\{\partial_{2}^{+,\oplus}(v)\}\#\{\partial_{2}^{+,\ominus}(v)\}$ reflects more the topology of the vector field $v=\alpha^{\ast}(\hat{v})$ than the topology of the surface $X$: in fact, any integral value of $J^{\alpha}$ can be realized by a traversally generic vector field $v$ on a 2disk $D$ which even admits a convex envelop! A portion of the boundary $\partial D$ looks like a snake (see Fig. 9) with respect to the vector field $\hat{v}$ of the envelop. For any $X$, the effect of deforming a portion of $\partial X$ into a snake is equivalent to adding several times a spike (an edge and a pair of univalent and trivalent verticies) to the graph $\mathcal{T}(v)$. Evidently, these operations do not affect $H_{1}(\mathcal{T}(v);\mathbb{Z})\approx H_{1}(X;\mathbb{Z})$.
In contrast, $\#\{\partial_{2}^{+,\oplus}(v)\}+\#\{\partial_{2}^{+,\ominus}(v)\}\geq 2\chi( X)$ has a topological significance for $X$.
For example, for $\alpha$ as in Fig. 8, $J^{\alpha}=0$. If we subject $\alpha$ to an isotopy that introduces a snakelike pattern of Fig. 9a, then for the new immersion $\alpha^{\prime}$, the invariant $J^{\alpha^{\prime}}=1$.
Remark 9.2.
Consider a connected oriented surface $X$ with a connected boundary. (We exclude the case of a 2disk.) It is a boundary connected sum of a few copies of $T^{\circ}$, the torus with a hole. A punctured torus admits an immersion $\alpha\colon T^{\circ}\to A$ in the annulus so that the cardinality of the fibers of $\theta:\alpha(\partial T^{\circ})\to S^{1}$ does not exceed 6 (see Fig. 8). Therefore, any connected oriented surface $X$ with boundary admits an immersion $\alpha\colon X\to A$ with the property $\#\{\theta^{1}(\theta_{\star})\cap\alpha(\partial X)\}\leq 6$ for all $\theta_{\star}\in S^{1}$.
Any $\hat{v}$generic immersion $\alpha\colon X\to A$ also produces a welldefined element $[K_{\alpha}]$ in the set of homotopy classes $[\mathcal{T}(v),\mathcal{F}_{\leq 2}]$ of maps from the trajectory graph $\mathcal{T}(v)$ to the functional space $\mathcal{F}_{\leq 2}$. Its construction is similar to the one of $J_{z(\alpha)}$. Consider the $\hat{v}$generated obvious map $Q_{\alpha}\colon\mathcal{T}(v)\to\mathcal{T}(\hat{v})\approx S^{1}$ (each $v$trajectory is contained in the unique $\hat{v}$trajectory). Put $K_{\alpha}=_{\mathsf{def}}J_{z(\alpha)}\circ Q_{\alpha}$.
Remark 9.3.
Note that, for some immersions $\alpha\colon X\to A$, the invariant $J^{\alpha}$ may be different from 0, but $[K_{\alpha}]$ may be trivial. For example, this is the case when $X$ is a disk with a snakelike boundary $\alpha(\partial X)$ with respect to $\hat{v}$. However, there exist immersions $\alpha$ with a nontrivial $[K_{\alpha}]$. For example, such is the immersion in [Katz2015], Figure 10, (1). At the same time, for $\alpha$ in [Katz2015], Figure 11, (3), $[K_{\alpha}]$ is trivial.
Since $\pi_{1}(\mathcal{F}_{\leq 2})\approx\mathbb{Z}$, it follows that $H_{1}(\mathcal{F}_{\leq 2};\mathbb{Z})\approx\mathbb{Z}$. In turn, this implies that the 1dimensional cohomology $H^{1}(\mathcal{F}_{\leq 2};\mathbb{Z})\approx\mathbb{Z}$. Thus $K_{\alpha}$ induces a homomorphism
$\displaystyle K_{\alpha}^{\ast}:H^{1}(\mathcal{F}_{\leq 2};\mathbb{Z})\to H^{1 }(\mathcal{T}(v);\mathbb{Z})\approx H^{1}(X;\mathbb{Z}).$ 
In particular, we get an element $K_{\alpha}^{\ast}(\kappa)\in H^{1}(X;\mathbb{Z})$, where $\kappa$ is a generator of $H^{1}(\mathcal{F}_{\leq 2};\mathbb{Z})\approx\mathbb{Z}$. $K_{\alpha}^{\ast}(\kappa)$ is a characteristic cohomology class of the given $\hat{v}$generic immersion $\alpha$. It is an invariant of the class $[\alpha]\in\mathbf{B}_{\mathsf{no\,\thetainflect.}}$.
Theorem 9.2 implies that, if two $\hat{v}$generic immersions $\alpha_{1},\alpha_{2}:X\to A$ are such that the pullbacks $\alpha_{1}^{\ast}(\hat{v})=\alpha^{\ast}_{2}(\hat{v})=v$, then $K_{\alpha_{1}}^{\ast}(\kappa)=K_{\alpha_{2}}^{\ast}(\kappa)$.
So the cohomology class $K_{\alpha}^{\ast}(\kappa)$ turns out to be a characteristic class of $v$! It is desirable to be able to reach this conclusion without relying on the cobordisms of curves’ patterns in $A$ with no $\theta$vertical inflections.
The space $\mathcal{F}_{\geq 3}$ has codimension 2 in $\mathcal{F}$; so a loop in $\mathcal{F}$ may be linked with $\mathcal{F}_{\geq 3}$. In fact, the value of $K_{\alpha}^{\ast}(\kappa)$ on any loop (1cycle) $\delta\colon S^{1}\to X$ equals to the linking number $\mathsf{lk}(K_{\alpha}(\delta),\mathcal{F}_{\geq 3})$. The validation of this claim and of its generalizations will appear in [Katz2016c]; [Katz2017].
In dimensions higher than two, similar considerations apply to produce characteristic classes of traversally generic flows. They are based on quite involved computations of homology of spaces of real monic polynomials with restricted combinatorics of their real divisors.
As a result of these computations, it turns out that the topology of highdimensional convex envelops is at least as intricate as the homotopy groups of spheres [Katz2016c].
10 A Glance at the $3D$Lands
Here is brief overview of how some elements of our program manifest themselves in dimension three.
Let $v$ be a traversally generic (see Definition 5.1 for the $2D$flows and its generalization in [Katz2009]; [Katz2014b] for the $n$dimensional flows) vector field on a compact connected 3dimensional manifold $X$ with boundary. The trajectories of such vector fields have the combinatorial tangency patterns from the poset $\Omega^{\bullet}_{{}^{\prime}\langle 2]}=\{(11);\,(2),(121);\,(1221),(31),(13)\}$.
Consider the number $gc(v)$ of $v$trajectories that are simply tangent (have tangency of order 2 in the sense of Definition 7.1) to the boundary $\partial X$ at two distinct points. We call this number $gc(v)$ the gradient complexity of $v$. It is an example of an invariant of the vector field $v$, considered up to the natural action of diffeomorphism group $\mathsf{Diff}(X)$ on vector fields. We minimize $gc(v)$ over all traversally generic vector fields on $X$ to produce the gradient complexity $gc(X)$, a new invariant of $X$.
Recall that any compact 3manifold can be assembled out of the standard building blocks (in technical terms, these blocks are “boundary irreducible” with no “essential annuli” compact 3manifolds) by gluing them along their boundary components [Jaco and Shalen1979]. It turns out that there exist only finitely many standard blocks whose gradient complexity $gc(\sim)$ is bounded from above ([Katz2009], Theorem 7.3).
Let $Y$ be a closed hyperbolic 3manifold, and let $X$ be obtained by removing a 3ball from $Y$. In [Katz2009], Theorem 7.5, we prove that the gradient complexity of $X$ is greater than or equal to the normalized hyperbolic volume of $Y$. In fact, this type of result is valid in any dimension ([Alpert and Katz2015], Theorem 1). In particular, compare this claim with Theorem 4.1, the 2dimensional analogue.
Our investigation of vector flows in Flatland has reached its conclusion. To find out how things flow in other lands—“the romances of many dimensions”—([Abbott1992]), the reader could consult with the references below.
Acknowledgements
We are grateful to the referees; their thoughtful and detailed recommendations helped to improve the quality of our presentation.References
 [Abbott1992] Abbott, E.: Flatland: A romance of many dimensions. Dover Publications, New York (1992)
 [Alpert and Katz2015] Alpert, H., Katz, G.: Using simplicial volume to count multitangent trajectories of traversing vector fields. Geometriae Dedicata 180(1), 323–338 (2015)
 [Arnold1989] Arnold, V.I.: Spaces of functions with moderate singularities, Func. Analysis & Its Applications, 23(3):1–10 (1989) (Russ)
 [Barannikov1994] Barannikov, S.A.: The framed Morse complex and its invariants, Advances in Soviet Mathematics, Vol. 21, V. I. Arnold Ed., 93115 (1994)
 [Cohen1991] Cohen R.L.: Topics in Morse theory Stanford University (1991)
 [Goresky and MacPherson1983a] Goresky, M., MacPherson, R.: Stratified Morse Theory. Proceedings of Symposia in Pure Mathematics 40(Part 1):517–533 (1983a)
 [Goresky and MacPherson1983b] Goresky, M., MacPherson, R.: Morse theory for the intersection homology groups, Analyse et Topologie sur les Espaces Singulieres, Astérisque #101 (1983b), 135–192, Société Mathématique de France
 [Goresky and MacPherson1989] Goresky, M., MacPherson, R.: Theory, Stratified Morse, Springer Verlag, N.Y.: Ergebnisse, vol. 14. Also translated into Russian and published by MIR Press, Moscow (1989). 1991
 [Gottlieb1996] Gottlieb, D.H.: All the way with GaussBonnet and the sociology of mathematics. Math. Monthly 103, 457–469 (1996)
 [Gromov1982] Gromov, M.: Volume and bounded cohomology. Publ. Math. I.H.E.S., tome 56:5–99 (1982)
 [Guth2009] Guth, L.: Minimal number of selfintersections of the boundary of an immersed surface in the plane. [math.DG] (2009)
 [Hardon2004] Hardon, D.G.C.: The Morse Complex for a Morse Function on a Manifold with Corners. (2004)
 [Hopf1937] Hopf, H.: Vectorfelder in $n$dimensionalen Mannigfaltigkeiten. Math. Annalen 96, 225–250 (1937)
 [Jaco and Shalen1979] Jaco, W.H., Shalen, P.B.: Seifert fibered spaces in 3manifolds. Mem. Am. Math. Soc. 21(220) (1979)
 [Katz2009] Katz, G.: Convexity of Morse Stratifications and gradient spines of 3manifolds. JP J. Geometry Topol. 9(1), 1–119 (2009)
 [Katz2014a] Katz,G.: Stratified Convexity and Concavity of Gradient Flows on Manifolds with Boundary. Appl. Math. 5, 2823–2848 (2014a). http://www.scirp.org/journal/am (also [mathGT]
 [Katz2014b] Katz, G.: Traversally Generic and Versal Flows: Semialgebraic Models of Tangency to the Boundary. Asian J. Math. (2014b, to appear). [mathGT]
 [Katz2014c] Katz, G.: Causal Holography of Traversing Flows (2014c). [mathGT]
 [Katz2015] Katz, G.: Flows in Flatland: a Romance of few Dimensions (2015). [math.GT]
 [Katz2016a] Katz, G.: The Stratified Spaces of Real Polynomials and Trajectory Spaces of Traversing Flows. JP J. Geometry Topol. 19(2), 95–160 (2016a). [math.GT]
 [Katz2016b] Katz, G.: Complexity of Shadows and Traversing Flows in Terms of the Simplicial Volume. J. Topol. Anal. 8(3), 501–543 (2016b). [mathGT]
 [Katz2016c] Katz, G.: Spaces of Functions with Constrained Zero Divisors, Convex Envelops of Traversing Flows, and their Characteristic Classes, manuscript (2016c)
 [Katz2017] Katz, G.: Morse Theory, Gradient Flows, Concavity, and Complexity on Manifolds with Boundary, monograph. World Sci. (2017)
 [Laudenbach2011] Laudenbach, F.: Morse Complex on Manifolds with Boundary. Geom. Dedicata 153, 47–57 (2011)
 [Malgrange1964] Malgrange, B.: The preparation theorem for differential functions. Differential Analysis (Papers presented at the Bombay Colloquium), pp. 203–208 (1964)
 [Milnor1965] Milnor, J.: Morse Theory. Princeton University Press, Princeton (1965)
 [Morse1929] Morse, M.: Singular points of vector fields under general boundary conditions. Am. J. Math. 51, 165–178 (1929)
 [Morse and Van Schaack1934] Morse, M., Van Schaack, G.B.: The critical point theory under general boundary conditions. Ann. Math. (2) 35(3):545571 (1934)
 [Teschi2012] Teschi, G.: Ordinary differential equations and dynamical systems. Graduate Studies in Math., 140, AMS publication (2012)
 [Vasiliev1994] Vasiliev, V.A.: Complements of Discriminants of Smooth Maps: Topology and Applications, Translations of Mathematical Monographs, vol. 98. Society publication, American Math (1994)
 [Whitney1937] Whitney, H.: On regular closed curves in the plane. Comp. Math. 4, 276–284 (1937)