The Ptolemy-Alhazen problem with source at infinity

1. Problem Definition

The paper addresses a specific variation of the classical Ptolemy-Alhazen problem, which seeks to find the point of reflection on a spherical mirror given a light source and an observer.

• Classical Context: Traditionally, both the source ($z_1$) and the observer ($z_2$) are at finite distances. The solution involves finding a point $w$ on the mirror boundary such that the angle of incidence equals the angle of reflection.
• The Specific Variation: This paper investigates the case where the light source is at infinity. In this scenario, light rays arrive as parallel beams (e.g., from a distant star or a satellite signal) rather than diverging from a finite point.
• Motivation: The problem was highlighted by V.H. Almeida Jr. in the context of astrophysics applications (specifically GPS signal transmission and radar scattering), where existing numerical methods showed instability or systematic discrepancies when the source distance is effectively infinite.

2. Mathematical Framework and Methodology

The authors utilize tools from geometric function theory, specifically the triangular ratio metric, and algebraic geometry to derive an exact analytical solution.

A. The Triangular Ratio Metric

The problem is framed using the triangular ratio metric $s_G(z_1, z_2)$ for a domain $G$. For a unit disk $B^2$, the reflection point $w$ on the boundary $\partial B^2$ maximizes the ratio:
$$s_{B^2}(z_1, z_2) = \frac{|z_1 - z_2|}{|z_1 - w| + |w - z_2|}$$
Geometrically, this corresponds to finding an ellipse with foci $z_1$ and $z_2$ that is internally tangent to the unit circle at $w$. The metric value equals the eccentricity of this maximal ellipse.

B. Derivation for Source at Infinity

When the source is at infinity, the light rays are parallel. The authors assume rays arrive parallel to the real axis from the right.

• Setup: Let the observer be at $f = re^{i\theta}$ (where $r > 1$) and the reflection point be $w$ on the unit circle ($|w|=1$).
• Geometric Property: The reflection point $w$ is the point of tangency between the unit circle and a parabola whose focus is the observer $f$ and whose directrix is perpendicular to the incoming light rays.
• Algebraic Derivation: By applying the reflection property (equality of angles) and utilizing complex conjugation properties ($w\bar{w}=1$), the authors derive a condition where a specific complex expression must be real. This leads to a quartic (fourth-degree) polynomial equation.

3. Key Results

A. The Quartic Equation (Theorem 3.1)

The primary result is an explicit fourth-degree polynomial equation that determines the reflection point $w$:
$$re^{-i\theta}w^4 - w^3 + w - re^{i\theta} = 0$$

• Solution: The roots of this equation can be solved using Cardano's formula (standard symbolic computation).
• Selection: Among the four roots, the physically relevant solution is the one lying on the unit circle within the specific angular range corresponding to the observer's position ($0 \le \phi \le \pi/2$).

B. Properties of the Roots (Theorem 3.4)

The authors prove that for $r > 1$, all four roots of the derived quartic equation lie strictly on the unit circle.

• Proof Strategy: They transform the complex quartic equation into a real quartic equation with real coefficients via a Möbius transformation ($z = i\frac{1+w}{1-w}$).
• Stability Analysis: By analyzing the discriminant ($\Delta$), the quantity $P$, and the quantity $D$ of the resulting real polynomial, they demonstrate that $\Delta > 0$, $P < 0$, and $D < 0$. These conditions guarantee that the equation has four distinct real roots, confirming that all four complex roots of the original equation lie on the unit circle.

C. The Envelope of Directrices (Section 4)

The paper further explores the geometric locus of the directrices of the parabolas associated with varying light directions.

• Result: The envelope of the family of directrices forms a Limaçon of Pascal.
• Equation: The envelope is given by $z^2\bar{z}^2 - 2(a^2+2)z\bar{z} + 4a(z+\bar{z}) + a^4 - 4a^2 = 0$.
• Shape: For the parameters considered, this curve contains an internal loop, described as a "tear drop". This provides a deeper geometric understanding of the reflection geometry when the source is at infinity.

4. Significance and Applications

• Numerical Stability: Previous methods for solving the Ptolemy-Alhazen problem (such as iterative algorithms or other formulations) suffered from numerical instability when the source distance was very large (the "one-finite" case). This paper provides a closed-form algebraic solution that avoids these instabilities.
• Astrophysics and Engineering: The solution is directly applicable to:
  • GPS Signal Transmission: Modeling signal reflection off the Earth's surface (modeled as a sphere) for satellite-to-satellite or satellite-to-ground links.
  • Radar and Acoustics: Solving scattering problems where sources are effectively at infinity.
  • Computer Vision: Enhancing drone vision systems that rely on mirror reflections.
• Validation: The authors note that their previous algorithms (for finite sources) have already been validated in astrophysics literature (e.g., DG1, DG2), where they outperformed competitors in accuracy and stability. This new formula extends that robustness to the infinite source case.

Conclusion

Fujimura and Vuorinen successfully extend the Ptolemy-Alhazen problem to the case of a source at infinity. By deriving a specific quartic equation and proving the location of its roots, they provide a mathematically rigorous and numerically stable method for calculating light reflection points on spherical mirrors. This work bridges geometric function theory with practical applications in signal processing and astrophysics, offering a superior alternative to iterative numerical methods for this specific geometric configuration.