On cusps of caustics by reflection: a billiard variation on Jacobi's Last Geometric Statement

Abstract: A point source of light is placed inside an oval. The $n$-th caustic by reflection is the envelope of the light rays emanating from the light source after $n$ reflections off the curve. We show that each of these caustics, for a generic point light source, has at least 4 cusps. This is a billiard variation on Jacobi's Last Geometric Statement, concerning the number of cusps of the conjugate locus of a point on a convex surface. We present various proofs, using different ideas, including the curve shortening flow and Legendrian knot theory.

\section{Introduction}
\subsection{Motivation and background}
The {\em conjugate locus} of a point on a surface is the locus of first conjugate points along  geodesics emanating from that point.
In his ``Lectures on Dynamics" \cite{Cl}, published posthumously,
Jacobi stated that  the conjugate locus of a generic point on an ellipsoid has exactly four cusps.
This {\it Last Geometric Statement of Jacobi} was proved only in this century, see \cite{IK}. Indeed, as
recently as the end of the 20th century, Marcel Berger wrote \cite{Ber}:
\begin{quote}
{\small
... this latter assumption depends on the scandalously unproved Jacobi ``statement": the conjugate locus of a non-umbilical point of an ellipsoid has exactly four cusps.}
\end{quote}

A related result is that the conjugate locus of a generic point on a convex surface has at least four cusps, see \cite{Wa} for a recent proof. This theorem was attributed to C. Carath\'eodory (1912) by W. Blaschke (sect. 103 of \cite{Bl}), who presented a sketch of the proof. This theorem belongs to a long list of results that stem from and are motivated by the celebrated 4-vertex theorem of S. Mukhopadhyaya. See \cite{Ar,GTT} for surveys.

The conjugate locus can be equivalently described as the locus of the first intersections of infinitesimally close geodesics emanating from a point. These geodesics may intersect more than once, and the loci of their intersections are known as second, third, etc., caustics of the point. Statements similar to Jacobi's statement and generalization to arbitrary convex surfaces about the  higher order caustics are still open. There is some experimental evidence that if the surface is an ellipsoid then, for a non-umbilic point, each such caustic has exactly four cusps, see \cite{Si} and Figure \ref{ellipsoids} from this paper (presented with permission).

\begin{figure}[ht]
\centering
\includegraphics[width=.3\textwidth]{./ellipsoids_left.png}\hspace{.2\textwidth}
\includegraphics[width=.3\textwidth]{./ellipsoids_right.png}
\caption{The first three caustics of an non-umbilic point on an ellipsoid.}
\label{ellipsoids}
\end{figure}

In this article we consider a billiard version of this  problem. Let $\gamma$ be an oval (a smooth strictly convex closed curve in ${\mathbb R}^2$), the boundary of a billiard table or, equivalently, an ideal mirror. Let $O$ be a point inside $\gamma$, a source of light. For $n=1,2.\ldots$, the 1-parameter family of rays that have undergone $n$ optical reflections in $\gamma$ envelopes a curve $\Gamma_n$, the {\em $n$-th caustic by reflection}. See Figure \ref{fig:caustics}.

\begin{figure}[ht]
\centering
\figsize{.4}{3caustics}
\caption{The   $n$-th caustic by reflection $\Gamma_n$ is the envelope of the family of rays emanating from $O$ that have undergone $n$ reflections by $\gamma$.}
\label{fig:caustics}
\end{figure}

These caustics may have singularities, generically, semi-cubical cusps. We always assume that the caustics $\Gamma_n$ are in general position in this sense. The singularities of caustics were thoroughly studied  by Bruce, Giblin, and Gibson; see \cite{BGG} and the references therein.

Figure \ref{fig:cc}a shows that a caustic by reflection may extend beyond the interior of $\gamma$, and furthermore,  it can be disconnected in the Euclidean plane; however, as the envelope of a 1-parameter family of lines, it is a connected curve in the projective plane ${\mathbb {RP}}^2$ (possibly, with singularities). Indeed, a 1-parameter family of lines is a curve in the space of lines, and the respective envelope is projectively dual to this curve.

\begin{figure}[ht]
\centering
\figsize{.45}{circle_caustics1c}
\hspace{.15\textwidth}
\figsize{.32}{circle_caustics2

\caption{(a)  The first three caustics by reflection in a circle, showing 4 cusps on each of them (red=1st caustic, blue=2nd, green=3rd). The small circle is the light  source. (b) The first caustic by reflection  in a circle with an external source of light (each ray optically reflects at both intersection points with the circle). The caustic in this case has only 2 cusps. In this paper we  consider only an internal light source.}
\label{fig:cc}
\end{figure}

\subsection{The main result and two conjectures}
Our main result is as follows.

\begin{theorem} \label{thm:main}
For every oval $\gamma\subset{\mathbb R}^2$, a generic light source  inside $\gamma$ and  $n\ge 1$, the $n$-th caustic by reflection   $\Gamma_n\subset{\mathbb {RP}}^2$ has at least four cusps.
\end{theorem}

We present  three proof sketches.

Let $\L$ be the space of directed lines in  ${\mathbb R}^2$. To each caustic  $\Gamma_n$ is associated its dual  curve  $C_n\subset\L$, corresponding to  the tangent lines along  $\Gamma_n$ (the rays of the $n$-th reflected beam). One can identify  $\L$ with the complement of the `north' and `south' pole of the unit sphere $S^2\subset{\mathbb R}^3$,  so that cusps of $\Gamma_n$ correspond to   inflection points of $C_n$ (points with vanishing   spherical geodesic curvature). Using  standard properties of convex billiards, we show that $C_n$ is a closed simple smooth curve in  $S^2$, intersecting every great circle. A theorem of B. Segre from 1968  \cite{Se,Wei}  states  that such a curve has at least four spherical inflection points, thus completing the proof of Theorem \ref{thm:main}.

Another approach, starts with a realization of $\L$ as the vertical cylinder circumscribing $S^2$ and the curve $C_n \subset \L$ representing  the tangent lines of $\Gamma_n$.
Following S. Angenent \cite{An}, apply
the curve shortening flow with respect to the flat metric on the cylinder to the curve $C_n$ to deform it  to the  graph of a function  $F:S^1\to{\mathbb R}$ with zero mean value. Spherical inflection points of $C_n$ correspond  to  the zeroes of $F''+F$, a function with vanishing constant  and first order Fourier terms. By the Sturm-Hurwitz theorem, it has at least four zeros.

Yet another approach is to use the relation between the cusps of the caustic $\Gamma_n$ and the vertices (critical points of the curvature) of its normal front $\Delta_n$, a closed planar curve whose normal lines, parametrized by $C_n$, are the lines tangent to $\Gamma_n$, see Figure \ref{fig:3curves}. The relation bewteen $\Gamma_n$ and $\Delta_n$ is the familiar relation between evolutes and involutes, see, e.g., \cite{GTT}.

 \begin{figure}[ht]
\centering
\fig{involute2}
\caption{The 2nd caustics $\Gamma_2$ (red) with an involute $\Delta_2$ (green) for an elliptical billiard table $\gamma$ (black). The rays (gray) correspond to $C_2$, are normal to the wave front $\Delta_2$ and are tangent to $\Gamma_2$.  }
\label{fig:3curves}
\end{figure}

We show that $\Delta_n$ exists as a closed curve, possibly with cusps (in fact, there is a 1-parameter equidistant family of such curves).  To show that $\Delta_n$ has at least four  vertices, we use a theorem of  Chekhanov and Pushkar \cite{ChP}, stating that a planar cooriented closed wave front has at least four vertices, provided its Legendrian lift to the space of cooriented contact element in the plane  is a `Legendrian unknot', that is, is Legendrian isotopic to the Legendrian lift of a circle.

The rest of the article provides background information and details of these arguments. In the last section we present some generalizations of Theorem \ref{thm:main} to spherical and hyperbolic geometry, as well as ``projective billiards.''

\medskip\noindent

We present two conjectures; the first one is supported by experimental evidence, the  second one might be over-optimistic.

\begin{conjecture} \label{conj:one}
If $\gamma$ is an ellipse, then the caustic by reflection $\Gamma_n$ for a  light source inside $\gamma$ and different from a focus has exactly four cusps for every $n\ge 1$ (see Figure \ref{ellipse1}).
\end{conjecture}

This conjecture is only known to hold in the case of $n=1$ (the `catacaustic', see next section).

\begin{figure}[ht]
\centering
\includegraphics[width=.3\textwidth]{./caustic2}\quad
\includegraphics[width=.3\textwidth]{./caustic5}\quad
\includegraphics[width=.3\textwidth]{./caustic8}
\caption{The 2nd, 5th and 8th caustics by reflection in an ellipse, each with 4 cusps (marked by gray disks).}
\label{ellipse1}
\end{figure}

\begin{conjecture} \label{conj:two}
If $\gamma$ is not an ellipse then, for some choice of light source  inside $\gamma$  and some $n\ge 1$, the caustic by reflection $\Gamma_n$ has more than four cusps.
\end{conjecture}

Figure \ref{fig:4plus} shows caustics  with $>4$ cusps for non-elliptical $\gamma$ .

\begin{figure}[ht]
\centering
\includegraphics[width=.4\textwidth]{./cata4.pdf}\hspace{.15\textwidth}
\includegraphics[width=.4\textwidth]{./cataegg.pdf}
\caption{First caustic by reflection with more than four cusps of non-elliptical ovals. Left: $x^4+y^4=1$, $O=(.6, .4)$. Right: $.5x^2+(1+.25x)y^2=1$, $O=(.5, .3)$ }
\label{fig:4plus}
\end{figure}

\subsection{Catacaustics}
The first caustics by reflections, called {\em catacaustics}, are well studied. We give a brief summary of what is known about them, referring to \cite{Bo,BGG,Hu} and the literature cited in these articles.

A version of the string construction that recovers a billiard curve from a billiard caustic (see, e.g., \cite{Ta}) makes it possible to reconstruct the curve $\gamma$ from its first caustic by reflection $\Gamma_1$. This construction involves a parameter, the length of the string. See Figure \ref{string} (left).

\begin{figure}[ht]
\centering
\figsize{.4}{string}
\hspace{.2\textwidth}
\figsize{.32}{osculating_kepler_conic1}
\caption{(a) The curve $\gamma$ is the locus of points $X$ such that $|OX|+ |XB| + \stackrel{\smile}{|BA|}=const$ (point $A$ is fixed on $\Gamma_1$). The point $Z$ is the reflection of $O$ in the tangency line to $\gamma$ at   $X$. The locus of points $Z$ is the orthotomic curve $\Delta_1$, orthogonal to $BZ$ and whose evolute is $\Gamma_1$. (b) The catacaustic $\Gamma_1$ (red) is the locus of 2nd foci $B$ of the osculating Kepler conic (gray) to the curve $\gamma$ (blue) with 1st foci at $O$.
}
\label{string}
\end{figure}

The orthotomic curve $\Delta_1$ is an involute of the catacaustic $\Gamma_1$,  see  Figure \ref{string} (left), and $\Gamma_1$ is the evolute of $\Delta_1$, that is, the envelope of its normals. The  cusps of $\Gamma_1$ correspond to the vertices of $\Delta_1$.
 It is known that when $\gamma$ is an ellipse and $O$ is not one of its foci then $\Gamma_1$ has 4 vertices \cite{BGG}.

 A Kepler conic is a conic with one focus fixed at the origin $O$. The curve $\gamma$ has a Kepler conic that has 3-point contact with it at every point, see \cite{BT}. The locus of the second foci of these osculating Kepler conics is the first caustic $\Gamma_1$ -- this follows from the optical properties of conics (a ray from one focus reflects to another focus). It follows that the cusps of the catacaustics $\Gamma_1$ correspond to the points where the Kepler conics hyperosculate the curve $\gamma$.

 \paragraph{Computer graphics and animations.} Most figures in this article were made
using the computer program Mathematica. They are complemented with some
animations on the web page   \url{https://www.cimat.mx/~gil/caustics/}.

\paragraph{Acknowledgements:}
GB acknowledges  support from the Shapiro visiting program in Penn State and   CONACYT Grant A1-S-4588. ST was supported by NSF grant DMS-2005444.

\section{Background material}

\subsection{The  phase cylinder and the billiard ball map}\label{sect:pc}

This section contains some standard material on mathematical billiards, see, e.g., \cite{Ta}.

Denote by $\L$ the space of oriented lines in ${\mathbb R}^2$. We use the `cylinder model' of $\L$, with  coordinates $(\alpha,p)$ defined as follows: $\alpha\in S^1={\mathbb R}/2\pi{\mathbb Z}$ is the direction of the line and $p\in {\mathbb R}$ is the signed distance from the oriented line to the origin $O$ (which we choose to be the center of the initial beam of light). The sign of $p$ is defined by the right-hand rule, see Figure \ref{lines}.  Thus $\L$ is an infinite cylinder.

\begin{figure}[ht]
\centering
\fig{lines}
\caption{The coordinates $(\alpha,p)$ of an oriented line in ${\mathbb R}^2$.}
\label{lines}
\end{figure}

The space $\L$ of oriented lines in ${\mathbb R}^2$ admits an area form, unique up to scale, invariant under the Euclidean group action. In coordinates, this area form is $\omega=d\alpha\wedge dp.$

The {\em phase cylinder} of the billiard system inside an oval $\gamma\subset{\mathbb R}^2$ is  the set $M\subset \L$ of oriented lines intersecting  $\gamma$. It is a bounded cylinder whose two  boundary components  correspond to the lines tangent to $\gamma$, one component for each orientation.  The ``equator'' $p=0$ corresponds to the lines through $O$, see Figure \ref{fig:pc}.

\begin{figure}[ht]
\centering
\figsize{.3}{phase_cylinder6}
\hspace{.2\textwidth}
\figsize{.4}{T}
\caption{Left: the phase cylinder $M\subset \L$. Right: the billiard ball map $T:M\to M$. }
\label{fig:pc}
\end{figure}

The billiard ball map $T:M\to M$, sending  an incoming ray to the reflected one, is an area preserving transformation, that is, $T^*\omega=\omega$.
Since $\omega=-d(pd\alpha)$, the differential 1-form $T^*(pd\alpha)-pd\alpha$ is closed. In fact, more is true: as we will now show, it is {\it exact}, that is, $T^*(pd\alpha)-pd\alpha=dF$ for some function $F:M\to{\mathbb R}$. An example of an area preserving, but non-exact, map   is   $(\alpha,p)\mapsto (\alpha, p+1).$

\begin{proposition}\label{prop:exact} The billiard ball map $T:M\to M$ is exact.
\end{proposition}

In order to  prove Proposition \ref{prop:exact},  consider another description of the phase cylinder, as the set of unit vectors with a foot point on $\gamma$, pointing inwards, the initial position and velocity of the billiard ball. These unit vectors are in one-to-one correspondence with the oriented lines that they generate.   Let $\gamma(t)$ be an arc length  counterclockwise parameterization and  $\varphi$ be the angle between the tangent $\gamma'(t)$ and the unit vector. See Figure \ref{fig:ellip}a.

\begin{figure}[ht]
\centering
\figsize{.4}{tphi}\hspace{.15\textwidth}
\figsize{.39}{chord}
\caption{(a) The coordinates $(t, \varphi)$ on $M$ and their relation to $(\alpha,p).$ (b) The generating function $L$ of the billiard ball map $T$.
}
\label{fig:ellip}
\end{figure}

Consider the differential 1-form $\cos\varphi\, dt$. Let $L=|\gamma(t_1)-\gamma(t)|$ be the distance between the intersection points of a line   with $\gamma$. See Figure \ref{fig:ellip}b.

\begin{lemma} \label{lm:exact}
$T^* (\cos\varphi\ dt)-\cos\varphi\ dt = dL$.
\end{lemma}

\begin{proof}
One has: $T(t,\varphi)=(t_1, \varphi_1)$ and
$$
\frac{\partial L(t,t_1)}{\partial t} = - \cos\varphi,\ \frac{\partial L(t,t_1)}{\partial t_1} =  \cos\varphi_1,
$$
that is,
$
dL= \cos\varphi_1\ dt_1 - \cos\varphi\ dt,
$
as needed. \end{proof}

Two differential 1-forms are cohomologous if their difference is the differential of a function.

\begin{lemma} \label{lm:cohom}
The 1-form $p\, d\alpha$ is cohomologous to $\cos\varphi\ dt$.
\end{lemma}

\begin{proof}
Using complex notation, see Figure \ref{fig:ellip}a, one has
$$
e^{i(\alpha-\varphi)}=\gamma'(t),\ p=\det( \gamma(t), e^{i\alpha} ).
$$
Differentiating  these equations, we get
$$
d\alpha\equiv d\varphi\ (\mbox{mod }dt),\ dp\equiv \sin\varphi\, dt\ (\mbox{mod }d\alpha),
$$ so
$$
d\alpha\wedge dp=\sin\varphi\, d\varphi \wedge dt=-d(\cos\varphi\ dt),
$$
hence $p d\alpha - \cos\varphi\ dt$ is closed.

To show that $p\, d\alpha - \cos\varphi\,dt$ is exact it suffices to show that its integral along a non-contractible closed curve in $M$ vanishes.
As such a curve take $C$, the boundary component of the phase cylinder $M$ given by $\varphi =0$. Clearly, $\int_C \cos\varphi\,dt$ equals the perimeter of $\gamma$.

On the other hand, by the Cauchy-Crofton formula, this perimeter equals $\pi$ times the average length of the orthogonal projection of $\gamma$ on a line, that is,  $\int_C p d\alpha$.
This implies the result.
\end{proof}

\noindent {\em  Proof  of Proposition \ref{prop:exact}.} Lemma \ref{lm:exact} shows that $T$  preserves the integral $\int_C\cos\varphi\ dt$, and Lemma \ref{lm:cohom} shows that $\int_C\cos\varphi\ dt=\int_C p\,d\alpha$, hence $T$ preserves the integral $\int_C p\,d\alpha$, i.e., it is exact.  \qed

\medskip\noindent

The {\em signed area} enclosed by an oriented  closed curve $C$ in $\L$ is the line integral $\int_C pd\alpha$. Clearly, the curve $C_0 $ representing the  initial beam of light encloses zero area. Since $T$ is an exact map, so does $C_n=T^n(C_0)$. See Figure \ref{fig:3beams}.

\begin{figure}[ht]
\centering
\fig{3beams}
\caption{The succesive iterates $C_n=T^n(C_0)$, for $n=1,2,3$,  drawn on the phase cylinder $\L$ (spread flat). Each is a simple closed smooth curve of 0 signed area (blue area = green area). The billiard table is the ellipse $4x^2/5 +y^2=1$ and the source is $(.6, .2)$.}
\label{fig:3beams}
\end{figure}

\subsection{Contact elements, Legendrian knots, and wave fronts}

A {\em contact element} in the plane is a pair consisting of a point and a line through it. A coorientation of a contact element is a choice of one of the sides of the line. More conceptually, the space of cooriented contact elements is the spherization of the cotangent bundle $ST^*{\mathbb R}^2$: assign to a covector $\eta$ its kernel, a tangent line, and define coorientation by choosing the side on which $\eta$ is positive.

The space of contact elements carries a contact structure, a 2-dimensional distribution defined by the ``skating condition": the foot point may move along the line, and the line may rotate about the foot point.
 Let $(x,y)$ be the standard  coordinates in ${\mathbb R}^2$ and $\theta$  the angle between the positive $x$-axis and the direction of the line;  then the contact distribution  is the kernel of the 1-form
$
\sin\theta\, dx-\cos\theta\, dy.
$

A smooth curve in $ST^*{\mathbb R}^2$ that is tangent to the contact distribution is called {\em Legendrian}. Its projection to the plane  is a {\em wave front}, a curve that may have singularities, generically semicubical cusps, but that has a tangent line at every point. Conversely, such a curve has a unique lift to the space of contact elements as a Legendrian curve.

We introduce   coordinates $(\alpha,p,z)$ in $ST^*{\mathbb R}^2$, where $(\alpha,p)\in T^*S^1$ are the coordinates of the orthogonal line at the foot point $A$, and $z$ is the (signed) distance  of the line to the origin $O$.  See Figure \ref{contact}.

\begin{figure}[ht]
\centering
\fig{contact}
\caption{Coordinates $(\alpha,p,z)$ on the space $ST^*{\mathbb R}^2$ of cooriented contact elements in ${\mathbb R}^2$. With a cooriented line $\ell$ through $A$ we associate the (signed) distance $z$ to $O$ and the coordinates $(\alpha,p)$ of the perpendicular oriented line at $A$, pointing  to the positive side of $\ell$.
}
\label{contact}
\end{figure}

This defines an identification of $ST^*{\mathbb R}^2$ with $J^1S^1$, the  space of 1-jets of functions $f:S^1\to{\mathbb R}$,  where $z=f(\alpha)$ and $p=f'(\alpha).$ On $J^1S^1$ there is a standard  contact form $dz-pd\alpha$ (the 1-jets of functions $z=f(\alpha)$ are Legendrian curves).

\begin{lemma} The identification $ST^*{\mathbb R}^2=J^1S^1$ is a contactomorphism.
\end{lemma}
\begin{proof} Using the notation of Figure \ref{contact},
$$
\alpha=\theta+\pi/2,\ p=\det(A, e^{i\alpha}),\ z=\det(A, e^{i\theta}),
$$
 so
 $$
 dz-pd\alpha=\det(dA, e^{i\theta})=\sin\theta dx-\cos\theta dy,
 $$
 as claimed.
 \end{proof}

We have two projections
\[
\begin{tikzcd}
& ST^*{\mathbb R}^2\arrow[dl,swap,"\pi_1"] \arrow[dr,"\pi_2"] \\
{\mathbb R}^2  && \L
\end{tikzcd}
\]
where  $\pi_1$ maps a cooriented contact element $(A,\ell)$ to $A$ and  $\pi_2$ maps it to the  line through $A$, orthogonal to $\ell$, oriented towards its positive side.

In coordinates, $\pi_2:
(\alpha,p,z)\mapsto (p,\alpha)$ (``forgetting $z$''). The fibers of  $\pi_1$ are Legendrian, spanned by the vector field $\partial_\alpha-z\partial_p+p\partial_z$, while the fibers of  $\pi_2$, spanned by the vector field $\partial_z$,  are transverse to the contact distribution. Hence $\pi_2$ projects a smooth Legendrian curve in $ST^*{\mathbb R}^2$ to a smooth curve in $\L$.

Conversely, let $C$ be a   closed  curve in $\L$. We want to lift it via $\pi_2:ST^*{\mathbb R}^2\to \L$ to a Legendrian curve $\widetilde C\subset ST^*{\mathbb R}^2$. Since the contact distribution on $ST^*{\mathbb R}^2$ is transverse to the fibers of $\pi_2$, once the initial point of the lifting is chosen, the lifting is uniquely determined, but it may fail to close up.

\begin{lemma} \label{lm:front}
The lifted curve $\widetilde C$ is closed if and only if $\int_C p d\alpha =0$: the curve $C$ encloses zero signed area.
\end{lemma}

\begin{proof} The curve
$\widetilde C$ is closed if and only if the value of the third coordinate $z$ is the same at the endpoints. Since $dz=p d\alpha$ along a Legendrian curve, the values of $z$ at the endpoints are equal if and only if  $\int_C p d\alpha =0$.
\end{proof}

The curve $C\subset \L$ defines a 1-parameter family of oriented lines.
The projection of the lifted Legendrian curve $\widetilde C\subset ST^*{\mathbb R}^2$ to ${\mathbb R}^2$ is a wave front $\Delta$ that is orthogonal to this family of lines and is cooriented by their directions. If a closed   front $\Delta$ exists, i.e., $C$ encloses zero signed area, then there exists a whole 1-parameter family of fronts that are equidistant from each other. This non-uniqueness corresponds to the choice of the initial point of the lifted curve $\widetilde C$.

The situation is the same as in the familiar  relation between evolutes and involutes: for an involute of a closed curve to close up it is necessary and sufficient for the curve to have zero signed length (the sign changes after each cusp), and the equidistant family of curves share their normals, and hence their evolutes.

\subsection{Vertices of wave fronts and Legendrian isotopies}

A {\em vertex} of a plane curve is an extremum of its curvature or, equivalently, a cusp of the evolute, the envelope of its normals. The notion of vertex extends to cooriented wave fronts: the curvature at cusps is infinite, changing from $-\infty$ to $\infty$ (so  cusps are not vertices).

The classical 4-vertex theorem asserts that a simple closed convex curve has at least four vertices.  Let $\Delta$ be a cooriented wave front whose Legendrian lift to $ST^*{\mathbb R}^2$  is embedded, i.e., is a {\em Legendrian knot}. V. Arnold conjectured \cite{Ar,Ar1} that if this Legendrian knot is  homotopic as a Legendrian knot to the Legendrian lift of a circle, then $\Delta$ has at least four vertices. This conjecture was proved by Chekanov and Pushkar \cite{ChP} using Legendrian knot theory.

A generic regular homotopy of a cooriented wave front is a composition of a number of moves, similar to the Reidemeister moves in knot theory, see Figure \ref{moves}, borrowed from \cite{ChP}. The first five moves are isotopies of the respective Legendrian knot, but the ``dangerous" self-tangency with coinciding coorientations correspond to self-intersection of the Legendrian lifted curve and changes the Legendrian knot type.

\begin{figure}[ht]
\centering
\includegraphics[width=\textwidth]{./moves}
\caption{Generic ``perestroikas" of cooriented wave fronts.
}
\label{moves}
\end{figure}

For example, the curve on the left of Figure \ref{curves} has only two vertices, but the curve on the right is Legendrian isotopic to a circle, therefore  it has at least four vertices no matter how one draws it. Thus these curves are not Legendrian isotopic. On the other hand,  the Whitney winding number of both curves is one, hence they are regularly isotopic.

\begin{figure}[ht]
\centering
\includegraphics[height=.23\textwidth]{./curve1}
\hspace{.2\textwidth}
\includegraphics[height=.23\textwidth]{./curve2}
\caption{Left: only two vertices; right: at least four vertices.
}
\label{curves}
\end{figure}

\subsection{Summary}
With each ``beam'' of light rays   (a 1-parameter family of oriented lines in ${\mathbb R}^2$) we have associated four curves,
$$ C\subset \L,\qquad  \widetilde C\subset ST^*{\mathbb R}^2,\qquad \Delta\subset{\mathbb R}^2,
\qquad\Gamma\subset {\mathbb {RP}}^2.
$$
The correspondences between the cusps, vertices and inflection points on these curves is depicted in Figure \ref{fig:co}.
\begin{figure}[ht]
\centering
\figsize{.6}{diagram1}
\caption
}
\label{fig:co}
\end{figure}

\section{Proofs of Theorem \ref{thm:main}}

\paragraph{First proof.}  Cusps of $\Gamma_n$ correspond, by projective duality, to inflection points of $C_n$. These  inflections points  are  3-point contacts of $C_n$ with the curves in $\L$ corresponding to the 1-parameter families of lines passing through a fixed point. If the point is $(a,b)\in{\mathbb R}^2$, then the respective curve in $\L$ is
the graph of the first harmonic
$$
p=a\sin \alpha - b \cos \alpha,
$$
that is, it is an ellipse obtained as the intersection of the cylinder $\L$ with a plane through the origin. Note that this graph encloses zero signed area.

Consider the central projection of this cylinder to the unit sphere. This projection sends the graphs of the first harmonics to great circles.

Since $C_n$ encloses zero signed area, it intersects every graph of the first harmonics. It follows that $\bar C_n$, the image of the curve $C_n$, is a smooth spherical curve that is not contained in any hemisphere. In particular, the convex hull of $\bar C_n$ contains the origin.

The (geodesic)  inflections of $\bar C_n$ in the standard metric of the sphere are its 3-point contacts with great circles. By the Segre theorem mentioned earlier, $\bar C_n$ has at least four inflections. Therefore so does $C_n$. \qed

\paragraph{Second proof.} Following \cite{An}, one can use the curve shortening flow to prove that the curve $C_n$ has at least four inflections.
Recall that under the curve shortening flow, each point of the curve moves in the normal direction with the speed equal to the curvature; see
\cite{GH,Gr} and the book \cite{CZ}.

Equip $\L$ with the flat Riemannian metric $d\alpha^2+dp^2$ and apply the curve shortening flow to $C_n$. Let $C_n(s)$ be arclength parametrization, then the flow is given by the partial differential equation $C_t=C_{ss}$.

A variation of the standard proof shows that the evolution is defined for all $t\geq 0$, deforming $C_n$  through embedded curves,  shrinking it to a horizontal curve $p=const$, which is a closed geodesic.
A version of the maximum principle implies that the number of inflections does not increase during this evolution, see \cite{An}. This is illustrated in Figure \ref{infl}.

\begin{figure}[ht]
\centering
\includegraphics[height=.35\textwidth]{./flow_vert.pdf}
\caption{In the curve shortening flow, two nearby inflection points may cancel each other, but they cannot appear on a convex arc.
}
\label{infl}
\end{figure}

Next, a version of Lemma 3.1.7 in \cite{GH}  or Lemma 1.10 of \cite{Gr} shows that the evolving curves enclose zero signed areas.

\begin{lemma} \label{lm:flow}
The curve shortening flow   $C_t=C_{ss}$  preserves the signed area $\int_C p d\alpha.$ \end{lemma}

\begin{proof}
Let $(\alpha(s),p(s))$ be arc length parameterization of a non-contractible  closed curve in $\L$, so that $\alpha_s^2+p_s^2=1$. Then its time evolution under the curve-shortening flow is given by
$
(\alpha_t,p_t)=kN=k(-p_s,\alpha_s),
$
where the subscript denotes the derivative, $k$ is the curvature, and $N$ is the unit normal. It follows that
$$
\frac{d \int p d\alpha}{dt} = \int [k \alpha_s^2 - p (kp_s)_s] ds = \int [k \alpha_s^2+kp_s^2] ds = \int k(s) ds,
$$
where the second equality is due to integration by parts.

It remains to note that the total curvature of a closed curve that goes around the cylinder equals zero.
\end{proof}

 As $C_n$ approaches  $C_0$, it is given by the graph of a function $p=F(\alpha)$. The inflection points of this graph are the points where it is tangent to 2nd order to  the graphs of functions of the form $h(\alpha)=a\cos(\alpha)+b\sin(\alpha) .$

 For each $\alpha\in S^1$, one can find unique $a,b$ such that
 $F(\alpha)=h(\alpha)$, $ F'(\alpha)=h'(\alpha).$ Since $h''+h=0,$ the equation $F''(\alpha)=h''(\alpha) $ holds if and only if  $F''(\alpha)+F(\alpha)=0.$ So the inflection points of $F$ are  the zeros  of the function  $G:=F''+F.$

 To conclude that $G$ has no less than four zeros, apply the Sturm-Hurwitz theorem that states that the number of zeros of a $2\pi$-periodic function is not less than the number of zeros of its first non-trivial harmonic, see, e.g., \cite{Ar,Ar1}.

 Since the differential operator $d^2+1$ preserves the order of Fourier terms and kills 1st order terms, $G$ has no first harmonics. Since the curve encloses zero signed area, $F$ has zero constant term, and so does $G$, as needed. \qed

 \begin{remark}
 {\rm The Sturm-Hurwitz theorem has many proofs, see Section 8.1 of \cite{OT}. Interestingly, one of them, due to G. Polya, makes use of the heat equation, a close relative  of the curve shortening flow.
 }
 \end{remark}

\paragraph{Third proof.}
This argument relies on the correspondence between the cusps of $\Gamma_n$ and the vertices of $\Delta_n$, its normal front.

Since $C_n$ encloses zero signed area, Lemma \ref{lm:front} implies that it admits a Legendrian lift $\widetilde C_n\subset ST^*{\mathbb R}^2$, and its  projection to ${\mathbb R}^2$ is a closed curve, possibly with cusps, which is normal to the rays of $C_n$.

A homotopy of the curve $C_n$ to $C_0$ in the class of smooth closed embedded curves that enclose zero signed area induces a Legendrian isotopy between the Legendrian knots   $\widetilde C_n$ and $\widetilde C_0$. (Such a homotopy is provided by the curve shortening flow but, unlike the second proof, one can use any other homotopy for this purpose).

Now our ``black box", the Pushkar-Chekanov  theorem \cite{ChP}, implies that $\Delta_n$ has at least four vertices. \qed

\section{Variations}

In this last section we briefly mention other variants of our result.

\paragraph{Spherical and hyperbolic geometry.}
One can extend  Theorem \ref{thm:main} to geodesically convex billiards in  spherical and hyperbolic geometries (the former lie in one hemisphere).

The space of oriented geodesics (great circles) in $S^2$ is identified with $S^2$ itself via the usual equator/pole correspondence. The phase cylinder $M\subset S^2$ parametrizes oriented geodesics intersecting $\gamma$. The billiard ball map   preserves the standard area form on $S^2$, and the curve $C_n=T^n(C_0)$ bisects the area of the sphere.

According to Arnold's ``tennis ball theorem" \cite{Ar} (which also follows from the theorem of Segre), a smooth closed embedded spherical curve that bisects the area has at least four inflections. As before,  the curve $C_n$ is dual to the caustic $\Gamma_n$, hence the latter has at least four cusps.

In the hyperbolic case, consider  the hyperboloid model  $H^2=\{x^2+y^2-z^2=-1, \ z>0\}$, with the Minkowski metric $dx^2+dy^2-dz^2$ restricted to $H^2$. An oriented geodesic is given by intersecting $H^2$  with  an oriented plane through the origin, the orthogonal complement (in the Minkowski sense) of a space-like unit vector.

These vectors comprise a hyperboloid of one sheet $H^{1,1}=\{x^2+y^2-z^2=1\}$, equipped with the area form induced from the ambient Minkowsky space (this area form is invariant under the group of motions  $SO(2,1)$).
The phase cylinder $M\subset H^{1,1}$ corresponds to geodesics intersecting $\gamma$, and the billiard ball map is exact area preserving.

To show that the curve $C_n$ has at least four inflections, we use the same argument as in our first proof of Theorem \ref{thm:main}: centrally project the hyperboloid to the sphere. This projection takes inflections to inflections, and the image of $C_n$ contains the origin in its convex hull. Then the Segre theorem implies the result.

\paragraph{Projective billiards.} For any convex curve $\gamma\subset{\mathbb R}^2$ with a transverse vector field $v$ along it one can define the projective billiard map $T:M\to M$ \cite{TaE,TaG}. The reflection law is as follows. Consider an incoming ray  at a point $x\in\gamma$ in the  direction $u$, decompose $u=u_1+u_2$,  where $u_1$ is tangent to $\gamma$ at $x$ and $u_2$ is a multiple  of $v(x)$. Then the  outgoing ray passes through $x$ in the direction $u_1-u_2$. Equivalently,  the tangent line, the transverse line, the incoming, and the outgoing ones, form a harmonic quadruple of lines.

If the transverse field consists of the normals, one has the usual  law ``the angle of incidence equals the angle of reflection".

If $\gamma$ is an origin-centered ellipse and the transverse field $v$ is given by the gradient of a homogeneous function of two variables, then the projective billiard ball map is again exact area preserving. The area form on the phase space $M$ is the same as the one on the space of oriented geodesics in the hyperbolic plane, but this time one considers the projective, or Cayley-Klein, model of hyperbolic geometry in the interior of the ellipse $\gamma$. The total area of $M$ is infinite in this case.

As before, Theorem \ref{thm:main} holds: see Figure \ref{fig:proj} for an illustration.

\begin{figure}[ht]
\centering
\includegraphics[width=.7\textwidth]{./proj1_crop}
\caption{The 1st  caustic by  reflection, showing 8 cusps,  in a projective billiard system with a circular table and the  exact transverse field $v= \nabla (x^4+y^4)$.
}
\label{fig:proj}
\end{figure}

\paragraph{Other initial beams.}
Finally, we note that Theorem \ref{thm:main}, along with its proofs, extends to some other initial beams of light. For example, one may consider a 1-dimensional source, an oval that lies inside $\gamma$ and that emanate rays of light in the outward normal directions.

Differential Geometry (math.DG); Dynamical Systems (math.DS); Symplectic Geometry (math.SG)


A new simple proof for Lum-Chua's conjecture

Abstract: The already proved Lum-Chua's conjecture says that a continuous planar piecewise linear differential system with two zones separated by a straight line has at most one limit cycle. In this paper, we provide a new proof by using a novel characterization for Poincaré half-maps in planar linear systems. This proof is very short and straightforward, because this characterization avoids the inherent flaws of the usual methods to study piecewise linear systems (the appearance of large case-by-case analysis due to the unnecessary discrimination between the different spectra of the involved matrices, to deal with transcendental equations forced by the implicit occurrence of flight time, ...). In addition, the application of the characterization allow us to prove that if a limit cycle exists, then it is hyperbolic and its stability is determined by a simple relationship between the parameters. To the best of our knowledge, the hyperbolicity of the limit cycle and this simple expression for its stability have not been pointed out before.

\section{Introduction}
\label{intro}
The study of limit cycles in planar piecewise linear differential systems dates back to Andronov et al. \cite{AndronovEtAl66} in 1937.  Since then, these systems have received a lot of attention by the scientific community mainly because of their wide range of application in applied science as idealization of nonlinear phenomenon (see, for instance, \cite{CarmonaEtAl02,FreireEtAl98,LlibreEtAl13} and the references therein).

The following continuous planar piecewise linear differential system with two zones separated by a straight line is the simplest possible configuration for a piecewise linear differential system,
\begin{equation}\label{s1}
\dot \mathbf{x} =
\left\{\begin{array}{l}
A_L\mathbf{x}+\mathbf{b}, \quad\textrm{if}\quad x_1\leqslant 0,\\
A_R\mathbf{x}+\mathbf{b}, \quad\textrm{if}\quad x_1\geqslant 0,
\end{array}\right.
\end{equation}
where $\mathbf{x}=(x_1,x_2)^T\in\ensuremath{\mathbb{R}}^2,$ $A_{L}=(a_{ij}^{L})_{2\times 2},$ $A_{R}=(a_{ij}^{R})_{2\times 2},$ with  $a_{12}^L=a_{12}^R=a_{12}$ and $a_{22}^L=a_{22}^R=a_{22},$ and $\mathbf{b}=(b_1,b_2)^T\in\ensuremath{\mathbb{R}}^2.$

A limit cycle, in the context of continuous planar piecewise linear systems, is a non-trivial closed orbit isolated of other closed orbits (see, for instance, Definition $1.6$ in work \cite{LumChua91}). In 1991, after computer experimentations, Lum and Chua stated the following conjecture:

\medskip

\noindent{\bf Lum--Chua's Conjecture.} (\cite{LumChua91}) {\it A continuous planar piecewise differential system with two zones separated by a straight line has at most one limit cycle. The limit cycle, if it exists, is either attracting or repelling.}

\medskip

This conjecture was proven in 1998 by Freire et al. \cite{FreireEtAl98}.  Their proof is performed on a large case-by-case study, which requires individual reasonings and techniques for every possible scenario depending on
the different combinations of configurations (center, focus, saddle, node, \ldots) given by the spectra of matrices $A_L$ and $A_R$. Moreover, their approach forces the appearance of the flight-time as a new implicit variable that hinders the study.
In 2013, Llibre et al. \cite{LlibreEtAl13} made use of Massera's approach \cite{Massera54} for proving a particular case of this conjecture, namely when the determinants
satisfy $\det(A_L)>0$ and $\det(A_R)\geqslant0$.

Our main goal in this paper is to provide a new and simple proof for Lum--Chua's conjecture. Our proof is based on the novel integral characterization for Poincar\'e half-maps for planar linear differential systems presented in \cite{CarmonaEtAl19}, which provides a unified way to deal with the problem, avoiding the large case-by-case study of the former proof and the unnecessary appearance of the superfluous flight-time. In addition, we prove that the limit cycle, if it exists, is hyperbolic and, consequently, either attracting or repelling.  Moreover,  the stability of the limit cycle is also given in terms of the sign of a simple relationship between the parameters (see Theorem \ref{main}). As far as we know, the hyperbolicity of the limit cycle and the easy characterization of its stability have not previously been stated in the literature.

Accordingly, Lum--Chua's Conjecture follows straightforwardly from the next theorem, which is the main result of this paper.

\begin{theorem}\label{main}
The continuous planar piecewise differential system \eqref{s1} has at most one limit cycle, which is hyperbolic if it exists. Moreover, in this case, $(a_{12}b_2-a_{22}b_1) \mathrm{tr}(A_{L})\neq0$ and the limit cycle is attracting $($resp. repelling$)$ provided that $(a_{12}b_2-a_{22}b_1) \mathrm{tr}(A_{L})<0$ $($resp. $(a_{12}b_2-a_{22}b_1) \mathrm{tr}(A_{L})>0),$ where $\mathrm{tr}$ stands for the trace of the matrix.
\end{theorem}

This paper is structured as follows. First, in Section \ref{sec:NFandIIF} we introduce the preliminary results needed to prove Theorem \ref{main}. More specifically, in Section \ref{sec:NF}, we introduce the {\it Li\'{e}nard Canonical Form} for continuous piecewise linear differential systems and, in Section \ref{sec:IIF}, we  present
a result, given in \cite{CarmonaEtAl19}, that allows to manage the graphs of Poincar\'{e} half-maps of piecewise linear differential systems as orbits of a specific cubic vector field. Section \ref{sec:proof} is completely devoted to the proof of Theorem \ref{main}.
 In order to illustrate the concepts and results shown in this paper, we include and comment, in Section \ref{sec:example}, several figures that correspond to a particular system. The choice of parameters for this system is not relevant to the ideas or conclusions developed in the paper and it is obvious that any other values of the parameters would also give scenarios consistent with the proved results. A last small section is dedicated to give some conclusions.

\section{Li\'enard Canonical form and Poincar\'{e} Half-Maps}
\label{sec:NFandIIF}
This section is devoted to introduce the preliminary results needed in the proof of Theorem \ref{main}.

\subsection{Li\'{e}nard Canonical Form} \label{sec:NF}
One can readily see that the assumption $a_{12}\neq0$ is a necessary condition for the existence of periodic solutions of system \eqref{s1}. In this case, from \cite{CarmonaEtAl02}, the linear change of variables $(x,y)=(x_1,a_{22} x_1-a_{12} x_2-b_1)$ transforms system \eqref{s1} into the following Li\'{e}nard canonical form,
\begin{equation}\label{cf}
(S_L) \left\{\begin{array}{l}
\dot x= T_L x-y\\
\dot y= D_L x-a
\end{array}\right.\quad \text{for}\quad x< 0,
\quad
(S_R) \left\{\begin{array}{l}
\dot x= T_R x-y\\
\dot y= D_R x-a
\end{array}\right.\quad \text{for}\quad x\geqslant 0,
\end{equation}
where $T_{L}=\mathrm{tr}(A_{L}),$ $T_{R}=\mathrm{tr}(A_{R}),$ $D_{L}=\det(A_{L}),$ $D_{R}=\det(A_{R}),$ and  $a=a_{12}b_2-a_{22}b_1.$ Notice that any limit cycle of system \eqref{cf} is anti-clockwise oriented and crosses the switching set
 $\Sigma=\{(x,y)\in\mathbb{R}^2: x=0\}$ twice.

\subsection{Poincar\'{e} Half-Maps }\label{sec:IIF}

Consider the following linear differential system
\begin{equation}\label{lienard}
\left\{\begin{array}{rl}
\dot x=&T x-y,\\
\dot y=&D x-a,
\end{array}\right.
\end{equation}
and take section
$\Sigma=\{(x,y)\in\mathbb{R}^2: x=0\}$.

Roughly speaking, the {\it Forward Poincar\'{e} Half-Map} $y_F$ of system \eqref{lienard} associated to section
$\Sigma$ maps a value $y_0\in[0,+\infty)$ to a value $y_1\in(-\infty,0]$, where $(0,y_1)$
is the first return to $\Sigma$ of the forward trajectory of system \eqref{lienard} starting at $(0,y_0).$ Analogously, we define the {\it Backward Poincar\'{e} Half-Map} $y_B$ for the backward trajectory. For a rigorous definition of both Poincar\'{e} half-maps, their domains and ranges, see \cite{CarmonaEtAl19}.

In the next section we see that $a\neq0$ is a necessary condition for the existence of limit cycles for piecewise linear system \eqref{cf}.
Under this condition, in \cite{CarmonaEtAl19}, it is proved that the graphs of the Poincar\'{e} half-map $y_1=y_F(y_0)$ and $y_1=y_B(y_0)$  are orbits of the following cubic vector field
\begin{equation}\label{dy1}
X(y_0,y_1)=-\big(y_1 \big(Dy_0^2-aTy_0+a^2\big),y_0 \big(Dy_1^2-aTy_1+a^2\big)\big),
\end{equation}
with the same orientation as $y_0>0$ increases.
Moreover, it is proved that the existence of the forward (resp. backward) Poincar\'{e} half-map for any value $y_0$ implies that the quadratic function $W(y)=Dy^2-aTy+a^2$ is strictly positive in the interval $[y_F(y_0),y_0]$ (resp. $[y_B(y_0),y_0]$). In particular, it is strictly positive for the domains and ranges of $y_F$ and $y_B$.

\section{Proof of Theorem \ref{main}}\label{sec:proof}
As discussed in Section \ref{sec:NF}, system \eqref{s1} is transformed into system \eqref{cf} through a linear change of variables, provided that $a_{12}\neq0,$ which is a necessary condition for the existence of periodic solutions.

Since systems $(S_L)$ and $(S_R)$, given in  \eqref{cf}, are linear, limit cycles of system \eqref{cf} cannot be completely contained in the closed half-planes $\{(x,y)\in\mathbb{R}^2: x\leqslant 0\}$ or $\{(x,y)\in\mathbb{R}^2: x\geqslant 0\}$. Then, it is a simple consequence of the Green's Theorem that no periodic orbit exists when $T_L T_R>0$ (see \cite{FreireEtAl98}). It is also easy to see that no limit cycle can exist when either the system is homogeneous, i.e. $a=0,$ or $T_L T_R=0.$ Thus, for the sake of our interest, it is sufficient to assume that $a\neq0$ and $T_L T_R<0.$

Let $y_{L}$ (resp. $y_{R}$) be the forward (resp. backward) Poincar\'{e} half-map of planar system $(S_L)$ (resp. $(S_R)$)
associated to section $\Sigma$,
and
let $I_{L}\subset [0,+\infty)$ (resp. $I_{R}\subset [0,+\infty)$) be its interval of definition.
Obviously,  no periodic solution can exist when $I_L\cap I_R=\emptyset.$ Thus, for $y_0\in I\defeq I_L\cap I_R,$ define the displacement function $\delta(y_0)= y_R(y_0)-y_L(y_0).$ Clearly, $\delta(y_0^*)=0$ if, and only if, there exists a periodic orbit passing through $(0,y_0^*)$ and $(0,y_1^*),$ $y_1^*=y_R(y_0^*)=y_L(y_0^*)$. Notice that if that periodic orbit is a limit cycle then it must be $y_0^*>0$ and $y_1^*<0$. Furthermore, it is a hyperbolic limit cycle if, and only if, $\delta'(y_0^*)\neq0.$ In this case, the limit cycle is attracting (resp. repelling) provided that $\delta'(y_0^*)<0$ (resp. $\delta'(y_0^*)>0$).

From now on, let us assume that $y_0^*\in I$, $y_0^*>0$, $\delta(y_0^*)=0,$ and $y_1^*=y_R(y_0^*)=y_L(y_0^*)<0$, i.e., there exist a periodic orbit of system \eqref{cf} that crosses transversally section $\Sigma$ in two points $(0,y_0^*)$ and $(0,y_1^*)$. Next, we are going to prove that the sign of the first derivative of the displacement function $\delta$ in $y_0^*$ is given by the sign of a quadratic function at point $(y_0^*,y_1^*)$.
From \eqref{dy1}, one gets
\begin{equation}\label{diffeq}
\dfrac{d y_{L}}{d y_0}(y_0^*)=\dfrac{y_0^* W_{L}(y_1^*)}{y_1^* W_{L}(y_0^*)}
\qquad \mbox{and} \qquad
\dfrac{d y_{R}}{d y_0}(y_0^*)=\dfrac{y_0^* W_{R}(y_1^*)}{y_1^* W_{R}(y_0^*)},
\end{equation}
being
\begin{equation}\label{iif}
\begin{array}{l}
W_L(y)=D_Ly^2-aT_Ly+a^2 \qquad \mbox{and} \qquad W_R(y)=D_Ry^2-aT_Ry+a^2.
\end{array}
\end{equation}
As mentioned in Section \ref{sec:IIF}, function $W_L$ (resp. $W_R$) must be strictly positive for every interval $[y_L(y_0),y_0] $ (resp. $[y_R(y_0),y_0] $) with $y_0\in I$.

From \eqref{diffeq}, one has
$\delta'(y_0^*)=C(y_0^*,y_1^*)F(y_0^*,y_1^*),$
where functions $C$ and $F$ are real functions defined by
\[
C(y_0,y_1)=\dfrac{-y_0(y_0-y_1)}{y_1W_R(y_0)W_L(y_0)}
\]
and
\begin{equation}\label{Fe1}
F(y_0,y_1)=\dfrac{W_L(y_1)W_R(y_0)-W_L(y_0)W_R(y_1)}{y_0-y_1}.
\end{equation}
Since $W_R(y_0^*)>0$ and $W_L(y_0^*)>0$,
then the sign of $\delta'(y_0^*)$ is determined by the sign of $F.$
Substituting \eqref{iif} into $F(y_0,y_1),$ we get
\begin{equation*}\label{Fe2}
F(y_0,y_1)=a^3(T_L-T_R)+a(D_L T_R-D_R T_L)y_0 y_1+a^2(D_R-D_L)(y_0+y_1),
\end{equation*}
which is a quadratic function. Moreover, curve $\gamma=\{(y_0,y_1)\in\ensuremath{\mathbb{R}}^2:F(y_0,y_1)=0\}$
describes a hyperbola, possibly degenerate. Denote
\begin{equation}\label{Q}
\begin{array}{l}
Q=\{(y_0,y_1)\in\ensuremath{\mathbb{R}}^2:\,
y_1\leqslant 0 \leqslant y_0\}.
\end{array}
\end{equation}
One can easily see that the curve $\gamma$ splits the set $Q\setminus\gamma$ into two disjoint connected set, $Q\setminus\gamma=R_+\cup R_-,$ where
\begin{equation}\label{Rpm}
R_{\pm}=\{(y_0,y_1)\in Q:\, \mathrm{sign}(F(y_0,y_1))=\pm\mathrm{sign}(aT_L)\}.
\end{equation}
Notice that $R_-$ could be the empty set, nevertheless $(0,0)\in R_+$ because $F(0,0)=a^3(T_L-T_R)$ and $T_L T_R<0$.

Now, we are going to see that  $(y_0^*,y_1^*)\in R_+$ and, therefore, the corresponding periodic orbit is a hyperbolic limit cycle whose stability is given by $\mathrm{sign}(aT_L)\neq0$. Notice that this implies the uniqueness of the limit cycle.

If $R_-=\emptyset,$ then it is trivial that $(y_0^{*},y_1^{*})\in R_+,$ so let us assume that $R_-$ is not empty.

Let us consider the graphs
$\gamma_{L}=\{(y_0,y_{L}(y_0)): \, y_0\in I_{L}\}$ and $\gamma_{R}=\{(y_0,y_{R}(y_0)): \, y_0\in I_{R}\}$, which are contained in $Q.$ Clearly, point $(y_0^*,y_1^*) \in\gamma_L\cap\gamma_R$.
From \eqref{dy1}, $\gamma_{L}$ and $\gamma_{R}$ are orbits of the following cubic vector fields
\begin{equation}\label{X}
X_{L}(y_0,y_1)=-\big(y_1W_{L}(y_0),y_0W_{L}(y_1)\big) \ \, \mbox{and}\ \,
X_{R}(y_0,y_1)=-\big(y_1W_{R}(y_0),y_0W_{R}(y_1)\big),
\end{equation}
respectively.
 From now on, orbits $\gamma_{L}$ and $\gamma_{R}$ will be endowed with the orientation given by the corresponding vector field. These objects are illustrated in Fig.~\ref{fields} of next section for a particular example.

In the sequel, we shall study the vector fields $X_{L}$ and $X_{R}$ along the curve $\gamma$ for $(y_0,y_1)\in\textrm{int}(Q).$ From \eqref{Fe1},  if $(y_0,y_1)\in\textrm{int}(Q),$ the equation $F(y_0,y_1)=0$ is equivalent to equality $W_L(y_1)W_R(y_0)=W_L(y_0)W_R(y_1).$
Thus, substituting this last equality into $\langle\nabla F(y_0,y_1),X_{L}(y_0,y_1)\rangle$ and $\langle\nabla F(y_0,y_1),X_{R}(y_0,y_1)\rangle$ and using expression \eqref{Fe1} of $F$ we get
\begin{equation*}\label{GL}
\begin{array}{rl}
G_{L}(y_0,y_1)=&\big\langle\nabla F(y_0,y_1),X_{L}(y_0,y_1)\big\rangle\big|_{(y_0,y_1)\in \gamma
}\vspace{0.2cm}\\
=&W_{L}(y_1)a\big(T_L W_R(y_0)-T_R W_L(y_0)\big)
\end{array}
\end{equation*}
and
\begin{equation*}\label{GR}
\begin{array}{rl}
G_{R}(y_0,y_1)=&\big\langle\nabla F(y_0,y_1),X_{R}(y_0,y_1)\big\rangle\big|_{(y_0,y_1)\in\gamma
}\vspace{0.2cm}\\
=&W_{R}(y_1)a\big(T_L W_R(y_0)-T_R W_L(y_0)\big),
\end{array}
\end{equation*}
respectively.
Since $T_L T_R<0,$ we conclude that
\begin{equation}\label{sgnGRL}
\mathrm{sign}(G_{L}(y_0,y_1))=\mathrm{sign}(G_{R}(y_0,y_1))=\mathrm{sign}(a T_L) \qquad \forall \ (y_0,y_1)\in\gamma.
\end{equation}
This means that each
curve $\gamma_{L}$ and $\gamma_{R}$
intersects $\gamma$ at most once  in $\textrm{int}(Q)$ and,  if this intersection exists, it must be crossing $\gamma$ from $R_-$ to $R_+$. This behavior is illustrated In Fig.~\ref{final} of the next section.

Now, since $a\neq0,$ the origin is a quadratic contact point of the continuous piecewise linear system \eqref{cf}. Thus, at least one of the Poincar\'{e} half-maps  is defined for $y_0> 0$ sufficiently small and can be continuously extended to $y_0=0$ with image $y_1=0$, see Remark 1 of  \cite{CarmonaEtAl19} (in the example considered in the next section, this is the case of the backward Poincar\'e half-map $y_R$, see Fig.~\ref{phase_planes}(b) and  Fig.~\ref{fields}(b)).
Consequently, the graph of such a Poincar\'{e} half-map contains the origin and, since $(0,0)\in R_+$, it cannot intersect the set $\gamma$ (see Fig.~\ref{final}). Hence, point $(y_0^*,y_1^*)\in \textrm{int}(Q)$  is contained in $R_+,$ that is $F(y_0^*,y_1^*)\neq0$ and $\mathrm{sign}(F(y_0^*,y_1^*))=\mathrm{sign}(a T_L).$ This implies that the corresponding periodic orbit is a hyperbolic limit cycle and its stability is determined by $\mathrm{sign}(aT_L),$ namely it is attracting (resp. repelling) provided that $a T_L<0$ (resp. $a T_L>0$).

The proof is concluded.

\section{Illustrative example}
\label{sec:example}

Let us consider the continuous planar piecewise linear system \eqref{cf} for values
\begin{equation}\label{values}
T_L=0.4, \ T_R=-0.3, \ D_L=3, \ D_R=0.1, \ \text{and} \ a=-0.2.
\end{equation}
Let us also consider the corresponding linear systems $(S_L)$ and $(S_R)$ involved in system \eqref{cf} for the same values given in \eqref{values}.
As it was said in last paragraph of Section \ref{intro}, the values of parameters given in \eqref{values} have been fixed just for the sake of illustrating, with nice and simple figures, some issues mentioned in the previous sections. Obviously, many other choices could also
be taken with the same purpose.

The phase planes of linear systems $(S_L)$ and $(S_R)$ are shown in Fig.~\ref{phase_planes}. Solid curves are real pieces of orbits of piecewise system \eqref{cf} but dotted ones are pieces of orbits of the linear systems that do not correspond to orbits of \eqref{cf}. The orbits of piecewise system \eqref{cf} are obtained by suitable concatenations of the solid curves.

Both linear systems correspond to focus configurations. However, while the equilibrium of system $(S_L)$ is a real unstable equilibrium of system \eqref{cf}, the stable equilibrium of system $(S_R)$ is not an equilibrium of system \eqref{cf}, because it is located in the open half-plane $\{(x,y)\in\mathbb{R}: x<0\}$. Besides that,
the domain of the forward Poincar\'e half-map $y_L$ is $[0,+\infty)$
and its range is $(-\infty,\bar{y}_1]$, where $\bar{y}_1=y_L(0)$
(see Fig. \ref{phase_planes}(a)), while the domain of the
backward Poincar\'e half-map $y_R$ is $[0,+\infty)$ and its range is $(-\infty,0]$,  taking into account that it can be continuously extended to the origin as $y_R(0)=0$ (see Fig. \ref{phase_planes}(b)). Note that polynomials $W_L(y)=D_L y^2-aT_L y+a^2$ and $W_R(y)=D_R y^2-aT_R y+a^2$ are strictly positive for every $y\in\mathbb{R}$. In particular, as mentioned  in Section \ref{sec:IIF}, they are strictly positive for the domains and ranges of both Poincar\'e half-maps, as stated in \cite{CarmonaEtAl19}.

\begin{figure}
\[
\begin{array}{cc}
\includegraphics[width=0.45\linewidth]{focoizquierda.eps}&
\includegraphics[width=0.45\linewidth]{focoderecha.eps}\\
\text{(a)}&\text{(b)}
\end{array}
\]
\caption{\label{phase_planes}
Phase planes of the linear systems involved in \eqref{cf} for values \eqref{values}: (a) system $(S_L)$, (b) system $(S_R)$.
Solid curves are real pieces of orbits of piecewise system \eqref{cf}, while dotted ones stand for pieces of orbits of the linear systems $(S_L)$ and $(S_R)$ that do not correspond to orbits of system \eqref{cf}.
Note that, in (a), the interval $(\bar{y}_1,0)$, where $\bar{y}_1=y_L(0)$, is not included in the range of $y_L$. However, in (b), the backward Poincar\'e half-map $y_R$ can be continuously extended to the origin as $y_R(0)=0$, due to the quadratic contact.}
\end{figure}

In Fig.~\ref{fields}, vector fields  $X_L$ and $X_R$, together with some of their orbits have been shown  in the fourth quadrant $Q$ given in \eqref{Q} (for the sake of clarity, vectors have been normalized in the figure). Among all the curves, the ones corresponding to the forward and backward Poincar\'e half-maps have been highlighted by their thickness. Taking into account the direction induced by the flow, we see that curve $\gamma_L$  begins at point $(0,\bar{y}_1)$ and $\gamma_R$  begins at the origin. This agrees with the previous descriptions of the domains and ranges of $y_L$ and $y_R$ obtained from Fig.~\ref{phase_planes}.

\begin{figure}
\[
\begin{array}{cc}
\includegraphics[width=0.45\linewidth]{izquierda.eps}&
\includegraphics[width=0.45\linewidth]{derecha.eps}\\
\text{(a)}&\text{(b)}
\end{array}
\]
\caption{\label{fields} (a) Sketch of vector field $X_L$ given in equation \eqref{X} and several orbits of its flow (note that the vectors have been normalized for a clearer view of their directions). The thickest curve corresponds to the graph $\gamma_L$ of the forward Poincar\'e half-map $y_L$. This graph contains point $(0,\bar{y}_1)$, where $\bar{y}_1=y_L(0)$.
(b)  Idem for vector field $X_R$ given in equation \eqref{X}. Note that, due to the second order contact, the origin lies in the graph $\gamma_R$ of the backward Poincar\'e half-map $y_R$.}
\end{figure}

A limit cycle for system \eqref{cf} with parameters given in \eqref{values} exists if  curves $\gamma_L$  and $\gamma_R$ have an isolated intersection point in the open set $\textrm{int}(Q)$.
Since in Fig.~\ref{final} we can observe such an intersection point, then the system has a limit cycle.

Some other important elements for the theoretical analysis developed in the previous section have been added in Fig.~\ref{final}. In fact, the dashed curve corresponds to hyperbola $\gamma$, given by equation $F(y_0,y_1)=0$  for the values of the parameters fixed in \eqref{values}, being $F$ the function defined in \eqref{Fe1}. The vector fields $X_L$ and $X_R$ (once they have been normalized) are sketched on curve $\gamma$ in order to show that they point to $R_+$ (defined in equation \eqref{Rpm}). Notice that vector fields $X_L$ and $X_R$ have the same direction and sense over $\gamma$ (this is a direct conclusion from the equalities $F(y_0,y_1)=0$ and \eqref{sgnGRL}). Therefore, their normalizations are equal over $\gamma$.

In this example, we can check that, as was theoretically deduced in previous section, at least one of the curves of the Poincar\'e half-maps ($\gamma_R$ in this case) is always contained in region $R_+$ because it begins at the origin and cannot cross curve $\gamma$. On the contrary, curve $\gamma_L$ begins from a point of region $R_-$ and must pass through $\gamma$ in order to intersect $\gamma_R$. Since a new intersection of $\gamma_L$ and $\gamma_R$ would require that both curves pass through curve $\gamma$ once more, which is impossible, the maximum number of limit cycles must be one.

\begin{figure}
\includegraphics[width=0.45\linewidth]{final.eps}\\
\caption{\label{final} Diagram, in the fourth quadrant $Q$, of the relative position between the graphs of the Poincar\'e half-maps $y_L$ and $y_R$ (that contains the origin) and the curve $\gamma=\{(y_0,y_1)\in Q: F(y_0,y_1)=0\}$, where $F$ is defined in \eqref{Fe1}. The graph of $y_L$ contains point $(0,\bar{y}_1)$, while the origin lies in the graph of $y_R$. Vector fields $X_L$ and $X_R$, normalized like in previous figures, are sketched over the curve $\gamma$ (dashed). Actually, both normalized vector fields coincide on $\gamma$. Note that the orbits of $X_L$ and $X_R$ can cross curve $\gamma$ from $R_-$ to $R_+$. }
\end{figure}

\section{Some conclusions and further directions}
\label{conclusion}
It has just been observed that the proof of Lum--Chua's conjecture can be significantly shorten and simplified, and additionally unified, by using the novel characterization of Poincar\'e half-maps given in \cite{CarmonaEtAl19}.

This fact motivates us to analyze other problems of piecewise linear systems by using the ideas, concepts, and techniques developed in this work.

On the one hand, we would like to look for better ways to solve some previously worked problems, those involving closing equations and Poincar\'e maps, with the additional purpose of developing a common method of study.
On the other hand, this approach allows us to consider open problems whose study with traditional techniques have turned out to be very hard to deal with.

Among the open problems, we can mention the uniqueness of limit cycles in planar piecewise linear differential systems without sliding region. This problem has been addressed by several authors, usually via large case-by-case analyses (see, for instance, \cite{FreireEtAl2013, HuanYang2013, HuanYang2014, LiLLibre2020, LlibreEtAl2008,  MedradoTorregrosa15}), some of them involving transcendental implicit equations due to the flight time.
However, in spite of the previous efforts, as far as we know, the uniqueness of limit cycles in this case has remained an open problem. Currently, we are working in this problem and the present technique has also been proved to be very effective  \cite{CarmonaEtAl20a}.

We can also mention the problem of the existence of a finite upper bound for the maximum number of limit cycles in discontinuous planar piecewise linear differential systems. One can find some partial results in the literature for some non generic families of piecewise linear differential systems (see, for instance, \cite{FreireEtAl12,LlibreEtAl15, Novaes17}). The technique and ideas developed in the present paper allowed us to easily obtain a finite upper bound, which we are currently  improving.

\section*{Acknowledgements}

VC and FFS are partially supported by the Ministerio de Ciencia e Innovación, Plan Nacional I+D+I cofinanced with FEDER funds, in the frame of the projects MTM2014-56272-C2-1-P, MTM2015-65608-P, MTM2017-87915-C2-1-P,  PGC2018-096265-B-I00 and by the Junta de Andaluc\'{i}a grants TIC-0130, P12-FQM-1658. DDN is partially supported by Funda\c{c}\~{a}o de Amparo \`{a} Pesquisa do Estado de S\~{a}o Paulo (FAPESP) grants 2018/16430-8, 2018/ 13481-0, and 2019/10269-3, and by Conselho Nacional de Desenvolvimento Cient\'{i}fico e Tecnol\'{ o}gico CNPq grants 306649/2018-7 and 438975/ 2018-9.