Entrada autorreferencial: Recopilación de los Ig Nobel de 2008 en la Mula Francis

Los Premios Nobel Ig de 2009 se anunciarán en la ceremonia que se celebrará el próximo jueves, 1 de octubre, a las 7:30 PM, hora de la costa este norteamericana. Dedicaré en este blog una entrada a los premios más interesantes, como ya hice el año pasado. Para abrir boca, os recuerdo las entradas del año pasado (quizás algún lector reciente no las haya leído o algún otro quiera releerlas).

2008 Ig Nobel, Paz: La dignidad y los sentimientos de nuestras amigas las plantas, 8 Octubre 2008

2008 Ig Nobel, Arqueología: El armadillo “arqueólogo” como pesadilla para el arqueólogo humano, 8 Octubre 2008

2008 Ig Nobel, Medicina: La importancia del precio de un medicamento en su eficacia terapéutica, 7 Octubre 2008

2008 Ig Nobel, Ciencias cognitivas: Hasta una protozoo puede encontrar el camino más corto en un laberinto, 6 Octubre 2008

2008 Ig Nobel, Economía: las bailarinas topless deberían evitar la píldora para lograr mayores propinas, 6 Octubre 2008

2008 Ig Nobel, Física: El enredo espontáneo de lo enredable o no le dejes el ovillo de lana al gato, 5 Octubre 2008

2008 Ig Nobel, Biología: Las pulgas de los perros son mejores saltadoras que las pulgas de los gatos, 5 Octubre 2008

2008 Ig Nobel, Literatura: “Eres un bastardo” o el acaso laboral en el trabajo, 5 Octubre 2008

2008 Ig Nobel, Química: Coca-Cola y Pepsi-Cola, los anticonceptivos del pobre, 4 Octubre 2008

2008 Ig Nobel, Nutrición: El sonido “crujiente” y el disfrute de una patata frita, 3 Octubre 2008

El problema de Thomson o como distribuir N cargas en una esfera con energía total mínima

El físico británico J. J. Thomson ganó el Premio Nobel en 1906 por el descubrimiento del electrón. En 1904 propuso un problema matemático muy difícil de resolver en general: ¿cuál es la configuración de mínima energía para N electrones (con una fuerza repulsiva 1/r2) en una superficie esférica? Para N pequeño obtener la solución óptima a este problema es fácil. Para 4, 6, y 12 corresponden a los vértices de un sólido platónico. Hasta N=400 se conocen las soluciones óptimas. Sin embargo, para N>400 sólo se conocen algunas pocas, el resto son sólo las mejores candidatos obtenidas por ordenador. Wales, McKay y Altschuler han obtenido por simulación las mejores configuraciones hasta el momento en el rango N de 400 a 4000. La vídeo muestra cinco de los nuevos resultados para N=400, 752, 1632, 3952, y 4352. Los tres primeros son configuraciones simétricas. Los dos últimos son configuraciones asimétricas ligeramente de menor energía que las simétricas observadas. ¿Serán óptimas? Nadie lo sabe pero la búsqueda de la demostración por ordenador continúa. Nos lo cuenta Tony Phillips, “Progress on the Thomson problem,” Take on Math in the Media, September, 2009.

Las 400 configuraciones óptimas fueron publicadas en David J. Wales, Sidika Ulker, “Structure and Dynamics of Spherical Crystals Characterised for the Thomson Problem,” Phys. Rev. B, 74, 212101 (2006) [Página con animaciones gif con los resultados óptimos]. El nuevo artículo con las soluciones cuasi-óptimas es David J. Wales, Hayley McKay, Eric L. Altschuler, “Defect motifs for spherical topologies,” Phys. Rev. B 79, 224115 (2009) [Página con animaciones gif con los resultados cuasi-óptimos]. 

En una configuración de mínima energía cada electrón está rodeado de 6 vecinos cercanos (hexágonos verdes en el vídeo) resultando en una carga efectiva nula, excepto ciertos electrones que están rodeados de 5 vecinos (pentágonos rojos) con una carga efectiva de +1, o de 7 vecinos (heptágonos azules) con una carga efectiva de -1. Siendo Ci el número de los vecinos cercanos al electrón i-ésimo, la red que conecta los electrones más cercanos entre sí define una triangulación de la superficie de la esfera con V=N vértices, E = (1/2) Σi Ci aristas y F = 2 E/3 caras. El teorema de Euler, V-E+F=2, aplicado a esta tringulación nos da Σi (6-Ci) = 12, es decir, la suma de las cargas efectivas debe ser igual a 12. El problema de optimización a resolver es dónde hay que colocar las cargas efectivas para minimizar la energía total.

Como se observa en el vídeo, para las configuraciones con N = 400, 752, 1632, 3952, y 4352, conforme N crece, el número de heptágonos también crece. Las tres primeras configuraciones son (aproximadamente) simétricas, con una simetría icosaédrica aproximada que en algunos casos, como para N=1632, es exacta. Lo más sorprendente es que en muchos casos, como los dos últimos ilustrados en el vídeo, las configuraciones simétricas no son siempre las de menor energía. La energía potencial se define geométricamente como P = Σi>j |ri – rj|-1, donde se ha representado el electrón i-ésimo con un vector unitario ri in R3. Por ejemplo, para N=4352 la configuración cuasi-óptima tiene energía potencial P = 9311276, mientras que la configuración con 12 rosetas colocadas simétricamente (algo parecido a la configuración con N=1632 del vídeo) sólo alcanza un valor de P = 9311299. Es decir, una configuración más simétrica cercana a la mejor es sólo ligeramente peor. Realmente sorprendente.