Showing posts with label Relatividad General. Show all posts
Showing posts with label Relatividad General. Show all posts

SUSY en el LHC, Cuerdas sin SUSY y lo que todo físico debería saber sobre Cuerdas

Hace poco encontré estos vídeos sobre supersimetría y teoría de cuerdas, dos temas íntimamente relacionados:

Supersimetría: el LHC regresa a la búsqueda de nueva física
Este vídeo está dirigido al público en general. La anfitriona entrevista a varios físicos, entre ellos a John Ellis, quien habla acerca del Run II del LHC y en particular sobre la búsqueda de nueva física y supersimetría. La explicación a nivel divulgativo es la usual, pero las entrevistas son interesantes y se presenta también el lado más escéptico de la propuesta de SUSY con Alvaro De Rújula que bromea un poco sobre viajar al pasado.

Entre otras cosas, si no se encuentra SUSY en el Run II del LHC, John Ellis señala que tendrá que limpiar su oficina.
El escritorio de John Ellis @CERN

Teoría de cuerdas en un universo sin supersimetría
Luego del Run I del LHC, los medios y "algunos físicos" comenzaron a decir que SUSY estaba muerto. En realidad, y como se explica en el vídeo anterior, SUSY podría estar en cualquier escala de energía y podría incluso ser inalcanzable por el LHC. De cualquier modo esto mete en problemas a muchos teóricos y en particular a los que se dedican a la construcción de modelos realistas de partículas (tipo modelo estándar) a partir de teoría de cuerdas (fenomenología de cuerdas).

Algunos recursos para abordar este tema conceptualmente son:
What if the LHC doesn't see SUSY? en Physics SE
Why string theory implies supersymmetry en The Reference Frame
Saúl Ramos, investigador del IFUNAM, aborda el tema de Cuerdas sin SUSY (y explica exactamente a qué se refiere esto en el tiempo 27:00 gracias a una pregunta) precisamente motivado por el hecho de que no ha habido evidencia alguna de SUSY en el LHC, entre otros.

El squark stop es la superpartícula asociada al top quark y era una de las primeras que se esperaba en el LHC 

La charla no es divulgativa (seminario de Altas Energías en la UNAM) pero es accesible para quien esté interesado sin ser experto. El audio es relativamente malo pero vale la pena intentar seguirlo (como sea se agradece la iniciativa del camarógrafo) y la acción en específico sobre la construcción de N=0 comienza en el tiempo 21:00.


Finalmente, no recuerdo si durante la presentación se cita el artículo en que está basada la presentación, pero parece que es este trabajo:
Non-supersymmetric heterotic model building
Michael Blaszczyk, Stefan Groot Nibbelink, Orestis Loukas, Saúl Ramos-Sánchez

We investigate orbifold and smooth Calabi-Yau compactifications of the non-supersymmetric heterotic SO(16)×SO(16) string. We focus on such Calabi-Yau backgrounds in order to recycle commonly employed techniques, like index theorems and cohomology theory, to determine both the fermionic and bosonic 4D spectra. We argue that the N=0 theory never leads to tachyons on smooth Calabi-Yaus in the large volume approximation. As twisted tachyons may arise on certain singular orbifolds, we conjecture that such tachyonic states are lifted in the full blow-up. We perform model searches on selected orbifold geometries. In particular, we construct an explicit example of a Standard Model-like theory with three generations and a single Higgs field.

Lo que todo físico debería saber sobre teoría de cuerdas
Al fin una charla de Eddie Witten que todos podemos entender... todos los físicos... todos los que conocen, al menos de lejitos, QFT y obviamente relatividad general.

Las diapositivas de la charla están disponibles aquí a través del sitio del Strings 2015.

A groso modo, en la primera parte, Witten comienza con la que sería una formulación 1-dimensional de gravedad cuántica, en la que el tensor de Riemann es idénticamente nulo y así la curvatura de Ricci y la acción de Einstein-Hilbert, con D campos escalares y llega a una teoría de campo en un espaciotiempo D-dimensional curvo. Luego incluso, con este modelo unidimensional, continúa calculando el propagador de una partícula de un punto a otro e introduce interacciones. Finalmente considera el caso dos dimensional, explica por qué Cuerdas es una teoría de gravedad cuántica (no esperen conocer al gravitón, aunque sólo quedan unas patadillas para describirlo) y discute cómo es que en cuerdas desaparecen las divergencias ultravioleta.


Hay además disponible en YouTube una charla de 2013 del mismo Witten llamada Feynman diagrams in String Theory que parece ser esencialmente la misma, aunque quizá con un poco más de detalle (en la primera diapositiva menciona "I will aim to explain the minimum about string perturbation theory that every quantum physicist should know"). Finalmente en algo con relación únicamente a Witten: How does Edward Witten know so much math?.

http://abstrusegoose.com/332

A cien años del nacimiento del teorema de Noether: su relación con la Relatividad General

Para noviembre de 1915, Albert Einstein había dado a la relatividad general su forma final a través de las ecuaciones que llevan su nombre. Unos días después, el poco menos conocido para la población general pero también brillantísimo David Hilbert, enviaría para publicación las mismas ecuaciones de campo obtenidas independientemente (aunque ya con la mesa conceptual puesta por Einstein) a través del irrazonablemente efectivo principio de mínima acción. Tanto menos conocida todavía es Amalie Emmy Noether, que por el mismo tiempo, e influenciada/motivada por Hilbert y el desarrollo de la relatividad general, tenía listo ya (aunque se publicaría en 1918) el bellísimo teorema que hoy en física se conoce por su nombre.

Imagen sobre Emmy Noether creada por el Perimeter Institute (publicada en su cuenta de Twitter)

Aunque desconozco su nivel de popularidad entre los matemáticos, Emmy Noether ciertamente es una rockstar entre los físicos y aún más entre los teóricos. Su vida es en sí misma bastante interesante, siendo desfavorables gran parte de las circunstancias de la época para una mujer que extrañamente estaba interesada en las matemáticas (concretamente se dice que la segunda mujer en dedicarse a las matemáticas en la historia reciente). Irónicamente, aunque Noether es ampliamente celebrada en la física, su teorema no fue más que un pequeño paseíllo en su interés por la teoría de grupos y el resto de su trabajo se restringió únicamente a las matemáticas.


Aunque el teorema de Noether no es del mismo calibre que el de una teoría (como la propia relatividad general), es un ingrediente esencial en el entendimiento de la realidad física. Aunque a muchas personas la palabra teorema les invoque algo escalofriante, la idea central del teorema de Noether es realmente simple: A cada simetría continua de las leyes de la física corresponde una cantidad conservada.


Usualmente se parafrasea de este modo para no meter la pata en algo esencial (como la palabra continua). Para los dummies como yo, básicamente se trata de que cuando un sistema físico permanece inalterado ante una alteración pequeñita dada, existirá una cantidad que no cambiará con el tiempo. El ejemplo más común es el que se enseña en la secundaria sobre la conservación de la energía (la energía no se crea ni se destruye) que corresponde a que las leyes de la física son las mismas sin importar el paso del tiempo (hay una simetría temporal continua: traslaciones temporales). De manera análoga, la conservación de momento (o cantidad de movimiento) corresponde a una simetría ante traslaciones espaciales y la conservación de momento angular corresponde a una simetría ante rotaciones espaciales. Estos son los casos más sencillos en la mecánica clásica, pero seguido las teorías físicas contienen otros tipos de simetría (a veces más evidentes, a veces menos evidentes) y aparecen nuevas cantidades conservadas: el mejor ejemplo es el de la carga eléctrica, que surge de la simetría de la propia teoría electromagnética. Así pues, el teorema de Noether es de una importancia crucial en la física aunque no se trate de una teoría en particular.

Vale, ya entrando en materia, aunque normalmente el teorema de Noether se introduce en mecánica clásica a los estudiantes de licenciatura, el teorema realmente germinó en Emmy Noether a través de la relatividad general cuando en 1915 fue invitada por David Hilbert y Felix Klein a trabajar a la Universidad de Gotinga. Una famosa anécdota*, que refleja claramente las ideas de género de la época, es la de algunos profesores de humanidades reclamando a la llegada de Noether
¿Qué pensarán nuestros soldados cuando regresen a la Universidad y encuentren que deben aprender a los pies de una mujer?
y Hilbert contestando
No veo que el género del candidato sea un argumento contra su admisión como Privatdozent. Después de todo, somos una Universidad, no un baño público.
A pesar de las objeciones, Emmy llegó a Gotinga, aunque sin recibir un salario (siendo mantenida por su familia, misma situación que tuvo, ya como profesora, desde hacía casi 7 años en su ciudad natal) por algún tiempo y tratada como mera ayudante de Hilbert, justo prácticamente al mismo tiempo en que Einstein fue invitado a dar una serie de charlas sobre sus ideas en relatividad general. Hilbert estaba ampliamente interesado en la teoría de Einstein y en la física en general, por lo que obtuvo por su cuenta (quizá poco antes o poco después de Einstein) las ecuaciones de campo de la teoría y las publicó apenas unos días después de que lo hiciera Einstein. Como sea, aunque como menciona Einstein en su artículo, la relatividad general estaba completa como construcción lógica, aún había un detalle que era desconcertante: la hasta entonces bien establecida conservación de la energía, que deja de ser cierta en esta teoría.

Desde Galileo (básicamente desde el inicio de la ciencia moderna), el concepto de simetría en física ha sido vital. Un grupo de simetría es una colección de un tipo particular de funciones matemáticas (como las traslaciones o las rotaciones) llamadas transformaciones, que al aplicarse a un objeto dado, lo dejan inalterado. Se dice entonces que un grupo de simetría X corresponde a una teoría física cuando las ecuaciones de la teoría no cambian al aplicarles las transformaciones de X. En particular, el grupo de simetría de la mecánica clásica es llamado grupo de Galileo y contiene las transformaciones de traslación, rotación y desplazamientos a velocidad constante. De manera análoga, el grupo de simetría de la relatividad especial es el llamado grupo de Poincaré, que contiene traslaciones, rotaciones y un tipo especial de rotaciones llamados boosts. En el caso de la relatividad general, el grupo de simetría es el de los 'cambios de coordenadas', de modo que la teoría es exactamente la misma sin importar con qué coordenadas se le vea (lo que los físicos llaman covariancia general y que hace a la teoría naturalmente formulable en términos tensoriales, lo que seguido es atemorizante cuando uno como estudiante la quiere aprender).

Una diferencia crucial entre los cambios de coordenadas y las traslaciones, las rotaciones o los 'boosts', es que de los primeros existe una infinidad, mientras que de los otros existe solamente una cantidad finita con la que se puede generar cualquier otro. Este detalle es básicamente el que hace que la conservación de energía no se cumpla en la relatividad general, pues la 'energía', o más precisamente el tensor de energía-momento, termina dependiendo de las coordenadas que se usan para calcularlo. Emmy Noether, pues, entendió esta diferencia esencial en lo que Hilbert llamaba 'teoremas de energía propios e impropios' y lo empleó para formular su teorema, que en realidad son dos, uno para cada caso.

Para entender con mayor profundidad la historia del desarrollo del teorema ciertamente se requiere al menos cierta familiaridad con la relatividad general; un recuento bastante interesante y que vale la pena leer es el siguiente:
E. Noether's Discovery of the Deep Connection Between Symmetries and Conservation Laws
Nina Byers
arXiv:physics/9807044 [physics.hist-ph]

Emmy Noether proved two deep theorems, and their converses, on the connection between symmetries and conservation laws. Because these theorems are not in the mainstream of her scholarly work, which was the development of modern abstract algebra, it is of some historical interest to examine how she came to make these discoveries. The present paper is an historical account of the circumstances in which she discovered and proved these theorems which physicists refer to collectively as Noether's Theorem. The work was done soon after Hilbert's discovery of the variational principle which gives the field equations of general relativity. The failure of local energy conservation in the general theory was a problem that concerned people at that time, among them David Hilbert, Felix Klein, and Albert Einstein. Noether's theorems solved this problem. With her characteristically deep insight and thorough analysis, in solving that problem she discovered very general theorems that have profoundly influenced modern physics.
Asimismo en el arXiv se encuentra una traducción al inglés del artículo original de Noether:
Invariant Variation Problems
Emmy Noether, M. A. Tavel

The problems in variation here concerned are such as to admit a continuous group (in Lie's sense); the conclusions that emerge from the corresponding differential equations find their most general expression in the theorems formulated in Section 1 and proved in following sections. Concerning these differential equations that arise from problems of variation, far more precise statements can be made than about arbitrary differential equations admitting of a group, which are the subject of Lie's researches. What is to follow, therefore, represents a combination of the methods of the formal calculus of variations with those of Lie's group theory. For special groups and problems in variation, this combination of methods is not new; I may cite Hamel and Herglotz for special finite groups, Lorentz and his pupils (for instance Fokker), Weyl and Klein for special infinite groups. Especially Klein's second Note and the present developments have been mutually influenced by each other, in which regard I may refer to the concluding remarks of Klein's Note.

Null geodesics and affine parameters

I've been away for a while preparing my IELTS test and I hope I'll be ok, whatever my fate is! So, concerning this entry, it came to mind this Physics SE answer I wrote some time ago, where the original poster commented he was only concerned with proper time. In turn I posted this link to another Physics SE answer which deals with what an affine parameter is.

So, the thing is pretty straightforward: the definition of proper time is just the arc length
\begin{equation}\tau=\int\sqrt{-ds^2}\end{equation} and for null geodesics, (by definition) $ds^2=0$, so that's why "proper time assigns the same value to all points (on a null) geodesic". For spacelike and timelike geodesics, $ds^2\neq0$ and $\dot{s}^2=\pm1$ when proper time is the parameter.

Now, recall the geodesic equations for a curve $x=x(\lambda)$ on an arbitrary parameter $\lambda$ with a tangent vector $U\equiv\frac{dx}{d\lambda}$,
\begin{equation}\left(\frac{dU}{d\lambda}\right)^\alpha=\frac{d^2x^\alpha}{d\lambda^2}+\frac{dx^\beta}{d\lambda}\frac{dx^\mu}{d\lambda}{\Gamma^\alpha}_{\mu\beta}=0\end{equation} usually written shorthand as
\begin{equation}(\nabla_{U}{U})^\alpha\equiv{U}^\beta\nabla_\beta{U}^\alpha=0\label{dagger}\end{equation}
So, let's now consider another arbitrary parameter ${\sigma=\sigma(\lambda)}$. Then,
\begin{align}\frac{d}{d\lambda}&=\frac{d\sigma}{d\lambda}\frac{d}{d\sigma}\nonumber\\[0.2in]
\frac{d^2}{d\lambda^2}&=\left(\frac{d\sigma}{d\lambda}\frac{d}{d\sigma}\right)\left(\frac{d\sigma}{d\lambda}\frac{d}{d\sigma}\right)\nonumber\\&=\frac{d\sigma}{d\lambda}\left[\frac{d}{d\sigma}\left(\frac{d\sigma}{d\lambda}\right)\frac{d}{d\sigma}+\frac{d\sigma}{d\lambda}\frac{d^2}{d\sigma^2}\right]\nonumber\\&=\frac{d\sigma}{d\lambda}\left[\frac{\left(\frac{d^2\sigma}{d\lambda^2}\right)}{\frac{d\sigma}{d\lambda}}\frac{d}{d\sigma}+\frac{d\sigma}{d\lambda}\frac{d^2}{d\sigma^2}\right]\end{align} and the geodesic equations on $\sigma$ are
\begin{equation}\frac{\left(\frac{d^2\sigma}{d\lambda^2}\right)}{\frac{d\sigma}{d\lambda}}\frac{dx^\alpha}{d\sigma}+\frac{d\sigma}{d\lambda}\frac{d^2x^\alpha}{d\sigma^2}+\frac{d\sigma}{d\lambda}\frac{dx^\beta}{d\sigma}\frac{dx^\mu}{d\sigma}{\Gamma^\alpha}_{\mu\beta}=0\end{equation} i.e., from (\ref{dagger}),
\begin{equation}U^\beta\nabla_\beta{U}^\alpha=-\frac{\left(\frac{d^2\sigma}{d\lambda^2}\right)}{\left(\frac{d\sigma}{d\lambda}\right)^2}U^\alpha\end{equation} so that from the start, in principle, one could've defined the geodesic equations as ${\nabla_UU\propto{U}}$; however, if
\begin{equation}\frac{d^2\sigma}{d\lambda^2}=0\,\,\Longrightarrow\,\,\sigma=a\lambda+b\label{ddagger}\end{equation} (\ref{dagger}) is recovered. These $\sigma$'s are called affine parameters, and this is why they have this simple form.

Now, (I think) this is where trouble comes. Some authors, e.g. Carroll, define an affine parameter as any parameter related to the proper time $\tau$ as $\sigma=a\tau+b$. The thing is that also, usually it is said that one should use an affine parameter for null geodesics (as in the previous Physics SE answer I mentioned before), but still with this last definition you end up with a constant parameter. So what is generally understood as an affine parameter is simply that with which (\ref{dagger}) is satisfied. For null geodesics, (\ref{dagger}) is trivially satisfied with $\tau$ (or any linear combination whatsoever) so that this parameter is indeed useless; however one is free to use any affine parameter as defined in (\ref{ddagger}) with non-constant $\lambda\in\mathbb{R}$ regardless of its physical relevance.

Now, the term $\nabla_{U}{U}=\frac{dU}{d\lambda}=\frac{d^2x}{d\lambda^2}$ is indeed the acceleration of an observer along the curve $x=x(\lambda)$, so that when using an affine parameter, the observer is not accelerating, while au contraire, with a non-affine parameter, the observer will be accelerated parallel to the direction of movement. This is the main relevance of this sort of parameter; in General Relativity it's the custom to deal with affine parameters where $\tau$, of course, is our favorite one for timelike geodesics, whereas with null geodesics people usually don't care what the parameter is as long as it works, at the end it's all about a goofy parameter ;-)

Gravedad de superficie en un espacio esféricamente simétrico y estático

Quizá la característica más sorprendente de los Agujeros Negros, es que exhiben propiedades termodinámicas. Cuenta la historia que alrededor del año 1970, el físico estadounidense John A. Wheeler (el responsable también de popularizar el término Agujero Negro o Black Hole) sugirió la pregunta: ¿Qué sucede si tiro una taza de té en un Agujero Negro? El meollo del asunto es qué es lo que pasa con la entropía del té caliente, ¿acaso simplemente se pierde para siempre?

Ya antes se había tenido indicaciones de que los Agujeros Negros presentaban propiedades termodinámicas únicamente a partir de propiedades clásicas empleando la Relatividad General, como el llamado teorema del área y otras propiedades dinámicas que tenían una perfecta analogía con las leyes de la Termodinámica Clásica. Uno de los alumnos de Wheeler, llamado Jacob Bekenstein (nacido en México, por cierto), pocos años después, junto con (de manera independiente) Stephen Hawking obtendrían el famoso resultado
\begin{equation}T_H\equiv\frac{\kappa}{2\pi},\hspace{0.5in}S_{BH}\equiv\frac{A}{4}\end{equation} (en unidades $\hbar=k_B=c=G=1$) donde $\kappa$ es la llamada gravedad de superficie y $A$ es el área superficial del Agujero Negro. Este resultado tiene gran relevancia, y vale la pena notar que estas cantidades nos dan información acerca del Agujero Negro entero únicamente conociendo su frontera. En cuanto a $\kappa$ nos bastará pensarla como la aceleración debida a la gravedad experimentada por una partícula de prueba (de masa despreciable) en la superficie del Agujero Negro (en general en la superficie de cualquier objeto).

Lo que no había podido lograr es calcular $T_H$ para Agujeros Negros tipo Schwarzschild, que es el resultado más sencillo que suele mostrarse en la literatura. El cálculo también es sencillo, aunque uno debe sacar un poco el colmillo para hacerlo ;-)

La definición en el contexto de la Relatividad General puede basarse en la idea de un horizonte de Killing, i.e. una hipersuperficie donde un vector de Killing $V$ de la métrica se anula. Para espaciotiempos estacionarios y asintóticamente planos (i.e. que en infinito se vuelven Minkowski), como Schwarzschild, el vector de Killing $t$-traslacional puede normalizarse tal que $V_\mu{V}^\mu=-1$ en infinito. Dado que $V^\mu{V}_\mu=0$ en el horizonte de Killing y que $V$ es normal a este mismo horizonte, entonces $\nabla_\nu(V^\mu{V}_\mu)$ también es normal al horizonte, de modo que $\nabla_\nu(V^\mu{V}_\mu)\propto{V}_\nu$ o bien, de manera equivalente (empleando $\nabla_\nu(V^\mu{V}_\mu)=0$ y la ec. de Killing $\nabla_{(\mu}{V}_{\nu)}=0$), existe algún $\kappa$ tal que
\begin{equation}V^\mu\nabla_\mu{V}_\nu=\kappa{V}_\nu\end{equation} que es la gravedad de superficie, además, por el llamado teorema de Fröbenius (véase e.g. el clásico libro de Robert M. Wald), se tiene que también en el horizonte,
\begin{equation}V_{[\mu}\nabla_\nu{V}_{\sigma]}=0\end{equation} y finalmente contrayendo esta última ecuación con $\nabla^\mu{V}^\nu$ y usando la ec. de Killing, se obtiene finalmente que
\begin{equation}\kappa^2=-\frac{1}{2}(\nabla^\mu{V}^\nu)(\nabla_\mu{V}_\nu)\end{equation}

En el caso de espacios esféricamente simétricos y estáticos es relativamente sencillo calcular la gravedad de superficie. Considérese en general toda métrica de la forma
\begin{equation}ds^2=g_{tt}(r)dt^2+g_{rr}(r)dr^2+r^2d\Omega^2\end{equation}

Sabemos que $V=\p_t=(1,0,0,0)$ es un vector de Killing (nótese que en efecto $V^\mu{V}_\mu=-1$ en $r=r_s$, i.e. en el horizonte de Killing). Primero véase que sólo los términos $\nabla^r{V}^t$ y $\nabla^tV^r$ (y la contraparte con índices covariantes) contribuyen; aquí es donde uno saca el colmillo ;-) ya que uno podría precipitarse a pensar que $\nabla^tV^r\stackrel{?}{=}0$ (incorrectamente, lo que yo hice), para ello se necesitan entonces las conexiones
\begin{align}
{\Gamma^t}_{rt}=\frac{1}{2}g^{tt}(\p_rg_{tt}),\hspace{0.5in}
{\Gamma^r}_{tt}=-\frac{1}{2}g^{rr}\left(\p_rg_{tt}\right)
\end{align} de donde se sigue que (usando que $g^{\alpha\alpha}=(g_{\alpha\alpha})^{-1}$ en este caso),
\begin{align}\nabla^rV^t&=g^{rr}{\Gamma^t}_{rt}=\frac{\p_rg_{tt}}{2g_{tt}g_{rr}}\\
\nabla^tV^r&=g^{tt}{\Gamma^r}_{tt}=-\frac{\p_rg_{tt}}{2g_{tt}g_{rr}}
\end{align} es decir, $\nabla^rV^t=-\nabla^tV^r$, lo que ya se sabía, siendo $V=\p_t$ un vector de Killing. Así entonces, finalmente
\begin{align}\kappa^2&=-\frac{1}{2}\left(\nabla^tV^r\nabla_tV_r+\nabla^rV^t\nabla_rV_t\right)=\frac{(-\p_rg_{tt})^2}{-4(g_{tt}g_{rr})^3}\label{gravedaddesup}\end{align} lo que da como resultado para Schwarzschild (siendo $r_s$ el radio del horizonte de Killing),
\begin{equation}\kappa^2=\frac{1}{4}\lim_{r\to{r}_s}\left(\frac{r_s}{r^2}\right)^2=\frac{1}{4r_s^2}\,\Longrightarrow\,\kappa=\frac{1}{2r_s}=\frac{1}{4M}\end{equation} lo que concuerda con el resultado $T_H=1/8\pi{M}$. Esto por supuesto es general, con ello se puede calcular e.g. que $\kappa=\frac{M}{R^2}$ para el caso Newtoniano (Schwarzschild con $r\to{R}$), $\kappa=\frac{\alpha^2+3r_+^2}{2\alpha^2r_+}$ en AdS-Schwarzschild (tomando en cuenta por supuesto que $r_+$ es el horizonte en este caso) o bien $\kappa=\frac{\sqrt{M^2-Q^2}}{\left(M+\sqrt{M^2-Q^2}\right)^2}$ para un Agujero Negro tipo Reissner–Nordström (no extremal y con $-g_{tt}=g_{rr}^{-1}=1-\frac{2M}{r}+\frac{Q^2}{r^2}$), lo que puede comprobarse fácilmente.

Geodésicas en AdS

Hace poco publiqué esta respuesta en Physics SE, misma que no fue elegida por el interrogador, que eligió en su lugar otra innecesariamente embrollada matemáticamente, como si el interrogador hubiese estado buscando cátedra en geometría diferencial cuando ignoraba siquiera lo que es un parámetro afín.

De cualquier modo, yo también me lié un poco con mi respuesta. En la parte de las geodésicas tipo luz básicamente copié el cálculo que hice en mi trabajo de servicio social, pero en el caso de las geodésicas tipo tiempo y tipo espacio me encontré con algunas dificultades (en realidad no pensé en la necesidad o plausibilidad de calcularlas explícitamente en el servicio social) y no pude ahondar más porque estaba en pleno periodo trimestral en la Universidad.

El objetivo de esta entrada entonces es básicamente reproducir analíticamente las geodésicas en el esquema del diagrama de Penrose de AdS que muestro acá:
Esquema hecho con TikZ: El espacio AdS
y graficar con algún software numérico, e.g. Mathematica, para ello voy a emplear la métrica (con variables angulares fijas)
\begin{equation}ds^2=-(1+r^2)dt^2+(1+r^2)dr^2\end{equation} con $r=\tan\mc{R}\in\mathbb{R}^+$, y para la cubierta universal, $t\in\mathbb{R}$.

Para las geodésicas de tipo tiempo, tomando $\lambda$ como el tiempo propio, ya que (tomando signatura +---) $\dot{s}^2=-1$, se tiene que
\begin{equation}-(1+r^2)\dot{t}^2+(1+r^2)\dot{r}^2=-1\end{equation} Ahora bien, $\p_t$ (i.e. el vector de componentes $V^t=1$, $V^i=0$) es un vector de Killing global, entonces
\begin{equation}V_t\dot{t}=g_{tt}\dot{t}=-(1+r^2)\,\dot{t}\equiv-E\label{da}\end{equation} es una constante de movimiento (la energía, siendo que proviene de una simetría $t$-traslacional); donde el signo negativo es simplemente tal que $\dot{t}>0$. De aquí entonces
\begin{equation}\dot{r}^2+(1+r^2)-E^2=0\end{equation} es decir,
\begin{equation}\frac{dr}{\pm\sqrt{(E^2-1)-r^2}}=d\lambda\end{equation} cuya solución más general, con la condición (arbitraria) $r(0)=0$, es
\begin{equation}r(\lambda)=\pm\sqrt{E^2-1}\tan\lambda\,|\cos\lambda|\end{equation} que de cualquier modo, restringiendo $r(\lambda)>0$, se tiene que reducir inevitablemente el dominio de $\lambda$ a intervalos de $\pi/2$ respecto a $\lambda=0$, es decir,
\begin{equation}r(\lambda)=\begin{cases}\sqrt{E^2-1}\tan\lambda\,|\cos\lambda|,&\text{si}&\lambda\in[n\pi,\frac{(2n+1)\pi}{2}]\\[0.1in]-\sqrt{E^2-1}\tan\lambda\,|\cos\lambda|,&\text{si}&\lambda\in[\frac{(2n+1)\pi}{2},(n+1)\pi]\end{cases},\,n\in\mathbb{Z}\end{equation} esto por supuesto de algún modo está relacionado con que AdS no está extendido a-priori, o que está enrollado en la coordenada temporal. Tomando cualquier caso, de la ec. (\ref{da}) con $t(0)=0$, se sigue que
\begin{equation}t(\lambda)=\arctan(E\,\tan\lambda)\end{equation} y así entonces escribiendo $\lambda=\lambda(t)$, se sigue finalmente que (asumiendo $E>0$; si se asume lo contrario los intervalos se invierten correspondientemente),
\begin{equation}r(t)=\begin{cases}\sqrt{E^2-1}\frac{\tan{t}}{E\sqrt{1+\frac{\tan^2t}{E^2}}},&\text{si}&t\in[n\pi,\frac{(2n+1)\pi}{2}]\\[0.1in]-\sqrt{E^2-1}\frac{\tan{t}}{E\sqrt{1+\frac{\tan^2t}{E^2}}},&\text{si}&t\in[\frac{(2n+1)\pi}{2},(n+1)\pi]\end{cases},\,n\in\mathbb{Z}\end{equation} y simplemente hay que usar $\mc{R}=\arctan{r}$. De aquí uno puede procurar simplificar cosas e.g. tomar intervalos de $\pi$ respecto a $t=0$ haciendo $\tan{t}\to|\tan{t}|$ (esto equivale a lo que escribí en la respuesta de Physics SE para $r(\lambda)$), lo que puede justificarse siendo removible la discontinuidad entre cada $\pi/2$ respecto a 0; de algún modo esto recupera la $2\pi$-periodicidad de AdS en la cubierta universal. Esta cuestión de los intervalos vuelve un poco delicado el tratamiento de este problema de obtener analíticamente las expresiones para las geodésicas; el caso de las geodésicas de tipo espacio es enteramente análogo tomando $\dot{s}^2=1$ en el tiempo propio $\lambda$, pero nuevamente uno debe ser cuidadoso al considerar en dónde están definidas las soluciones físicamente relevantes.

Graficando las soluciones entonces para diversos valores de $E$, se obtiene (extendiendo a la cubierta universal)
como se esperaba, si se quiere rotar el gráfico, simplemente hay que graficar $t=t(\mc{R})$.

Finalmente una cuestión relevante es que, independientemente del signo que tenga, en el caso de geodésicas tipo tiempo, la cantidad $E$ debe satisfacer $|E|\geq1$. De esto no he encontrado nada al respecto en otros lados y no lo he podido discutir ampliamente con mi asesor, sin embargo no es descabellado que estas geodésicas tengan una restricción en la energía pues implican movimiento masivo y la cantidad 1 está asociada a un cambio de longitud de arco en el tiempo propio, de modo que para que una partícula se mueva en una de estas geodésicas es coherente que su energía deba exceder o igualar en proporción esta razón de cambio (en unidades correspondientes), si la iguala, simplemente no se mueve (i.e. se queda en el polo $\mc{R}=0$). También debe notarse que cuando $E\to\infty$, las geodésicas se tornan de tipo luz, como también se esperaba.

Agujeros Negros en AdS

He estado ausente algún tiempo intentando terminar la licenciatura, y como mostré en la entrada anterior, he estado trabajando en la geometría (clásica) del espacio de Anti-de Sitter. Ahora quiero compartir algunos comentarios sobre Agujeros Negros en AdS, basándome en este ensayo de Peng Zhao.

El caso más sencillo es (A)dS-Schwarzschild, y la idea es que partiendo de una métrica esféricamente simétrica en forma de Schwarzschild,
\begin{equation}ds^2=-f(r)\,dt^2+f^{-1}(r)\,dr^2+r^2d\Omega^2\end{equation} uno puede tomar $f$ tal que cerca del origen se tenga un espacio de Schwarzschild y asintóticamente se tenga (A)dS (para una obtención más metódica de esta métrica véase e.g. la p.823 de las notas de M. Blau). Las métricas de AdS y de Schwarzschild, respectivamente, pueden escribirse como
\begin{align}ds_\text{AdS}^2&=-\left(1+\frac{r^2}{\alpha^2}\right)\,dt^2+\left(1+\frac{r^2}{\alpha^2}\right)^{-1}dr^2+r^2d\Omega^2\\
ds_\text{Sch}^2&=-\left(1-\frac{r_s}{r}\right)\,dt^2+\left(1-\frac{r_s}{r}\right)^{-1}dr^2+r^2d\Omega^2\end{align} donde $\alpha$ es una constante positiva (el radio de AdS) y $r_s$ es el radio de Schwarzschild. De aquí, se define entonces la métrica de AdS-Schwarzschild como
\begin{equation}ds^2=-\left(1+\frac{r^2}{\alpha^2}-\frac{r_s}{r}\right)\,dt^2+\left(1+\frac{r^2}{\alpha^2}-\frac{r_s}{r}\right)^{-1}dr^2+r^2d\Omega^2\end{equation} que en efecto se comporta como un agujero negro de Schwarzschild para $r$ pequeña y como AdS para $r$ grande. Lo propio puede hacerse empleando de Sitter.

***(Actualización)
La primera cuestión es a AdS se le puede asociar una energía finita y por tanto el análogo a una pared de potencial infinita en su infinito asintótico. Por teorema de Noether, en AdS la energía es una carga conservada, $E=-P_\mu{V}^\mu$ (el signo negativo es simplemente por consistencia, de cualquier modo siempre se le puede cargar a $E$ que es constante); en este caso nos interesa la energía $E_\text{loc}$ medida por un observador local (es decir, que localmente está en un sistema plano) estático, en donde tendríamos simplemente que $P=(E,\vec{0})$ y $V=(1,\vec{0})$. Lo que nos interesa es encontrar esta relación para este observador en el sistema global de coordenadas de AdS (esféricas). Considérense los vectores transformados $p$ y $v$. Primero empleemos la condición de normalización $v^\mu{v}_\mu=-1$ y el hecho de que el observador es estático también en el sistema global de AdS, i.e. $v=(v^t,\vec{0})$, esto implica que $v^t=1/\sqrt{-g_{tt}}$, entonces en este sistema coordenado,
\begin{equation}E_\text{loc}=-p_\mu{v}^\mu=-g_{\mu\nu}p^\mu{v}^\nu=\sqrt{-g_{tt}}p^t\end{equation} Ahora bien, en el sistema del observador local estático simplemente $p^t=m\dot{t}$ con $m$ la masa de la partícula (el observador), que podemos llevar a la forma en coordenadas de AdS sabiendo que $P^t=g_{tt}\dot{t}=-E$ en AdS, de modo que $p^t=-mE/g_{tt}\equiv-{E}_0/g_{tt}$ y así,
\begin{equation}E_\text{loc}(r)=\frac{E_0}{\sqrt{-g_{tt}(r)}}\label{1}\end{equation} de donde se lee que la energía medida por un observador local en algún punto $r$ tiene un corrimiento al rojo (debido al campo gravitacional) respecto a la energía $E_0$. En el ensayo de Zhao se escribe $E_0$ como $E_\infty$, sin embargo la notación se antoja un poco desafortunada, ya que ésta es la energía total medida por un observador local estático en $r=0$; la notación probablemente proviene de la interpretación de ser la energía medida por un observador alejado infinitamente del campo gravitacional del espacio, lo que resulta útil e.g. en Schwarzschild, donde $E$ en efecto es la energía (por unidad de masa) medida por un observador local estático en $r\to\infty$, aquí de hecho $E_\text{loc}\to0$ cuando $r\to\infty$, contrario a Schwarzschild, en donde $E_\text{loc}\to{mE}$ cuando $r\to\infty$ y $E_\text{loc}\to\infty$ cuando $r\to{r}_s$, esto hace la gran diferencia para ver a AdS como una caja. Nótese también que los vectores de momento y de velocidad (i.e. el vector de Killing transformado) se recuperan para el caso local en AdS con $r=0$ y en Schwarzschild con $r\to\infty$. Como se ve adelante, pueden construirse estados térmicos en AdS identificando periódicamente una coordenada de tiempo imaginaria con periodo $\beta_0=T^{-1}$. Esto significa entonces que (\ref{1}) puede escribirse también como $T_\text{loc}=\frac{T}{\sqrt{-g_{tt}}}$ y que la energía total de la radiación térmica en AdS es finita sin necesidad de encerrarla en una caja.
***

Luego entonces, tomando $r=r_+$ el horizonte de eventos de AdS-Schwarzschild, i.e. la raíz más grande de
\begin{equation}V(r)\equiv1-\frac{r_s}{r}+\frac{r^2}{\alpha^2}\end{equation} se fija el valor de $r_s$ a
\begin{equation}V(r_+)=0\,\Longrightarrow\,r_s=\frac{r_+(\alpha^2+r_+^2)}{\alpha^2}\end{equation} y entonces en efecto al ver la métrica cerca del horizonte, i.e. con $r=r_++\epsilon^2$ con $|\epsilon|\ll1$, para el orden dominante en $\epsilon$,
\begin{equation}V(r_++\epsilon^2)\approx\frac{\alpha^2+3r_+^2}{\alpha^2r_+}\epsilon^2\end{equation} y haciendo una rotación de Wick, $\tau=it$, se sigue para la métrica AdS-Schwarzschild,
\begin{equation}ds^2=\frac{4\alpha^2r_+}{\alpha^2+3r_+^2}\left[d\epsilon^2+\left(\frac{\alpha^2+3r_+^2}{2\alpha^2r_+}\right)^2\epsilon^2d\tau^2\right]+(r_++\epsilon^2)^2d\Omega^2\end{equation} Aquí puede resultar poco claro cómo ver que necesariamente $\tau$ es periódica. El punto es entender lo que significa que se tenga una singularidad cónica. La idea es básicamente que si $\tau$ tiene un cierto periodo $\beta_0$, el espacio se cubrirá completamente, lo que no sucederá necesariamente con cualquier otro valor; la analogía que utiliza Zhao es la de las coordenadas polares $ds^2=dr^2+r^2d\theta^2$, donde $\theta$ debe tener un periodo de $2\pi$ para corresponder a la métrica del plano, de otro modo se tiene un hueco del tipo de una rebanada de pastel, lo que lleva a algo así como un pico de un cono, lo que es la singularidad cónica.
En este caso es particularmente fácil verificar cuál es el periodo $\beta_0$ de $\tau$, dado que el término entre corchetes de la métrica AdS-Schwarzschild cerca del horizonte es análoga a la del plano en coordenadas polares. Se debe satisfacer entonces que
\begin{equation}\frac{\alpha^2+3r_+^2}{2\alpha^2r_+}\int_0^{\beta_0}d\tau\int_0^{r_+}\epsilon\,{d}\epsilon=\pi{r}_+^2\,\Longrightarrow\,\beta_0=\frac{4\pi\alpha^2r_+}{\alpha^2+3r_+^2}\end{equation}
En adelante la cosa es bastante clara, al menos antes de llegar a la parte de AdS/CFT, y el resultado $T=\beta^{-1}$ creo que es bastante majo, pues entonces la temperatura del agujero negro puede obtenerse únicamente con argumentos geométricos, sin invocar directamente la gravedad de superficie. Sólo hay detalles cuando se menciona la distribución canónica, aunque al final igual lo que se calcula es $S=\beta\langle{E}\rangle+\ln{Z}$. La intención es llegar a la dualidad fluido/gravedad, pero de eso tal vez escriba después  ;-)

La geometría clásica del espacio de (Anti)-de Sitter

He pasado poco más de seis meses entrándole a la geometría (clásica) de los espacios (A)dS en mi proyecto de servicio social en la licenciatura. El producto son estas notas de casi 200 páginas que cubren los temas de Relatividad Especial, Relatividad General y los espacios de (Anti)-de Sitter. La razón de haberlo hecho así es que en realidad tuve que aprender todos estos temas a lo largo del servicio. Las notas me parecen bien digeridas para cualquier estudiante de licenciatura interesado en temas de gravitación y no creo que haya ningún problema para cualquiera que las quisiera leer. El trabajo presuntamente está libre de errores conceptuales mayores y a lo más podrían existir errores estéticos o de dedo, mismos que procuraré señalar aquí al encontrarlos.

Resumen
Los espacios de de Sitter y de Anti-de Sitter son básicamente espaciotiempos de curvatura constante con ciertas propiedades peculiares que resultan físicamente relevantes. Estos espacios pueden verse como una solución de las ecuaciones de campo de la Relatividad General cuando se considera una constante cosmológica ya sea positiva o negativa. Si bien el espacio ${AdS}$ juega un papel limitado en la Relatividad General, contrario al espacio ${dS}$, éste puede resultar de gran relevancia en otras teorías que incorporan gravedad como lo es teoría de cuerdas dentro del marco de la dualidad AdS/CFT. Aunque este trabajo no es exhaustivo, se busca introducir primeramente los elementos físicos y matemáticos básicos de la Relatividad Especial y General, dado que estos temas actualmente pertenecen únicamente a cursos optativos en la UAM-I. En general se cubren tres secciones básicas: 1) Relatividad Especial, 2) Relatividad General, y 3) Geometría y el espacio $(A)dS$.

En la sección 1 se introduce la notación tensorial y algunos conceptos básicos sobre tensores, luego se hace un resumen de los conceptos básicos de la Relatividad Especial y como complemento se tratan brevemente la electrodinámica relativista y las ecuaciones cuántico-relativistas. Finalmente se muestran diversos problemas resueltos, cuya formulación ha sido traducida de las diversas fuentes bibliográficas citadas.

En la sección 2 se presentan de manera más sólida los conceptos acerca de tensores y se introduce el tema de variedades diferenciables y curvatura a manera de formular la teoría de manera tensorial. Finalmente se obtienen las ecuaciones de campo por formulación Lagrangiana y asimismo se muestran algunos problemas resueltos, con el principal objetivo de complementar el contenido de la sección.

En la sección 3 se tratan finalmente los espacios clásicos de ${(A)dS_n}$ introduciendo diversos conceptos geométricos que aquí se dirigen principalmente a los campos de Killing y los espacios máximamente simétricos. A partir de aquí se estudian diversos aspectos de ${(A)dS_n}$ en varios sistemas coordenados y se ilustran los correspondientes diagramas de Penrose a modo de estudiar la estructura causal de cada espacio. Ambos espacios en los distintos sistemas coordenados se logran a través de diversas parametrizaciónes de hipersuperficies que luego se encajan dentro de un espacio coordenado de dimensión mayor. El tratamiento en ambos casos se hace en un número de dimensiones arbitrarias, de modo que el lector pueda simplemente considerar cualquier caso de interés. Finalmente se tratan los espacios ${(A)dS_n}$ en el contexto de la Relatividad General, considerando el papel que juegan desde un punto de vista cosmológico y se revisa la solución de $AdS$-Schwarzschild.

Enlace al documento en la Colección de Tesis Electrónicas de la UAM-I:
http://tesiuami.izt.uam.mx/uam/aspuam/presentatesis.php?recno=16576&docs=UAMI16576.pdf

Esquema hecho con TikZ: El espacio AdS

Atractores Cosmológicos

En la UAM no existe tal cosa como la titulación por tesis a nivel licenciatura (que yo sepa), sin embargo existe algo análogo llamado Proyecto Terminal en donde uno básicamente tiene que desarrollar un tema en el que esté interesado y enfrentarse a las dificultades que implica la investigación científica.

Soy de la idea de que sería algo verdaderamente bueno el que hubiera acceso libre al trabajo generado en estos proyectos, al menos por parte de la misma comunidad; hay compañeros que parecen estudiar cosas bastante interesantes y generar contenido que puede resultar útil, e.g. código fuente, además de proporcionar algo así como una guía para quien apenas está en vías de decidir qué camino tomar. Aunque hay algo parecido, parece estar desactualizado, además de estar poco promovido y -a mi parecer- resultar poco atractivo. Como sea, lo siguiente es mi trabajo del primer proyecto (de dos), que cubre temas que no son vistos regularmente en la licenciatura (dos de ellos están escritos aquí: Sistemas Hamiltonianos Singulares y El modelo FLRW) y que afortunadamente entra en armonía con las noticias recientes acerca de inflación cosmológica.

Resumen
La naturaleza de las soluciones de tipo atractor en cosmología inflacionaria ha sido controversial, ya que este aparente comportamiento en sistemas cosmológicos con campo escalar y potencial inflacionario contradiría el teorema de Liouville, que afirma que el volumen de cualquier región dada en el espacio fase permanece invariante. En este trabajo se presentan algunos de los elementos físicos y matemáticos básicos necesarios para entender el concepto de atractores en modelos cosmológicos con campo escalar, así como para poder abordar algunos resultados recientes en el tema, como son [1] y [2]. Los principales temas tratados son: i) Sistemas Hamiltonianos Singulares, ii) El modelo de FLRW, iii) Las ecuaciones de Einstein acopladas a un campo escalar y iv) Inflación, Medidas y Atractores. Este trabajo cubre la primera parte de lo que se planea como dos proyectos terminales que se cubrirán en los dos trimestres consecutivos 14-I y 14-P. Por facilidad y plausibilidad se cita eventualmente código escrito para el software Wolfram Mathematica 8 al realizar diversos cálculos.
El enlace: Atractores Cosmológicos I

El modelo de Friedmann-Lemaître-Robertson-Walker

La esencia de la relatividad general es que el espaciotiempo es una variedad cuatrodimensional dotada de una métrica ${g_{\alpha\beta}}$ de signatura de Lorentz (-+++) relacionada con la densidad de materia vía las ecuaciones de Einstein ${G_{\alpha\beta}=8\pi{G}T_{\alpha\beta}}$. Una cuestión vital de la relatividad general es conocer qué solución a las ecuaciones de Einstein corresponde a nuestro universo, o al menos a un modelo idealizado de nuestro universo.

El principio cosmológico declara que en gran escala el universo es isótropo y homogéneo. La principal evidencia para considerar este principio es la estabilidad de la temperatura de la radiación cósmica de fondo o radiación de fondo de microondas (CMB por sus siglas en inglés: cosmic microwave background) medida por el explorador COBE y la sonda WMAP. La forma de obtener el modelo de Friedmann-Lemaître-Robertson-Walker o simplemente FLRW, es precisamente aplicar las simetrías de homogeneidad e isotropía a las ecuaciones de Einstein. En seguida se investiga entonces la estructura del espaciotiempo bajo estas condiciones.

Homogeneidad e Isotropía
Otro principio que germinó desde el tiempo de Copérnico es el que dice que no ocupamos (los humanos) una posición privilegiada en nuestro universo, y que si nos moviéramos a alguna otra región del mismo, las características básicas de los alrededores parecerían ser las mismas. Asimismo suele asumirse que no importa en qué dirección miremos, las observaciones a gran escala serán básicamente las mismas. Esto constituye el llamado principio de Copérnico.


Estas ideas son de algún modo hipótesis que se hacen sobre la estructura del universo y corresponden al concepto de homogeneidad e isotropía, respectivamente. Como se menciona en la introducción, estos principios han recibido apoyo por distintas observaciones, lo que las hace un buen candidato de estudio.

Un espaciotiempo ${(\mathscr{M},g)}$ se dice homogéneo si existe una foliación ${\sigma_t}$ de hipersuperficies tipo espacio en un parámetro $t$ tal que para cada $t$ y cualesquiera puntos ${p,q\in\sigma_t}$ existe una isometría $\phi$ tal que ${\phi(p)=q}$. Esto significa que en cualquier instante de tiempo todo punto del espacio debe verse como cualquier otro.

El espaciotiempo ${(\mathscr{M},g)}$ será isótropo si existe una foliación $\alpha$ de curvas congruentes tipo tiempo con vector tangente $u$ tales que para cualquier ${p\in\alpha}$ y cualesquiera vectores tangente unitarios tipo espacio ${v,w\in{T}_p(\mathscr{M})}$ existe una isometría $\varphi$ tal que ${\varphi(p)=p}$, ${\varphi^*(u)=u}$ y ${\varphi^*(v)=w}$. Esto significa que en cualquier punto es imposible construir vectores espaciales preferidos.

Las hipersuperficies homogéneas ${\sigma_t}$ son ortogonales a cada vector $u$, es decir, a las líneas de mundo de los observadores isotrópicos. Véase que isotropía para todos los observadores implica homogeneidad para todos los observadores. Es posible construir universos homogéneos y anisotrópicos pero no universos inhomognéneos e isotrópicos. Un universo isótropo implica que no hay un centro para el universo. Esto es relevante al considerar el origen del universo o Big Bang, pues debido a la isotropía, no hay un punto privilegiado o centro en el cual éste haya ocurrido.

Por supuesto el universo no es exactamente homogéneo; hay irregularidades locales como las estrellas y las galaxias, al referirse a que lo es en gran escala, uno se refiere en un orden de millones de años luz. El principio de Copérnico entonces puede traducirse en que a gran escala el universo es esféricamente simétrico en cada punto.

La métrica FLRW
En 1944 Arthur Geoffrey Walker mostraría que una simetría esférica exacta alrededor de cada punto implicaría que el universo es espacialmente homogéneo y que admitiría un grupo de seis parámetros de isometrías cuyas superficies de transitividad son 3-superficies de curvatura constante ([1]). Tal espacio es el llamado de Friedmann-Lemaître-Robertson-Walker o FLRW (algunas personas suelen llamarlo simplemente con un nombre o dos; lo más conveniente parece ser llamarlo simplemente FLRW). Lo que se quiere entonces es determinar las geometrías de las hipersuperficies ${\sigma_t}$ para todo valor $\kappa$ de curvatura. Uno sólo debe preocuparse por curvatura positiva, nula y negativa, que puede describirse por múltiplos de ${\kappa=-1,0,1}$.

Para los respectivos $\kappa$ se tiene una 3-esfera, un 3-plano y un 3-hiperboloide, cuyas métricas (elementos de línea), en coordenadas esféricas, cartesianas e hiperbólicas, respectivamente, son
\begin{equation}d\sigma^2=\begin{cases}d\psi^2+\sin^2\psi\left(d\theta^2+\sin^2\theta\,d\phi^2\right),&\text{si}\,\kappa=1,\,\psi\in[0,2\pi)\\
dx^2+dy^2+dz^2,&\text{si}\,\kappa=0\\
d\zeta^2+\sinh^2\zeta\left(d\theta^2+\sin^2\theta\,d\phi^2\right),&\text{si}\,\kappa=-1,\,\zeta\in[0,\infty)\end{cases}\end{equation} o bien, en términos de los ángulos polar $\theta$ y azimutal $\phi$,
\begin{equation}d\sigma^2=d\psi^2+\left(d\theta^2+\sin^2\theta\,d\phi^2\right)\begin{cases}\sin^2\psi,&\text{si}\,\kappa=1,\,\psi\in[0,2\pi)\\
\psi^2,&\text{si}\,\kappa=0,\,\psi\in[0,\infty)\\
\sinh^2\psi,&\text{si}\,\kappa=-1,\,\psi\in[0,\infty)\end{cases}\end{equation} que puede escribirse de manera sucinta en términos de la curvatura $\kappa$ como ([2])
\begin{equation}ds^2=\frac{dr^2}{1-\kappa{r}^2}+r^2\left(d\theta^2+\sin^2\theta\,d\phi^2\right)\label{1}\end{equation} con los cambios de variable ${r=\sin\psi,\,\psi,\,\sinh\psi}$ para ${\kappa=1,0,-1}$, respectivamente, únicamente restringida en el primer caso con $|r|<1$. En los casos del 3-plano y el 3-hiperboloide, los hiperespacios serán difeomorfos (un difeomorfismo es un equivalente de un isomorfismo pero en variedades diferenciables) a ${\mathbb{R}^3}$ mientras que para las 3-esferas, éstas serán difeomorfas a ${\mathcal{S}^3}$. En el primer caso los espacios son infinitos y las posibilidades para el universo son llamadas universos abiertos, mientras que en el segundo los espacios son compactos, que son finitos y sin frontera, y naturalmente los universos son llamados universos cerrados. De esta manera, al saber que estas soluciones son buenos candidatos para modelar aproximadamente nuestro universo, se abre la pregunta acerca de si el universo es abierto o cerrado.

De aquí entonces para construir la métrica FLRW, se llevan las coordenadas de cada ${\sigma_t}$ a cada una de las otras hipersuperficies por medio de cada observador isótropo (se asignan coordenadas fijas a cada observador) y se dota de un reloj a cada observador isótropo, i.e. se etiqueta cada hipersuperficie por un tiempo propio $\tau$ (también llamado tiempo cósmico), de modo que todo evento en el universo está dado por $\tau$ y las coordenadas espaciales. De este modo entonces la métrica FLRW es
\begin{equation}ds^2=-d\tau^2+a^2(\tau)d\sigma^2\end{equation} donde el factor $a$ es llamado factor de escala (determina la escala total de la métrica espacial) y puede depender de $\tau$.

Las ecuaciones de Friedmann
La métrica entonces es conocida y se puede conocer fácilmente el tensor de energía-momento dada la propiedad de isotropía del universo FLRW, entonces las ecuaciones de Einstein
\begin{equation}G_{\mu\nu}\equiv{R}_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R=8\pi{G}T_{\mu\nu}\end{equation} serán útiles para determinar la evolución dinámica espacial del universo, i.e. el factor de escala ${a(\tau)}$.

La distribución de materia en este caso debe ser también homogénea e isótropa, de modo que el tensor de energía-momento puede modelarse por un fluido perfecto, i.e. un fluido que está completamente especificado por la densidad de energía propia (i.e. en un marco de reposo) $\rho$ y una presión isotrópica propia $p$, i.e.
\begin{equation}T_{\mu\nu}=(p+\rho)u_\mu{u}_\nu+pg_{\mu\nu}\end{equation} donde ${u_\alpha}$ es la 4-velocidad del fluido. Por isotropía, el fluido debe estar en reposo en las coordenadas comóviles (las coordenadas propias de cada punto del fluido), entonces
\begin{equation}u^\mu=(1,\vec{0})\end{equation} y el tensor de energía-momento es de la forma
\begin{equation}T_{\mu\nu}=\begin{pmatrix}\rho&\vec{0}^\mathrm{\,T}\\
\vec{0}&g_{ij}p\end{pmatrix}\label{4}\end{equation} Para calcular las componentes $G_{\alpha\beta}$ lo más sencillo es recurrir a la computadora, siempre que se entienda qué es lo que se está haciendo. Aquí emplearé nuevamente Mathematica. El siguiente código tiene la ventaja de que funcionará para cualquier métrica. Recuérdese que para los símbolos de Christoffel,
\begin{equation}{\Gamma^\rho}_{\mu\nu}=\frac{1}{2}g^{\rho\sigma}\left(\p_{\mu}g_{\sigma\nu}+\p_{\nu}g_{\sigma\mu}-\p_{\sigma}g_{\mu\nu}\right)\end{equation} lo que en Mathematica puede escribirse como una función para la métrica y que trabaje con las coordenadas apropiadas


Luego para el tensor de Riemann, a partir de
\begin{equation}{R^\rho}_{\sigma\mu\nu}=\p_\mu{\Gamma^\rho}_{\sigma\nu}-\p_\nu{\Gamma^\rho}_{\sigma\mu}
+{\Gamma^\rho}_{\mu\alpha}{\Gamma^\alpha}_{\sigma\nu}-{\Gamma^\rho}_{\nu\alpha}{\Gamma^\alpha}_{\sigma\mu}\end{equation} de manera análoga,


Y así también, para el tensor de Ricci, sabiendo que
\begin{equation}R_{\mu\nu}={R^\alpha}_{\mu\alpha\nu}\end{equation} de manera análoga,


y finalmente para el escalar de Ricci, ya que
\begin{equation}R={R^\alpha}_\alpha=g^{\alpha\beta}R_{\alpha\beta}\end{equation} de manera análoga


y de este modo se puede calcular fácilmente el tensor ${G_{\alpha\beta}}$. Empleando entonces la métrica FLRW general del elemento (\ref{1}),
\begin{equation}g_{\alpha\beta}=\text{diag}\left(-1,\frac{a^2}{1-\kappa{r}^2},a^2r^2,a^2r^2\sin^2\theta\right)\end{equation} en las coordenadas ${(t,r,\theta,\phi)}$, lo que en Mathematica puede escribirse como


De aquí entonces simplemente puede definirse el tensor de Einstein y evaluarse con la métrica y las coordenadas


de modo que se obtiene que
\begin{align}
G_{00}&=3\frac{\kappa+\dot{a}^2}{a^2}\\
G_{11}&=-\frac{\kappa+\dot{a}^2+2a\ddot{a}}{1-\kappa{r}^2}\\
G_{22}&=-r^2\left(\kappa+\dot{a}^2+2a\ddot{a}\right)\\
G_{33}&=-r^2\sin^2\theta\left(\kappa+\dot{a}^2+2a\ddot{a}\right)
\end{align} con todas las demás componentes nulas. De aquí es evidente que
\begin{align}
{G^0}_0&=g^{0\alpha}G_{\alpha0}=-3\frac{\kappa+\dot{a}^2}{a^2}\\
{G^i}_j\delta^j_i&=g^{i\alpha}G_{\alpha{i}}\delta^j_i=-\frac{1}{a^2}\left(\kappa+\dot{a}^2+2a\ddot{a}\right)
\end{align} y además de (\ref{4}),
\begin{equation}{T^\mu}_\nu=\text{diag}(-\rho,p,p,p)\end{equation} de modo que resulta conveniente (de forma equivalente, simplemente es más evidente) escribir las ecuaciones de movimiento de Einstein como
\begin{equation}{G^\mu}_\nu=8\pi{G}{T^\mu}_\nu\end{equation} y las ecuaciones de Einstein se reducen a
\begin{align}
\left(\frac{\dot{a}}{a}\right)^2+\frac{\kappa}{a^2}&=\frac{8\pi{G}}{3}\rho\label{2}\\
2\frac{\ddot{a}}{a}+\left(\frac{\dot{a}}{a}\right)^2+\frac{\kappa}{a^2}&=-8\pi{G}p
\end{align} que a su vez pueden combinarse para producir la ecuación independiente de la curvatura
\begin{equation}\frac{\ddot{a}}{a}=-\frac{4\pi{G}}{3}(\rho+3p)\label{3}\end{equation} A (\ref{2}) junto con (\ref{3}) suele llamársele (y en adelante así se hará aquí) las ecuaciones de Friedmann, o algunos autores le llaman ecuación de Friedmann a (\ref{2}) y a (\ref{3}) ecuación de aceleración.

Recuérdese ahora la interpretación de $a$ como un factor que determina la escala total de la métrica espacial. Un resultado impresionante que se sigue de (\ref{3}), es que el universo no es estático dados ${\rho,p>0}$ para la materia ordinaria. Esto significa que el universo está expandiéndose si ${\dot{a}>0}$ o contrayéndose si ${\dot{a}<0}$ y de modo tal que $a$ en el pasado sea más pequeña hasta que en algún momento ${a=0}$ en donde las componentes del tensor de curvatura divergen y se tiene una singularidad. De esto último se sigue la interpretación del Big Bang, misma que es vigente, sin embargo desde hace algunos años se sabe que el universo de hecho se está expandiendo y no contrayendo ([3]).

Es bien conocida la anécdota de que Einstein se incomodó con el resultado de un universo que no es estático e introdujo la constante cosmológica $\Lambda$, misma que luego llamaría su más grande metida de pata. La idea se deshecharía hasta las observaciones hechas en [3]. Las ecuaciones de Einstein con la modificación $\Lambda$ pueden escribirse como
\begin{equation}G_{\mu\nu}+\Lambda{g}_{\mu\nu}=8\pi{G}T_{\mu\nu}\label{5}\end{equation} de donde se siguen soluciones análogas con el término $\Lambda$, lo que puede cambiar el carácter de (\ref{3}). La constante cosmológica seguido se interpreta como una presión negativa (aunque esto disgusta a más de uno, incluyendo cosmólogos) y de hecho puede incluirse por medio de transformaciones de $\rho$ y $p$, por lo que en adelante salvo que se mencione explícitamente, omitiré el carácter positivo de $p$ y $\rho$ y mantendré el tratamiento general.

Véase luego que (\ref{2}) puede escribirse como,
\begin{equation}\rho{a}^3=\frac{3a}{8\pi{G}}\left(\dot{a}^2+\kappa\right)\end{equation} cuya derivada es
\begin{equation}\dot{\rho}a^3+3a^2\dot{a}\rho=\frac{3\dot{a}}{8\pi{G}}\left(\dot{a}^2+\kappa+2a\ddot{a}\right)\end{equation} y también por (\ref{3}) se tiene que
\begin{equation}-3a^2\dot{a}p=\frac{3\dot{a}}{8\pi{G}}\left(\dot{a}^2+\kappa+2a\ddot{a}\right)\end{equation} y entonces se sigue que
\begin{equation}\dot{\rho}+3\frac{\dot{a}}{a}(\rho+p)=0\label{6}\end{equation} que es la misma ecuación que puede obtenerse por conservación de energía-momento, i.e. resolviendo
\begin{align}\nabla_\mu{T}^{\mu\nu}=\p_\mu{T}^{\mu\nu}+T^{\alpha\nu}\Gamma^{\mu}_{\alpha\mu}=0\end{align} que nuevamente empleando el código escrito en Mathematica puede calcularse fácilmente.

Dos soluciones particulares
En la escala cósmica que aquí se ha referido como gran escala, las galaxias pueden idealizarse como granitos de polvo cuyas velocidades aleatorias son tan pequeñas que la presión de este polvo es despreciable, de modo que
\begin{equation}p=0\end{equation} Sin embargo, además de materia no relativista, el universo contiene radiación electromagnética que ejerce una presión no nula
\begin{equation}p=\frac{\rho}{3}\end{equation} con $\rho$ la densidad de energía (véase e.g. [4]).

Por la ecuación de conservación de energía-momento (\ref{6}), se sigue que para el caso de polvo ${p=0}$,
\begin{equation}\underbrace{\dot{\rho}a^3+3\rho\dot{a}a^2}_{\frac{d}{dt}(\rho{a}^3)}=0\end{equation} es decir
\begin{equation}\rho{a}^3=\text{cte}\end{equation} mientras que para el caso de radiación ${p=\rho/3}$, se sigue de manera análoga, que
\begin{equation}\rho{a}^4=\text{cte}\end{equation} Con esto entonces para el caso de polvo, la ecuación de Friedmann (\ref{2}) se escribe
\begin{equation}\dot{a}^2+\kappa-\frac{\mc{C}_1}{a}=0\label{7}\end{equation} con ${\mc{C}_1\equiv8\pi\rho{a}^3/3}$ una constante. De manera análoga para el caso de radiación se tiene \begin{equation}\dot{a}^2+\kappa-\frac{\mc{C}_2}{a^2}=0\label{8}\end{equation} con ${\mc{C}_2\equiv8\pi\rho{a}^4/3}$ constante.

Consideremos primero el caso plano ${\kappa=0}$. Para el caso de polvo se tiene
\begin{equation}\left(\frac{da}{d\tau}\right)^2=\frac{\mc{C}_1}{a}\end{equation} cuya solución con ${a(0)=0}$ es
\begin{equation}a=\left(\frac{3}{2}\sqrt{\mc{C}_1}\,\tau\right)^{2/3}\end{equation} mientras que para el caso de radiación, de manera análoga, la solución a la respectiva ecuación es

Factor de escala ${a(\tau)}$ para ${\kappa=0}$

Para el caso esférico ${\kappa=1}$, se tiene, en el caso de polvo ${p=0}$, la ec. de Friedmann,
\begin{equation}\left(\frac{da}{d\tau}\right)^2=\frac{\mc{C}_1}{a}-1\end{equation} es decir,
\begin{equation}\frac{da}{\sqrt{\frac{\mc{C}_1}{a}-1}}=d\tau\end{equation} que puede resolverse introduciendo el tiempo conforme $\eta$ (intento de traducción de conformal time)
\begin{equation}d\eta\equiv\frac{d\tau}{a}\end{equation} cuya interpretación a un tiempo ${t_0}$ es la cantidad de tiempo que le tomaría a un fotón viajar desde nuestra localización hasta la distancia más lejana observable si en ese instante dejara de expandirse el universo en ${t_0}$. Así, se tiene que
\begin{align}
\int{d\eta}&=\int\frac{da}{\sqrt{\mc{C}_1a-a^2}}
\end{align} cuya solución con ${a(\eta=0)=0}$ es
\begin{equation}\eta=\arcsin\left(\frac{2a-\mc{C}_1}{\mc{C}_1}\right)+\frac{1}{2}\pi\end{equation} es decir
\begin{equation}a=\frac{\mc{C}_1}{2}\left(1-\cos\eta\right)\end{equation} y ya que también ${a\,d\eta=d\tau}$, se sigue que con ${\eta(\tau=0)=0}$,
\begin{equation}\tau=\frac{\mc{C}_1}{2}\left(\eta-\sin\eta\right)\end{equation} de modo que se tiene la solución en términos del parámetro de tiempo conforme $\eta$ y puede graficarse con una curva paramétrica. Para el caso de radiación ${p=\rho/3}$, de manera análoga se tienen las soluciones
\begin{align}
a&=\sqrt{\mc{C}_2}\sin\eta\\
\tau&=\sqrt{\mc{C}_2}(1-\cos\eta)
\end{align}

Factor de escala ${a(\tau)}$ para ${\kappa=1}$

Finalmente en el caso hiperbólico ${\kappa=-1}$, de manera análoga se encuentra la solución en términos del parámetro $\eta$, para el caso de polvo ${p=0}$,
\begin{align}
a&=\frac{\mc{C}_1}{2}(\cosh\eta-1)\\
\tau&=\frac{\mc{C}_1}{2}(\sinh\eta-\eta)
\end{align} y para el caso de radiación ${p=\rho/3}$,
\begin{align}
a&=\mc{C}_2\sinh\eta\\
\tau&=\mc{C}_2(\cosh\eta-1)
\end{align}

Factor de escala para k=-1

En conjunto se tiene el siguiente gráfico para ${a(\tau)}$,

Factor de escala ${a(\tau)}$ para polvo y radiación para cada ${\kappa}$

en donde se puede ver cualitativamente el comportamiento de las soluciones para cada $\kappa$. Resalta el hecho de que en ${\kappa=0}$ existen tiempos ${\tau\neq0}$ para los cuales la rapidez de expansión ${\dot{a}}$ cambia de signo y lleva a otro $\tau$ para el cual el factor de escala se anula, llevando al llamado Big Crunch, mientras que para ${\kappa=0,-1}$ el universo sigue expandiéndose para todo $\tau$. De (\ref{7}) y de (\ref{8}), puede graficarse el comportamiento del valor absoluto de la velocidad, que en el caso de ${\kappa=1}$, se tiene

Rapidez de expansión ${|\dot{a}(\tau)|}$ para polvo y radiación para ${\kappa=1}$

mientras que para ${\kappa=0,-1}$,

Rapidez de expansión |\dot{a}|} para polvo y radiación para k=0,-1, donde las curvas punteadas corresponden a k=-1

donde puede verse que $\displaystyle{\lim_{\tau\to\infty}|\dot{a}|=0}$ para ${\kappa=0}$, mientras que $\displaystyle{\lim_{\tau\to\infty}|\dot{a}|=1}$ para ${\kappa=-1}$.

El universo de Friedmann sigue siendo de los mejores candidatos para modelar el universo en el que vivimos, y sin embargo quedan cuestiones que tomar en cuenta como la constante cosmológica, y otras más que entender, como la materia oscura para poder tener un panorama mucho más completo. Para un universo de Friedmann, en general, dada una ecuación de estado termodinámico ${p=p(\rho)}$ (que podría incluir a la constante cosmológica) y las ecuaciones de Friedmann, queda completamente determinada la evolución del universo.

Parámetro de Hubble $H$
La tasa de expansión del universo está caracterizada por el llamado parámetro de Hubble, definido como
\begin{equation}H=\frac{\dot{a}}{a}\end{equation} en honor a Edwin Hubble. Alrededor del año 1920, Hubble descubrió una relación lineal entre las distancias de las galaxias y su corrimiento al rojo debido a su alejamiento relativo, llamado {velocidad de recesión}. Esta relación velocidad-distancia, seguido llamada la ley de Hubble, se expresa como
\begin{equation}\vec{v}=H\vec{r}\end{equation} con $\vec{v}$ la velocidad de recesión y $\vec{r}=a\vec{x}$ la distancia a una galaxia dada, con $\vec{x}$ la coordenada comóvil. Recientemente ([6]) se ha medido el valor de $H$, encontrándose un valor positivo, implicando que en efecto el universo se está expandiendo. El valor de $H$ en el tiempo presente es llamado constante de Hubble y se denota como ${H_0}$.

Parámetro de densidad $\Omega$
La ecuación de Friedmann (\ref{2}) puede escribirse con el parámetro de Hubble como
\begin{equation}H^2=\frac{8\pi{G}}{3}\rho-\frac{\kappa}{a^2}\end{equation} A partir de esta ecuación, se define la densidad crítica ${\rho_c}$ como aquella que vuelve nula la curvatura, i.e.
\begin{equation}\rho_c\equiv\frac{3H^2}{8\pi{G}}\end{equation} de aquí entonces se define el parámetro de densidad $\Omega$ como
\begin{equation}\Omega\equiv\frac{\rho}{\rho_c}\end{equation} y de este modo se reescribe la ecuación de Friedmann como
\begin{equation}H^2=\frac{8\pi{G}}{3}\rho_c\Omega-\frac{\kappa}{a^2}=H^2\Omega-\frac{\kappa}{a^2}\end{equation} y entonces se sigue que
\begin{equation}\Omega=1+\frac{\kappa}{a^2{H}^2}\end{equation} y de este modo si ${0<\Omega<1}$, se tiene curvatura negativa y un universo abierto, si ${\Omega=1}$, la curvatura es nula, el universo es plano y abierto y finalmente si ${\Omega>1}$ la curvatura es positiva y se tiene un universo cerrado.

Tómese en cuenta que se deben considerar los distintos tipos de materia para un valor total de $\Omega$. Recientemente se ha encontrado un valor aproximado de 1 para $\Omega$ total ([6]), por lo que se piensa que actualmente el universo es aproximadamente plano.

Parámetro de desaceleración $q$
El llamado parámetro de desaceleración es una manera de cuantificar la rapidez de la razón de expansión del universo y se define como
\begin{equation}q\equiv-\frac{a\ddot{a}}{\dot{a}^2}\end{equation} Puede verse que esta cantidad surge de la siguiente manera. Considérese una expansión del factor de escala en serie de Taylor alrededor de un tiempo ${t}$,
\begin{equation}a(\tau)=a(t)+\dot{a}(t)(\tau-t)+\frac{1}{2}\ddot{a}(t)(\tau-t)^2+\ldots\end{equation} y dividiendo sobre ${a(t)}$,
\begin{equation}\frac{a(\tau)}{a(t)}=1+H(\tau-t)+\frac{1}{2}qH^2(\tau-t)^2+\ldots\end{equation} donde $H$ es el parámetro de Hubble introducido anteriormente.

Constante Cosmológica $\Lambda$
Un hecho característico de la relatividad general es que relaciona vía las ecuaciones de Einstein el tensor de energía-momento con el campo gravitacional. En gravitación, a diferencia de la física que no incorpora la gravedad, los valores puntuales de la energía importan, no sólo las diferencias de energía de uno a otro estado, lo que abre la posibilidad de introducir una energía de vacío, i.e. una densidad de energía característica del espacio vacío ([5]).

Para construir el tensor de energía-momento del vacío, lo primero que se pediría es que sea un invariante de Lorentz, de modo que no tenga direcciones privilegiadas. La forma de que esto ocurra es que ${T_{\mu\nu}}$ sea proporcional a la métrica,
\begin{equation}T_{\mu\nu}^{(\text{vac})}\equiv-\rho_{\text{vac}}{g}_{\mu\nu}\end{equation} donde mañosamente se ha escrito la densidad de energía ${-\rho_{\text{vac}}}$ (constante), dado que al comparar con (\ref{4}), puede verse que tenemos que el vacío es un fluido perfecto tal que la presión es igual a la densidad de energía con signo opuesto,
\begin{equation}p_{\text{vac}}=-\rho_{\text{vac}}\end{equation} De aquí entonces al considerar el tensor de energía-momento del vacío, las ecuaciones de Einstein se escriben
\begin{equation}G_{\mu\nu}=8\pi{G}\left(T_{\mu\nu}-\rho_{\text{vac}}g_{\mu\nu}\right)\end{equation} que comparando con la ecuación que incluye la histórica metida de pata $\Lambda$ de Einstein (\ref{5}), se sigue que
\begin{equation}\Lambda\equiv8\pi{G}\rho_{\text{vac}}\end{equation} de modo que los términos constante cosmológica y densidad de energía del vacío son intercambiables. En mecánica cuántica, siendo una teoría sin gravedad, surge este tipo de densidad de energía del vacío como la llamada energía de punto cero, que existe aunque no haya partículas presentes. Aunque a veces se trata a la constante cosmológica como nula, como se ha mencionado, ésta es una de las mejores candidatas para representar a la energía oscura y/o la actual expansión acelerada del universo.

[1] Stephen Hawking & George Ellis,The Large Scale Structure of Space-Time, Cambridge University Press, 1975.
[2] José A. Vázquez, Sistemas dinámicos en cosmologías para campo escalar, Tesis de maestría, CINVESTAV, 2007.
[3] Adam G. Riess et. al., Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant, arXiv:astro-ph/9805201 (1998).
[4] Leopoldo García-Colín, Introducción a la Termodinámica Clásica, Trillas, 4a ed., 2008.
[5] Sean M. Carroll, Spacetime and Geometry: An Introduction to General Relativity, Addison Wesley, 2004.
[6] Planck Collaboration, Planck 2013 results. XVI. Cosmological parameters, arXiv:1303.5076 [astro-ph.CO] (2013).

Aniversario luctuoso de Bernhard Riemann

Imagen (recarga la página)
Georg Friedrich Bernhard Riemann nació en un poblado en el Reino de Hanóver, en Alemania el 17 de septiembre de 1826 y murió de tuberculosis en Italia el 20 de julio de 1866, i.e. hace 147 años. Riemann es de los personajes que más admiro, y me he enterado un poco más de su vida por vez primera en el libro Hyperspace de Michio Kaku, que recomiendo bastante, aunque se haga énfasis en solo una parte de su trabajo.

Riemann fue un jóven tímido, introvertido y temeroso de hablar en público de una manera enfermiza. Provenía de una familia pobre, donde era el segundo de seis hijos, cuyo padre era un pastor luterano y veterano de guerra, y la madre había muerto cuando él era niño. Riemann inclusive se dedicó a estudiar ampliamente la biblia, con el plan de convertirse en pastor y apoyar económicamente a su familia, tanto que su padre y él se las arreglaron para que llegara a la Universidad de Göttingen a estudiar teología.

Pero el genio de Riemann había sobresalido desde edades tempranas y la Universidad de Göttingen era un lugar donde las matemáticas fluían como el aire, de modo que pronto Riemann estaba aprendiendo del que probablemente sea el matemático más grande de la historia: Carl Friedrich Gauss. Con apoyo de su padre, Riemann terminó en la Universidad de Berlin, con gente como Peter Gustav Dirichlet, quien se dice sería una gran influencia para Riemann.

Las contribuciones de Riemann fueron bastantes y su influencia en la física es indudable, por ejemplo, con la invención de la geometría riemanniana que es base de la teoría de la relatividad general y probablemente la primer sugerencia de que la realidad física puede estar compuesta de dimensiones extra. Entre las palabras más usadas que llevan su nombre están las integrales de Riemann, las superficies de Riemann, la función zeta-Riemann, la hipótesis de Riemann, etc.

Bernhard Riemann definitivamente fue un matemático adelantado a su época por varias décadas, y hasta hoy se mantiene vigente tanto su legado como varias de sus interrogantes. Es evidente la razón por la que me resulta fascinante la vida de este gran personaje, siendo que en sus condiciones y en poco menos de 40 años lograra cosas tan asombrosas. Puedes leer un poco más en The MacTutor History of Mathematics.

Componentes independientes del tensor de curvatura deRiemann-Christoffel

Sigo trabajando por mi cuenta con temas que me interesan, en este caso tensores. En general me he guiado con "Mathematical Methods for Physicists" de Arfken & Weber para introducirme al cálculo tensorial. El siguiente enunciado corresponde al ejercicio 2.6.6 de la 5a edición.

Se pide demostrar que el número de componentes independientes del tensor de curvatura Riemann-Christoffel (relatividad general), 4-dimensional de grado cuatro ${R_{\rho\sigma\mu\nu}}$, que en general tiene 256 componentes, se reduce consecuentemente a 20 a través de las siguientes relaciones de simetría con índices que corren de 0 a 3,
\begin{align}R_{\rho\sigma\mu\nu}&=-R_{\rho\sigma\nu\mu}=-R_{\sigma\rho\mu\nu}\,\Rightarrow\,36\text{ componentes}\\[0.1in]R_{\rho\sigma\mu\nu}&=R_{\mu\nu\rho\sigma}\,\Rightarrow\,21\text{componentes}\\[0.1in]R_{\rho\sigma\mu\nu}&+R_{\rho\mu\nu\sigma}+R_{\rho\nu\sigma\mu}=0\,\Rightarrow\,20\text{componentes}\end{align} En general, un tensor $\boldsymbol{\chi}$ de grado n en el espacio N-dimensional tiene ${N^n\equiv\xi_{\boldsymbol\chi}}$ componentes, por ello
\begin{equation}\xi_{\mathbf{R}}=4^4=256\end{equation} La primer relación mostrada es de antisimetría. Pensemos en cualquier tensor $\mathbf{A}$ de grado ${n=2}$ tal que
\begin{equation}A_{ij}=-A_{ji}\end{equation} por lo que necesariamente se cumple, para k en el rango de índices,
\begin{equation}A_{kk}=-A_{kk}=0\end{equation} entonces se sigue que el número de componentes independientes de $\mathbf{A}$ se reduce a *
\begin{equation}\xi_{\boldsymbol{A}}=\frac{1}{2}\left[N^2-N\right]=\frac{1}{2}N(N-1)\end{equation} En este caso podemos hacer analogía con el tensor de grado dos tomando parejas de índices como si se tratase de una sola, respectivamente. Existen ${N^2}$ modos de ordenar cada pareja de índices en el espacio N-dimensional **, entonces
\begin{equation}\xi_\mathbf{R}=\left[\frac{1}{2}N(N-1)\right]^2=6^2=36\end{equation} Continuando, la relación que sigue implica la simetría que nos es familiar en una matriz cuadrada o un tensor grado dos
\begin{equation}A_{ij}=A_{ji}\end{equation} pero esta relación implica que en general ${A_{kk}\neq{0}}$, por lo que tendríamos que agregar los N elementos de índices iguales a la relación anterior, entonces tenemos que
\begin{equation}\xi_{\boldsymbol{A}}=\frac{1}{2}N(N-1)+N=\frac{1}{2}N(N+1)\end{equation} Tengamos en cuenta que para la pareja de parejas individuales de índices ${[\rho\sigma][\mu\nu]}$, el tensor es simétrico para la pareja y antisimétrico para cada pareja individual, por tanto
\begin{align}\xi_\mathbf{R}&=\frac{1}{2}\left[\frac{1}{2}N(N-1)\right]\left[\frac{1}{2}N(N-1)+1\right]\nonumber\\[0.1in]&=\frac{1}{8}(N^4-2N^3+3N^2-2N)\nonumber\\[0.1in]&=21\end{align} Para la última relación, nótese que para cualquier tensor $\boldsymbol{\Psi}$, por definición
\begin{equation}0=-\Psi_{\sigma\mu\nu}-\Psi_{\mu\nu\sigma}-\Psi_{\nu\sigma\mu}\end{equation} de donde se sigue precisamente antisimetría, en este caso nula.

De modo análogo al caso anterior, de agrupar términos, aquí se sigue que
\begin{equation}R_{\rho[\sigma\mu\nu]}=0\hspace{0.25in}\Longrightarrow\hspace{0.25in}-R_{[\sigma\mu\nu]\rho}=0\hspace{0.25in}\Longrightarrow\hspace{0.25in}R_{[\rho\sigma\mu\nu]}=0\end{equation} es decir, que la parte antisimétrica de $\mathbf{R}$ [entre corchetes], es nula, lo que significa que se trata de un tensor completamente antisimétrico.

Para relacionar esta propiedad con las dos anteriores, véase que para cualquier tensor ${\boldsymbol\Psi}$ de grado cuatro completamente antisimétrico,
\begin{equation}\Psi_{\alpha\beta\gamma\delta}=-\Psi_{\alpha\beta\delta\gamma}=-\Psi_{\beta\alpha\gamma\delta}=\Psi_{\beta\alpha\delta\gamma}\end{equation} se tiene un tensor antisimétrico para cada índice y simétrico para cada pareja de índices, se sigue entonces que esta propiedad de $\mathbf{R}$ es independiente de las dos anteriores.

Un tensor completamente antisimétrico, por ejemplo, es el símbolo Levi-Civita
\begin{equation}\epsilon_{\alpha\beta\gamma\delta}=\left\{\begin{array}{lll} 0, & \text{si cualquier indice es identico a otro} \\ 1, & \text{si }\alpha\beta\gamma\delta \text{ son permutaciones pares de }0123 \\ -1, & \text{si } \alpha\beta\gamma\delta \text{ son permutaciones impares de }0123 \end{array}\right.\end{equation} y es evidente que sus componentes independientes son exactamente 1.

En general para un tensor completamente antisimétrico de grado n en el espacio N-dimensional, las componentes independientes serán las combinaciones de n en N. ***

Así pues, se concluye que
\begin{align}\xi_\mathbf{R}&=\frac{1}{8}(N^4-2N^3+3N^2-2N)-\binom{N}{n}\nonumber\\[0.1in]&=\frac{1}{8}(N^4-2N^3+3N^2-2N)-\frac{N!}{4}\nonumber\\[0.1in]&=\frac{1}{8}(N^4-2N^3+3N^2-2N)-\frac{(N-3)(N-2)(N-1)N}{4}\nonumber\\[0.1in]&=\frac{1}{12}N^2(N^2-1)=20\end{align} * Pruébense los casos más sencillos N=2,3,4

** ${N=4\,\Longrightarrow\,\{00,01,10,02,20,\ldots,13,31,23,32\}}$ (se le llama regla del producto de pares ordenados)

*** Por ejemplo es claro que si n>N todas las componentes serán nulas, pues siempre se repetirá uno u otro índice. De igual modo, véase que la fórmula obtenida antes es simplemente el caso particular ${\binom{4}{2}=6}$, pero cada pareja puede ordenarse en ${4^2}$ modos, por lo que se tienen ${\binom{4}{2}^2=6^2=36}$ componentes independientes.

La conjetura de Rhie y el Teorema Fundamental del Álgebra

Fue el año 1919 cuando Arthur Eddington confirmaba la desviación de la luz de las estrellas frente al campo gravitatorio del Sol y así entonces también confirmaba uno de los aspectos claves de la relatividad general de Einstein.

La predicción de la relatividad general confirmada por Eddington, es que el camino de un rayo de luz entre un observador y la fuente se verá desviado con la curvatura del espacio-tiempo producida por cualquier campo gravitatorio (planetas, galaxias, u hoyos negros que producen fuertes campos gravitatorios) entre ambos. Una consecuencia de esto es el efecto de lente gravitacional, y el resultado sería que el observador vería la luz desde direcciones distintas, como si hubiera más de una fuente de luz ahí.

Imagen (si ves este texto recarga la pag)


Imagen (si ves este texto recarga la pag)


Imaginemos n masas puntuales interviniendo entre la fuente y el observador con un centro de masa que es precisamente nuestro lente gravitacional y consideremos el caso en que la distancia entre estas masas es pequeña (comparada con la distancia entre las masas y el observador) de modo que podamos modelar esta distribución con un plano, -llamado plano de lente-. La llamada ecuación de lente -sabiendo que trabajamos con un plano y por tanto podemos usar números complejos- proyecta el k-ésimo punto a la posición ${z_k}$ del plano de lente mediante
$$w=z-\sum_{k=1}^n\frac{\sigma_k}{\bar{z}-\bar{z}_k}$$ con ${\sigma_k\in\mathbb{R}}$ relacionada a la k-ésima masa. Un caso particular y espectacular se da cuando n=1, es decir, cuando hay una sola masa y ésta se encuentra directamente enfrente de la fuente (desde el punto de vista del observador), para ${w=0}$ se tiene una circunferencia con centro en la lente gravitacional. Esta circunferencia es llamada Anillo de Einstein.

Imagen (si ves este texto recarga la pag)

Recientemente Dmitry Khavinson y Genevra Neumann decidieron estudiar el teorema fundamental del álgebra -sin idea alguna de astrofísica- para polinomios complejos armónicos (polinomios en 2D, ${P(x,y)}$, que satisfacen la ecuación de Laplace: ${\nabla^2P=0}$), de la forma ${f(z)=p(z)/q(z)-\bar{z}}$, mientras Sun Hong Rhie trabajaba con la ecuación de lente para casos n>1, haciendo pública la conjetura de que el mayor número de imágenes producidas por un n-lente gravitacional es 5n-5, que no demostró rigurosamente, sin embargo Khavinson y Neumann demostraron que para la ecuación
$$r(z)=\sum_{k=1}^n\frac{\sigma_k}{z-z_k}+\bar{w}$$ siendo de la misma forma de la cual partieron sus investigaciones, se tenían precisamente, como cota superior, 5n-5 soluciones, y que el encontrar los ceros en la ecuación anterior significa encontrar las imágenes producidas por el n-lente gravitacional para n>1.

Así, detrás de las ideas que surgieron únicamente de la astrofísica se encontraba una demostración puramente matemática, y detrás del teorema fundamental del álgebra en el caso de polinomios complejos armónicos estaba la respuesta a un caso del efecto de la lente gravitacional.

En la siguiente liga se puede encontrar uno de los artículos de Khavinson y Neumann sobre el tema: From the fundamental theorem of algebra to astrophysics; a "harmonious path"