Carl Størmer and his Numbers

Basic Information

• Paper ID: 2511.03030
• Title: Carl Størmer and his Numbers
• Authors: Matthew Kroesche, Lance L. Littlejohn, Graeme Reinhart
• Classification: math.HO (History of Mathematics), math.NT (Number Theory)
• Publication Date: November 4, 2025
• Paper Link: https://arxiv.org/abs/2511.03030

Abstract

In many proofs of Fermat's theorem on sums of two squares, the minimal positive residual solution $x_0$ of the quadratic congruence equation $x^2 \equiv -1 \bmod p$ plays a crucial role, where $p$ is a prime with $p \equiv 1 \bmod 4$. Such an $x_0$ is called a Størmer number, named after Norwegian mathematician and astronomer Carl Størmer (1874-1957). This paper establishes necessary and sufficient conditions for $x_0 \in \mathbb{N}$ to be a Størmer number for some prime $p \equiv 1 \bmod 4$. Størmer's investigation of these numbers originated primarily from his research on identities expressing $\pi$ as finite linear combinations of certain values of the Gregory-MacLaurin series $\arctan(1/x)$. Since fewer than 600 digits of $\pi$ were known by 1900, approximating $\pi$ was an important problem. An identity discovered by Størmer in 1896 was used by Kanada and his team in 2002 to compute 1.24 trillion digits of $\pi$.

Research Background and Motivation

Problem Background

1. Importance of Fermat's Two-Square Theorem: This theorem, called by G.H. Hardy "one of the most beautiful theorems in arithmetic," states that a prime $p$ can be expressed as the sum of two positive integer squares if and only if $p \equiv 1 \bmod 4$.
2. Central Role of Quadratic Congruences: Most proofs of Fermat's two-square theorem involve studying the quadratic congruence equation $x^2 \equiv -1 \bmod p$. For example, the proof independently given by Hermite and Serret in 1848 assumes the minimal residual solution $x_0$ is known, then uses the Euclidean algorithm to generate an algorithm for solving positive integers $a, b$ satisfying $p = a^2 + b^2$.
3. Historical Significance: Although Girard stated the theorem in 1632 and Fermat announced a proof in 1640, the first recorded existence proof was given by Euler in 1749, and the uniqueness proof was completed by Gauss in 1801.

Research Motivation

1. Theoretical Necessity: When the prime $p$ is large, finding $x_0$ is practically difficult, requiring systematic determination conditions.
2. Historical Value: Carl Størmer's work connects number theory and the approximation of $\pi$, possessing significant mathematical historical importance.
3. Applied Value: Størmer's identities continue to play a role in high-precision $\pi$ computation today.

Core Contributions

1. Established complete characterization of Størmer numbers: Provided necessary and sufficient conditions for a positive integer $x_0$ to be a Størmer number of some prime (Theorem 2.3)
2. Proved injectivity of the Størmer function: Demonstrated that each Størmer number corresponds to a unique prime (Theorem 2.1)
3. Provided heuristic proof of natural density of Størmer numbers: Supporting the Everest-Harman conjecture that the natural density of Størmer numbers is $\ln 2$
4. Systematically elucidated connections between Størmer numbers and Gregory numbers: Demonstrated Størmer's theory of $\pi$ approximation and its modern applications
5. Provided abundant computational examples: Illustrated practical applications of the theory through concrete examples

Methodology in Detail

Task Definition

Input: Positive integer $x_0$Output: Determine whether $x_0$ is a Størmer number for some prime $p \equiv 1 \bmod 4$; if so, find the corresponding prime $p$Constraints: $p$ must be a prime of the form $4n+1$, and $1 < x_0 \leq (p-1)/2$, $x_0^2 \equiv -1 \bmod p$

Core Theoretical Framework

Definition 2.1 (Størmer Number)

Let $p = 4n+1$ be a prime. A positive integer $x_0$ satisfying:

• (i) $1 < x_0 \leq \frac{p-1}{2}$
• (ii) $x_0^2 \equiv -1 \bmod p$

is called a Størmer number of $p$, denoted $S(p) = x_0$.

Key Theorems

Theorem 2.1 (Injectivity): The Størmer function $S: P_{4n+1} \to \mathbb{N}$ is injective.

Proof Sketch: Assume $S(p_1) = S(p_2)$ but $p_1 \neq p_2$. Without loss of generality, let $p_1 > p_2$. Since $p_j | (x_0^2 + 1)$, we have $p_1 p_2 | (x_0^2 + 1)$, thus: $x_0^2 + 1 \geq p_1 p_2$

On the other hand, since $x_0 \leq (p_2-1)/2$ and $p_2 \leq p_1 - 2$, we can derive: $x_0^2 + 1 < p_1 p_2$

This produces a contradiction, therefore $p_1 = p_2$.

Theorem 2.2 (Prime Properties): If $x_0^2 + 1 = 2^r p_1^{r_1} \cdots p_m^{r_m}$, then each $p_j \equiv 1 \bmod 4$.

Theorem 2.3 (Necessary and Sufficient Conditions): Let the prime factorization of $n = x_0^2 + 1$ be $2^r p_1^{r_1} \cdots p_m^{r_m}$, where $p_1 < \cdots < p_m$. Then there exists a prime $p$ such that $x_0$ is its Størmer number if and only if $2x_0 + 1 \leq p_m$. If the condition is satisfied, then $p = p_m$ and $S(p_m) = x_0$.

Technical Innovations

1. Unified determination framework: Transformed the problem of determining Størmer numbers into analyzing the prime factorization of $x_0^2 + 1$
2. Constructive proof: Not only provided determination conditions but also offered concrete methods for finding the corresponding prime
3. Connecting number theory and analysis: Established deep connections between Størmer numbers, Gregory series, and the approximation of $\pi$

Experimental Setup

Data Presentation

The paper provides two important data tables:

Table 1: Ordered pairs of primes and corresponding Størmer numbers

(5,2), (13,5), (17,4), (29,12), (37,6), (41,9), (53,23)
(61,11), (73,27), (89,14), (97,22), (101,10), (109,33), (113,15)
...

Table 2: The first several Størmer numbers in increasing order

1, 2, 4, 5, 6, 9, 10, 11, 12, 14, 15, 16, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 33, 34, 35, 36, 37, 39, 40, 42, ...

Verification Method

The paper verifies theoretical results through concrete examples:

Example 1: $x_0 = 3$

• $x_0^2 + 1 = 10 = 2 \times 5$
• $2x_0 + 1 = 7 > 5$
• Therefore, 3 is not a Størmer number

Example 2: $x_0 = 15$

• $x_0^2 + 1 = 226 = 2 \times 113$
• $2x_0 + 1 = 31 \leq 113$
• Therefore, $S(113) = 15$

Experimental Results

Natural Density Analysis

The paper provides numerical evidence for the natural density of Størmer numbers:

These data support the Everest-Harman conjecture: $\lim_{n \to \infty} \frac{|\{S(p) | p \leq n\}|}{n} = \ln 2 \approx 0.693$

Heuristic Density Proof

The paper provides a heuristic proof based on probability theory:

For $x_0 \in \mathbb{N}$, if $x_0 = S(p)$, then $2x_0 + 1 \leq p \leq x_0^2 + 1$. Assuming each integer has a "probability" of $2/(p-1)$ of being a Størmer number, the probability that $x_0$ is a Størmer number is approximately:

$\sum_{2x_0+1 \leq p \leq x_0^2+1, p \equiv 1 \bmod 4} \frac{2}{p-1} \sim \sum_{2x_0 \leq p \leq x_0^2} \frac{1}{p} \sim \ln 2$

Using Mertens' estimate yields the final result.

Gregory Numbers and π Approximation

Theoretical Foundation

Størmer established connections between Størmer numbers and Gregory numbers $t_x := \arctan(1/x)$:

Theorem 5.1 (Størmer's Theorem): If $x \in \mathbb{N}$ is not a Størmer number, then $t_x$ can be uniquely expressed as a finite linear combination of $t_n$, where each $n$ is a Størmer number with $n < x$.

Historical Applications

1. Machin's Formula (1706): $\frac{\pi}{4} = 4\arctan\frac{1}{5} - \arctan\frac{1}{239}$
2. Euler's Formula: $\frac{\pi}{4} = 5\arctan\frac{1}{7} + 2\arctan\frac{79}{3}$
3. Størmer's Formula (1896): $\frac{\pi}{4} = 44\arctan\frac{1}{57} + 7\arctan\frac{1}{239} - 12\arctan\frac{1}{682} + 24\arctan\frac{1}{12943}$

Modern Significance

Størmer's 1896 identity was used by Kanada's team in 2002 to compute 1.24 trillion digits of π, demonstrating the enduring value of classical mathematical theory.

Related Work

Historical Development

1. Fermat's Two-Square Theorem: Development from Girard (1632) through Fermat (1640) to Euler (1749)
2. Smith's Method (1855): Elegant existence proof using Euclidean algorithm and continued fraction theory
3. Hermite-Serret Algorithm (1848): Algorithm for constructing sum-of-squares decomposition based on known Størmer numbers

Related Research

1. Conway-Guy Definition: Slightly different definition of Størmer numbers, including the case n=1
2. Lehmer's Work: Research on identities of arccot functions
3. Todd's Theorem: Characterization of irreducibility of arctan functions

Conclusions and Discussion

Main Conclusions

1. Complete Characterization: Established necessary and sufficient determination conditions for Størmer numbers, resolving a long-standing problem
2. Theoretical Unification: Unified Størmer numbers, Fermat's two-square theorem, Gregory series, and π approximation within a single theoretical framework
3. Historical Value: Systematically reviewed Carl Størmer's mathematical contributions and their modern significance

Limitations

1. Density Proof: The proof of natural density remains heuristic, lacking mathematical rigor
2. Computational Complexity: For large $x_0$, determining whether it is a Størmer number requires factoring $x_0^2 + 1$, which may be computationally difficult
3. Application Scope: Primarily limited to theoretical mathematics with relatively limited practical applications

Future Directions

1. Rigorous Density Proof: Seeking rigorous mathematical proof of the Everest-Harman conjecture
2. Algorithm Optimization: Developing more efficient algorithms for determining and computing Størmer numbers
3. Generalization Research: Exploring applications of similar structures in other mathematical branches

In-Depth Evaluation

Strengths

1. Theoretical Completeness: Provides complete mathematical characterization of Størmer numbers, filling theoretical gaps
2. Historical Value: Systematically introduces Carl Størmer's life and contributions, possessing significant mathematical historical importance
3. Strong Connectivity: Skillfully connects number theory, analysis, and π computation across multiple mathematical branches
4. Rich Examples: Illustrates theoretical applications through abundant concrete examples, enhancing readability
5. Clear Writing: Well-structured paper with logical progression from basic definitions to advanced applications

Weaknesses

1. Computational Efficiency: For practical applications, the computational complexity of prime factorization may become a bottleneck
2. Density Proof: The proof of natural density lacks mathematical rigor, being merely heuristic
3. Limited Applications: While theoretically complete, modern practical value is relatively limited
4. Originality: Primarily systematizes existing theory; original contributions are relatively modest

Impact

1. Theoretical Contribution: Provides new tools and perspectives for number theory research
2. Educational Value: Suitable as teaching material for number theory and history of mathematics
3. Inspirational Significance: May inspire further research into related mathematical structures
4. Reproducibility: Theoretical results are clear and computational examples are easily verifiable

Applicable Scenarios

1. Number Theory Research: Provides new methods for studying quadratic congruences and prime properties
2. Mathematics Education: Serves as an excellent teaching case connecting different mathematical branches
3. History of Mathematics Research: Provides important material for understanding 19th-20th century mathematical development
4. Algorithm Design: Provides theoretical foundation for related computational problems

References

The paper cites 21 important references spanning from classical number theory textbooks to modern computational mathematics, reflecting the depth and breadth of the research. Key references include Andrews' Number Theory, Conway and Guy's The Book of Numbers, and Størmer's original papers.