Una de las predicciones más intrigantes de la cosmología nacida del Big Bang son los defectos topológicos del espacio-tiempo, en particular las cuerdas cósmicas, capaces de dejar una huella en la radiación de fondo cósmico. En un universo en expansión, estos objetos, surgidos durante transiciones de fase, se convierten en fuentes potentísimas de ondas gravitacionales, creando un fondo estocástico que detectores como LIGO, construido entre otros por Rainer Weiss y Kip Thorne, intentan registrar. Sin embargo, en las teorías de gran unificación las cuerdas pueden ser metaestables: se desintegran produciendo espontáneamente pares de monopolos magnéticos. Este proceso — el efecto túnel cuántico en el marco de la teoría cuántica de campos — suprime la señal de ondas gravitacionales a bajas frecuencias y depende críticamente del valor exacto de la acción de rebote, que antes se calculaba solo en la aproximación de cuerda infinitesimalmente delgada.
Los autores abandonaron las simplificaciones analíticas y acudieron a métodos numéricos. Consideraron un modelo mínimo con dos etapas de ruptura de simetría SU(2) → U(1) → nada, donde primero nacen los monopolos y luego las cuerdas. Para calcular la acción de rebote —la cantidad central que determina la tasa de desintegración— se empleó el método del flujo de gradiente en una red espaciotemporal. En lugar de buscar directamente el punto de silla del funcional de acción, lo cual se complica por la presencia de un modo negativo, las simulaciones buscaban el mínimo de la acción con la posición del centro del monopolo fijada (ρ_E = R) y luego maximizaban la acción respecto a R. Los perfiles de los campos de bosones gauge y escalares de Higgs se describieron mediante un ansatz general restringido por simetría, sin recurrir a suposiciones adicionales. Las simulaciones en la red se realizaron con distintos pasos para extrapolar los resultados al límite continuo.
Los experimentos numéricos revelaron una reducción significativa de la acción de rebote en comparación con la fórmula clásica S_B^(PV) = π m_M^2/μ. Para un ejemplo ilustrativo con masas iguales de partículas de Higgs y de gauge en cada uno de los dos sectores, la acción de rebote cae a valores por debajo de 60 (en unidades de m_W^2/(g^2 m_γ^2)) para una razón de masas m_W/m_γ ≈ 3, mientras que el límite de cuerda delgada da alrededor de 60.5. La extrapolación al límite continuo permitió construir una función de ajuste fiable que describe S_B con alta precisión en todo el rango de parámetros. Para jerarquías grandes (m_W/m_γ ≫ 1) los resultados convergen suavemente al límite de cuerda delgada, confirmando la autoconsistencia del método. Esto significa que la tasa de desintegración Γ ∝ exp(-S_B) aumenta en órdenes de magnitud, y las cuerdas cósmicas viven menos de lo que se suponía. Como consecuencia, la frecuencia característica f_low, por debajo de la cual el espectro de ondas gravitacionales se corta, se desplaza hacia arriba: de los nanohercios a los micro- y milihercios.
En la práctica, esto se traduce en una revisión del fondo de ondas gravitacionales esperado. Regiones de parámetros que antes se consideraban seguras para los detectores terrestres ahora pueden caer dentro de la zona de exclusión de LIGO/Virgo. Por otro lado, la señal detectada por las colaboraciones de cronometraje de púlsares (PTA), como NANOGrav, requiere una interpretación más precisa: para explicarla con cuerdas metaestables, es necesario aumentar la razón de escalas de ruptura de simetría, lo cual facilita la construcción de modelos consistentes de gran unificación con inflación intermedia. En otras palabras, tener en cuenta el grosor finito de las cuerdas hace que la física a la escala de gran unificación sea más accesible a la comprobación experimental.
El método propuesto de simulación en red puede generalizarse a relaciones arbitrarias de masas de las partículas que surgen en la ruptura de simetría. El siguiente paso será incluir los efectos relacionados con inhomogeneidades en las cuerdas, así como la creación térmica de monopolos. Esto permitirá predecir con mayor precisión la forma del espectro de ondas gravitacionales y, posiblemente, explicar el reciente indicio de señal en PTA sin conflictos con las restricciones de LIGO.
Los resultados tendrán impacto en la astronomía de ondas gravitacionales, la cosmología observacional y la física de partículas elementales, ya que vinculan datos experimentales concretos con los parámetros de las teorías de gran unificación.
La tarea inmediata es aplicar este enfoque a grupos de gauge más realistas, como SU(2)×U(1), e incluir las correcciones relativistas al movimiento de las cuerdas.
El trabajo aborda cuestiones fundamentales en la intersección de la teoría cuántica de campos y la cosmología, incluyendo el problema de la desintegración del falso vacío en un espacio-tiempo curvado. Además, proporciona una nueva herramienta para interpretar la enigmática señal de los conjuntos de cronometraje de púlsares, basados en la observación de estrellas de neutrones: si se trata del primer eco de la física de cuerdas o si tiene un origen astrofísico, es algo que se deberá dilucidar en los próximos años.
🎯 Una cuerda cósmica con un grosor de milmillonésimas del tamaño de un protón puede tener una tensión equivalente a la masa de un millón de toneladas por cada centímetro de longitud; es comparable a la masa de un rascacielos comprimida en una regla.