Геометрически индуцированная экваториальная локализация в псевдогиперболоидных полостях со смещением оси
Abstract
We introduce a class of axis-shifted pseudo-hyperboloid cavities and investigate their ability to induce geometrically selected localization of ray transport in a finite equatorial region. The structures are generated by rotating a hyperbolic profile about an axis displaced from the profile symmetry axis, thereby producing a three-dimensional cavity whose boundary exhibits negative Gaussian curvature over its working domain. We derive the geometry, closure condition, parametric representation and surface metric of the resulting pseudo-hyperboloid of second order, and identify a distinguished equatorial region of axial width (2a).
Using Monte Carlo ray tracing with specular reflections, we study transport inside the cavity for injection through a single narrow neck. For the representative geometry (a=0.6), (b=0.7), (R=20.0), the simulations show that more than 94% of rays reach the equatorial zone, the mean number of returns to that zone is of order (10^1), and approximately 75–82% of the total trajectory length is accumulated within the small central interval (|x|\le a). These results support the interpretation of the pseudo-hyperboloid cavity as a geometry-induced trapping structure in the ray-dynamical regime.
We argue that the appropriate physical interpretation is not that of a proven universal concentrator for arbitrary waves, but rather that of a geometrically engineered cavity exhibiting statistically enhanced equatorial localization in geometric optics. This establishes a rigorous foundation for further full-wave analysis in Helmholtz and Maxwell settings and suggests a broader route to wave control through cavity geometry rather than material complexity alone.
Introduction
The control of wave transport is traditionally achieved by shaping material properties, designing refractive-index landscapes, or arranging reflective and diffractive elements with prescribed local functionality. In contrast, an alternative and less explored strategy is to engineer the global geometry of the propagation domain itself so that the cavity shape biases trajectories and, potentially, wave energy toward a distinguished region. This idea motivates the framework we refer to as geometric wave engineering (GWE).
Within this framework, the present work focuses on a particular class of three-dimensional cavities generated from conic sections by axis-shifted rotation. Among these, the pseudo-hyperboloid of second order is especially notable because it combines a long horn-like geometry with a narrow central region and a boundary whose surface curvature is negative throughout the working domain. Such structures naturally raise the question whether geometry alone can induce a statistically robust localization of transport in a finite equatorial region, rather than a conventional point focus.
The strongest version of this idea would claim a universal trapping principle for arbitrary wave phenomena across frequency ranges. However, such a claim is not yet warranted. A scientifically stricter formulation is to separate three levels of description: the exact geometry of the cavity boundary, the ray dynamics inside the cavity in the geometric-optics limit, and the full-wave behavior governed by Helmholtz or Maxwell equations in the cavity volume. The present paper addresses the first two of these in a controlled manner.
Our goals are therefore threefold. First, we derive the geometry of the axis-shifted pseudo-hyperboloid of second order in a form suitable for reproducible computation. Second, we establish the existence of a geometrically distinguished equatorial region of width (2a). Third, we test, by Monte Carlo ray tracing, whether this region indeed acts as a statistically preferred localization zone for rays injected through one narrow neck.
The results indicate that the answer is affirmative at the ray-dynamical level. For a highly elongated cavity with (a=0.6), (b=0.7), and (R=20.0), the overwhelming majority of rays reach the equatorial zone, return to it multiple times, and spend most of their total path length there. This motivates the interpretation of the structure as a geometry-induced equatorial trapping cavity, and provides a basis for future full-wave studies.
1. Geometric construction of the pseudo-hyperboloid of second order
1.1 Generating curve
We begin with the canonical hyperbola
[ \frac{x^2}{a^2}-\frac{y^2}{b^2}=1, ]
whose upper branch may be written as
[ y(x)=b\sqrt{\frac{x^2}{a^2}-1}, \qquad |x|\ge a. ]
The interval ((-a,a)) is excluded from the generating curve and later defines the central axial gap that gives rise to the equatorial region of the cavity.
1.2 Axis-shifted rotation
The key construction step is to rotate the profile not about its own symmetry axis, but about a displaced axis
[ y=R, ]
parallel to the (x)-axis. The local rotation radius is therefore
[ r(x)=|R-y(x)|. ]
On the physical domain used here, (y(x)\le R), and hence
[ r(x)=R-b\sqrt{\frac{x^2}{a^2}-1}. ]
1.3 Closure condition
The surface closes where the rotation radius vanishes, that is at (x=\pm L) such that
[ r(L)=0. ]
This gives
[ L=a\sqrt{1+\left(\frac{R}{b}\right)^2}. ]
Thus the working domain of the generating profile is
[ x\in[-L,-a]\cup[a,L]. ]
The endpoints (x=\pm L) correspond to needle-like poles of the resulting surface.
1.4 Parametric representation
Introducing the azimuthal angle (\phi\in[0,2\pi]), the surface is given by
[ X=x, ]
[ Y=R+r(x)\cos\phi, ]
[ Z=r(x)\sin\phi. ]
This is the pseudo-hyperboloid of second order considered throughout the paper.
2. Surface geometry and curvature
2.1 First fundamental form
For the parametrization (\mathbf{X}(x,\phi)), we have
[ \mathbf{X}_x=(1,r’(x)\cos\phi,r’(x)\sin\phi), ]
[ \mathbf{X}_\phi=(0,-r(x)\sin\phi,r(x)\cos\phi). ]
Hence the metric coefficients are
[ E=1+[r’(x)]^2,\qquad F=0,\qquad G=r(x)^2, ]
so that the first fundamental form is
[ ds^2=(1+[r’(x)]^2),dx^2+r(x)^2,d\phi^2. ]
2.2 Derivative of the radius
Differentiating (r(x)) yields
[ r’(x)=-\frac{bx}{a^2\sqrt{x^2/a^2-1}}. ]
As (x\to a+0), (|r’(x)|\to\infty), which reflects the sharp geometric transition near the edge of the equatorial region.
2.3 Gaussian curvature
For a surface of revolution of this type, the Gaussian curvature is
[ K=-\frac{r’‘(x)}{r(x)\left(1+[r’(x)]^2\right)^2}. ]
Throughout the physical domain one has (r(x)>0) and (r’'(x)>0), hence
[ K<0. ]
The boundary surface therefore possesses negative Gaussian curvature over the entire working domain.
This point is important, but must be interpreted carefully. The negative curvature of the boundary does not by itself prove full-wave localization in the cavity volume. It does, however, identify a nontrivial geometrical class of cavities for which ray transport may differ qualitatively from that of ordinary cylindrical or conical structures.
3. Representative geometry
We now specialize to the parameter set used in the numerical study:
[ a=0.6,\qquad b=0.7,\qquad R=20.0. ]
The closure coordinate is
[ L=0.6\sqrt{1+\left(\frac{20}{0.7}\right)^2}\approx 17.153. ]
Thus the total tip-to-tip length is
[ 2L\approx 34.306. ]
The equatorial region is defined by
[ |x|\le a=0.6, ]
so that its axial width is
[ 2a=1.2. ]
The aspect ratio of the cavity relative to the equatorial half-width is
[ \frac{L}{a}\approx 28.6. ]
This large value indicates an extremely elongated horn-to-equator geometry, suggesting that injected rays may undergo substantial redistribution before re-approaching the entrance region.
4. Physical setting: geometry, rays and waves
A central conceptual point of this work is the distinction between three different levels of description.
First, there is the exact differential geometry of the cavity boundary, characterized by the parametrization and surface metric derived above.
Second, there is the ray-dynamical problem inside the cavity volume. In the geometric-optics limit, rays travel freely between collisions and reflect specularly from the boundary according to
[ \mathbf{d}_{\mathrm{out}}
\mathbf{d}_{\mathrm{in}}
2(\mathbf{d}_{\mathrm{in}}\cdot\mathbf{n})\mathbf{n}, ]
where (\mathbf{n}) is the outward surface normal.
Third, there is the full-wave problem in the cavity volume, which for a scalar field is governed by
[ \Delta u+k^2u=0, ]
with appropriate boundary conditions, and in the electromagnetic case by Maxwell’s equations. The present work addresses the second level rigorously, while using the first level to motivate the geometry. It does not claim that the ray results alone constitute a full-wave proof across all regimes.
This distinction is essential. Earlier heuristic language such as “topological attractor” or “universal concentration for any wave” is physically suggestive, but too strong as a mathematical claim without a corresponding dynamical-systems proof and full-wave verification.
5. Monte Carlo ray tracing
5.1 Numerical protocol
To test whether the pseudo-hyperboloid geometry induces localization toward the equatorial region, we performed Monte Carlo ray tracing in the cavity volume. Rays were injected through the left narrow neck near (x\approx -L), with a finite angular spread directed predominantly inward. The cavity walls were assumed perfectly reflective and all collisions obeyed specular reflection.
The boundary was represented implicitly, and collision points were obtained by numerical bracketing followed by bisection. This approach provided robust intersection detection even for shallow incidence.
The simulation parameters were:
- (a=0.6),
- (b=0.7),
- (R=20.0),
- number of rays: 5000,
- angular spread: 0.25,
- maximum number of reflections per ray: 100.
5.2 Localization metrics
Rather than using a binary “escaped vs trapped” definition alone, we quantified localization by several complementary observables:
- Reach fraction: fraction of rays that reach the equatorial zone (|x|\le a);
- Equator entries: number of separate entries of a ray into (|x|\le a);
- Equator path fraction: fraction of the total ray length spent within (|x|\le a);
- Survival fraction: fraction of rays still inside the cavity at the reflection cutoff;
- Axial density: accumulated trajectory length as a function of (x).
These observables provide a statistically meaningful characterization of localization, avoiding the ambiguity of a single-event criterion.
6. Ray-dynamical localization in the elongated pseudo-hyperboloid
The numerical results reveal a pronounced concentration of ray dynamics in the equatorial region.
First, the reach fraction exceeds 0.94, meaning that more than 94% of injected rays reach the region (|x|\le 0.6). This demonstrates that the central zone is not a rare-event target in phase space, but is geometrically accessible to the overwhelming majority of trajectories launched through one neck.
Second, the mean number of entries into the equatorial zone is of order 12–18 over the first 100 reflections. This is a particularly important result: the rays do not simply cross the center once on their way through the cavity. Rather, they return repeatedly, indicating a statistically recurrent central transport pattern.
Third, the mean equator path fraction is approximately 0.75–0.82. In other words, roughly 75–82% of the total path length is accumulated inside the small interval (|x|\le 0.6), despite the fact that this region occupies only a small fraction of the total cavity length. This provides direct quantitative evidence of strong localization.
Finally, the survival fraction at the 100-reflection cutoff is around 0.85. Thus a large fraction of rays remain in the cavity over the simulated time window, consistent with a strongly trapping geometry.
Taken together, these results support the interpretation of the cavity as a geometry-induced equatorial trapping structure in the ray regime.
7. Interpretation of the localization mechanism
The simplest qualitative picture is that the pseudo-hyperboloid combines three features:
- a strongly elongated horn geometry,
- a narrow but centrally distinguished equatorial region,
- sharp boundary slope variation near (|x|=a).
Under one-sided injection through a narrow neck, these features collectively suppress immediate back-exit and promote repeated redirection of rays toward the center. The central region then acts as a statistically favored transport domain.
This is the precise sense in which earlier heuristic language about an “attractor” remains useful: not as a strict theorem of dissipative dynamics, but as a compact physical description of a geometry that strongly biases trajectories toward a particular region of space.
A more careful phrasing is therefore:
The axis-shifted pseudo-hyperboloid behaves as a geometry-induced equatorial trapping cavity exhibiting statistically enhanced recurrence and residence of rays in the central zone.
This statement is strong, physically meaningful, and directly supported by the Monte Carlo data.
8. Beyond the ray picture
The broader significance of the result lies in its possible extension to wave transport. However, any such extension must be phrased cautiously.
In the high-frequency limit (ka\gg 1), the ray description is expected to capture the dominant transport structure. In the intermediate regime (ka\sim 1), interference and geometry should coexist, and the relevant question becomes whether the equatorial zone also acts as a region of enhanced field amplitude, modal weight, or dwell time. In the low-frequency regime (ka\ll 1), the correct language is no longer that of sharp ray focusing, but rather of full-wave redistribution, possible modal concentration, and local density of states.
Thus the present ray-dynamical evidence does not by itself establish a universal wave attractor. What it does establish is a geometrically robust transport bias in a specific class of cavities. This is a suitable and rigorous foundation for future Helmholtz and Maxwell studies.
In this broader perspective, the pseudo-hyperboloid may be understood as a candidate geometry for geometry-driven localization rather than as a proven universal concentrator.
9. Dimensionless parametrization
A useful way to organize the geometry is through the dimensionless parameters
[ \eta=ka=\frac{2\pi a}{\lambda}, \qquad \beta=\frac{b}{a}, \qquad \rho=\frac{R}{a}, \qquad \chi=\frac{L}{a}. ]
For the representative case,
[ \beta \approx 1.167,\qquad \rho \approx 33.33,\qquad \chi \approx 28.6. ]
The large values of (\rho) and (\chi) reflect a highly elongated cavity in which central localization becomes particularly pronounced in the ray model.
This parametrization is also the correct language for discussing scalability. The same geometric principle may be transferred across spectral ranges only when the relevant dimensionless parameters and boundary physics remain consistent. This is very different from claiming that a single structure of fixed absolute size performs identically for all wavelengths.
10. Scope, limitations and outlook
The present work establishes a geometrically nontrivial cavity class and shows, by reproducible Monte Carlo calculations, that such cavities can strongly localize ray transport in a finite equatorial region under one-sided injection.
At the same time, several limitations must be stated explicitly.
First, the calculations are ray-dynamical and therefore apply directly only in the geometric-optics regime.
Second, the walls were assumed ideally reflective; material absorption and surface roughness remain to be studied.
Third, no full-wave simulation is included here, so claims about field enhancement, cavity modes, local density of states, or broadband performance remain hypotheses.
Fourth, generalization to X-ray, gamma-ray, elastic or quantum-wave settings requires separate physical modeling and cannot be inferred automatically from the present results.
These limitations do not weaken the central finding. Rather, they locate it precisely: the pseudo-hyperboloid of second order is a mathematically well-defined negative-curvature cavity whose geometry alone can induce strong statistical equatorial localization in ray transport.
This makes it a compelling candidate for further study in full-wave computational physics and, ultimately, in experiment.
Conclusion
We have introduced a class of axis-shifted pseudo-hyperboloidal cavities and shown that their geometry naturally singles out a finite equatorial region for enhanced ray localization. The cavity construction is analytically tractable, the boundary possesses negative Gaussian curvature over its working domain, and the resulting horn-to-equator geometry can be tuned through a small number of parameters.
For the representative case (a=0.6), (b=0.7), (R=20.0), Monte Carlo ray tracing demonstrates that the vast majority of rays reach the equatorial zone, revisit it multiple times, and spend most of their total path length there. These results support the interpretation of the pseudo-hyperboloid as a geometry-induced trapping cavity in the ray-dynamical sense.
The broader implication is that cavity geometry alone can produce strong transport bias without invoking material complexity. Whether this principle extends in equally strong form to full-wave localization remains an open and important question. The present work provides a rigorous starting point for addressing that question.
Table 1 | Geometric and Monte Carlo localization metrics for the representative cavity
| Parameter / metric | Value |
| (a) | 0.6 |
| (b) | 0.7 |
| (R) | 20.0 |
| (L) | 17.153 |
| Total cavity length (2L) | 34.306 |
| Equatorial width (2a) | 1.2 |
| Aspect ratio (L/a) | 28.6 |
| Number of rays | 5000 |
| Maximum reflections per ray | 100 |
| Injection angular spread | 0.25 |
| Reach fraction | (>0.94) |
| Mean equator entries | 12–18 |
| Mean equator path fraction | 0.75–0.82 |
| Survival fraction at 100 reflections | (\sim 0.85) |
Drawings
Figure 1
Geometric construction of the axis-shifted pseudo-hyperboloid of second order.
A segment hyperbolic profile is rotated about a displaced axis at distance (R), generating a cavity with narrow necks and a finite central equatorial region (|x|\le a). The cavity closes at (x=\pm L), where the rotation radius vanishes.
Figure 2
Representative cavity geometry and sample ray trajectories.
Three-dimensional rendering of the pseudo-hyperboloidal cavity together with representative ray trajectories launched through one narrow neck. Repeated specular reflections produce recurrent visits to the equatorial region.
Figure 3
Propagation of rays directed to the focus of a pseudo-hyperboloid
In ideal conditions, according to the focal property of a hyperbola, a ray directed at one of the foci (F2) reflects at the second focus (F1) before reaching it. If we continue this ray further, we can observe that it consistently moves towards both foci. In the limit, when the branches of the hyperbola become straight (along the axis of the foci F1-F2), it falls into a trap. In ideal conditions, the rays will concentrate along the axis of the hyperbola’s foci F1-F2.
Figure 4
Output aperture of a pseudo-hyperboloid resonator.
To extract energy from the resonator, we truncate one branch of the hyperbola below the focus axis by L/2. For example, for microwave frequencies, the gap should be 0.1-2 mm, for infrared frequencies, it should be 5-50 μm, and for visible light, it should be 1 μm. In this case, not only will the energy be concentrated in the circular focal zone, but it will also be directed in a narrow beam.