23 February 2011

Ulam spiral

The Ulam spiral, or prime spiral (in other languages also called the Ulam Cloth) is a simple method of visualizing the prime numbers that reveals the apparent tendency of certain quadratic polynomials to generate unusually large numbers of primes. It was discovered by the mathematicianStanislaw Ulam in 1963, while he was doodling during the presentation of a “long and very boring paper”[1] at a scientific meeting. Shortly afterwards, in an early application of computer graphics, Ulam with collaborators Myron Stein and Mark Wells used MANIAC II at Los Alamos Scientific Laboratory to produce pictures of the spiral for numbers up to 65,000.[2][1][3] In March of the following year, Martin Gardner wrote about the Ulam spiral in his Mathematical Games column;[1] the Ulam spiral featured on the front cover of the issue of Scientific American in which the column appeared.

Ulam constructed the spiral by writing down a regular rectangulargrid of numbers, starting with 1 at the center, and spiraling out:

He then circled all of the prime numbers and he got the following picture:

To his surprise, the circled numbers tended to line up along diagonal lines. The image below is a 200×200 Ulam spiral, where primes are black. Diagonal lines are clearly visible, confirming the pattern. Horizontal and vertical lines, while less prominent, are also evident.

Ulam spiral of size 200×200

All prime numbers, except for the number 2, are odd numbers. 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 of the Ulam spiral. What is startling is the tendency of prime numbers to lie on some diagonals more than others.

Tests so far confirm that there are diagonal lines even when many numbers are plotted. The pattern also seems to appear even if the number at the center is not 1 (and can, in fact, be much larger than 1). This implies that there are many integer constants b and c such that the function:

f(n) = 4n2 + bn + c

generates, as n counts up {1, 2, 3, ...}, a number of primes that is large by comparison with the proportion of primes among numbers of similar magnitude.

Remarkably, in a passage from his 1956 novel The City and the Stars, author Arthur C. Clarke describes the prime spiral seven years before it was discovered by Ulam. Apparently, Clarke did not notice the pattern revealed by the Prime Spiral because he "never actually performed this thought experiment."[4]

Contents

In their famous 1923 paper on the Goldbach Conjecture, Hardy and Littlewood stated a series of conjectures, one of which, if true, would explain some of the striking features of the Ulam spiral. This conjecture, which Hardy and Littlewood called “Conjecture F”, is a special case of the Bateman–Horn conjecture and asserts an asymptotic formula for the number of primes of the form ax2+bx+c. Rays emanating from the central region of the Ulam spiral making angles of 45° with the horizontal and vertical correspond to numbers of the form 4x2+bx+c with b even; horizontal and vertical rays correspond to numbers of the same form with b odd. Conjecture F provides a formula that can be used to estimate the density of primes along such rays. It implies that there will be considerable variability in the density along different rays. In particular, the density is highly sensitive to the discriminant of the polynomial, b2−16c.

The primes of the form 4x2-2x+41 with x=0, 1, 2, ... have been highlighted. The prominent parallel line in the lower half of the figure corresponds to 4x2+2x+41 or, equivalently, to negative values of x.

Conjecture F is concerned with polynomials of the form ax2+bx+c where a, b, and c are integers and a is positive. If the coefficients contain a common factor greater than 1 or if the discriminant Δ=b2−4ac is a perfect square, the polynomial factorizes and therefore produces composite numbers as x takes the values 0, 1, 2, ... (except possibly for one or two values of x where one of the factors equals 1). Moreover, if a+b and c are both even, the polynomial produces only even values, and is therefore composite except possibly for the value 2. Hardy and Littlewood assert that, apart from these situations, ax2+bx+c takes prime values infinitely often as x takes the values 0, 1, 2, ... This statement is a special case of an earlier conjecture of Bunyakovsky and remains open. Hardy and Littlewood further assert that, asymptotically, the number P(n) of primes of the form ax2+bx+c and less than n is given by

where A depends on a, b, and c but not on n. By the prime number theorem, this formula with A set equal to one is the asymptotic number of primes less than n expected in a random set of numbers having the same density as the set of numbers of the form ax2+bx+c. But since A can take values bigger or smaller than 1, some polynomials, according to the conjecture, will be especially rich in primes, and others especially poor. An unusually rich polynomial is 4x2−2x+41 which forms a visible line in the Ulam spiral. The constant A for this polynomial is approximately 6.6, meaning that the numbers it generates are almost seven times as likely to be prime as random numbers of comparable size, according to the conjecture. This particular polynomial is related to Euler's prime-generating polynomialx2−x+41 by replacing x with 2x, or equivalently, by restricting x to the even numbers. Hardy and Littlewood's formula for the constant A is

In the first product, p is a prime dividing both a and b; in the second product, is an odd prime not dividing a. The quantity ε is defined to be 1 if a+b is odd and 2 if a+b is even. The symbol is the Legendre symbol. A quadratic polynomial with A ≈ 11.3, currently the highest known value, has been discovered by Jacobson and Williams.[6][7]

Robert Sacks devised a variant of the Ulam spiral in 1994. In the Sacks spiral the non-negative integers are plotted on an Archimedean spiral rather than the square spiral used by Ulam, and are spaced so that one perfect square occurs in each full rotation. (In the Ulam spiral, two squares occur in each rotation.) Euler's prime-generating polynomial, x2−x+41, now appears as a single curve as x takes the values 0, 1, 2, ... This curve asymptotically approaches a horizontal line in the left half of the figure. (In the Ulam spiral, Euler's polynomial forms two diagonal lines, one in the top half of the figure, corresponding to even values of x in the sequence, the other in the bottom half of the figure corresponding to odd values of x in the sequence.)