For nearly a decade, astronomers believed they had solved one of the great mysteries of cosmic alchemy, attributing the production of the universe’s heaviest elements, such as gold and platinum, primarily to kilonovae from colliding neutron stars. This consensus was further solidified by the landmark 2017 detection of gravitational waves from such a merger.

 

However, new research from Los Alamos National Laboratory, published in The Astrophysical Journal Letters, challenges this narrative. A team led by Marko Ristić demonstrates through advanced supercomputing simulations that long-duration gamma-ray bursts, some of the most energetic explosions in existence, can produce kilonova-like signatures without requiring a neutron-star merger. Instead, the team proposes that collapsing massive stars, or "collapsars," can generate these characteristic optical and infrared emissions via a previously overlooked nucleosynthesis mechanism within their relativistic jets. This discovery is more than a mere astrophysical curiosity; it highlights how modern high-performance computing is fundamentally transforming our understanding of the cosmic origins of the elements.

Gamma-ray bursts (GRBs) are brief flashes of extraordinarily energetic radiation capable of releasing more energy in seconds than the Sun will emit during its entire lifetime. For decades, astronomers divided GRBs into two categories: short bursts produced by compact-object mergers and long bursts generated by collapsing massive stars. Observations of GRB 211211A and GRB 230307A complicated that picture. Although both events lasted roughly 40 seconds, far longer than typical merger-driven bursts, they exhibited infrared signatures resembling kilonovae, leading many researchers to conclude they originated from neutron-star mergers. The Los Alamos team questioned that assumption.

 

Their work proposes that the observed emission can be reproduced by a collapsar, a rapidly rotating massive star collapsing into a black hole while launching relativistic jets through its interior. Rather than producing heavy elements through tidal disruption of neutron stars, the model creates neutron-rich conditions within the jet and the surrounding cocoon.

Testing that idea required a computational effort spanning multiple physics domains. Researchers combined nuclear reaction networks, magnetohydrodynamic simulations, Bayesian inference techniques, radiative transfer calculations, and large-scale parameter exploration. The project united scientists from Los Alamos' Theoretical Division, Computational Division, and Center for Theoretical Astrophysics.

 

At the heart of the study was the Portable Routines for Integrated nucleoSynthesis Modeling (PRISM) framework, which simulated the creation of heavy elements through rapid neutron-capture processes. The calculations explored how intense photon fields inside collapsar jets generate neutrons that subsequently seed nucleosynthesis in the surrounding cocoon. The researchers modeled both weak and full r-process scenarios and calculated the resulting radioactive heating over time. Those outputs became inputs for another computationally intensive stage: radiation transport simulations.

To predict what astronomers should observe, the team employed Los Alamos' SuperNu code, a sophisticated Monte Carlo radiation transport framework widely used in transient astrophysics. SuperNu follows millions of photon packets as they interact with expanding ejecta, accounting for absorption, emission, scattering, radioactive heating, and detailed atomic opacities. The simulations used high-fidelity opacity tables generated by Los Alamos atomic physics codes and modeled spectra across wavelengths ranging from ultraviolet through infrared. The computational workflow produced synthetic observations that could be directly compared with telescope data from GRB 211211A and GRB 230307A.

 

Rather than relying on simplified analytical approximations, the researchers simulated the detailed microphysics governing how light emerges from expanding stellar debris. The result was remarkable. A single weak-r-process ejecta component reproduced both the optical and infrared observations associated with the two gamma-ray bursts.

The study also demonstrates how artificial intelligence and statistical computing are becoming essential tools for modern astrophysics. The team generated dozens of high-fidelity SuperNu simulations across a wide range of ejecta masses and velocities. Because running radiation transport calculations for every possible parameter combination would be prohibitively expensive, researchers trained Gaussian Process surrogate models to emulate the simulation outputs. These surrogate models enabled rapid Bayesian inference using the RIFT parameter-estimation framework, allowing the team to explore vast parameter spaces while preserving the accuracy of the underlying simulations.

 

This combination of supercomputing and machine-learning acceleration represents an increasingly common pattern across computational science, where advanced simulations generate data and AI-driven techniques help scientists navigate the resulting complexity.

One of the study's most innovative computational components appears in its appendix. The researchers developed a three-fluid, two-dimensional magnetohydrodynamic simulation that tracks protons, neutrons, and alpha particles inside collapsar jets. The simulation investigated a phenomenon the team calls a magnetic "sieve." Strong magnetic fields trap charged particles near the jet core while allowing neutral neutrons to migrate into the surrounding cocoon. Under sufficiently intense magnetic fields, approaching 10¹² gauss, the resulting environment becomes neutron-rich enough to sustain rapid neutron-capture nucleosynthesis.

 

The simulation solved coupled continuity, momentum, energy, transport, and magnetic induction equations using an implicit-explicit numerical approach designed for stiff plasma systems. Without modern high-performance computing resources, such calculations would be effectively impossible.

Perhaps the most surprising conclusion is that a red kilonova does not necessarily imply the creation of the heaviest r-process elements. The team's simulations show that a weak r-process producing elements only up to approximately mass number 130 can generate the observed red evolution traditionally interpreted as evidence for lanthanide production. In other words, astronomers may need to reconsider some long-standing assumptions about how heavy elements are identified in explosive transients. If confirmed, the discovery could reshape our understanding of how the universe manufactures many of its elements.

The broader significance of the work extends well beyond gamma-ray bursts. The study exemplifies how supercomputing has become a primary scientific instrument. Researchers combined nuclear physics, plasma physics, radiation transport, machine learning, Bayesian statistics, and astrophysical modeling into a single computational framework capable of testing ideas that cannot be reproduced in any terrestrial laboratory. The calculations relied on resources provided through Los Alamos National Laboratory's Institutional Computing Program, underscoring the increasingly central role of advanced computing infrastructure in modern astrophysical discovery.

 

A generation ago, astronomers could only observe the universe. Today, they can recreate it. And as this Los Alamos research demonstrates, the next breakthrough in understanding the origin of the elements may emerge not from a telescope alone, but from the convergence of supercomputing, artificial intelligence, and computational physics operating at unprecedented scale. The universe still holds many secrets. Increasingly, supercomputers are becoming the tools that allow us to uncover them.