Global Convergence of the Return Dynamics in the Class $\mathcal{O}_{C}$

Let $C\subset\mathbb{R}^{N}$ be a compact convex set with non-empty interior and boundary $\partial C$ of class $C^{2}$. We denote by

the outward unit normal vector to $\partial C$ at $c\in\partial C$.

Let $\Omega\subset\mathbb{R}^{N}$ be an open set containing $C$. We consider domains $\Omega\in\mathcal{O}_{C}$ as introduced in the Introduction [1]: for almost every point $x\in\partial\Omega$ where the inward unit normal $n(x)$ exists, the half-line

intersects the convex set $C$.

This condition guarantees that inward normal rays from the boundary $\partial\Omega$ reach the convex core $C$.

We refer to domains $\Omega\in\mathcal{O}_{C}$ as admissible domains.

The geometry of the region between $\partial C$ and $\partial\Omega$ is encoded by the thickness function

defined by

Geometrically, $d(c)$ measures the distance from the boundary $\partial C$ to $\partial\Omega$ along the outward normal direction.

The thickness function provides a parametrization of the outer boundary through the radial map

defined by

Let $x\in\partial\Omega$ and denote by $n(x)$ the inward unit normal to $\partial\Omega$ at $x$. Let $t(x)$ be the smallest non-negative number such that

The reciprocal map is then defined by

The composition of the radial and reciprocal maps defines the return map

Starting from a point $c\in\partial C$, the map first moves to the outer boundary $\partial\Omega$ along the outward normal direction via $\Phi$, and then returns to $\partial C$ along the inward normal to $\partial\Omega$ via $\pi$. The resulting transformation generates a discrete dynamical system

on the hypersurface $\partial C$. The purpose of the present paper is to analyze the global asymptotic behavior of this dynamical system.

Let the thickness function

be of class $C^{2}$. The return map

defined by $F=\pi\circ\Phi$ admits the expansion

where the remainder term satisfies

Equation (3.3) shows that the leading displacement of a point under the return map is proportional to the negative gradient of the thickness function.

The expansion (3.3) reveals that the return dynamics behaves, to leading order, like a gradient descent iteration on the hypersurface $\partial C$. Indeed, neglecting the higher-order remainder term, the iteration

formally reduces to

Thus the return map may be interpreted as a variable-step gradient descent for the thickness function, with step size

This gradient structure immediately suggests the existence of a natural Lyapunov function for the dynamics. Define

As will be shown in the next section, the function $V$ decreases along trajectories of the return map. This observation provides the key tool for analyzing the global asymptotic behavior of the iteration.

Let $d:\partial C\to\mathbb{R}_{+}$ be the thickness function introduced in Section 2, assumed of class $C^{2}$. We define the energy functional

This quantity measures the squared thickness between the two boundaries along the outward normal direction.

Let $(c_{k})$ be a trajectory of the return map

Using the first-order expansion (3.3), we obtain the displacement

A Taylor expansion of the thickness function gives

where $\xi_{k}$ lies on the geodesic segment between $c_{k}$ and $c_{k+1}$.

Substituting (4.3) into (4.4) and using the remainder estimate (3.4), we obtain, after squaring and collecting terms,

More precisely, under the assumption that $\nabla_{\partial C}d$ is Lipschitz and that the curvature of $\partial C$ is bounded, one obtains the following uniform estimate.

Thus the Lyapunov function decreases whenever the gradient of the thickness function is non-zero, provided the gradient is sufficiently small.

The inequality (4.6) implies that the Lyapunov function $V(c_{k})$ is nonincreasing to leading order, and the energy cannot increase when the gradient is nonzero.

Since $\partial C$ is compact, the trajectory $(c_{k})$ admits accumulation points. Moreover, any accumulation point must satisfy

We assume throughout this section that the thickness function

satisfies the following properties.

1. 1.

   $d\in C^{2}(\partial C)$.

2. 2.

   There exist constants $d_{\min},d_{\max}>0$ such that

   $$0<d_{\min}\leq d(c)\leq d_{\max}\quad\text{for all }c\in\partial C.$$

   (5.2)

3. 3.

   The tangential gradient $\nabla_{\partial C}d$ is Lipschitz.

4. 4.

   The critical points of $d$ are nondegenerate (Morse condition) [3].

5. 5.

   The gradient of the thickness function is uniformly small:

   $$\|\nabla_{\partial C}d\|_{L^{\infty}(\partial C)}\leq\varepsilon_{0},$$

   (5.3)

   for $\varepsilon_{0}>0$ sufficiently small (to be chosen below).

From Proposition 4.1 we have

with constants $a,b>0$.

Under the small-gradient assumption (5.3), we have

Choosing $\varepsilon_{0}$ sufficiently small so that $b\varepsilon_{0}^{2}<a$, we obtain

for some constant $\eta>0$ and for all $c\in\partial C$.

Let $(c_{k})$ be a trajectory of the return map $c_{k+1}=F(c_{k})$. Summing (5.6) over $k$ yields

Therefore

Let

Since $d$ is Morse on the compact manifold $\partial C$, this set is finite [3].

Let $U$ be an open neighborhood of $\mathrm{Crit}(d)$. Because $|\nabla_{\partial C}d|$ is continuous and nonzero on the compact set $\partial C\setminus U$, there exists $\varepsilon_{U}>0$ such that

Using (5.6), we obtain a uniform decrease: there exists $\gamma_{U}>0$ such that

Since $V$ is bounded below and strictly decreases by at least $\gamma_{U}$ whenever the trajectory lies in $\partial C\setminus U$, the sequence $(c_{k})$ can visit $\partial C\setminus U$ only finitely many times. Therefore every trajectory eventually enters the neighborhood $U$ and remains in $U$ for all sufficiently large $k$.

Since $\partial C$ is compact, the sequence $(c_{k})$ admits accumulation points. Let $c^{*}$ be such a point. By continuity of $\nabla_{\partial C}d$ and (5.8), we obtain

Thus every accumulation point is a critical point of the thickness function.

We now prove that the entire sequence converges.

Recall that the return map is defined by $F=\pi\circ\Phi$. A point $c^{*}\in\partial C$ is a fixed point of the dynamics if $F(c^{*})=c^{*}$.

Thus the equilibrium configurations of the dynamical system coincide with the critical points of the thickness function.

We now study the behavior of the return map near a critical point. Let $c^{*}$ be a critical point of $d$, i.e.,

Consider a small perturbation

Expanding the gradient of $d$ gives

Moreover, since (6.2) holds, the thickness function satisfies

Using the first-order expansion (3.3), we obtain

Substituting (6.4) and (6.5) into (6.6), we obtain

Therefore,

Hence the Jacobian matrix of the return map at $c^{*}$ is

Let $\lambda_{1},\ldots,\lambda_{N-1}$ be the eigenvalues of the tangential Hessian $\nabla_{\partial C}^{2}d(c^{*})$. The eigenvalues of the linearized return map are therefore

The classification follows from the linearization (6.9) and standard results in dynamical systems. For hyperbolic fixed points, the stable manifold theorem [5] guarantees the existence of local stable and unstable manifolds.

Let $c^{*}$ be a critical point of the thickness function. If $c^{*}$ is attracting, the linearization satisfies

where $\rho$ denotes the spectral radius. Since the Hessian of $d$ at a nondegenerate critical point is invertible, the linearization has no eigenvalue equal to $1$, and the fixed point is isolated.

Moreover, if $|\mu_{i}|\neq 1$ for all $i$, the fixed point is hyperbolic and the stable manifold theorem [5] applies. In that case there exists a neighborhood $U\subset\partial C$ such that every trajectory starting in $U$ converges to $c^{*}$.

The basin of attraction of $c^{*}$ is defined as

Thus the return dynamics induces a decomposition of $\partial C$ into basins of attraction associated with the critical points of the thickness function.

The geometry of the region between the convex core $C$ and the outer boundary $\partial\Omega$ is encoded by the thickness function

This function defines a scalar landscape on the manifold $\partial C$. The gradient $\nabla_{\partial C}d$ describes the directions along which the thickness increases most rapidly, while the critical points

correspond to locations where the thickness is locally extremal.

The return map arises from the geometric round-trip

Starting from a point $c\in\partial C$, the radial map $\Phi$ moves along the outward normal direction until reaching the boundary $\partial\Omega$. The reciprocal map $\pi$ then follows the inward normal of $\partial\Omega$ until returning to the convex core.

Although the resulting transformation acts entirely on $\partial C$, its generator lies in the surrounding geometry of the domain $\Omega$. The dynamics on $\partial C$ is therefore induced by the geometry of the region between the two boundaries.

The first-order expansion (3.3) of the return map shows that

where the remainder term is of higher order in $\|\nabla_{\partial C}d(c)\|$.

Thus the leading behavior of the dynamics is governed by the negative gradient of the thickness function. The return dynamics therefore acts as a variable-step gradient descent whose step size is proportional to the local thickness $d(c)$ on the manifold $\partial C$.

Moreover, the Lyapunov estimate (5.6) implies that the energy

decreases along trajectories, confirming that the dynamics is gradient-like.

The global convergence result (Theorem 5.2) shows that, under suitable regularity and nondegeneracy assumptions, every trajectory converges to a critical point of the thickness function.

The stability analysis (Theorem 6.2) further shows that the nature of each equilibrium is determined by the Hessian of $d$ at the corresponding point. In particular:

• •

  local minima of $d$ correspond to attracting configurations,

• •

  saddle points generate unstable directions,

• •

  local maxima correspond to repelling configurations.

In this way, the geometry of the domain $\Omega$ organizes the phase space $\partial C$ through the structure of the thickness landscape. The asymptotic behavior of the return dynamics is therefore fully determined by the geometric properties of the region between $\partial C$ and $\partial\Omega$.

To illustrate that the hypotheses are nonempty, consider a simple setting where $C$ is the unit sphere in $\mathbb{R}^{3}$ and $\Omega$ is a slightly deformed sphere given in radial coordinates by

with $\varepsilon>0$ sufficiently small and $f:S^{2}\to\mathbb{R}$ a smooth Morse function with nondegenerate critical points [3].

In this setting, the thickness function satisfies

uniformly on $S^{2}$. In particular,

so that the gradient of $d$ is of order $\varepsilon$. Therefore the small-gradient condition (5.3) is satisfied for sufficiently small $\varepsilon$.

Moreover, the critical points of $d$ coincide with those of $f$, and the nondegeneracy of $f$ implies that $d$ satisfies the Morse condition.

By choosing $f$ with a single nondegenerate minimum (for instance, a localized bump), the thickness function admits a unique global minimizer. In this case, Theorem 5.2 implies that every trajectory converges to this point, yielding a global attractor.

Furthermore, the return map $F$ is a small perturbation of the identity, with leading behavior governed by the gradient structure (3.3). This example shows that the hypotheses of the convergence theorem are realizable and that the dynamics indeed drives the system toward the thinnest region of the shell.

The geometric mechanism uncovered here suggests several directions for further investigation.

• •

  The classification of dynamical systems generated by domains in the class $\mathcal{O}_{C}$ [1].

• •

  The study of inverse problems relating the geometry of $\Omega$ to the structure of the induced return map.

• •

  Extensions of the round-trip construction to other geometric settings where boundary interactions generate hidden dynamical systems.

These directions indicate that the interaction between geometry and induced dynamics provides a natural framework for analyzing hidden geometric transformations through their associated dynamical systems.