The reason this phase had remained hidden for so long lies in the architecture of previous X-ray missions. Most narrow-field X-ray telescopes could not find transient events on their own; they relied on gamma-ray detections to provide an initial localization before repointing their instruments toward the source. By the time they arrived, the earliest and softest X-ray emission had already passed. Einstein Probe, with its wide-field soft X-ray monitoring capability, was watching the sky at the right moment and in the right energy band, providing a direct observation of the earliest stage of the short gamma-ray burst explosion. The finding, published in Science Bulletin, suggests that the standard picture of these mergers, in which the action is essentially over within a couple of seconds, is incomplete.

An Li, a PhD student at Beijing Normal University who serves as a Transient Advocate for Einstein Probe, was on shift when the event unfolded. According to Li, the event initially appeared to be an ordinary short gamma-ray burst, producing a bright flash lasting less than half a second that was detected simultaneously in gamma rays by SVOM-GRM and Insight-HXMT and in X-rays by EP-WXT. He swiftly responded to the onboard alerts and conducted a preliminary analysis. Rather than disappearing, the source kept producing episodes of soft X-ray emission for nearly ten minutes, an behavior that immediately set the event apart from the hundreds of short bursts recorded by earlier missions.

The energy budget of this extended emission was substantial, but its spectral character made it extraordinarily difficult to detect with conventional instrumentation. Professor Bin-Bin Zhang of Nanjing University, a co-corresponding author of the paper who initiated the in-depth study, explained that the spectrum was so soft that, for a burst at a typical cosmological distance, it would have fallen below the detection threshold of conventional gamma-ray instruments such as Swift’s Burst Alert Telescope. Previous missions would therefore have recorded only the brief gamma-ray flash and missed the prolonged activity that Einstein Probe revealed. As Zhang put it, the observations show that what appears to be a typical short gamma-ray burst can actually conceal a much longer and richer episode of activity at soft X-ray energies.

Determining the physical nature of the explosion required more than a single instrument. The team mounted an extensive international follow-up campaign spanning X-ray, optical, and radio wavelengths, coordinating facilities across the globe and in orbit. Professor Eleonora Troja of the University of Rome Tor Vergata, a co-corresponding author whose group obtained the key redshift information from spectroscopic analysis, emphasized that the coordinated multiwavelength observations were essential. They allowed the team to identify and study the burst’s host galaxy, measure its distance, and, critically, rule out an accompanying supernova. That exclusion provided strong evidence linking the extraordinary X-ray emission to a compact object merger rather than the death of a massive star, which produces long gamma-ray bursts and supernovae.

The deeper analysis delivered perhaps the most consequential result of the study: the long-lasting X-ray emission appeared to be powered directly by the merger remnant itself, not by the expanding blast wave slamming into surrounding material. Yi-Han Iris Yin, a PhD student in the Department of Physics and the Hong Kong Institute of Astronomy and Astrophysics at The University of Hong Kong, led the analysis of the high-energy emission. As a co-corresponding author, she found that the event’s rapid variability, its spectral evolution, and the subsequent behavior of both the X-ray and optical afterglows all pointed toward sustained activity from a central engine operating long after the initial short gamma-ray burst had faded. In standard afterglow models, emission is dominated by the external shock; here, the data told a different story.

One plausible explanation advanced by the team is that the merger produced a rapidly rotating, highly magnetized neutron star, known as a magnetar, which powered the extended X-ray emission through continued injection of energy. A stable magnetar remnant of this kind would rotate hundreds of times per second, winding up extreme magnetic fields that can tap the enormous rotational energy of the newborn object and channel it into electromagnetic radiation. Such a remnant has major implications for fundamental physics, because whether a merger produces a long-lived neutron star or collapses promptly into a black hole depends sensitively on the true maximum mass of neutron stars, which in turn depends on the poorly constrained equation of state of matter at nuclear densities.

The discovery also resonates strongly with the field of multi-messenger astronomy. Since the first joint detection of electromagnetic signals and gravitational waves from merging neutron stars in 2017, astronomers have been searching for electromagnetic counterparts to gravitational-wave sources that can reveal what happens during and after these violent cosmic collisions. Troja noted that the newly discovered soft X-ray component provides a new probe, indicating that fast X-ray transients are also electromagnetic counterparts to gravitational-wave sources and may originate from compact object mergers. If confirmed by future joint detections, fast X-ray transients could join kilonovae as standard signposts of mergers, helping telescopes and gravitational-wave detectors triangulate events across the sky.

Importantly, the researchers argue that this phenomenon may not be rare at all. Similar soft X-ray emission could accompany many more short gamma-ray bursts, but it may have escaped detection simply because previous missions lacked the capability to capture prompt emission below gamma-ray energies. In that sense, EP250704a may be the first clearly observed member of a population that has been hiding in plain sight, its soft X-ray glow too faint and too long-lived for gamma-ray monitors, and too prompt for repointing X-ray telescopes. Wide-field, sensitive soft X-ray monitoring changes the observational calculus entirely, potentially multiplying the number of merger events whose central engines can be studied directly.

The broader stakes of the finding extend to some of the deepest questions in astrophysics. Yin observed that the discovery extends our view of neutron star mergers beyond the brief gamma-ray flash, and that by revealing this previously hidden soft X-ray phase, Einstein Probe opens a new window for studying neutron star merger remnants and may ultimately help constrain the neutron star equation of state. Zhang added that the findings demonstrate Einstein Probe’s unique capability to uncover new classes of transient phenomena and strengthen its role in the era of multi-messenger astronomy, in which gravitational waves and electromagnetic radiation are studied together to understand some of the most extreme events in the Universe. As more wide-field X-ray monitors come online and gravitational-wave detectors grow more sensitive, events like EP250704a are likely to become a routine part of the astronomer’s toolkit, transforming short gamma-ray bursts from split-second flashes into extended laboratories for probing matter under the most extreme conditions nature allows.

Subject of Research: Soft X-ray emission from a neutron star merger associated with a short gamma-ray burst

Article Title: A hidden X-ray chapter of cosmic collisions comes to light

Article References: A hidden X-ray chapter of cosmic collisions comes to light. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: Einstein Probe, short gamma-ray burst, neutron star merger, magnetar, soft X-ray transient, gravitational waves, multi-messenger astronomy, central engine, kilonova, SVOM, Insight-HXMT, equation of state

The new work, published in The European Physical Journal C by Kazem Bitaghsir Fadafan of Shahrood University of Technology, Jesús Cruz Rojas of the National Autonomous University of Mexico, and Jonas Mager of the Vienna University of Technology, applies holography, the correspondence that links strongly coupled quantum field theories to weakly coupled gravitational theories in higher dimensions. The central obstacle in describing quark matter inside neutron stars is that the relevant regime of quantum chromodynamics, or QCD, is strongly coupled. First-principles lattice simulations fail there because of the notorious fermion sign problem at large densities, while traditional phenomenological approaches such as nucleon effective field theories struggle under the extreme conditions found in stellar cores. Holographic models offer a way around this impasse: instead of computing the equation of state of dense quark matter directly, one solves a comparatively tractable problem in a higher-dimensional gravitational dual and reads off the properties of the quantum field theory from the geometry.

The researchers built their model on the well-studied D3/D7-brane configuration of type IIB string theory, a bottom-up construction in which D7-branes embedded in the curved background of D3-branes describe quarks moving in a gauge theory. The key innovation lies in the dilaton, a scalar field that controls the running of the gauge coupling in the dual theory. Rather than keeping the dilaton constant, as in the simplest version of the model, the team adopted a phenomenological profile that interpolates smoothly and monotonically from its ultraviolet value to a finite, regular value in the infrared. Three parameters, denoted A, lambda and kappa, govern the precise shape of this profile, controlling the infrared value of the dilaton, the scale at which conformal symmetry is broken, and the steepness of the transition. This tailored profile captures essential features of QCD such as chiral symmetry breaking and produces a deconfined yet massive quark phase at finite density.

In the holographic setup, the embedding function of the D7-brane encodes the quark mass and condensate, while a gauge field living on the brane describes the quark chemical potential and density. Two distinct phases emerge: a Minkowski, or vacuum, phase with broken chiral symmetry and zero baryon density, and a quark phase in which chiral symmetry is restored and the density is nonzero. By computing the on-shell action of the brane, which the holographic dictionary identifies with the grand canonical potential, the researchers obtained the pressure as a function of chemical potential and hence the equation of state of the quark phase. At asymptotically large densities the model correctly reproduces the perturbative QCD behavior, with the pressure scaling as the fourth power of the chemical potential.

The team also attempted to describe the baryonic phase within the same holographic framework, modeling nucleons as D5-branes wrapped on a five-sphere. The nontrivial dilaton profile stabilizes these wrapped branes against collapse, which is itself a notable achievement. However, in the homogeneous, smeared approximation used here, the resulting nuclear equation of state turned out to be unrealistic, predicting pressures in sharp conflict with phenomenology just above the onset of the baryonic phase. The authors showed analytically that near the transition the chemical potential scales as the one-third power of the density, forcing a fourth-order transition that cannot reproduce the first-order behavior expected of isospin-symmetric QCD. They therefore discarded the holographic baryon phase and instead adopted the phenomenological equations of state of Hebeler and collaborators, derived from chiral effective field theory and constrained by nuclear physics and observation.

With the baryonic phase fixed phenomenologically and the quark phase supplied by holography, the researchers scanned the parameter space of the model, varying the AdS radius and the t’Hooft coupling. They found a region of parameters, with the AdS radius between roughly 0.015 and 0.02 inverse megaelectronvolts and the t’Hooft coupling between 1.9 and 3, that produces rather stiff equations of state with only weakly first-order transitions between baryonic and quark matter. This softness of the transition is crucial: it allows a smooth conversion from nuclear matter to quark matter inside a star rather than a violent discontinuity that would destabilize the configuration. When combined with the stiff Hebeler equation of state, three parameter choices produced equations of state supporting stable quark cores, with phase transitions occurring at chemical potentials around 400 megaelectronvolts.

To connect the microscopic equation of state to observable stellar properties, the team solved the Tolman-Oppenheimer-Volkoff equations of general relativity, which determine the structure of non-rotating compact stars. The resulting mass-radius curves show stars that begin their lives as ordinary nucleonic objects and develop quark cores as the central density increases. Stability requires that the mass increase with central energy density, and the stars with quark cores satisfy this criterion all the way up to the maximum of the mass-radius curve, beyond which a radial mode becomes unstable and the star collapses into a black hole. The maximum masses obtained for the three quark-star-supporting parameter choices are 2.17, 2.13 and 1.90 solar masses, the largest of which is broadly consistent with the heaviest precisely measured neutron star mass of about 2.35 solar masses, albeit somewhat favoring lower values.

The authors also computed the tidal deformability, the quantity that measures how easily a star is deformed by the gravitational field of a companion and that leaves a characteristic imprint on gravitational wave signals from binary neutron star inspirals. Once a star becomes heavy enough to develop a quark core, the tidal deformability drops rapidly, falling by roughly 80 units between the onset of the quark phase and the maximum allowed mass. This rapid decrease is a potentially observable signature of quark matter formation. However, the very stiffness of the baryonic phase required to obtain stable quark cores causes the model to overshoot the tidal deformability constraint from the GW170817 event at 1.4 solar masses, predicting a value near 950 compared with the observational bound of about 190 with large uncertainties. The authors note that a baryonic equation of state intermediate between the stiff and medium phenomenological cases would likely ease this tension while still permitting stable quark cores.

Several other results add nuance to the picture. The speed of sound in the quark phase remains below the speed of light, although the stiff baryonic phase violates a recently derived transport bound. More strikingly, the polytropic index of the holographic quark phase exceeds the value of 1.75 that has been proposed in the literature as a criterion for the onset of quark matter, reaching values as high as 2.5 just after the transition. This suggests that the criterion may be too restrictive, and it contrasts with predictions from the V-QCD holographic model, which places transitions at much larger energy densities with smaller polytropic indices. The authors emphasize that their model cannot definitively establish whether stable quark stars exist in nature, since the answer depends on parameter choices, but it demonstrates that such objects can in principle emerge from holographic models, opening a computational window onto their phenomenology. Future work will explore localized D5-brane configurations, which may improve the low-density baryonic description, along with neutrino transport and rotating stars with quark cores, as multi-messenger observations continue to tighten the constraints on the densest matter in the universe.

Subject of Research: Holographic modeling of stable massive quark stars and the equation of state of dense quark matter in compact star cores

Article Title: Properties of stable massive quark stars in holography

Article References: Bitaghsir Fadafan, K., Cruz Rojas, J., & Mager, J. (2026). Properties of stable massive quark stars in holography. The European Physical Journal C, 86(9), Article 1094. https://doi.org/10.1140/epjc/s10052-026-16338-z

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16338-z

Keywords: quark stars, holographic QCD, neutron stars, D3/D7-branes, equation of state, chiral symmetry breaking, tidal deformability, compact stars, gauge-gravity duality, quark matter, TOV equations, dense nuclear matter

The motivation stems from a familiar frustration. The Lambda-CDM model — the concordance cosmology in which about 25 percent of the Universe’s energy budget is cold dark matter and roughly 70 percent is a cosmological constant — fits an impressive range of observations, from the cosmic microwave background to galaxy clustering. Yet it offers no fundamental explanation for either dark sector, and persistent tensions in the data, including disagreements over the Hubble constant and the clustering amplitude S8, hint that the framework may be incomplete. Decaying dark matter has emerged as a compelling alternative because a parent particle that slowly transforms into lighter daughters could simultaneously relieve the small-scale structure problems of cold dark matter and ease the Hubble and S8 tensions. There is, the authors argue, no compelling reason why dark matter should be perfectly stable in the first place.

The team’s framework is deliberately called semi-cosmographic, a middle path between fully model-dependent fits and completely agnostic data reconstructions. Pure cosmography expands observable quantities such as the luminosity distance as series in redshift, with the expansion coefficients constrained directly by data. But simple Taylor expansions break down above redshift one, precisely where much modern data lies. The researchers instead adopt a Padé rational approximant — a ratio of two polynomials in the variable xi, defined as the square root of one plus redshift — to describe the luminosity distance, which improves convergence at high redshift and reproduces sensible asymptotic behavior in both the low- and high-redshift limits. The three Padé parameters are treated as free quantities to be pinned down by observation rather than tied to any preselected dark energy model.

Into this flexible expansion history the authors embed a physically motivated two-body decaying dark matter sector. In their scenario, a non-relativistic parent dark matter particle decays with a constant rate Gamma, equal to the inverse of its lifetime tau, into two daughters: a massless relativistic particle, interpretable as dark radiation, and a massive daughter particle. A single dimensionless parameter epsilon controls how the parent’s rest energy is shared between the two channels. When epsilon is small, the massive daughter receives only a tiny recoil kick and behaves essentially like cold matter; when epsilon approaches one half, the daughter is born nearly relativistic. Crucially, because each decay injects daughter particles with a fixed physical momentum that then redshifts as the Universe expands, the massive daughter population is a superposition of particles produced at different epochs — early-produced daughters cool into effectively cold matter while late-produced ones remain warm. The result is a time-dependent effective equation of state for the massive daughter, which the team computes self-consistently by integrating the coupled continuity equations that govern how the parent, massless daughter and massive daughter densities evolve.

The heart of the analysis is the reconstruction of an effective residual dark energy. Once the Padé expansion history is specified and the decaying dark matter, baryon and radiation densities are known, whatever remains under the standard Friedmann equation for a spatially flat Universe must be attributed to a residual dark-energy component. Its effective equation of state, w_phi of z, is therefore not an independent model prediction but a conditional reconstruction — a diagnostic of what kind of dark energy behavior is required if the decaying dark matter picture holds. The authors are careful to stress this interpretive caveat: different parameterizations of the background expansion or different matter-sector assumptions could yield quantitatively different residual reconstructions even when the fits to data are comparably good.

For the observational constraints, the team turned to the DESI DR2 baryon acoustic oscillation data set, comprising thirteen measurements of transverse, radial, volume-averaged and anisotropic distance combinations across redshifts from about 0.295 to 2.33, along with their full covariance matrix. These geometric observables probe the smooth expansion history with exquisite precision, and the analysis recovered few-percent-level constraints on the Padé parameters and the background density parameters. However, the decay parameters — the kick strength epsilon and the lifetime tau — remained only weakly constrained by BAO data alone. The reason is subtle but important: the Padé ansatz is flexible enough that changes in the decay sector can be absorbed into adjustments of the expansion history, so many combinations of the two parameter sets reproduce nearly the same distance-redshift relation. Geometry alone cannot break that degeneracy.

This is where the 21 centimeter line enters the story. Rather than detecting individual galaxies, intensity mapping measures the collective radio emission from neutral hydrogen across large swaths of sky, tracing the large-scale distribution of matter in the post-reionization Universe between redshifts one and three. Unlike baryon acoustic oscillations, the 21 centimeter power spectrum is sensitive not only to the background expansion but also to the growth rate and scale dependence of matter clustering — exactly the quantities that decaying dark matter modifies. Decay reduces the late-time matter density and free-streaming of the kicked daughters suppresses clustering below a characteristic scale, effects that cannot be mimicked by simply reshaping the smooth background. The researchers constructed an optimistic mock power-spectrum data set at redshift 1.75 for an SKA1-MID-like interferometer with 197 dishes of fifteen-meter diameter, a system temperature of 60 Kelvin, and 4000 hours of observing time, using a simulation-based emulator calibrated on N-body simulations of two-body decaying dark matter to model the nonlinear suppression of the matter power spectrum.

The forecast results are striking. When the mock 21 centimeter likelihood is combined with the DESI DR2 BAO data, the decay parameters collapse from a broad, poorly constrained region to tight posteriors: a kick parameter of roughly 1.1 percent and a lifetime of approximately 29 gigayears at 68 percent credibility — values that place the dark matter comfortably in a regime where it behaves almost, but not exactly, like the standard cold variety. The reconstructed equation of state of the massive daughter is driven to values of order ten to the minus five over the redshift range zero to 2.3, meaning the daughter behaves effectively as cold matter for late-time structure formation. Meanwhile, the residual dark-energy equation of state, which was broadly consistent with a cosmological constant in the BAO-only analysis, becomes noticeably tighter in the joint forecast, with the band overlapping the value minus one considerably less — a weak tension that the authors interpret cautiously rather than as evidence against the cosmological constant, especially since the uncertainty grows with redshift.

The team also verified that their approach is internally consistent. Because the sound horizon used to calibrate the BAO measurements was fixed to its Planck CMB value, one might worry that decaying dark matter with a relativistic daughter could alter pre-recombination physics and invalidate that calibration. The researchers checked every posterior sample against a no-decay reference, quantifying shifts in the pre-recombination expansion rate, the matter-radiation equality redshift and the sound horizon itself. The verdict: almost no decay occurs before the drag epoch, and the maximum sound-horizon shift falls below the CMB uncertainty, confirming that the fixed calibration is self-consistent within the allowed parameter regions. Comparisons with earlier literature, including the work of Abellán and collaborators and of Fuß and Garny, show broad agreement in the low-epsilon, long-lifetime regime, though the inclusion of clustering information nudges the preferred lifetime somewhat shorter.

The authors are candid about the limitations of their forecast. The neutral hydrogen bias used in the mock analysis is calibrated on standard cold dark matter simulations and may not fully capture how decaying dark matter reshapes the halo population; foreground contamination from Galactic synchrotron emission — many orders of magnitude brighter than the cosmological signal — and the chromatic foreground wedge of the interferometer could degrade the constraints by factors of two to three, and the analysis deliberately excludes such systematics as an optimistic proof of concept. Even so, the central lesson stands with unusual clarity: geometric probes like BAO, however precise, cannot by themselves disentangle a changing expansion history from genuine particle decay. Growth-sensitive information — whether from 21 centimeter intensity mapping, weak gravitational lensing or the Lyman-alpha forest — is essential to isolate the physics of an unstable dark sector. As the Square Kilometre Array moves toward operational maturity, the crossroads of cosmography, decaying dark matter and radio cosmology may offer one of the sharpest tests yet of whether the Universe’s most abundant matter is truly as permanent as we have assumed.

Subject of Research: Semi-cosmographic reconstruction of residual dark energy in a two-body decaying dark matter cosmology using DESI DR2 BAO data and 21 cm intensity-mapping forecasts

Article Title: Padé semi-cosmographic reconstruction of residual dark energy in decaying dark matter cosmology: DESI DR2 BAO constraints and mock 21 cm forecasts

Article References: Yadav, M., Chavan, P., & Sarkar, T. G. (2026). Padé semi-cosmographic reconstruction of residual dark energy in decaying dark matter cosmology: DESI DR2 BAO constraints and mock 21 cm forecasts. The European Physical Journal C, 86(9), Article 1093. https://doi.org/10.1140/epjc/s10052-026-16335-2

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16335-2

Keywords: dark matter, dark energy, decaying dark matter, Padé approximant, cosmography, DESI DR2, baryon acoustic oscillations, 21 cm intensity mapping, SKA, equation of state, cosmological constant, large-scale structure

The problem begins with the sheer range of scales involved. In the final orbit before two neutron stars merge, the stellar cores slam into one another at a substantial fraction of the speed of light, generating a shear layer roughly a kilometer wide that becomes Kelvin-Helmholtz unstable. This instability, the same mechanism that shapes wind-blown clouds on Earth, shreds the interface between the stars into vortices that fragment into ever smaller eddies, producing a turbulent cascade that spans from kilometer scales down to about a nanometer, where viscosity finally converts kinetic energy into heat. The Reynolds number of this flow, a measure of the ratio of inertial to viscous forces, is a staggering ten to the power of sixteen. Simulating every eddy directly, the approach known as direct numerical simulation, would require computational resources that scale as the Reynolds number cubed, making it utterly impossible for the foreseeable future.

Turbulence is not confined to the moment of contact. Once the stars have merged, the remnant, whether a massive neutron star or a newly formed black hole, is typically encircled by a hot, dense accretion disk. There, the magnetorotational instability stirs the plasma, redistributing angular momentum and governing how matter spirals inward or is flung outward. The way turbulence transports angular momentum determines whether the remnant neutron star collapses promptly to a black hole or survives as a long-lived object, and it controls the mass ejection that powers the kilonova flashes and the nucleosynthesis of r-process elements. It may also amplify magnetic fields to magnetar levels, potentially launching the relativistic jets that produce short gamma-ray bursts. In short, nearly every observable signature of a neutron star merger is touched by turbulence somewhere along the way.

The mathematical machinery for handling unresolved turbulence has a long history in Newtonian fluid dynamics. The classic approach, Reynolds averaging, splits the flow into a mean component and fluctuations, yielding equations for the mean motion that contain an extra term, the Reynolds stress, which encapsulates the momentum carried by the unresolved eddies. A more practical alternative for simulations is large-eddy simulation, or LES, in which the equations are filtered over a length scale comparable to the numerical grid. The filtered equations resemble the original ones but include subgrid-scale stresses that must be modeled. The central difficulty, known as the closure problem, is that these stresses depend on information about the small scales that the simulation does not compute, so modelers must supply approximate relations, or closures, that capture the net effect of the missing physics using only the resolved quantities.

Extending this framework to general relativity introduces subtleties that have no Newtonian counterpart. Radice and Hawke review how averaging or filtering the equations of relativistic hydrodynamics produces effective stresses even when the underlying fluid is ideal, and how the nonlinear structure of the fluxes demands additional closure relations, including one for turbulent mass diffusion. More troubling still is the question of covariance: the averaging operations used in practice are tied to a particular slicing of spacetime, which breaks the four-dimensional symmetry of Einstein’s theory. Recent work has explored building the averaging procedure around a physical observer rather than a coordinate slice, showing that the coarse-grained equations then take the form of a non-ideal relativistic fluid, complete with bulk viscosity, shear stresses, and heat transport terms that arise purely from the turbulence. Even the equation of state, the relation linking pressure, density, and energy, is modified by averaging, since fluctuations in density generate corrections that behave like an additional pressure.

In practice, most published neutron star merger simulations to date have used the simplest possible strategy: implicit large-eddy simulation, which sets the subgrid stresses to zero and relies on the intrinsic numerical dissipation of shock-capturing schemes to mimic the effect of unresolved turbulence. This approach has been remarkably successful in other fields, but the review is blunt about its limitations in this context. The modified equation analysis shows that numerical dissipation can indeed act like an effective viscosity, but implicit methods require a significant fraction of the inertial range to be resolved before results converge, and no neutron star merger simulation has yet been demonstrated to be in that regime. The alternative is explicit modeling. Radice’s own general-relativistic large-eddy simulations employ a relativistic version of the Smagorinsky closure, in which the turbulent viscosity is estimated from a mixing length set by the local scale of the flow and the speed of sound. A third family of methods, gradient or approximate-deconvolution models, reconstructs the effect of the filter algebraically and has the advantage of introducing no tunable parameters beyond the filter width itself.

The payoff of these techniques is already visible in the study of magnetic field amplification. Early Newtonian simulations suggested that the Kelvin-Helmholtz instability could amplify even weak seed fields to magnetar strengths of around ten to the fifteenth gauss, but general-relativistic calculations initially failed to reproduce this, simply because their grids were too coarse. Later, extraordinarily high-resolution simulations by Kenta Kiuchi and collaborators showed that the saturated field strength kept climbing with resolution, with no sign of convergence, precisely because the magnetic back-reaction only halts the cascade at centimeter scales, far below anything a global simulation can resolve. When subgrid models were introduced, the picture changed dramatically. Gradient-model simulations by Ricard Aguilera-Miret, Carlos Palenzuela, and colleagues achieved converged results, confirming that weak fields are indeed amplified to ten to the sixteenth gauss and that the statistical properties of the resulting turbulence are remarkably insensitive to the unknown initial magnetic configuration inside the stars, a reassuring result for predictive modeling.

The same simulations revealed tantalizing evidence of an inverse cascade, in which the characteristic scale of the magnetic field grows from roughly half a kilometer immediately after merger to several kilometers a hundred milliseconds later, as turbulent resistivity rearranges field lines into larger structures. Meanwhile, measurements of the effective viscosity generated by magnetic stresses suggest it remains modest in the dense core of the remnant, implying that turbulence is unlikely to distort the post-merger gravitational-wave signal enough to compromise plans to probe the equation of state of nuclear matter with next-generation detectors such as the Einstein Telescope and Cosmic Explorer. On the other hand, turbulence and dynamo action are expected to leave a significant imprint on the long-term evolution of the remnant, its mass ejection, and its multi-messenger emission, from kilonova light curves to the engines of short gamma-ray bursts.

Much remains to be done. The review highlights open questions about whether angular momentum transport accelerates or delays the collapse of the remnant neutron star, about the topology of the amplified magnetic fields and whether tangled configurations better explain the energetics of gamma-ray bursts, and about the formidable challenge of uncertainty quantification in a parameter space already crowded with uncertain inputs. Validation is particularly thorny: unlike wind tunnels, neutron stars offer no laboratory tests, so models must be calibrated against resolved simulations whose own fidelity is uncertain, and tuned to observables, such as gravitational waves and neutrino signals, that differ from the statistical quantities conventionally used in closure validation. The authors anticipate rapid progress on three fronts: simulations that combine sophisticated microphysics, magnetohydrodynamics, and large-eddy closures; improved phenomenological subgrid models tested in local calculations; and data-driven, machine-learned closures that learn the missing physics directly from high-resolution data. As gravitational-wave astronomy enters its next generation, taming turbulence may prove the key to turning collisions of dead stars into precision measurements of matter at its densest.

Subject of Research: Turbulence modelling in general-relativistic simulations of binary neutron star mergers

Article Title: Turbulence modelling in neutron star merger simulations

Article References: Radice, D., & Hawke, I. (2024). Turbulence modelling in neutron star merger simulations. Living Reviews in Computational Astrophysics, 10(1), Article 1. https://doi.org/10.1007/s41115-023-00019-9

Image Credits: AI Generated

DOI: 10.1007/s41115-023-00019-9

Keywords: neutron star mergers, turbulence, large-eddy simulation, gravitational waves, magnetohydrodynamics, Kelvin-Helmholtz instability, magnetic field amplification, subgrid models, numerical relativity, accretion disks, r-process nucleosynthesis, equation of state

The scientific stakes have risen dramatically since 1995, when the first exoplanet around a main sequence star was detected. More than 6,000 exoplanets are now catalogued, spanning rocky Super-Earths, puffy Mini-Neptunes, and gas giants unlike anything in our own solar system. To interpret their masses, radii, gravity fields, and atmospheric spectra, scientists need equations of state, phase diagrams, transport coefficients, and optical properties for the dominant planetary ingredients across an enormous range of conditions, from dilute molecular clouds to searingly hot planetary cores. Traditional wide-range data tables such as the Sesame tables and the chemical-model equation of state of Saumon, Chabrier, and van Horn served a generation of modelers well, but the precision now delivered by missions such as Juno and Cassini exposes their limits. Ab initio simulations, anchored only in the fundamental laws of quantum mechanics and electromagnetism, have become the gold standard for this data.

The methodological core of the review is a careful walkthrough of density functional theory molecular dynamics, or DFT-MD. The many-body Schrödinger equation is untamable in its full form, so the Born-Oppenheimer approximation separates the fast electrons from the comparatively sluggish ions. The electrons are then treated quantum statistically through finite-temperature density functional theory, following the Hohenberg-Kohn theorem and the Kohn-Sham formalism introduced by Mermin’s finite-temperature extension, while the ions move classically under forces derived from the electronic structure at every molecular dynamics step. The quality of the results hinges on the exchange-correlation functional, the one genuinely approximate ingredient. Simple local density approximations and generalized gradient approximations such as PBE remain workhorses, but thermal functionals, range-separated hybrids, and van der Waals corrected functionals are increasingly important in the warm dense regime, where pressure ionization and molecular dissociation occur.

Extracting thermodynamic data from these simulations is itself a technical tour de force. The thermal equation of state combines ideal ionic motion, Coulomb interactions among the ions, and electronic pressures derived through Hellmann-Feynman forces. Caloric quantities follow from ensemble-averaged energies. To close the loop on entropies, which planetary cooling models absolutely require, the authors describe thermodynamic integration, coupling-constant integration, and the two-phase thermodynamic method based on velocity autocorrelation functions, including an advanced memory-function variant that handles superionic water. Nuclear quantum effects, often neglected because classical ions are computationally cheap, can be reinstated through harmonic post-processing corrections derived from the vibrational power spectrum. From the electronic structure, the Kubo-Greenwood formalism then yields DC electrical and thermal conductivities, viscosities, diffusion coefficients, reflectivities, and Rosseland mean opacities, all quantities that feed directly into interior and dynamo models.

The payoffs are clearest for the gas giants. Shock-compression experiments on deuterium, from gas guns to the National Ignition Facility’s lasers and Sandia’s pulsed power Z machine, have reached pressures exceeding a terapascal, and these Hugoniot data now discriminate sharply between theoretical approaches. Chemical models struggle to capture the insulator-to-metal transition of hydrogen, whereas DFT-MD with appropriate functionals reproduces the observed compression maximum. The transition itself, predicted to occur at roughly 1 to 3 megabar and with a critical point between about 1,300 and 1,900 kelvin, remains under active refinement, but theory and experiment are converging. For helium, a noble gas that resists metallization to far greater pressures, high-precision static and gas-gun data on pre-compressed samples agree with several modern functionals, although the highest-pressure laser-driven measurements still carry uncertainties too large to pick a winner.

Perhaps the most consequential phenomenon in hydrogen-helium mixtures is demixing. As the mixture cools and the hydrogen metallizes, the two elements separate into helium-rich and helium-poor phases, and the denser helium-rich droplets sink toward the core in a process poetically known as helium rain. This rain releases gravitational energy as heat, reshapes the composition profile, and slows the planet’s cooling. Ab initio calculations of the miscibility diagram, first by Lorenzen, Holst, and Redmer and later refined by Schöttler and Redmer with non-ideal mixing entropies and a van der Waals functional, indicate that both Jupiter and Saturn cross into the immiscibility regime. A single experimental study by Brygoo and colleagues in 2021 found evidence of hydrogen-helium immiscibility at Jupiter-interior conditions, and hints at demixing temperatures even higher than theory predicts, keeping this debate alive.

The interplay between these material data and spacecraft observations is rewriting our picture of the giants. Saturn’s interior path runs close to the critical point of the hydrogen insulator-to-metal transition, producing negative thermal expansion coefficients and a stably stratified layer between roughly 0.71 and 0.75 Saturn radii, and helium rain promotes further stratification between about 0.37 and 0.55 radii. Both stratified zones fall silent in the dynamo sense, helping explain Saturn’s remarkably axisymmetric magnetic field. Juno’s precision gravity moments for Jupiter, combined with first-principles equations of state, now support models featuring dilute, fuzzy cores enriched in heavy elements rather than compact rocky centers. For Saturn, ring seismology and the Cassini Grand Finale gravity data point to differential rotation and a diluted core as well, possibly topped by a helium layer born of demixing.

For the ice giants, the story turns to the so-called planetary ices of water, ammonia, and methane. First-principles equations of state for water, pioneered by French and colleagues in 2009 and continually refined, underpin the state-of-the-art AQUA equation of state used in planetary models, and recent coupling-constant integration work by Militzer suggests Uranus and Neptune may be 15 to 30 percent colder in their interiors than previously thought. Water’s crowning glory is the superionic phase, in which protons diffuse freely through a lattice of oxygen ions, delivering high ionic conductivity. Laser-driven shock experiments with nanosecond X-ray diffraction confirmed superionic ice in 2019, and diamond anvil cell studies have since mapped both face-centered cubic and body-centered cubic superionic lattices, although their stability fields still disagree substantially. Recent work even hints at close-packed mixed structures and doubly superionic compounds in which two ionic species drift through the lattice of a third.

Carbon chemistry adds a flourish worthy of the wildest planetary speculation. Methane under gigabar pressures dissociates into heavier hydrocarbons and ultimately separates into diamond and hydrogen, and machine-learning-accelerated simulations by Cheng, Hamel, and Bethkenhagen now suggest demixing into carbon occupies a far larger region of the methane phase diagram than earlier estimates. Diamond rain inside Uranus and Neptune has thus moved from thought experiment to physically grounded prediction. Hydrogen-water demixing, meanwhile, could establish stably stratified thermal boundary layers that impede heat flow and help explain why Uranus radiates so little internal heat compared with Neptune. For mixtures, the linear mixing approximation proves adequate away from phase transitions, deviating only a few percent in density, but genuine multi-component simulations increasingly powered by machine-learning potentials are essential wherever demixing, superionicity, or carbon networking dominate.

The review closes with a synthesis of where the field must go. Future missions, notably ESA’s PLATO for exoplanet characterization and NASA’s planned Uranus Orbiter and Probe, will shower planetary modelers with fresh gravity, atmospheric, and seismological data that only first-principles material properties can decode. On the experimental side, next-generation X-ray free-electron lasers and new dynamic diamond anvil cell techniques will probe ever more extreme states, while machine-learning force fields trained on DFT-MD extend simulations to sizes and timescales once unthinkable. Entropy determination for complex multi-component mixtures remains the central technical hurdle, because it controls the demixing diagrams and phase boundaries on which every interior, evolution, and dynamo model ultimately rests. What is no longer in doubt is that simulating matter at planetary conditions from first principles has transformed planetary science from informed guesswork into a predictive, quantitative discipline, one that can peer through thousands of kilometers of opaque fluid and read the hidden architecture of the giant worlds that share our cosmic neighborhood.

Subject of Research: Ab initio simulation of warm dense matter to model the interiors, thermal evolution, and magnetic fields of giant planets

Article Title: Ab initio description of warm dense matter applied to the physics of giant planets

Article References: Bethkenhagen, M., Preising, M., & Redmer, R. (2026). Ab initio description of warm dense matter applied to the physics of giant planets. Living Reviews in Computational Astrophysics, 12(1), Article 4. https://doi.org/10.1007/s41115-026-00030-w

Image Credits: AI Generated

DOI: 10.1007/s41115-026-00030-w

Keywords: warm dense matter, ab initio simulations, DFT-MD, giant planets, planetary interiors, helium rain, superionic water, hydrogen metallization, ice giants, equation of state, magnetic dynamo, diamond rain