Equidistribution of points on a sphere

  • Context: Graduate 
  • Thread starter Thread starter iironiic
  • Start date Start date
  • Tags Tags
    Points Sphere
Click For Summary
SUMMARY

The discussion focuses on the recursive construction of equidistributed points on a sphere using irrational numbers α and β. The points are defined by u_n = {nα} and v_n = {nβ}, with spherical coordinates θ_n = 2πu_n and ψ_n = arccos(2v_n - 1). The goal is to demonstrate that the limit of the ratio of points in a spherical cap S approaches h/2 as n approaches infinity. The use of Mathematica is mentioned for calculating the number of points k that satisfy the conditions for equidistribution.

PREREQUISITES
  • Understanding of irrational numbers and their properties
  • Familiarity with spherical coordinates and their mathematical representation
  • Knowledge of limits and asymptotic behavior in calculus
  • Experience with Mathematica for computational analysis
NEXT STEPS
  • Research the properties of equidistribution in number theory
  • Learn about the application of spherical coordinates in geometry
  • Explore the use of Mathematica for numerical simulations and calculations
  • Investigate the concept of simulated annealing for optimization problems
USEFUL FOR

Mathematicians, physicists, and computer scientists interested in geometric distributions, optimization algorithms, and numerical analysis techniques.

iironiic
Messages
9
Reaction score
0
I am struggling to find a way to count the number of irrational points defined recursively which satisfies specific conditions. This is the only hurdle I need to overcome, but the problem is proving itself extremely difficult.

Here is my construction of a recursion. My goal is to essentially prove that this recursion will guarantee a set of equidistributed points.

Let [itex]α[/itex] and [itex]β[/itex] be any positive irrational value. Define [itex]u_n = \{n α\} = n α \mod 1[/itex] and [itex]v_n = \{n β\} = n β \mod 1[/itex].

Consider the spherical coordinates, defined as follows:

[itex]θ_n = 2 π u_n[/itex]

[itex]ψ_n = \arccos{2 v_n - 1}[/itex]


This construction guarentees uniform distribution of points on the surface of the sphere. To show equidistribution, consider a spherical cap [itex]S[/itex] with height [itex]h[/itex] on the unit sphere (sphere with radius one). Given a particular [itex]n \in \mathbb{N}[/itex], my goal is to find the number of points that lie in [itex]S[/itex]. Inevitably, what I want to show is as follows:

[itex]\displaystyle \lim_{n \rightarrow \infty} \frac{\# p_n \in S}{n} = \frac{2 π h}{4π} = \frac{h}{2}[/itex].

This is equivalent to finding the number of [itex]n \in \mathbb{N}[/itex] (we'll call this number [itex]k[/itex]) satisfying,

[itex]\displaystyle u_n \in \left(0, \frac{\arccos{\frac{1-h}{2\sqrt{v_n(1-v_n)}}}}{2 \pi} \right) \cup \left(1- \frac{\arccos{\frac{1-h}{2\sqrt{v_n(1-v_n)}}}}{2 \pi}, 1 \right)[/itex]

[itex]\displaystyle v_n \in \left( \frac{1-\sqrt{1-(1-h)^2}}{2} , \frac{1+\sqrt{1-(1-h)^2}}{2}\right)[/itex]

For all [itex]0 \leq h \leq 1[/itex]. I used Mathematica to find [itex]k[/itex] and it seems that [itex]\frac{k}{n} \rightarrow \frac{h}{2}[/itex] where [itex]n[/itex] denotes the number of points on the surface of the entire sphere.

My question is, how do you use the definition of [itex]u_n[/itex] and [itex]v_n[/itex] to determine at least, the bounds for [itex]k[/itex]? Thank you!
 
Physics news on Phys.org
I'll follow this with interest. We had to equally distribute points on a sphere (to do with cell towers). The only way of doing it that we could find was by simulated annealing.
 
iironiic said:
[itex]ψ_n = \arccos{2 v_n - 1}[/itex]
You mean this, right?
[itex]ψ_n = \arccos(2 v_n - 1)/itex]<br /> <blockquote data-attributes="" data-quote="" data-source="" class="bbCodeBlock bbCodeBlock--expandable bbCodeBlock--quote js-expandWatch"> <div class="bbCodeBlock-content"> <div class="bbCodeBlock-expandContent js-expandContent "> Given a particular [itex]n \in \mathbb{N}[/itex], my goal is to find the number of points that lie in [itex]S[/itex]. </div> </div> </blockquote>I think you mean, for a given [itex]m \in \mathbb{N}[/itex] you want the number of points n <= m in S.<br /> Your expressions beyond that look much too complicated. Have you oriented the cap in the most convenient way? I get cos(2πu<sub>n</sub>) = cos(θ<sub>n</sub>) > 1 - h as the only constraint.[/itex]
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 21 ·
Replies
21
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
Replies
1
Views
1K
  • · Replies 16 ·
Replies
16
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K