Converse theorem on a contraction metric for a periodic orbit

Abstract: Contraction analysis uses a local criterion to prove the long-term behaviour of a dynamical system. A contraction metric is a Riemannian metric with respect to which the distance between adjacent solutions contracts. If adjacent solutions in all directions perpendicular to the flow are contracted, then there exists a unique periodic orbit, which is exponentially stable and we obtain a bound on the rate of exponential attraction. In this paper we study the converse question and show that, given an exponentially stable periodic orbit, a contraction metric exists on its basin of attraction and we can recover the bound on the rate of exponential attraction.

\section{Introduction}

The stability and basin of attraction of periodic orbits is an important problem in many applications. Already the determination of a periodic orbit is a non-trivial task as it involves solving the differential equation. The classical definition of stability, as well as its study using a Lyapunov function require the knowledge of the position of the periodic orbit which in many applications can only be approximated. An alternative way to study the stability and basin of attraction is contraction analysis, which is a local criterion and does not require us to know the location of the periodic orbit.

Throughout the paper we will study the autonomous ODE
\begin{eqnarray}
\dot{{\bf x}}&=&{\bf f}({\bf x})\label{ODE}
\end{eqnarray}
where ${\bf f}\in C^\sigma(\mathbb R^n,\mathbb R^n)$ with $\sigma\ge 1$. We denote the solution ${\bf x}(t)$ with initial condition ${\bf x}(0)={\bf x}_0$ by $S_t{\bf x}_0={\bf x}(t)$ and  assume that it exist for all $t\ge 0$.

In the next definition we will define a contraction metric on $\mathbb R^n$.
Note that $M({\bf x})$ defines a point-dependent scalar product through $\langle {\bf v},{\bf w}\rangle_{M({\bf x})}={\bf v}^TM({\bf x}){\bf w}$ for all ${\bf v},{\bf w}\in\mathbb R^n$.

\begin{definition}[Contraction metric]\label{Riemannian}
	A Riemannian metric is a function $M\in C^0(G,\mathbb S^n)$, where $G\subset \mathbb R^n$ is open and $\mathbb S^n$ denotes the symmetric $n\times n$ matrices, $M({\bf x})$ is positive definite for all ${\bf x}\in G$ and  the orbital derivative of $M$ exists for all ${\bf x}\in G$ and is continuous, i.e.
$$M'({\bf x})=\frac{d}{dt} M(S_t{\bf x})\big|_{t=0}$$
exists and is continuous. A sufficient condition for the latter is that $M\in C^1(G,\mathbb S^n)$; then $M_{ij}'({\bf x})=\nabla M_{ij}({\bf x})\cdot {\bf f}({\bf x})$ for all $i,j\in\{1,\ldots,n\}$.

Define \begin{eqnarray}
L_M({\bf x};{\bf v})&:=&\frac{1}{2}{\bf v}^T\left(M({\bf x})D{\bf f}({\bf x})+D{\bf f}({\bf x})^TM({\bf x})+M'({\bf x})\right){\bf v}. \label{LM}
\end{eqnarray}

The Riemannian metric $M$ is called {\bf  contraction metric in $K\subset  G$ with exponent $-\nu<0 $} if $L_M({\bf x})\le -\nu$ for all ${\bf x}\in K$, where
\begin{eqnarray}
L_M({\bf x})&:=&\max_{{\bf v}^TM({\bf x}){\bf v}=1,{\bf v}^T M({\bf x}){\bf f}({\bf x})=0}L_M({\bf x};{\bf v}).\label{L_M}
\end{eqnarray}
\end{definition}

The following theorem shows the implications of the existence of such a contraction metric on a certain set in the phase space.

\begin{theorem}
	Let $\varnothing\not= K \subset \mathbb R^n$ be
	a compact, connected and positively invariant set which contains no equilibrium.
	Let $M$ be a   contraction metric in $K$ with exponent $-\nu<0 $, see Definition \ref{Riemannian}.

Then there exists one and only one periodic orbit $\Omega\subset K$. This periodic orbit
is exponentially asymptotically stable, and the real parts of all Floquet exponents --
except the trivial one -- are less than or equal to $-\nu$.
 Moreover, the basin of attraction $A(\Omega)$ contains $K$.
 \end{theorem}

This theorem goes back to Borg \cite{borg} with $M({\bf x})=I$, and has been extended to a general Riemannian metric \cite{stenstroem}. For more results on contraction analysis for a periodic orbit see \cite{hartman,hartmanbook,leonov3,leonov96}.

 Note that a similar result holds with an equilibrium if the contraction takes place in all directions ${\bf v}$, i.e. if $L_M({\bf x})\le -\nu$ in \eqref{L_M} is replaced by
 ${\mathcal L}_M({\bf x}):=\max_{{\bf v}^TM({\bf x}){\bf v}=1}L_M({\bf x};{\bf v})\le -\nu$. For more references on contraction analysis see \cite{lohmiller}, and for the relation to Finsler-Lyapunov functions see \cite{finsler}.

 Note that $L_M({\bf x})$ is a continuous with respect to ${\bf x}$ and, as we will show in the paper, also locally Lipschitz-continuous. Due to the maximum, however, it is not differentiable in general.

In this paper we are interested in converse results, i.e. given an exponentially stable periodic orbit, does a Riemannian contraction metric as in Definition \ref{Riemannian} exist?
\cite{lohmiller} gives a converse theorem, but here $M(t,{\bf x})$ depends on $t$ and will, in general, become unbounded as $t\to\infty$.  In \cite{giesl04} the existence of such a contraction metric was shown on a given compact subset of $A(\Omega)$, first on the periodic orbit, using Floquet theory, and then on $K$, using a Lyapunov function. The local construction, however, neglected the fact that the Floquet representation of solutions of the first variation equation along the periodic orbit is in general not real, but complex. We will show in this paper, that, by choosing the complex Floquet representation appropriately, the constructed Riemannian metric is real-valued, thus justifying the arguments in \cite{giesl04}. Moreover, we will show the existence of a Riemannian metric on the whole, possibly unbounded basin of attraction by using a new construction. The Riemannian metric will be arbitrarily close to the true rate of exponential attraction. Let us summarize the main result of the paper in the following theorem.

\begin{theorem}\label{main}
	Let $\Omega$ be an exponentially stable periodic orbit of $\dot{{\bf x}}={\bf f}({\bf x})$, let $-\nu$ be the largest real part of all its non-trivial Floquet exponents and ${\bf f}\in C^\sigma(\mathbb R^n,\mathbb R^n)$ with $\sigma\ge 3$.

Then for all $\epsilon\in (0,\nu/2)$ there exists a contraction metric $M\in C^{\sigma-1}( A(\Omega), \mathbb S^{n})$ in $A(\Omega)$ as in Definition \ref{Riemannian}  with exponent $-\nu+\epsilon<0 $, i.e.
\begin{eqnarray}L_M({\bf x})&=&\frac{1}{2}\max_{{\bf v}^TM({\bf x}){\bf v}=1,{\bf v}^T M({\bf x}){\bf f}({\bf x})=0}{\bf v}^T\left(M({\bf x})D{\bf f}({\bf x})+D{\bf f}({\bf x})^TM({\bf x})+M'({\bf x})\right){\bf v}
	\nonumber\\
	&\le& -\nu+\epsilon\label{esti}
\end{eqnarray}
holds for all ${\bf x}\in A(\Omega)$.
\end{theorem}

The metric is constructed in several steps: first on the periodic orbit, then in a neighborhood, and finally in the whole basin of attraction. In the proof, we define a projection of points ${\bf x}$ in a neighborhood of the periodic orbit onto the periodic orbit, namely onto ${\bf p}\in \Omega$, such that  to $({\bf x}-{\bf p} )^TM({\bf p}){\bf f}({\bf p})=0$. This is then used to synchronize the times of solutions through ${\bf x}$ and ${\bf p}$, and to define a time-dependent distance between these solutions, which decreases exponentially.

Let us compare our result with  converse theorems for a contraction metric for an equilibrium. In \cite{converse}, three converse theorems were obtained: Theorem 4.1 constructs a metric on a given compact subset of the basin of attraction (see \cite{giesl04} for the case of a periodic orbit), Theorem 4.2 constructs a metric on the whole basin of attraction (see this paper for the case of a periodic orbit), while Theorem 4.4 constructs a metric as solution of a linear matrix-valued PDE (see \cite{other} for the case of a periodic orbit). The latter construction is beneficial for its computation by solving the PDE, and it also constructs a smooth function; however, the exponential rate of attraction cannot be recovered, which is an advantage of the approach in this paper.

Let us give an overview over the paper: In Section \ref{Floquet} we prove a special Floquet normal form to ensure that the contraction metric that we later construct on the periodic orbit is real-valued. In Section \ref{conv} we prove the main result of the paper, Theorem \ref{main}, showing the existence of a Riemannian metric on the whole basin of attraction. The section also contains Corollary \ref{help_other}, defining a projection onto the periodic orbit and related estimates. In the appendix we prove that $L_M$ is locally Lipschitz-continuous.

\section{Floquet normal form}
\label{Floquet}

Before we consider the Floquet normal form, we will prove a lemma which calculates  $L_M({\bf x})$ for the Riemannian metric $M({\bf x})=e^{2V({\bf x})}N({\bf x})$.

\begin{lemma}\label{lem}
	Let $N\colon \mathbb R^n\to \mathbb S^n$ be a Riemannian metric and $V\colon \mathbb R^n\to \mathbb R$ a continuous and orbitally continuously differentiable function.

	Then $M({\bf x})=e^{2V({\bf x})}N({\bf x})$ is a Riemannian metric and
	\begin{eqnarray*}
		L_M({\bf x})&=&L_N({\bf x})+V'({\bf x}).
	\end{eqnarray*}
\end{lemma}
\begin{proof}
	It is clear that $M({\bf x})$ is a positive definite for all ${\bf x}$ since $e^{2V({\bf x})}>0$.
	We have
	\begin{eqnarray*}
		L_M({\bf x};{\bf v})&=&
		\frac{1}{2}{\bf v}^T\left(M({\bf x})D{\bf f}({\bf x})+D{\bf f}({\bf x})^TM({\bf x})+M'({\bf x})\right){\bf v}\\
		&=&
		\frac{1}{2}{\bf v}^T\bigg(e^{2V({\bf x})}N({\bf x})D{\bf f}({\bf x})+e^{2V({\bf x})}D{\bf f}({\bf x})^TN({\bf x})\\
		&&\hspace{1.2cm}+
		e^{2V({\bf x})}(2V'({\bf x})N({\bf x})+N'({\bf x}))\bigg){\bf v}\\
		&=&
		\frac{1}{2}{\bf w}^T\left(N({\bf x})D{\bf f}({\bf x})+D{\bf f}({\bf x})^TN({\bf x})+
		N'({\bf x})\right){\bf w}
		+{\bf w}^TN({\bf x}){\bf w} \, V'({\bf x})
	\end{eqnarray*}
	with ${\bf w}=e^{V({\bf x})}{\bf v}$, so $L_M({\bf x};{\bf v})=L_N({\bf x}; {\bf w})+{\bf w}^TN({\bf x}){\bf w} \, V'({\bf x})$.
	Thus,
	\begin{eqnarray*}

		L_M({\bf x})&=&\max_{{\bf v}^TM({\bf x}){\bf v}=1,{\bf v}^TM({\bf x}){\bf f}({\bf x})=0}L_M({\bf x};{\bf v})\\
		&=&\max_{{\bf w}^TN({\bf x}){\bf w}=1,{\bf w}^T N({\bf x}){\bf f}({\bf x})=0}\left[L_N({\bf x};{\bf w})
		+{\bf w}^TN({\bf x}){\bf w} V'({\bf x})\right]\\
		&=&L_N({\bf x})+ V'({\bf x}).
	\end{eqnarray*}
	This shows the lemma.
\end{proof}

In order to show later that our constructed Riemannian metric $M$ is real-valued, we will construct a special Floquet normal form in Proposition \ref{prop} such that the matrix in \eqref{real} is real-valued. In Corollary \ref{coro} we will show estimates in the case that \eqref{var} is the first variation equation of a periodic orbit.
The proof of the following proposition is inspired by \cite{chicone}.

\begin{proposition}\label{prop}
	Consider the periodic differential equation
	\begin{eqnarray}\dot{{\bf y}}&=&F(t){\bf y}\label{var}
	\end{eqnarray} where $F\in C^s(\mathbb R,\mathbb R^{n\times n})$ is $T$-periodic, $s\ge 1$ and denote by $\Phi\in C^s(\mathbb R,\mathbb R^{n\times n})$ its principal fundamental matrix solution with $\Phi(0)=I$.

	Then there exists a $T$-periodic function $P\in C^s(\mathbb R,\mathbb C^{n\times n})$ with $P(0)=P(T)=I$ and a  matrix $B\in \mathbb C^{n\times n}$ such that for all $t\in\mathbb R$
	$$\Phi(t)=P(t)e^{Bt}.$$

	Denote by $\lambda_1,\ldots,\lambda_r\in \mathbb R\setminus \{0\}$ the pairwise distinct real eigenvalues and by $\lambda_{r+1},\overline{\lambda_{r+1}},\ldots,$
	$\lambda_{r+c},\overline{\lambda_{r+c}}\in\mathbb C\setminus \mathbb R$ the pairwise distinct pairs of complex conjugate complex eigenvalues of $\Phi(T)$ with algebraic multiplicity $m_j$ of $\lambda_j$.
	For $\epsilon>0$ there exists a non-singular matrix $S\in \mathbb R^{n\times n}$ such that $B=SAS^{-1}$ with $A=\operatorname{blockdiag}(K_1,K_2,\ldots,K_{r+c})$ and $K_j\in \mathbb C^{m_j\times m_j}$ for $j=1,\ldots,r$ and   $K_j\in \mathbb R^{2m_j\times 2m_j}$ for $j=r+1,\ldots,r+c$ as well as
	$$\frac{1}{2}{\bf w}^*(A^*+A){\bf w}\le \sum_{j=1}^{r+c}c_j\sum_{i=1}^{m_j}|w_{i+\sum_{k=1}^{j-1}m_k}|^2
	\text{	for all }{\bf w}\in \mathbb C^n,$$
	where $c_j=\left(\frac{\ln |\lambda_j|}{T}+\epsilon\right)$ if $m_j\ge 2$ and  $c_j=\frac{\ln |\lambda_j|}{T}$ if $m_j=1$.

 Moreover, we have
\begin{eqnarray}
(P^{-1}(t))^*(S^{-1})^*S^{-1}P^{-1}(t)&\in& \mathbb R^{n\times n}\label{real}
\end{eqnarray}
	for all $t\in \mathbb R$.
\end{proposition}
\begin{proof}
	Since $F\in C^s$, we also have $\Phi\in C^s(\mathbb R,\mathbb R^{n\times n})$.
	Noting that $\Psi(t):=\Phi(t+T)$ solves \eqref{var} with $\Psi(0)=\Phi(T)$, we obtain from the uniqueness of solutions that
	\begin{eqnarray}
		\Phi(t+T)&=&\Psi(t)=\Phi(t)\Phi(T)\text{ for all }t\in \mathbb R.\label{uni}
	\end{eqnarray}

	Consider $C:=\Phi(T)\in \mathbb R^{n\times n}$ which is non-singular and hence all eigenvalues of $\Phi(T)$ are non-zero. Let $\epsilon':=\frac{1}{2}\min\left( \frac{\epsilon T}{2},1\right)$ and $S\in \mathbb R^{n\times n}$ be such that
	$S^{-1}CS=:J$ is in real Jordan normal form with the $1$ replaced by $\epsilon'  |\lambda_j|$ for each eigenvalue $\lambda_j$, i.e. $J$ is a block-diagonal matrix with blocks $J_j$ of the form
	$J_j=\left(\begin{array}{lllll}\lambda_j&\epsilon' |\lambda_j|&&&\\
	&\lambda_j&\epsilon' |\lambda_j|&&\\
	&&\ddots&\ddots&\\
	&&&\lambda_j&\epsilon' |\lambda_j|\\
	&&&&\lambda_j\end{array}\right)\in\mathbb R^{m_j\times m_j}$ for real eigenvalues $\lambda_j$ of $C$ and
	$J_j=\left(\begin{array}{ccccccc}\alpha_j&-\beta_j&\epsilon' r_j&&&&\\
	\beta_j&\alpha_j&&\epsilon' r_j&&&\\
	&&\ddots&&\ddots&&\\
&&&\alpha_j&-\beta_j&\epsilon' r_j&\\
	&&&\beta_j&\alpha_j&&\epsilon' r_j\\
	&&&&&\alpha_j&-\beta_j\\
	&&&&&
	\beta_j&\alpha_j\end{array}\right)\in\mathbb R^{2m_j\times 2m_j}$ for each pair of complex eigenvalues $\alpha_j\pm i\beta_j$ of $C$, where $r_j=\sqrt{\alpha_j^2+\beta_j^2}$ and $m_j$ denotes the dimension of the generalized eigenspace of one of them; note we have pairs of complex conjugate eigenvalues since $C$ is real.

	This can be achieved by letting  $S_1\in\mathbb R^{n\times n}$ be an invertible matrix such that
	$S^{-1}_1CS_1$ is the standard real Jordan Normal Form with  $1$ on the super diagonal. Then define $S_2$ to be a matrix of blocks
	$$\operatorname{diag}(1,\epsilon' |\lambda_j|,
	(\epsilon')^2|\lambda_j|^2,\ldots,(\epsilon')^{m_j-1}|\lambda_j|^{m_j-1})$$
	for real $\lambda_j$ and

	$$\operatorname{diag}(1,1,\epsilon' |\lambda_j|,\epsilon' |\lambda_j|, \ldots,(\epsilon')^{m_j-1}|\lambda_j|^{m_j-1},(\epsilon')^{m_j-1}|\lambda_j|^{m_j-1})$$ for a pair of complex conjugate eigenvalues $\lambda_j$ and $\overline{\lambda_j}$. Setting $S=S_1S_2$ yields the result.

	For each of the blocks, we will now construct a matrix $K_j\in\mathbb C^{m_j\times m_j}$ for real eigenvalues $\lambda_j$ and $K_j\in\mathbb R^{2m_j\times 2m_j}$ for each pair of complex eigenvalues $\alpha_j\pm i\beta_j$ such that
	$$e^{K_j T}=J_j,$$
	which shows with $B=SAS^{-1}$, where $A:=\operatorname{blockdiag}(K_1,\ldots,K_r)$,
\begin{eqnarray}
	e^{BT}&=&Se^{AT}S^{-1}=S\operatorname{blockdiag}(e^{K_1T},\ldots,e^{K_rT})S^{-1}\nonumber\\
	&=&SJS^{-1}=C=\Phi(T).\label{*}
	\end{eqnarray}

	We distinguish between three cases:  $\lambda_j$ being real positive, real negative or complex.
	Using the series expansion of $
	\ln (1+x)$ we obtain for a nilpotent matrix $M\in \mathbb R^{n\times n}$
\begin{eqnarray}
	\exp\left(\sum_{k=1}^\infty \frac{(-1)^{k+1}}{k}M^k\right)&=&I+M;\label{log}
	\end{eqnarray}
	note  that the sum is actually finite.

	\vspace{0.3cm}

\noindent	\underline{\bf Case 1: $\lambda_j\in \mathbb R^+$}

	Writing $J_j=\lambda_j(I+\epsilon'  N)$  with the nilpotent matrix $N=
	\left(\begin{array}{llll}0&1&&\\
	&\ddots&\ddots&\\
	&&0&1\\
	&&&0\end{array}\right)\in\mathbb R^{m_j\times m_j}$, we define
	$$K_j=\frac{1}{T}\left( (\ln \lambda_j)I+\sum_{k=1}^{m_j-1} \frac{(-1)^{k+1}}{k} (\epsilon')^kN^k\right)\in \mathbb R^{m_j\times m_j}.$$
Since $I$ and $N$ commute, we have with \eqref{log} and $N^k=0$ for $k\ge m_j$
	\begin{eqnarray*}
		\exp(K_j T)=\lambda_j \left(I+\epsilon' N\right)=J_j.
	\end{eqnarray*}

\vspace{0.5cm}

\noindent	\underline{\bf Case 2: $\lambda_j\in \mathbb R^-$}

 With the nilpotent matrix $N=
 \left(\begin{array}{llll}0&1&&\\
 &\ddots&\ddots&\\
 &&0&1\\
 &&&0\end{array}\right)\in\mathbb R^{m_j\times m_j}$ we	write $J_j=-|\lambda_j| (I-\epsilon' N)$ and define
	$$K_j=\frac{1}{T}\left( (i\pi +\ln |\lambda_j|)I+\sum_{k=1}^{m_j-1}  \frac{(-1)^{k+1}}{k}(-\epsilon')^kN^k\right)\in \mathbb C^{m_j\times m_j}.$$
	Since $I$ and $N$ commute, and $N^k=0$ for $k\ge m_j$ we have with \eqref{log}
	\begin{eqnarray*}
		\exp(K_j T)=-|\lambda_j | \left(I-\epsilon' N\right)=J_j.
	\end{eqnarray*}

	\vspace{0.3cm}

	\noindent
	\underline{\bf Case 3: $\lambda_j=\alpha_j+i\beta_j$ with $\beta_j\not=0$ }

	We only consider one of the two complex conjugate eigenvalues $\lambda_j$ and $\overline{\lambda_j}$ of $\Phi(T)$.
		Writing $\lambda_j$ in polar coordinates gives
	$\lambda_j=\alpha_j+i\beta_j=r_je^{i\theta_j}=r_j\cos\theta_j + i r_j\sin \theta_j$ with $r_j>0$ and $\theta_j \in (0,2\pi)$. Then,
	defining $R_j=r_j\left(\begin{array}{cc}\cos \theta_j &-\sin  \theta_j\\
	\sin \theta_j&\cos \theta_j\end{array}\right)$,  ${\mathcal R}=\operatorname{blockdiag} (R_j,R_j,\ldots,R_j)\in \mathbb R^{2m_j\times 2m_j}$ and the nilpotent matrix $\mathcal{N}
	\in\mathbb R^{2m_j\times 2m_j}$ having $2\times 2$ blocks of $ \left(\begin{array}{cc}\cos \theta_j &\sin  \theta_j\\
	-\sin \theta_j&\cos \theta_j\end{array}\right)=\left(\begin{array}{cc}\cos \theta_j &-\sin  \theta_j\\
	\sin \theta_j&\cos \theta_j\end{array}\right)^{-1}$ above its diagonal, we have
	$J_j={\mathcal R}(I+\epsilon' \mathcal{N})$.
	We define $\Theta=\left(\begin{array}{cc}0 &-\theta_j\\
	\theta_j&0\end{array}\right)$ and
	$$K_j=\frac{1}{T}\left( (\ln r_j)I+\operatorname{blockdiag}(\Theta,\Theta,\ldots,\Theta)+ \sum_{k=1}^{2m_j-2}\frac{(-1)^{k+1}}{k} (\epsilon') ^k\mathcal{N}^k\right)\in \mathbb R^{2m_j\times 2m_j}.$$
Since $I$, $\operatorname{blockdiag}(\Theta,\Theta,\ldots,\Theta)$ and $\mathcal{N}$ commute, we have, using $\mathcal{N}^k=0$ for $k\ge 2m_j-1$ and \eqref{log}
	\begin{eqnarray*}
		\lefteqn{
			\exp(K_j T)}\\&=&r_j \operatorname{blockdiag}\left( \left(\begin{array}{cc}\cos \theta_j &-\sin  \theta_j\\
			\sin \theta_j&\cos \theta_j\end{array}\right),\ldots,\left(\begin{array}{cc}\cos \theta_j &-\sin  \theta_j\\
			\sin \theta_j&\cos \theta_j\end{array}\right)\right) (I+\epsilon'\mathcal{N})\\
		&=&J_j.
	\end{eqnarray*}

	We can now define $P\in C^s(\mathbb R, \mathbb C^{n\times n})$ by $P(t)=\Phi(t)e^{-Bt}$, which satisfies $P(0)=I$ and
	\begin{eqnarray*}
		P(t+T)&=&\Phi(t+T)e^{-BT}e^{-Bt}\\
		&=&\Phi(t)\Phi(T)e^{-BT}e^{-Bt}\text{ by }\eqref{uni}\\
		&=&P(t)\text{ by }\eqref{*}
		\end{eqnarray*}
		for all $t\ge 0$, so in particular $P(T)=P(0)=I$. We can now write
			$$\Phi(t)=P(t)e^{Bt}.$$
			This shows the first statement of the proposition.

\vspace{0.3cm}
	We now evaluate $A^*+A=\operatorname{blockdiag}(K_1^*+K_1,\ldots,K_r^*+K_r)$. Let us consider $K_j$ as in the three cases above. If $m_j=1$, then $K_j$ below does not contain the last sum with $\epsilon'$ and the form of $c_j$ is immediately clear.

	\vspace{0.3cm}

	\noindent
	\underline{\bf Case 1: $\lambda_j\in \mathbb R^+$}

	$$K_j=\frac{1}{T}\left( (\ln \lambda_j)I+\sum_{k=1}^{m_j-1} \frac{(-1)^{k+1}}{k}(\epsilon')^kN^k\right) \in\mathbb R^{m_j\times m_j};$$
	hence, for ${\bf w}\in \mathbb C^{m_j}$
	\begin{eqnarray*}
		\lefteqn{
		\frac{1}{2}{\bf w}^*(K_j^*+K_j){\bf w}}\\
		&=&
		\frac{\ln \lambda_j}{T}\sum_{i=1}^{m_j}|w_i|^2\\
		&&+\epsilon'
		\frac{1}{2T}
		\left(\overline{w_1}w_2+w_1\overline{w_2}+\overline{w_2}w_3+w_2\overline{w_3}+\ldots+
		\overline{w_{m_j-1}}w_{m_j}+w_{m_j-1}\overline{w_{m_j}}\right)\\
		&&
		-\frac{(\epsilon')^2}{2}
		\frac{1}{2T}
		\left(\overline{w_1}w_3+w_1\overline{w_3}+\overline{w_2}w_4+w_2\overline{w_4}+\ldots+
		\overline{w_{m_j-2}}w_{m_j}+w_{m_j-2}\overline{w_{m_j}}\right)\\
		&&+\ldots\\
		&&
		+(-1)^{m_j}\frac{(\epsilon')^{m_j-1}}{m_j-1}
		\frac{1}{2T}
		\left(\overline{w_1}w_{m_j}+w_1 \overline{w_{m_j}}\right)\,.
	\end{eqnarray*}

	Note that the Cauchy--Schwarz inequality implies $ \mathbb {R}\ni\overline{\xi}\eta+\xi\overline{\eta}\le |\xi|^2+|\eta|^2$, which yields that, using $\epsilon'=\frac{1}{2}\min\left(\frac{\epsilon T}{2},1\right)$
	\begin{eqnarray*}
		\frac{1}{2}{\bf w}^*(K_j^*+K_j){\bf w}
		&\le &		\frac{\ln \lambda_j}{T}\sum_{i=1}^{m_j}|w_i|^2\\
		&&+\frac{\epsilon'+( \epsilon')^2 + \ldots+(\epsilon')^{m_j-1} }{T}\sum_{i=1}^{m_j}|w_i|^2\\
		&\le &
		\left(
		\frac{\ln \lambda_j}{T} +\epsilon\left(\frac{1}{2}+\frac{1}{4}+\frac{1}{8}+\ldots \right)\right)\sum_{i=1}^{m_j}|w_i|^2\\
		&\le&
		\left(
		\frac{\ln \lambda_j}{T} +\epsilon\right)\sum_{i=1}^{m_j}|w_i|^2.
	\end{eqnarray*}

	\vspace{0.3cm}
	\newpage

	\noindent
	\underline{\bf Case 2: $\lambda_j\in \mathbb R^-$}

	$$K_j=\frac{1}{T}\left( (i\pi +\ln |\lambda_j|)I+\sum_{k=1}^{m_j-1}  \frac{(-1)^{k+1}}{k}(-\epsilon')^kN^k\right)\in \mathbb C^{m_j\times m_j};$$
	hence, for ${\bf w}\in \mathbb C^{m_j}$
	\begin{eqnarray*}
		\lefteqn{
		\frac{1}{2}{\bf w}^*(K_j^*+K_j){\bf w}}\\
		&=&
		\frac{\ln |\lambda_j|}{T}\sum_{i=1}^{m_j}|w_i|^2\\
		&&+\epsilon'
		\frac{1}{2T}
		\left(\overline{w_1}w_2+w_1\overline{w_2}+\overline{w_2}w_3+w_2\overline{w_3}+\ldots+
		\overline{w_{m_j-1}}w_{m_j}+w_{m_j-1}\overline{w_{m_j}}\right)\\
		&&
		-\frac{(\epsilon')^2}{2}
		\frac{1}{2T}
		\left(\overline{w_1}w_3+w_1\overline{w_3}+\overline{w_2}w_4+w_2\overline{w_4}+\ldots+
		\overline{w_{m_j-2}}w_{m_j}+w_{m_j-2}\overline{w_{m_j}}\right)\\
		&&+\ldots\\
		&&
		+(-1)^{m_j}\frac{(\epsilon')^{m_j-1}}{m_j-1}
		\frac{1}{2T}
		\left(\overline{w_1}w_{m_j}+w_1 \overline{w_{m_j}}\right)\\
		&\le&
		\left(
		\frac{\ln |\lambda_j|}{T} +\epsilon\right)\sum_{i=1}^{m_j}|w_i|^2
	\end{eqnarray*}
	similarly to case 1.

	\vspace{0.3cm}

	\noindent
	\underline{\bf Case 3: $\lambda_j=\alpha_j+i\beta_j$ with $\beta_j\not=0$ }

	Recall that
	$$K_j=\frac{1}{T}\left( (\ln r_j)I+\operatorname{blockdiag}(\Theta,\Theta,\ldots,\Theta)+ \sum_{k=1}^{2m_j-2} \frac{(-1)^{k+1}}{k}(\epsilon' )^k\mathcal{N}^k\right)\in \mathbb R^{2m_j\times 2m_j};$$
	where $\Theta=\left(\begin{array}{cc}0 &-\theta_j\\
	\theta_j&0\end{array}\right)$ and  the nilpotent matrix $\mathcal{N}$ has $2\times 2$ blocks of $ \left(\begin{array}{cc}\cos \theta_j &\sin  \theta_j\\
	-\sin \theta_j&\cos \theta_j\end{array}\right) $ on its super diagonal.
	Note that all entries of $\mathcal{N}^k$, $k\in \mathbb N$ are real and have an absolute value of $\le 1$ as they are of the form
	$\cos(k\theta_j)$ and $\pm \sin(k\theta_j)$ for $k=1,2,\ldots$.
	Hence, for ${\bf w}\in \mathbb C^{2m_j}$
		\begin{eqnarray*}
		\frac{1}{2}{\bf w}^*(K_j^*+K_j){\bf w}
		&=&
		\frac{\ln r_j}{T}\sum_{i=1}^{2m_j}|w_i|^2\\
		&&+\epsilon'
		\frac{1}{2T}
		\bigg(\cos\theta_j (\overline{w_1}w_3+w_1\overline{w_3})+\sin\theta_j (\overline{w_1}w_4+w_1\overline{w_4})\\
		&&
		-\sin\theta_j (\overline{w_2}w_3+w_2\overline{w_3})+\cos\theta_j (\overline{w_2}w_4+w_2\overline{w_4})+\ldots
		\bigg)+\ldots
		\end{eqnarray*}
		\begin{eqnarray*}
		&\le&
		\frac{\ln r_j}{T}\sum_{i=1}^{2m_j}|w_i|^2\\
		&&+2
		\frac{\epsilon'+(\epsilon')^2+\ldots+(\epsilon')^{2m_j-1}}{T}\sum_{i=1}^{2m_j}|w_i|^2
		\\
		&\le &
		\left(
		\frac{\ln r_j}{T} +\epsilon\left(\frac{1}{2}+\frac{1}{4}+\frac{1}{8}+\ldots \right)\right)\sum_{i=1}^{2m_j}|w_i|^2\\
		&\le&
		\left(
		\frac{\ln r_j}{T} +\epsilon\right)\sum_{i=1}^{2m_j}|w_i|^2
	\end{eqnarray*}
	since $\epsilon'= \min \left(\frac{\epsilon T}{2},1\right)$.
	This shows the second statement of the proposition.

	To show that $(P^{-1}(t))^*(S^{-1})^*S^{-1}P^{-1}(t)$ has real entries, note that
	\begin{eqnarray*}
		P^{-1}(t)&=&e^{Bt}\Phi^{-1}(t)\\
		&=&Se^{At}S^{-1}\Phi^{-1}(t)
		\end{eqnarray*}
		so that
	\begin{eqnarray*}
(P^{-1}(t))^*(S^{-1})^*S^{-1}P^{-1}(t)
&=&(\Phi^{-1}(t))^*(S^{-1})^*
(e^{At})^*e^{At}S^{-1}\Phi^{-1}(t).
\end{eqnarray*}
It is thus sufficient to show that $(e^{At})^*e^{At}$ is real-valued, since all other matrices are real-valued.
Note that since
	$A=\operatorname{blockdiag}(K_1,\ldots,K_r)$, we have
	\begin{eqnarray*}
		e^{At}&=&\operatorname{blockdiag}(e^{tK_1},\ldots,e^{tK_r}),\\
		(e^{tA})^*
		e^{tA}&=&\operatorname{blockdiag}((e^{tK_1})^*e^{tK_1},\ldots,(e^{tK_r})^*e^{tK_r})
	\end{eqnarray*} and the
	blocks where $K_j$ have only real entries are trivially real-valued (cases 1 and 3).
	In case 2,   $K_j=\frac{1}{T}((i\pi+ \ln |\lambda_j|)I+N')$, where $N'\in\mathbb R^{m_j\times m_j}$ is a nilpotent, upper triangular matrix. Then, noting that $I$ and $N'$ commute,
	\begin{eqnarray*}
		e^{tK_j}&=&e^{\frac{  t}{T}(i\pi+\ln |\lambda_j|)} \exp\left(\frac{t}{T}N'\right)\\
		(e^{tK_j})^*&=&e^{\frac{  t}{T}(-i\pi+\ln |\lambda_j|)} \exp\left(\frac{t}{T}(N')^T\right)\\
		(e^{tK_j})^*e^{tK_j}&=&e^{\frac{  2t}{T}\ln |\lambda_j|}\exp\left(\frac{t}{T}(N')^T\right)\exp\left(\frac{t}{T}N'\right),
	\end{eqnarray*}
	which has real entries.
\end{proof}

\begin{corollary}\label{coro}
	Consider the ODE $\dot{{\bf x}}={\bf f}({\bf x})$ with $f\in C^\sigma(\mathbb R^n,\mathbb R^n)$, $\sigma\ge 2$ and let $S_t{\bf q}$ be an exponentially stable periodic solution with period $T$ and ${\bf q}\in \mathbb R^n$.
	 Then the first variation equation $\dot{{\bf y}}=D{\bf f}(S_t {\bf q}){\bf y}$ is  of the form
	as in the previous proposition with $s=\sigma-1$; 1 is a single eigenvalue of $\Phi(T)$ with eigenvector ${\bf f}({\bf q})$ and all other eigenvalues of $\Phi(T)$ satisfy $|\lambda|<1$.
	More precisely, if $-\nu<0$ is the maximal real part of all non-trivial Floquet exponents, we have $\frac{\ln |\lambda|}{T}\le -\nu$.
	With the notations of Proposition \ref{prop} we can assume that $\lambda_1=1$ and $S{\bf e}_1={\bf f}({\bf q})$.

	Then we have for all $\epsilon> 0$
	\begin{eqnarray*}
		{\bf f}(S_t{\bf q})&=&P(t)S{\bf e}_1\text{ for all }t\in\mathbb R\\
	\text{ and }\frac{1}{2}{\bf w}^*(A^*+A){\bf w}&\le& \left(-\nu+\epsilon\right)(\|{\bf w}\|^2-|w_1|^2).
	\end{eqnarray*}
	for all ${\bf w}\in \mathbb C^n$, where $\|{\bf w}\|=\sqrt{{\bf w}^* {\bf w}}$.
\end{corollary}
\begin{proof}
	Since ${\bf f}(S_t{\bf q})$ solves \eqref{var}, we have ${\bf f}(S_t{\bf q})=P(t)e^{Bt}{\bf f}({\bf q})$ and, in particular for $t=T$,  ${\bf f}({\bf q})={\bf f}(S_T{\bf q})=e^{BT}{\bf f}({\bf q})$.
	Hence, \begin{eqnarray*}
		{\bf f}(S_t{\bf q})&=&P(t)
Se^{At}S^{-1}{\bf f}({\bf q})\\
&=&P(t)Se^{At}{\bf e}_1\\
&=&P(t)	S{\bf e}_1
\end{eqnarray*}
since $K_1=0$ in the definition of $A$.
Proposition \ref{prop}  shows the result taking $\lambda_1=1$ and $m_1=1$ into account.
\end{proof}

\section{Converse theorem}\label{conv}

We will prove Theorem \ref{main}, showing that a contraction metric exists for an exponentially stable periodic orbit in the whole basin of attraction. Moreover, we can achieve the bound  $-\nu+\epsilon$ for $L_M$  for any fixed $\epsilon>0$, where $-\nu$ denotes the largest real part of all non-trivial Floquet exponents.

Note that we consider contraction in directions ${\bf v}$ perpendicular to ${\bf f}({\bf x})$ with respect to the  metric $M$, i.e. ${\bf v}^TM({\bf x}){\bf f}({\bf x})=0$. One could alternatively consider directions perpendicular to ${\bf f}({\bf x})$ with respect to the Euclidean metric, i.e. ${\bf v}^T{\bf f}({\bf x})=0$, but then the function $L_M$ needs to reflect this, see \cite{other,leonov1}.

In the proof we will first construct $M=M_0$ on the periodic orbit $\Omega$ using Floquet theory. Then, we define a projection $\pi$ of points in a neighborhood $U$ of $\Omega$ onto $\Omega$ such that $({\bf x}-\pi({\bf x}))^TM_0(\pi({\bf x})){\bf f}(\pi({\bf x}))=0$, which will be used to synchronize the time of solutions such that
$\pi(S_\tau {\bf x})=S_{\theta_{\bf x}(\tau)}\pi({\bf x})$. Finally, $M$ will be defined through a scalar-valued function $V$ by
$M({\bf x})=M_1({\bf x})e^{2V({\bf x})}$, where $M_1=M_0$ on the periodic orbit.

\begin{proof} [of Theorem \ref{main}]
Note that we assume ${\bf f}\in C^\sigma(\mathbb R^n,\mathbb R^n)$ to achieve more detailed results concerning the smoothness and assume lower bounds on $\sigma$ as appropriate for each result; we always assume $\sigma\ge 2$.

\vspace{0.3cm}
\newpage
\noindent \underline{\bf I. Definition and properties of $M_0$ on $\Omega$}

\noindent We fix a point ${\bf q}\in\Omega$ and consider the first variation equation
\begin{eqnarray}
\dot{{\bf y}}&=&D{\bf f}(S_t{\bf q}){\bf y}\label{variation}
\end{eqnarray}
which is a $T$-periodic, linear equation for ${\bf y}$, and $D{\bf f}\in C^{\sigma-1}$. By Proposition \ref{prop} and Corollary \ref{coro}  the principal fundamental matrix  solution $\Phi\in C^{\sigma-1}(\mathbb R,\mathbb R^{n\times n})$ of (\ref{variation}) with $\Phi(0)=I$ can be written as
$$\Phi(t)=P(S_t{\bf q})e^{Bt},$$
where $B\in\mathbb C^{n\times n}$;
note that $P\in C^{\sigma-1}(\mathbb R,  \mathbb C^{n\times n})$ can be defined on the periodic orbit as it is $T$-periodic. By the assumptions on $\Omega$, the eigenvalues of $B$ are $0$ with algebraic multiplicity one and the others have a real part  $\le -\nu<0$.

  We define $S$ as in Proposition \ref{prop} and define the $C^{\sigma-1}  $-function
  \begin{eqnarray}
M_0(S_t {\bf q})&=&P^{-1}(S_t{\bf q})^*(S^{-1})^*S^{-1}P^{-1}(S_t{\bf q})\in\mathbb R^{n\times n}.\label{M_0}
\end{eqnarray}
Note that $M_0(S_t{\bf q})$ is real by Proposition \ref{prop}, symmetric, since it is Hermitian and real, and positive definite  by
\begin{eqnarray}
	{\bf v}^TM_0(S_t{\bf q}){\bf v}	&=&\|S^{-1}P^{-1}(S_t{\bf q}){\bf v}\|^2 \mbox{ for all } {\bf v}\in \mathbb {R}^n\label{N=M}
\end{eqnarray}
and since $S^{-1}P^{-1}(S_t{\bf q})$ is non-singular.

  We will now calculate
   $L_{M_0}(S_t{\bf q};{\bf v})$.    First, we have for the orbital derivative
     $$
  M_0'(S_t{\bf q})=(({P}^{-1}(S_t{\bf q}))')^\ast (S^{-1})^\ast S^{-1}P^{-1}(S_t{\bf q}) +P^{-1}(S_t{\bf q})^\ast (S^{-1})^\ast S^{-1}(P^{-1}(S_t{\bf q}))'\,.
$$

Furthermore, by using $(P^{-1}(S_t{\bf q})P(S_t{\bf q} ))'=0$, we obtain $$(P^{-1}(S_t{\bf q}))'=-P^{-1}(S_t{\bf q})(P(S_t{\bf q}))' P^{-1}(S_t{\bf q}).$$ In addition, since $t\mapsto P(S_t{\bf q} )e^{Bt }$ is a solution of \eqref{variation}, we have
  $(P(S_t{\bf q}))'=D{\bf f}(S_t{\bf q})P(S_t{\bf q})- P(S_t{\bf q} )B$. Altogether, we get
  \begin{eqnarray}
 (P^{-1}(S_t{\bf q} ))'&=&-P^{-1}(S_t{\bf q})D{\bf f}(S_t{\bf q}) +B P^{-1}(S_t{\bf q})\,.\label{alto}
  \end{eqnarray}
  Hence,
  \begin{align*}
  M_0'(S_t{\bf q}) &=-D{\bf f}(S_t{\bf q})^TP^{-1}(S_t{\bf q})^\ast  (S^{-1})^\ast S^{-1}P^{-1}(S_t{\bf q})\\
   &\quad \,+P^{-1}(S_t{\bf q})^\ast B^\ast (S^{-1})^\ast S^{-1}P^{-1}(S_t{\bf q})
    \\
    &\quad\, -P^{-1}(S_t{\bf q})^\ast (S^{-1})^\ast S^{-1}P^{-1}(S_t{\bf q})D{\bf f}(S_t{\bf q})   \\
    &\quad\, +P^{-1}(S_t{\bf q})^\ast (S^{-1})^\ast S^{-1}B P^{-1}(S_t{\bf q}).
  \end{align*}

  Thus, we obtain
  \begin{eqnarray*}
    \lefteqn{M_0(S_t{\bf q}) D{\bf f}(S_t{\bf q})+ D{\bf f}(S_t{\bf q})^TM_0(S_t{\bf q})+ M_0'(S_t{\bf q})}
    \\
    &= &    P^{-1}(S_t{\bf q})^\ast B^\ast (S^{-1})^\ast S^{-1}P^{-1}(S_t{\bf q}) +P^{-1}(S_t{\bf q})^\ast (S^{-1})^\ast S^{-1}B P^{-1}(S_t{\bf q})\,.
  \end{eqnarray*}
  Furthermore, we have for ${\bf v}\in  \mathbb {R}^n$

  \begin{eqnarray}
  L_{M_0}(S_t{\bf q};{\bf v})&=&
  \frac{1}{2}  \lefteqn{ {\bf v}^T\left(M_0(S_t{\bf q}) D{\bf f}(S_t{\bf q})+ D{\bf f}(S_t{\bf q})^TM_0(S_t{\bf q})+ M_0'(S_t{\bf q})\right){\bf v}}\nonumber
    \\
     &=& {\bf v}^T P^{-1}(S_t{\bf q})^\ast (S^{-1})^\ast \left( \frac{1}{2}\left(S^\ast B^\ast (S^{-1})^\ast+S^{-1}BS\right)\right) S^{-1}P^{-1}(S_t{\bf q}){\bf v}\nonumber
    \\
    &=&{\bf w}^\ast\left(\frac{1}{2}\left(A^\ast+A\right)\right){\bf w}\,,\label{**}
  \end{eqnarray}
  where ${\bf w}:=S^{-1}P^{-1}(S_t{\bf q}){\bf v}\in  \mathbb {C}^n$ and $A=S^{-1}BS$.

  For ${\bf v}\in\mathbb R^n$ with ${\bf v}^TM_0(S_t{\bf q}){\bf v}=1$ and ${\bf f}(S_t{\bf q})^TM_0(S_t{\bf q}){\bf v}=0$  we have
  \begin{eqnarray*}
  {\bf w}^*{\bf w}&=&{\bf v}^T P^{-1}(S_t{\bf q})^* (S^{-1})^*
 S^{-1}P^{-1}(S_t{\bf q}){\bf v}\\
  &=&{\bf v}^TM_0(S_t{\bf q}){\bf v}\\
  &=&1
  \end{eqnarray*}
  and, using ${\bf e}_1=S^{-1}P^{-1}(S_t{\bf q}){\bf f}(S_t{\bf q})$ from Corollary \ref{coro} \begin{eqnarray*}
  w_1&=&{\bf e}_1^*  {\bf w}\\
  &=&{\bf f}(S_t{\bf q})^T P^{-1}(S_t{\bf q})^* (S^{-1})^*
  S^{-1}P^{-1}(S_t{\bf q}){\bf v}\\
  &=&{\bf f}(S_t{\bf q})^TM_0(S_t{\bf q}){\bf v}\\
  &=&0.
  \end{eqnarray*}

This shows with Corollary \ref{coro} and \eqref{**}
\begin{eqnarray}
	L_{M_0}(S_t{\bf q})&=&
	\max_{{\bf v}^TM_0(S_t{\bf q}){\bf v}=1,{\bf v}^TM_0(S_t{\bf q}){\bf f}(S_t{\bf q})=0}
	L_{M_0}(S_t {\bf q};{\bf v})\nonumber \\
	&\le&
	\max_{{\bf w}\in \mathbb C^n,w_1=0,\|{\bf w}\|=1}
	(-\nu+\epsilon)(\|{\bf w}\|^2-|w_1|^2)\nonumber\\
	& \le&-\nu+\epsilon.\label{LMv}
\end{eqnarray}

\vspace{0.3cm}
\noindent \underline{\bf II. Projection}

\noindent
Fix a point ${\bf q}\in \Omega$ on the periodic orbit.
For ${\bf x}$ near the periodic orbit we define the projection
$\pi({\bf x})=S_\theta {\bf q}$ on the periodic orbit orthogonal
to ${\bf f}(S_\theta {\bf q})$ with respect to the scalar product $\langle {\bf v},{\bf w}\rangle_{M_0(S_\theta {\bf q})}={\bf v}^T M_0(S_\theta {\bf q}){\bf w}$ implicitly by \eqref{proj} below.
The following lemma is based on the implicit function theorem and shows that the projection can be defined in a neighborhood of the periodic orbit, not just locally.

\begin{lemma}\label{imp}
Let $\Omega$ be an exponentially stable periodic orbit of $\dot{{\bf x}}={\bf f}({\bf x})$ where ${\bf f}\in C^\sigma( \mathbb R^n,\mathbb R^n)$ with $\sigma\ge 2$.

Then there is a compact, positively invariant neighborhood $U$ of $\Omega$ with $U\subset A(\Omega)$ and a function $\pi\in C^{\sigma-1}(U,\Omega)$ such that $\pi({\bf x})={\bf x}$ if and only if ${\bf x}\in \Omega$. Moreover, for all ${\bf x}\in U$ we have
\begin{eqnarray}
	( {\bf x}-\pi({\bf x}))^TM_0(\pi({\bf x})){\bf f}(\pi({\bf x}))&=&0.\label{proj}
	\end{eqnarray}
\end{lemma}
\begin{proof} Fix a point ${\bf q}\in \Omega$ and define $M_0$ by \eqref{M_0}.
Define the $C^{\sigma-1}$ function
$$G({\bf x},\theta)=({\bf x}-S_\theta {\bf q})^TM_0(S_\theta {\bf q}){\bf f}(S_\theta {\bf q})$$
 for ${\bf x}\in \mathbb R^n$, $\theta\in \mathbb R$.

Define the following constants:
\begin{eqnarray*}
	\min_{{\bf p}\in \Omega}\|{\bf f}({\bf p})\|&= & c_1>0\\
	\max_{{\bf p}\in \Omega}\|{\bf f}({\bf p})\|&= & c_2\\
	\max_{{\bf p}\in \Omega}\|D{\bf f}({\bf p})\|&=&c_3\\
	\max_{{\bf p}\in \Omega}\|P({\bf p})\|&=&p_1\\
	\max_{{\bf p}\in \Omega}\|P^{-1}({\bf p})\|&=&p_2\\
	\min_{{\bf p}\in \Omega}\|M_0({\bf p})\|&= & m_1>0\\
	\max_{{\bf p}\in \Omega}\|M_0({\bf p})\|&= & m_2\\
	\max_{{\bf p}\in \Omega}\|M_0'({\bf p})\|&= & m_3,
\end{eqnarray*}
with the matrix norm $\|\cdot\|=\|\cdot\|_2$, which is induced by the vector norm $\|\cdot\|=\|\cdot\|_2$ and is sub-multiplicative.
We will first prove the following quantitative version of the local implicit function theorem, using that $\theta$ is one-dimensional.

\begin{lemma}\label{implicit}
There are constants $\delta,\epsilon>0$ such that
for each point ${\bf x}_0=S_{\theta_0}{\bf q} \in \Omega$, there is a function $p_{{\bf x}_0}\in C^{\sigma-1}( B_\delta({\bf x}_0), B_\epsilon(\theta_0))$ such that
  for all $({\bf x},\theta)\in B_\delta({\bf x}_0)\times B_\epsilon(\theta_0)$
$$G({\bf x},\theta)=0\Longleftrightarrow \theta=p_{{\bf x}_0}({\bf x}).$$
If ${\bf x}\in B_{\delta/2}({\bf x}_0)$, then $p_{{\bf x}_0}({\bf x})\in B_{\epsilon/2}(\theta_0).$
\end{lemma}
\begin{proof}
We have
\begin{eqnarray*}
G_\theta({\bf x},\theta)&=&\frac{d}{d\theta}( {\bf x}-S_\theta {\bf q})^TM_0(S_\theta {\bf q}){\bf f}(S_\theta {\bf q})\\
&=&-{\bf f}(S_\theta {\bf q})^TM_0(S_\theta{\bf q}){\bf f}(S_\theta {\bf q})\\
&&+
( {\bf x}-S_\theta {\bf q})^TM_0'(S_\theta {\bf q}){\bf f}(S_\theta {\bf q})\\
&&
+( {\bf x}-S_\theta {\bf q})^TM_0(S_\theta {\bf q})D{\bf f}(S_\theta
{\bf q}){\bf f}(S_\theta {\bf q})\,.
\end{eqnarray*}
With $\min_{\theta \in
[0,T]}{\bf f}(S_\theta {\bf q})^TM_0(S_\theta{\bf q}){\bf f}(S_\theta {\bf q})\ge c_1^2m_1>0$ we have for all
$ \|{\bf x}-S_\theta {\bf q}\|<\delta_2:= \frac{c_1^2m_1}{2c_2 (m_3+m_1c_3)}$
\begin{eqnarray*}
G_\theta({\bf x},\theta)&<& -c_1^2m_1
+\delta_2c_2 (m_3+m_1c_3)\
= \ -\frac{c_1^2m_1}{2}<0.
\end{eqnarray*}

Let $\delta_1:=\frac{\delta_2}{2} $ and $\epsilon_1:=\frac{\delta_2}{2c_2}$.
For
 any ${\bf x}_0=S_{\theta_0}{\bf q}\in \Omega$ we have for all  ${\bf x}\in \mathbb R^n$ with
$\|{\bf x}-{\bf x}_0\|<\delta_1$ and all $\theta\in \mathbb R$ with  $|\theta-\theta_0|<\epsilon_1$
\begin{eqnarray}
G_\theta({\bf x},\theta)&<& -\frac{c_1^2m_1}{2}<0 \label{Ftheta}
\end{eqnarray} since
$$\|{\bf x}-S_\theta{\bf q}\|\le \|{\bf x}-{\bf x}_0\|+\|S_{\theta_0}{\bf q}-S_\theta {\bf q}\|<
\delta_1+|\theta_0-\theta| c_2<\delta_2.$$

Since $G({\bf x}_0,\theta_0)=0$ we have with $\epsilon:=\epsilon_1/2$ by \eqref{Ftheta}  \begin{eqnarray*}
G({\bf x}_0,\theta_0+\epsilon)&<&-\frac{c_1^2m_1}{2}\epsilon,\\
G({\bf x}_0,\theta_0-\epsilon)&>&\frac{c_1^2m_1}{2}\epsilon.
\end{eqnarray*}

Furthermore, we have
\begin{eqnarray*}
\nabla_{\bf x} G({\bf x},\theta)&=&{\bf f}(S_\theta {\bf q})^TM_0(S_\theta {\bf q}),\\
\|\nabla_{\bf x} G({\bf x},\theta)\|&\le&c_2m_1
\end{eqnarray*}
for all ${\bf x}\in \mathbb R^n$ and $\theta\in \mathbb R$.

Define $\delta:=\min\left(\delta_1,\frac{c_1^2m_1}{4c_2m_1}\epsilon\right)$.
For $\|{\bf x}-{\bf x}_0\|< \delta$ we have
\begin{eqnarray*}G({\bf x},\theta_0+\epsilon)&\le&G({\bf x}_0,\theta_0+\epsilon)\\
	&&+\int_0^1
\nabla_{\bf x} G({\bf x}_0+\lambda ({\bf x}-{\bf x}_0),\theta_0+\epsilon)\,d\lambda \cdot ({\bf x}-{\bf x}_0) \\
&<&-\frac{c_1^2m_1}{2}\epsilon+c_2m_1 \delta \\
&\le &-\frac{c_1^2m_1}{4}\epsilon  \ < \ 0\\
\text{ and }G({\bf x},\theta_0-\epsilon)&>&\frac{c_1^2m_1}{4}\epsilon  \ >\ 0.
\end{eqnarray*}
Since $G({\bf x},\theta)$ is strictly decreasing with respect to $\theta$ in $B_\epsilon(\theta_0)$ by \eqref{Ftheta}  , the intermediate value theorem implies that there is a unique $\theta^*\in (\theta_0-\epsilon,\theta_0+\epsilon)$ such that $G({\bf x},\theta^*)=0$, which defines a function $p_{{\bf x}_0}({\bf x})=\theta^*$. The statement for $\epsilon/2$ and $\delta/2$ follows similarly.
The smoothness of $p_{{\bf x}_0}$ follows by the classical implicit function theorem, since $G\in C^{\sigma-1}$.
\end{proof}

Now we want to show the uniqueness of the function in a suitable neighborhood $\widetilde{U}$ of $\Omega$. Denote the minimal period of the periodic orbit by $T$; we can assume that $\epsilon<T$.
Define
$$c:=\min_{{\bf p}\in \Omega}\min_{\theta\in [-T/2,T/2]\setminus (-\epsilon/2,\epsilon/2)}\|S_\theta{\bf p}-{\bf p}\|>0.$$
We can conclude that if $\|S_\theta{\bf p}-{\bf p}\|\le c/2$ with ${\bf p}\in \Omega$ and $|\theta|\le  T/2$, then $|\theta|<\epsilon/2$.

Let $\delta'=\min(\delta/2,c/4)$. Since $\Omega$ is compact and $\Omega\subset \bigcup_{{\bf x}_0\in \Omega}B_{\delta'}({\bf x}_0)$, there is a finite number of ${\bf x}_i=S_{\theta_i}{\bf q}\in \Omega$, $i=1\ldots,N$, with
\begin{eqnarray}
	\Omega&\subset& \bigcup_{i=1}^N B_{\delta'}({\bf x}_i)=:\widetilde{U},\label{cover}
	\end{eqnarray}
such that $\widetilde{U}$ is an open neighborhood of $\Omega$.
We want to show that the $p_{{\bf x}_i}=p_i$ define a unique function $p\colon \widetilde{U}\to  S^1_T$, where $S^1_T$ are the reals modulo $T$ such that $p=p_i$ on $B_{\delta'}({\bf x}_i)$. We need to show that if ${\bf x}\in B_{\delta'}({\bf x}_i)\cap B_{\delta'}({\bf x}_j)$, then $p_i({\bf x})=p_j({\bf x})$.

Let  ${\bf x}\in B_{\delta'}({\bf x}_i)\cap B_{\delta'}({\bf x}_j)$ and, without loss of generality,  $|\theta_j-\theta_i|\le T/2$ since the $\theta_i$ and $\theta_j$ are uniquely defined only modulo $T$. Then
\begin{eqnarray*}
\|{\bf x}_i-S_{\theta_j-\theta_i}{\bf x}_i\|&=&
\|{\bf x}_i-{\bf x}_j\|\\&\le &\|{\bf x}_i-{\bf x}\|+\|{\bf x}-{\bf x}_j\|\\
&<&2\delta'\ \le\  \min(\delta,c/2).
\end{eqnarray*}
Hence, $|\theta_j-\theta_i|<\epsilon/2$.

Since ${\bf x}\in B_{\delta/2}({\bf x}_i)\cap B_{\delta/2}({\bf x}_j)$, we have $p_i({\bf x})\in B_{\epsilon/2}(\theta_i)$ and  $p_j({\bf x})\in B_{\epsilon/2}(\theta_j)$ by Lemma \ref{implicit}.
 Then
\begin{eqnarray*}
|p_i({\bf x})-\theta_j|
&\le&|p_i({\bf x})-\theta_i| +|\theta_i-\theta_j| \
< \ \epsilon
\end{eqnarray*}
and similarly $p_j({\bf x})\in B_\epsilon(\theta_i)$.
Moreover, ${\bf x}\in B_{\delta}({\bf x}_j)\cap B_{\delta}({\bf x}_j)$.   Lemma \ref{implicit} implies that $\theta=p_i({\bf x})$ if and only if $G({\bf x},\theta)=0$ if and only if $\theta=p_j({\bf x})$, which shows $p_i({\bf x})=p_j({\bf x})$.

Since $\Omega$ is stable, we can choose $\Omega\subset U^\circ\subset U\subset \widetilde{U}$ such that $U$ is compact and positively invariant. For ${\bf x}\in U$ define $\pi({\bf x})=S_{p({\bf x})}{\bf q}$. Since $p$ is defined by $p_{{\bf x}_i}$, we have by Lemma \ref{implicit} $0=G({\bf x},p({\bf x}))=({\bf x}-\pi({\bf x}))^TM_0(\pi({\bf x})){\bf f}(\pi({\bf x}))$.

If ${\bf x}=S_\theta{\bf q} \in \Omega$, then there is a ${\bf x}_i=S_{\theta_i}{\bf q}\in \Omega$ by \eqref{cover} such that ${\bf x}\in B_{\delta'}({\bf x}_i)$ and thus, as above, $|\theta-\theta_i|<\epsilon/2$. Hence, by Lemma \ref{implicit}, as ${\bf x}\in B_\delta({\bf x}_i)$ and $\theta\in B_\epsilon(\theta_i)$, $p_i({\bf x})=\theta$ and thus $\pi({\bf x})={\bf x}$, as this satisfies $0=G({\bf x},\theta)$. If ${\bf x}\not\in \Omega$, then, since $\pi({\bf x})\in \Omega$, ${\bf x}\not=\pi({\bf x})$. This shows the lemma.
\end{proof}

\noindent \underline{\bf III. Synchronization}

\noindent In this step we synchronize the time between the solution $S_t{\bf x}$ and the solution on the periodic orbit $S_\theta \pi({\bf x})$ such that \eqref{first} holds. This will enable us to later define a distance between $S_t{\bf x}$ and $\Omega$ in Step IV.
\begin{definition}
		For ${\bf x}\in U$ we can define $\theta_{\bf x}\in C^{\sigma-1}(\mathbb R^+_0,\mathbb R)$ by $\theta_{\bf x}(0)=0$ and
	\begin{eqnarray}
		S_{\theta_{\bf x}(t)}\pi({\bf x})&=&\pi(S_t {\bf x})\label{first}
		\end{eqnarray}
	for all $t\ge 0$.

	We have \begin{eqnarray}
	\dot{\theta}_{\bf x} (t)
		&=&({\bf f}(S_t{\bf x})^TM_0(S_{\theta_{\bf x}(t)} \pi({\bf x})){\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x})))\nonumber\\
		&&\bigg({\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x}))^TM_0(S_{\theta_{\bf x}(t)} \pi({\bf x})){\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x}))\nonumber\\
		&& -
		( S_t{\bf x}-S_{\theta_{\bf x}(t)} \pi({\bf x}))^T
		\big[M_0'(S_{\theta_{\bf x}(t)} \pi({\bf x})) {\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x}))\nonumber\\
	&&\hspace{0.5cm}	+M_0(S_{\theta_{\bf x}(t)} \pi({\bf x}))D{\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x})) {\bf f}(S_{\theta_{\bf x}(t)} \pi({\bf x}))\big]\bigg)^{-1}.\label{eqbruch}
	\end{eqnarray}
	The denominator of \eqref{eqbruch} is strictly positive for all $t\ge 0$ and ${\bf x} \in U$.
\end{definition}
\begin{proof}
Denote $\pi({\bf x})=:{\bf p}\in \Omega$. Observe, that both sides of \eqref{first} equal for $t=0$.
 For any $t\ge 0$, $S_t{\bf x}\in U$ and  $\pi(S_t {\bf x})$ denotes a point on the periodic orbit, so we can write it as
$\pi(S_t {\bf x})=S_{\theta_{\bf x}(t)}{\bf p}$. Note that $\theta_{\bf x}(t)$ is only uniquely defined modulo $T$, however, it is uniquely defined by the requirement that $\theta_{\bf x}$  is a continuous function.

By \eqref{proj}, we have
$$(S_t {\bf x}-S_{\theta_{\bf x}(t)}{\bf p})^TM_0(S_{\theta_{\bf x}(t)} {\bf p}){\bf f}(S_{\theta_{\bf x}(t)}{\bf p})=0.$$
Hence, $\theta_{\bf x}(t)$ is implicitly defined by
 \begin{eqnarray}
 	Q(t,\theta)&=&(S_t{\bf x}-S_\theta {\bf p})^TM_0(S_\theta {\bf p}){\bf f}(S_\theta {\bf p})=0.
 	\label{theta}
 \end{eqnarray}
 Note that $\theta_{\bf x}\in C^{\sigma-1}(\mathbb R^+_0,\mathbb R)$ by the Implicit
 Function Theorem which implies
  \begin{eqnarray*}
 	\frac{ d\theta_{\bf x}}{dt}
 	&=&-\frac{\partial_t Q(t,\theta)}{\partial_\theta Q(t,\theta)}\bigg|_{\theta=\theta_{\bf x}(t)}
	\nonumber\\
 	&=&({\bf f}(S_t{\bf x})^TM_0(S_{\theta_{\bf x}(t)} {\bf p}){\bf f}(S_{\theta_{\bf x}(t)} {\bf p}))\nonumber\\
 	&&\bigg({\bf f}(S_{\theta_{\bf x}(t)} {\bf p})^TM_0(S_{\theta_{\bf x}(t)} {\bf p}){\bf f}(S_{\theta_{\bf x}(t)} {\bf p})\nonumber\\
 	&&\hspace{0.3cm}-
 	( S_t{\bf x}-S_{\theta_{\bf x}(t)} {\bf p})^TM_0'(S_{\theta_{\bf x}(t)} {\bf p}) {\bf f}(S_{\theta_{\bf x}(t)} {\bf p})\nonumber\\
 	&&\hspace{0.3cm}-( S_t{\bf x}-S_{\theta_{\bf x}(t)} {\bf p})^TM_0(S_{\theta_{\bf x}(t)} {\bf p})D{\bf f}(S_{\theta_{\bf x}(t)} {\bf p}) {\bf f}(S_{\theta_{\bf x}(t)} {\bf p})\bigg)^{-1}.
 	 \end{eqnarray*}

	With the notations of the proof of Lemma \ref{imp}, for $S_t{\bf x}\in U$ there is a point ${\bf x}_i=S_{\theta_i} {\bf q} \in \Omega$ such that $S_t{\bf x}\in B_{\delta'}({\bf x}_i)$. We have $S_t{\bf x}\in B_\delta({\bf x}_i)$ and, modulo $T$, we have $p_i(S_t{\bf x})=\theta_{\bf x}(t)\in B_\epsilon(\theta_i)$. Hence, the denominator is $>\frac{c_1^2m_1}{2}$ by \eqref{Ftheta}.
 \end{proof}

\begin{lemma}\label{change}
For ${\bf x}\in U$ we have
$$S_{\theta_{S_\tau {\bf x}}(t)}\pi(S_\tau {\bf x})=S_{\theta_{\bf x}(t+\tau)}\pi({\bf x})$$
for all $t,\tau\ge 0$.
\end{lemma}
\begin{proof}
We apply \eqref{first} to the point $S_\tau {\bf x}$ and the time $t$, obtaining
$$S_{\theta_{S_\tau {\bf x}}(t)}\pi(S_\tau {\bf x})=\pi(S_t S_\tau {\bf x}).$$
Now we apply \eqref{first} to the point $ {\bf x}$ and the time $t+\tau$, obtaining
$$S_{\theta_{\bf x}(t+\tau)}\pi({\bf x})=\pi(S_{t+\tau} {\bf x}).$$
As both right-hand sides are equal by the semi-flow property, this proves the
statement.
\end{proof}

\vspace{0.3cm}
\noindent \underline{\bf IV. Distance to the periodic orbit}

\noindent
In the following lemma we define a distance of points in $U$ to the periodic orbit, and we show that it decreases exponentially.

\begin{lemma}\label{def_d}
Let $\epsilon<\min(1,\nu/2)$ and $\sigma\ge 2$. Then there is a positively invariant, compact neighborhood $U$ of the periodic orbit $\Omega$ such that
	the function $d\in C^{\sigma-1}(U,\mathbb R^+_0)$, defined by
$$d({\bf x})=({\bf x}-\pi({\bf x}))^TM_0(\pi({\bf x}))({\bf x}-\pi({\bf x}))$$
 satisfies $d({\bf x})=0$ if and only if ${\bf x}\in \Omega$. Moreover, $d'({\bf x})<0$ for all ${\bf x}\in U\setminus \Omega$
and
\begin{eqnarray*}
d(S_t{\bf x})&\le& e^{2(-\nu+2\epsilon) t}d({\bf x})\text{    for all ${\bf x}\in U$ and all $t\ge 0$,}\\
1-\epsilon&\le& \dot{\theta}_{\bf x}(t)\ \le  \ 1+\epsilon\text{  for all $t\ge 0$}.
\end{eqnarray*}
\end{lemma}

\begin{proof}
Note that $d$ is $C^{\sigma-1}$ as all of its terms are. As $M_0({\bf x})$ is positive definite, $d({\bf x})=0$ if and only if ${\bf x}=\pi({\bf x})$, i.e. ${\bf x}\in \Omega$ by Lemma \ref{imp}.
Define
\begin{eqnarray}
c^*&:=&\frac{ \epsilon}{2 p_1\, p_2 \|S^{-1}\| \,  \|S\| }>0,
\label{defc}\\
c_4&=&2c_2\frac{2c_3m_2+m_3+c^*m_2}{c_1^2m_1},\label{c4def}\\
v^*&:=&\frac{\epsilon}{2 p_1\, p_2 \|S^{-1}\| \,  \|S\| \, c_4(c_3+\|B\|)}\,.\label{defv}
\end{eqnarray}
where the constants where defined in Step II, proof of Lemma \ref{imp}.

For ${\bf y}\in U$ we use the Taylor expansion around $ \pi({\bf y})\in \Omega$. Hence, there is a function ${\boldsymbol \psi}({\bf y})$ satisfying
\begin{eqnarray}
{\bf f}({\bf y})&=&{\bf f}(\pi({\bf y}))+D{\bf f}(\pi({\bf y}))({\bf y}-\pi({\bf y}))+{\boldsymbol \psi}({\bf y})\label{Taylor}
\end{eqnarray}
with $\|{\boldsymbol \psi}({\bf y})\|\le c^* \|{\bf y}-\pi({\bf y})\|$ for all ${\bf y} \in U$, noting that $\Omega$ is compact,
 where we choose $U$ still to be a positively invariant, compact neighborhood of $\Omega$, possibly smaller than before and such that also have
\begin{eqnarray}
\|{\bf y}-\pi({\bf y})\|&\le&\delta'=\min\left(v^*, \frac{c_1^2m_1}{2c_2[m_3+ m_2c_3] },\frac{\epsilon}{c_4},1\right)\text{ for all ${\bf y}\in U$.}\label{vest}
\end{eqnarray}

Recall that, due to the definition of $M_0$ and \eqref{first} we have
\begin{eqnarray*}
d({\bf x})&=&({\bf x}-\pi({\bf x}))^T(P^{-1}(\pi({\bf x})))^*(S^{-1})^*S^{-1}P^{-1}(\pi({\bf x}))({\bf x}-\pi({\bf x}))\\
d(S_t{\bf x})&=&(S_t{\bf x}-S_{\theta_{\bf x}(t)}\pi({\bf x}))^T(P^{-1}(S_{\theta_{\bf x}(t)}\pi({\bf x})))^*(S^{-1})^*\\
&&S^{-1}P^{-1}(S_{\theta_{\bf x}(t)}\pi({\bf x}))(S_t{\bf x}-S_{\theta_{\bf x}(t)}\pi({\bf x})).
\end{eqnarray*}

Now let us calculate the orbital derivative, denoting
$\theta(t):=\theta_{\bf x}(t)$.
\begin{eqnarray*}
d'(S_t{\bf x})&=&\bigg[\frac{d}{dt}\left(P^{-1}(S_{\theta(t)}\pi({\bf x}))\right)(S_t{\bf x}-S_{\theta(t)}\pi({\bf x}))\\
&&+P^{-1}(S_{\theta(t)}\pi({\bf x}))[{\bf f}(S_t{\bf x})-\dot{\theta}(t){\bf f}(S_{\theta(t)}\pi({\bf x}))]\bigg]^*\\
&&(S^{-1})^*S^{-1}P^{-1}(S_{\theta(t)}\pi({\bf x}))(S_t{\bf x}-S_{\theta(t)}\pi({\bf x}))\\
&&+
(S_t{\bf x}-S_{\theta(t)}\pi({\bf x}))^T(P^{-1}(S_{\theta(t)}\pi({\bf x})))^*(S^{-1})^*S^{-1}\\
&&
\bigg[\frac{d}{dt}\left(P^{-1}(S_{\theta(t)}\pi({\bf x}))\right)(S_t{\bf x}-S_{\theta(t)}\pi({\bf x}))\\
&&+P^{-1}(S_{\theta(t)}\pi({\bf x}))[{\bf f}(S_t{\bf x})-\dot{\theta}(t){\bf f}(S_{\theta(t)}\pi({\bf x}))]\bigg].
\end{eqnarray*}
We denote ${\bf p}:=\pi({\bf x})$ and  ${\bf v}(t):=S_t{\bf x}-S_{\theta(t)}\pi({\bf x})=S_t{\bf x}-\pi(S_t{\bf x})$ by  \eqref{first}.
Hence, using \eqref{vest} for ${\bf y}=S_t{\bf x}\in U$ since ${\bf x}\in U$, which is positively invariant, we have
\begin{eqnarray}
\|{\bf v}(t)\|&\le&\delta'=\min\left(v^*, \frac{c_1^2m_1}{2c_2[m_3+ m_2c_3] },\frac{\epsilon}{c_4},1\right)\label{vest2}
\end{eqnarray}
for all $t\ge 0$.
 We have
 $\frac{d}{dt}\left(P^{-1}(S_{\theta(t)}\pi({\bf x}))\right)=\dot{\theta}(t)(-P^{-1}(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p})+B P^{-1}(S_{\theta(t)}{\bf p}))$ by \eqref{alto}.
 Thus,
\begin{eqnarray}
d'(S_t{\bf x})&=&
\bigg[\dot{\theta}(t)(-P^{-1}(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p})+B P^{-1}(S_{\theta(t)}{\bf p})){\bf v}(t)\nonumber\\
&&+P^{-1}(S_{\theta(t)}{\bf p})[{\bf f}(S_t{\bf x})-\dot{\theta}(t){\bf f}(S_{\theta(t)}{\bf p})]\bigg]^*\nonumber\\
&&(S^{-1})^*S^{-1}P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)\nonumber\\
&&+{\bf v}(t)^T(P^{-1}(S_{\theta(t)}{\bf p}))^*(S^{-1})^*S^{-1}\nonumber\\
&&
\bigg[\dot{\theta}(t)(-P^{-1}(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p})+B\nonumber P^{-1}(S_{\theta(t)}{\bf p})){\bf v}(t)\nonumber\\
&&+P^{-1}(S_{\theta(t)}{\bf p})[{\bf f}(S_t{\bf x})-\dot{\theta}(t){\bf f}(S_{\theta(t)}{\bf p})]\bigg].\label{eq33}
\end{eqnarray}

Using the Taylor expansion \eqref{Taylor} for ${\bf y}=S_t{\bf x}$, we obtain  with $\pi(S_t{\bf x})=S_{\theta(t)}{\bf p}$,
\begin{eqnarray}
{\bf f}(S_t{\bf x})&=&{\bf f}(S_{\theta(t)}{\bf p})+D{\bf f}(S_{\theta(t)}{\bf p}) {\bf v}(t)+{\boldsymbol \psi}(S_t{\bf x})
\label{Taylor2}
\end{eqnarray}
and thus with \eqref{eqbruch}
 \begin{eqnarray*}
 	\lefteqn{
 \dot{\theta}(t)-1}\nonumber\\
&=&\bigg( {\bf f}(S_t{\bf x})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})-{\bf f}(S_{\theta(t)}{\bf p})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\\
&&\hspace{0.5cm}+{\bf v}(t)^TM_0'(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})+{\bf v}(t)^TM_0(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\bigg)\\
&&\bigg({\bf f}(S_{\theta(t)}{\bf p})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})-{\bf v}(t)^TM_0'(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\\
&&\hspace{0.5cm}-{\bf v}(t)^TM_0(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\bigg)^{-1}\\
 &=&\bigg( {\bf v}(t)^TD{\bf f}(S_{\theta(t)}{\bf p})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})+{\boldsymbol \psi}(S_t{\bf x})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\\
  &&+\hspace{0.5cm}  {\bf v}(t)^TM_0'(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})+{\bf v}(t)^TM_0(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\bigg)\\
&&\bigg({\bf f}(S_{\theta(t)}{\bf p})^TM_0(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})-{\bf v}(t)^TM_0'(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\\
&&\hspace{0.5cm}-{\bf v}(t)^TM_0(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})\bigg)^{-1}
 \end{eqnarray*}
 which shows, using \eqref{vest2} and \eqref{c4def},
 \begin{eqnarray}
  | \dot{\theta}(t)-1|&\le&\frac{\|{\bf v}(t)\| c_2[2c_3 m_2 +m_3]+
  \|{\boldsymbol \psi}(S_t{\bf x})\| m_2c_2  }{c_1^2m_1-\|{\bf v}(t)\|c_2[m_3+ m_2c_3] }\nonumber\\
   &\le&2c_2\frac{ 2c_3 m_2 +m_3+
  c^* m_2}{c_1^2m_1 }\|{\bf v}(t)\|=c_4\|{\bf v}(t)\|\le \epsilon\,.\label{est}
  \end{eqnarray}
 In particular, we have
  $1-\epsilon\le  \dot{\theta}(t)\le 1+\epsilon$, which shows the existence of $\theta(t)$ for all $t\ge 0$, $\dot{\theta}(t)>0$ for all $t\ge 0$, that $\theta(t)$ is a bijective function from $[0,\infty)$ to $[0,\infty)$ and $\lim_{t\to\infty}\theta(t)=\infty$.

  Hence, we have from \eqref{eq33} and \eqref{Taylor2}
    \begin{eqnarray*}
d'(S_t{\bf x})&=&\big[(1-\dot{\theta}(t))P^{-1}(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf v}(t)
+B P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)\\
&&\hspace{0.5cm}-(1-\dot{\theta}(t))B P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)+(1-\dot{\theta}(t))
P^{-1}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})
\\
&&\hspace{0.5cm}+P^{-1}(S_{\theta(t)}{\bf p}){\boldsymbol \psi}(S_t{\bf x})\big]^*
(S^{-1})^*S^{-1}P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)\\
&&+{\bf v}(t)^T(P^{-1}(S_{\theta(t)}{\bf p}))^*(S^{-1})^*S^{-1}\big[(1-\dot{\theta}(t))P^{-1}(S_{\theta(t)}{\bf p})D{\bf f}(S_{\theta(t)}{\bf p}){\bf v}(t)\\
&&\hspace{0.5cm}+B P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)-(1-\dot{\theta}(t))B P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)\\
&&\hspace{0.5cm}
+(1-\dot{\theta}(t))
P^{-1}(S_{\theta(t)}{\bf p}){\bf f}(S_{\theta(t)}{\bf p})+P^{-1}(S_{\theta(t)}{\bf p}){\boldsymbol \psi}(S_t{\bf x})\big]\\
&\le&2
\|S^{-1}P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)\|\, \|S^{-1}\| \,  \|P^{-1}(S_{\theta(t)}{\bf p})\|\,\\
&&\hspace{0.5cm}
\left[|1-\dot{\theta}(t)|(\|D{\bf f}(S_{\theta(t)}{\bf p})\| +\|B\|)\|{\bf v}(t)\|+\|{\boldsymbol \psi}(S_t{\bf x})\|\right]\\
&&+{\bf v}(t)^*(P^{-1}(S_{\theta(t)}{\bf p}))^*\left[(S^{-1})^*S^{-1}B +B^*(S^{-1})^*S^{-1}
\right]P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)
\end{eqnarray*}
using $$0={\bf f}(S_{\theta(t)}{\bf p})^*M_0(S_{\theta(t)}{\bf p}){\bf v}(t)={\bf f}(S_{\theta(t)}{\bf p})^*(P^{-1}(S_{\theta(t)}{\bf p}))^*(S^{-1})^*S^{-1}P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)$$ by \eqref{theta}.

Setting ${\bf w}(t)=S^{-1}P(S_{\theta(t)}{\bf p})^{-1}{\bf v}(t)$, we obtain, using \eqref{est}
and \eqref{vest2}
\begin{eqnarray*}
d'(S_t{\bf x})
&\le&2\, p_2
\|{\bf w}(t)\|\, \|S^{-1}\|  \|{\bf v}(t)\|
\left[c_4  (c_3+\|B\|) \|{\bf v}(t)\|+c^*\right]\\
&&+{\bf w}(t)^* \left[S^{-1}B S+S^*B^*(S^{-1})^*\right]{\bf w}(t)\\
&\le&2  p_1 p_2 \|S\| \,
 \|S^{-1}\|  \left[c_4 (c_3+\|B\|)\|{\bf v}(t)\|+c^*\right]\,\|{\bf w}(t)\|^2\\
 &&+{\bf w}(t)^* \left[A+A^*\right]{\bf w}(t)\\
&\le&2\epsilon\,\|{\bf w}(t)\|^2+{\bf w}(t)^* \left[A+A^*\right]{\bf w}(t)
\end{eqnarray*}
by \eqref{defv} and \eqref{defc}.
 Noting that
$$w_1(t)={\bf e}_1^*{\bf w}(t)={\bf f}(S_{\theta(t)}{\bf p})^*(P^{-1}(S_{\theta(t)}{\bf p}))^*(S^{-1})^*S^{-1}P^{-1}(S_{\theta(t)}{\bf p}){\bf v}(t)=0$$
we have   with Corollary \ref{coro}
$${\bf w}(t)^* \left[A+A^*\right]{\bf w}(t) \le   2 (-\nu+\epsilon)
\|{\bf w}(t)\|^2.$$

Altogether, we have
\begin{eqnarray*}
d'(S_t{\bf x})
&\le&
\left[  2\epsilon-2\nu+2\epsilon\right]\|{\bf w}(t)\|^2\\
&=&2(-\nu+2\epsilon) d(S_t{\bf x}),
\end{eqnarray*}which shows  $d(S_t{\bf x})\le e^{2(-\nu+2\epsilon)  t}d({\bf x})$ for all ${\bf x}\in U$ and $t\ge 0$.
\end{proof}

Let us summarize the results obtained so far in the following corollary.

\begin{corollary}\label{help_other}
Let $\Omega$ be an exponentially stable periodic orbit of $\dot{{\bf x}}={\bf f}({\bf x})$ with ${\bf f}\in C^\sigma(\mathbb R^n,\mathbb R^n)$ and $\sigma\ge 2$, such that
$-\nu<0$ is the maximal real part of all non-trivial Floquet exponents.

	For   $\epsilon_0\in (0,\min(\nu,1))>0$  there is a compact, positively invariant neighborhood $U$ of $\Omega$ with $\Omega\subset U^\circ$ and $U\subset A(\Omega)$, and a map $\pi\in C^{\sigma-1}( U, \Omega)$ with $\pi({\bf x})={\bf x}$ if and only if ${\bf x}\in\Omega$.

	Furthermore, for a fixed ${\bf x}\in U$, there is a bijective $C^{\sigma-1}$ map
	$\theta_{\bf x}\colon [0,\infty)\to[0,\infty)$ with inverse
		$t_{\bf x}=\theta_{\bf x}^{-1}\in C^{\sigma-1}( [0,\infty),[0,\infty))$
	such that $\theta_{\bf x}(0)=0$ and
	$$\pi(S_t{\bf x})=S_{\theta_{\bf x}(t)}\pi({\bf x})$$
	for all $t\in [0,\infty)$.
	We have $\dot{\theta}_{\bf x}(t)\in \left[1-\epsilon_0,1+\epsilon_0\right]$ for all $t\ge 0$ and
	 $\dot{t}_{\bf x}(\theta)\in \left[1-\epsilon_0,1+\epsilon_0\right]$ for all $\theta\ge 0$.

	 Finally, there is a constant $C>0$ such that
		\begin{eqnarray}|\dot{t}_{\bf x}(\theta)-1|&\le& Ce^{(-\nu+\epsilon_0) \theta}\label{res1}\\
	\|S_{t_{\bf x}(\theta)}{\bf x}-S_\theta \pi({\bf x})\|&\le& Ce^{(-\nu+\epsilon_0)  \theta}\|{\bf x}-\pi({\bf x})\|\label{res2}
	\end{eqnarray}
	for all $\theta\ge 0$ and all ${\bf x}\in U$.
	\end{corollary}
\begin{proof}
Setting $\epsilon:=\frac{\epsilon_0}{2(1+\nu)}\le\min\left( \frac{\epsilon_0}{2},\frac{1}{2}\right)\le  \min\left( \frac{\nu}{2},1\right)$, all results  follow directly from Lemma \ref{def_d}   by using the inverse $t(\theta)$ of $\theta(t)$. Indeed,  we have\begin{eqnarray*}
|\dot{t}_{\bf x}(\theta)-1|&=&\left|\frac{1-\dot{\theta}_{\bf x}(t(\theta))}{\dot{\theta}_{\bf x}(t(\theta))}\right|\\
&\le&\frac{\epsilon}{1-\epsilon}\\
&\le&2\epsilon\le\epsilon_0\,.
\end{eqnarray*}
Furthermore, we have by  \eqref{est} and noting that
$ m_1 \|S_{t_{\bf x}(\theta)}{\bf x}-S_\theta \pi({\bf x})\|^2\le d(S_{t_{\bf x}(\theta)}{\bf x}) \le m_2 \|S_{t_{\bf x}(\theta)}{\bf x}-S_\theta \pi({\bf x})\|^2$
\begin{eqnarray*}
|\dot{t}_{\bf x}(\theta)-1|&\le&
\left|\frac{1-\dot{\theta}_{\bf x}(t(\theta))}{1/2}\right|\\
&\le&2c_4\|{\bf v}(t(\theta))\|\\
&\le&\frac{2c_4}{\sqrt{m_1}}\sqrt{d(S_{t(\theta)}{\bf x})}\\
&\le&Ce^{(-\nu+2\epsilon)t(\theta)}\sqrt{d({\bf x})}\\
&\le&Ce^{(-\nu+2\epsilon)(1-2\epsilon)\theta}\\
&\le&Ce^{(-\nu+2\epsilon(1+\nu)-4\epsilon^2)\theta}\\
&\le&Ce^{(-\nu+\epsilon_0)\theta},
\end{eqnarray*}
using $t(\theta)=\int_0^\theta \dot{t}(\tau)\,d\tau\ge \theta(1-2\epsilon)$ and that $d({\bf x})$ is bounded in $U$.
Similarly, we can prove \eqref{res2} from Lemma \ref{def_d}.
\end{proof}

\vspace{0.3cm}
\noindent \underline{\bf V. Definition of $M_1$ and $M$ in $A(\Omega)$}

\noindent
For all ${\bf x}\in U$ we have defined the distance
$$d({\bf x})=({\bf x}-\pi({\bf x}))^TM_0(\pi({\bf x}))({\bf x}-\pi({\bf x}))$$
in Lemma \ref{def_d}
which is $C^{\sigma-1}$. Let $\iota>0$ be so small that
the set $\Omega_{2\iota}:=\{{\bf x}\in U\colon d({\bf x})\le 2\iota\}$ satisfies $\Omega_{2\iota}\subset U^\circ$.  Define
the $C^\infty$ functions $h_1\colon \Omega_\iota \to [0,1]$, $h_2\colon \Omega_{2\iota} \to [0,1]$ such that $h_1({\bf x})=1$ for all $d({\bf x})\le \frac{\iota}{3}$ and $h_1({\bf x})=0$ for all $d({\bf x})\ge \frac{2}{3} \iota$, and $h_2({\bf x})=1$ for all $d({\bf x})\le \frac{4}{3}\iota$ and $h_2({\bf x})=0$ for all $d({\bf x})\ge \frac{5}{3} \iota$.  Set
$$M_1({\bf x}):=\left\{\begin{array}{ll}
I&\mbox{ if }{\bf x}\not \in  \Omega_{2\iota},\\
(1-h_2({\bf x}))I+h_2({\bf x})M_0(\pi({\bf x}))&\mbox{ if }{\bf x} \in  \Omega_{2\iota}.\end{array}\right.$$
It is clear that $M_1({\bf x})$ is positive definite for all ${\bf x}\in \mathbb R^n$, $M_1$ is $C^{\sigma-1}$  and $M_1(\pi({\bf x}))=M_0(\pi({\bf x}))$ for all ${\bf x}\in \Omega_{\frac{4}{3}\iota}$.

We will define the Riemannian metric $M$ through $M_1$ and a scalar-valued function $V\colon A(\Omega)\to \mathbb R$, which will be defined later.
Let us denote $\mu:=\nu-\epsilon>0$. The function $V$ will be continuous and continuously orbitally differentiable and satisfy
\begin{eqnarray}
V'({\bf x})&=&-L_{M_1}({\bf x})+r({\bf x}), \text{ where }\label{eqV}\\
r({\bf x})&=&\left\{\begin{array}{ll}
-\mu&\mbox{ if }{\bf x}\not \in  \Omega_\iota,\\
-\mu(1-h_1({\bf x}))+h_1({\bf x})L_{M_1}(\pi({\bf x}))&\mbox{ if }{\bf x} \in  \Omega_\iota.\end{array}\right.
\end{eqnarray}
Note that $$r({\bf x})\le -\mu$$
for all ${\bf x}\in \mathbb R^n$. Indeed, for ${\bf x}\in\Omega_\iota$ we have   $L_{M_1}(\pi({\bf x}))=L_{M_0}(\pi({\bf x}))\le -\mu$ as $\pi({\bf x})\in \Omega$, see \eqref{LMv}, and thus
\begin{eqnarray*}
r({\bf x})&=& -\mu+\underbrace{h_1({\bf x})}_{\ge 0}\underbrace{(\mu+L_{M_1}(\pi({\bf x}))}_{\le 0}
\ \le\  -\mu\,.
\end{eqnarray*}

Then we define
$$M({\bf x})=e^{2V({\bf x})}M_1({\bf x}).$$
We obtain by Lemma \ref{lem}
$$L_M({\bf x})=L_{M_1}({\bf x})+V'({\bf x})
=L_{M_1}({\bf x})-L_{M_1}({\bf x})+r({\bf x})\le-\mu.$$

This shows the theorem. In the last steps we will define the function $V$ and prove the properties stated above.

\vspace{0.3cm}
\noindent \underline{\bf VI. Definition of $V_{loc}$}

\noindent
We define $V_{loc}({\bf x})$ for $ {\bf x}\in \Omega_\iota$. Note that $\Omega_\iota$ is positively invariant by Lemma \ref{def_d}, so
$S_t {\bf x}\in \Omega_\iota$ for all $t\ge 0$. We define
\begin{eqnarray}
V_{loc}({\bf x})&=&\int_0^\infty [L_{M_1}(S_t {\bf x})-L_{M_1}(S_{\theta_{\bf x}(t)}
\pi({\bf x}))]\,dt.\label{Vdefi}
\end{eqnarray}
We have $V_{loc}({\bf x})=0$ for all ${\bf x}\in\Omega$.
We will show that the $V_{loc}$ is well-defined, continuous and orbitally continuously differentiable for all ${\bf x}\in \Omega_\iota$ and that
 (\ref{eqV}) holds for all ${\bf x}\in \overline{\Omega_{\iota/3}}$.

For ${\bf x}\in U$, define
\begin{eqnarray*}
g_T(\tau,{\bf x})&=&\int_\tau^{T+\tau} [L_{M_1}(S_t {\bf x})-L_{M_1}(S_{\theta_{\bf x}(t)}
\pi({\bf x}))]\,dt.
\end{eqnarray*}
  By Lemma \ref{def_d} there is a constant $C>0$ such that, defining ${\bf p}:=\pi({\bf x})\in \Omega$,
 \begin{eqnarray}
 \label{expsta}
 \|S_t{\bf x}-S_{\theta_{\bf x}(t)}{\bf p}\|&\le& Ce^{-\mu_0 t}
 \end{eqnarray}
 for all $t\ge 0$ and all ${\bf x}\in U$ with $\mu_0:=\nu-2\epsilon>0$; note that $S_{\theta_{\bf x}(t)}{\bf p}=\pi(S_t{\bf x})$ by \eqref{first}.

 Now, we use  Lemma \ref{Lipschitz} and $\sigma\ge 3$, showing that $L_{M_1}$ is Lipschitz-continuous on the compact set $U$; note that $\sigma-1\ge 2$.

 Hence,
  \begin{eqnarray*}
 \left|L_{M_1}(S_t {\bf x})-L_{M_1}(S_{\theta_{\bf x}(t)}
\pi({\bf x}))\right|
&\le&L  C_1\left\|S_t {\bf x}-S_{\theta_{{\bf x}}(t)}{\bf p}\right\|\\
&\le&L C_2 e^{-\mu_0 t}
\end{eqnarray*}
by \eqref{expsta},
which is integrable over $[0,\infty)$. Hence, by Lebesgue's dominated convergence theorem, the function $g_T(\tau,{\bf x})$ converges point-wise for $T\to \infty$ for all $\tau\ge 0$  and ${\bf x}\in U$.

Choose $\theta_0>0$ so small that $S_{-\theta_0}\Omega_{\iota}\subset U$.
We have that
\begin{eqnarray*}
\lefteqn{\frac{\partial }{\partial \tau}g_T(\tau,{\bf x})}\\&=&
[L_{M_1}(S_{T+\tau} {\bf x})-L_{M_1}(S_{\theta_{\bf x}(T+\tau)}
\pi({\bf x}))]
- \left(L_{M_1}(S_{\tau}{\bf x})-L_{M_1}(S_{\theta_{{\bf x}}(\tau)}{\bf p})\right)\\
&=&
[L_{M_1}(S_T(S_\tau {\bf x}))-L_{M_1}(S_{\theta_{S_\tau{\bf x}}(T)}
\pi(S_\tau{\bf x}))]
- \left(L_{M_1}(S_{\tau}{\bf x})-L_{M_1}(S_{\theta_{{\bf x}}(\tau)}{\bf p})\right)
\end{eqnarray*}
by Lemma \ref{change}. For ${\bf x}\in\Omega_\iota$, the right-hand side
converges uniformly in $\tau \in  (-\theta_0,\theta_0)$
as $T\to \infty$ to
$- \left(L_{M_1}(S_{\tau}{\bf x})-L_{M_1}(S_{\theta_{{\bf x}}(\tau)}{\bf p})\right)$ by the same estimate as above. Hence, we can exchange $\frac{d}{d\tau}$ and $\lim_{T\to\infty}$.
Altogether, we thus have for all ${\bf x}\in\Omega_\iota$, using Lemma \ref{change}
\begin{eqnarray*}
V_{loc}'({\bf x})&=&\frac{d}{d\tau}V_{loc}(S_\tau {\bf x})\bigg|_{\tau=0}\\
&=&\frac{d}{d\tau}\int_0^\infty [L_{M_1}(S_{t+\tau} {\bf x})-L_{M_1}(S_{\theta_{S_\tau{\bf x}}
(t)}
\pi(S_\tau{\bf x}))]\,dt\bigg|_{\tau=0}\\
&=&\frac{d}{d\tau}\lim_{T\to\infty}\int_0^T  [L_{M_1}(S_{t+\tau} {\bf x})-L_{M_1}(S_{\theta_{\bf x}(t+\tau)}
\pi({\bf x}))]\,dt\bigg|_{\tau=0}\\
&=&\frac{d}{d\tau}\lim_{T\to\infty}\int_\tau^{T+\tau} [L_{M_1}(S_{t} {\bf x})-L_{M_1}(S_{\theta_{\bf x} (t)}
\pi({\bf x}))]\,dt\bigg|_{\tau=0}\\
&=&\frac{d}{d\tau}\lim_{T\to \infty}g_T(\tau,{\bf x})\bigg|_{\tau=0}\\
&=&\lim_{T\to \infty}\frac{d}{d\tau}g_T(\tau,{\bf x})\bigg|_{\tau=0}\\
&=&- L_{M_1}({\bf x})+L_{M_1}({\bf p})
\end{eqnarray*}
and in particular, that $V_{loc}$ is continuously orbitally differentiable.
Note that $V_{loc}'({\bf x})=-L_{M_1}({\bf x})+r({\bf x})$
for all ${\bf x}\in\overline{\Omega_{\iota/3}}$.

\noindent \underline{\bf VII. Definition of $V_{glob}$ in $A(\Omega)$}

\noindent
For the global part note that
$V_{loc}$ is defined and smooth in $\Omega_{\iota}$ and we have
$V_{loc}'({\bf x})=-L_{M_1}({\bf x})+r({\bf x})$ for all ${\bf x}\in \overline{\Omega_{\iota/3}}$.
The global function $V_{glob}\colon A(\Omega)\setminus \Omega\to \mathbb R$ is defined as the solution of the non-characteristic Cauchy problem
\begin{eqnarray}
\begin{array}{lcl}
  \nabla V_{glob}({\bf x})\cdot {\bf f}({\bf x})&=&-L_{M_1}({\bf x})+r({\bf x})
\mbox{  for }{\bf x}\in A(\Omega)\setminus
\Omega\\
V_{glob}({\bf x})&=&V({\bf x}) \mbox{ for }
{\bf x}\in\Gamma,\end{array}\label{cauchy}
\end{eqnarray} where
$\Gamma=\{{\bf x}\in U\mid d({\bf x})=\iota/3\}$.

 In particular, we can construct the solution by first defining the function $\tau\in C^\sigma(A(\Omega)\setminus \Omega,\mathbb R)$ implicitly by
$$d(S_{\tau}{\bf x})=\iota/3.$$
Since ${\bf x}\in A(\Omega)\setminus \Omega$,  there exists a $\tau$ satisfying the equation, and since $d'({\bf x})<0$ for all ${\bf x}\in \Gamma$, $\tau({\bf x})$ is unique. The function $\tau$ is $C^{\sigma-1}$, since $d$ and $S_\tau$ are. We have
$\tau'({\bf x})
-1$. Then the function
$$V_{glob}({\bf x})=\int_0^{\tau({\bf x})} q(S_t{\bf x})\,dt+V_{loc}(S_{\tau({\bf x})}({\bf x}))$$
with $q({\bf x}):=L_{M_1}({\bf x})-r({\bf x})$ is continuous and orbitally continuously differentiable and satisfies \eqref{cauchy}, noting that $S_{\tau({\bf x})}({\bf x})=S_{\tau(S_\theta{\bf x})}(S_\theta{\bf x})$ for all $\theta\ge 0$. Indeed, for ${\bf x}\in \Gamma$ we have $ V_{glob}({\bf x})=V_{loc}({\bf x})$ and we have
\begin{eqnarray*}
V_{glob}'({\bf x})&=&\frac{d}{d\theta}\left(\int_0^{\tau(S_\theta {\bf x})} q(S_{t+\theta}{\bf x})\,dt+V(S_{\tau(S_\theta {\bf x})}(S_\theta {\bf x}))\right)\bigg|_{\theta=0}\\
&=&\frac{d}{d\theta}\left(\int_\theta^{\tau(S_\theta {\bf x})+\theta} q(S_{t}{\bf x})\,dt+V(S_{\tau( {\bf x})}( {\bf x}))\right)\bigg|_{\theta=0}\\
&=&\left( q(S_{\tau(S_\theta {\bf x})+\theta}{\bf x})(\tau'( {\bf x})+1)
-q(S_\theta {\bf x})\right)\bigg|_{\theta=0}\\
&=&-q({\bf x})
	\end{eqnarray*}
	since $\tau'({\bf x})=-1$.

	Note that we have $V_{glob}({\bf x})=V_{loc}({\bf x})$ for ${\bf x}\in \overline{\Omega_{\iota/3}}\setminus \Omega$, and hence $V_{glob}$ can be extended to a continuous and orbitally continuously differentiable function $V$ on  $A(\Omega)$ satisfying \eqref{eqV} by setting $V_{glob}({\bf x}):=V_{loc}({\bf x})=0$ for all ${\bf x} \in \Omega$.
This proves the theorem.
\end{proof}

\section*{Conclusions}

In this paper we have proven a converse theorem, showing the existence of a contraction metric for an exponentially stable periodic orbit. The metric is defined in its basin of attraction and the bound  on the function $L_M$ is arbitrarily close to the true  exponential rate of attraction.

\begin{appendix}

	\section{Local Lipschitz-continuity of $L_M$}

	In the appendix we prove that the function $L_M$ is locally Lipschitz continuous.

\begin{lemma}\label{Lipschitz} Let ${\bf f}\in C^2(\mathbb R^n,\mathbb R^n)$ and $M\in C^2(\mathbb R^n,\mathbb S^n)$ such that $M({\bf x})$ is positive definite for all ${\bf x}\in \mathbb R^n$.

 Then $L_M$ is locally Lipschitz continuous on $D=\{{\bf x}\in \mathbb R^n\mid {\bf f}({\bf x})\not={\boldsymbol 0}\}$.

\end{lemma}
\begin{proof}
For ${\bf y}\in D$ we define a projection $P_{\bf y}\colon \mathbb R^n\to\mathbb R^n$  onto the $(n-1)$-dimensional space of vectors ${\bf w}\in \mathbb R^n$ with ${\bf f}({\bf y})^TM({\bf y}){\bf w}=0$  by
$$P_{\bf y} {\bf v} ={\bf v} -\frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}{\bf f}({\bf y})$$
for all ${\bf y}\in D$ and all ${\bf v}\in \mathbb R^n$.
Note that indeed
\begin{eqnarray*}
{\bf f}({\bf y})^TM({\bf y})P_{\bf y} {\bf v}
&=&{\bf f}({\bf y})^TM({\bf y}){\bf v} -\frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}
{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})\\
&=&{\boldsymbol 0}.
\end{eqnarray*}

Fix ${\bf x}\in D$ and choose a basis ${\bf v}_1={\bf f}({\bf x}),{\bf v}_2,\ldots,{\bf v}_n$ of $\mathbb R^n$ such that ${\bf v}_i^TM({\bf x}){\bf v}_j=0$ for $i\not=j$.
Choose $\epsilon>0$ such that
\begin{eqnarray}
{\bf f}({\bf y})^TM({\bf x}){\bf f}({\bf x})&\not=&0\label{not0}
\end{eqnarray} holds for all ${\bf y}\in B_\epsilon({\bf x})$; note that for ${\bf y}={\bf x}$ we have ${\bf f}({\bf x})^TM({\bf x}){\bf f}({\bf x})\not=0$.

For ${\bf y}\in B_\epsilon({\bf x})$ we define ${\bf w}_1={\bf f}({\bf y})$ and ${\bf w}_i=P_{\bf y}{\bf v}_i$ for $i=2,\ldots,n$. We show that $({\bf w}_1,\ldots,{\bf w}_n)$ is a basis of $\mathbb R^n$.

Let us first show that ${\bf w}_i\not={\boldsymbol 0}$ for $i=2,\ldots,n$. Assuming the opposite, we have
\begin{eqnarray}
{\bf v}_i&=&\frac{{\bf f}({\bf y})M({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}{\bf f}({\bf y})\label{first_eq}\\
0&=&\frac{{\bf f}({\bf y})M({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}\nonumber
\end{eqnarray}
 multiplying by ${\bf f}({\bf x})^TM({\bf x})$ from the left
as ${\bf f}({\bf x})^TM({\bf x}){\bf f}({\bf y})\not=0$ by \eqref{not0}. This, however, implies by \eqref{first_eq} that ${\bf v}_i={\boldsymbol 0}$ which is a contradiction.
${\bf w}_1\not={\boldsymbol 0}$ follows directly from \eqref{not0}.

We express ${\bf f}({\bf y})=\sum_{j=1}^n\beta_j{\bf v}_j$ and note that multiplying this equation by ${\bf f}({\bf x})^TM({\bf x})$ from the left gives $$0\not={\bf f}({\bf x})^TM({\bf x}){\bf f}({\bf y})=\beta_1 {\bf f}({\bf x})^TM({\bf x}){\bf f}({\bf x})$$
by \eqref{not0}, i.e. in particular $\beta_1\not=0$.

To show that the ${\bf w}_i$ form a basis, we  assume $\sum_{i=1}^n\alpha_i{\bf w}_i={\boldsymbol 0}$. Multiplying this equation by ${\bf f}({\bf y})^TM({\bf y})$ from the left gives
$\alpha_1 {\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})=0$ by the projection property, hence $\alpha_1=0$.

Hence,
\begin{eqnarray*}
{\boldsymbol 0}&=&\sum_{i=2}^n\alpha_i\left[{\bf v}_i-\frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}{\bf f}({\bf y})\right]\\
&=&\sum_{i=2}^n\alpha_i{\bf v}_i-
\sum_{i=2}^n\sum_{j=1}^n \frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}\beta_j{\bf v}_j.
\end{eqnarray*}
Using that ${\bf v}_j$ is a basis, we can conclude that the coefficient in front of ${\bf v}_1$ is zero, namely
$$
\sum_{i=2}^n \frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}\beta_1=0.$$
Since $\beta_1\not=0$, we have $
\sum_{i=2}^n \frac{{\bf f}({\bf y})^TM({\bf y}){\bf v}_i}{{\bf f}({\bf y})^TM({\bf y}){\bf f}({\bf y})}=0$.
Plugging this back in, we obtain $\sum_{i=2}^n\alpha_i{\bf v}_i={\boldsymbol 0}$, which shows $\alpha_2=\ldots=\alpha_n=0$ as the ${\bf v}_i$ are linearly independent.

Now define the matrix-valued function $Q\colon B_\epsilon({\bf x})\to \mathbb R^{n\times n}$ by the columns
$$Q({\bf y})=({\bf w}_1({\bf y}),\ldots,{\bf w}_n({\bf y})).$$
Note that
$Q\in C^2(B_\epsilon({\bf x}),\mathbb R^{n\times n})$ due to the smoothness of ${\bf f}$ and $M$, and $Q$ is invertible for every ${\bf y}$. We have ${\bf w}^TM({\bf y}){\bf f}({\bf y})=0$ if and only if ${\bf w}\in \operatorname{span} ({\bf w}_2({\bf y}),\ldots,{\bf w}_n({\bf y})$, which in turn is equivalent to
${\bf u}\in \operatorname{span}({\bf e}_2,\ldots,{\bf e}_n)=:E_{n-1}$, where ${\bf u}=Q({\bf y})^{-1}{\bf w}$ and ${\bf e}_1,\ldots,{\bf e}_n$ denotes the standard basis in $\mathbb R^n$.

Now we write
\begin{eqnarray*}
\lefteqn{L_M({\bf y})}\\
&=&\max_{{\bf w}^TM({\bf y}){\bf w}=1,{\bf w}^TM({\bf y}){\bf f}({\bf y})=0}
\frac{1}{2}{\bf w}^T\left[M({\bf y})D{\bf f}({\bf y})+D{\bf f}({\bf y})^TM({\bf y})+M'({\bf y})\right]{\bf w}\\
&=&
\max_{{\bf u}^TQ({\bf y})^TM({\bf y})Q({\bf y}){\bf u}=1,{\bf u}\in E_{n-1}}
\frac{1}{2}{\bf u}^TQ({\bf y})^T\\
&&\hspace{1cm}\left[M({\bf y})D{\bf f}({\bf y})+D{\bf f}({\bf y})^TM({\bf y})+M'({\bf y})\right]Q({\bf y}){\bf u}.
\end{eqnarray*}
Denoting by $[A]_{n-1}\in \mathbb S^{n-1}$ the lower-right square $(n-1)$ matrix of $A\in \mathbb S^{n}$ and
 with ${\bf u}=\left(\begin{array}{l}0\\\widetilde{\bf u} \end{array}\right)$, where $\widetilde{\bf u}\in\mathbb R^{n-1}$ we get
\begin{eqnarray*}
L_M({\bf y})
&=&
\max_{\widetilde{\bf u}^T[Q({\bf y})^TM({\bf y})Q({\bf y})]_{n-1}\widetilde{\bf u}=1,\widetilde{\bf u}\in \mathbb R^{n-1}}
\frac{1}{2}\widetilde{\bf u}^T\bigg[Q({\bf y})^T\big[M({\bf y})D{\bf f}({\bf y})\\
&&\hspace{1cm}
+D{\bf f}({\bf y})^TM({\bf y})+M'({\bf y})\big]Q({\bf y})\bigg]_{n-1}\widetilde{\bf u}.
\end{eqnarray*}
Now denote by $\operatorname{Chol}(A)$ the unique Cholesky decomposition of the symmetric, positive definite matrix $A\in \mathbb S^{n-1}$, such that $\operatorname{Chol}(A)$ is an invertible, upper triangular matrix with $\operatorname{Chol}(A)^T\operatorname{Chol}(A)=A$. Denoting $C({\bf y}):=\operatorname{Chol}([Q({\bf y})^TM({\bf y})Q({\bf y})]_{n-1})\in \mathbb R^{(n-1)\times (n-1)}$ and $\widetilde{\bf v}=C({\bf y})\widetilde{\bf u}\in\mathbb R^{n-1}$ we have
\begin{eqnarray*}
L_M({\bf y})
&=&
\max_{\|\widetilde{\bf v}\|=1,\widetilde{\bf v}\in \mathbb R^{n-1}}
\frac{1}{2}\widetilde{\bf v}^T(C^{-1}({\bf y}))^T\\
&&\hspace{0.8cm}\left[Q({\bf y})^T\left[M({\bf y})D{\bf f}({\bf y})+D{\bf f}({\bf y})^TM({\bf y})+M'({\bf y})\right]Q({\bf y})\right]_{n-1}C^{-1}({\bf y})\widetilde{\bf v}\\
&=&
\max_{\|\widetilde{\bf v}\|=1,\widetilde{\bf v}\in \mathbb R^{n-1}}
\widetilde{\bf v}^TH({\bf y})\widetilde{\bf v}\\
&=&\lambda_{max}(H({\bf y}))
\end{eqnarray*}
where $H({\bf y})\in \mathbb S^{n-1}$ is defined by
\begin{eqnarray*}H({\bf y})&=&\frac{1}{2}(C^{-1}({\bf y}))^T\left[Q({\bf y})^T\left[M({\bf y})D{\bf f}({\bf y})+D{\bf f}({\bf y})^TM({\bf y})+M'({\bf y})\right]Q({\bf y})\right]_{n-1}\\
&&C^{-1}({\bf y}).
\end{eqnarray*}

The function ${\bf y}\to H({\bf y})$ is continuously differentiable as the Cholesky decomposition, the inverse, the operation $[\cdot]_{n-1}$, $Q$, $M$, $D{\bf f}$ and $M'$ are continuously differentiable by the assumptions. Hence, the function $H({\bf y})$ is locally Lipschitz-continuous. The function $\lambda_{max}$ is globally Lipschitz-continuous, hence, $L_M$ is locally Lipschitz-continuous.
\end{proof}
\end{appendix}

{\small
	\bibliographystyle{plain}