Demystify the Prime Number Sequence

"My general view of mathematics is that most of the complicated things we learn have their origins in very simple examples and phenomena." – Dr. Richard Evan Schwartz, Chancellor's Professor of Mathematics, Brown University

"Everything should be made as simple as possible, but not simpler." – Albert Einstein

"Seek simplicity and distrust it." – Alfred North Whitehead

This work is licensed by its author, Gary William Croft, under a Creative Commons Attribution-ShareAlike 3.0 Unported License.

This site explores a deterministic algorithm and geometry in the form of a spiral sieve encompassing eight factorization progressions that intertwine like an octal helix and ultimately determine the distribution of prime numbers greater than 5 (prime numbers defined as whole numbers greater than 1 that are only evenly divisible by 1 and themselves).

The sequence populating the Prime Spiral Sieve can be variously defined as:

♦ Natural numbers not divisible by 2, 3 or 5 (and, given no prime number > 5 is divisible by 2, 3 or 5, it's axiomatic that the domain contains all prime numbers > 5, starting with 7 ... and their multiplicative multiples ...). It follows that all members of our domain are relatively prime (aka coprime or mutually prime) to 2, 3 and 5.

♦ Natural numbers ≡ to {1, 7, 11, 13, 17, 19, 23, 29} modulo 30.

♦ 1 {+6 + 4 + 2 + 4 + 2 + 4 + 6 + 2} {repeat ... ∞}.

♦ 30n+1, 30n+7, 30n+11, 30n+13, 30n+17, 30n+19, 30n+23, 30n+29.

♦ Natural numbers modulo 30 that distribute to the following 8 angles: 12° ... 84° ... 132° ... 156° ... 204° ... 228° ... 276° ... 348°.

♦ The number 1 and all odd numbers >3 with digital root of 1, 2, 4, 5, 7 or 8 and terminating digit of 1, 3, 7 or 9.

♦ It's important to note that the beautiful symmetries encountered within this domain, both numeric and geometric, are largely a consequence of patterns rooted in its period-24 digital root, which has the following repetition cycle: {1, 7, 2, 4, 8, 1, 5, 2, 4, 1, 5, 7, 2, 4, 8, 5, 7, 4, 8, 1, 5, 7, 2, 8} {repeat ...}. This is equivalent to three rotations around the Prime Spiral Sieve.The first 24 members of our domain are especially important in relation to these cycles. For convenience, and to emphasize the periodicity of 24 from whence most of the symmetries we explore emanate, we occasionally frame them as: Numbers ≡ to {1, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89} modulo 90.

Below is a picture of the Prime Spiral Sieve showing rotations spanning numbers 1 thru 259 of the infinite sequence populating the sieve, as defined above:

Hidden deep within this sieve, a radical spin on modulo 30 wheel factorization, are mysteriously beautiful symmetries and geometries, profound in their implications. This sieve not only forms the basis of an extremely efficient and widely used prime number sieving algorithm of the 'non-probabilistic' kind, but it demonstrates how prime numbers are distributed within a radial geometry that effectively defragments the Ulam Spiral and ultimately leads us to "the theory of everything."

This domain, once fathomed, reveals itself to be a beautiful mathematical object in and of itself. It can be conceived as both an infinite spiral and–when matrix factorized at the digital root level–an ever-expanding 4-sided pyramid structured at the deepest level by palindromic sequences. Regardless, the real power of this geometry becomes evident when we triple it dimensionally and explore modulo 90 factorization symmetries at the digital root level, culminating in the
Magic Mirror Matrix, a 'calculatory geometry' that serves as a prime factorization sequencer accounting for the first 1000 prime numbers, *exactly*, and ultimately all prime numbers >5.

Most profoundly, we'll discover how the eight spiraling factorization algorithm's generating composite numbers within the Prime Spiral Sieve's orbits, i.e., the multiplicative multiples of the domain's sequence starting with 7^{2}, are structured by rotational symmetry groups in the shape of equilateral triangles ( {1,4,7} {2,5,8} {3,6,9} ) that form 3 x 3 Latin Squares that in turn extrapolate into beautiful complex polygons.

To help visualize these pattern generating algorithms, imagine a musical instrument with 24 strings possessing a neck of infinite length upon which are played a score of period-24 chord progressions comprised exclusively of composite numbers ... the unplayed 'notes' being all prime numbers greater than 5 ...]

Only elementary arithmetic and geometry are required to understand what you find here. Our approach, rooted in empirical observation, is more expository than technical, analogous to "writing through the curriculum." Testimony supporting this 'languaging' of math comes from Stephen Hawking in his *Brief History of Time: A Reader's Companion* (1992): "Equations are necessary if you're doing accountancy, but they are the boring part of mathematics. Most of the interesting ideas can be conveyed by words or pictures." [Some equations, of course, are beautiful in their own right, as Dr. Hawking would no doubt acknowledge.] Though our method requires minimal mathematical knowledge, it can nonetheless deliver insight and aesthetic pleasure.

What follows is a mixture of well known fact and conjecture (the latter so labeled) informed by more than two decades of heuristic
experimentation complemented by research. With one exception, our Proof by Construction for the Digital Root Sequencing of Twin Primes, you'll be disappointed if you're seeking rigorous proofs here. Yes, we're well aware of the mantra
"Observations Aren't Proofs." However, we hope you'll concede that "observations precede proofs," the Reimann Hypothesis being the most famous example. Speaking of which, Karl Sabbagh, writing in *The Riemann Hypothesis: The Greatest Unsolved Problem in Mathematics*, quotes highly regarded mathematician Hugh Montgomery as saying: "Philosophically, I'm a firm believer in trying to form conjectures far beyond what one can prove ..." Or, as Marcus Du Sautoy put it, prefacing a quote by Richard Feynman, "[A]s a physicist I hold Feynman's dictum dear: 'A great deal more is known than has been proved.'" In that spirit, if your intent is to better *understand* how prime numbers are sequenced, read on. Through words, illustrations and elementary arithmetic, we'll show how prime numbers are sequenced deterministically; their seeming randomness an illusion.

If this site "proves" anything, it's that those with only an elementary education in mathematics can have direct access to the beautiful symmetries encompassing algorithmic order embedded within the seeming chaos of prime number and related patterns. It's hard to imagine a presumably enlightened "Spirit of the Universe" barring the vast majority of intelligent beings who haven't pursued higher mathematics a measure of insight into the deep mystery "curled up" in the roots of the number universe. After all, the experience is exhilarating, if not spiritualizing.

To date this site has attracted more than 58,900 "absolute unique visitors" (Google Analytics speak) from 194 countries. Scores have returned repeatedly, and none have offered counter-examples. From a hermeneutic perspecitve, the patterns explored herein speak for themselves ...*Res ipsa loquitur ...*

[Note: Twelve of the sequences discussed on this site have been published by the author on the *On-Line Encyclopedia of Integer Sequences*:

A227863: Numbers congruent to {1,49} mod 120.

A227896: Digital root of Fibonacci numbers indexed by natural numbers not divisible by 2, 3 or 5.

A230462: Numbers congruent to {1, 11, 13, 17, 19, 29} mod 30.

A232878: Twin prime pairs which sum to perfect squares.

A232880: Twin primes with digital root 2 or 4.

A232881: Twin primes with digital root 5 or 7.

A232882: Twin primes with digital root 8 or 1.

A233766: Digital root of Lucas numbers indexed by natural numbers not divisible by 2, 3 or 5.

A230113: Digital root of Fibonacci plus Lucas digital root indexed by numbers not divisible by 2, 3 or 5.

A240924: Digital root of squares of numbers not divisible by 2, 3 or 5.

A246508: Digital root of numbers ≡ to {1, 7, 11, 13, 17, 19, 23, 29} mod 30.

A261371: Prime numbers that sum to a prime when added to their digital root.

Also note that an M.I.T. licensed (Python) Prime Factorization Tool compared seven non-probabilistic prime number sieving algorithms, and the programmer deemed the Croft Spiral (aka Prime Spiral Sieve) the 'fastest and most efficient' of those tested. Quoting the programmer, "The fastest method, Croft, is over 1000 times faster than the slowest."]

The genesis of most if not all repeating prime number patterns described in the mathematics literature, e.g., Twin Primes, Cullen Primes,
Chen Primes, Sexy Primes, Cousin Primes, Sophie Germain Primes, Siamese Primes, Cunningham Chains ... the list goes on and on ... can be readily deciphered
using the Prime Spiral Sieve as an analytical tool employing modular arithmetic
(and specifically, modulo 30 relationships). Here are two examples supporting this claim, i.e. using this sieve to analyze and predict
Siamese Primes (n^{2} -2 and n^{2} +2 are primes) and
Sophie Germain Primes (p and 2p+1 are primes), keeping in mind that these will make more
sense after you've read what follows.

The most obvious of these repeating prime patterns are the three Twin Prime Distribution Channels, described at length on this site. These and all other such repeating–albeit intermittent and seemingly random–patterns are fundamentally sub-patterns of the set of natural numbers not divisible by 2, 3 or 5 when arrayed in 8 dimensions, whether in a matrix or spiral form. The illusion of randomness results from the overlapping sequences of the eight algorithmic "chord progressions" that factorize the domain. We will be discussing these progressions when we get to prime factorization.

The first rotation of the sieve, comprised of 8 members, (1, 7, 11, 13, 17, 19, 23 and 29), is the deterministic key to everything that follows. These are the first 8 counting numbers not divisible by 2, 3 or 5, a sequence which by definition includes (and *only* includes) 1 and all primes ≥ 7 and their multiplicative multiples (and, as you'll see below, it's conjectured that the entire set can be generated by a simple expression involving 2, 3 and 5).

[Note: Given our domain is limited to numbers ≡ to {1,7,11,13,17,19,23,29} modulo 30," only ϕ(m)/m = 8/30 or 26.66% of natural numbers need be sieved. Also note that if you plug the number 30 into Euler's totient function, phi(n): phi(30)= 8, with the 8 integers (known as totatives) smaller than and having no factors in common with 30 being: 1, 7, 11, 13, 17, 19, 23 and 29, i.e., what are called "prime roots" above. Thirty is the largest integer with this property.]

The integer **30**, product of the first three prime numbers (2, 3 and 5), and thus a primorial, plays a powerful role organizing the array's perfect symmetry, viz., in the case of the 8 prime roots:

**1+29=30; 7+23=30; 11+19=30;** and **13+17=30.**

In *The Number Mysteries* well-known physicist and mathematics popularizer Marcus Du Sautoy writes: "In the world of mathematics, the numbers 2, 3, and 5 are like hydrogen, helium, and lithium. That's what makes them the most important numbers in mathematics." Although 2, 3 and 5 are the only prime numbers not included in the domain under discussion, they are nonetheless integral to it: First of all, they sieve out roughly 3/4ths of all natural numbers, leaving only those nominally necessary to construct a geometry within which prime numbers can be optimally arrayed. The remaining 26.66% (to be a bit more precise) constituting the array can be constructed with an elegantly simple interchangeable expression that incorporates the first three primes. It's conjectured that this expression can be configured (albeit by trial-and-error) to produce **all** (*and only*) the numbers in the array (and their negatives):** x ^{n}y^{n} ± z^{n}** where x=2, y=3 and z=5. Thus:

Given that all prime numbers > 5 are in the array, it is conjectured that this expression can be configured to generate all primes >5. What is critical to understand, is that the invisible hand of 2, 3 and 5, and their factorial 30, create the *structure* within which the balance of the prime numbers, i.e., all those greater than 5, are arrayed algorithmically–as we shall demonstrate. Primes 2, 3 and 5 play out in modulo 30-60-90 cycles (decomposing to {3,6,9} sequencing at the digital root level). Once the role of 2, 3 and 5 is properly understood, all else falls beautifully into place.

The Prime Spiral Sieve possesses remarkable structural and numeric symmetries. For starters, the intervals between the prime roots (and every subsequent row or rotation of the sieve) are perfectly balanced, with a period 8 difference sequence of: {6, 4, 2, 4, 2, 4, 6, 2}. The entire domain can thus be defined as **1 {+6 +4 +2 +4 +2 +4 +6 +2} {repeat ... ∞}**. As we've already suggested, the number 30 figures large in our modulo 30 domain. The Prime Spiral Sieve is Archimedean in that the separation distance between turns equals 30, ad infinitum. The first two rotations increment as follows:

Interestingly, the sum of the 2nd rotation = 360. Is it coincidental that the product of the first three primorials, 2, 6 and 30 = 360? Or is it coincidental that when you multiply the first five Fibonacci numbers in sequence, you produce 1, 2, 6 and 30? And, speaking of the Fibonacci number sequence, there is symmetry mirroring the above in the relationship between the terminating digits of Fibonacci numbers and their index numbers equating to members of the array populating the Prime Spiral Sieve:

Remarkably, the sequence of Fibonacci terminating digits indexed to the prime roots,13,937,179 (see graphic, above), is a prime number and a member of a prime pair (with 13,937,177), though, if you're curious, not a reversible prime (although the reversal is a semi-prime: 9,461 * 10,271 = 97,173,931, and you'll note that both its prime factors have two combinations summing to 10). In addition, 13,937,179 when added to its reversal 97,173,931 = 111,111,110 (in strict digital root terms, the sum is 11,111,111) and the entire repeating (and palindromic) Fibo sequence end-to-end (equivalent to two rotations around the sieve) gives you this equivalency: 1,393,717,997,173,931 ≡ 11,111,111 (mod 111,111,110)... (and interestingly, 11,111,111 * 111,111,110 = 1234567876543210 and 111,111,110/11,111,111 = 10). Also 1,393,717,997,173,931 is divisible by the repunits 11 and 1,111 and 11,111,111.

Echoing the Fibonacci patterns just described, the terminating digits of the *prime roots* (17,137,939), when added to *their*
reversal (93,973,171) = 111,111,110. And, when you connect the prime root terminating digit sequence to its reversal,
the entire palindromic sequence end-to-end produces this: 1,713,793,993,973,171 ≡ 111,111,111 (mod 111,111,110)
[And in this case, 111,111,111 * 111,111,110 = 12345678876543210.]. And if that isn't enough, 1,713,793,993,973,171 is *also* divisible
by the repunits 11 and 1,111 and 11,111,111.

Well, not quite enough, because there's yet another related dimension of symmetry: The terminating digits of the prime root
*angles* (24,264,868; see illustration of Prime Spiral Sieve) when added to *their*
reversal (86,846,242) = 111,111,110, not to mention this sequence possesses symmetries that dovetail perfectly with the prime root and
Fibo sequences, including the fact that when it is connected to its reversal (giving us 2,426,486,886,846,242), it's divisible by the
repunits 11 and 1,111 and 11,111,111.

And when you combine the terminating digit symmetries described above, capturing three rotations around the sieve in their actual sequences, you produce the ultimate combinatorial symmetry:

Here's yet another fascinating dimension of symmetry: the pattern of 9's created by decomposing and summing either the digits of
Fibonacci numbers indexed to the first two rotations of the spiral (a palindromic pattern {1393717997173931} that repeats every 16 Fibo
index numbers) or, similarly, decomposing and summing the prime root angles. The decomposition works as follows (in digit sum arithmetic
this would be termed summing to the digital root): F_{17} (the 17th Fibonacci number) = 1597 = 1 + 5 + 9 + 7 = 22 = 2 + 2 = 4:

Another dimension of symmetry involves the terminating digits of the prime roots and their angles: those paired with like terminating digits being separated by 120°: 1(12°) and 11(132°) ... 13(156°) and 23(276°) ... 7(84°) and 17(204°) ... 19(228°) and 29(348°). Another consideration with regard to terminating digits, is that one can easily construct, by combining all numbers with the same terminating digits, a four-fold arithmetic progression in increments of +10 and +20, starting with 1, 7, 13 and 19. Thus, combining 1(12°) and 11(132°) gives us: 1, 11, 31, 41, 61, 71, 91, [+10+20] ... n; combining 7(84°) and 17(204°) gives us 7, 17, 37, 47, 67, 77, 97, [+10+20] ... n; combining 13(156°) and 23(276°) gives us 13, 23, 43, 53, 73, 83, 103, [+10+20] ... n; and, combining 19(228°) and 29(348°) gives us 19, 29, 49, 59, 79, 89, 109, [+10+20] ...n. Looking at the array in this configuration, however, has borne no fruit.

As fascinating as the symmetries examined above may be, they are but a prelude to the beautiful patterns we'll explore when we discuss digital root sequencing and the Trinity of Triangles and Magic Squares rooted in Vedic Arithmetic that drive factorization algorithms within this domain. And, finally, if you want to jump ahead and view the most stunning symmetrical object found on this site, check out the Magic Mirror Matrix that maps factorizatons at the digital root level and accounts for the first 1000 prime numbers–exactly.

Around the perimeter of the spiral sieve (pictured below) you'll note that the 8 radii are labeled in relation to their modulo 30 prime roots, i.e., 1(12°); 7(84°); 11(132°); 13(156°); 17(204°); 19(228°); 23(276°) and 29(348°). These relate to the fact that the circle is segmented into 30 equal sectors or radii separated by 12° (30*12°=360°), although only the eight radials that are the focus of this study are shown.

This sieve "exposes" the twin primes, aligning as they do along three distinct "distribution channels." One obvious implication, is that those numbers in the array congruent to {7} modulo 30 (radial angle 84°) and {23} modulo 30 (radial angle 276°) can be excluded as twin prime candidates (and, by definition, all prime numbers distributed along these two diagonals. with the exception of 7, which is twinned with 5, will be what are known as "isolated primes"). Later we explain how twin prime candidates can be segregated from all other positive integers and be partitioned into three columnar sets covertly aligned by the first three prime numbers (encoded in angles).

The array is rooted in the first three prime numbers: 2, 3 and 5 and their product, 30, the 3rd primorial. This array reveals that the first three primes play a very special role in creating the symmetrical geometries that align the distribution of all subsequent prime numbers, thus distinguishing them from all other primes. Primes 2, 3 and 5 are like 8-legged spiders assigned to spin the beautiful spiraling web in which the remaining prime numbers are arrayed along assigned threads. (For a detailed listing of Number 30's attributes, plus reference links click here: The Number 30).

It is conjectured that all (and only) the numbers in this array (and their negatives) can be derived using the
interchangeable expression incorporating the first three prime numbers, 2, 3 and 5, where x=2, y=3 and z=5.
Thus: ** x ^{n}y^{n} ± z^{n}**,

All prime numbers (with the exception of 2, 3 and 5) are distributed along 8 diagonals in intervals of 30, starting with "prime roots": 1, 7, 11, 13, 17, 19, 23 and 29 (thus: 1...31...61...91...n; 7...37...67...97...n; etc.).

The products of *any* combination of factors in the array = a number in the array, e.g.,
7*11 = 77; 7*11*13 = 1001; etc. Conversely, all factors for composite numbers in the array can be found in the array.

Every composite number in our modulo 30 domain can be derived from the product of two terms in the domain multiplied together, and these multipliers need not necessarily be prime themselves. For example, 5831, which is congruent to 11, modulo 30, and therefore in the array, is the product of 49 x 119 = 5831. In this example, neither 49 (7 x 7) nor 119 (7 x 17) are prime, though both are members of the array.

The sum of *any* sequential odd number of addends in the array = a number in the array, e.g.,
1+7+11 = 19; 1+7+11+13+17 = 49; etc.

Because the digital roots of all prime root angles are either 3 or 6, any prime root angle times another will produce a product whose digital root = 9, e.g., PR7 (84°) x PR29 (348°) = 84 x 348 = 29232 = dr(9).

Any number in the array x 30 + 1 = a number in the array.

The sum of the angles for 2(24°), 3(36°) and 5(60°) = 120°, and the sum of the prime roots (1+7+11+13+17+19+23+29) also = 120. This is because the prime roots are an arithmetic anagram for the angles of the first three primes, thus: 11+13 = 24; 17+19 = 36; and 1+7+23+29 = 60. The sum of the second rotation = 360 ... 3(2[24°] + 3[36°] + 5[60°]) = 30[360°]

The array reveals beautifully symmetrical relationships:

1[12°] + 29[348°] = 30[360°]

7[84°] + 23[276°] = 30[360°]

11[132°] + 19[228°] = 30[360°]

13[156°] + 17[204°] = 30[360°]

Mod 30 of all numbers in this array (and thus all primes other than 2, 3 and 5) must be 1, 7, 11, 13, 17, 19, 23 or 29.

The sum of the digital root sums of the prime roots (1, 7, 11, 13, 17, 19, 23, 29) = 1+7+2+4+8+1+5+2 = 30.

This sieve reveals why all primes >5 are adjacent to a multiple of six, as the prime root radii are adjacent to 6(72°); 12(144°); 18(216°); 24(288°); and 30(360°). [And you'll note that the digital root sums of all adjacent angles equal 9.]

The modulo 90 congruence of any member of this domain can be determined by digital root and terminating digit configured in an xy matrix: arrayed by digital root {1, 2, 4, 5, 7, 8} on the vertical axis and by terminating digits {1, 3, 7, 9} on the horizontal axis.

By definition prime factors for all composite numbers within this domain must originate from within it, which is why all composite numbers are reducible to one or more modulo 90 'congruency dyads' traceable to the first 24 members of this domain:

For a detailed discussion of efficient factorization and prime number sieving algorithms, as well as an in-depth analysis of
the 8-chord progression and deterministic modulo 90 digital root dyad sequences underlying all factorizations employing this sieve,
click here: **Prime Number Sieving Algorithms**.

If you use Python (the computer programming language designed for the development of scientific, engineering and
mathematics applications) and want to cut to the chase, check out the MIT licensed Python module dubbed "pyprimes" designed to run and compare *non-probabilistic*
prime number sieving algorithms, including the Sieve of Eratosthenes and the Prime Spiral Sieve (referred to by its alternative name, "Croft Spiral Sieve"),
here or here at
code.google.com.
The programmer, we're pleased to report, rates the Prime Spiral Sieve "fastest/best" and recommends it as "the preferred way of generating
prime numbers" compared to the several sieves tested.

Although "pyprimes" is an efficient algorithm, we're waiting for a savvy programmer to appreciate that processing speed can be increased *exponentially* by plotting loci within a 24-wide matrix generated by 24 each period-24 deterministic modulo 90 progressions (which can run in series-parallel, for even greater speed). Once the 'loci chord' progression algorithm(s) are set up and initialized, *the need for multiplication and/or division is eliminated*. This method involves *pattern generation*–not number crunching. This is neither factorization nor primality testing, per se, but rather an extremely efficient deterministic (not probabilistic) way to index composite numbers within this domain, i.e., numbers not divisible by 2, 3 or 5 when framed congruent to {1, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89} modulo 90). The matrix can be populated numerically *after-the-fact*, all composite numbers having been located by the algorithmically generated loci patterns which constitute an array of points on a modulo 90 matrix. When the process of pattern superimposition is complete, all primes greater than 5 within the domain up to a given n are completely differentiated from composite numbers.

The world of twin primes is truly magical. For a detailed discussion of the factorization algorithms and symmetry groups (tensor matrices, orthogonal Latin squares, palindromes, equilateral triangles, star polygons, etc.) that ultimately determine the distribution of twin primes along three modulo 30 'channels' go to this page: **Twin Primes Demystified: Prime Pairing Distribution Algorithms and Symmetries**.

There you will be introduced to the beautiful 'Palindromagon' (pictured below), a complex polygon generated by tiered digital root dyad cyles central to the twin prime distribution channels (as well as modulo 90 factorization sequences). Its name comes from the fact that the triangulations generating it sum to a period-18 palindrome consisting of the six possible permutations of {3,6,9}, which in turn can be permutated to produce two 3 x 3 Latin squares with rows, columns and principal diagonals all summing to 18:

Although not particularly elegant, there is a method whereby one can calculate the number of prime numbers in a given range with considerable precision using this Sieve. To illustrate, below is a step-by-step procedure to determine the number of primes in the range 1 thru 10,000 [Note: Please consider this section to be more conceptual outline than final product in as much as the author will not be satisfied until the process is drastically improved.]:

- First, calculate the number of integers within the range not divisible by 2, 3 and 5. Since 8 out of every 30 integers fall into this category, we have 8/30 = .2666666666... or 26.66%. Multiplying 10,000 times .266666... gives us 2666. (This is our starting point. The steps that follow will identify the number of non-primes to subtract until we've boiled 2666 down to reflect only prime numbers.)
- Next, given that the square root of 10,000 is 100, we know that we need to perform chord factorizations from 7chord to 97chord to generate products (i.e., non-primes) less than or equal to 10,000 (reference "The 8-Chord Progression," above). Doing so will index all composite numbers in this array (aside from those divisible by 2, 3 and 5, already removed in step 1).
- We then count the number of calculations required to complete step 2, above. This gives us 1847. (To open an Excel spreadsheet showing in detail how this number is derived, click here. This may take several seconds!)
- Because the factorization process generates a number of duplicate products, we must count and subtract these from the total in step 3. [Note: These duplications were described earlier, and the example given is repeated here: 7*77 = 539 is equivalent to 11*49 = 539; both expressions being equivalent to 7*7*11 = 539.]. In our example, we find that there are 411 such duplications; subtracting these from the number of factorizations (1847-411) gives us 1436. To open an Excel spreadsheet listing the duplicates and their count, click here.]
- Next, we subtract the result in step 4 (1436) from the result in step 1 (2666), i.e., 2666-1436 = 1230.
- The final step is to subtract 1 from the total in step 5 to account for the fact that the number 1 is not a prime; this
leaves us with a balance of 1229, which is
*exactly*the number of primes between 1 and 10,000.

No doubt this procedure can be simplified. For example, one could probably devise formulas to determine the number of factorizations and duplicates, eliminating the need for a spreadsheet count.

All perfect squares within our domain (numbers not divisible by 2, 3 or 5) possess a digital root of 1, 4 or 7 and are congruent to either {1} or {19} modulo 30. By definition, this includes the squares of all prime numbers greater than 5. We can easily explain this from a digital root perspective given that the digital roots of members of our domain are restricted to 1, 2, 4, 5, 7 or 8 (Numbers with digital root 3, 6, or 9 can't be members because they are divisible by 3.). Thus the digital root of squares is likewise restricted, as follows (and note the palindrome):

1 x 1 = 1

2 x 2 = 4

4 x 4 = 7

5 x 5 = 7

7 x 7 = 4

8 x 8 = 1

By arithmetic law, perfect squares can only have terminating digits of 1, 4, 5, 6, 9 or 0. Only two of these final digits (1 and 9) apply to our domain, i.e., for numbers congruent to {1, 11, 19 or 29} modulo 30. In turn, numbers congruent to {11, 29} sequence digital roots 2, 5 or 8, and therefore – as we demonstrated above – there can be no perfect squares among them. And so it is that the distribution of squares is narrowed to numbers congruent to {1, 19} modulo 30, which is to say they distribute along two – and only two – radii of the Prime Spiral Sieve: 12° (numbers congruent to {1} modulo 30) and 228° (numbers congruent to {19} modulo 30). (This is also consistent with the fact that the quadratic residues for modulo 30 (making them congruent with perfect squares) are 1 and 19.)

[And it follows that all squares in this series distribute evenly to two of the three twin prime distribution channels, described above, negating a significant percentage of potential twin prime pairs.]

The matrix below illustrates the distribution of squares from 1*1 thru 59*59 (squares hi-lited in blue):

Summarizing the above relationships in mathematical terms (and in the knowledge that these modular relationships apply to the squares of all prime numbers ≥ 7) we get:

for all n where n mod 30 = 1, n^{2} mod 30 = 1

for all n where n mod 30 = 29, n^{2} mod 30 = 1

for all n where n mod 30 = 11, n^{2} mod 30 = 1

for all n where n mod 30 = 19, n^{2} mod 30 = 1

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

for all n where n mod 30 = 7, n^{2} mod 30 = 19

for all n where n mod 30 = 23, n^{2} mod 30 = 19

for all n where n mod 30 = 13, n^{2} mod 30 = 19

for all n where n mod 30 = 17, n^{2} mod 30 = 19

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

When the digital root of perfect squares is sequenced within a modulo 30 x 3 = modulo 90 horizon, beautiful symmetries in the form of 24 period repeating palindromes are revealed, which the author has documented on the *On-Line Encyclopedia of Integer Sequences* as Digital root of squares of numbers not divisible by 2, 3 or 5 (A24092):

1, 4, 4, 7, 1, 1, 7, 4, 7, 1, 7, 4, 4, 7, 1, 7, 4, 7, 1, 1, 7, 4, 4, 1

In the matrix pictured below, we list the first 24 elements of our domain, take their squares, calculate the modulo 90 congruency and digital roots of each square, and display the digital root factorization dyad for each square (and map their collective bilateral 9 sum symmetry):

The Ulam Spiral arrays prime numbers in fragmented spiral and diagonal formations. Quoting from Wikipedia: "Since in the Ulam spiral adjacent diagonals are alternatively odd and even numbers, it is no surprise that all prime numbers lie in alternate diagonals ... What is startling is the tendency of prime numbers to lie on some diagonals more than others." From this one might deduce that the Ulam Spiral is very likely a scrambled version of the Prime Spiral Sieve as the latter demonstrates how all prime numbers (except 2, 3 and 5) are fundamentally arrayed along eight (and only eight) diagonals.

It would appear, circumstantially, that the Prime Spiral Sieve is mathematically harmonious and perhaps isomorphic with the most complex and visually arresting Lie group, named E_{8}, which–like the Prime Spiral Sieve–is 8-dimensional (E_{8} is pictured below superimposed with a star polygon and the 8 radii of the modulo 30 factorization wheel). This group was recently in the news as possibly being a key to unifying theories in gravity and particle physics to create the proverbial "theory of everything." The number 30–integral to the Prime Spiral Sieve–is the Coxeter Group number *h*, dual Coxeter number and the highest degree of fundamental invariance of E_{8}. You'll note, looking at the graphical representation of E_{8} below, that the perimeters of every one of its multiple concentric circles possesses 30 points. And, not surprisingly, E_{8} has 2-, 3- and 5-torsion and its exponents are the co-primes up to 30, i.e., 1, 7, 11, 13, 17, 19, 23, and 29–numbers you're very familiar with if you've read to this point ... which brings us
full circle Ο:

We'll close with a graphic showing E_{24} superimposed with the 24 radials of a modulo 90 factorization wheel and the 15 points of a 15-point star represented with red dots, each point separated by 24°:

Your feedback welcome! Email: gwc@hemiboso.com

This work is licensed by its author Gary William Croft under a Creative Commons Attribution-ShareAlike 3.0 Unported License.