The scale separation in space is staggering. Reynolds numbers, which measure the ratio of inertial to viscous effects, routinely reach 10 billion or larger in astrophysical flows, meaning turbulence is essentially unavoidable. Unlike water or air, astrophysical plasmas are almost perfectly conducting, so their dynamics are governed by the MHD equations, coupling currents, magnetic fields and the Lorentz force. In our Galaxy, magnetic fields of around 5 microgauss reach rough equipartition with turbulent kinetic energy; in galaxy clusters, fields of 1 to 3 microgauss sit near a twentieth of equipartition. These fields are not decorative: they accelerate particles, fill the Universe with non-thermal radiation, and even enable the accretion of matter onto black holes through the magnetorotational instability, a process estimated to be the most potent energy source in the cosmos, exceeding thermonuclear burning in stars.

The mathematical backbone of the subject remains the Kolmogorov cascade picture. Energy injected at large scales passes, without dissipation, through an inertial range of scales until viscosity finally wins at the Kolmogorov scale. Dimensional analysis then yields the famous spectrum in which energy content scales as the wavenumber to the minus five-thirds power, with a dimensionless Kolmogorov constant that experiments and simulations place near 1.6. Beryesnyak emphasizes a rigorous tool called scaling convergence: because the underlying equations contain no preferred scale, small-scale statistics from simulations of different resolutions should collapse onto a single universal curve when properly normalized, a technique that drives statistical error virtually to zero and has even resolved the tiny intermittency correction of about 0.04 to the spectral slope.

MHD turbulence, however, refuses to behave like its hydrodynamic cousin. A large-scale magnetic field cannot be transformed away, and it remains dynamically important on every scale. Linearizing the MHD equations reveals four wave modes, of which the transverse Alfvén mode dominates the nonlinear cascade. Early theorists Iroshnikov and Kraichnan imagined weak, wave-like interactions, but later work showed that turbulence becomes stronger, not weaker, as it cascades. Energy transfer proceeds preferentially perpendicular to the field, producing extreme anisotropy. Goldreich and Sridhar proposed that this anisotropy is capped by critical balance, where the cascade time matches the wave period, yielding a perpendicular spectrum of minus five-thirds and the relation that the parallel wavenumber scales as the perpendicular wavenumber to the two-thirds power.

High-resolution direct numerical simulations have now put these ideas to demanding tests. Using the scaling convergence method on simulations up to 4096 cubed grid points, Beresnyak found that the perpendicular spectrum converges best near a slope of minus 1.7, consistent with the Kolmogorov picture rather than competing minus three-halves models. The Alfvénic Kolmogorov constant was measured as 3.3 with a total value near 4.2 when the passively advected slow mode is included. Remarkably, the residual energy, the difference between magnetic and kinetic energy, turns out to be a constant fraction, about 15 percent, of the total energy throughout the inertial range, with a corresponding Alfvén ratio of roughly 0.74, resolving earlier conceptual difficulties with theories that predicted scale-dependent behavior.

An elegant theoretical result connects the parallel spectrum to Lagrangian statistics. Because oppositely directed Alfvén wave packets propagate along magnetic field lines at a fixed speed, measuring fluctuations along the field is mathematically equivalent to following a fluid element in time. This argument yields a parallel spectrum proportional to the wavenumber to the minus second power, scaled by the inverse of the Alfvén speed, without ever invoking critical balance. Numerical measurements along the local magnetic field overwhelmingly confirm this minus-two law, matching observations from the solar wind, where spacecraft such as Helios 2 have recorded clean power-law spectra over decades of frequency.

Imbalanced turbulence, where waves traveling one direction dominate, presents a harder puzzle, and it is the norm in the solar wind and near astrophysical jets. Because critical balance cannot hold simultaneously for counter-propagating waves of unequal amplitude, several competing models were proposed. Simulations with systematically varied imbalance show that the Lithwick-Goldreich-Sridhar model captures the spectra at small imbalances, while the Beresnyak-Lazarian model, which relaxes locality for the dominant component, best matches the energy ratios and the diverging anisotropies of the two populations at strong imbalance.

Perhaps the most consequential result concerns the small-scale dynamo, the process by which turbulence amplifies weak magnetic fields. Once the kinematic, exponential phase ends, magnetic energy grows linearly in time as turbulence converts a fixed fraction of cascade power into magnetism. The measured efficiency constant is small, about 0.05, but the implications are enormous: in galaxy clusters, this keeps the ratio of magnetic to thermal energy constant at roughly 40 over the past 10 billion years, in agreement with Faraday rotation observations. The review also issues a caution to simulators. Because numerical Reynolds numbers are vastly smaller than astrophysical ones, starting a simulated young object such as a collapsing cloud with a vanishing field can artificially delay magnetization, and implicit large-eddy codes with zero initial field produce no field at all, in gross contradiction with nature, where the dynamo always jump-starts itself.

The frontier now extends to supersonic turbulence in molecular clouds, where Mach numbers near 10 produce densities varying by orders of magnitude and log-normal probability distributions sculpted by slow-mode shocks and sheared by Alfvénic motions, with direct consequences for star formation theory. It extends, too, to magnetic reconnection, where current sheets tear and spawn their own strong, critically balanced turbulence, producing reconnection rates of about 1.5 percent of the Alfvén speed that are independent of resistivity. From solar flares to black hole jets to the magnetization of the cosmic web, MHD turbulence has become the connective tissue of modern astrophysics, and the convergence of theory, simulation and spacecraft measurement suggests the field is finally converging on universal answers.

Subject of Research: Theory and numerical simulation of magnetohydrodynamic turbulence and its astrophysical applications

Article Title: MHD turbulence

Article References: Beresnyak, A. (2019). MHD turbulence. Living Reviews in Computational Astrophysics, 5(1), Article 2. https://doi.org/10.1007/s41115-019-0005-8

Image Credits: AI Generated

DOI: 10.1007/s41115-019-0005-8

Keywords: MHD turbulence, magnetohydrodynamics, astrophysical plasmas, Kolmogorov cascade, Alfvén waves, turbulent dynamo, solar wind, galaxy clusters, magnetic reconnection, interstellar medium, numerical simulations, star formation

The core obstacle is captured by a single number, the Reynolds number, which measures the ratio of non-linear inertial forces to viscous damping in a flow. When this number exceeds roughly a few thousand, flows become turbulent. Astrophysical flows routinely reach values that are astronomically larger; the turbulent convection zone of the Sun, for instance, is estimated to have a Reynolds number of order ten to the fourteenth. Direct numerical simulations that resolve every scale down to the microscopic dissipation length would require computing resources that scale with the Reynolds number cubed. For solar convection alone, the review notes, a single dynamical time would demand roughly ten to the forty-second floating point operations, far beyond even exascale machines.

Large eddy simulations sidestep this impossibility through a mathematical procedure called spatial filtering. A low-pass filter, defined by a convolution with a kernel such as a box, Gaussian, or sharp cut-off filter, separates the flow into smoothed, large-scale variables and fluctuating components on scales smaller than the filter length, which is usually identified with the grid spacing. Because the equations of fluid dynamics are non-linear, filtering them generates new terms: a subgrid-scale turbulence stress tensor that describes the momentum exchange between resolved and unresolved eddies, a turbulent pressure proportional to the subgrid-scale kinetic energy, and additional fluxes for energy and chemical species. The trace of the stress tensor defines the subgrid-scale turbulence energy, the kinetic energy hidden in eddies too small for the grid to see. Modelling these terms is known as the closure problem, since the filtered equations form an infinite hierarchy of moments that must be truncated.

Several families of closures dominate astrophysical practice. The simplest is the Smagorinsky model, which assumes a local balance between turbulence production and dissipation and yields an eddy viscosity proportional to the square of the grid scale times the local rate of strain. More sophisticated is the one-equation model, which solves a transport equation for the subgrid-scale turbulence energy itself, including production by shear, dissipation, turbulent diffusion, and pressure-dilatation effects. A third approach, structural modelling, reconstructs the unresolved stresses directly from gradients of the resolved velocity and magnetic fields, without tuneable coefficients or extra equations. For magnetohydrodynamics, the structural model also provides a closure for the subgrid-scale electromotive force, the term that governs how unresolved turbulent motions amplify magnetic fields through dynamo action.

A crucial subtlety is that most astrophysical simulation codes already dissipate energy through their numerical schemes. Finite-volume and finite-difference truncation errors behave like diffusion terms, producing an effective numerical viscosity that mimics turbulent viscosity on the grid scale. This observation underlies implicit large eddy simulation, or ILES, in which no explicit model is used and the numerics themselves play the role of the closure. The review explains why ILES often suffices: because inertial-range scaling is largely independent of the dissipation mechanism, statistics such as energy spectra can be reproduced as long as the dynamical range is adequate. However, the so-called bottleneck effect distorts spectra near the grid scale, and explicit subgrid-scale models add only a marginal effect for strongly diffusive second-order schemes. The picture changes for higher-order methods and mesh-free Lagrangian codes, which lack intrinsic diffusion and benefit far more from explicit modelling.

Validating these models without experiments or direct numerical simulations is a distinctive challenge of astrophysics. The standard tools are a priori tests, in which turbulence data from simulations are explicitly filtered and the correlations between modelled and true subgrid-scale terms are measured, and a posteriori comparisons of turbulence statistics between LES and ILES. Such analyses have calibrated the coefficients of the eddy-viscosity and gradient-diffusion closures, revealing, for example, that the turbulent Prandtl number is around ten rather than unity, and that purely linear eddy-diffusivity closures fail badly for magnetic stresses, matching the sign of the energy cascade only about half the time. Non-linear structural closures, by contrast, achieve correlations close to unity for the kinetic and magnetic stress tensors and the electromotive force, especially when a compressibility correction is included.

The payoff of this machinery appears across a remarkable range of cosmic phenomena. In thermonuclear supernovae, where a white dwarf detonates in a deflagration that no grid can resolve, the subgrid-scale turbulence energy sets the effective flame propagation speed, replacing the microscopic laminar burning velocity with a turbulent one. In galaxy simulations, the subgrid-scale turbulence energy can be fed by supernova feedback as an additional source term, and it allows a locally variable, turbulence-regulated star formation efficiency to be computed from the turbulent Mach number and virial parameter, reproducing observed star formation laws without assuming them. In cosmological simulations, filtered equations formulated in co-moving coordinates show that unresolved turbulent pressure can contribute a non-negligible fraction of the support against gravity in galaxy clusters, while shear-improved models cleanly separate genuine turbulence from gravity-driven bulk flows.

Perhaps the most striking recent success concerns magnetic field amplification in binary neutron star mergers. Kelvin-Helmholtz instabilities in the shear layer between the merging cores trigger a small-scale dynamo, and the growth rate of the magnetic field normally increases with resolution because the smallest resolved eddies dominate the amplification. By combining a relativistic generalization of the structural subgrid-scale model with high-order schemes, researchers achieved, for the first time, a numerically converged magnetic field amplification in merger simulations: the large eddy simulation reached the saturated field strength at half the resolution required without the model, a substantial saving in computational cost. The resulting magnetic energy spectrum even displays the expected Kazantsev scaling, and the post-merger field appears universal, insensitive to the initial magnetic configuration of the two stars.

Subgrid-scale modelling also matters for the chemical evolution of galaxies. In mesh-free codes without numerical diffusion, metals expelled by supernovae would never mix with their surroundings without an explicit turbulent diffusivity, and studies show that such mixing significantly reshapes the metallicity distributions of the circumgalactic and warm-hot intergalactic media. Looking forward, the review argues that the main utility of subgrid-scale models lies precisely in this treatment of complex sub-resolution physics, from turbulent burning and star formation to dynamo action and metal transport. While statistical properties of well-resolved turbulence remain largely insensitive to the choice of model, the coupling between resolved and unresolved scales is real and physically motivated closures capture it better than numerical truncation errors alone. As exascale computing opens new regimes, large eddy simulations are poised to remain an indispensable bridge between the eddies we can compute and the turbulent universe we observe.

Subject of Research: The methodology of large eddy simulations and subgrid-scale modelling for turbulent astrophysical flows

Article Title: Large eddy simulations in astrophysics

Article References: Schmidt-Brückner, W. (2025). Large eddy simulations in astrophysics. Living Reviews in Computational Astrophysics, 11(1), Article 2. https://doi.org/10.1007/s41115-025-00023-1

Image Credits: AI Generated

DOI: 10.1007/s41115-025-00023-1

Keywords: large eddy simulation, subgrid-scale modelling, astrophysical turbulence, magnetohydrodynamics, computational astrophysics, turbulent dynamo, type Ia supernovae, star formation, galaxy clusters, neutron star mergers, turbulent mixing, cosmological simulations

The key observable is the dispersion measure, a quantity derived from the tiny frequency-dependent delay imposed on a radio burst as it traverses ionized gas. Because free electrons slow lower-frequency radio waves slightly more than higher-frequency ones, each burst arrives smeared across a characteristic sweep, and the total delay encodes the column density of electrons along the line of sight. After subtracting the contribution of our own Milky Way and the burst’s host galaxy, the remaining extragalactic dispersion measure is a direct tally of the ionized baryons lying between the source and the observer — gas threaded through the intergalactic medium, the circumgalactic halos of intervening galaxies, and the hot atmospheres of groups and clusters. This is precisely the baryonic material whose distribution has been sculpted by supernova explosions, jets from supermassive black holes, and other feedback processes that blow gas out of galaxies and redistribute it across megaparsec scales.

That redistribution matters enormously for cosmology. Surveys of weak gravitational lensing — the subtle distortion of galaxy images by intervening matter — measure the clustering of all matter, dark and luminous alike, on these same scales. But the theoretical predictions that lensing measurements are compared against must account for how feedback expels gas from dark matter halos, smoothing the matter distribution and suppressing the amplitude of the matter power spectrum on small scales. If that suppression is mis-modeled, inferred cosmological parameters such as the clumpiness of matter, the sum of neutrino masses, and the properties of dark energy can all be biased. Until now, the primary probes of this baryonic physics have been X-ray observations of hot gas and the Sunyaev–Zel’dovich effects, in which hot electrons leave signatures in the cosmic microwave background. Both approaches, however, have faced tensions and systematics, particularly in galaxy groups and lower-mass clusters.

The new analysis demonstrates that fast radio bursts are already competitive with the legacy measurements from the Atacama Cosmology Telescope and the eROSITA X-ray telescope. Using a halo-model inference framework calibrated on hydrodynamical simulations, the team converted the observed dispersion measures of 114 bursts into constraints on the gas content of dark matter halos more massive than about 10¹³ solar masses, and on the degree to which feedback suppresses the matter power spectrum. The framework builds on the well-established Macquart relation, the observed linear trend between a burst’s dispersion measure and its redshift, which was previously used to account for the Universe’s so-called missing baryons. Here the authors went a step further, treating the scatter and fluctuations around that relation as a signal of how baryons are distributed within and around halos.

Technically, the inference proceeds by modeling gas profiles within halos using flexible analytic prescriptions whose parameters — including a characteristic feedback mass scale and the radial extent of gas ejected from halos — are constrained through a Markov Chain Monte Carlo fit to the burst sample. A critical systematic is the dispersion measure contributed by each burst’s own host galaxy, which the team treated as a nuisance parameter, allowing its mean and scatter to float in the fit. Robustness tests splitting the sample into lower- and higher-redshift subsets, and allowing the host contribution to evolve with the cosmic star formation history, showed that the inferred feedback constraints are stable across all configurations, indicating that uncertainties about host galaxies do not drive the result. Lower-redshift bursts in the sample anchor the dispersion measure–redshift relation and pin down the host distribution, while higher-redshift bursts sharpen the sensitivity to feedback physics.

The measurements reveal clear signatures of efficient gas expulsion from massive halos, quantified as a suppression of the matter power spectrum at wavenumbers between roughly 0.1 and 3 h Mpc⁻¹ — the regime where weak lensing surveys are most sensitive to baryonic effects and where existing X-ray and thermal Sunyaev–Zel’dovich measurements have disagreed. By constraining the gas mass fraction within group- and cluster-scale halos, the burst data provide an independent check on the scaling relations that underpin cluster cosmology, and they help arbitrate the tensions that have emerged between different baryon surveys. The result establishes fast radio bursts as a genuinely new probe of feedback-regulated structure formation, complementary to lensing, X-ray and microwave-background methods because it is sensitive to all ionized gas regardless of its temperature.

The implications extend beyond astrophysics into fundamental physics. Because baryonic feedback and cosmology are partially degenerate in lensing measurements, independently constraining feedback breaks those degeneracies and sharpens cosmological inference. Forecast analyses accompanying the study show that combining fast radio burst dispersion statistics with weak lensing from a survey like the Vera C. Rubin Observatory’s Legacy Survey of Space and Time — through dispersion-galaxy and dispersion-shear cross-correlations — would tighten constraints on the feedback mass scale dramatically and propagate into substantially improved limits on the sum of neutrino masses and on dynamical dark energy parameters. In an era when percent-level control of baryonic physics is a prerequisite for precision cosmology, a probe that measures the baryons directly and independently is a valuable asset.

The field is poised for rapid growth. Instruments now coming online, including the Deep Synoptic Array, the Canadian Hydrogen Observatory and Radio-transient Detector, and CHIME/FRB with its outrigger stations, are expected to deliver localized bursts at rates hundreds of times higher than current samples, extending the redshift baseline and improving the statistical power of dispersion-based probes. The authors’ forecasts indicate that within the next decade, fast radio bursts could deliver leading constraints on baryonic physics, rivalling or exceeding the multi-probe combinations that currently define the field. As sample sizes grow, the same data will also refine our understanding of the bursts’ own progenitors and host galaxy populations, closing the loop between astrophysics and cosmology in a single dataset.

For now, the result stands as a striking demonstration of scientific serendipity: signals once dismissed as inexplicable flashes of radio noise have become a census of the invisible matter that threads the Universe, revealing not only where cosmic matter resides, but how the explosive feedback of galaxies has smoothed it away. In revealing the suppressed clustering of matter through the electrons it left behind, 114 fleeting radio flashes have delivered a measurement that decades of X-ray and microwave observations have struggled to pin down — and they promise much more to come.

Subject of Research: Using fast radio burst dispersion measures to quantify baryonic feedback and the suppression of the matter power spectrum in structure formation.

Article Title: Signatures of suppressed matter clustering revealed by fast radio bursts

Article References: Sharma, K., Krause, E., Ravi, V., Connor, L., Anbajagane, D., & Rajendra Singh, P. (2026). Signatures of suppressed matter clustering revealed by fast radio bursts. Nature Astronomy. https://doi.org/10.1038/s41550-026-02957-9

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02957-9

Keywords: fast radio bursts, cosmology, baryonic feedback, matter power spectrum, weak lensing, galaxy clusters, dispersion measure, dark matter halos, missing baryons, Sunyaev–Zel'dovich effect, large-scale structure, Nature Astronomy