Quantum mechancis in your android device

mayo 23, 2013

El título es una paráfrasis del libro Quantum Mechanics Using Computer Algebra: Includes Sample Programs in C++, SymbolicC++, Maxima, Maple, and Mathematica .

Esta es la segunda edición. La primera tenía el mismo título, pero cambiaban los entornos en los que se ofrecía software. En particular -de cara al propósito de ésta entrada- en vez de máxima usaba reduce. Cualquiera que haya seguido recientemente el blog adivinará que la parte de máxima esta justificada por la reciente versión de máxima para android. Si menciono la primera edición y reduce es porque resulta que también existe una versión de reduce para android. Aviso, cómo ahí explica hay dos reduce para android, el que enlazo y uno que se llama Android reduce. El primero es un entorno que incluye toda la funcionalidad pero que debe ejecutarse desde la línea de comandos. El segundo es una aplicación visual que permite hacer una serie de operaciones mediante una interfaz visual que, por debajo, llaman al core. Ah, y, por supuesto, en android hay varios compiladores de C++ (de symbolic C++ no lo sé, acabo de enterarme de que existe eso xD) así que también puede usarse ese código con uno de los compiladores.

En definitiva, que en ese libro vamos a tener un montón de código para resolver problemas típicos de mecánica cuántica, desdealgunos muy elementales hasta tópicos mas avanzados como teoría de perturbaciones, scattering, etc. También incluye algunos temas que podrían considerarse mecánica cuántica relativista, cómo las matrices de Dirac. Un fallo de la primera edicion es que no contiene ningún diskette/CD-ROM (es del 1995, por esa época no había nada mas sofisticado, y bueno, tampoco hace falta un DVD o blue ray para meter ese código y si acaso las versiones de los programas), no sé si la segunda versio´n lo incluirá, esperemos que sí.

Por cierto, cómo habréis podido ver en la web es tienen algunos pdf com parte de los primeros capítulos que incluye algo de software. Un código de máxima sacado de ese capítulos es el siguiente: /* wavepacket.mac */
depends (psi,x,psi,t);
depends (f1,x,f1,t);
depends (f2,x,f2,t);
f1: B/(sqrt(1+%i*hb*t/(m*a2)));
f2: exp(-(x*x)/(2*a2*(1+%i*hb*t/(m*a2))));
psi: f1*f2;
res1: hb*hb*diff(psi,x,2)/(2*m);
res2: %i*hb*diff(psi,t);
result1: (res1+res2)/f2;
result2: expand(result1);
print(result2);

He probado a ejecutar ese código en el maxima para windows y el de android y en ambos se ejecuta correctamente. Hay alguna diferencia eso sí. En el de windows se puede hacer copy & paste de todo y darle a ejecutar y funciona, haciendo todos los cálculos de un tirón. En el de Android eso no va y he tenido que ir copiando el código línea por línea, lo cuál es un poco pesado. En la práctica eso no sería problema porque lo suyo es poner este código en un archivo de texto (extensión txt, max o mac) y cargarlo con el comando load( path del archivo). Aviso, lo de las rutas en linux/android son un poco lio par ala gente que está acostumbrada a windows. Si ponéis los scripts de máxima en una carpeta llamada “maximaScripts” de la microSD que viene con el dispositivo (microSD0 cómo suelen llamarla los exploradores de archivos) tenéis que escribir como ruta: “/mnt/microsd/maximascrips.nombrearchivo.mac” (espero no haberme equivocado que lo estoy escribiendo de memoria).

También probé a ejecutar algo de código de reduce e igualmente funciona correctamente en la versión para android. Para cosas tan sencillas (a nivel de proceso) cómo estas el tiempo de ejecución (probado en un note 10.1 y en un note II – ambos con cuatro núcleos a 1.4 y 16 Ghz respectivamente, y 2 GB de RRAM) es despreciable y no se puede notar ninguna diferencia al tiempo que tarda en ejecutarse en un PC (AMD fusion de 6 núcleos a 3.2 GHZ, y 14 GB de RAM). Por supuesto para cálculos complejos la potencia es importante, pero al fin y al cabo la primera versión del libro se escribió para ordenadores 486, o si acaso Pentium, con velocidades de reloj del orden de los 100 MGHZ y con una memoria que podría ir entre uno 8 y 32 MB de RAM. Vale que los micros de ARM para móviles no son tan potentes cómo sus correspondientes actuales de PC, pero sin duda le dan mil vueltas a los ordenadores de esa época.

Aparte de esté código de mecánica cuántica también testeé el maxima de android con algo de cálculo tensorial, en particular le pedí que me sacase el tensor de Ricci a partir de una métrica. Una vez más el cálculo lo hizo de manera casi instantánea y no hubo diferencia apreciable al tiempo que requiere el mismo cálculo en el PC de sobremesa. Eso sí, hacer el mismo cálculo a mano podría requerir horas, así que ahí se ve lo fantástico de tener este tipo de potencia de cálculo en un dispositivo que puedas llevar a cualquier lado, en particular porque éste tipo de cálculos son los que suele necesitar un físico teórico. Incluso si alguien trabaja en teoría de cuerdas se va encontrar con la necesidad de hacer cálculo tensorial, no digamos ya si trabaja en agujero negros en teoría de cuerdas. Con estos packages de cálculo, los libros y artículos que puedes llevar en el tablet, y un stylus para escribir prácticamente hay todo lo que se puede necesitar para hacer física teórica cómodamente. Por supuesto en un tablet windows 8 con stylus aparte de estos programas tienes mathematica y maple, que aún no están portados para android, y que tienen mas scripts en la red para más tareas, pero pese a ello uno se puede apañar mas que dignamente con lo que ya hay. En el peor de los casos si uno quiere un mathemática puede pagar por ese invento extraño del wolphram alpha pro que da una funcionalidad muy similar a mthematica, pero en versión en la nube (que es un tipo de informática que detesto, y no digamos ya de pago, as´ique conmigo que no cuenten xD).

Bien, esto de poder comprar un libro es fantástico, pero sería interesante tener -cómo anticipé- material para descargar desde la web. Hay varias páginas. Por ahora la más interseante que he encontrado (par aun teórico) es MAXIMA BY EXAMPLE: DETAILED EXAMPLES OF THE USE OF THE MAXIMA COMPUTER ALGEBRA SYSTEM escrita por el físico Edwin L. (Ted) Woollett.

Por supuesto si uno sube de nivel sería interesante tener software para hacer cálculo de cosas cómo el espectro de una teoría supersimétrica a partir de su Lagrangiano, o algo para hacer diagramas de Feynman, o vaya usted a saber qué. Cómo dije antes mucho de ese software está para mathematica y mapple, que por ahora no están disponibles en Android. Pero si hay mucho para matlab, y si está para matlab está para Octave, que si está disponible en android. En fín, no nos vamos a engañar, para windows sigue habiendo mas material, pero comparado a lo que había hace un año ahora tenemos Maxima, Octave, mathstudio y reduce (por citar los mas importantes) y teniendo en cuenta que no hay móviles con windows 8 y sí con android y que lo mejor es poder compartir todo el software en el móvil, en el tablet y en el ordenador, y que eso sólo puede hacerse con Android, pues está claro cuál es el camino. Bueno, en el futuro próximo los de Ubuntu van a sacar una versión de Linux para móviles y tablets, y par linux ya sí que hay prácticamente lo mismo que para windows.

Entiendo que hay mucha gente que tiene la idea de que lo suyo es dedicarse a la física en un despacho, con el ordenador al lado. O si acaso en una biblioteca, con un portátil. Pero eso limita mucho, yo quiero disponer de todo el material posble en cualquier lado, porque, a ver, ¿quien no ha quedado con amigos, estos se van a la discoteca, se pierden por ahí, y uno se queda aburrido (eso puede ocurrir incluso si los amigos no se pierden xD) sin nada interesante que hacer? Pues en esos caso se saca del bolsillo el note II y se pone a hacer física y si uno necesita calcular un tensor de Ricci que mejor que tener un software que lo haga y luego ya hacer a mano los cálculos que corresponda ;) .

Por supuesto, si se es del grupo de gente alérgico a la tecnología siempre puede optar por opciones mas…clásicas xD.

Caos espectral en mecánica cuántica (o perturbación crítica)

mayo 15, 2013

Actualización: He modificado el título de la entrada para llamar a la teoría propuesta por un nombre mas descriptivo, y posiblemente mas llamativo también ;) .

Hace unas fechas escribí una breve reseña sobre caos cuántico. La ligereza de la misma puede haber servido cómo pista de que no era un tema en el que hubiese profundizado demasiado, por decirlo suavemente. El motivo es que no era un tema que me hubiese resultado especialmente atractivo hasta ahora. Cómo dije ahí esa gente tiente a ocuparse de mecánica cuántica muy sencilla y -además- en última instancia no existe un análogo cuántico del caos clásico así que no terminaba de ver un buen motivo para profundizar en ese campo.

Con los precedentes anteriores uno podría plantearse cómo es que ahora dedico una segunda entrada al tema. El motivo proviene de algo aparentemente muy alejado de la mecánica cuántica: el cálculo numérico. Cómo parte de la titulación de matemáticas es necesario cursar al menos dos asignaturas sobre cálculo numérico. En el primer año se ven varios temas relativamente sencillos cómo cálculo de ceros de funciones, interpolación de funciones, integración/derivación numérica, resolución de sistemas lineales y, la parte mas extensa, álgebra lineal numérica.

En los cursos de álgebra lineal se aprende a calcular autovalores y autovectores de manera analítica. Cómo suele ser habitual el método que es eficiente para hacer algo con “lápiz y papel” dista mucho de ser el mejor para hacerlo mediante algoritmos en un ordenador. En el caso del álgebra numérica sucede lo mismo, pero, además, se aprenden algunos hechos tan interesantes cómo sorprendentes.

Recordemos que para calcular los autovalores de una matriz M uno debe hacer el determinante de \mid M -\lambda \mathbb{I} \mid . Esto da cómo resultado un polinomio en \lambda . Una vez se tienen las raíces de ese polinomio se sustituyen en la expresión \mid M -\lambda \mathbb{I} \mid y se calculan el/los autovectores correspondientes a ese autovalor. En los ejercicios típicos las matrices están “cocinadas” para que el polinomio tenga raíces enteras (o al menos alguna de ellas entera y que el resto puedan obtenerse a partir de ahí por el método de Rufini). Por supuesto en problemas reales no sucede eso casi nunca (no estoy 100% seguro pero creo recordar que en el conjunto de todos los polinomios los que cumplen eso tendrían medida nula -definiendo una medida más o menos natural en el espacio de polinomios, claro-). Por ese motivo, en la práctica, uno debería resolver el polinomio por métodos numéricos y luego implementar un algoritmo que calculase el autovector a partir de ese autovalor.

El caso es que esa tarea tan sencilla de calcular las raíces de un polinomio tiene sutilezas inesperadas. La clave del asunto es que polinomios muy similares pueden tener raíces no tan similares. Es decir, una pequeña incertidumbre en el valor de un coeficiente del polinomio se traduce en una gran diferencia entre los valores de las raíces, osea, cómo exclamaría Malcom: ¡CAOS!. Antes de seguir con consideraciones teóricas sobre esto dejo un ejemplo de matriz cuyo polinomio característico tiene esa sensibilidad:

A=\left(\begin{array}{ccc}    -149&-50&-154\\    537&180&546\\    -27&-9&-5    \end{array}\right)

Uno puede hacer el polinomio característico y verificar que sus raíces son {1,2,3}. Si uno modifica muy ligeramente esa matriz, por ejemplo modificando el segundo elemento de la diagonal de 180 a 180.01 podría comprobar (recomiendo usar algún programa informático) que las nuevas raíces el polinomio característico (vamos, los autovalores de la matriz), son { 0.207, 2.3008, 3.50} Es decir, una modificación de uno de los elmentos de la matriz del orden de 10^{ -5} modifica todos los autovalores en magnitudes entorno al 50% lo cuál es algo realmente impresionante. En los libros o manuales elementales sobre cálculo numérico no se suele comentar mucho más al respecto y se pasa directamente a enseñar métodos para el cálculo de esos autovalores (según en que manuales se limita al método de potencias para el cálculo del autovalor dominante) y autovectores (normalmente el método QR y variantes). Tampoco suelen hacer las cuentas de cómo cambian los autovectores así que me hice el cálculo para los dos casos anteriores. Para la matriz original los autovectores también son bastante distintos, con variaciones incluso mayores que las de los autovalores.

El caso es que para alguien con una base de física una matriz es, sobre todo, un operador cuántico (para un matemático una matriz puede ser un montón de cosas, y dependiendo del caso se la estudia de muchas maneras diferentes xD), los autovalores las autoenergías (si el operador es el hamiltoniano) y los autovectores las autofunciones de onda cuánticas. Y claro, inmediatamente (al menos yo es lo que pensé nada mas leer ese resultado) es que si un operador cuántico tiene un comportamiento tan exótico uno podría pensar que algo extraño podría pasar con Mary…digo la cuántica ;) . El caso es que no pude dedicarme inmediatamente a profundizar en ese hecho chocante, pero siempre estuvo ahí en segundo plano, cómo una inquietud, y en cuanto se dió la oportunidad analicé más a fondo el asunto. La primera duda que surge es sí esa incertidumbre en el cálculo de autovalores es debida a errores de redondeo en algoritmos numéricos o si obedece a una causa mas fundamental. Rastrear en extensos libros sobre álgebra lineal numérica no me llevó a ningún lado pero una búsqueda en google me llevó a un foro dónde se trataban esos temas y ahí daban un ejemplo muy sencillo que -sí no interpreto mal- resuelve la cuestión. En concreto plantean el caso del poinomio \lambda^3 -\epsilon=0 cuya solución \lambda=\sqrt[3] \epsilon que no es derivable en el entorno de 0 y ese es el origen de la sensibilidad del polinomio. En la misma web mencionaban un ejemplo mas complejo de una matriz, dependiente de un parámetro, que originaba polinomios “sensibles” a variaciones de ese parámetro. Esa matriz era además simétrica (autoadjunta) lo cuál es bueno pues los operadores cuánticos deben ser autoadjuntos. Con eso ya se tiene bastante información relevante, el problema es “fundamental” y no de redondeo, se identifica el problema (o al menos un factor del mismo) y se pueden analizar familias de matrices, no una sóla.

La siguiente reflexión que a uno se le ocurre es plantearse cómo puede suceder esto con sistemas lineales. Al fin y al cabo la cuántica es lineal, y uno aprende, estudiando Sturn-Liouville (bien sea mediante análisis clásico o, de modo riguroso, en cálculo funcional) que la clase de operadores lineales autoadjuntos son buenos y maravillosos y nos dan una base del espacio de Hilbert de soluciones del problema ¿que más se puede pedir?. Bien, el caso es que bajo esa aparente inocencia los operadores lineales (incluidos los autoadjuntos) ocultan muchas sorpresas y uno, a poco que los estudie, se da cuenta de que son unos grandes desconocidos que guardan en su interior muchas pautas insospechadas. Pero, volviendo al principio, si la cuántica es lineal ¿de dónde surgen todas estas “sensibilidades”? Bien, los operadores son lineales, sí, pero las operaciones para extraer información de ellos (sus autovalores y autovectores en el caso de la cuántica) implican formar expresiones no lineales. La tarea de obtener el determinante es no lineal, y para calcular autovalores el resultado de la misma es un polinomio no lineal. Digamos que hay mucha no-linealidad escondida.

Bien, esa es la matemática, pero queda analizar un poco como afecta eso a la física. La idea básica -en un esquema meramente formal- es muy simple. Tenemos un sistema cuántico, todo lo particular que haga falta, cuyo hamiltoniano podemos considerar que es, en alguna base apropiada, una matriz finita. Los elementos de esa matriz en general van contener términos que se deben obtener de manera experimental (por ejemplo si es el hamiltoniano de un electrón en un campo eléctrico el valor de ese campo podría ser un dato experimental). Entonces eso significa que tenemos dos Hamiltonianos, H y H’ que difieren por una pequeña cantidad. Siguiendo la costumbre de teoría cuántica de perturbaciones podríamos escribir H'=H_0 + \epsilon H_1 aunque, en este caso, no nos importa (necesariamente) que H_0 sea resoluble analíticamente. La idea es que uno podría esperar que los autovalores de H y H’ fuesen muy similares (es el fundamento de la teoría de perturbaciones, en particular ahí se exige, cómo prueba de consistencia, que la diferencia entre un autovalor del sistema sin perturbar y el perturbado sea menor que la diferencia entre dos autovalores del sistema sin perturbar). Pero, cómo acabamos de ver, esto no siempre tiene porque suceder. Yo estudié esto por mi cuenta y elaboré un poco algunas consecuencias sencillas. Más adelante descubrí que hay una línea muy reciente de investigación, liderada por Michel Berry (el de la famosa fase de Berry) y llaman a esto “perturbación crítica”. Aún tengo que explorar mas el tema de lo que hace esa gente y cuanto se parece a lo que yo estoy considerando.

Pero sigamos con el quid de la cuestión. La idea es que dos sistemas con hamiltonianos muy similares pueden tener energías muy diferentes. Podría darse el caso medir el campo eléctrico del hamiltoniano con precisión de varios decimales y que pudiésemos resolver el problema y aún así los resultados no nos servirían para predecir, en la práctica los valores posibles de las energías. Pero puede ser peor aún, cómo los autovectores, que son las funciones de onda, también cambian mucho. Imaginemos que el campo eléctrico fluctua en el tiempo. Si colocásemos el sistema en un estado inicial de superposición y midiéramos las frecuencias con las que se da cada autovalor de la energía estas no tendrían una distribución probababilística. El motivo es que al fluctuar el campo fluctuan los autovalores y no siempre estamos trabajando con autoenergías similares. Y, cómo además varían los autovalores, las probabilidades de ocupar cada autovalor también fluctúan. Es decir, podríamos tener un sistema del que sabemos el Hamiltoniano con mucha precisión, poder resolverlo analíticamente, y aún así, en la práctica, no poder obtener ninguna información útil respecto a que nos vamos a encontrar.

Para tratar esos sistemas- creo yo, habría que optar por una descripción en términos de ecuaciones diferenciales estocásticas (para una introducción ver por ejemplo este pdf) en la que aparte del término determinsta (la ecuación de Schröedinger) habría un término de “ruido”. Eso sí, ese término no tendría porque ser browniano sino que su naturaleza dependería de la naturaleza analítica del parámetro del hamiltoniano que dicta la “sensibilidad” del mismo y podría bautizarse algo así cómo “ruido espectral”.

Para ir concluyendo hago una reflexión importante. Esto no es caos cuántico. En un sistema caótico clásico tenemos que la dinámica (el hamiltoniano) es fijo y hay sensibilidad en las condiciones iniciales (que no pueden medirse con precisión infinita). En cuántica el observable fundamental es la función de onda y no las posiciones/momentos. Y la unitariedad de la evolución cuántica implica que si las funciones de onda en un instante dado difieren por una cantidad pequeña esa diferencia se mantendrá constante en el tiempo. Esto plantea una duda conceptual de cómo si la realidad es cuántica en sistemas clásicos, que son el límite de los cuánticos (teorema de Erenfest) puede haber caos. Por supuesto mi argumento no implica que haya caos cuántico porque aquí lo que tenemos es algo distinto. Tenemos que el propio hamiltoniano (lo que dicta la dinámica) es el que está sujeto a una incertidumbre experimental y cómo consecuencia de la misma los observables cuánticos (autovalores, autofunciones) son muy sensibles a variaciones de esa incertidumbre. Por supuesto, y esto sería curioso de analizar en comparativas, en sistemas clásicos también hay esa incertidumbre en el valor exacto del hamiltoniano, y también hay operaciones de obtener autovalores y autovectores para obtener soluciones en algunos de esos sistemas (por ejemplo osciladores armónicos acoplados) así que esta sensiblidad extra, esta “perturbación crítica” afectaría por igual al mundo clásico y al cuántico, y tal vez (o tal vez no, vaya usted a saber xD) seria interesante comparar las diferencias entre ambos mundos para esos sistemas.

Para finalizar algunas palabras sobre lo que hace la gente de caos cuántico, que está relacionada con la naturaleza de los autovalores de los operadores audoadjuntos. Resulta, por ejemplo, que los hamiltonianos que presentan simetrías tienen una distribución de los autovalores muy diferente de los que no tienen simetrías (si se quieren buscar detalles usar los términos “quantum chaos, random matrix). Es un tema curioso, sobre el que tal vez lea mas, o tal vez no. Pero tras ver esto de la sensibilidad de los autovalores y que la distribución de los mismos depende de las simetrías del hamiltoniano está claro que bajo su inocente apariencia los operadores autoadjuntos tal vez puedan ocultar auténticos “animales patológicos” en su interior y que posiblemente la mecánica cuántica mas elemental guarde aún muchas sorpresas importantes en contradicción con la idea de que es un “animal doméstico y conocido”. Y eso si nos restringimos a cuántica elemental, y estudiando matrices finitas (a saber que pasa con las infinitas que son lo común en mecánica cuántica). Pero el caso es que las teorías cuánticas de campos también son, en el fondo, teorías cuánticas “normales” y, no sé ¿cómo podría ser el grupo de renormalización de un hamiltoniano de campos que fuera el análogo de uno de partículas “sensible”? ¿Tal vez el flujo de renormalización hiciese evolucionar el valor de las constantes de acoplo de manera caótica según nos movemos hacia energías mas altas? O, si la estructura de los autovalores depende de la simetría ¿que pasa en los fenómenos de ruptura espontánea de simetría?

En definitiva, que me da la impresión de que el formalismo de la mecánica cuántica convencional puede ser mucho, mucho mas rico de lo esperado, y que, por ejemplo, tal vez algunos fenómenos que por argumentos de “naturalidad” podrían parecer muy improbables a lo mejor no lo sean debido a que en algún punto hay alguna “sensiblidad” oculta en algún punto.

Preones, una buena idea con muchos problemas

mayo 5, 2013

Actualización: Veo que nadie ha hecho ninguna propuesta en la línea que pido en el último párrafo. Dejo esto unos días mas para ver si alguien se anima y escribiré la actualización después.

Antes de ir con la entrada en sí voy a hacer un comentario general sobre el panorama de la blogosfera de la física en España. De la lista de blogs que tengo enlazados me encuentro con que hay varios blogs que empezaron bien pero llevan años sin actualizar. Otros que han tenido un ritmo de publicación elevado están muy parados. Kanijo, por ejemplo, lleva sin actualizar desde Febrero y el blog de Migui parece haber sido abducido por alguno de los magufos que tanto le critican. Otro blog magnífico, el de física en la ciencia ficción, directamente ha sido borrado (y no me consta que tenga sustituto) . Sólo Francis continúa infatigable. Me pregunto cuál será el motivo -si hay alguno común- o si es sólo alguna coincidencia. A nivel particular debo decir que he disminuido mucho mi ritmo de publicación por lo que podríamos llamar estrés post traumático tras habernos encontrado con que el LHC nos ha traído el peor escenario posible.

Bien, vamos con el tema de los preones (no confundir con los priones, las proteínas con pliegue anómalo que están asociadas a la enfermedad de Creutzfeldt-Jakob, o síndrome de las vacas locas) , que, además, está bastante relacionado con el LHC. Debo decir que había oído hablar varias veces de los preones en diversos sitios. Hablando de manera rápida los preones serían unas partículas hipotéticas mas pequeñas que los quarks y los leptones y tales que ambos (quarks y leptones) estarían formados de estos preones. La referencia mas extensa sobre los mismos la leí en el blog de Lubos( Preons probably can’t exist ) en el que argumenta porque son una idea muy simple y porque hay motivos para pensar que presenta bastantes inconsistencias.Por cierto, esa entrada ha sido editada desde que yo la leí y le ha añadido bastantes fórmulas (el blog de lubos -que siempre ha sido muy bueno- ha mejorado mucho desde que tiene el sistema de Mathajax para escribir fórmulas) Aparte de esa entrada ha publicado en otras la ausencia de indicios de los mismos en el LHC (el peor escenario abarca también a los preones).

Bien, cómo quiera que ahora mismo se me hace tarde, y que el post de Lubos está muy bien y creo que todo el mundo debería leerla dejo a medias esta entrada y la editaré más adelante. Mientras tanto los que hayan leído la entrada de Lubos pueden probar a ver si adivinan cómo se podría (al menos en principio) mantener viva la idea de los preones pese a los sólidos argumentos en contra de Lubos. Si nadie da con nada yo comentaré algunos de los que hay por ahí propuestos al respecto, tanto los medio sensatos cómo los propios de la santísima trinidad ;-) .

Software científico en Android: Máxima y Octave

mayo 1, 2013

En su momento, cuando hablé de la saga note (phablet y tablet) y del HTC flyer, todos equipos con sistema operativo android, mencioné algún software científico que podía usarse en esos dispositivos.

Resalté programas cómo mathstudio, que tiene una gran cantidad de funciones incorporadas y un lenguaje de script. Es un buen programa que, desafortunadamente, parece haber paralizado su desarrollo y lleva congelado en la versión 5.3 alrededor de un año. También comenté la existencia de programas cómo Addi y mathmatiz que eran unos clones restringidos de matlab. Posteriormente llegó a android una versión de Sage que permite conectar con el servidor web a través de una interfaz dedicada. También llegó matlab, pero. lo mismo que Sage, en una versión “para la nube” que permite conectar bien con los servidores del fabricante o bien con una versión con una versión de matlab en un PC local. Habrá a quien esto de la nube le parezca bien pero para mí en particular me parece un invento inútil.

Antes de seguir con el software una pequeña consideración sobre hardware y sistemas operativos. A finales de verano apareció windows 8. Yo había probado en su momento la versión beta en un tablet archos que inicialmente venía con windows 7. En la época en que tuve el archos android estaba casi recién nacido y había muy poco software serio para ese sistema. Además, el hardware no era muy potente. Pero ahora la situación a nivel hardware ha cambiado. El note II y el note 10.1 tienen un procesador con cuatro núcleos y 2 GB de RAM. Ambos traen 16 GB de almacenamiento interno (yo quería haber pillado las versiones de 32 o 64 GB, pero no es fácil encontrarlas) y una microSD de 32 y 64 GB respectivamente en cada uno. Cierto que esos micros no son tan potentes y la memoria no es tan veloz cómo la de PC, pero aún así siguen siendo unas especificaciones bastante potentes. Mi último – y único, que nunca me han entusiasmado esos chismes- portátil (que tengo prestado con carácter indefinido) era un dual core a 1.8 GHZ con 2 GB de RAM y 128 GB de disco duro. Debo decir, que, en general la sensación de uso con los note es mucho mejor que con el portátil. Y dado el mayor tamaño de los programas windows el almacenamiento es equivalente.

Vale, en el último párrafo he comparado un hardware de éste año (o finales del pasado) en android con otro de hace unos 6 años en windows, tal vez es injusto, aunque, todo sea dicho, con ese portátil podía usar todo el sotware matemático que necesitaba y buena parte del musical (que es el software mas exigente que suelo usar en ordenadores). Ahora hay tablets windows 8(versión pro – y RT para arquitecturas ARM, aunque no le veo yo ningún futuro a ésta versión-) cómo el surface de microsoft y Samsung Ativ. Probé este último en un corte inglés y no me convenció nada. Era demasiado grande cómo para llevarlo cómodamente por ahí y la suite del S-pen no estaba a la altura de la de Android (aunque haya alternativas, claro). La potencia y fluidez no estaban mal, pero, así, en general, el windows 8 me parece un híbrido un tanto fallido. Las aplicaciones de escritorio -que son las potentes- apenas pueden optimizarse para uso táctil y las hechas para la interfaz táctil metro tienen demasiadas limitaciones (casi tantas cómo las de IOS y desde luego mucho más que las de Android) así que no confío yo en ver ninguna especialmente estupenda. Y, desde luego, ahora mismo no hay ningún software de los importantes portados a esa interfaz. Total, que al menos por ahora el camino en talbets/móviles es android, y por eso esta entrada.

Y vamos ya por fín con los dos programas protagonistas. Seguiré el orden cronológico de aparición, que es el contrario del de la entrada.

Octave

Para quien no lo conozca decir que octave es un clon muy, muy completo de matlab para linux (aunque por supuesto hay versión para PC, y ahora android). En algunas ingenierías se usa éste en las clases en vez de matlab así que su uso está bastante extendido.

El responsable del port es el autor del anteriormente mencionado programa Addi. Para instalar el programa hace falta instalar también el emulador de terminal para android, que es desde dónde se va a teclear el código de octave. Para poder usar gráficos hay que instalar otro programa, droidplot. Y aparte de los packages que vienen por defecto pueden instalarse otros adicionales. En ésta descripción sobre el modo en que se instala Octave se ve una gran ventaja de los programas Android sobre los que hay para windows 8 (interfaz metro) y, desde luego, IOS. Hablo de la posibilidad de que los programas se comuniquen unos con otros, lo cuál permite crear una infraestructura de algo similar a los plugins típicos de windows (y MAC y creo que también en Linux, aunque mi experiencia con linux es limitada), que son la verdadera potencia a la mayoría de sus programas mas famosos.

Coincidió que cuando instalé octave -allá por navidades- tuve que dar clases de octave (apoyándome en mi conocimiento de matlab obtenido en las asignaturas de cálculo numérico) a un chaval, estudiante de ingeniería, que lo usaba en linux. Por la misma época también estuve ayudando a alguien a hacer un modelo sobre temas de ecología matemática. Para estudiar ese modelo hubo que elaborar un programa en matlab. El caso es que probé el código matlab de las clases y del programa en la versión de octave para android y fuí capaz de ejecutarlo (una vez me hice con el modo de usar el programa) en android. Un poco mas adelatne me tocó hacer un programa en matlab (basándome en algo ya hecho) para un tema mucho mas complejo (ecuaciones en derivadas parciales no lineales con valores de frontera) y una vez más el código que iba en matlab iba en octave.

También probé algún código típico de cálculo numérico (runge kutas y cosas así) y no me dió problemas. Cómo no uso matlab (ni hago cálculo numérico habitualmente) tengo claro que par mis necesidades en ese campo octave cumple de sobra, aunque la interfaz de matlab sea mas amigable.

Máxima on Android

Máxima es un software para linux, con versión para PC, orientado al cálculo simbólico. Sólo lo conocía hasta ahora por referencia. Me enteré de su existencia cuando buscaba software en android que pudiese hacer cálculo tensorial (en una época en que- además- necesitaba hacer ese tipo de cálculos sí o sí y cómo no encontré alternativa al final hice parte de ellos a mano y la otra parte en el mathemática del PC. con lo cuál se perdía el “do everything everywhere” que es la esencia de los tablets). En la wiki sobre software científico informaban de la lista de programas que hacían cálculo simbólico con tensores y aparecía mencionado máxima, pero el problema es que por ese entonces no estaba disponible para Android. Hoy, en una de esas búsquedas que hago periódicamente, he descubierto que al fín lo han portado. Y, claro, lo he instalado ipso facto. Ocupa bastante, 90 MB (30 el programa base y 60 de librerías extra que se bajan posteriormente a instalar la aplicación). Cómo no lo conocía previamente y acabo de instalarlo no puedo decir mucho sobre él. He hecho las pruebas básicas (usando el teclado de sagemath para álgebra, que es bastante cómodo para introducir expresiones matemáticas, aunque hay otros teclados similares e iré probando cuál es mejor). Veo que la sintáxis no es la misma que la de mathemática en algunos aspectos (por ejemplo, los parámetros de las funciones van entre paréntesis en vez de entre corchetes) aunque sospecho que no es demasiado distinta. Mirando el manual veo que trae soporte integrado para cálculo tensorial y no depende de packages externos, cómo mathemática, lo cuál parece una ventaja. He leído, buscando información, que en algunas universidades en vez de mathematica o maple se usa este programa en algunas asignaturas así que asumo que debe ser bastante potente. En todo caso está claro que tanto octave cómo éste máxima son programas profesionales- y gratuitos- que pueden hacer un trabajo cuanto menos similar a matlab y mathematica. Cómo quiera que están disponibles en tablets android – y matlab y mathematica no- sospecho que van a tener un boom. Yo desde luego voy a instalar también la versión windows de máxima y en la medida de lo posible haré todo en este entorno para así tener compatibilidad entre tablet y PC. Por cierto, Octave está disponible para IOS sólo en versión “en la nube” (vamos, cómo si no existiese por tanto xD) y desde luego máxima no está disponible (ni creo que lo llegue a estar). Los usuarios de la manzana -hasta dónde yo sé- deben conformarse con mathstudio (ok no está mal) y mathomatic. De éste último hhay versión también para android, pero no me he molestado en instalarla porque no parece ofrecer nada que no tengan ninguno de los programas anteriores.

Aparte de programas de corte tan profesional cómo Octave y Máxima, y el resto de mencionados, en android hay bastante mas software científico de diverso pelaje (por ejemplo uno específico para estadística, statistic calculator). Desde luego hay mucho más para windows, pero la brecha se va reduciendo. Dado que el 9x de los móviles a día de hoy son Android y en tablets ya se está en torno al 45% sin duda es de esperar que la situación siga mejorando. Por ejemplo, hay otro programa “famosillo” que también ha sido portado a Android, reduce. Yo no lo he usado nunca, pero sé de su existencia por un libro titulado “quantum mechanics thorught computers” que trae mucho software de mecánica cuántica en diversos entornos, entre ellos éste reduce.

Y no sólo a nivel software mejora Android, a nivel hardware también. El note III se rumorea que vendrá con 3 GB de RAM y con un mínimo de 32 GB de almacenamiento interno. Las microSD de 128 GB se supone que ya están a la venta (aunque aún están muy caras y no he intentando pillar ninguna) y se rumorea que las de 256 no tardarán en salir. Otro aspecto interesante, ya fuera del ámbito científico, es que android ya tiene soporte para teclados midi (vía el USB on the go). Yo lo he probado con el note 10.1, con el programa caustic, y funciona bien. El problema es que la latencia aún es algo alta (unos 80 ms) -al menos en le note- pese a que se suponía que con la versión jelly bean iba a mejorar.

Ciertamente IOS y sus drivers coreaudio siguen estando por encima de Android en asuntos musicales, y, en consecuencia, hay mucho mas software. Pero no podía cerrar esta entrada sin mencionar que uno de los programas mas famosos de la informática musical, Fruity Loops (vale que no es de los realmente grandes cómo cubase, cakewalk, abbleton live o reason, pero sigue siendo uno muy usado por mucha gente)también acaba de ser portado para Android. Aún no lo he probado, pero no creo que tarde mucho en hacerlo ;) .

La singularidad tecnológica ya está aquí

abril 25, 2013

Al menos los primeros pasos. Estamos en los albores de la llegada al gran público de una serie de avances que pueden llegar a replantear muchos aspectos de la vida cotidiana. Voy a analizar algunos de ellos y algunas posibles consecuencias en la vida cotidiana que pueden plantear.

Uno de ellos, posiblemente uno de los que más publicidad tendrá, son las google glass (y sus variantes de otros fabricantes). Para presentarlas dejo un vídeo promocional:

En relación a estas google glass recomiendo a la gente que intente ver el tercer episodio de la primera temporada de la serie de TV “black mirror”. Ahí plantea una sociedad dónde la gente lleva un chip de memoria que graba todo lo que ve y oye (y posiblemente más aspectos, pero eso no se explicita en la serie).En la serie se plantea las consecuencias de eso para los problemas de pareja, en particular las infidelidades, que quedan grabadas.

Bien, las google glass hacen fotos, e imagino que podrán grabar vídeo (posiblemente puedan conectarse por cable con un móvil que vaya en el bolsillo para tener mas espacio de almacenamiento). Por ley los dispositivos de grabación deben avisar cuando hacen fotos (con un sonido) o cuando están grabando (y los mecanismos de autofocus y similares contribuyen a ello en todo caso). Pero, desde luego, hay por ahí patches para saltarse esas limitaciones. Con un dispositivo al uso uno puede intuir cuando le están grabando, pero con las glass no hay manera. Total, que al que no le guste que le hagan fotos (cómo es mi caso) lo lleva crudo.

Bien, eso es uno de los muchos aspectos de esas gafas (que no señalan en el vídeo promocional) a considerar. Hay otro que me resulta muy interesante. A poco que tengan conexiones adecuadas se podría adjuntar a esas gafas una pequeña “webcam”. En particular se podría colocar ea webcam en la parte trasera de la cabeza y podría activarse siempre que se quisiera para que la imagen de la webcam se proyectase en las google glass. Eso significaría que se podría ver lo que está delante de uno y lo que está detrás. ¿Es eso útil? En principio podría parecer que sí, que es muy útil, pero a lo mejor luego uno se pone a ello y resulta que tampoco hace tanta falta. Lo interesante es que la posibilidad estará ahí dentro de nada y podremos comprobarlo. Supongo que a alguna gente le será útil y a otra no, pero ya digo, lo veremos.

Una variante de las gafas y la webcam vendría con los dispositivos voladores pequeños de diverso pelaje, cómo por ejemplo éste.

Por ahora están diseñados para ser controlados desde un móvil o un tablet (ok, o desde un portatil, pero eso ya no es para llevar a todos lados). Pero sin duda podrán controlarse en algún momento desde las google glass. Esto nos lleva a que por las gafas podamos controlar y ver lo que hace el bochocopter. Creo que pocas soluciones mejores que esa hay para poder tener una sensación de primera mano de lo que es volar (al menos mientras no haya algún dispositivo de antigravedad xD).

Bien, esto de las gafas está bien, pero es sólo el inicio. Vamos con un paso más, por ahora me refiero a tecnología destinada a uso comercial. Estoy hablando de un proyecto de Samsung para controlar “directamente con la mente” sus dispositivos. Aquí tenéis una breve reseña del invento. Usar eso para manejar un tablet no tiene mucha historia (salvo para gente con serios problemas de movilidad cómo tetraplejicos o gente con los brazos amputados que no dispongan de prótesis robóticas, de las que luego hablaré). Pero quien maneja un tablet puede manejar otro dispositivo. Por ejemplo, el note 10.1 tiene un sensor de control remoto. Con ese sensor universal (vía software) se puede controlar cualquier dispositivo con interfaz de infrarojos (televisor, cadenas de sonido, etc, etc). Vale, con las google glass se podría hacer lo mismo si las gafas llevasen sensor de infrarojos (y si las google glass no lo llevan ya saldrán las “Samsung glass” que seguramente lo lleven).

El caso es que con las gafas y el sensor éste de ondas cerebrales podremos tener una interfaz con la tecnología que no dependerá de las manos. El abanico de posibilidades es obviamente inmenso, y a saber dónde lleva. Por supuesto para discapacitados será una bendición, pero para gente normal representa posibles ventajas y posibles peligros (por ejemplo, aún más sedentarismo).

La mayoría de la gente recordará cómo en la segunda película de star wars (el imperio contraataca), a Luke Skywalker le corta la mano Darth Vader, pero no es un gran problema porque enseguida se la sustituyen por una artificial. Aún no hemos llegado a ese nivel de perfección, pero las prótesis robóticas ya tienen mucha funcionalidad. Dejo un ejemplo entre los muchos disponibles en el mercado.

En ese vídeo no se muestra la mano insertada en el brazo de una persona amputada. Por supuesto el principal objetivo de esos miembros son esas personas, pero hay mas opciones. Esos miembros biónicos están desarrollados para que sean controlados mediante señales de los miembros amputados, Pero, en realidad, nada impediría que se pudiesen controlar con otros dispositivos, cómo por ejemplo el que está desarrollando Samsung para “control directo con la mente”. Imagino que en los primeros modelos la sensibilidad y “el ancho de banda” será pequeño, pero según vaya mejorando ese ancho de banda aumentará y habrá un control mas preciso. Y entonces uno podría plantearse si es o no útil manejar algún brazo robótico extra, aparte de los dos brazos naturales. Habrá gente que esto le parezca innecesario, pero a otra no, por ejemplo, el doctor octopus sin duda es favorable a esa idea ;) .

Tal vez los brazos extra no sean del gusto de todo el mundo. Pero ahora que están a punto de estrenar la tercera entrega de Iron Mán seguro que más de uno se plantear si podría disponer de una armadura cómo la de Tony Stark. Bien, pués, hasta cierto punto, sí. Hay muchos proyectos de exoesqueletos. Estos detectan el movimiento del cuerpo humano y lo acompañan, aumentando su potencia. Por cierto, este exoesqueleto lo produce una empresa llamada cyberdyne ¿a alguien le suena ese nombre? ;) :

El uso de estos exoesqueletos es muy amplio. Los militares lo quieren para que los soldados puedan cargar mas material lo cuál proporciona grandes ventajas estratégicas. Esto no es nuevo en la ciencia ficción. Posiblemente una de las novelas mas famosas en las que se plantean esos trajes para uso militar es “tropas del espacio” de Robert A. Heinlein. Pues bien, cómo el resto de casos que he mencionado, ya no es ciencia ficción, esa tecnología la tenemos ya aquí.

Otros candidatos obvios al uso de esos trajes son los parapléjicos que gracias a ellos pueden recuperar la capacidad de caminar. Por ahora estos aparatos tiene autonomía limitada, pero si hay algo en lo que se está invirtiendo mucho dinero y recursos es el tema de las baterías así que, poco a poco (o rápido si hay alguna sorpresa) mejorará mucho esa faceta.

La cuestión de si un exoesqueleto es útil al público en general es abierta en mi opinión. Posiblemente sean demasiado caros cómo para adquirirlos por mero capricho, pero tal vez se vayan generalizando a diversos usos profesionales. Tenemos, por ejemplo, su uso en el manejo de cargas pesadas, cómo el que usa Ripley para pelear con el Alien. Cómo el resto de tecnologías de ciencia ficción que he mencionado ya es una realidad. En el vídeo que dejo se muestra uno militar, pero su aspecto recuerda mucho al de Ripley. Hay otros modelos no militares, pero esto es lo primero que he encontrado.

Podría seguir con el tema de la robótica yendo por el tema de los terminators, que también están aquí (incluso hay ya terminators de aspecto humanoide, no sólo los aviones no trippulados y dirigidoos por inteligencia artificial tipo predator, sus versiones en helicóptero, lancha, etc, y los diveros trastos que caminan vuelan y nadan imitando a animales diversos). Pero no lo voy a hacer. Voy con otra tecnología.

En la versión de cine de “tropas del espacio” se veía de modo muy gráfico una tecnología muy comúnmente tratada en la CF, los órganos artificiales. En concreto era cuando al protagonista, Tico, le curan una horrible herida en la pierna mediante una máquina que va creando a toda velocidad tejido mientras el soldado está en un tanque regenerador.

Una vez más no estamos aún ahí. Pero ya hay tecnologías para imprimir órganos artificiales (véase el enlace para los detalles, que no es exactamente cómo suena, aunque sí muy similar). Aparte de órganos complejos se pueden imprimir huesos. Esto es mas barato, y preciso, que las prótesis metálicas. No sé exactamente los detalles pero adivino que puede ser muy útil en muy diversos casos. Hay mencionan cosas cómo reconstruir mandíbulas o espinas dorsales. Pero imagino que tal vez también pueda ser útil para reponer huesos grandes (tibias, fémur, peroné)que hayan quedado tan dañados en accidentes de tráfico que hagan inviable una recuperación completa de la movilidad.

Toda esta miríada de tecnologías eran ciencia ficción hace pocas décadas ya es realidad. En la mayoría de las novelas se analizaban las consecuencias de uno o unos pocos de estos factores. Ahora tenemos casi todos ellos a la vez, y es difícil adivinar cómo van a interaccionar entre sí. Desde luego cuando entre de manera masiva, añadiéndose a lo que ya hay (móviles y tablets, que son una auténtica maravilla, a la cabeza) van a hacer que el mundo de aquí a 5 años pueda llegar a no parecerse demasiado al mundo actual. Y todo ello en medio de una crisis de modelo social (más que económica) por el empeño de convertir el estado de bienestar en una pesadilla cyberpunk. Y en cuanto esas tecnologías empiecen a aterrizar tal vez ya habrá otras de similar impacto potencial, e incluso tal vez estarán al acecho otras aún mas sorprendentes. Por ejemplo, la nanotecnología tarde o temprano (espero que temprano) posibilitará la construcción del ascensor espacial (y muchas, muchas otras cosas, claro). Hay investigadores que afirman que tienen bastante avanzado un reactor de fusión de pequeño tamaño para añadir a una nave espacial que pueda llevar la nave a Marte en entre uno y tres meses: http://www.abc.es/ciencia/20130408/abci-planean-enviar-humanos-marte-201304081041.html.

En la novela (muy sobrevalorada en mi opinión) “la estrella de Pandora” de Peter F. Hamilton uno de los protagonistas llega a Marte en una nave, tras un largo viaje, sólo para encontrarse con que a lo largo de su viaje se ha desarrollado la tecnología de agujeros de gusano y que le están esperando allí unos universitarios que llegaron antes y sin ningún sacrificio. La creación de agujeros de gusano -que según la relatividad general son posibles- tiene un serio talón de aquiles, que requieren materia exótica. Hasta hace nada eso parecía un problema insalvable. Pero ¿podría haber alguna novedad inesperada al respecto? Bueno, se puede conseguir violar la condición de energía positiva débil (la clave de la materia exótica)
mediante el efecto de Casimir. Y, curiosamente, puede que otro trabajo del casimiro ese pueda ser la clave para que esa sorpresa se haga realidad ;) .

Quantum chaos: las aproximaciones habituales

marzo 21, 2013

Cuando en los 80 hubo la revolución del caos (Ian Malcom’s dixit) inmediatamente se empezó a tratar de ver si podría haber algo similar en mecánica cuántica. El caso es que si bien ha habido progreso en ese área la verdad es que no se ha obtenido nada que pueda considerarse revolucionario.

Los motivos son varios. Uno de ellos es que la gente de quantum chaos no suelen ser físicos teóricos habituados a la teoría cuántica de campos sino gente de física mas aplicada mas habituados a trabajar con la ecuación de Schröedinger que es lineal y no encaja con el paradigma del caos clásico dónde el comportamiento caótico surge de ecuaciones no lineales. En realidad en algunas áreas se trabaja con alguna variante de la ecuación de Schroedinger no lineal. Cómo podéis leer en la wiki la ecuación no lineal surge en física clásica y cuando se pasa a la versión cuántica de esa ecuación se obtiene una ecuación lineal.

No voy a extenderme mucho en lo que se ha hecho en esta línea de trabajo porque no lo conozco muy a fondo. remito a la gente a la entrada de la wiki para el punto de vista standard.

El caso es que en principio la cuántica y el caos son antagónicos. El caos esta ligado a ecuaciones no lineales y la cuántica es la teoría de operadores lineales sobre espacios de Hilbert (equipados si a alguien le importa el rigor), o espacios de Fock (productos tensoriales infinitos de espacios de Hilbert por así decirlo) si vamos a teorías de campos o quizás a construcciones con *-álgebras si se trabaja en las versiones mas rigurosas de cuantización en espacios curvos. En cualquiera de estos casos uno esta trabajando con teorías lineales. Si uno tiene una teoría clásica y la cuantiza la no linealidad de la teoría clásica se pierde y obtiene una teoría lineal. En ese sentido alguna gente se dedica a estudiar cómo la cuantización “suaviza” el caos clásico.

El caso es que uno podría plantearse si en esa linealización no pudiera perderse algo. Cojamos un ejemplo típico, y muy importante, una teoría gauge no abeliana. En esa las ecuaciones clásicas de los campos gauge son no lineales. En la interpretación habitual se considera que un campo nos da la probabilidad de crear una partícula. Para que esa interpretación de partículas tenga sentido lo que se hace es cuantizar la teoría entorno a un vacío concreto de la teoría cuántica. Si el vacío cambia (cómo pasa con el fenómeno de ruptura espontanea de simetría asociada al Higgs) la interpretación de partículas clásica, y de hecho la esencia de la física, cambia. Lo que yo no tengo claro es que la mecánica cuántica, al menos tal cómo la conocemos, sea capaz de describir esa transición entre vacíos. Bueno, sí, tenemos los instantones que pueden ser tratados por integrales de caminos pero no sé hasta que punto no omiten algún aspecto importante.

Bien, entonces, si la cuántica, al menos tal cuál la entendemos, es una teoría de operadores lineales, todo el “fenómeno caos” de evolución sensible a las condiciones iniciales debe ser olvidado ¿o no?. El punto clave es plantearse que es caos. Si uno se plantea esto cómo partir de un hamiltoniano de un sistema y considerar que la no linealidad de este hamiltoniano se traduzca en que una pequeña diferencia entre las condiciones iniciales en el tiempo de las condiciones iniciales evolucione a una gran diferencia de las condiciones finales al cabo de un lapso relativamente corto de tiempo(efecto mariposa) no creo que pueda obtenerse realmente algo así cómo un quantum chaos que sea de una importancia tan esencial cómo el caos clásico.

Pero si nos permitimos una ligera extensión del concepto de caos a algo cómo que una pequeña ignorancia experimental en el sistema pueda conducir a que la evolución temporal de ese sistema amplifique de manera no lineal esas pequeñas ignorancias cuando el sistema evoluciona entonces la historia es muy diferente. Os daría referencias sobre eso pero hasta dónde sé no hay ninguna, y espero que no las haya, al menos hasta que no escriba yo algo al respecto ;) .

¿Realmente los alemanes son buenos en ciencia?

marzo 8, 2013

La respuesta a la pregunta debería ser obviamente afirmativa. ¿Por que plantearla entonces? El motivo es muy sencillo, ese engendro conocido cómo escuela austriaca de economía. Resulta muy extraño que en un país más o menos civilizado se opte por darle carta blanca a las majaderías acientíficas de esa gente. En la época en que se empezó a pergreñar esa abominación intelectual los alemanes del momento lo trataron correspondía, con desprecio. ¿Que ha pasado para que ahora se les de crédito. Por supuesto parte de la respuesta proviene de otros factores (por ejemplo, el diseño antidemocrático de la unión (de empresaurios) europea. Pero aquí voy a ocuparme del aspecto que reza el título, la parte científica.

Negar que Alemania ha sido grande en ciencia es imposible. En matemáticas tenemos así, a bote pronto, a Gaus y a Riemman. Pero desde luego hay más, Jacobi, hilbert, Klein, Weierrstras, Bessel, etc. Podéis ver un listado general en la wiki, uno de los del siglo XIX y otro de los del siglo XX.

En física teórica tenemos a Einstein, a Planck (por esa época hablar de física teórica es un poco aventurado, pero bueno), a Sommerfeld, a Heisenberg, a Schröedinger (mec, no, ese es austricao, pero seguro que por el nombre mucha gente piensa que es alemán. Su formación la hizo en Austria principalmente pero su época mas productiva la realizó en universidades Alemanas) y, en general, la universidad de Gottinga que fué una de las cunas (sino la cuna de la física cuántica.

Y bueno, si seguimos por la rama de la relatividad general tenemos a Karl Schwasrschil, el autor de la famosa solución de las ecuaciones de Einstein estáticas con simetría radial que llevaría al descubrimiento de los agujeros negros y luego otro nombre muy sonoro, la soución de Reissener-Nordstrom para agujeros negros cargados (Reissner era alemán, nordstron sin embargo era finlandés) y luego la solución de Kerr-Newman para un agujero negro rotatorios que, ah, no, estos no son alemanes, Kerr es neozelandés y Newman es americano, pero vaya, el nombre podría dar el pego así que por eso lo aclaro.

Bien, queda claro, en el siglo XIX y en la primera mitadl del XX en matemática había muchos y muy buenos matemáticos, y a finales del XIX y principios (hasta la 2ª guerra mundial) del siglo XX tenían a muchos de los mejores físicos teóricos. Hasta aquí más o menos lo más conocido y lo que da fama a los alemanes como científicos. Pero vamos a intentar ir un paso mas allá en el análisis.

Leyendo un libro de divulgación sobre la conjetura de poincaré me enteré de cosas muy curiosas sobre la historia de las matemáticas. Cosas cómo que por ejemplo los USA habían alcanzado su máximo número de matemáticos en los 70 y desde entonces han ido cayendo y que en Alemania sólo empezó a haber matemática (y sospecho que física) en serio a raíz de que Inglaterra empezase a ser la primera potencia del mundo, empujada por la revolución que supuso la ciencia de Newton y, secundariamente, la revolución industrial. Puede verse la poca relevancia matemática de los alemanes en el siglo XVII consultando el listado de matemáticos alemanes de ese siglo. Realmente conocidos están Gauss, que es más del XIX que del XVIII, Leibniz y si acaso Goldbach.

A esa nueva ciencia se habían ido sumando otros países, cómo por ejemplo Francia, o Suiza y Holanda, y les iba muy bien. Los alemanes, viendo que se quedaban atrás, terminaron por optar por invertir en ese área y lograron que un ser tan singular cómo Gauss se quedase en tierras bávaras en vez de irse a una universidad Suiza (para lo cuál tuvieron que pagarle muy bien). Bien, ese idílio de Alemania con la matemática y la física se prolongó hasta que llegaron los nazis al poder y apostaron por una ciencia “práctica” (que enseguida se vió que no valía un pimiento) y la mayoría de los grandes físicos, con Einstein a la cabeza, huyeron del país.

¿Que fué de la física en Alemania después de los Nazis? Bien, Alemania durante la postguerra tuvo serios problemas, aunque económicamente se recupero más o menos bien- y relativamente rápido- a raíz del plan Marshall. Ahora lo que es en física teórica nunca se recuperaron. Siendo físico teórico debo decir que no recuerdo ahora mismo (y he estado pensando en ello un rato) ningún nombre de teórico alemán que haya hecho alguna contribución importante desde la 2ª guerra mundial. Ahora mismo, en teoría de cuerdas, me suenan dos nombres, Lüest y Theisen, autores de uno de los primeros libros sobre teoría de cuerdas que está escrito e a raíz de un curso en una universidad alemana lo cuál junto a la sonoridad de los nombres, me lleva a pensar que los autores son alemanes. Luest, aparte de ese libro, me resulta conocido por sus trabajos en la construcción de modelos fenomenológicos a partir de D-branas, en particular sus artículos “the LHC string hunter Iy II” que, cómo el resto de modelos fenomenológicos basados en cuerdas han sido vapuleados por el LHC (que, por ahora, no ha cazado nada, usando la terminología de los artículos).

Si dejamos la teoría de cuerdas y aceptamos a día de hoy a la LQG cómo una alternativa válida en gravitación cuántica (algo que mucha gente discutiría enérgicamente) podriamos señalar dos nombres, Martin Bojowald, que trabaja en loop quantum cosmology y en asymptotic safe gravity.

Y luego tenemos a Sabine Hossenfander que es famosa por su blog, back-reaction, y tal vez pueda decirse que es relativamente buena en fenomenología de la gravitación cuántica. Y ya, no soy capaz de recordar ni un sólo nombre de físico teórico alemán relevante.

Puede parecer algo el indicar estos nombres. Pero, sinceramente, si ese es el nivel de la física teórica alemana no está mejor que la Española. Podria señalar nombres cómo Luis Ibañez, importante fenomenólogo de cuerdas, por ejemplo, los modelos de compactificación LARGE, de los que es creador junto al también muy reconocido J.C. Quevedo , Luis Alvarez-Gaumé, autor español cuyo nombre aparece en cualquier libro de teoría de cuerdas por un trabajo que hizo en colaboración con Edward Witten sobre anomalías en teoría de cuerdas, Enrique Alvarez que posiblemente sean mejores que los físicos Alemanes mencionados (aunque ya sabemos que las comparaciones son odiosas).

Pero vale, podría ser que me fallase la memoria o algo a´si. Vamos pues a consultar la lista de premios nobel de física. Observamos que entre el 1932, año en que le dan el nobel a Heisenberg y el 54, en el que se lo dan a un tal Bothe (de cuya existencia me acabo de enterar y que el tema por el que le dieron el nobel me resulta absolutamente desconocido) no hay ningún Alemán. Osea, que en 32 años no hay ni un sólo premio noble Alemán y el que rompe la sequía es de una talla infinitamente menor que Heisenberg. Pero- vale, estaban recuperándose de la guerra y de la fuga de sus grandes teóricos. Sigamos, el siguiente premiado llega 7 años después, en el 1961, Mößbauer. No es un precisamente teórico pero al menos su trabajo es relativamente conocido. Podríamos pensar que ahora, una vez recuperados económicamente, y tras juntar dos nobel en 7 años, volverían a la gran época. Pero no, el siguiente nobel se hace esperar hasta el 1985 ¡24 años sin un nobel cuando ya sí son de nuevo una potencia económica!. El galardonado es un tal Klaus von Klitzing por el efecto Hall cuántico. dado que el tema por el que fué galardonado lleva el nombre de otro podemos sospechar que no fué el que mas aporto al tema. Si entramos a leer los detalles vemos que en realidad lo que hizo fué estudiar cómo varía a saltos la corriente en conductores en determinadas circunstancias. Soy teórico y no físico del estado sólido, pero aún así puede afirmar que su trabajo no tiene ni de lejos la relevancia de temas de ese área 8estado sólido/materia condensada) cómo los trabajos de Landau, el descubrimiento de los transistores de estado sólido por parte de Cooper o la teoría BCS (Bardeen-Cooper, Schaefer) de la wsuperconductividad.

Ya, mas recientes, tenemos en el 87 a Johannes Georg Bednorz galardonado por el descubrimiento de la superconductividad en materiales cerámicos, que hizo abrigar la esperanza de tener superconductores a temperatura ambiente. Dado que incluso hoy no tenemos una teoría medianamente fundamental (y ni siquiera si tenemos una teoría digna del nombre, pero vamos, no es mi campo de estudio así que ni afirmo ni niego) está claro que es un descubrimiento experimental en un tema experimental. Es decir. 0% teórico. Para no ir caso por caso decir que desde entonces han tenido algún nobel más, pero todos en temas de física del estado sólido u otros temas aplicados, cómo optronica. Es decir que desde la huida de físicos teóricos durante el nazismo nunca han vuelto a tener un nobel en ese área. Y hasta hace relativamente poco los dos nobel intermedios que tuvieron no se puede decir que hayan sido de los mas prestigiosos precisamente.

Bien, vamos con el otro gran apartado de la ciencia de élite, las matemáticas. En concreto vamos al premio mas prestigioso de las matemáticas, la medalla fields y sus ganadores. Aqí el trabajo es muy sencillo. Desde que se instaura el premio, en el año 1936, ¡sólo un matemático alemán ha obtenido el galardón!. En concreto se trata de Gerd Faltings que obtuvo el premio por demostrar una conjetura de geometría algebraica, la conjetura de Mordell (o conjetura de Morder-weyl que es lo que le oigo decir siempre a un amigo mío, doctorado en teoría de números y geometría algebraica). Desde luego el bagaje en la medalla fields de los matemáticos alemanes es absolutamente lamentable. Mas patético resulta si lo comparamos con otro gran país Europeo cómo Francia que tiene 11 medallas fields.

Bien, esta claro que las medallas fields dejan a la altura del betún a Alemania. Vamos a ver si por algún otro lado levantan algo la cabeza. Veamos de nuevo la lista de matemáticos alemanes del siglo XX.

Tenemos algunos nombres ilustres algunos y cuanto menos destacados otros. EStán Cantor, Hilbert, Klein, Dedekind, Frobenius, Hartogs, Hólder, Klein, Kutta, Maric, Minkowsky, Neumann y Noether. Si miramos sus biografias observamos una pauta común. Son gente nacia en el XIX o principios del XX y que mueren antes de mediados de siglo. Los otros nombres yo no los conozco de nada. Tal vez sea porque sus contribuciones son muy recientes o porque no son terriblemente relevantes (o ambas cosas si hemos de juzgar por su no aparición en las medallas field).

Resumiendo. Alemania llegó algo mas tarde que otros países a la ciencia moderna (la que vino a raíz de Newton) pero cuando lo hizo obtuvo grandes resultados con matemáticos y físicos de primer nivel. Su época dorada puede situarse entre mediados del XIX y la segunda guerra mundial Después de eso en física teórica nunca volvieron a tener un nivel digno de comentar y sólo en ramas mas aplicadas de la física han remontado un poco cabeza, en especial a partir de los 80. En matemáticas la época de gloria también va en la misma época, pero con una diferencia, tras la diáspora de científicos de primera fila después del nazismo nunca han vuelto a ser nada en matemáticas.

Y ahí estamos, ahora en Alemania no hay físicos teóricos ni matemáticos especialmente dignos de mención. Con ese panorama tan patético se entiende un poco mejor que hayan podido caer en manos de los trastornados de la economía austriaca. O a lo mejor es al revés y su aquiescencia con la descerebrada economía austriaca y su anarcopapialismo les ha llevado a la indigencia científica en la que están actualmente.

Si nos vamos a la tecnología…nos alargaríamos mucho, y tampoco es el terreno que mas me entusiasme. Pero vamos, no todo es malo en la Alemania reciente, hicieron algo realmente bueno, traer el wing tsun de china y darle una excelente difusión en Europa. Y bueno, tienen buena música electrónica (sería interesante ver que paso con su música clásica, que así de primeras me parece que también tuvo un gran bajón después del nazismo del que no se llegó a recuperar).

¿Algo más digno de destacar de la Alemania actual? yo diría que no :P .

Y para cerrar un dato muy revelador. El instituto científico Alemán mas prestigioso es el Max Planck. Pues bien, el director del mismo es un Español, Ignacio Cirac, que está en las quinielas de los nobel año sí y año también.

Los cocodrilos y la segunda ley de Snell. Ataque desde el ángulo límite.

marzo 7, 2013

Últimamente estoy -por variados motivos- perezoso a la hora de escribir en este blog y en los otros. Hoy salgo de esa atonía para escribir sobre un tema ligero de física elemental. Podría haber hablado de la conferencia Moriond dónde se están publicando análisis de los datos del año pasado del LHC, pero cómo hasta ahora casi todo sigue en la tediosa tónica de no dar indicios decentes de física de mas allá del modelo standard dejo a otros ocuparse de ello. Podría hablar de los últimos resultados sobre aspectos cuánticos de agujeros negros obtenidos a partir de teorías de cuerdas (o teorías conformes) pero preferiría hacerlo cuando haya algún resultado experimental, aunque sea en análogos de materia condensada. En definitiva, a la espera de la publicación en breve de nuevos datos sobre observación experimental de la materia oscura la física teórica parece un poco estancada a la espera de alguna revolución inesperada, que desde luego es algo que no puede descartarse ;) .

Bien, pues vamos ya con los cocodrilos, que cómo dice una canción (“tu por mí”) de Cristina Rosenvingue hay que tener cuidado con ellos, que vienen despacio, y nunca los ves. Admito que soy muy aficionado a los documentales “de bichos” (y sin son dinosaurios mejor) y también a las películas sobre el mismo tema. El mes pasado en cuatro pusieron un ciclo de pelis de esa temática con el termino comercial “animalator” y aunque ya había visto varias las volví a ver, y las que no había visto, pues con mas motivo también las ví.

El caso es que dando clases de biofísica he tenido que repasar varias veces la ley de snell (que en su momento ví en óptica de tercero) y haciendo con un alumno uno de los problemas se me ocurrió una relación entre ella y las pelis (y documentales) de cocodrilos.

Cómo se observa en el vídeo del inicio una de las tácticas de los cocodrilos es acercarse sigilosamente a las orillas de los ríos dónde esta abrevando algún ñú o gacela. El caso es que los animales de la orillas están muy atentos al agua pero aún así no ven llegar al cocodrilo. En los documentales los biólogos echan la culpa a que el agua está muy sucia. Sin duda ese se un factor, pero ¿es el único factor? Yo diría que no. Me parece que tal vez la segunda ley de Snell y el ángulo límite jueguen un papel mas importante.

Recordemos cuál es esta ley:

n_i.sin (\alpha_i) = n_r. sin (\alpha_r)

En la fórmula i hace referencia al ángulo incidente y r al refractado. Los ángulos se miden respecto a la normal que separa los medios. Un dibujo dónde se puede apreciar la geometría implicada es este.

Para los propósitos de tema tratado debemos imaginar que el rayo incidente proviene de el agua y el refractado sale al aire.

Lo interesante es que cuando nr>ni va a existir un ángulo incidente, denominado ángulo límite, para el cuál el refractado sale con ángulo de 90%. Es decir “sale” paralelo a la superficie y por tanto no llega a salir realmente.

Si ponemos en la fórmula que el ángulo refractado es 90 y despejamos obtenemos la expresión del ángulo límite:

\alpha_{lim}= arcsin( n_r/m_i) .

Para mas detalles ver la entrada de la wiki

Con datos más o menos típicos (ni=1.33, nr=1) tenemos que el ángulo límte es del orden de uos 50º. Entonces lo que tenemos es que si el cocodrilo va bajo la superficie del agua, pero muy cerca de la superficie el animal de la orilla. con la cabeza pegada al agua- tendría que el cocodrilo se acerca a él con un ángulo mucho menor de esos 50º incluso cuando está casi al lado. Por tanto debido a la ley de Snell el cocodrilo es totalmente invisible para el animal, incluso si el agua estuviese totalmente limpia. Le basta pues al reptil con ser sigiloso y no agitar el agua para poder atrapar a su presa por muy buena vista que tenga ésta y muy atenta que esté a la superficie del agua. La opción que tendría el animal es levantar de vez en cuando la cabeza a una altura mayor pués así el cocodrilo estaría respecto a la cabeza a un ángulo mucho mayor y sería visible.

Ya, cómo idea peregrina, me pregunto si podría haber algún motivo para que el espacio de algún modo tuviese algún ángulo de refracción por algún exótico motivo, y que eso se diese de manera istótropa. En ese caso desde cualquier ángulo a la tierra sólo llegaría una cierta fracción de luz, ya que el resto quedaría fuera del ángulo límite, y habría potencialmente mucha materia que sería invisible incluso aunque emitiese luz. Eso daría un escenario para explicar lo de la dichosa materia oscura, pero vamos, me resulta complejo imaginar cómo reunir las condiciones que pido para esa refracción isótropa en el vacío.

Cómo conclusión, haced caso a Cristina y mucho cuidado con los cocodrilos ;) .

P.S. Ahora entiendo porque Cristina se dió cuenta de que venían los cocodrilos. Si se sumerge la cabeza en el agua no hay fenómeno de refracción y se ve llegar al cocodrilo desde cualquier ángulo ;.

Por cierto, nunca imaginé a Cristina usando electrónica en la base de sus temas. Imagino que ahora que ha aprendido a tocar el piano también ha empezado a usar sintetizadores y de ahí el cambio. Sea como sea muy mono el tema, que desconocía hasta hoy mismo.

Introducción a la LQG canónica II

diciembre 24, 2012

Continúo la anterior entrada y concluyo esta introducción a la “vieja” LQG. Insisto, a estas alturas la LQG tiene poco prestigio y hay buenos motivos para no tomársela muy en serio, pero. cóm dije antes, desde un punto de vista académico, con vistas a otros asuntos, tal vez puedan ser útiles algunas técnicas matemáticas de la LQG.

——————————————————————————————

Habíamos visto por encima el formalismo ADM. Vamos a mantener la misma foliación del espacio-tiempo que para evitar confusiones de notación aquí denoto por Et=MxR.

Voy a poner algunas expresiones explicitas de las fórmulas propias de ese formalismo porque en la LQG aparecerán de distinta manera. Si alguien se marea con las fórmulas tranquilidad, que intentaré discutir el significado intuitivo de las mismas para que se pueda seguir la argumentación sin depender (demasiado) de ellas:

V_a = -2 q_{ac} D_b P^{bc}

Aquí D es una diferencial covariante compatible con qab, la métrica en el slice (corte) M. Y los P´s son los momentos canónicos definidos por:

P^{ab}=\frac{\kappa\delta S}{\delta q_{ab}}=\sqrt{\det(q)}  [q^{ac} q^{bd}-q^{ab} q^{cd}] K_{cd}

Esta es la ligadura de difeomorfismos espaciales. Representa la invariancia de teoría bajo cambios de coordenadas.

C=\frac{1}{\sqrt{\det(q)}}  [q_{ac} q_{bd}-\frac{1}{2} q_{ab} q_{cd}] P^{ab} P^{cd}-  \sqrt{\det(q)} R

En términos de estas dos ligaduras el hamiltoniano queda:

H=\int d^3x [\lambda P+\lambda^a P_a+U^a V_a+N C]

Aquí los varios P´s que aparecen son los momentos conjugados a qab, N y Ua. Las λ son multiplicadores de Legendre que sirven para implementar las ligaduras (es algo similar a los multiplicadores de legendre que aparecen en problemas de minimización de funciones sometidas a una condición). Se puede ver que las Ua y N juegan también un papel cómo ligaduras. De hecho debo decir que esta fórmula esta tomada de uno de los reviews clásicos de la LQG, el de Thieman, y en este punto resulta algo confuso, la relación con formas mas habituales de este hamiltoninao puede verse del siguiente modo: lo que ocurre es que esas P’s están relacionadas con las ligaduras C y V, de hecho son ellas mismas cambiadas de signo, con lo cual se eliminan del analisis efectuado sobre la base de las ligaduras de primera clase, y al final solo te queda el Hamiltoniano dependiendo de C y V, siendo la función lapso y el vector desplazamiento multiplicadores de Lagrange, y así se obtiene el Hamiltoniano empleado.

A modo de resumen rescato lo más importante de lo introducido en este formalismo. Tenemos las variables que juegan el papel de “coordenadas”, las qab y sus momentos conjugados Pab. Cómo no se pueden “despejar” unos en función de otros en el formalismo hamiltoniano aparecen dos ligaduras. Una, V, que representa la invariancia bajo difeomorfismos y otra, C, que representa la dinámica del sistema y por eso se conce cómo ligadura hamiltoniana. Esta ligadura es la que, cómo dije en la entrega anterior, se conoce cómo ecuación de Wheeler-de Witt, bueno, para ser mas exactos la “versión cuántica” de esta ligadura, ya explicaré al final este matiz.

Dejo una imagen que pueda aclarar un poco esto de las foliaciones, que es importante visualizarlo para entender la parte clásica de estas teorías:

Cómo avancé la LQG era una especie de “cambio de variables” para el formalismo hamiltoniano. Para entrar en LQG necesitamos introducir el formalismo del vielbein:

Este formalismo (también conocido cómo el de la tétrada) fue introducido por Elie Cartán. Entre otras cosas permite introducir fermiones en la relatividad general, algo muy útil, pero aquí no haré hincapié en ese aspecto.

La tétrada eui, o vielbein, es (así le gusta decir a la gente) una especie de raíz cuadrada la métrica (en realidad es un sistema de referencia orto-normal, bueno algo así), las ecs. que cumple son:

g^{\nu\mu}=e_\nu^ie_\mu^j\eta^{ij}

\eta_{ij}=e_\nu^ie_\mu^jg^{\nu\mu}
η= métrica de Lorentz.

Esta construcción vale también para la métrica tridimensional q definida en los cortes M de la foliación del espacio-tiempo en MxR.

Asociada al vielbein de esta la conexión de spin:

\Gamma_{\nu}^{i_j} que cumple

\partial_a e_b^j-\Gamma^c_{ab} e_c^j+\epsilon_{jkl} \Gamma_a^k e_b^l=0
Dónde \Gamma_{ab}^c son los símbolos de Christoffel asociados a la métrica qab

Aprovecho para indicar que en el formalismo de la tétrada hay un aspecto interesante, la posibilidad de considerar la conexión, y no la métrica cómo el campo básico que describe la gravitación, formalmente esto quiere decir que aplicamos el principio de mínima acción a este campo y no a la métrica. Es este detalle el que permite que en LQG haya la posibilidad de tener gij=0, es decir, que no haya espacio, mientras que en el formalismo normal la ausencia de gravedad indica que tenemos η cómo vacío.

Pero sigamos, vamos a introducir una nueva cantidad A. Esta cantidad va a jugar un papel similar a un campo gauge SU(2), pero con matizaciones que ya iremos viendo. Uno se puede preguntar, con toda la razón del mundo, ¿que pinta aquí un grupo SU(2)? Voy a intentar aclararlo un poco.

Es importante de que no detallé mucho a cómo se introducía la conexión, bien, esa conexión se puede asociar con un grupo, el grupo de la fibre bundle. El caso es que aquí pasamos a trabajar en un slice tipo ADM.
Por tanto tenemos una 3-variedad, y con esa 3 variedad construimos un fibrado. ¿Que grupo tiene ese fibrado?, pués el grupo de invarianzas locales. Puesto que estamos en un slice riemaniano ese grupo es el SO(3), es decir, localmente el espacio es euclideo y tiene invarianza bajo rotaciones euclideas. SO(3) tiene de grupo recubirdor SU(2), así que ya hemos llegado a puerto. Si no trabajáramos con el slice tendríamos una 4-variedad y el grupo sería SO(3,1) (Correspondencia al la invarianza lorentz local).

Para quien no sepa nada de fibrados intentaré explicarlo de otro modo. A muchos os sonará el siguiente “experimento mental” de Einstein. Una persona en caída libre dentro de un campo gravitatorio no experimenta ninguna fuerza. Sin embargo si esa misma persona se ve parada en su caída por algún tipo de superficie notará que el campo gravitatorio le empuja contra la superficie. Es decir en ese punto tenemos dos posibilidades, estar en caída libre y no sentir ningún campo gravitatorio y estar sometido a una fuerza gravitatoria. En términos matemáticos eso significa que localmente siempre podemos encontrar un sistema de referencia, el de caída libre, en el que el campo gravitatorio no existe. Pero un sistema sin campo gravitatorio se describe por la relatividad espacial, que tiene el grupo SO(3,1) como invariante). Ahora bien, esto puede hacerse sólo localmente. Para dos puntos cercanos el sistema “libre” de uno no lo será para otro. Lo que se puede, sin embargo, es asociar a cada punto del espacio tiempo un sistema de referencia en caída libre (sistema de referencia inercial). Bien, eso es una fibración. Y en este caso cómo cada punto de la fibra (cada sistema de referencia inercial) tiene un grupo de invarianza tenemos un fibrado principal en el que el grupo actua sobre cada fibra. Sí, vale, imagino que este intento de explicación “fácil” no es muy comprensible, pero es que los fibrados no son la cosa más sencilla del mundo.

Bien, tras esta explicación pongo la expresión de A.

A_a^j=\Gamma_a^j+\beta K_{ab} e^b_j,\;\;E^a_j=\sqrt{\det(q)} e^a_j/\beta

Aquí Kij es la curvatura extrínseca, β es un parámetro libre de la teoría conocido cómo parámetro de barbero-Inmirizzi y ya iremos viendo que juega un papel importantísimo en la teoría. Puede tomar cualquier valor complejo distinto de 0. En la formulación inicial de Astekhar este parámetro tomaba el valor i (raíz de -1). Con este valor la forma de la ligadura hamiltoniana (que ya veremos) tomaba una forma muy sencilla. Para este valor el SU(2) asociado es un SU(2) complejo y por tanto no compacto. Sin embargo si se toman valores reales el SU(2) es real y por tanto compacto. Esto facilita muchas cosas, pero no en particular la forma de la ligadura hamiltoniana que en general tendrá una forma complicada, pero ya digo que iremos viendo cosas varias de este parámetro a lo largo del tema.

Ya tenemos unas nueva variables, cómo aduje que era nuestra intención buscar, Las A´s (“campo gauge”) y las E´s (tétrada densitizada, i.e. dividida por la raíz de la métrica). Se puede expresar el lagrangiano de la RG en estas variables y se pueden calcular los momentos canónicos asociados. Cómo era de esperar aparecen ligaduras. Y supongo que no es de extrañar que las ligaduras en lsa nuevas variables guarden similitud en su significado a las ligaduras en las otras variables, así pués las denotamos igual y tenemos:

V_a = F_{ab}^j E^b_j

Aquí F es el análogo al tensor que en los campos gauge juega el papel de tensor campo electromagnético y tiene la expresión:

F=2(dA+[A,A])

La “ligadura hamiltoninana” tiene la forma:
C = \frac{F_{ab}^j \epsilon_{jkl} E^a_j E^b_l}{\sqrt{|\det(E)|}}\;\;+\;\;  \mbox{ Mas }

Dónde “Más” indica términos adicionales de similar factura.

En este formalismo sin embargo aparece una nueva ligadura. Cómo he dicho La A juega el papel de un campo gauge y se sabe que los campos gauge tiene una ligadura, la ley de Gauss:

G_j:=\partial_a E^a_j+\epsilon_{jkl} A_a^k E^b_l

Lo interesante es que ahora todas las ligaduras son polinómicas en los campos, en contraste con le caso ADM en que no lo son. Por tanto son más sencillas de resolver. Más sencillas, evidentemente, no quiere decir sencillas, de hecho son dificilísimas de resolver, pero al menos se puede hacer algo con ellas. Todo el programa de la LQG canónica consiste, en realidad, en “pelearse” con estos malos bichos. Ya iremos viendo cómo se hace, pero voy a dar algunos matices más necesarios para entender por dónde iremos.

Lo primero es señalar que puede que no os hayáis dado cuenta pero hasta ahora todo lo que se ha hecho en este post es RG clásica. Todos los campos son clásicos. Las ecuaciones de ligadura pueden tratarse y resolverse de manera puramente clásica mediante, por ejemplo, el programas de Dirac, que consiste en cambiar los corchetes de Poisson por “corchetes de Dirac” que seleccionan los estados físicos y tal y cuál. Pero en LQG se va a intentar hacer relatividad cuántica. Así pués, siguiendo los pasos habituales, se van a sustituir las “posiciones” y los “momentos” por operadores cuánticos. Aquí las “posiciones” eran las qab y los “momentos” Las E´s. Pués bien, siguiendo el procedimiento habitual en QFT se convierten estos campos clásicos en operadores que actúan sobre un espacio de Hilbert. Ahora bien, las constraints se expresan en términos de estos campos ¿cómo se van a expresar en la versión cuántica? Pues del mismo modo que se hace con las teorías gauge normales. Se impone que los estados físicos sean anulados por las ligaduras, es decir:

L|φ>=0

Donde estoy denotando por L una ligadura arbitraria y |φ> representa el estado cuántico. Bien, ahora ya he explicado con más detalle matemático el programa de la LQG y he introducido su parte clásica. No voy a analizar en este post, ya lo dije antes, cómo se resuelven las ligaduras, pero al menos si voy a explicar de dónde viene el nombre de loop quantum gravity.

En las teorías gauge ordinarias, basadas en un campo gauge A, se construyen los “wilson loops” correspondientes. ¿Que es un Wilson loop? la definición matemática es esta:

dónde C es un contorno (curva cerrada) y P es el “path-ordering” que no definiré pués no me voy a meter a fondo con los wilson loops. En teorías gauge ordinarias los wilson loops permiten reconstruir, a nivel clásico, la conexion gauge. A nivel cuántico son operadores que crean excitaciones del campo gauge localizadas en el lopp (el contorno C).

Eso en teorías gauge “ordinarias” ¿que papel juegan en LQG?. Bien, más o menos la idea es que también se convierten en operadores. Recordemos que los operadores en mecánica cuántica tienen autoestados. Pués bien, los autoestados asociados a los Wilson loops generan estados cuánticos que satisfacen dos de las ligaduras anteriores, la ley de gauge y la V de los difeomorfismos. Aparte de eso tienen otros significados intuitivos relacionados con áreas y demás. Pero ya anticipé que no son la mejor base posible para la LQG y han sido abandonados en favor de unos nuevos tipos de estados construidos a partir de grafos (con ciertas relaciones de equivalencia metidas por ahí en medio), las spin networks que también resuelven las dos ecuaciones de ligadura y que nos dan el espacio de hilbert “kinemático”. Pero eso para una próxima entrega.

Por cierto ¿habéis visto alguna pega hasta aquí? No, ¿verdad? Pués veréis, resulta que se aduce que este cambio de variables es de algún modo local (imagino que se refieren a que hay que hacerlo carta por carta -hablando en términos de variedades diferenciables-). El caso es que, anticipo ya, uno de los resultados de la LQG es que el operador área actuando sobre los estados de las redes de spin tienen un espectro discreto, es decir, que el área esta discretizada en LQG. Pué bien, se critica que esa discretizacion se debe exclusivamente a ese carácter local del cambio de variables. Yo particularmente es algo que no veo nada claro (no he visto el desarrollo pormenorizado de esa crítica). Pero la verdad es que es algo absolutamente crucial para todo el programa de la LQG, tanto el canónico cómo el resto. En la LQG canónica porque varios de sus éxitos dependen de ello. Muy asociada a la LQG están la teorías de double special relativity que imponiendo la la relatividad especial de Einstein la existencia de una longitud mínima (algo razonable si existe un área mínima) implican que la velocidad de la luz en el vacío depende de su frecuencia (algo que podría ser constatado experimentalmente dentro de no mucho si es cierto). Las DSR tiene sus propios problemillas, pero sin ese resultado de cuantizacion del área pierden su más fuerte motivación. Pero es más grave aún, algunos puntos de partida para teorías tipo LQG (triangulaciones dinámicas, cálculo cuántico de Regge) parten de un espacio-tiempo discretizado. Aunque hay motivos heurísticos para considerar plausible esa discretización la verdad es que ese resultado de que el operador área tiene espectro discreto es un resultado mucho más potente que da mucho más crédito a todos esos programas. Si ese resultado esta mal, es decir, es un producto de la elección de variables y no refleja ninguna verdad “profunda” la LQG pierde bastante de su sentido.

Introducción a la LQG canónica I

diciembre 24, 2012

Había escrito esto, hace ya años, antes en otros blogs, pero cómo quiera que esos no tiene soporte nativo para Latex resulta que ahora no se ven las fórmulas. Al ponerlo aquí me aseguro de que no tengo que andar retocando nada en el futuro. Por cierto, aunque hubo una época que parecía una buena opción actualmente la LQG está bastante perdida cómo una opción válida para tener una teoría de la gravedad cuántica. Eso sí, puede resultar interesante cómo una herramienta matemática para analizar cierto tipo de cuestiones de un interés posiblemente académico y por eso me he molestado en rescatar esos antiguos posts. Sin mas dilación voy con ello.

——————————————————————————————–

He hablado mucho por aquí de la loop quantum gravity pero hasta ahora no había explicado nada de ella. A raíz de empezar a participar en un foro de CF, http://www.sedice.com, me encontré con que me pidieron un post sobre el tema. He aquí la primera parte del mismo (ligeramente retocada). Aviso, en el post se repiten cosas que ya he explicado, y mas extensamente en este journal, pero creo que no viene mal tenerlo aquí todo junto.

Respecto a los conocimientos previos voy a intentar abarcar varios espectros de lectores. a todos ellos les asumiré una cierta familiaridad con la cuántica y con la relatividad general cuanto menos a un nivel elemental. Pero también tendré en mente a quien tenga una sólida base en física y no conozca este campo concreto. Imagino que cada cuál podrá sacar de párrafos concretos una idea que estará en función de su nivel. Y tras estos prolegómenos vamos con el tema.

Primero explico porque se investiga en estos temas. En la física establecida tenemos dos grupos de teorías, unas que explican todos los datos experimentales relacionados con las interacciones en las que interviene las fuerzas electromagnéticas y las nucleares débiles y nucleares fuertes. Estas son teorías cuánticas. Para ser concretos son teorías cuánticas relativistas. Están son bastante distintas de las teorías cuánticas “ordinarias”. Paso a explicar porque pues es relevante y me sirve par introducir unos conceptos que usaré más adelante.

En cuántica no relativista el problema tipo es resolver la ecuación de Schröedinger independiente del tiempo. Esta es una ecuación en derivadas parciales para una función de onda. Una vez resuelta se tiene la expresión matemática de esa función de onda que nos da información sobre el sistema, los observables (energía, momento angular total, momento angular en el “eje z” y spin normalmente). La solución a esa ecuación no va a ser casi nunca única y va a haber una familia de funciones de onda. Cada una tendrá sus propios valores de los observables.

He usado la expresión “ecuación de Schröedinger”. Esto es correcto, pero conviene resaltar un punto. No hay una sola ec de Schrëdinger. Par cada sistema físico hay una. Por “sistema físico” me refiero a situación física bajo estudio. Si estoy estudiando el átomo de hidrógeno esa ecuación tendrá un término que haga referencia al potencial culombiano. Si estudio el oscilador armónico habrá un término que haga referencia al potencial armónico. Esto es un poco lioso si sólo vemos la ecuación como una ecuación en derivadas parciales. Las cosas se aclaran más introduciendo algunos conceptillos.

En mecánica clásica newtoniana tenemos la famosísima ecuación de movimiento F=m.a. (fuerza igual a masa por aceleración) que es una ecuación diferencial. Ahí al igual que en cuántica, tenemos que F varía de un sistema físico a otro. Aunque útil la ec. de Newton es un poco latosa de usar en muchos casos prácticos. Por eso surgieron otros formalismos matemáticamente equivalentes. Hablaré del lagrangiano y el Hamiltoniano.

La idea del Lagrangiano es escribir una expresión que representa el sistema físico. Y a partir de esa expresión dar un “algoritmo” que permita obtener las ecuaciones de movimiento. Para sistemas clásicos este lagrangiano normalmente va a ser de la forma:

1. L= T-V

Aquí T es la energía cinética y tiene la forma T=1/2mv 2 y V es el potencial en el que se mueve la partícula. Cómo sabréis, espero, la velocidad es la derivada respecto al tiempo de la expresión que nos da la posición de la partícula en un momento dado.

Bien, este lagrangiano es clave en mecánica clásica y en teoría cuántica relativista. En teoría clásica porque se extraen de el las ecuaciones de movimiento (y otra información) En teoría cuántica relativista (a partir de ahora QFT, i.e. quantum field theory) aparte de las ecs. clásicas se obtiene a partir de el Lagrangiano muchas otras cosas (cantidades conservadas, elementos de matriz S, etc) mediante otros procedimientos. Aprovecho ya para indicar una diferencia. Hay un lagrangiano para una partícula y un lagrangiano para un campo. El lagrangiano de una partícula da, como dije, las ecuaciones de Newton. El lagrangiano para un campo da las ecuaciones de campo. Por ejemplo se puede escribir un lagrangiano cuyas ecuaciones clásicas son las ecuaciones de Maxwell del electromagnetismo. Para ser precisos se habla de lagrangiano para una partícula y de densidad lagrangiana para un campo, pero no siempre es necesario precisar tanto. Es importante señalar que el lagrangiano de un campo no siempre tiene la estructura de la ecuación 1. Por ejemplo el del campo electromagnético no tiene esa estructura.

Bien, dije lagrangianos y hamiltonianos. Aún no he hablado del hamiltoniano. En realidad todos habéis visto, en cierto modo, el hamiltoniano. El hamiltoniano s una expresión para la energía. En el caso sencillo de una partícula moviéndose en un potencial el análogo de la expresión 1 será:

2. H= T +V

Voy rápidamente con matizaciones muy importantes. El hamiltoniano es más que la expresión de la energía. Tanto en mecánica newtonina cómo en Lagrangiana se pueden obtener expresiones para la energía que son de la forma 2. Así pues hay más truco. Primero aclaro que es la energía. En formalismo lagrangiano es “una integral primera de las ecuaciones de movimiento”. Intuitivamente puede verse que ciertos sistemas, los no disipativos, conservan una magnitud en su evolución temporal, esta magnitud conservada es la energía.

La clave que diferencia el formalismo hamiltoniano del lagrangiano y que es fundamental en LQG es que en el formalismo hamiltoniano se expresa la energía no en términos de las velocidades sino en términos de los momentos. Explico esto para la partícula libre (i.e. V=0). En este caso L= 1/2 m.v2 . Bien, todos habréis oído hablar del momento lineal, denotado p. Y recordareis que tiene la expresión P=m.v. Así pues podemos despejar v en función de p, i.e. v=v/m. Si sustituimos en el lagrangiano nos da H=(1/2m). P2.

Bien, esto tan sencillo aparentemente oculta muchísimas sutilezas que son claves tanto en teorías gauge como en LQG. Para empezar todos sabéis que es el momento he dicho. ¿Seguro? En realidad no es tan sencillo. El momento tiene una definición técnica, una derivada parcial del Lagrangiano respecto a la velocidad (con una variante adecuada para densidades lagrangianas relativistas y no relativistas).

Técnicamente el Hamiltoniano surge de una “transformada de Legenddre” del lagrangiano, pero a efectos prácticos puede verse cómo una expresión obtenible del hamiltoniano mediante una operación matemática sencilla (que en el caso de la partícula clásica se traduce en el cambio de la expresión 1 por la expresión 2) y la sustitución de las “velocidades” por “los momentos”.

El primer paso puede implementarse sin (excesivos) problemas en cualquier teoría. El segundo ya no. ¿Por qué? Pues porque la sustitución de las velocidades por los momentos puede no estar definida de manera única. Aparecen lo que se llaman ligaduras. Aclaro, algunos habrán estudiado ya lagrangianos y hamiltonianos en mecánica clásica (y a estas alturas estarán aburridos cómo una ostra). Allí hablan de ligaduras y de que el lagrangiano permite escribir las ecs de Newton en “coordenadas generalizadas”. Bien, esas ligaduras no son las mismas conceptualmente a las que acabo de introducir.

Vamos un poco más con las ligaduras. En mecánica clásica casi no hay sistemas con ligaduras. En cuántica relativista las hay siempre en los casos más interesantes. Y esto es porque las famosas (espero) teorías gauge son teorías con ligaduras. El electromagnetismo clásico es una teoría gauge. Y tiene una ligadura. Esta ligadura os sonará, se suele conocer cómo la ley de Gauss y se conocía mucho antes de que Dirac introdujera el concepto de ligadura para sistemas cuánticos (con intención, precisamente, de tratar la gravedad).

Volveremos más tarde con las ligaduras, pero retomemos la cuántica. Ahí deje pendiente ver que era la ecuación de Schröedinger. Pues bien, la ecuación de Schröedinger independiente del tiempo es lo que se conoce cómo una ecuación de auto-valores si se expresa en términos del Hamiltoniano, en concreto:

3 .Hφ = λ.φ ( φ denota la función de onda)

Por supuesto ahora el hamiltoniano es el hamiltoniano cuántico. Este se obtiene del clásico sustituyendo el momento clásico p por el operador momento (básicamente la derivada con respecto a la posición) y la coordenada de la partícula (que suele aparecer en el potencial) por el operador momento (en realidad ahí casi no hay cambio en la notación, aunque sí en el concepto, ahora la coordenada multiplica a la función de onda).

Esta es la ec. independiente del tiempo. La dependiente del tiempo nos dice que la derivada (parcial) de la función de onda respecto al tiempo es igual al hamiltoniano actuando sobre la función de onda.

Dije al principio que esto es para cuántica no relativista. ¿Que pasa con la cuántica no relativista? Ya anticipe algo, pero voy ahora con ello.

Por culpa del famoso E=m.c2 tenemos que si hay suficiente energía se pueden formar partículas. De hecho como existen partículas sin masa en reposo, por ejemplo los fotones, siempre pueden formarse partículas. Es decir, no hay un número fijo de partículas en el sistema. La cuántica “ordinaria” esta bien para tratar con un número fijo de partículas (en muchos casos con una sola) pero no con un número indeterminado. Por esto, y por otros motivos que no trataré, en cuántica relativsta uno ya no se ocupa de calcular funciones de onda. Ahí lo que se busca es calcular probabilidades de transición (elementos de matriz S). Es decir, la probabilidad de que si en un instante de tiempo (formalmente t en el infinito pasado) hay una serie de partículas incidentes que chocan (interaccionan) en un tiempo posterior (formalmente en infinito futuro) dan por resultado una serie de posibles partículas resultantes. Hay, normalmente, varios resultados posibles para las partículas resultantes, cada uno con su probabilidad. También hay que considerar las distribuciones angulares en que salen despedidas esas partículas. Estrictamente el estudio de choques no es exclusivamente relativista. Existe matriz S para cuántica no relativista. Pero el caso es que el formalismo matriz S es mucho mas importante en QFT. Por cierto, experimentalmente no se obtiene la matriz S sino una magnitud relacionada llamada sección eficaz de dispersión.

Bien, vamos con un aspecto muy importante. Las partículas incidentes y resultantes se corresponden con estados descritos por sistemas de partículas libres (en el infinito pasado e infinito futuro las partículas no interaccionan). En la región de transición “reina” el término de interacción. Los cálculos de las interacciones son complicados. Feynman, para ayudar a entenderlos, y también efectuar los cálculos, introdujo el concepto de partícula virtual. Este es muy importante. La idea es que un fotón y un positrón (por ejemplo) interactúan intercambiando un “fotón virtual” que es la partícula mediadora.

Bueno. un largo post y aún no he dicho nada de LQG. Muchos ya conocéis lo que he expuesto (supongo). Pero necesitaba introducir en un contexto mas sencillo una serie de conceptos que os permitirán entender mejor la LQG. Vamos ya con ello, que ya es hora.

La relatividad general, que describe la gravitación, es una teoría clásica, i.e. no cuántica. De hecho es la unica interacción que hay que no es cuántica. Eso es raro. Uno podría pensar que “el universo es así”. Y además los fenomenos cuanticos que involucren a la gravitacion son muy irreleantes a efectos prácitos (bueno, eso dicen xD). Pero tratar a la vez fenomenos cuánticos y fenomenos gravitatorios de forma clásica es inconsistente. Así que definitivamente hay que buscar una gravitación cuántica (o alguna alternativa, la que sea).

Aquí ya vamos a hacer una importante escisión. De un lado hay físicos expertos en QFT. De otro expertos en RG (relatividad general), y cada uno “barre para casa”. Los QFT, empezando por Feynman, intentaron plantear la gravedad cómo una QFT. Al igual que le fotón media la interaccion lelectromagnética (y los diversos bosones intermedios las otras interacciones) debía haber un mediador de la gravedad. El famoso gravitón. ¿Cómo obtener el gravitón?

Hay dos caminos. Por consideraciones técnicas se sabe que debe ser una partícula de spin 2 (los bosones de las otras interacciones, fotón, gluones bosones W+- y Z) son de spin 1. Si se parte de un Lagrangiano para una partícula de spin 2 en un espacio plano (sin gravedad) y se aplican ciertas condiciones de consistencia se llega a una recursión que termina llevando a las ecuaciones de Einstein (hay ahí muchas sutilezas, claro).

Otro camino es partir de las ecuaciones de Einstein. Voy a tratar este con calma porque es clave para la LQG. La relatividad especial se describe (en el formalismo de Minkowsky) por una métrica en un espacio de 4 dimensiones (la cuarta es c.t, dónde c es la velocidad de la luz yt el tiempo). Esta métrica es “plana”. Para entendernos, constante. Siempre puede reducirse (mediante un cambio de coordenadas), para cualquier punto a la forma de una matriz diagonal (-1,1,1,1). El campo gravitatorio modifica esa métrica y aunque localmente siempre puede reducirse a esa forma globalmente no. Digamos que vamos a tener una métrica dónde cada elemento de matriz es una función de las coordenadas y el tiempo.

Aclaro que es una métrica. es un objeto matemático (matriz, tensor, forma diferencial, depende del formalismo, pero para entendernos matriz nos vale) que permite medir longitudes de vectores. La longitud de un vector ordinario es |v|=sqrt(vx2 + vy 2 +vz2). Pués bien, si el vector tiene cuatro dimensiones su métrica será algo que generalizará esa expresión. y nos permite calcular su longitud. A partir de la longitud de dos vectores podemos calcular el ángulo que forman. Y midiendo ángulos podemos comprobar si vivismo en un espacio plano o en uno curvo. Vale, hasta aquí una minidivulgación de relatividad. Seguimos.

Las ecuaciones de Einstein son unas ecuaciones para la métrica. unas ecuaciones diferenciales. nos relacionan unas derivadas de la métrica, que representan la curvatura del espacio-tiempo con un término que representa la energía (y por tanto masa) de la materia. Formalmente:

4. Curvatura= Tensor energía momento.

Einstein llegó a esta ecuación mediante algo llamado “desviación geotética”, pero no nos ocupamos ahora de eso. Son, las ecuaciones de Einstein, unas ecuaciones no lineales complicadas, que implican muchos conceptos resbaladizos y etc. Pero si recordáis el principio de este post recordareis que os dije que la ec de Schröedinger vista sólo como ec perdía. Se ganaba mucho viéndola como proveniente de un hamiltoniano. Y este hamiltoniano proviene a su vez de un lagrangaino. ¿Existen lagrangianos y hamiltonianos para las ecs. de Einstein? Ahí quería yo llegar (de hecho podría haber empezado el post por aquí).

Lagrangianos si hay. Es fácil llegar a ellos. De hecho a mi me explicaron la RG partiendo de un Lagrangiano cuya expresión se obtenía por consideraciones de simetría (y “un poquito” de geometría diferencial xD). A partir de ese Laggrangiano por las ecs de Euler-Lagrange se llegaba a las ecuaciones de Einstein. Vale, lagrangianos sí, ¿y el hamiltoniano? Uff, eso es otra historia, y el corazón de parte de la LQG. Pero antes otro inciso.

Os dije que los QFT querían gravedad a partir de partículas virtuales, el gravitón. Bien, el gravitón ellos lo identificaban, en el formalismo geométrico, cómo una perturbación de la métrica. Formalmente gμν=nsub>μν + hsub>μν (g es la métrica general, n la de minkowski y h la perturbación o gravitón. La idea es expresar el lagrangiano de la RG en términos de esta división de la métrica y considerar la v cómo un campo. Este es un campo que se corresponde a una partícula de spin 2. En realidad se puede usar este formalismo considerando perturbaciones respecto a métricas curvas. El primero en introducir esta idea fué Feynman y el primero en desarrollarla formalmente fué Bryce de Witt con el método del campo de fondo.. En QFT se pueden usar tanto el lagrangiano como el hamiltoniano para obtener matriz S. De Witt obtuvo mediante esta descomposición elementos de matriz S para el campo gravitatorio, eso sí de un modo muy lioso. Mas adelante t´hoof y Veltman (y luego otros) usaron los métodos de las por entonces triunfantes teorías gauge para expresar lo mismo en formas ligeramente mas sencillas.

Un momento, ¿he dicho que de Witt obtuvo expresiones para la matriz S? Entonces ¿no acaba ahí la historia? Pues no, claro. El problema es que esas expresiones llevaban a resultados infinitos. En realidad nada nuevo, las QFT siempre llevan a resultados infinitos. Lo que pasa es que algunas QFT permiten que se eliminen esos infinitos (son renormalizables) otras no. Las no renormalizables se consideran inútiles. Y la RG cuantizada de este modo salió no renormlalizable. La gente de QFT se encontraron, trabajando con la cromodinámica cuantica, con la necesidad (o mas bien conveniencia) de cuantizar objetos unidimensinonales, digamos cuerdas. Tras una serie de vueltas rocambolescas se vió que esas cuerdas tenían por ahí un graviton y ya se lió. Treinta años de teoría de cuerdas sin ningún resultado práctico, pero eso es otra historia, que díria el cronista de Conán el bárbaro.

Los puristas de RG vieron que esa cuantización de perturbaciones a la métrica desde su perspectiva era muy poco fiel al espíritu de la RG. Así que no les sorprendía que diera una teoría inconsistente. Así pues buscaron una teoría no perturbativa. Pero antes, mucho antes, buscaron un hamiltoniano. Ahí una vez mas tenemos a Bryce de Witt (el mismo que dirigió la tesis a Everet de la que surgió la interpretación de los muchos mundos de la mecánica cuántica). Sin embargo este formalismo se conoce como formalismo ADM (Arnowitz, Desser Misner). Ellos trabajaron en este tema motivados por un aspecto interesante de la RG. Enel resto de teorías la energía es un concepto bien definido. En RG sin embargo no lo es. Una definición tentativa usa lo que se conocen como “pseudotensores” (expresiones que son tensoriales únicamente bajo un grupo restringido de transformaciones) como por ejemplo el de Landau-lipshitz, pero no son aplicables en general, y además no son nada elegantes. Se supone que el hamiltoniano debería dar una definición correcta de la energía aplicable en cualquier caso. En realidad no tuviern un éxito total, pero para casos asintóticamente minkowskianos funciono y dió lugar a lo que se llama masa ADM de un espacio-tiempo

Bien, obtener un hamiltoniano es muy compilado para la RG. Se tiene que seguir una serie de pasos un tanto sorprendentes. Lo primero es descomponer el espacio-tiempo en un producto (foliación) de espacio y tiempo (no entraré en conceptos de topología), formalmente pondré T=ExR.

Bien, lo siguiente es introducir una métrica para el espacio tridimensional curvo E (o más bien una familia de espacios, uno para cada valor de la coordenada “temporal”). Aparte hay una función de paso entre el espacio en un “instante” t1 y uno t2 que se descompone en:

5. T= N +uS

Aquí T ya no es la energía cinética. Y hay subíndices varios de por medio que omito poner. Indico simplemente que N es la función de lapso y que S es el vector de desplazamiento (de una hipersuperficie a otra).

Bien, el truco consiste introducir una métrica intrínseca y una curvatura intrínseca (referente al espacio tridimensional E) y expresar la métrica y curvatura en 4 dimensiones en términos de esta métrica y esta curvatura. La curvatura intrínseca se expresará en términos de las funciones lapso y vector desplazamiento indicados anteriormente. Pues bien, una vez hecho esto podeos plantear un lagrangiano en términos de estas nuevas variables. El lagrangiano “normal” no tiene nada claramente identificable con “velocidades” o “momentos” que nos guíen a buscar un hamiltoniano. El lagrangiano expresado en estas variables sí (no puedo entrar en detalles de que es cada cosa por motivos de las limitaciones para escribir fórmulas, y además tampoco aportaría mucho teniendo en cuanta que lo que estoy viendo ahora es un paso intermedio hacia la LQG)

Pero claro, no todo iba a ser tan fácil (cómo si hasta aquí hubiera sido fácil xD). Este lagrangiano es singular. Es decir, que los momentos no pueden despejarse a partir de las velocidades. Hay ligaduras. Pero con ligaduras o sin ellas puede haber un hamiltoniano. el hamiltoniano contienen ligaduras.

Ahora cuento una cosa que no conté en su momento. Un hamiltoniano con ligaduras tiene dos tipos de ecuaciones. Las normales hamilton y las de ligadura. Como, eso sí, dije las ecuaciones de ligadura para el campo electromagnético es la ley de gauss. Esas ecuaciones de ligadura relacionan estados físicos equivalentes. La otra parte del hamiltoninao que no es ligadura, dicta la evolución temporal de los estados físicos.

Pues bien, en RG tenemos un “pequeño problemilla”. Resulta que el hamiltoniano no tiene ni una sola parte que no sea de ligadura. Por tanto no hay evolución temporal.

Bien, ya casi tenemos planteado el quid de la cuestión. El hamiltoniano este tiene dos ecuaciones de ligadura. una viene a representar la invariancia bajo cambio de coordenadas de la RG (invariancia bajo difeomorfismos dicho de otro modo).Es una ecuación relativamente sencilla. La otra ec. de ligadura es mucho mas compleja, es no polinomial. Y tiene nombre propio. Se llama ecuación de Wheeler de Whitt. Viene a representar la “física” del sistema.

Esta ecuación es dificilísima de resolver. Sólo en casos muy sencillos hay soluciones aproximadas. Son lo que se conocen como modelos de “miniespacio” o “minisuperespacio” (nada que ver con la supersimetría pese al nombre).

Bien, antes de exponer mas cosas voy con el Lagrangiano. Para ser exactos con el lagrangiano “normal”. No he dicho nada pero el tratamiento de la QFT basado en lagrangianos parte de la integral de caminos de Feynman. La idea intuitiva es que una partícula cuántica recorre todas las trayectorias entre un punto y otro. Cada trayectoria contribuye a la integral de camino con un factor de exponencial imaginaria, e-iS dónde S es la acción (la integral del lagrangiano a lo largo de la trayectoria). De ese modo los estados que cumplen las ecuaciones clásicas son los que más contribuyen-los siguientes que mas contribuyen, si los hay, son los instantones-). Bien, esto para partículas puntuales. Un campo pasa por todas sus configuraciones posibles y etc. Los detalles son terribles, claro. Aunque conceptualmente delicada y matemáticamente no muy bien definida la integral de caminos es muy potente para generar elementos de matriz S. Las teorías de campos gauge se cuantizan normalmente en formalismo lagrangiano (aunque también pueden cuantizarse con hamiltonianos) + integral de caminos.

Bien ¿puede aplicarse la integral de caminos a la RG? Ya dije que en forma perturbativa es lo que se había hecho. Pero Hawking usó una pequeña variante, lo que se conoce como euclidean quantum gravity en que rota a un tiempo imaginario. La integral de caminos es casi incalculable para la RG, pero para ciertas configuraciones sencillas puede calcularse en una aproximación conocida cómo aproximación WKB (wenzel krammers brillouin). Y, curiosamente, esta aproximación nos lleva a una ecuación que, “curiosamente”, es la ecuación de Wheller de Witt.

Nota, hasta ahora no he dado ninguna bibliografía, creo que es hora de recomendar algo. Para el tema de la euclidean quantum gravity hay libros enteros, pero no me parece que merezca la pena leerlos pués no parece tener un gran futuro. Pero hay un librito que recoge unas conferencias de Hawking y Penrose que tratan en parte estos temas y que es fino, i.e. no muy extenso, y muy agradable de leer. Se llama “la naturaleza del espacio y el tiempo”. La formulación ADM de la gravedad viene explicad en los apéndices del famoso Wald De relatividad general.

Bien, pués hasta aquí por ahora. He introducido una serie de conceptos que sitúan el origen de la LQG, plantean las dificultades que había en un momento dado. En otro futuro post indicaré cómo desde la via hamiltoniana se introducen unas nuevas variables (en vez de la métrica la conexión) que dan un hamiltoninao de una estructura de ligaduras más sencilla que puede resolverse en términos de funciones cilíndricas la primera y de spin networks la segunda (y explicaré como los lazos, willson loops, también resolvían esas ligaduras en versiones mas antiguas de la LQG). Explicaré cómo ese formalismo permite deducir que el área esta cuantizada. Explicaré que esa cuantizacion permite calcular la entropía de un agujero negro.

También comentaré que en casos sencillos (simetrías) hay variantes mas simples de la LQG que se aplican a la cosmología (loop quantum cosmology) y a las singularidades de los agujeros negros. Explicaré el problema del tiempo en LQG canónica. Y explicaré que ese problema, entre otros motivos, lleva a otro punto de partida, las spin foams, que son una especie de versión lagrangian y de integral de caminos de la LQG canónica (basada en hamiltonianos). Pero comentaré que los lagrangianos allí usados no son exactamente el de Einstein sino “más o menos equivalentes” y que hay varias aproximaciones.

También explicaré que hay el problema del límite clásico. No se puede a partir del hamiltoninao llegar a un gravitón a un espacio de Minkowsky,y por tanto conectar con la RG clásica. En spin foams hay un cálculo de un propagador de un gravitón, que sí daría esa conexión.

Y por supuesto explicaré las críticas que les hacen los de cuerdas a la LQG.


Seguir

Recibe cada nueva publicación en tu buzón de correo electrónico.