Showing posts with label primes. Show all posts
Showing posts with label primes. Show all posts

Monday, 27 July 2026

Finding Numbers Within The Copeland-Erdos Constant

I watched this video recently about the Copeland-Erdos constant.

There following text was included below the video:

Write the primes in a row — 2, 3, 5, 7, 11, 13, and on forever — put a decimal point in front, and you get a single number: 0.235711131719232931… the Copeland–Erdős constant. In 1946, Copeland and Erdős proved it is normal in base ten: every finite string of digits appears in its expansion, and appears exactly as often as pure chance would predict — each digit one tenth of the time, each pair one hundredth, every block its fair share. So the opening of pi (314159) is in there, at digit 6,903. So is 42, at digit 214. So is every birthday, phone number, and message ever written.

What makes this remarkable: normality is the rule, not the exception — a randomly chosen real number is normal with probability one (Borel, 1909) — yet we still cannot prove that pi, e, or the square root of two are normal. The first explicit example was Champernowne's constant (0.123456789101112…), normal because it was built by pasting in every integer on purpose. The primes do no such favor: they thin out, follow no pattern, and were never designed to be tidy — and their concatenation is normal anyway. This video builds the number, defines normality, shows why the famous constants resist proof, and walks the intuition (via the prime number theorem) for why the patternless primes still even out in the limit.

I was interested in determining if and where a particular number, let's say my diurnal age today, might occur within this constant. Gemini created a program to do this, although I've written similar programs in the past myself for finding numbers within other constants like \( \pi\), \(e\) and \( \sqrt{2}\). Here is a permalink for the Gemini program. I've incorporated this program into my daily number analysis. Here is the output for the number associated with my diurnal age today: 28239.

The position in the Copeland-Erdos constant

Target '28239' first appears at position: 61014

What's of interest is whether a number can be equal to its index. I got Gemini to investigate and here is what it came up with:

While no fixed point where the 1-based index $P(N)$ equals $N$ itself has been found among smaller integers in the Copeland–Erdős constant, number-theoretic density principles dictate that such numbers must exist among larger composite numbers, even though they are mathematically impossible among prime numbers.

To see why this dichotomy exists, we have to look at the growth rate of the string indices for primes versus composite numbers.

Why Primes Never Equal Their Index

For any prime number $p$, its "natural" appearance occurs when $p$ itself is concatenated into the decimal stream. The 1-based index where $p$ begins is equal to $1$ plus the total number of digits of all primes strictly smaller than $p$:

$$P_{\text{natural}}(p) = 1 + \sum_{q < p} \lfloor \log_{10}(q) + 1 \rfloor$$

Using Chebyshev's function $\theta(p) = \sum_{q \le p} \ln(q)$, which by the Prime Number Theorem is asymptotically equivalent to $p$, we can evaluate the sum of the base-10 logarithms of all primes up to $p$:

$$\sum_{q < p} \log_{10}(q) = \frac{1}{\ln 10} \sum_{q < p} \ln(q) \sim \frac{p}{\ln 10} \approx 0.4343 \times p$$

Because $\frac{1}{\ln 10} \approx 0.4343 < 1$, the cumulative number of digits written into the stream grows at less than half the rate of the prime numbers themselves. By the time we concatenate a large prime $p$, its starting position is only around the $0.4343p$ mark. If $p$ happens to appear even earlier as an accidental substring across prime boundaries, its index is smaller still.

Therefore, for all primes, $P(p)$ is strictly less than $p$, making a fixed point $P(p) = p$ impossible.

Why Composites Can Form Fixed Points

Composite numbers behave completely differently because they never get a "natural" turn in the prime concatenation. They can only appear in two ways:

  1. As an accidental boundary string formed across the junction of two or more adjacent primes (for example, 4 first appears at index 21, spanning across the end of $3\mathbf{7}$ and the start of $\mathbf{4}1$).
  2. As a substring entirely within a much larger prime $q > N$.

Because the Copeland–Erdős constant is normal in base 10 (proven by Copeland and Erdős in 1946), every finite decimal string is guaranteed to appear eventually. However, the expected waiting time to see a specific $k$-digit composite string by chance in a normal sequence is roughly $10^k \ln 10 \approx 2.3 \times 10^k$.

Because $N$ itself sits between $10^{k-1}$ and $10^k$, a typical composite number that does not get an early boundary match will appear at an index much larger than $N$.

The Oscillation of $P(N) - N$

We can see this sharp split by examining the first few integers in the constant (OEIS A229190):

Integer $N$ Index $P(N)$ Type Behavior
2 1 Prime $P(N) < N$
3 2 Prime $P(N) < N$
4 21 Composite $P(N) > N$
5 3 Prime $P(N) < N$
6 31 Composite $P(N) > N$
7 4 Prime $P(N) < N$
8 41 Composite $P(N) > N$
9 12 Composite $P(N) > N$
10 47 Composite $P(N) > N$

Notice how consistently the primes lag behind their own values ($P(N) < N$), while the composite numbers overshoot ($P(N) > N$).

As we examine larger numbers, "lucky" composite numbers will occasionally appear unusually early due to boundary concatenations, while "unlucky" ones will appear late. Because the difference function $f(N) = P(N) - N$ oscillates wildly between negative values (for primes and early-matching composites) and large positive values (for late-matching composites), the graph of $P(N)$ repeatedly jumps across the line $Y = N$.

While finding a specific integer where $P(N) = N$ requires extensive string-searching through millions of digits, this continuous crossing confirms that fixed points exist among composite numbers.

Friday, 17 July 2026

Artin's Constant

My diurnal age today is 28229, a prime, with a reciprocal that has a period of 28228. It is thus one of those primes \(p\) with a period that is given by \( \dfrac{1}{p-1}\).

I got to wondering what proportion of primes have reciprocals that are one less than the prime itself. I asked Gemini. It's response follows and comprises the rest of this post.

Approximately 37.3956% of all prime numbers have reciprocals with a repeating decimal period of maximum length, \(p - 1\). In number theory, primes with this property are known as full reptend primes (or long primes) in base 10.

The Number Theory Behind Maximal Periods

The number of digits in the repeating period of the reciprocal \(1/p\) is determined by the multiplicative order of \(10\) modulo \(p\), denoted as \(\text{ord}_p(10)\).

By Fermat's Little Theorem, for any prime \(p \neq 2, 5\), the period length \(k\) is the smallest positive integer such that:

\[10^k \equiv 1 \pmod{p}\]

Since the multiplicative group \((\mathbb{Z}/p\mathbb{Z})^\times\) has order \(\phi(p) = p - 1\), Lagrange's theorem dictates that the period length \(k\) must always divide \(p - 1\).

For the period to achieve the absolute maximum length of \(p - 1\), the integer \(10\) must be a primitive root modulo \(p\). When \(10\) is a primitive root, the successive powers of \(10 \pmod{p}\) cycle through every possible non-zero remainder from \(1\) to \(p - 1\) before repeating, generating a cyclic period of length \(p - 1\).


Artin's Constant and Asymptotic Density

The exact asymptotic density of primes for which a given integer is a primitive root is governed by Artin's conjecture on primitive roots.

Because the base \(a = 10\) is neither a perfect square nor \(-1\), Artin's conjecture predicts that the proportion of full reptend primes among all primes up to \(N\) converges to a universal mathematical constant as \(N \to \infty\). This value is known as Artin's constant (\(C_{\text{Artin}}\)), defined by the infinite product over all prime numbers \(q\):

\[C_{\text{Artin}} = \prod_{q \text{ prime}} \left(1 - \frac{1}{q(q-1)}\right) \approx 0.3739558136\dots\]

While Emil Artin's conjecture remains unconditionally unproven for specific bases like 10, Christopher Hooley proved in 1967 that the asymptotic density is exactly \(C_{\text{Artin}}\) under the assumption of the Generalized Riemann Hypothesis (GRH). Extensive computational verification confirms that the distribution stabilizes at this ~37.3956% threshold (cataloged as decimal expansion A005596 in the OEIS).


Period Lengths of Small Primes

Comparing the first few primes highlights how the multiplicative order behaves when \(10\) is (and is not) a primitive root:

Prime (\(p\)) Reciprocal (\(1/p\)) Period Length (\(\text{ord}_p(10)\)) Full Reptend? (\(k = p - 1\))
3 \(0.\overline{3}\) \(1\) No
7 \(0.\overline{142857}\) \(6\) Yes
11 \(0.\overline{09}\) \(2\) No
13 \(0.\overline{076923}\) \(6\) No
17 \(0.\overline{0588235294117647}\) \(16\) Yes
19 \(0.\overline{052631578947368421}\) \(18\) Yes
23 \(0.\overline{0434782608695652173913}\) \(22\) Yes

The sequence of base-10 full reptend primes continues: \(7, 17, 19, 23, 29, 47, 59, 61, 97, 109, 113, 131, 149, \dots\) (OEIS sequence A001122). Noticeably, primes like \(13\) (where \(10^6 \equiv 1 \pmod{13}\)) and \(37\) (where \(10^3 \equiv 1 \pmod{37}\)) fall into sub-cyclic periods because their multiplicative orders are proper divisors of \(p - 1\).

Saturday, 27 June 2026

Prime Magic Squares

The number associated with my diurnal age today, \( \textbf{28209} \), is a member of OEIS A269324:


A269324
: magic sums of 3 x 3 magic squares composed of primes in an arithmetic progression.

I was interested in finding out what these primes were and how they could be arranged so I sought help from Gemini. Here was what it came up with:

To construct a 3×3 magic square of primes in an arithmetic progression for the magic sum 28209 (a term from OEIS A269324), we need to determine the central prime, find the common difference of the progression, and then map those primes to the cells of the grid. Here is the complete step-by-step breakdown of the mathematics and the final arrangement.

1. Finding the Central Prime

In any 3×3 magic square, the central cell is always the average of the magic sum. Because the magic sum (\( S \)) is the total of any given row, column, or diagonal, the center prime is exactly \( S / 3 \).

$$ \text{Center} = \frac{28209}{3} = 9403 $$

Because the 9 primes must form an arithmetic progression, the central cell of the magic square must also be the median (the 5th term) of the sorted sequence. Therefore, our arithmetic progression of primes is centered at 9403.

2. Determining the Common Difference

Let \( d \) be the common difference. Our 9 primes will take the form of \( 9403 + kd \), where \( k \) ranges from \( -4 \) to \( 4 \).

According to prime arithmetic progression properties, an arithmetic progression of \( k \) primes must have a common difference divisible by every prime \( q \le k \) (unless the progression actually begins with \( q \)). Since our sequence contains 9 primes and does not start with 2, 3, 5, or 7, the common difference \( d \) must be a multiple of their primorial:

$$ d = 2 \times 3 \times 5 \times 7 = 210 $$

By testing the multiples of 210, we can find the exact difference where all 9 resulting terms are prime:

  • \( d = 210 \) fails because the term \( 9403 - 210 = 9193 \), which is \( 29 \times 317 \).
  • \( d = 420 \) fails because \( 9403 + 420 = 9823 \), which is \( 11 \times 19 \times 47 \).
  • \( d = 630 \) fails because \( 9403 - 3(630) = 7513 \), which is \( 11 \times 683 \).
  • \( d = 840 \) succeeds.

3. The Nine Primes

Using \( d = 840 \) and our center prime of 9403, we can generate the 9 primes for the magic square by stepping out by 840 in both directions:

  • \( 9403 - 4(840) = \) 6043
  • \( 9403 - 3(840) = \) 6883
  • \( 9403 - 2(840) = \) 7723
  • \( 9403 - 1(840) = \) 8563
  • 9403 (Center)
  • \( 9403 + 1(840) = \) 10243
  • \( 9403 + 2(840) = \) 11083
  • \( 9403 + 3(840) = \) 11923
  • \( 9403 + 4(840) = \) 12763

4. Arranging the Magic Square

To arrange these into a valid 3×3 magic square, we map our ordered sequence to the classic Lo Shu magic square structure. In the traditional 1–9 Lo Shu square, the numbers are placed so that opposites balance out across the center (5). Substituting our 1st through 9th primes directly into the equivalent positions of the Lo Shu configuration yields the following layout:

11923 6043 10243
7723 9403 11083
8563 12763 6883

Every row, column, and main diagonal in this grid perfectly sums to 28209. Keep in mind that any of the 8 standard symmetries (rotations or reflections) of this exact grid will also produce a valid magic square for this progression.

Friday, 19 June 2026

Sierpinski Numbers

I watched this video on YouTube video about Sierpiński numbers. It's quite informative but I thought I'd find out more.


I asked Gemini to find out more. Here is its report:

The Sierpiński Problem: An Exhaustive Analysis of Covering Systems, Primality, and Distributed Computational Number Theory

The intersection of abstract Diophantine number theory and large-scale computational mathematics has produced some of the most rigorous and collaborative problem-solving efforts in modern scientific history. At the center of this nexus lies the Sierpiński problem, a profound question that challenges mathematicians to determine the smallest odd natural number $k$ for which the expression $k \cdot 2^n + 1$ never yields a prime number for any positive integer $n$. Such numbers, known formally as Sierpiński numbers, represent a striking anomaly in the asymptotic distribution of primes, guaranteeing an infinite sequence of composite integers.

This comprehensive report explores the mathematical properties of Sierpiński numbers, the life and theoretical contributions of their discoverer Wacław Sierpiński, the foundational proof mechanisms involving cyclotomic covering sets, and the ongoing, decades-long distributed computing efforts—most notably the "Seventeen or Bust" initiative—that seek to resolve the Sierpiński problem and its associated mathematical variants.

Historical Context: The Life and Legacy of Wacław Sierpiński

To fully grasp the significance of Sierpiński numbers, it is essential to examine the historical and intellectual environment of their namesake, Wacław Franciszek Sierpiński (1882–1969). An eminent Polish mathematician, Sierpiński made outstanding and foundational contributions to set theory, point set topology, the theory of functions of a real variable, and number theory [cite: 1, 2, 3].

Born in Warsaw during the Russian Empire's occupation of Poland, Sierpiński’s early education was severely constrained by the systematic suppression of Polish culture [cite: 4, 5]. Following sweeping educational mandates implemented between 1869 and 1874, the Russian authorities forced all secondary schools and universities in the territory to operate entirely in the Russian language [cite: 3, 6]. The implicit geopolitical aim was to suppress Polish intellectual advancement; however, this hostile environment inadvertently forged Sierpiński's resilience and fierce intellectual independence [cite: 5].

Despite these formidable barriers, Sierpiński enrolled in the Department of Mathematics and Physics at the Czar's University (the official, Russian-controlled iteration of the University of Warsaw) in 1899 [cite: 1, 7]. His prodigious talent in number theory was recognized almost immediately. In 1903, the department announced a competition for an essay on the number theory contributions of the distinguished Russian mathematician Georgy Voronoy [cite: 3]. Sierpiński submitted a dissertation that won the university's gold medal in 1904 [cite: 1, 7].

His prize-winning work significantly advanced the famous Gauss circle problem. Specifically, if $R(r)$ denotes the number of integer lattice points $(m, n)$ contained within a circle of radius $r$ centered at the origin, there exists a constant $C$ and a number $k$ such that the error term is bounded by $|R(r) - \pi r^2| < C r^k$ [cite: 1, 7]. Carl Friedrich Gauss had originally proved in 1837 that the minimal value $d$ for $k$ is $d \le 1$. Sierpiński's major contribution was proving that this inequality could be tightened to $d \le 2/3$ [cite: 1, 7]. In 1913, the prominent German mathematician Edmund Landau shortened Sierpiński's proof, publicly describing the underlying mathematics as exceptionally profound [cite: 1, 7].

Despite his academic triumph, Sierpiński's Polish nationalism led him to withdraw his gold-medal essay from the university's Russian-language Izvestia journal, choosing instead to publish it years later in a Polish mathematical magazine [cite: 7]. Following his graduation in 1904, his academic journey pivoted toward set theory and topology. He became deeply fascinated by the continuum hypothesis, the axiom of choice, and fractal curves, eventually introducing geometric concepts that bear his name, such as the Sierpiński triangle (or gasket), the Sierpiński carpet, and the Sierpiński space-filling curve [cite: 4]. He was also the first to provide a concrete example of an absolutely normal number in 1916—a number whose digits occur with equal asymptotic frequency in any given base [cite: 6].

During World War II, Sierpiński courageously continued his academic work within the "Underground Warsaw University" while officially employed as a municipal clerk, smuggling his mathematical papers out of occupied Poland to be published in Italy [cite: 1]. However, his enduring fascination with the structural properties of integers eventually led his focus back to pure number theory. In 1960, utilizing advanced techniques in modular congruences, Sierpiński published a seminal proof establishing the existence of the numbers that now bear his name [cite: 8, 9, 10].

Mathematical Foundations and Polignac's Conjecture

The expression $k \cdot 2^n + 1$, where $k$ is an odd positive integer and $n$ is a positive integer, generates a specific class of integers. When the exponentiated term overtakes the multiplier (specifically, when $2^n > k$), the resulting integer is known as a Proth number, named after the 19th-century French mathematician François Proth [cite: 1, 3]. Proth numbers are highly significant in computational number theory because their primality can be tested with exceptional algorithmic efficiency using Proth's theorem, circumventing the need for exhaustive trial division [cite: 1, 11].

A Sierpiński number is defined specifically as an odd natural number $k$ such that the Proth expression $k \cdot 2^n + 1$ is composite for every integer $n \ge 1$ [cite: 10, 12]. For the vast majority of $k$ values, evaluating successive values of $n$ quickly yields a prime number. For instance, if $k = 3$, setting $n = 2$ produces $3 \cdot 2^2 + 1 = 13$, which is prime; thus, 3 is immediately disqualified as a Sierpiński number [cite: 2, 3]. The core mathematical anomaly of a Sierpiński number is that its sequence completely evades the infinitely many prime numbers scattered throughout the integers.

The conceptual architecture of Sierpiński's 1960 proof relies heavily on addressing a flawed hypothesis known as Polignac's conjecture. In 1849, the French mathematician Alphonse de Polignac conjectured that every odd integer greater than 1 could be expressed as the sum of a prime number and a power of two (i.e., $p + 2^k$) [cite: 13, 14]. While Leonhard Euler had previously noted isolated counterexamples like 127 and 959, Polignac's conjecture remained a subject of intense debate until 1950 [cite: 14].

In 1950, the prolific Hungarian mathematician Paul Erdős and the Dutch mathematician J. G. van der Corput independently disproved Polignac's conjecture with rigorous finality [cite: 15]. Erdős proved that there exists an infinite arithmetic progression of odd numbers that cannot under any circumstances be represented as a sum of a prime and a power of two [cite: 9]. To achieve this, Erdős utilized a mathematical construct known as a "covering set" (or covering system) of congruences [cite: 13]. Erdős observed that by carefully selecting a finite set of prime numbers—specifically $\{3, 5, 7, 13, 17, 241\}$—one could create a modular interlocking system where every possible exponent $n$ aligns with a modulus that forces the resulting expression to be divisible by at least one of those selected primes [cite: 9, 16].

Sierpiński’s 1960 Proof and the Mechanics of Covering Sets

Building directly upon Erdős's foundational concept of covering systems, Sierpiński proved his eponymous theorem: there exist infinitely many odd positive integers $k$ such that $k \cdot 2^n + 1$ is strictly composite for all positive integers $n$ [cite: 17].

The proof mechanism relies on a finite covering set of congruences. A covering set $P = \{p_1, p_2, \dots, p_m\}$ is a finite set of prime numbers constructed so that for every positive integer $n$, the expression $k \cdot 2^n + 1$ is divisible by at least one $p_i \in P$ [cite: 18, 19]. Sierpiński constructed a precise covering of exponents $n$ using the moduli $\{2, 4, 8, 16, 32, 64\}$. This specific covering ensures that every integer $n$ satisfies at least one modular condition, triggering divisibility by a corresponding prime.

The mechanics of Sierpiński’s 1960 proof can be modeled by mapping the exponent conditions to their required prime divisors and the resulting congruences for $k$.

Index | Modulus for n | Exponent Congruence | Prime Divisor (p_i) | Required Congruence for k
1 | 2 | n = 1 mod 2 | 3 | k = 1 mod 3
2 | 4 | n = 2 mod 4 | 5 | k = 1 mod 5
3 | 8 | n = 4 mod 8 | 17 | k = 1 mod 17
4 | 16 | n = 8 mod 16 | 257 | k = 1 mod 257
5 | 32 | n = 16 mod 32 | 65537 | k = 1 mod 65537
6 | 64 | n = 32 mod 64 | 641 | k = 1 mod 641
7 | 64 | n = 0 mod 64 | 6700417 | k = -1 mod 6700417

Table 1: The covering system utilized in Wacław Sierpiński's 1960 proof [cite: 18].

The primes $\{3, 5, 17, 257, 641, 65537, 6700417\}$ form the covering set $C$ [cite: 1]. This specific set is highly delicate and demonstrates a deep structural connectivity to Fermat numbers. It functions explicitly because the fifth Fermat number, $F_5 = 2^{32} + 1$, factors cleanly into two distinct prime divisors: 641 and 6,700,417 [cite: 1].

By invoking the Chinese Remainder Theorem, Sierpiński proved that there exists a simultaneous solution $k$ to all the modular constraints listed in Table 1 [cite: 1]. Furthermore, because the prime moduli are pairwise coprime, the solutions for $k$ repeat periodically, forming an infinite arithmetic progression [cite: 1]. Any $k$ falling within this progression ensures that $k \cdot 2^n + 1$ is divisible by at least one prime in the covering set $C$ for any conceivable $n$, thus guaranteeing absolute compositeness [cite: 17, 19]. Solving this system of congruences yields the smallest $k$ native to Sierpiński's original method—an unwieldy 20-digit integer: 15,511,380,746,462,593,381 [cite: 17, 20].

John Selfridge and the Formulation of the Sierpiński Problem

While Sierpiński brilliantly proved the existence of these infinite progressions, he did not attempt to hunt for the absolute smallest example of such a number. In 1962, the American mathematician John Selfridge—renowned for his work on Fermat numbers, primality testing, and his proof alongside Erdős that the product of consecutive integers is never a power—addressed this gap [cite: 8, 21].

Selfridge discovered a far more elegant and numerically compact covering set [cite: 1, 22]. He demonstrated that $k = 78557$ is a Sierpiński number by employing a covering set consisting of only seven small primes: $\{3, 5, 7, 13, 19, 37, 73\}$ [cite: 3, 8]. For $k = 78557$, the covering functions beautifully over the modulus 36 [cite: 19]. Selfridge proved that every number of the form $78557 \cdot 2^n + 1$ is perfectly divisible by at least one prime in this covering set, thereby ensuring that no prime number can ever be generated by the sequence [cite: 8, 18, 21].

A critical nuance in the definition of a Sierpiński number is the strict requirement that $k$ must be an odd integer [cite: 8]. If this constraint were removed, the number $k = 65536$ ($2^{16}$) would likely be a smaller candidate [cite: 23]. For $k = 65536$, the expression becomes $2^{16} \cdot 2^n + 1 = 2^{n+16} + 1$ [cite: 24]. This sequence only yields primes if $n+16$ is a power of 2, generating Fermat numbers. Since it is widely believed (though unproven) that $2^{16} + 1$ is the largest prime Fermat number, $65536 \cdot 2^n + 1$ would technically be composite for all $n > 0$ [cite: 23, 25]. To avoid entangling the Sierpiński problem with the unproven finiteness of Fermat primes, the definition strictly excludes even numbers, solidifying 78,557 as the premier candidate [cite: 3, 24].

Following his discovery, Selfridge engaged in private correspondence with Paul Erdős, conjecturing that 78,557 was not merely a Sierpiński number, but the absolute smallest possible Sierpiński number [cite: 8, 26]. Determining the mathematical truth of this conjecture is known as the Sierpiński Problem [cite: 8, 10].

Algebraic Compositeness: Aurifeuillean Factorizations

While covering sets represent the dominant mathematical paradigm for proving that a number is a Sierpiński number, they are not the sole mechanism. In 1995, the mathematician A. S. Izotov demonstrated that certain fourth powers could be proven to be Sierpiński numbers without establishing a comprehensive covering set for all possible values of $n$ [cite: 8, 27].

Izotov's proof relied on an advanced algebraic identity known as Aurifeuillean factorization [cite: 1]. He showed that expressions of the form $t^4 \cdot 2^{4m+2} + 1$ can be factored purely algebraically, independent of modular prime divisibility:

$$t^4 \cdot 2^{4m+2} + 1 = (t^2 \cdot 2^{2m+1} + t \cdot 2^{m+1} + 1) \cdot (t^2 \cdot 2^{2m+1} - t \cdot 2^{m+1} + 1)$$ [cite: 1, 28].

Because this polynomial factorization always splits the expression into two distinct integers strictly greater than 1, it proves that any exponent of the form $n \equiv 2 \pmod 4$ natively gives rise to a composite number [cite: 8]. Therefore, a covering set is only required to eliminate the remaining exponent classes: $n \equiv 0, 1, \text{ and } 3 \pmod 4$ [cite: 8]. Izotov's technique provides a profound second-order insight: the compositeness of the $k \cdot 2^n + 1$ sequence can arise from inherent algebraic geometry and polynomial expansion, rather than relying solely on the cyclical, interlocking gears of prime modular arithmetic.

The Distributed Computation Era: Solving the Sierpiński Problem

To rigorously prove Selfridge's conjecture that 78,557 is indeed the smallest Sierpiński number, mathematicians face a daunting task of exhaustive elimination. They must prove that every single odd integer $k < 78557$ is not a Sierpiński number [cite: 3, 8]. A candidate $k$ is successfully eliminated if and only if a positive integer $n$ can be found such that $k \cdot 2^n + 1$ is prime [cite: 8].

Pre-Internet Computational Efforts

In the late 1970s and early 1980s, testing these numbers was an arduous process constrained by the hardware limitations of early academic mainframes. By 1983, computational searches had determined that there were 70 values of $k < 78557$ for which no prime had been found for exponents $n \le 8000$ [cite: 10, 16]. Over the subsequent 14 years, targeted algorithmic testing slowly eliminated 48 of those multipliers [cite: 10].

A major paradigm shift occurred in August 1997 with the introduction of the Proth.exe software developed by Yves Gallot [cite: 2]. This program heavily optimized Proth's theorem for standard consumer personal computers, allowing amateur mathematicians to participate in the search. Utilizing Proth.exe, a decentralized group of enthusiasts (including Lew Baxter, Marc Thibeault, and Janusz Szmidt) systematically eliminated more candidates [cite: 22]. By the end of 2001, the list of unsolved multipliers below 78,557 had been reduced to exactly 17 [cite: 1, 10].

The "Seventeen or Bust" Project (2002–2016)

Recognizing that the exponent $n$ required to find a prime for the remaining 17 values was reaching into the hundreds of thousands—necessitating vast, decentralized amounts of computational power—two college undergraduates, Louis Helm and David Norris, along with Michael Garrison, conceived the distributed computing project "Seventeen or Bust" in March 2002 [cite: 16, 29]. The project's public client was released on April 1, 2002, allowing volunteers worldwide to donate their idle CPU cycles to test specific $k$ and $n$ pairings [cite: 29, 30].

The project was immensely successful and heralded the golden era of crowdsourced mathematics. Within its first year, it eliminated five $k$ values. Over its 14-year lifespan as an independent project, Seventeen or Bust found massive primes that eliminated 11 of the 17 candidates [cite: 16]. The primes discovered were mathematically significant; for instance, the prime found for $k = 19249$ by Konstantin Agafonov in 2007 contained nearly 4 million digits, making it the largest known non-Mersenne prime at the time [cite: 2, 30].

In 2010, the global volunteer platform PrimeGrid officially partnered with Seventeen or Bust to accelerate the search using the BOINC (Berkeley Open Infrastructure for Network Computing) infrastructure [cite: 31, 32]. Unfortunately, in April 2016, the Seventeen or Bust project suffered a catastrophic datacenter failure resulting in the loss of critical un-backed-up data [cite: 10, 33]. Consequently, PrimeGrid assumed full administration of the project and absorbed its mission [cite: 16, 29, 33].

Shortly after this transition, PrimeGrid facilitated the discovery of a 12th prime, eliminating $k = 10223$. Discovered by Szabolcs Péter of Hungary on October 31, 2016, the prime $10223 \cdot 2^{31172165} + 1$ contains a staggering 9,383,761 decimal digits [cite: 1, 16, 33]. The calculation required nearly 9 days of continuous processing on an Intel i7 CPU [cite: 16, 33]. At the time of its discovery, it was the 7th largest known prime overall, the largest known Proth prime, and the largest known Colbert number [cite: 16, 33, 34].

k Value | Exponent (n) | Decimal Digits | Date Discovered | Discoverer
46157 | 698207 | 210186 | Nov 27, 2002 | Stephen Gibson
65567 | 1013803 | 305190 | Dec 2, 2002 | James Burt
44131 | 995972 | 299823 | Dec 5, 2002 | Anonymous
69109 | 1157446 | 348431 | Dec 6, 2002 | Sean DiMichele
54767 | 1337287 | 402569 | Dec 23, 2002 | Peter Coels
5359 | 5054502 | 1521561 | Dec 6, 2003 | Randy Sundquist
28433 | 7830457 | 2357207 | Dec 30, 2004 | Team Prime Rib
27653 | 9167433 | 2759677 | Jun 8, 2005 | Derek Gordon
4847 | 3321063 | 999744 | Oct 15, 2005 | Richard Hassler
19249 | 13018586 | 3918990 | Mar 26, 2007 | Konstantin Agafonov
33661 | 7031232 | 2116617 | Oct 17, 2007 | Sturle Sunde
10223 | 31172165 | 9383761 | Oct 31, 2016 | Szabolcs Péter

Table 2: The 12 mega-primes discovered by the Seventeen or Bust and PrimeGrid collaborations, successfully eliminating their respective $k$ values from the Sierpiński problem [cite: 31, 32].

Current Status: "Five or Bust"

As of the mid-2020s, decades of concerted computational effort have distilled the original Sierpiński problem down to exactly five remaining unproven multipliers [cite: 8, 35]. While the project name could accurately be updated to "Five or Bust," PrimeGrid retains the original moniker for historical continuity [cite: 29, 35].

To definitively prove that 78,557 is the smallest Sierpiński number, a prime must be found for each of the following $k$ values:

Remaining Candidates for the Sierpinski Problem
k = 21181
k = 22699
k = 24737
k = 55459
k = 67607

Table 3: The five remaining odd multipliers below 78,557 [cite: 8, 10, 35].

Computational Infrastructure and Hardware Dynamics

The search continues under PrimeGrid utilizing a global network of hundreds of thousands of volunteers, supplying thousands of TeraFLOPS of processing power to test exceedingly large exponents [cite: 31]. The search space has long surpassed $n = 36420000$, venturing into algorithmic territory where any prime discovered will be well over 10 million decimal digits long [cite: 10, 30].

The underlying mathematics of the search relies on Fast Fourier Transform (FFT) algorithms to multiply ultra-large integers efficiently. Applications like LLR (Lucas-Lehmer-Riesel), PRST, and Genefer are heavily optimized for both CPU and GPU execution [cite: 16, 35, 36]. The workflow is strictly bifurcated into two phases: sieving and primality testing. Sieving acts as a highly efficient filter, removing candidate exponents that have small algebraic factors. However, the deeper the sieve goes, the slower the rate of candidate removal becomes, eventually hitting an "optimal depth" where sieving takes just as much computational time as running a full primality test [cite: 28, 37].

Because computing large $n$ values requires intense, uninterrupted CPU time, multithreading is frequently employed. A single task running on one CPU core can take upwards of two to four weeks on older machines, making the tuning of parameters like max cpus critical for participants [cite: 35]. Interestingly, deep hardware analysis indicates that hyperthreading (or SMT) frequently decreases overall throughput for LLR primality testing. This negative performance scaling occurs because FFT calculations rely heavily on cache continuity and memory bandwidth rather than instruction pipeline multiplexing, prompting project organizers to advise volunteers to disable hyperthreading for LLR tasks [cite: 35, 38].

PrimeGrid also conducts rigorous "double checking" of legacy data. Because Seventeen or Bust suffered data loss, PrimeGrid systematically retests specific ranges of $n$ (comparing mathematical residues) to ensure no primes were missed due to hardware calculation errors or incomplete validation [cite: 30, 39, 40]. Two of the original twelve primes (for $k = 4847$ and $k = 33661$) were discovered exclusively because volunteers were double-checking previously computed ranges, highlighting the necessity of cryptographic-level verification in these searches [cite: 29, 39].

Generalizations and Extended Number Theory Problems

The theoretical proofs and computational methodologies developed for the Sierpiński problem have laid the groundwork for an entire taxonomy of associated mathematical conjectures. By altering the parameters, constraints, or signs of the $k \cdot 2^n + 1$ sequence, number theorists have generated a suite of derivative problems.

The Prime Sierpiński Problem

In 1976, Nathan Mendelsohn proved that $k = 271129$ is a Sierpiński number, utilizing the covering set $\{3, 5, 7, 13, 17, 241\}$ [cite: 1, 2]. Because 271,129 is itself a prime number, this discovery initiated the Prime Sierpiński Problem: determining the smallest Sierpiński number that is also a prime [cite: 8, 22].

To establish 271,129 as the smallest prime Sierpiński number, all prime integers less than 271,129 must be tested. Currently, there are nine prime multipliers for which no prime of the form $k \cdot 2^n + 1$ has been found: 22699, 67607, 79309, 79817, 152267, 156511, 222113, 225931, and 237019 [cite: 8, 22]. Notably, the first two values (22,699 and 67,607) are less than 78,557, meaning they are shared candidates with the original Seventeen or Bust project [cite: 10, 28]. PrimeGrid operates the "Prime Sierpinski Problem" (PSP) subproject specifically to eliminate these remaining values [cite: 3, 8, 28].

The Extended Sierpiński Problem

Assuming Selfridge's conjecture holds and 78,557 is the absolute smallest Sierpiński number, and assuming Mendelsohn's 271,129 is the smallest prime Sierpiński number, a structural gap exists. The Extended Sierpiński Problem asks: Is 271,129 the second smallest Sierpiński number overall? [cite: 8, 22].

Answering this requires searching the interval $78557 < k < 271129$ for any composite Sierpiński numbers [cite: 8, 22]. Eliminating the candidates in this range requires finding a prime for each composite $k$. As of the most recent computational updates, 11 composite multipliers remain in this interval: 21181, 24737, 55459, 91549, 131179, 163187, 200749, 209611, 227723, 229673, and 238411 [cite: 22]. (Again, 21181, 24737, and 55459 are shared with the original Five or Bust list, demonstrating the interwoven nature of these conjectures [cite: 22].)

Riesel Numbers and Dual Compositeness

The algebraic inverse of the Sierpiński function is the sequence $k \cdot 2^n - 1$. An odd natural number $k$ for which this expression is composite for all positive integers $n$ is known as a Riesel number [cite: 13, 15]. In 1956—four years prior to Sierpiński's proof—the Swedish mathematician Hans Riesel proved the existence of such numbers and demonstrated that $k = 509203$ is a Riesel number [cite: 15, 41].

Like Sierpiński numbers, Riesel numbers rely heavily on covering systems. However, the required primes often differ significantly due to the $-1$ altering the cyclotomic properties and algebraic factorization geometry [cite: 13, 19]. Unlike the $+1$ form, which allows flexible use of small primes like 5 and 7, Riesel coverings often require more complex congruences to avoid overlapping failures [cite: 8]. The quest to prove 509,203 is the absolute smallest Riesel number continues through PrimeGrid's "The Riesel Problem" subproject [cite: 6, 11].

Brier Numbers: Simultaneously Sierpiński and Riesel

A profound synthesis of these two sequences was theorized by Eric Brier, who asked whether an integer could be simultaneously a Sierpiński number and a Riesel number [cite: 8, 23, 42]. Such an integer, now known as a Brier number, must satisfy the condition that both $k \cdot 2^n + 1$ and $k \cdot 2^n - 1$ are composite for all positive integers $n$ [cite: 8, 17, 43].

Because a Brier number must satisfy the congruences of two independent covering sets simultaneously, the values of $k$ escalate dramatically. The smallest known Brier number, discovered by Christophe Clavier in 2013, is a 41-digit behemoth: 3,316,923,598,096,294,713,661 [cite: 1, 19, 24]. By leveraging a theorem of D. Shiu on the distribution of primes in arithmetic progressions, it has been mathematically proven that for every positive integer $r$, there exist $r$ consecutive primes that are also Brier numbers, hinting at the deeply patterned, yet elusive, nature of prime distributions [cite: 23, 42]. A 2021 Dartmouth College study by Pomerance et al. further proved that infinitely many numbers are simultaneously Sierpiński, Riesel, and Carmichael numbers, creating an intersection of pseudo-primality and forced compositeness [cite: 1].

Base Extensions, Repdigits, and the First Kind

While the standard definition operates in Base 2, the mathematics generalizes readily to any integer Base $b \ge 2$, where one searches for $k$ such that $k \cdot b^n + 1$ is composite for all $n$ [cite: 25, 44]. PrimeGrid operates projects like the Sierpinski/Riesel Base 5 problem to eliminate candidates where $b=5$ [cite: 6, 44]. In 2025, PrimeGrid successfully eliminated $k = 67612$ in Base 5 by discovering a 3.8-million-digit prime (which took only 1 hour and 24 minutes of PRP testing on an AMD Ryzen 9 7945HX), bringing the remaining Base 5 candidates down to 27 [cite: 6, 36].

Furthermore, Sierpiński numbers of the first kind refer to primes of the form $k^k + 1$. Sierpiński himself proved that for $k^k + 1$ to be prime for any $k>1$, $k$ must take the highly restrictive form $2^{2^j}$, tying the sequence directly back to Fermat numbers [cite: 45, 46]. Mathematical investigations have also extended to detecting Sierpiński and Riesel numbers among restrictive integer subsets, such as Fibonacci sequences, Lucas numbers, and repdigits (numbers consisting solely of a single repeated digit in a given base) [cite: 15, 47].

Broader Implications and Conjecture Outlook

The quest to resolve the Sierpiński problem acts as a crucible for both pure mathematical theory and applied computer science. At the theoretical level, the existence of Sierpiński numbers directly links to Paul Erdős's combinatorial number theory. Erdős conjectured that for any Sierpiński number $k$, the smallest prime divisor of $k \cdot 2^n + 1$ must be bounded as $n$ approaches infinity [cite: 15, 27]. If this is proven true, it solidifies the hypothesis that all Sierpiński numbers are strictly the product of covering systems, ruling out any undiscovered algebraic anomalies that might cause infinite compositeness without a finite, periodic covering set [cite: 13, 27].

At the applied level, the search drives massive innovations in high-performance computing (HPC). The algorithms used by PrimeGrid to multiply integers spanning millions of digits—relying on discrete Fourier transforms and advanced memory bandwidth management—routinely push the physical boundaries of hardware stability [cite: 31, 35]. The massive prime numbers discovered in these searches are frequently utilized to stress-test modern CPU and GPU architectures, successfully identifying microscopic, silicon-level calculation errors before those processors are deployed in commercial, scientific, or cryptographic environments.

Conclusion

The Sierpiński problem stands as a towering testament to the enduring complexity of prime number distribution. What began in 1960 as a brilliant theoretical proof utilizing interlocking covering sets of prime congruences has evolved into one of the largest, most sustained distributed computing efforts in human history.

John Selfridge’s 1962 discovery of $k = 78557$ established the absolute baseline, and despite decades of aggressive computational searching by thousands of global volunteers, his conjecture that it is the smallest Sierpiński number remains flawlessly intact. The field has systematically stripped away candidates, transforming the original "Seventeen or Bust" initiative into the current "Five or Bust" reality.

As algorithmic efficiency improves and global computational power scales, the mathematical community steadily marches toward a definitive resolution of the five remaining candidates: 21181, 22699, 24737, 55459, and 67607. Whether a prime is eventually found for each, finally proving Selfridge right, or whether one of these five numbers stubbornly remains composite for all $n$, thereby rewriting number theory textbooks, the journey itself has fundamentally advanced our understanding of prime dynamics, modular coverings, and the very limits of computational mathematics.

References

  • [1] https://en.wikipedia.org/wiki/Sierpi%C5%84ski_number
  • [2] http://www.primegrid.com/forum_thread.php?id=13544
  • [3] https://primegrid.fandom.com/wiki/Extended_Sierpinski_Problem
  • [4] https://www.primegrid.com/
  • [5] https://www.primegrid.com/forum_thread.php?id=972
  • [6] https://www.primegrid.com/forum_thread.php?id=1647
  • [7] http://www.prothsearch.com/sierp.html
  • [8] https://primegrid.fandom.com/wiki/Prime_Sierpinski_Problem
  • [10] https://t5k.org/glossary/page.php?sort=SierpinskiNumber
  • [11] https://primegrid.fandom.com/wiki/Seventeen_or_Bust
  • [12] https://arxiv.org/abs/2106.07376
  • [13] https://math.colgate.edu/~integers/aa12/aa12.pdf
  • [14] https://math.dartmouth.edu/~carlp/BFLPST_AMS121001.pdf
  • [15] https://www.fq.math.ca/Papers1/49-4/BaczkowskiFasorantiFinch.pdf
  • [16] https://www.primegrid.com/forum_thread.php?id=1647
  • [18] http://www.primegrid.com/forum_thread.php?id=13544
  • [19] https://www.primegrid.com/forum_thread.php?id=1750
  • [20] https://www.primegrid.com/forum_thread.php?id=972
  • [21] https://www.primegrid.com/forum_thread.php?id=1538
  • [22] https://www.primegrid.com/forum_thread.php?id=9614
  • [23] https://www.primegrid.com/forum_thread.php?id=8360
  • [24] https://www.primegrid.com/forum_thread.php?id=13913
  • [27] https://grokipedia.com/page/primegrid
  • [28] https://grokipedia.com/page/Safe_and_Sophie_Germain_primes
  • [30] https://www.primegrid.com/forum_thread.php?id=13913
  • [33] https://www.primegrid.com/forum_thread.php?id=9614
  • [35] https://cs.uwaterloo.ca/journals/JIS/VOL10/Jones/jones67.pdf
  • [37] https://math.dartmouth.edu/~carlp/BFLPST_AMS121001.pdf
  • [38] https://en.wikipedia.org/wiki/Sierpi%C5%84ski_number
  • [39] https://grokipedia.com/page/covering_set
  • [40] https://en.wikipedia.org/wiki/John_Selfridge
  • [41] https://encyclopedia.pub/entry/37744
  • [46] https://boincsynergy.ca/wiki/PrimeGrid
  • [47] https://cs.uwaterloo.ca/journals/JIS/VOL10/Jones/jones67.pdf
  • [48] https://en.wikipedia.org/wiki/Sierpi%C5%84ski_number
  • [49] https://primegrid.fandom.com/wiki/Seventeen_or_Bust
  • [50] http://www.prothsearch.com/sierp.html
  • [51] https://www.primegrid.com/forum_thread.php?id=1647
  • [53] https://www.primegrid.com/forum_thread.php?id=972
  • [54] https://www.primegrid.com/forum_thread.php?id=8360
  • [55] https://www.primegrid.com/forum_thread.php?id=7767
  • [61] https://grokipedia.com/page/primegrid
  • [62] https://www.primegrid.com/forum_thread.php?id=13913
  • [63] https://www.primegrid.com/forum_thread.php?id=9614
  • [64] http://www.primegrid.com/old_news.php
  • [66] https://www.primegrid.com/forum_thread.php?id=1538
  • [67] https://www.primegrid.com/forum_thread.php?id=9614
  • [71] https://www.primegrid.com/forum_thread.php?id=1538
  • [72] http://www.primegrid.com/forum_thread.php?id=1304&nowrap=true
  • [78] https://www.epcc.ed.ac.uk/whats-happening/articles/how-do-you-solve-problem-sierpinski
  • [79] http://www.primegrid.com/old_news.php
  • [81] https://www.primegrid.com/forum_thread.php?id=7116
  • [82] Unknown Source / Custom Execution
  • [84] https://mathshistory.st-andrews.ac.uk/Biographies/Sierpinski/
  • [86] https://ns1.almerja.com/more.php?idm=80076
  • [89] https://www.primepuzzles.net/problems/prob_049.htm
  • [92] https://poetcommons.whittier.edu/cgi/viewcontent.cgi?article=1009&context=math
  • [94] http://www.prothsearch.com/sierp.html
  • [96] https://encyclopedia.pub/entry/37744
  • [98] https://www.primegrid.com/forum_thread.php?id=9614
  • [103] http://www.primegrid.com/forum_thread.php?id=7356
  • [104] https://www.primegrid.com/forum_thread.php?id=8513&nowrap=true&menu=top
  • [108] https://www.primegrid.com/forum_thread.php?id=9614
  • [109] https://boincsynergy.ca/wiki/PrimeGrid
  • [111] https://www.primegrid.com/forum_thread.php?id=1647
  • [113] https://math.colgate.edu/~integers/aa12/aa12.pdf
  • [116] https://poetcommons.whittier.edu/cgi/viewcontent.cgi?article=1009&context=math
  • [117] https://people.math.sc.edu/filaseta/papers/SierpinskiEtCoPap5.pdf
  • [118] https://oeis.org/A121270
  • [119] https://oeis.org/A121270/internal
  • [125] https://scispace.com/pdf/the-mathematics-of-paul-erdos-i-25xah963v8.pdf
  • [126] https://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1639-08.pdf
  • [128] https://en.wikipedia.org/wiki/Riesel_number
  • [129] https://www.researchgate.net/publication/256775044_On_Cullen_numbers_which_are_both_Riesel_and_Sierpinski_numbers
  • [131] https://people.math.sc.edu/filaseta/ConsecutiveWDDBrierNumbers.pdf
  • [132] https://arxiv.org/pdf/2101.08898

Sunday, 31 May 2026

Palindromic Day 28182

Palindromic properties of 28182 (showing only sequence members up to 40000):


A098834: palindromic Smith numbers.

4, 22, 121, 202, 454, 535, 636, 666, 1111, 1881, 3663, 7227, 7447, 9229, 10201, 17271, 22522, 24142, 28182, 33633, 38283

A Smith number is a composite number where the sum of its digits equals the sum of the digits of its prime factors. For 28182:$$ \begin{align} 28182 &\rightarrow 2+8+1+8+2 = 21 \\ 28182 &= 2 \times 3 \times 7 \times 11 \times 61\\ &\rightarrow 2 + 3 + 7 + 1+1+6 +1 =21 \end{align}$$


A046395: palindromes that are the product of 5 distinct primes.

6006, 8778, 20202, 28182

Here \(28182 = 2 \times 3 \times 7 \times 11 \times 61 \)


A099052: all palindromes of length > 1 in the decimal expansion of \(e\).

\(e\) = 2.71828182845904523536028747135266 ...


A045571: numbers that are palindromic, divisible by 11 and have an odd number of digits.

121, 242, 363, 484, 616, 737, 858, 979, 10901, 11011, 12221, 13431, 14641, 15851, 17171, 18381, 19591, 20702, 21912, 22022, 23232, 24442, 25652, 26862, 28182, 29392, 30503, 31713, 32923, 33033, 34243, 35453, 36663, 37873, 39193

All the palindromic numbers with an even number of digits are divisible by 11. The number of palindromic numbers with \(2k+1\) digits that are divisible by 11 is \((10^{k+1} + (-1)^k)/11\), and their asymptotic relative density within the set of all palindromic numbers with an odd number of digits is 1/11 (from OEIS comments).


A113838
: palindromes sandwiched between twin primes.

4, 6, 282, 828, 858, 2112, 21012, 21612, 23832, 26262, 26862, 28182

Here of course the twin primes are 28181 and 28183.


A032751
: palindromic Super-3 Numbers.

4554, 6776, 17471, 22322, 22722, 28182

Super-3 numbers \(n\) are of the form \(3 \times n^3 \) and contain three consecutive 3's.

Here \(3 \times 28182^3 = 67148557\textbf{333}704\)

Friday, 29 May 2026

Wieferich Primes

This video on a special sort of prime piqued my curiosity. It's very well explained. Here is the video:


I asked Gemini to find more about these primes:

To understand Wieferich primes, we first need to look at Fermat’s Little Theorem, which states that for any prime $p$ and any integer $a$ not divisible by $p$:

$$a^{p-1} \equiv 1 \pmod{p}$$

A Wieferich prime takes this a step further. It is a prime number $p$ that satisfies this congruence not just modulo $p$, but modulo $p^2$.

The defining congruence for a generalized Wieferich prime (to base $a$) is:

$$a^{p-1} \equiv 1 \pmod{p^2}$$

Multiplying both sides by $a$, this is sometimes written as $a^p \equiv a \pmod{p^2}$. When mathematicians refer to a "Wieferich prime" without specifying a base, they are universally referring to base 2.

History and Fermat's Last Theorem

The concept was introduced in 1909 by German mathematician Arthur Wieferich. At the time, mathematicians were still fiercely trying to prove Fermat's Last Theorem (FLT), which states that $x^n + y^n = z^n$ has no non-zero integer solutions for $n > 2$.

Wieferich proved a monumental theorem regarding the "First Case" of FLT (where the exponent $p$ does not divide $x$, $y$, or $z$). He proved that if the First Case of Fermat's Last Theorem fails for a prime exponent $p$ (meaning a solution actually exists), then that prime $p$ must satisfy the condition $2^{p-1} \equiv 1 \pmod{p^2}$.

Because Wieferich primes are incredibly rare, this theorem effectively proved the First Case of FLT for almost all prime numbers, long before Andrew Wiles proved the entire theorem in 1994.

Current Status

Wieferich primes are astonishingly rare. Despite massive computational efforts, we only know of exactly two standard (base 2) Wieferich primes:

    • 1093 (discovered by Waldemar Meissner in 1913)
    • 3511 (discovered by Nicolaas Beeger in 1922)

Distributed computing projects, notably PrimeGrid, have systematically searched for new Wieferich primes. As of recent searches, it has been computationally verified that no other Wieferich primes exist below $6.7 \times 10^{15}$.

Mathematically, their theoretical status remains an open question. It is strongly conjectured that there are infinitely many Wieferich primes, but it has not even been rigorously proven that there is a finite or infinite number of them.

Generalized Wieferich Primes ($a \ge 2$)

We can generalize the definition to check for primes $p$ that satisfy $a^{p-1} \equiv 1 \pmod{p^2}$ for any integer base $a \ge 2$.

When tracking these through integer sequence catalogs like the OEIS, specific bases have their own dedicated entries. Here is a comprehensive list of the known generalized Wieferich primes for the first several bases:

Base ($a$) Known Wieferich Primes ($p$) Notes / OEIS Sequence
2 1093, 3511 The "standard" Wieferich primes. (A001220)
3 11, 1006003 Known as Mirimanoff primes. (A014127)
4 1093, 3511 Same as base 2, because $4 = 2^2$.
5 2, 20771, 40487, 53471161, 1645333507, 6692367337, 188748146801 (A014128)
6 66161, 534851, 3152573, 741584281375, 491531
8 1093, 3511 Same as base 2, because $8 = 2^3$.
9 11, 1006003 Same as base 3, because $9 = 3^2$.
10 3, 487, 56598313 (A014129)
11 71
12 2693, 11429, 81537269
13 2, 863, 1747591
14 29, 353, 5834244249152913116
16 1093, 3511 Same as base 2, because $16 = 2^4$.

If a prime $p$ is a Wieferich prime in base $a$, it will also be a Wieferich prime in base $a^k$ for any integer $k$. This is why bases 4, 8, and 16 perfectly mirror base 2, and base 9 mirrors base 3. The computational heavy lifting to find these values is immense, relying on optimized algorithms for modular exponentiation to scan through billions of candidate primes.

I was surprised not to find Arthur Wieferich listed in the MacTutor history of mathematicians but then I came across this interesting information about him on Wikipedia:

Arthur Josef Alwin Wieferich (April 27, 1884 – September 15, 1954) was a German mathematician and teacher, remembered for his work on number theory, as exemplified by a type of prime numbers named after him. 
 
He was born in Münster, attended the University of Münster (1903–1909) and then worked as a school teacher and tutor until his retirement in 1949. He married in 1916 and had no children. 
 
Wieferich abandoned his studies after his graduation and did not publish any paper after 1909. His mathematical reputation is founded on five papers he published while a student at Münster.

The first three papers are related to Waring's problem. His fourth paper led to the term Wieferich prime, which are p such that \(p^2\) divides \(2^(p-1) - 1\).

I got Gemini to research why he never published anything after 1909: 

The historical record does not contain a definitive letter, diary entry, or professional correspondence where Arthur Wieferich explicitly states why he walked away from research. Biographical information on him is extremely sparse, largely compiled from the archives of the University of Münster and state records in Hannover. However, piecing together his timeline and the structure of the German academic system at the time provides a clear picture of why a mathematician of his caliber would vanish from the literature.

The German Academic Bottleneck

In the early 20th century, publishing groundbreaking papers—as Wieferich did between \( 1908 \) and \( 1909 \)—was not a guaranteed ticket to a university professorship. To secure a permanent academic post in Germany, a scholar had to complete a Habilitation and then work for years as a Privatdozent. A Privatdozent was an external lecturer who received no salary from the university, relying solely on the meager fees paid by the specific students attending their lectures.

Unless a young mathematician had independent wealth or substantial family backing, this path was financially ruinous. Wieferich was the son of a businessman, but evidently chose the alternative: the highly respected, financially secure civil service path of a Gymnasiallehrer (secondary school teacher).

The Trajectory of a Teacher

Wieferich's mathematical burst occurred entirely while he was a student, likely heavily influenced by a number theory course taught by Max Dehn in \( 1907 \). His famous paper on Fermat's Last Theorem, which defined what are now called Wieferich primes—primes \( p \) such that \( p^2 \mid 2^{p-1} - 1 \)—was published when he was just \( 25 \).

Upon graduating in \( 1909 \), he passed the state examination and immediately entered the school system. He was not a researcher who slowly faded away due to lack of ideas; he simply changed careers upon graduation. He taught in a succession of towns: Konitz, Elbing, Zoppot, Neustadt, Jülich, Stade, and finally Meppen.

He did not sever ties with the mathematical community immediately. He remained a member of the Deutsche Mathematiker-Vereinigung (the German Mathematical Society) until about \( 1929 \), twenty years after his last publication. This suggests he maintained a personal interest in the field, even if his profession left no time or incentive for active research.

A Minor Footnote: The Kempner Correction

There is one mathematical footnote regarding his departure. His first paper in \( 1908 \) tackled Waring's problem, attempting to prove that every integer can be represented as the sum of at most \( 9 \) positive cubes. The formidable Edmund Landau publicly praised the paper as a massive advance in elementary number theory.

However, Wieferich's proof contained a logical gap. It was not fatal to the overall idea, but it required an external correction, which was published by Aubrey J. Kempner in \( 1912 \). There is no historical evidence suggesting this correction discouraged Wieferich—he was already three years into his teaching career by the time Kempner published the fix—but it is the only major post-publication event related to his work during his lifetime.

Following World War II, from \( 1945 \) until his retirement in \( 1949 \), Wieferich supported himself as a private tutor. He died in \( 1954 \) in Meppen, leaving behind no children and only five papers—yet enough of a legacy to ensure his name remains permanently attached to number theory.

Thursday, 28 May 2026

Full Report on Brun's Constants

 Here is Gemini's full report on Brun's Constants:

Comprehensive Analysis of Brun's Constants: Twin Primes (B2) and Prime Quadruplets (B4)

Introduction to Prime Constellations and Asymptotic Densities

Within the discipline of analytic number theory, the distribution of prime numbers across the integer continuum represents one of the most profound and heavily investigated landscapes for mathematical inquiry. While the Prime Number Theorem provides a robust asymptotic framework for the frequency of primes generally, dictating that the number of primes less than or equal to a given real number $x$ is approximately $x / \ln(x)$, the behavior of prime constellations remains an area of intense theoretical and computational research. Prime constellations are specific groupings of prime numbers separated by fixed, predetermined intervals. The simplest and most historically scrutinized of these constellations is the twin prime pair, consisting of two primes separated by a distance of two, taking the form $(p, p+2)$. Extending this logic to denser topologies yields higher-order constellations, such as prime quadruplets, which represent the closest permissible clustering of four distinct prime numbers.

Among the most significant mathematical discoveries regarding these constellations is Viggo Brun's 1919 theorem. Utilizing revolutionary combinatorial techniques that would later become known as Brun's sieve, he successfully established that the sum of the reciprocals of all twin primes converges to a distinct, finite limit. This limit is universally recognized today as Brun's Constant for twin primes, denoted as $B_2$. The convergence of this series fundamentally differentiated the set of twin primes from the set of all prime numbers. It is a foundational proof in mathematics, dating back to Leonhard Euler and later formalized by Mertens, that the sum of the reciprocals of all prime numbers diverges to infinity. By proving that the harmonic series of twin primes converges, Brun demonstrated that twin primes are exceptionally sparse compared to primes in general, laying the groundwork for modern sieve theory.

Beyond the twin prime pairs, the architecture of prime constellations extends to prime quadruplets. The sum of the reciprocals of these specific prime quadruplets similarly converges to a distinct finite value, identified as Brun's Constant for prime quadruplets, or $B_4$. Calculating the precise numerical values of these twin and quadruplet constants has historically functioned as a premier benchmark for computational mathematics. The extraction of $B_2$ and $B_4$ relies not merely on brute-force enumeration, which is mathematically inadequate due to the agonizingly slow convergence of the series, but on sophisticated asymptotic extrapolation models utilizing the Hardy-Littlewood conjectures, advanced prime sieving algorithms, and conditional bounds relying on the Generalized Riemann Hypothesis (GRH). This report provides an exhaustive, peer-level analysis of Brun's Constants, investigating the theoretical foundations of their convergence, the historical trajectory of their numerical estimation, modern algorithmic paradigms for tabulating prime patterns up to immense thresholds, and the rigorous unconditional and conditional upper bounds established by contemporary mathematical literature.

Theoretical Foundations of Viggo Brun's Theorem

The Density of Twin Primes and Sieve Theory

A twin prime is strictly defined as a pair of prime numbers $(p, p+2)$. To evaluate the distribution and overall density of these pairs, mathematicians utilize the twin prime counting function, denoted mathematically as $\pi_2(x)$, which quantifies the absolute number of primes $p \leq x$ for which $p+2$ is also a prime number. The convergence of the sum of the reciprocals of twin primes is a direct mathematical consequence of rigorous upper bounds established on the density of the twin prime sequence.

Brun's pioneering application of pure sieve methods allowed him to demonstrate a strict upper bound on this counting function. He proved that the number of twin primes less than $x$ obeys the following asymptotic limitation:

$$\pi_2(x) = O\left( \frac{x (\log \log x)^2}{(\log x)^2} \right)$$

This explicit bound is profoundly significant because it indicates that twin primes occur far less frequently than general primes, being separated by nearly an additional logarithmic factor. For general primes, the density is governed by $x / \log x$, whereas the twin prime density is suppressed by an additional divisor of $\log x$, alongside the highly dampening double-logarithmic numerator.

From this strict density bound, the mathematical intuition for convergence naturally emerges. By applying the techniques of partial summation to the density bound, the infinite sum of the reciprocals of the twin primes can be evaluated. The sum for Brun's Constant $B_2$ is defined as the addition of both elements of the twin prime pair:

$$B_2 = \left(\frac{1}{3} + \frac{1}{5}\right) + \left(\frac{1}{5} + \frac{1}{7}\right) + \left(\frac{1}{11} + \frac{1}{13}\right) + \left(\frac{1}{17} + \frac{1}{19}\right) + \dots$$

In explicit mathematical terms, this infinite series must either possess finitely many terms, or it must possess infinitely many terms that undergo strict and absolute convergence to a specific value. The determination of this value, $B_2$, defines the scope of twin prime research.

The Informational Paradox of Convergence

Brun's theorem holds immense historical weight within analytic number theory, yet it generates a profound informational paradox regarding the unresolved Twin Prime Conjecture. The Twin Prime Conjecture postulates that there are infinitely many twin primes spanning the integer continuum. If the sum of the reciprocals of twin primes diverged to infinity—exactly as it does for the set of all primes—that mathematical divergence would serve as undeniable, conclusive proof that there are infinitely many twin prime pairs.

Because the sum strictly converges to a finite value, the theorem remains entirely agnostic on the finiteness of the set itself. It is mathematically impossible to conclude from Brun's theorem alone whether the set of twin primes is finite or infinite; the convergence merely establishes that the twin primes form a relatively "small set" within the broader integer continuum. However, the arithmetic nature of Brun's Constant provides a conditional diagnostic tool that mathematicians continue to probe. If $B_2$ were ever formally proven to be an irrational number, it would definitively imply the existence of infinitely many twin primes. This is because any finite sum of rational numbers (which the reciprocals of primes are) must necessarily result in a rational number. Therefore, irrationality would guarantee an infinite series of terms. Conversely, if $B_2$ is ever proven to be a rational number, the Twin Prime Conjecture would remain an unsolved mystery, as infinite sums of rational numbers can readily converge to rational limits. Because the rationality of $B_2$ remains unknown, the evaluation of the constant relies heavily on expanding the numerical precision of its calculation.

The Hardy-Littlewood Conjecture and Asymptotic Extrapolation

The Computational Impossibility of Direct Summation

The summation of twin prime reciprocals poses an extreme computational hurdle for computer scientists and mathematicians: the series converges at an agonizingly slow and irregular pace. Directly summing the reciprocals of twin primes up to a specific prime threshold $p$, a value denoted as $B_2(p)$, fails to provide an accurate estimate of the true infinite sum $B_2$ due to massive truncation errors. Direct calculations reveal that $B_2(p)$ grows with high volatility in the lower integer registers and flattens out into imperceptible microscopic increments as $p$ scales.

To put the agonizing rate of this convergence into perspective, number theorists have calculated that if one were to rely strictly on direct, brute-force summation without predictive mathematical modeling, the physical value of the sum $B_2(p)$ would not even reach the integer threshold of 1.9 until $p$ approximated $10^{530}$. A computation of this colossal magnitude surpasses the theoretical limits of any physical computing architecture, demanding an alternative mathematical approach to circumvent the limitations of brute force.

Applying Conjecture B to Prime Pairs

To successfully bypass this computational ceiling, number theorists utilize asymptotic extrapolation derived directly from the Hardy-Littlewood Twin Prime Conjecture, widely referenced in mathematical literature as Conjecture B. Proposed in 1922, Conjecture B provides a highly accurate probabilistic approximation for the twin prime counting function $\pi_2(n)$ for substantially large values of $n$. The formulation relies on the logarithmic integral $Li_2(n)$ and is mathematically expressed as:

$$\pi_2(n) \sim 2C_2 Li_2(n) = 2C_2 \int_{2}^{n} \frac{dt}{\log^2(t)} \quad \text{or} \quad \pi_2(n) \sim 2C_2 \frac{n}{\log^2(n)}$$

In this formulation, $C_2$ represents the foundational twin prime constant. The necessity of $C_2$ arises from the fact that the probability of a number $n$ being prime is roughly $1/\ln(n)$, but the events of $n$ and $n+2$ being prime are not statistically independent due to shared divisibility traits (e.g., if $n$ is odd, $n+2$ is guaranteed to be odd). $C_2$ corrects for this arithmetic dependence and is defined as an infinite product spanning all primes $p \geq 3$:

$$C_2 = \prod_{p \geq 3} \left(1 - \frac{1}{(p-1)^2}\right) \approx 0.6601618158...$$

An elegant analytical consequence of this conjecture allows mathematicians to algebraically re-express the partial direct sums of Brun's constant, $B_2(p)$, incorporating an asymptotic error term. Using the Hardy-Littlewood framework, the relationship is strictly defined as:

$$B_2(p) = B_2 - \frac{4C_2}{\log(p)} + O\left(\frac{1}{\sqrt{p}\log(p)}\right)$$

The Extrapolated Sum Function $B^*_2(p)$

By isolating the true limit $B_2$ from the asymptotic expansion, mathematicians can define a newly extrapolated function, denoted as $B^*_2(p)$, which converges to the true limit of Brun's Constant at an exponentially faster rate than the raw direct sum. This extrapolated value is calculated dynamically as:

$$B^*_2(p) = B_2(p) + \frac{4C_2}{\log(p)}$$

Plotting the direct sum $B_2(p)$ as a function of the inverse logarithmic parameter $1/\log(p)$ generates a strictly linear relationship. The resulting line possesses a negative slope that mathematically approximates $-4C_2$, which evaluates numerically to -2.64064726. The theoretical point at which this extrapolated line intersects the vertical y-axis represents the infinite limit $p = \infty$, yielding the true theoretical value of Brun's Constant assuming the validity of the Hardy-Littlewood heuristics.

The drastic and immediate improvement provided by the Hardy-Littlewood extrapolation method is highly evident when comparing the direct summation against the extrapolated sums across ascending magnitudes of the prime threshold $p$. The data below illustrates this profound divergence in convergence efficiency:

Threshold $p$ Direct Sum $B_2(p)$ Extrapolated Sum $B^*_2(p)$
$10^2$ $1.330990365719...$ $1.904399633290...$
$10^4$ $1.616893557432...$ $1.903598191217...$
$10^6$ $1.710776930804...$ $1.901913353327...$
$10^8$ $1.758815621067...$ $1.902167937960...$
$10^{10}$ $1.787478502719...$ $1.902160356233...$
$10^{12}$ $1.806592419175...$ $1.902160630437...$
$10^{14}$ $1.820244968130...$ $1.902160577783...$
$10^{15}$ $1.825706013240...$ $1.902160582249...$
$10^{16}$ $1.830484424658...$ $1.902160583104...$

As demonstrated by the structured numerical data, at the extreme computational threshold of $p = 10^{16}$, the direct sum has merely reached 1.83048, representing an unacceptable relative truncation error. Conversely, the extrapolated sum has stabilized with extreme precision, providing the most robust modern estimates of the true value of $B_2$.

Chronology of Computational Estimation for $B_2$

The pursuit of an increasingly accurate numerical value for Brun's Constant has driven some of the most intensive computational mathematics efforts of the last century. Because estimating $B_2$ relies on summing billions of minuscule floating-point fractions, these efforts continuously operate at the absolute limits of contemporary microprocessor architecture.

Pre-Digital Enumerations to Early Supercomputing

The foundational efforts to estimate twin prime densities occurred entirely in the pre-computer era, relying on exhaustive manual arithmetic. In 1878, the mathematician J. W. L. Glaisher successfully enumerated $\pi_2(10^5)$ utilizing hand-calculated integer tables to study the localized distribution of prime pairs. Decades after Brun's 1919 proof of strict convergence, the advent of digital supercomputing allowed for massive expansions in algorithmic capability.

In 1974, Shanks and Wrench pushed the envelope by computing the twin primes up to 2 million, evaluating the constant's early trajectory. In 1976, Richard P. Brent utilized mainframe architecture to push this boundary to an unprecedented $1 \times 10^{11}$. His evaluation of irregularities in twin prime distributions produced a highly accurate heuristic estimate of $B_2 \approx 1.902160540$.

Thomas Nicely and the Pentium FDIV Bug

The most historically famous and consequential computation of Brun's Constant occurred in the mid-1990s through the independent work of Dr. Thomas R. Nicely. Attempting to generate a high-precision estimate, Nicely meticulously expanded the twin prime enumeration up to $10^{14}$. Because the computation of Brun's constant requires the consecutive addition of incredibly small fractional reciprocals floating at the extreme edge of machine precision, it acts as a severe mathematical stress-test for floating-point arithmetic hardware. Nicely noted that even after summing the first billion ($10^9$) terms, the relative error of the raw sum was still greater than 5%, necessitating absolute hardware reliability for the ensuing trillions of operations.

In 1994, during these exhausting calculations, Nicely noticed persistent inconsistencies in his extrapolated sums that defied mathematical logic. Subsequent rigorous debugging revealed that the discrepancy was not algorithmic, but rather physically embedded in the silicon: Nicely had discovered a critical flaw in the Floating-Point Unit (FPU) of the newly released Intel Pentium microprocessor. The hardware error, famously immortalized as the Pentium FDIV bug, caused the processor to return erroneous binary floating-point results when dividing certain numbers. Intel subsequently attributed the error to missing entries in the microcode lookup table utilized by the FPU for accelerated division algorithms. Specific hardware calculations, such as 0001 / 824633702441.0, failed silently, generating errors beyond the eighth significant digit which compounded disastrously in Nicely's massive summation loops.

Nicely's mathematical discovery forced a massive, multi-million-dollar hardware recall, costing Intel an estimated 475 million to 1 billion in write-offs, illustrating a unique historical intersection where pure analytic number theory served as an inadvertent diagnostic tool for industrial microprocessor architecture. Accounting for and correcting the hardware limitations, Nicely heuristically estimated Brun's constant at $B_2 \approx 1.902160578$ for the $10^{14}$ threshold.

The $10^{16}$ Threshold and Modern Refinements

The early 2000s marked the transition toward distributed, memory-optimized computations. In 2002, mathematicians Pascal Sebah and Patrick Demichel conducted an exhaustive evaluation of all twin primes up to $10^{16}$. Sebah verified Nicely's earlier results independently and expanded the raw sum of twin prime reciprocals strictly below $10^{16}$ to exactly $1.830484424658...$. Applying the Hardy-Littlewood extrapolation algorithm to this direct sum, Sebah and Demichel derived the current, globally accepted consensus estimate for Brun's Constant:

$$B_2 \approx 1.902160583104$$

Nicely himself continued to stretch these boundaries over the following decade, subsequently expanding his own prime pair computations to $1.6 \times 10^{15}$ by January 2010. These subsequent runs continually reinforced the mathematical stability of the 1.902160583 estimation, confirming that no anomalous density fluctuations existed in the higher prime registers.

The cultural footprint of this specific numerical constant even extended into corporate finance. When Google placed bids during the Nortel patent auction in 2011, they utilized famous mathematical constants to signal their strategic bidding architecture. One of their exact bids was precisely $1,902,160,540, mirroring the first ten digits of Brun's constant as derived by Richard Brent in 1976.

The Twin Prime Counting Function $\pi_2(n)$ at Extreme Scale

To achieve the precise evaluation of $B_2$, the absolute count of twin primes up to massive integers must be flawlessly tabulated. The behavior of the twin prime counting function $\pi_2(n)$ verifies the density limitations predicted by Hardy-Littlewood. Analytical tabulations show that as $n$ increases by orders of magnitude, the absolute count of twin primes grows, yet their relative density to the overall integer baseline continuously rarefies.

The absolute evaluation of twin primes up to $10^{18}$ showcases the computational magnitude required for these modern estimates:

Threshold $n$ Twin Prime Count $\pi_2(n)$
$10^{11}$ $224,376,048$
$10^{12}$ $1,870,585,220$
$10^{13}$ $15,834,664,872$
$10^{14}$ $135,780,321,665$
$10^{15}$ $1,177,209,242,304$
$10^{16}$ $10,304,195,697,298$
$10^{17}$ $90,948,839,353,159$
$10^{18}$ $808,675,888,577,436$

As the threshold escalates to $10^{18}$, observing over 808 trillion twin prime pairs validates the structural integrity of the $2C_2 n / (\ln \ln(n))^2$ predictive bounds, ensuring that no sudden divergences skew the calculation of the reciprocal sums.

Rigorous Analytical Bounds on Brun's Constant $B_2$

While the extrapolated heuristic value of 1.902160583104 is universally accepted as highly accurate by the computational mathematics community, proving rigorous unconditional bounds for the infinite sum remains an exceptionally complex analytical challenge. Direct enumeration only provides a floor for the sum; determining the strict mathematical ceiling requires bounding the error terms of infinite series. Significant advancements have been made in both unconditional mathematical frameworks and in bounds conditionally dependent on the Generalized Riemann Hypothesis (GRH).

Unconditional Upper Bounds

The first modern rigorous unconditional upper bound on Brun's Constant was established by Richard Crandall and Carl Pomerance. Relying strictly on proven density theorems and avoiding any unproven conjectures, Crandall and Pomerance definitively proved that $B_2 < 2.347$.

This bound remained the gold standard until a significant mathematical breakthrough was achieved in 2018. Mathematicians Dave Platt and Tim Trudgian published a comprehensive analysis that improved the unconditional bounds by approximately 13% over the prior published results. Utilizing Riesel and Vaughan's bounds for the counting function $\pi_2(x)$, Platt and Trudgian integrated improvements on older sieving inequalities originally developed by Montgomery and Vaughan. Furthermore, they utilized explicit analytical bounds on the sums of divisors to control the error variables at infinite scales.

Their integration techniques heavily utilized bounding equations of the form $\int t^{6/5} \log t$ and incorporated the Exponential Integral $Ei(x)$ to manage the convergence decay. Through this rigorous arithmetic interval analysis, Platt and Trudgian successfully established a tightened unconditional interval:

$$1.840503 < B_2 < 2.288513$$

(An intermediate optimized calculation within their framework also references a slightly tighter bound of $B_2 < 2.288490$). This tightened mathematical corridor limits the maximum theoretical ceiling of $B_2$ to 2.288513 without requiring the assumption of complex analytical conjectures, marking a triumph of unconditional sieve theory.

Bounds Conditional on the Generalized Riemann Hypothesis (GRH)

Considerably sharper bounds can be achieved by assuming the truth of the Generalized Riemann Hypothesis (GRH). GRH postulates that the non-trivial zeros of all Dirichlet L-functions lie precisely on the critical line where the real part equals 1/2. Because GRH allows for deeply restricted error terms regarding the distribution of primes in specific arithmetic progressions, it provides highly constrained, noise-free density bounds for twin primes.

In 2007, Dominic Klyve produced an unpublished doctoral thesis at Dartmouth College utilizing GRH to derive a substantially tighter upper bound. Klyve relied heavily on advanced numerical integration methods over constrained intervals to compute the subsequent bounds, ultimately showing conditionally that if GRH holds, $B_2 < 2.1754$.

The bounds were pushed further in 2025 by Lachlan Dunn, who published an analysis further constraining this conditional upper limit. Dunn achieved a rigorously proven upper bound of $B_2 < 2.1594$ under the explicit assumption of GRH. Dunn achieved this by developing a specialized sieve framework based on an intermediate sum architecture. Rather than attempting to bound the entire infinite sum directly as a monolithic entity, Dunn decomposed the mathematical evaluation into distinct, manageable regions.

These regions were represented as intervals $B(m_i, m_{i+1})$, where $m_1 = x_0$ up to an arbitrary boundary $k$. This intermediate sum approach provided critical analytical flexibility. Because different mathematical estimation methods yield tighter bounds over varying intervals of $x$, Dunn could apply highly specialized local estimations to each regional block. For instance, setting boundaries at $L_1 = 10^8$ and $L_2 = 10^{10}$, Dunn computed values exactly below $L_1$, utilized rounded checkpoints between $L_1$ and $L_2$, and applied a complex asymptotic approximation for values exceeding $L_2$, utilizing functions such as $D(L_2, z) = 12.6244 \sqrt{z} - 9 \sqrt{L_2}$. By combining these highly localized results and controlling the explicit error variables $E(x; [d_1, d_2])$, Dunn constructed a sharper, completely rigorous global bound under GRH, representing the definitive state-of-the-art constraint on $B_2$.

Advanced Sieving Algorithms and Computational Moduli

To truly appreciate the massive computational bounds reached for calculating $B_2$ at $10^{16}$, one must understand the underlying sieve optimizations utilized by researchers like Pascal Sebah to strip away computational redundancy and memory latency.

Direct software applications of the classical Sieve of Eratosthenes are highly inefficient at this scale because they force the processor to evaluate numbers irrespective of obvious composite traits. To exponentially accelerate the process of locating twin primes, modern sieve algorithms operate entirely within modular arithmetic frameworks utilizing a base $m$ (representing integers in the form $mk + r$, where $0 \leq r < m$). This mathematical representation drastically shrinks the physical subset of integers the CPU must actively load into memory to sieve.

Modulo Representations and the Primorial Function

If the base $m$ is chosen as the product of the first consecutive prime numbers—a value yielded by the primorial function $p\#$—algorithms can pre-eliminate vast swaths of composite numbers without executing a single division operation. The efficiency of this process is dictated by the primorial modulo function $\phi_2(m)$, which outputs the exact number of integer pairs $r$ and $r+2$ that remain relatively prime with $m$.

  • Modulo 6: Because all prime numbers greater than 3 must take the form $6k+1$ or $6k+5$, any search for twin primes can be strictly restricted. The computational sieve restricts itself purely to pairs of the form $(6k+5, 6k+7)$. This elementary optimization instantly drops the proportion of integer space required to process down to $\phi(m)/m = 2/6$, or exactly 33.3%.
  • Modulo 30: Scaling up the primorial base, primes greater than 5 exist solely in the modular forms $30k + r$ (where $r$ is restricted to $\{1, 7, 11, 13, 17, 19, 23, 29\}$). Restricting candidate couples specifically to $(30k+11, 30k+13)$, $(30k+17, 30k+19)$, and $(30k+29, 30k+31)$ drops the actively sieved proportion to $6/30$, or precisely 20%.

This optimization cascades dramatically as the primorial base expands, heavily reducing the strain on processor cache lines and memory bandwidth:

  • $m = 210$ ($2 \times 3 \times 5 \times 7$): The required sieve proportion drops to 14.3%.
  • $m = 2310$ ($210 \times 11$): The required sieve proportion drops to 11.7%.
  • $m = 30030$ ($2310 \times 13$): The required sieve proportion drops to 9.9%.
  • $m = 510510$ ($30030 \times 17$): The required sieve proportion drops to 8.7%.

In the landmark 2002 evaluation of Brun's Constant $B_2$ that successfully breached the $10^{16}$ integer threshold, Pascal Sebah structured his sieve algorithms to execute exactly at modulo 30030 and 510510. By mathematically forcing the algorithm to consider less than 10% of the integer continuum, the computation bypassed crippling floating-point bottlenecks, demonstrating how abstract algebraic structures directly catalyze hardware performance limits.

Prime Quadruplets and $B_4$: Structural and Theoretical Characteristics

Prime Constellations and Quadruplet Topology

While twin primes represent the simplest generalized prime constellation, number theorists dedicate equal scrutiny to denser, highly localized clusters. A prime quadruplet is formally defined as a sequence of four distinct prime numbers separated by the absolute minimum possible distances allowable by arithmetic.

Due to the rigid divisibility rules of integers modulo 3, the shortest possible distance separating four prime numbers cannot be a continuous span of 6 (such as the sequence $p, p+2, p+4, p+6$). In any such sequence, exactly one of the numbers is mathematically guaranteed to be divisible by 3, instantly rendering it composite. The only exception in all of mathematics to this rule is the singular cluster $(3, 5, 7)$, which forms a prime triplet but cannot extend to a quadruplet. Consequently, the absolute minimum distance capable of spanning four prime numbers is 8, forcing the constellation to adopt the strict structural form:

$$(p, p+2, p+6, p+8)$$

With the singular exception of the initial anomalous quadruplet $(5, 7, 11, 13)$, every prime quadruplet in existence must strictly follow the arithmetic progression form $\{30n + 11, 30n + 13, 30n + 17, 30n + 19\}$ for some integer $n$. This inflexible structural constraint guarantees that none of the four prime components are divisible by 2, 3, or 5. In standard base 10 representations, this arithmetic property means the first prime of any quadruplet always terminates in the digit 1, and the final prime of the quadruplet always terminates in the digit 9, a unique formation that is commonly referred to in number theory as a "prime decade". Intriguingly, the sequence of prime quadruplets $\{11, 13, 17, 19\}$ is alleged by some paleoanthropologists to appear carved into the ancient Ishango bone, marking a deep historical fascination with these specific numerical clusters, though this interpretation remains academically disputed.

Distinction from Cousin Primes and Fractal Topologies

It is critical in mathematical literature to strictly disambiguate Brun's Constant for prime quadruplets ($B_4$) from "Brun's Constant for cousin primes." Cousin primes are simply pairs of prime numbers separated by a gap of four, taking the form $(p, p+4)$. While the sum of the reciprocals of cousin primes is also convergant and occasionally designated ambiguously as $B_4$ in overlapping literature, it evaluates entirely differently from the quadruplet sum and represents a distinctly separate topological subset of prime pairs.

The behavior of cousin primes was extensively analyzed by mathematician Marek Wolf, who heuristically estimated the cousin prime constant at roughly 1.19704 utilizing primes up to $2^{42}$. Wolf's analysis revealed extraordinary structural properties when comparing the distribution of cousin primes against twin primes. By plotting the function $W(x)$, which represents the absolute numerical difference between the total count of twin primes and cousin primes up to $x$, Wolf demonstrated that the distribution difference possesses intense fractal properties.

The fractal dimension of the $W(x)$ plot is approximately 1.48. This value is mathematically significant because it is remarkably close to the fractal dimension of a standard Brownian motion path (random walk), suggesting that the localized dominance of twin primes over cousin primes—or vice versa—behaves as a stochastic, random process. Furthermore, Wolf utilized box-counting methods to show that the exact threshold points where the counts of twin primes and cousin primes are perfectly identical also form a fractal set with a dimension of roughly 0.51. These fractal dynamics underscore the profound geometrical complexity hidden within prime gaps, differentiating the analysis of cousin primes entirely from the rigid structure of prime quadruplets.

The Calculation and Density of $B_4$

Evaluating the Reciprocal Sum

Just as Viggo Brun utilized pure sieving techniques to prove the convergence of twin prime reciprocals in 1919, he simultaneously proved that the sum of the reciprocals of prime quadruplets strictly converges to a finite value. This limit defines Brun's Constant for prime quadruplets ($B_4$).

The value $B_4$ is defined mathematically as the infinite sum over all valid prime quadruplets:

$$B_4 = \sum \left(\frac{1}{p} + \frac{1}{p+2} + \frac{1}{p+6} + \frac{1}{p+8}\right)$$

Through massive computational effort mirroring his work on $B_2$, Thomas Nicely evaluated the sum of reciprocals for all prime quadruplets up to extreme hardware thresholds. By carefully managing the floating-point arithmetic to prevent truncation bleed, Nicely determined the constant with an exceptionally tight 99% confidence interval. The value is universally accepted as:

$$B_4 \approx 0.8705883800 \pm 0.0000000005$$

This figure establishes that the total density and inverse sum of prime quadruplets is substantially smaller than that of twin primes (1.902), an expected outcome given the extreme combinatorial rarity required for four primes to cluster sequentially within an integer span of eight.

Hardy-Littlewood Density for Quadruplets

Similar to the modeling of twin primes, the asymptotic density and absolute frequency of prime quadruplets can be accurately predicted using the generalized first Hardy-Littlewood conjecture. The frequency formula for a prime quadruple constellation $P_x(p, p+2, p+6, p+8)$ involves a highly complex multiplicative constant designed to correct for the probabilistic dependence of four localized prime events.

The Hardy-Littlewood constant tailored specifically for prime quadruplets is denoted as $c$. Unlike $C_2$, the constant $c$ is derived from a much more complex infinite product over all primes $p \geq 5$:

$$c = \frac{27}{2} \prod_{p \geq 5} \frac{p^3(p-4)}{(p-1)^4} \approx 4.151180864...$$

Using this exact corrective constant, the asymptotic frequency of prime quadruplets up to a boundary $x$ behaves proportionally to the logarithmic integral of the fourth degree:

$$P_x \sim c \int_{2}^{x} \frac{dt}{(\ln t)^4}$$

This mathematical model accurately mirrors the empirical sparsity of quadruplets as they trend toward infinity, providing the framework needed to extrapolate the trailing tail of the $B_4$ summation far beyond what computers can manually sieve.

Modern Algorithmic Paradigms for Prime Patterns

While Thomas Nicely originally drove enumeration algorithms up to $2 \cdot 10^{16}$ to ascertain the accepted value of $B_4$, computational demands inherently required severe algorithmic breakthroughs to push prime quadruplet tabulation beyond that physical memory threshold. In 2019 and 2020, researchers Jonathan P. Sorenson and Jonathan Webster published highly optimized algorithms designed explicitly for finding primes in generalized linear patterns.

Sorenson and Webster mathematically modeled the search for an admissible linear pattern $P = (f_1(x), \dots, f_k(x))$ of length $k$. By abstracting the sieve to a generalized polynomial format, they could dynamically detect generalized prime sequences—including twin primes where $k=2$, and prime quadruplets where $k=4$. The researchers introduced two primary algorithmic strategies that achieved unprecedented computational complexity limits, bypassing the traditional limitations of raw prime sieves.

  • Algorithm 1 (Time-Optimized): This approach requires $O_P(n / (\log \log n)^k)$ arithmetic operations while utilizing $O(k\sqrt{n})$ space. It achieves this highly accelerated state by extending the Atkin-Bernstein prime sieve and integrating it seamlessly with a space-saving wheel sieve framework, heavily minimizing unnecessary modulus checks.
  • Algorithm 2 (Space-Optimized): Recognizing that physical Random Access Memory (RAM) often throttles prime enumeration long before processor speed does, their second algorithm utilizes slightly more clock time—$O_P(n k / (\log \log n)^{k-1})$ arithmetic operations—but drastically reduces the memory requirements. It operates utilizing only $O(n^{1/c})$ space, where $c > 2$ is a fixed constant.

Deploying parallel implementations of the space-optimized Algorithm 2 across highly distributed computing hardware, Sorenson and Webster successfully found and tabulated every single prime quadruplet up to the massive threshold of $10^{17}$. In doing so, they successfully extended Nicely's original reciprocal sum data for $B_4$, evaluating the partial sums $S_4(X)$ across far higher intervals to verify the stability of the constant.

The results of their tabulation showcase the extreme precision of the $B_4$ convergence as it approaches the $10^{17}$ limit:

Threshold $X$ Quadruplet Count $\pi_4(X)$ Sum of Reciprocals $S_4(X)$
$1 \cdot 10^{16}$ $25,379,433,651$ $0.870477691234...$
$2 \cdot 10^{16}$ $46,998,268,431$ $0.870483710948...$
$3 \cdot 10^{16}$ $67,439,513,530$ $0.870487031043...$
$4 \cdot 10^{16}$ $87,160,212,807$ $0.870489302002...$

These algorithmic advancements establish 0.870588 as a mathematically rigid asymptote, empirically verifying that as the boundary reaches $10^{17}$, the reciprocal sum remains exceptionally tight and predictable, confirming the theoretical dampening modeled by the Hardy-Littlewood integral.

Analytical Implications and Future Trajectories

The ongoing, multi-generational research surrounding the precise calculation of $B_2$ and $B_4$ yields profound secondary insights into the intertwined disciplines of analytic number theory and computer science.

First, the historical interplay between prime number extrapolations and hardware capabilities highlights pure mathematics as an ultimate, unyielding stress-test for micro-architecture. The discovery of the Pentium FDIV bug by Thomas Nicely during his reciprocal sum tabulations serves as a vital historical benchmark. It confirmed that calculating reciprocal aggregates over highly sparse datasets exposes microscopic precision faults within arithmetic logic units that typical software validations and floating-point benchmarks fail to detect. The reciprocal sums of prime constellations will continue to serve as boundary-condition diagnostics for next-generation quantum hardware and highly parallelized computational processors.

Second, the structural reliance on the Generalized Riemann Hypothesis to formulate tighter bounds—such as Lachlan Dunn's 2025 refinement proving $B_2 < 2.1594$—reveals an explicit analytical bridge between the zero-distribution of the Riemann zeta function and the topological density of prime groupings. The fact that the most optimized unconditional upper bound (2.288513) established by Platt and Trudgian remains markedly looser than the GRH-conditional bound demonstrates the specific, quantifiable "shadow" that the unproven Riemann Hypothesis casts over prime distribution geometry. The numerical delta between these two values physically maps the exact degree of statistical noise that proving GRH would permanently eliminate from prime counting functions.

Finally, the shift toward algorithmically generalizing prime detection, as demonstrated by Sorenson and Webster's pattern sieve generating the $B_4$ quadruplet data up to $10^{17}$, signals a departure from single-constellation targeting. As modern space-optimized $O(n^{1/c})$ algorithms become universally adopted, the computational bottleneck for evaluating mathematical constants shifts away from physical memory limitations and strict computational overhead, moving toward fundamental integer factorization boundaries and combinatorial logic limits.

Conclusion

Brun's Constants occupy a uniquely critical junction in modern mathematics, acting simultaneously as a definitive analytical solution and an enduring informational enigma. Viggo Brun's 1919 proof that the reciprocals of twin primes and prime quadruplets strictly converge was instrumental in the maturation of combinatorial sieve theory, firmly establishing that localized prime groupings exist with vastly lower relative densities than the general prime sequence.

However, this strict mathematical convergence effectively blinds mathematicians from deducing the infinitude of these constellations. To extract their true numerical values, researchers cannot rely on brute computational force, as direct sum evaluations for $B_2 \approx 1.902160583$ would theoretically require prime enumerations approaching $10^{530}$ to achieve sufficient precision. Instead, analytic extrapolation through the sophisticated Hardy-Littlewood hypotheses acts as an obligatory mathematical bridge over the chasm of physical impossibility, successfully yielding the modern asymptotic constants $B_2 \approx 1.902160583104$ and $B_4 \approx 0.8705883800$.

The evolution of estimating these constants maps the parallel evolution of computational power. From Glaisher's rudimentary hand calculations in the 19th century, to Nicely's inadvertent auditing of microprocessor logic arrays, to Sorenson and Webster's highly parallel $10^{17}$ space-optimized tabulations, computational rigor has continually scaled. Simultaneously, theoretical constraints have tightened, reaching rigorous unconditional bounds set by Platt and Trudgian and highly optimized GRH-conditional limits derived by Dunn. The sustained pursuit of Brun's Constants continues to define the absolute frontiers of computational endurance, ensuring that prime constellations remain a central pillar of analytic number theory.