1.06. Skewes' Numbers
"Even if we had a googol googol googol googol googol ... ... ... ... ... ... ... ... ... ... ... ... googol googol googol googol (say 10^31 times) observable universes you still wouldn't have enough particles to write out Skewes' Number!!! Now THAT is a proper description." ~ Sbiis Saibian (source)
Introduction
We've looked a variety of large numbers so far: numbers that exist within the real world, numbers encountered in science, and numbers used to tell the world about large numbers. But now, it's time to look at two famous numbers that appeared in mathematics: the Skewes' numbers.
The Skewes' numbers are large upperbounds to the solution of a problem whose answer is still not known, and they were named after Stanley Skewes who proved them to be upperbounds. The smaller bound assumed the Riemann hypothesis to be true, and the larger one about 20 years later did not assume it. Although the upperbounds have been improved, the numbers have gotten quite some fame due to their size. In this article I'll talk about how Skewes' numbers came to be, then the numbers themselves and how big they are, and finally the further history of the problem in the years following Skewes' proof.
Who Is Stanley Skewes?
Stanley Skewes was a South African mathematician (born 1899, died 1988) who, according to Wikipedia,^{[}^{1}^{]} is not known for much other than the numbers we'll talk about in this article. I find it interesting that such a humble figure played a notable part in the world of large numbers and even broke a record related to them; it goes to show how niche of a field googology is. Skewes was a student of a more famous mathematician named John Edensor Littlewood, who did a lot of important research in the mysterious workings of prime numbers, and whose work we will further discuss in this article. You may already know that prime numbers are notorious for not having a discernible pattern of occurrence, which causes allure and frustration for mathematicians alike. Stanley Skewes also worked with the study of prime numbers quite a lot, which leads us to how his eponymous numbers came to be.
Skewes' Problem
The problem Skewes was working on was this question:
What is the smallest natural number x such that π(x) > li(x), where π(x) is the prime counting function and li(x) is the logarithmic integral function?
Both π(x) and li(x) are simple functions which serve as approximations of each other, and we'll take a closer look at what each of them are.
The PrimeCounting Function
π(x) is the prime counting function, a function represented with the Greek letter pi (π). Here, the letter pi is not used to represent the number equal to about 3.14159. π(x) is defined as the number of prime numbers less than or equal to x. For example, π(5) = 3, since there are 3 numbers less than or equal to 5 which are prime, and they are 2, 3, and 5. Here are the first 20 values of π(x):
π(1) = 0
π(2) = 1
π(3) = 2
π(4) = 2
π(5) = 3
π(6) = 3
π(7) = 4
π(8) = 4
π(9) = 4
π(10) = 4
π(11) = 5
π(12) = 5
π(13) = 6
π(14) = 6
π(15) = 6
π(16) = 6
π(17) = 7
π(18) = 7
π(19) = 8
π(20) = 8
And here are some larger values of π(x), according to Wikipedia^{[2]}:
π(1000) = 168
π(1,000,000) = 78,498
π(1,000,000,000) = 50,847,534
π(10^12) = 37,607,912,018
π(10^15) = 29,844,570,422,669
π(10^18) = 24,739,954,287,740,860
π(10^21) = 21,127,269,486,018,731,928
π(10^24) = 18,435,599,767,349,200,867,866
Hold on now. How is it possible to calculate π(x) for such large inputs?! Don't you need to manually find every prime number below the value you plug in? As it turns out, you don't need to do it by brute force (i.e. checking every prime); there are methods for calculating π(x) using some mathematical tricks. We won't worry about how exactly you can do that, just that it is possible to calculate such values.
Besides the algorithms for calculating π(x) (which do have a limit for which inputs you can realistically use), there are functions to approximate π(x), which is where the second part of Skewes' problem comes into play.
Approximating π(x)
The distribution of prime numbers has fascinated mathematicians since ancient times, largely since that distribution has no apparent regular patterns. However, over the years mathematicians have proved some useful facts relating to it. The most important one is the prime number theorem, which states that the limit of π(x)/(x/ln(x)) as x approaches infinity is 1. In other words, π(x) and x/ln(x), as x approaches infinity, will be approximately the same value. ln(x), as you probably know, is the function defined as the number a such that e^{x} = a, where e is the important mathematical constant equal to about 2.781828, most commonly defined as the limit of (1+1/x)^{x} as x approaches infinity.
To see that π(x) and x/ln(x) are generally about the same value, let's compare π(x) with x/ln(x) for some small inputs:
With these results, it looks like π(x) will generally be larger than x/ln(x). But these small examples aren't nearly enough to make good conclusions, so let's get some larger examples:
Notice that x/ln(x) always seems to lag a bit behind π(x), though the ratio between x/ln(x) and π(x) becomes closer to 1 as x gets larger. It turns out that when x ≥ 11, x/ln(x) will always be a little smaller than π(x).
There are better approximations to π(x) than x/ln(x). The approximation we're focusing on here is the logarithmic integral function, noted as li(x). If you know your calculus, think of the function as the integral:
x
/\

 1/ln(t) dt

\/
0
Otherwise, just think of li(x) as the area between the graph of y = 1/ln(x) from 0 to x and the xaxis. The graph of y=1/ln(x) looks like this:
and the graph of li(x) looks like this:
The reason why at some points the area is negative is because area of a function below the xaxis is considered negative, while area above it is considered positive.
The logarithmic integral function is a much better approximation for π(x) than x/ln(x) is. Observe some values of each of those functions for large values of x:
(the values for li(x) in the table came from source [3]
li(x), as you can see, is quite a good approximation! Take a look at π(x) (black), li(x) (red), and x/ln(x) (blue) graphed together, and see how these functions stack up:
A(graph taken from [4] but modified)
However, li(x) seems to always be a little bigger then π(x). For quite a while, most mathematicians thought this is always the case no matter how big x gets. However, if you think about the erratic behavior of the distribution of prime numbers, making such a claim seems a lot like jumping to conclusions. Therefore, it became an intriguing mission to see if there is a point where li(x) is less than π(x).
In 1914, John Edensor Littlewood proved that there does indeed exist a number x such that li(x) is less than π(x)^{[5A]}—in fact, there exists an infinite number of x's where li(x) is less than π(x). Since there exists such a number x, that means that there must be a smallest value of x such that li(x) is less than π(x). Finding that smallest x is where Stanley Skewes comes into play.
Skewes' Proof
In 1933, Stanley Skewes attempted to solve the problem of the smallest number x where π(x) > li(x). He did not fully solve the problem, but he did give an upperbound to the solution, i.e. a number that the solution can be shown to be no greater than. He arrived at the upperbound:
a number which can be approximated in decimal as:
and it has an 8.85decilliondigit number of digits. I'll say more on its size later; for now, it's worth noting that the number surpasses a googolplex but falls short of a googolduplex, and it easily broke the record of the largest number used in serious mathematics.
However, there's a little twist: Skewes solved the problem assuming the Riemann hypothesis to be true. What exactly is the Riemann hypothesis? It is a hypothesis that's somewhat more complicated than the other things we've learned about so far, and that's what we'll examine in detail next.
The Riemann Hypothesis
The Riemann hypothesis is one of the most important conjectures in mathematics. It can be thought of as a hypothesis so notable that it has an indirect relevance in the field of large numbers, as it makes the difference between Skewes' number and the second Skewes' number, which are two of the "classic" large numbers alongside the googol, the googolplex, Graham's number, and a few others. The hypothesis was proposed by Bernhard Riemann in 1859, and to this day it remains yet to be found out whether it is true or not. According to Wikipedia,^{[6]} the hypothesis is the subject of several books, and finding whether or not it's true is part of both Hilbert's 23 problems and the Millennium prize problems. Suffice to say, the Riemann hypothesis is a pretty big deal.
But what exactly is the Riemann hypothesis? It has to do with a wellknown function called the Riemann zeta function, a function denoted with the Greek letter zeta (ζ). It's a function that can input any real or complex number. A complex number is any number of the form a+bi, where a and b are real numbers, and i is the imaginary unit defined as the square root of 1. In a complex number, a is the real part and b is the imaginary part, and real numbers are a subset of the complex numbers where the imaginary part is 1.
Here is how the Riemann zeta function is defined. When the real part of the number inputted is greater than 1, ζ(x) is defined as the infinite sum, 1/1^{x} + 1/2^{x} + 1/3^{x} ..., which can be noted with the sum notation as:
∞
Σ 1/n^{x}^{}n=1
It can be extended to other complex numbers using analytic continuation of that series.
The zeta function is known for having fascinating behavior. Here are some values of the function:
ζ(1) = 1/1^{1} + 1/2^{1} + 1/3^{1} .... = 1+2+3+4+5 ... Although this sum diverges to infinity with the usual calculations, using infinite series it can be derived to be equal to 1/12. We won't go into detail on how exactly that paradoxical sum can be derived, but this serves as an example of an unusual value of the function.
ζ(0) = 1/1^{0} + 1/2^{0} + 1/3^{0} + 1/4^{0} ... = 1/1 + 1/1 + 1/1 ... = 1+1+1+1 ..., which can be strangely derived to equal 1/2. Once again we won't go into detail about this strange sum, other than that it is strange. N.H. Abel famously said about such series that "divergent series are the invention of the devil". I can't say I disagree with him.
ζ(1) = 1/1^{1} + 1/2^{1} + 1/3^{1} + 1/4^{1} ... = 1 + 1/2 + 1/3 + 1/4 ... This sum, like the previous two, diverges to infinity, which may be surprising. Since the terms of this sum get smaller and smaller, you may expect it to converge to a finite number, but in fact it once again reaches infinity. The difference between the first x terms of this sum and ln(x), as x approaches infinity, is a mathematical constant known as the EulerMascheroni constant, denoted with the Greek letter gamma (γ), which is equal to about 0.57721566. Even stranger, unlike ζ(1) and ζ(0), a finite value for this sum cannot be derived with infinite series. However, if it is interpreted as a limit:
then that limit takes on the value of the EulerMascheroni constant, which is just as strange.
ζ(2) = 1/1^{2} + 1/2^{2} + 1/3^{2} ... ≈ 1.64493007.... = π^{2}/6. This sum is notable because it gives an example of pi showing up in a place you wouldn't expect it to. The sum was proven to be equal to π^{2}/6 by Euler in 1734.
ζ(3) = 1/1^{3} + 1/2^{3} + 1/3^{3} ... ≈ 1.2020561... This number is known as Apéry's constant, and unlike ζ(2) it has no known compact expression like π^{2}/6 for ζ(2). However there are many other infinite series that sum up to this constant.
ζ(4) = 1/1^{4} + 1/2^{4} + 1/3^{4} ... ≈ 1.082323... = π^{4}/90. Like ζ(2), this is an example of pi cropping up in a place you don't expect it to, and it has a nice compact expression.
However, plugging in real numbers into the zeta function only scratches the surface of what the zeta function has to offer. If you plug complex numbers into the function, things get wild, and that's where the Riemann hypothesis comes into play. What the hypothesis asks is best explained visually. To start off, here's a graph of the zeta function with any real or complex inputs:
This graph is a somewhat unusual graph. Here, the x and y axes represent the inputs of the function, while the colors represent the outputs. The xcoordinate of a point represents the real part of the input number, while the ycoordinate represents the imaginary axis of the number. For example, the point (1, 1) which is one unit up and one unit right of the origin would represent 1+i, whose real part and imaginary parts are both 1.
Now the colors are a bit more complicated. The coloring is based on the complex plane, which is the 2dimensional coordinate plane mapping each complex number to a point on that plane: the xcoordinate is the real part and the ycoordinate is the imaginary part. If the imaginary part is 0, then the number in the complex plane is a real number. The coloring is defined based on the output value's position in the complex plane: the darker a color is, the closer it is to 0, and light colors denote values far from 0. Black is zero, and white is infinity whether it be positive, negative, or whatever else (such as the limit of a*2+a*i when a approaches infinity). The hue of the color denotes which angle the number is from 0 in the complex plane. Red denotes positive real numbers (such as 1), cyan denotes negative real numbers (such as 1), yellowgreen denotes positive multiples of i (such as i), and violet denotes negative multiples of i (such as i). And that's what all those colors mean.
The Riemann hypothesis focuses on what the zeros of the function are: points where ζ(x) = 0. In this graph, the zeros are point of pure black. The dark zone in the middle of the graph isn't quite black, but rather represents numbers quite close to 0. The exceptions are negative multiples of two—those are the zeta function's trivial zeros. The only known nontrivial zeros of the zeta function are black points in this graph which all line up on a vertical line in the middle of the graph. Those zeros in this graph are all complex numbers with the real part equal to 1/2. But are there any nontrivial zeros of the zeta function that do not have the real part equal to 1/2? That's where we learn what the Riemann hypothesis is—it says that the only nontrivial zeros of the Riemann zeta function are where the real part of the input is equal to 1/2. Though this seems very likely to be true, nobody has been able to prove for sure whether it's true. To disprove this hypothesis, all you'd need to do is find one nontrivial zero of the function that doesn't meet these criteria. But proving it true is an ordeal that has stumped mathematicians to this day.
Now, how exactly would the Riemann hypothesis relate to Skewes' number? As it turns out, many of the formulas for the prime counting function make use of the zeros of the zeta function. The reasons for this relation are rather complicated, but they relate to this curious alternate formula for the zeta function, discovered by Euler:
ζ(x) = 1/(12^{x}) * 1/(13^{x}) * 1/(15^{x}) * 1/(17^{x}) * 1/(111^{x}) ...
With the notation for infinite products this can be written as:
∞
Π 1/(1p^{x})
p prime
This alternate formula for the zeta function, as you can see, uses all the prime numbers in its definition. It's a big reason why the Riemann hypothesis relates the distribution of prime numbers. Hopefully, you now have some idea why the Riemann hypothesis is notable to Skewes' problem.
The Second Skewes' Number
In 1955, Stanley Skewes revisited the problem that led to his famous number—the first point where π(x) is first greater than li(x)—but this time he did not assume the Riemann hypothesis to hold true. This constraint made his problem more difficult, and he arrived at an upperbound bigger than Skewes' number, equal to:
This number has come to be known as the second Skewes' number. It's bigger than a googolduplex, and it has roughly a 3.3*10^{963}digitnumber of digits. In decimal, it can be approximated as:
The number dethroned the original Skewes' number as the largest number used in serious mathematics. The two Skewes' numbers became honored as the largest numbers used in mathematics until the much larger Graham's Number came along in 1977, and in subsequent years Graham's number has itself lost that title, but let's not get ahead of ourselves.
Now that we've reviewed the two Skewes' numbers, we're now ready to talk about how big they are.
How big is Skewes' number?
To get an idea of how big either Skewes' number is, we must first convert to decimal. Doing this isn't as hard as you may think. You can make use of the common logarithm log(x), defined as the number a with the property that 10^{a} = x. We start with the expression for the first Skewes' number:
e^{e^e^79}
Then, by converting e to 10^{log(e)}, since by definition 10^{log(x)} = x, we can change this expresion to:
(10^{log(e)})^{e^e^79}
To simplify the expression we use one of the laws of exponents, (a^{b})^{c} = a^{b*c}, and so we obtain:
10^{log(e)*e^e^79}
Next we convert another e to 10^{log(e)}:
10^{log(e)*(10^log(e))^e^79}
and use the laws of exponents again:
10^{log(e)*10^(log(e)*e^79)}
To get rid of the leftmost log(e), we need to convert it to 10^{log(log(e)):}
10^{10^log(log(e))*10^(log(e)*e^79)}
Then we use another law of exponents, this time a^{b}*a^{c} = a^{b+c}:
10^{10^(l}^{og(log(e))+log(e)*e^79}^{)}
The bold expression, log(log(e))+log(e)*e^{79} is an expression that evaluates to a small enough number for a scientific calculator to handle. Plugging that expression into the Keisan online calculator^{[3]} with 50 digits of accuracy gives us the result:
8,852,142,197,543,270,606,106,100,452,735,038.5506795094512173
Thus, we obtain the estimate for Skewes' number:
10^{10^8,852,142,197,543,270,606,106,100,452,735,038.5506795094512173}
which can be compactly approximated as:
10^{10^8.852142197*10^33}
Just how big is that number?! It's certainly hard to give an analogy of its size. Let's consider just the difficulty of writing Skewes' number, or more specifically, the integer part of Skewes' number. The integer part of any real number is just the number with everything after the decimal part cut off. For example, the integer part of 441.93 is 441, and the integer part of 3.16457 is 3.
To find the number of digits in a large number, just take its decimal logarithm, add 1, and round it up to the nearest integer. Finding the logarithm of a large number when written in decimal exponent form is easy: just chop off the bottom 10. In our case the base10 logarithm of Skewes' number is about:
10^{8.852142197*10^33}
And adding 1 makes basically no difference to the number anyway, ditto for rounding up or down. So there we have the number of digits in Skewes' number, and now we can consider how hard it is to write it out.
First off, imagine you could fit a digit into each Planck volume in the observable universe (from here on out I'll just say "universe" instead of "observable universe"). Not that you could actually do that, so this is an EXTREMELY generous estimate of how many digits could fit in the universe. Fitting a digit in each Planck volume would give you:
3*10^{185} digits — not even close!
So we need to take things up a notch. How about we make a secondorder universe, which is a sphere of as many universes as there are Planck volumes in the universe? Another way to put this is shrinking the universe down to the size of a Planck volume, and fitting as many of those shrunkendown universes as you can into our universe, then expanding all those universes to normal size—that's pretty insane alright! And fitting a digit into each Planck volume in the secondorder universe gives you ...
9*10^{370} digits — still not even close!!
Well, let's then make a thirdorder universe, which would be like shrinking down our secondorder universe to the size of a Planck volume, and then filling us our universe with all those shrunkendown secondorder universes. That would be equal to the number of Planck volumes in the universe cubed, which gives you:
2.7*10^{556} digits — >:(
How about we make a sphere of thirdorder universes with as many spheres as there are Planck volumes in the universe, which is a fourth order universe? Or a sphere of the same amount of fourthorder universes, which is a fifthorder universe? Or a sixth, seventh, eighth, ninth, or tenthorder universe? Wait, how about a millionthorder universe, how many digits will that give us? No conventional calculator could calculate that number of digits, so we'll need to use Robert Munafo's online calculator called Hypercalc.^{[8]} Hypercalc can handle numbers far larger than any ordinary calculator can. It can handle any number we'll encounter in section 1 of this site, and many of the numbers we'll encounter in section 2! Raising 3*10^{185} to the millionth power, which is how many digits the millionth order universe would give us, equals:
1.8*10^{185,477,421} digits
This number has 185 million digits, but the number of digits in Skewes' number has 8.85 decillion digits. In terms of more familiar illions, that is 8.85 billion trillion trillion digits!! So we're still far from being able to write Skewes' number!
How about a trillionth order universe? With Hypercalc would give us about:
10^{1.85477*10^14} digits
We're definitely getting somewhere, but this still doesn't cut it. Even by multiplying this number by itself a quintillion times you wouldn't get to the number of digits in Skewes' number! So what order universe would we need to write Skewes' number?
We would need a 48 nonillionth order universe to make it happen. A nonillion is a million trillion trillions! It's near impossible to even comprehend that, but we can try:
Imagine starting with the plain old observable universe, and scaling it up by a factor of 3*10^{185} (the number of Planck volumes in the universe) every second. In 5 seconds we'd already reach past 10^{1000} digits, in a day we reach about 10^{15,951} digits, and in a year we reach a 31millionth order universe which could fit about 10^{5,800,000,000} digits! But we wouldn't reach a universe vast enough to write Skewes' number until 1.52 septillion years! That's about a quintillion times the age of the universe. To get an idea of how long of a time that is, consider that even a quintillion seconds is twice the age of the universe! Could things get any more insane?!
As it turns out, we've only gotten to the difficulty of writing Skewes' number! We haven't even gone close to getting a picture of how big Skewes' number itself is—we've only been blown away by its number of digits!!
To get a feel of how big Skewes' number really is, we need to take things up a notch. First off, let's get back to the 48nonillionthorder universe that could fit the digits of the integer part of Skewes' number. Let's call that universe a Skewes' universe. But now, let's not just have as many of those universes as there are Planck volumes in the universe. That's only going one extra step in universe orders from the 48 nonillion steps we took to get to the Skewes' universe. Instead, let's make as many Skewes' universes as there are Planck volumes in the Skewes' universe! Call that a secondorder Skewes' universe. Likewise, as many secondorder Skewes' universes as there are Planck volumes in the Skewes' universe would make a thirdorder Skewes' universe. Continue with a fourthorder, fifth, sixth, seventh ... tenth, 100th, 1000th ... millionth, billionth, trillionth ... decillionth ... centillionth ... it takes about as many steps up as there are Planck volumes in the Skewes' universe to get to a universe with Skewes' number of Planck volumes!
The Curious Case of Skewes' Approxima
It's so easy to get utterly completely blown away by numbers this big ... and yet, this isn't the end of the insanity. Consider the decimal representation of Skewes' number:
10^{10^8.852142197*10^33}
This can also be approximated as:
10^{10^10^33.94704}
Now consider a common approximation for Skewes' number, which Sbiis Saibian calls "Skewes' Approxima":^{[5}^{B}^{]}
10^{10^10^34}
This looks like quite a good approximation, right? After all, 34 is a very good approximation for 33.94704. You can see how close they are through this image:
If you look very closely at the rectangle on the right, you can tell that these lines aren't exactly the same length. Still, they're very close: if you plan to build a bridge that's 34 meters long and it ends up 33.94704 meters instead, that really wouldn't make much of a difference. Similarly, if you plan on making some calculations based on a city of 33,947 people, 34,000 is a reasonable approximation because cities' sizes tend to be in flux anyway.
But as we'll see, exponents make a BIG difference. Let's compare 10^{33.94704 }to 10^{34}:
Now this is a much more noticeable difference. 10^{33.94704} is about 8.852*10^{33}, or 8.852 decillion, whereas 10^{34} is ten decillion. Imagine devising some Rubik's Cubeesque puzzle and finding out that it has 8.852 decillion possible arrangements. You wouldn't want to approximate it as 10 decillion, would you? That's the difference an exponent makes... but it's nothing compared to the difference that two exponents make.
Comparing 10^{10^33.94704 }to 10^{1}^{0^34} is where things really start to mess with your head. Look at this image:
The decimal expansion of 10^{10^33.94704 }is 10^{33.94704} digits long; it's an 8.852 decilliondigit number. 10^{10^34}, on the other hand, is a 10 decilliondigit number; that's 1.148 decillion more digits. This does not mean that 10^{10^34} is 1.148 decillion times greater than 10^{10^33.94704}. When you add two zeros to one followed by a big wall of zeros, that doesn't just multiply it by two, but multiples it by 100! Similarly, 10^{10^34} is a 1.148 decilliondigit number of times greater than 10^{10^33.94704}. You can visualize this through the image above: imagine taking a cube of 10^{10^33.94704 }particles, then take a thousand of those cubes, a thousand of those cubes, and repeat this process 382.68 nonillion times. Only then will you get to 10^{10^34}! Even this approximation itself uses approximations; if you overshoot the number of times to iterate, you'll get a number at least a thousand times bigger than 10^{10^34}. And that's no good, is it? 33 too many iterations will get you about a googol times bigger than 10^{10^34}. Does your head hurt yet?!
Well, if it makes you feel any better, squaring 10^{10^33.94704}—which is to say, multiplying it by itself—will get you past 10^{10^34}. Specifically, it will get you up to 10^{10^34}^{.24807}, a 17.704 decilliondigit number. This is sure to make the difference seem modest relative to our final comparison: Skewes' number and Skewes' Approxima.
To calculate how much bigger 10^{10^10^34} is than Skewes' number, we'll represent it using a somewhat closer approximation, namely 10^{10^33.94704}.
Remember the Skewes' universes we went over when exploring the size of Skewes' number? Let's redefine a Skewes' universe as not just a universe big enough to fit the digits of the integer part of Skewes' number, but as a universe that is made of a Skewes' number of Planck volumes! Now, let's not just have as many Skewes' universes as there are Planck volumes in the universe, or as many as there are digits in Skewes' number. Instead, let's have Skewes' number of Skewes' universes (under our new definition)! Call that a secondorder Skewes' universe. Then of course, Skewes' number of secondorder Skewes' universes would be a thirdorder Skewes' universe. Then we can have a fourthorder, fifthorder, sixthorder, seventh, eighth, ninth, tenth ... 100th ... 1000th ... millionthorder universes ... and so on ...
So then does it takes 48 nonillion scalings up from Skewes' universe to reach 10^{10^10^34} Planck volumes? Nope, it's much much worse. With 48 nonillion scalings up we would not even reach a universe with 10^{10^10^33.94705} Planck volumes (in contrast with 10^{10^10^33.94704}, which is approximately Skewes' number)! Hold on now, how can that be right?! Let's take a closer look at what's going on:
Each time we scale up our nthorder Skewes' universe to a (n+1)th order Skewes' universe, we multiply the number of Planck volumes by Skewes' number. A nthorder Skewes' universe would therefore have S^{n} particles, where S is Skewes' number. Let's see what happens with a 48nonillionth order Skewes' universe:
(10^{10^10^33.94704})^{4.8*10^31} — that's Skewes' number raised to the 48nonillionth power.
= 10^{(10^10^33.94704)*(4.8*10^31)} — just using the laws of exponents here.
≈ 10^{(10^10^33.94704)*(10^31.6812)}
= 10^{10^(10^33.94704+31.6812)}
Now let's consider 10^{33.94704} + 31.6812. This is itself much smaller than 10^{33.94705} which is about 10^{33.94704} + 3*10^{28}. All this is completely counterintuitive—and yet we haven't gotten close to the sphere of 10^{10^10^34} Planck volumes!
Let's continue the nthorder Skewes' universe idea now. After the 48nonillionth order, let's continue through the orders, with a decillionthorder Skewes' universe, a vigintillionthorder Skewes' universe ... googolth, centillionth, 10^{1000}th, 10^{10,000}th ... 10^{1,000,000}th ... 10^{10^12}th order ... and so on ... when we reach a 10^{1.15*10^33}th order Skewes' universe, that would be made up of about 10^{10^10^34} Planck volumes. Absolutely mindblowing!!! I think an image won't do these insane levels of scaling up justice; instead, I'll let you imagine it in your head, if you haven't already.
As you can see, getting to the number of digits in Skewes' number is already utter insanity, and then getting to Skewes' number itself is far more insane. Unbelievably, it also takes a massive climb from Skewes' number to what looks like a good approximation for Skewes' number!
As you can see, Skewes' number is an utterly unfathomable number, with some pretty mindbending properties. However, it's nothing compared to the even more monstrous second Skewes' number, which is the next number for us to examine.
How big is the second Skewes' number?
The second Skewes' number, which is less known than the original Skewes' number, is equal to:
e^{e^e^e^7.705}
Converting that to decimalgives us:
(10^{log(e)})^{e^e^e^7.705}
= 10^{log(e)*e^e^e^7.705}
= 10^{(10^log(log(e)))*(10^log(e))^e^e^7.705}
= 10^{10^log(log(e))*10^(log(e)*e^e^7.705)}
= 10^{10^(log(log(e))+log(e)*e^e^7.705)}
Now, what to do with that log(log(e))? Adding it to log(e)*e^{e7.705} will have almost no effect on it, since log(e)*e^{e7.705} is so big. Therefore we can chop it off and still be able to approximate the second Skewes' number. Thus we continue with:
10^{10^(log(e)*e^e^7.705)}
= 10^{10^10^(log(log(e))*(10^log(e))^e^7.705)}
= 10^{10^(10^log(log(e))*10^(log(e)*e^7.705))}
= 10^{10^10^(log(log(e))+log(e)*e^7.705)}
Now log(log(e))+log(e)*e^{7.705} is actually small enough to calculate with an ordinary calculator. Plugging it into the Keisan online calculator gives us 963.5185064675691205137. Therefore the second Skewes' number is about:
10^{10^10^963.5185064675691205137}
This number is not only much much larger than a googolplex, it's also larger than a googolduplex which is:
10^{10^10^100}
Hopefully you are starting to get the idea that the second Skewes' number is HUGE ... but how much larger is it than Skewes' number? Let's return to the idea of Skewes' universes. First there is the Skewes' universe of Skewes' number of Planck volumes, then the secondorder Skewes' universe made of Skewes' number of Skewes' universes, then thirdorder, fourthorder, and so on. But let's not stop at the 10^{1.15*10^33}th order, but go to however many scalings up it takes us to get to the second Skewes' number.
Even after a googolplex scalings, we wouldn't reach the second Skewes' number. It would take about 10^{10^963.518 }scalings to reach it. To get a feel of how many scalings that is, let's use a time analogy.
Imagine a universe much much more vast than ours, whose volume in Planck volumes is, at any moment, denoted by a number N. In the beginning of that universe, its volume in Planck volumes is Skewes' number. Now imagine that every Planck time from there on out the universe's volume multiplies by Skewes' number. This means that one Planck time into the universe's time it would be equivalent to a secondorder Skewes' universe, and in the next Planck time after that it would be a thirdorder Skewes' universe, and so on. After just one second, it would be equivalent to roughly a 10^{43}th order Skewes' universe!
Now imagine going forward in time from this beginning, up to 13.7 billion years (the age of our universe). At that point that universe would be as old as our universe is. Call that time T1. But now, imagine zooming waaaaaaaay past T1, to the point where the span between the Big Bang and T1 is just one span amongst as many spans as Planck times since the Big Bang! In other words, we've now dwarfed the age of the universe by a factor of the age in Planck times! Call that point in time T2, the point where the universe is that old.
So at this mysterious strange time T2, will that universe be made of at least 10^{10^10^34} Planck volumes? Nope, not even close. The age of the universe in Planck times is about 8.03*10^{60},^{[9A]} and multiplying that by itself gives you an age of 6.45*10^{121} Planck times for the point in time, T2. This means that that universe at that point in time is equivalent to a 6.45*10^{121}th order Skewes' universe, aka Skewes' number to the power of 6.45*10^{121} Planck volumes. As with the 48nonillionthorder Skewes' sphere, the number is much less than something like 10^{10^10^33.94705}. Simply observe:
(10^{10^10^33.94704})^{6.45*10^121}
= 10^{(10^10^33.94704)*(6.45*10^121)}
~ 10^{(10^10^33.94704)*(10^121.81)}
= 10^{10^(10^33.94704+121.81)}
So at what point in time would we reach 10^{10^10^34} Planck volumes? We'll need to scale things up more. Imagine scaling the age of the T2 universe up by a factor of the age of the T1 universe in Planck times. The point in time where we've reached that point, T3, is about 5.18*10^{182} Planck times after the start of the universe. Then we can have T4, which is scaling up the time from the start of the universe to T3 by a factor of 8.03*10^{60}, then T5, T6, T7 ... T10 ... T100 ... T1000 ... ... T1,000,000 ... ... T10^{12}. Turns out that we finally reach 10^{101034} Planck volumes at T1.9*10^{31}, or after 19 nonillion scalings up. Now think about that for a minute. You may think that a nonillion scalings up is not entirely unimaginable, but think again. If you do one scaling up per second starting from the Big Bang of our universe, in the present day you'd be under a trillionth of the way there. And guess what ... we aren't even close to the second Skewes' number!
After a googol scalings up (i.e. T10^{100}) we'd get a universe with about 10^{10^(6*10^101)} Planck volumes. That's over a googolduplex. Now a googol is MUCH LARGER than the age of our universe in Planck times. That means that if you go start with T1 at the Big Bang, then do one scaling up (e.g. from T35 to T36) every Planck time, and continued to the present day, you still wouldn't reach a universe even close to large enough to fit a googolduplex Planck volumes, let alone the second Skewes' number of Planck volumes.
At T5.15*10^{961}, we FINALLY REACH a second Skewes' number of Planck volumes. Now 5.15*10^{961} (that's more than a centillion centillion centillions) is a number that is itself far far far beyond human comprehension. And that's the number of scalings up of a device that was itself used to get an idea of all the scalings up, needed to reach the second Skewes' number!
Is that mindblowing? Good, hopefully you realize now that the Skewes' numbers are exotically huge. They exhibit utterly counterintuitive behavior we don't see with numbers we use in daytoday life at all; and yet, they can be defined quite easily using standard mathematics! These large numbers are nothing like the puny millions, billions, and trillions we encounter in daytoday life—they're utterly exotic numbers that they can't even compare to!!!
The funny thing is that the most commonly given example of the size of numbers this big is that all of their digits wouldn't fit in the universe; however, this doesn't come near capturing a proper idea of their size, since even a number like 10^{10^200} (which is much smaller than either Skewes' number) would be GUARANTEED not to fit, even with fitting each digit in a Planck volume! To actually get an idea of how amazingly big these numbers must be, we need to say goodbye to even considering realistically writing all their digits, and take much more radical methods! This is just what I did, a rare example of an attempt to get a true feel of how big Skewes' numbers really are (Sbiis Saibian's article on Skewes' numbers^{[5A]} is another example). It's actually a little annoying that it is so popular to simply say that any number bigger than a googolplex has so many digits that they wouldn't fit in the universe; in a sense, this represents shying away from the enormity of these monsters, when it isn't too hard to use something much much more mindblowing.
But before we get carried away into a discussion of people not even wanting to get a full picture of the magnitude of these giants, let's switch gears and go into the further history of Skewes' numbers after Skewes' proofs.
Better bounds for Skewes' problem
After Skewes' proof in 1955, better upperbounds for the solution to Skewes' problem were found. According to Wikipedia's article on Skewes' number, the first proof of a better upperbound was by R. Sherman Lehman in 1966, and it proved that the first crossing was below 1.65*10^{1165}.^{[10]} Now think about the drastic decrease in upperbound from either Skewes' number. The digits of Skewes' number could not be stored in the observable universe even if you dwarf it by a factor of its volume in Planck volumes a nonillion times, and you couldn't store the digits of the second Skewes' number even if your dwarf it by the same factor a centillion centillion centillion times. However, this number, the next best upperbound found, could easily have all its digits stored on a sheet of paper. This represents a drastic improvement in the upperbound from the Skewes' numbers. This better upperbound is not even close to being an "exotic number" like the Skewes' numbers are! And yet, it's still an insane number in comparison to our everyday millions, billions, and trillions, or even numbers like a googol or a centillion.
In addition, Lehman's proof took away Skewes' numbers' status as practical upperbounds, and from there on out they became recognized instead as footnotes or mathematical curiosities, as former upperbounds that ended up surprisingly huge. Nonetheless Skewes' numbers retain their recognition as "legitimate" numbers due to their use other than just for the sake of largeness. They're important pieces of the history of both the study of prime numbers and of googology, just like how we still recognize scientists who weren't quite right about what they studied but still paved the way for our current understanding.
For the next five years, Skewes' numbers were still recognized as the largest numbers in mathematics. But they were dethroned from that title in 1971 by a far far larger number used as an upperbound in an unrelated problem in Ramsey theory by mathematicians Ronald Graham and Bruce Rothschild. That number is a whole different kind of number from Skewes' number, as it can't be represented with mainstream mathematical notation and needs a whole new notation to express! An alternate, larger version of that number was devised in 1977 by mathematics writer Martin Gardner as an easiertoexplain version of Graham and Rothschild's number, and that number is today known as Graham's number. To this day, Graham's number is still regularly honored as the largest number used in serious mathematics, although numbers that are themselves FAR LARGER than Graham's number have since then taken that title, such as TREE(3) and SCG(13), which aren't upperbounds at all but rather extremely large solutions to problems in mathematics! But let's not get ahead of ourselves. We'll discuss Graham's number in detail in section 2, and for now let's get back to the further history of Skewes' problem.
H.J.J. te Riele in 1987 proved a still better upperbound to the problem, equal to e^{e^27/4}, or about 8.17*10^{370}.^{[11][9B]} This bound is equal to roughly eighty thousand vigintillion centillion. It's quite a lot smaller than Lehman's bound from earlier, having about a third as many digits. Then in 2000, the bound was improved to 1.39822*10^{316} by Bays and Hudson.^{[12]} It may be tricky to remember that this bound is about a septendecillion (that's a quintillion quintillion quintillions) times smaller than Riele's. After 2000, a few small improvements were made to the upperbound. The most recent such improvement was in 2010 by Stefanie Zegowitz,^{[13]} reducing the upperbound to e^{727.951332973}, which is about 1.39716*10^{316}. This is only 0.07% smaller than the previous bound 1.39822*10^{316}. Robert Munafo incorrectly claims that this value is the solution to Skewes' problem.^{[9B]}
To this day, the solution to Skewes' problem remains unknown and an open question. However, it seems likely that the solution is around 1.39716*10^{316}, considering that nobody has been able to improve the bound very much from the slightly larger 1.39822*10^{316}. But we still don't know the solution for sure.
In addition, there is a current lower bound for the solution to Skewes' problem, a number that the solution is known to be no smaller than. If a number X is a lower bound for the solution to a problem, that means that we know that the solution must be a number greater than X. For Skewes' problem there must be some lower bound: we can check each individual number up to some point for whether or not π(n) > li(n). However, after a certain point doing that checking will become difficult. Therefore, proving a lower bound for this problem is harder than finding an upper bound! To find an upper bound, you just need to prove that some number n less than a number x has the property π(n) > li(n), but to prove a lower bound, you need to prove that π(n) > li(n) does not hold true for all n less than your lowerbound x. Since 2015, the best known lowerbound for Skewes' problem has been 10^{19}, or ten quintillion, proven by Jan Büthe.^{[11][14]} While I can't really comment on how he proved that lower bound, that allows us to say that, we currently know that:
10^{14} < N < 1.39716*10^{316}
where N is the smallest positive integer such that π(N) > li(N). As you can see, Skewes' problem is still yet to be solved, though we're certainly closer to an answer than we were in Skewes' time.
Conclusion
Today, Skewes' numbers continue to have a wellearned spot among the classic large numbers, and are given an honorable mention in the large number discussion. However, I doubt you've heard their real size explained to such detail as I did! These numbers are pretty insane and headspinning numbers, and they continue to be noted as mathematical curiosities, and as examples of really really big numbers. But as it turns out, these numbers aren't even close to being the biggest numbers ever, not even the biggest numbers in the popular imagination. They are just barely starting the endless numberscape of googology!
In any case, the Skewes' numbers are among the most famous large numbers these days, and therefore are worth discussing in detail in an article like this. Up next, we'll look at extensions to the illions that zoom way past both of the Skewes' numbers!