En una habitación oscura suena música. En el suelo se ha esparcido polvo fino, casi invisible. Las bajas frecuencias hacen que se agrupe en patrones caprichosos: filamentos, cúmulos, vacíos. De manera similar, la gravedad ha ido tejiendo durante miles de millones de años las galaxias en la telaraña cósmica, solo que en lugar de ondas sonoras actúa la atracción. Para entender con qué "melodía" comenzó todo, hay que retroceder en el tiempo: volver del patrón ya formado al puñado inicial de polvo. Esta es la tarea de reconstruir las condiciones iniciales del universo.
Las simulaciones informáticas convencionales son demasiado pesadas para ese viaje hacia atrás. Calculan honestamente cada empujón gravitatorio, pero una iteración puede llevar horas. Los métodos aproximados son más rápidos, pero recortan esquinas: pierden esa dinámica de encaje donde nacen los halos. Los autores del trabajo proponen un compromiso: un autómata celular con cerebro de red neuronal. El modelo LNCA (autómata neurocelular lagrangiano) no sigue una rejilla fija de densidad, sino las propias partículas de materia y se mueve con ellas. Solo aprende a corregir una conjetura burda pero físicamente significativa: la aproximación de Zeldóvich. Sobre esta base, el modelo añade desplazamientos residuales cuidadosos. Así se mantiene la precisión a gran escala, y el aprendizaje dibuja los detalles finos.
El secreto del ahorro está en la localidad. Cada célula del autómata actualiza su estado mirando solo a sus vecinos más cercanos, y lo hace paso a paso. El resultado no es solo una instantánea final, sino toda la historia del movimiento, como una filmación cuadro a cuadro del baile de las partículas de polvo. Las reglas de este baile se rotan y se desplazan junto con la materia: la arquitectura posee equivariancia rotacional y traslacional. La precisión final es impresionante. Hasta escalas de unos 12 megapársecs (aproximadamente 40 millones de años luz), el modelo no se equivoca más que en unos pocos por cien. Es un péndulo fiable para retroceder en el tiempo, desde la telaraña actual hasta un comienzo casi homogéneo.
tal como lo escribieron los autores, en inglés arXiv:2607.27320
Field-level inference of cosmological initial conditions from galaxy surveys requires a forward model that is simultaneously accurate in the non-linear regime, computationally efficient, and fully differentiable. Traditional N-body simulations are accurate but computationally prohibitive for iterative inference, while approximate solvers like Lagrangian Perturbation Theory (LPT) fail to capture the knotty halo-forming dynamics of the cosmic web at late times. We introduce the \textit{Lagrangian Neural Cellular Automaton} (LNCA), a hybrid deep learning framework that can be applied to emulate structure formation as a local, iterative dynamical process on a comoving lattice. Unlike standard Eulerian Convolutional Neural Networks (CNNs) which map fixed density fields, the LNCA operates in the Lagrangian frame, advecting the computational graph itself to follow the flow of mass. By training the network to learn only the \textit{residual} displacement corrections to the Zeldovich approximation, we achieve high-fidelity emulation of the non-linear physics while guaranteeing accuracy at large scales. We further constrain our model to produce complete trajectories, not just final states, by adopting an equivariant cellular automaton architecture, which recurrently iterates on its internal states to yield a dynamic history. The resulting model is strictly local, translationally and rotationally equivariant, and naturally supports continuous time integration, making it a reliable differentiable forward model for reconstructing the initial conditions of the universe from lightcone data. Our trained model supports percent-level precision in the power and cross spectra well into the non-linear regime ($$k \lesssim 0.5 \, h \text{Mpc}^{-1}$$), while requiring $$\sim10^4$$ times fewer learned parameters than comparable models which take the form of an interpretable internal dynamic rule set.
La aproximación de Zeldóvich, que el modelo toma como base, fue ideada en 1970 para describir el crecimiento de las perturbaciones en un universo en expansión. Es tan acertada que todavía sirve como "aproximación de orden cero" en muchos cálculos cosmológicos.
En la novela "El fin de la eternidad", Isaac Asimov describe una organización que corrige la historia calculando las consecuencias de los cambios más mínimos. LNCA es un paso hacia ese "computador de la realidad", capaz de reconstruir las condiciones iniciales a partir del patrón actual de galaxias.