Circumference to radius ratio

In Texas kids learn about the number “pi” in 6th grade. I do not know how it is in Kansas. But we are not in Kansas anyway, so we can learn it now.

We are not in Texas, and we are not in Kansas. We are on the unit disk with hyperbolic geometry defined by the transitive action of SU(1,1) – as it was discussed in the recent series of posts. We know the line element:

(1)   \begin{equation*}ds^2=\frac{4(dx^2+dy^2)}{(1-x^2-y^2)^2}.\end{equation*}

We can now ask this pressing question: what is the ratio of the circumference of a circle to its radius?

Any point of the disk is as good as any other, but in the standard coordinates (x,y) that we are using the origin of the disk is most convenient. Let us introduce polar coordinates

(2)   \begin{eqnarray*} x&=&\rho\cos\phi,\\ y&=&\rho\sin\phi. \end{eqnarray*}

We need to express our line element in terms of (\rho,\phi). We calculate

    \[ dx =d\rho\cos \phi-\rho\sin\phi d\phi,\]

    \[ dy =d\rho \sin\phi+\rho\cos\phi d\phi,\]

    \[dx^2= d\rho^2\cos^2\phi-2\rho\sin\phi\cos\phi d\rho d\phi+\rho^2\sin^2\phi d\phi^2,\]

    \[dy^2= d\rho^2\sin^2\phi+2\rho\sin\phi\cos\phi d\rho d\phi+\rho^2\cos^2\phi d\phi^2,\]

    \[dx^2+dy^2=d\rho^2+\rho^2 d\phi^2,\]

    \[ds^2=4\frac{d\rho^2+\rho^2d\phi^2}{(1-\rho^2)^2}.\]

Consider now a circle with \rho=\rho_0 and \phi running from 0 to 2\pi. What is the length C of its circumference? We have to integrate ds. Going around circumference \rho is constant, so d\rho=0. Therefore

    \[C=\int_0^{2\pi}ds=\int_0^{2\pi}\frac{2\rho_0 d\phi}{1-\rho_0^2}=\frac{4\pi \rho_0}{1-\rho_0^2}.\]

What is the length r of its radius? Along the radius \phi is constant, so d\phi=0, so

    \[r=\int_0^{\rho_0} ds=\int_0^{\rho_0}\frac{2d\rho}{1-\rho^2}=2\mathrm{arctanh}\rho_0.\]

From the last equation we get

    \[\rho_0=\tanh\frac{r}{2}.\]

Substituting into the formula for C:

    \[C=4\pi\frac{\sinh(r/2)/\cosh(r/2)}{1-\sinh^2(r/2)/\cosh^2(r/2)}=2\pi\frac{2\sinh(r/2)\cosh(r/2)}{\cosh^2(r/2)-\sinh^2(r/2)}.\]

Thus

(3)   \begin{equation*}C=2\pi \sinh r,\end{equation*}

and

(4)   \begin{equation*}\frac{\mathrm{circumference}}{\mathrm{radius}}=\frac{C}{r}=2\pi\frac{\sinh r}{r}.\end{equation*}

We can find this formula in Wikipedia article Hyperbolic geometry, in the section Circles and disks (putting there the curvature R=-1).

Looking at the graph of the function \sinh r/r

we see that the ratio circumference to radius on the disk is always greater than 2\pi. In fact the ratio is growing faster and faster with increasing r.

Hyperbolic angels and demons

In SU(1,1) preserves angles we have seen this picture:

Fig. 1

These are straight lines on the disk obtained from the interval [-0.5,0.5] on the real axis by the action of a one-parameter family A_1(t) of left shifts from SU(1,1). Since our geometry is, by definition, invariant under such transformations, all these lines are of equal length.

I thought it would be interesting to overlay these lines, obtained numerically, with Escher’s Angels and Demons picture to see how well Escher’s drawing corresponds to calculations. Here is the result:

Escher’s Angels and Demons and hyperbolic geometry

As far as I can see angels and demons do not know how to use computers. They work using intuition and aesthetic feelings. And yet Escher’s drawing corresponds to calculations almost perfectly.

In the last post, Following Einstein: deriving Riemannian metric on the Poincaré disk, \url{http://arkadiusz-jadczyk.eu/blog/2017/04/einstein-deriving-riemannian-metric-poincare-disk/} we have obtained the formula for the distance on the disk:

(1)   \begin{equation*}ds^2=\frac{4(dx^2+dy^2)}{(1-x^2-y^2)^2}.\end{equation*}

We can use it to find, say, the distance of a point (0,y_0) on the vertical axis from the origin. Tha path is paramatrized by (0,y) where y varies from 0 to y_0. Therefore the length of the path
is

    \[ s=\int_0^{y_0} ds= \int_0^{y_0} \frac{2 dy}{1-y^2}=2|\mathrm{arctanh} (y_0)|.\]

So, we have

(2)   \begin{equation*} s(y)=2|\mathrm{arctanh} (y)|.\end{equation*}

For producing Fig. 1 we used the family A_1(t) of SU(1,1) matrices;

    \[A_1(t)=\begin{bmatrix}\cosh(t)&i\sinh(t)\\-i\sinh(t)&\cosh(t).\end{bmatrix}\]

A_1(t) implements fractional linear transformation

    \[A_1(t):z\mapsto \frac{\cosh(t)z-i\sinh(t)}{i\sinh(t)z  +\cosh(t)}.\]

Thus

(3)   \begin{equation*}A_1(t):0\mapsto z(t)=-i\tanh(t).\end{equation*}

The distance from the origin to z(t), along the vertical axis, is therefore

    \[s(t)=2|\mathrm{arctanh}(-\tanh(t))|=2|t|.\]

It follows that the parallel lines in Fig. 1 are all equidistant, if the distance between them is measured along the vertical axis.

That was our first use of the formula for the line element of the non-Euclidean hyperbolic geometry of the Poincaré disk. We will have to play with it a little bit more, so that we will not be afraid that it will bite us.

Following Einstein: deriving Riemannian metric on the Poincaré disk

There is no geometry without metry. That is, without metric. It was Bernhard Riemann who gave the foundations to what we call today “Riemannian geometry”. When Einstein learned Rimennian geometry, he was able to compete with smart mathematicians and lay the foundations of the theory of gravitation – the General Relativity Theory.

“[In 1912] I suddenly realized that Gauss’s theory of surfaces holds the key for unlocking this mystery. I realized that Gauss’s surface coordinates had a profound significance. However, I did not know at that time that Riemann had studied the foundations of geometry in an even more profound way. I suddenly remembered that Gauss’s theory was contained in the geometry course given by Geiser when I was a student… I realized that the foundations of geometry have physical significance. My dear friend the mathematician Grossmann was there when I returned from Prague to Zürich. From him I learned for the first time about Ricci and later about Riemann. So I asked my friend whether my problem could be solved by Riemann’s theory [Pais’s italics], namely, whether the invariants of the line element could completely determine the quantities I had been looking for.”

Albert Einstein, as quoted by Abraham Pais in Subtle is the Lord, Pais’s scientific biography of Einstein.

So, it is time for us to get familiar with this “line element“. Here it is, written with the hand of Albert Einstein himself, probably in 1913, in his Zurich notebook:

From A Peek into Einstein’s Zurich Notebook“” by John D. Norton

We want to find our “line element” on the disk. Equivalently: we want to find the metric tensor G_{\mu\nu}.

That is, when from a point with disk coordinates (x,y) we make a small step, and arrive at the point with coordinates (x+dx,y+dy) what will be the length ds of this step? If the geometry were Euclidean, and our coordinates were Cartesian, rectangular, the answer would be simple: ds^2=dx^2+dy^2 from the Pythagoras formula. But our geometry is non-Euclidean. And we want our step length to be invariant under SU(1,1) transformations, which are non-linear on the disk. What to do?

We do what Einstein did in 1912

Einstein1912

We calculate. By the way: as we can learn from History of General Relativity and Gravitation, on his way to General Theory of Relativity Einstein was making errors

Two fundamental errors led Einstein to reject generally covariant gravitational field equations for over two years as he was developing his general theory of relativity. The first is now well known. It was the presumption that weak, static gravitational fields must be spatially flat and a corresponding assumption about his weak field equations. I conjecture that a second hitherto unrecognized error also defeated Einstein’s efforts: he unwittingly reified his spacetime coordinate systems. The same error, months later, allowed the hole argument to convince Einstein that all generally covariant gravitational field equations would be physically uninteresting.

While doing our calculations we will try to avoid errors, at least those computational. We will listen to Gottfried Leibniz who, as mentioned in the last note, has made this advice:

“It is unworthy of excellent men to lose hours like slaves in the labor of calculation which could be relegated to anyone else if machines were used.”

Again we will use a machine, though a different one than in From SU(1,1) to the Lorentz group.

As we have learned in From SU(1,1) to the Lorentz group, a point (x,y) on the disk determines a point (X,Y,T) on the hyperboloid in three-dimensional space-time

Space-time hyperboloid and the Poincare disk models

When we are making a small step (dx,dy) on the disk, we are making also a small step (dX,dY,dT) on the hyperboloid. Well, I am joking. When our step on the disk is very close to the boundary, the corresponding step on the hyperboloid can be “huge”. What I really mean is: when we make an infinitesimally small step …. The concept of an infinitesimally small step belongs to the calculus, and I will assume that we know how to use the elementary calculus concepts.

In the the previous post we used the formulas expressing (X,Y,T) as functions of (x,y):

(1)   \begin{eqnarray*} X&=&\frac{2x}{1-x^2-y^2},\\ Y&=&\frac{2y}{1-x^2-y^2},\\ T&=&\frac{1+x^2+y^2}{1-x^2-y^2}. \end{eqnarray*}

From these formulas we can calculate (dX,dY,dT) in terms of (dx,dy) according to the following general prescription

(2)   \begin{eqnarray*} dX&=&\frac{\partial X}{\partial x}dx+\frac{\partial X}{\partial y}dy,\\ dY&=&\frac{\partial Y}{\partial x}dx+\frac{\partial Y}{\partial y}dy,\\ dT&=&\frac{\partial T}{\partial x}dx+\frac{\partial T}{\partial y}dy. \end{eqnarray*}

We can do these calculations by hand, but that would be unworthy of excellent men. Machine will do it for us. Once we know (dX,dY,dT) we can calculate the corresponding line element using the formula of special relativity, a generalization of the Pythagoras law from space to space-time. I our case it is now convenient to chose the signature (+1,+1,-1) then

(3)   \begin{equation*}ds^2=dX^2+dY^2-dT^2.\end{equation*}

I let Maxima to do the calculations using this code:

X(x,y):=2*x/(1-x^2-y^2);
Y(x,y):=2*y/(1-x^2-y^2);
T(x,y):=(1+x^2+y^2)/(1-x^2-y^2);
factor(ratsimp(
(diff(X(x,y),x)*dx+diff(X(x,y),y)*dy)^2+
(diff(Y(x,y),x)*dx+diff(Y(x,y),y)*dy)^2-
(diff(T(x,y),x)*dx+diff(T(x,y),y)*dy)^2);

If you want, save this code as a file, say, “ds2.mac”. Install and run wxMaxima, use File/Batch File … to load the code. The answer is produced immediately:

(4)   \begin{equation*}ds^2=\frac{4(dx^2+dy^2)}{(1-x^2-y^2)^2}.\end{equation*}

And that is the famous line element of hyperbolic geometry – see Metric and volume element on the Poincaré disk on Wikipedia.

We will yet have to find out what can we do with it? Lot of things. We can find geodesics (we have already guessed what they are), curvature …. Yes, our disk is curved even if it does not look so.