I wrote a little bit about the last digits of consecutive twin primes yesterday in this blog post:

The distribution of the last digits of consecutive twin primes is strange

My data is in this spreadsheet:

My google doc with twin prime counts – see the “Improved Twin Prime Sheet” tab

I’m looking at these patterns because of the new result about prime numbers discussed in this Evelyn Lamb article:

Evelyn Lamb’s article about the new result about last digits of prime number

Last night – just by coincidence – I noticed that the value of one of the error terms I was wondering about was surprisingly close to 1/(100*e). I was playing with the numbers up to the 5 billionth prime yesterday. The numbers up to the 6 billionth prime made this error term even a bit closer to 1/(100*e).

I have a hard time believing that the error term I was noticing can be represented in such a simple way, but for now . . . here is my guess at the count for the various digit sequences of consecutive twin primes up to a number X.

Let be the number you are counting up to.

Let be the number of twin primes up to

Let be the number

Let be shorthand for

The the number of consecutive twin primes with the following last digit sequences is approximately:

{(1,3) (1,3)} : (N/9)*( 1 – – 9* )

(1,3) (7,9) : (N/9)*( 1 + +9* )

(1,3) (9,1) : (N/9)*( 1 )

(7,9) (1,3) : (N/9)*( 1 + )

(7,9) (7,9) : (N/9)*( 1 – – 9* )

(7,9) (9,1) : (N/9)*( 1 + 9* )

(9,1) (1,3) : (N/9)*( 1 + 9* )

(9,1) (7,9) : (N/9)*( 1 )

(9,1) (9,1) : (N/9)*( 1 – 9* )

Comparing the actual counts to the formula up to 1 billionth, 2 billionth, and up to the 6 billionth prime you get these results (from the google doc)

Actual Counts:

The formulas above give these estimates:

and the errors between the formula and the actual count are:

The strange thing about this result is that it appears that there are more pairs of consecutive twin primes with last digitsĀ (9,1) (9,1) than there are with either last digits (1,3) (1,3) or (7,9) (7,9).

The difference is that strange factor of 1 – 1/(100*e).