PAHs are the workhorses of the interstellar medium in galaxies like our own. These carbonaceous nanoparticles, consisting of linked aromatic rings, absorb ultraviolet photons from young massive stars and re-emit the energy in a set of characteristic mid-infrared bands at 3.3, 6.3, 7.7, 8.6 and 11.3 micrometres, among others. In metal-rich star-forming galaxies, these features account for 5 to 20 percent of the total infrared luminosity, making PAHs a widely used tracer of both star formation rates and molecular gas. They also play a physical role: in neutral gas, PAHs heat the medium through the photoelectric effect and regulate its ionization balance by recombining ions and electrons.

For decades, however, observations with Spitzer and other infrared telescopes have shown that PAH emission collapses as metallicity falls. The fraction of the dust budget locked up in PAHs drops sharply around a metallicity of roughly 30 percent of the solar value, and below about 10 percent solar, PAHs had remained stubbornly undetected even by JWST. Two competing explanations have been debated: either the fragile molecules are being destroyed more rapidly by harsh radiation fields, shocks and electron collisions in the poorly shielded, dust-poor gas, or their formation is being suppressed in the first place because of carbon scarcity and a lack of dense gas where grains can grow.

The team, led by Elizabeth Tarantino of the Space Telescope Science Institute, targeted Sextans A because it offered a rare combination of properties. The galaxy sits 1.4 megaparsecs away, contains about 6.2 times ten to the seventh solar masses of atomic gas, and is still forming stars at a modest rate of roughly 0.012 solar masses per year. Crucially, Herschel had previously detected cold dust there, signalling that a JWST search for PAHs had a realistic chance of success. The observations used NIRCam and MIRI imaging in filters centred on the 3.3, 7.7 and 11.3 micrometre PAH features, flanked by continuum filters that allowed the team to subtract the overwhelming glow of stars and hot dust.

The detection was unambiguous. In the brightest star-forming region of the galaxy, the spectral energy distribution of a compact clump showed clear excess flux in all three PAH filters relative to the bracketing continuum bands, confirming that the signal came from aromatic molecules rather than a spurious background source. What surprised the team most was the morphology. Instead of the extended, diffuse PAH emission seen in metal-rich galaxies, where PAH-emitting regions can span 500 parsecs or more, the PAHs in Sextans A are confined to tiny clumps just 3 to 10 parsecs across, with average radii of about 2 parsecs. Many are unresolved even at JWST’s exquisite resolution.

That compactness explains why previous instruments missed the emission entirely. Spitzer, with its coarser resolution, would have smeared these point-like clumps into invisibility against the bright continuum. The team quantified the effect by comparing the PAH fraction measured over the full far-infrared beam with the fraction measured on the clumps themselves, finding that the clumps are dramatically richer in PAHs relative to small dust grains than their surroundings. The measured ratio of PAH luminosity to total infrared luminosity in Sextans A, about 0.031 percent, extends the well-known metallicity trend to its lowest point yet, but the authors caution that beam dilution means the true PAH fraction within the clumps is substantially higher than the galaxy-wide average suggests.

The band ratios told an equally important story. The 3.3 to 11.3 micrometre ratio, which traces grain size, indicated that the PAH population in Sextans A consists of small and neutral grains, consistent with model predictions for inhibited grain growth. More strikingly, the 3.3 to 7.7 micrometre ratios were so elevated that neither the Draine-Li model grids nor the alternative THEMIS dust models could reproduce them, implying the PAHs are even more neutral than the models’ lowest-ionization cases. Because the 7.7 micrometre feature arises primarily from charged PAHs, its weakness suggests the molecules live in dense, cold, well-shielded gas where photoionization is suppressed.

Perhaps the most counterintuitive result concerns the radiation field. If enhanced photodestruction were the main cause of PAH poverty at low metallicity, the smallest and most fragile grains should be scarcest where ultraviolet radiation is strongest, driving the 3.3 micrometre feature down relative to longer wavelengths. Instead, the team found the opposite: the 3.3 micrometre ratios correlate positively with H-alpha and ultraviolet flux, likely because harder radiation fields preferentially boost the 3.3 micrometre feature itself. The absence of evidence for radiation processing undermines the destruction-only scenario and instead favours inhibited grain growth as the dominant explanation for the PAH deficit in Sextans A.

The authors argue that the compact clumps are active sites of in situ PAH formation within a dense, shielded phase of the interstellar medium. Because the dust-to-gas ratio scales with metallicity, a much larger gas column is needed in Sextans A to achieve the same visual extinction as in the Milky Way, shrinking the physical volume where shielding is effective and confining PAH growth to molecular cores. Moreover, PAH growth timescales scale roughly with density and inversely with metallicity, meaning growth in Sextans A proceeds at least thirty times more slowly than at solar abundance, naturally producing a population dominated by small grains.

The implications reach far beyond one nearby dwarf. Sextans A demonstrates that PAHs can form and survive in environments resembling those of galaxies at high redshift, where JWST is now routinely detecting aromatic emission at metallicities well below solar. If the PAH life cycle in such systems is governed by the balance between grain growth in dense clumps and shielding from hard radiation fields, then interpreting PAH-based star formation and molecular gas tracers in the early Universe will require accounting for where, and in what physical conditions, these remarkable molecules manage to build themselves.

Subject of Research: Detection and characterization of polycyclic aromatic hydrocarbon dust grains in the extremely metal-poor dwarf galaxy Sextans A using JWST imaging

Article Title: Growth of aromatic hydrocarbon dust particles in the extremely metal-poor galaxy Sextans A

Article References: Tarantino, E. J., Roman-Duval, J., Sandstrom, K. M., Whitcomb, C. M., Smith, J.-D. T., Draine, B. T., Boyer, M. L., Chastenet, J., Chown, R., Clark, C. J. R., Elyajouri, M., Gordon, K. D., Hensley, B. S., Lai, T. S.-Y., Lindberg, C. W., McQuinn, K. B. W., Newman, M. J. B., Telford, O. G., Putte, D. V. D., & Williams, B. F. (2026). Growth of aromatic hydrocarbon dust particles in the extremely metal-poor galaxy Sextans A. Nature Astronomy. https://doi.org/10.1038/s41550-026-02969-5

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02969-5

Keywords: polycyclic aromatic hydrocarbons, Sextans A, JWST, metal-poor galaxies, interstellar dust, dwarf galaxies, mid-infrared emission, interstellar medium, grain growth, Local Group, dust shielding, Nature Astronomy

The research, published in the Monthly Notices of the Royal Astronomical Society, was led by Anmol Aggarwal, a PhD student at UCL’s Mullard Space Science Laboratory, together with Dr Ralph Schoenrich, also based at the same laboratory. Their central insight is deceptively simple: when a star explodes, its material is not flung evenly into space. Different elements are ejected in different directions, with oxygen potentially streaming one way and sulphur another. New stars that condense from this debris may therefore inherit a lopsided sample of the supernova’s ingredients rather than its overall chemical blend.

Previous models, the researchers argue, assumed that the gas ejected by a supernova was thoroughly mixed before it formed new stars. Under that assumption, the peculiar chemical signatures observed in certain stars could not be produced by an ordinary supernova, and something far more energetic was required. Aggarwal and Schoenrich built mathematical models that instead accounted for the uneven, asymmetric way material is expelled, fitting how much material from each layer or region of an exploding star would be needed to match the composition of the observed stars.

Across the two studies, the team examined four stars that had previously been judged to have formed from material expelled by a hypernova. In the first study, they found that for three of these stars, formation from an ordinary supernova was statistically favoured over a hypernova. The unusual abundance patterns that had seemed to demand an ultra-powerful explosion could instead arise naturally if the stars formed from incomplete, patchy samples of supernova debris. As Aggarwal put it, stars with a very unusual mix of ingredients were most likely formed from ordinary supernovae, and all the evidence for hypernovae suddenly disappears.

To understand why this matters, it helps to consider how a supernova actually works. When a massive star exhausts its nuclear fuel, its core collapses under gravity, triggering an immense explosion that blasts the star’s outer layers into space. Just before this happens, the star is structured like an onion, with the heaviest elements concentrated at the centre and progressively lighter elements in the outer layers. A hypernova is theorised to occur when the same collapse happens in an even more massive, rapidly spinning star, producing an explosion of vastly greater energy.

These explosions seed nearby star-forming clouds with heavy elements, and the process gives astronomers a powerful investigative tool. By analysing the composition of a star, researchers can infer its family history, effectively reading the chemical record of the stellar explosions that contributed to the gas from which it formed. But this forensic method depends on assumptions about how well ejected material mixes with the surrounding interstellar gas before new stars are born. Schoenrich emphasised that, despite decades of work in this area, scientists still do not know how thoroughly supernova material gets mixed before forming new stars. The new research suggests that some mixing occurs, but that it is incomplete, a conclusion with consequences far beyond the handful of stars studied.

The second study focused on a single red giant in the Milky Way’s halo, the vast cloud of stars surrounding the galaxy’s disc. This star has an exceptionally peculiar chemical make-up, including an abundance of certain heavy elements such as silver and uranium. Elements this heavy can only be produced by neutron star mergers, the collisions of ultra-dense stellar remnants, or, in theory, by the hypernova of a highly magnetised massive star capable of forging both these heavy elements and the standard supernova elements no heavier than iron and nickel.

Earlier work had ruled out a neutron star merger as the origin of this star’s material, on the grounds that the star is extremely metal-poor, meaning it contains a very low proportion of heavier elements created in stellar explosions. Such metal-poor stars are understood to have formed from gas barely enriched by earlier generations of stars. The previous analysis relied on a model of the Milky Way in which the galaxy’s gas was enriched with heavy elements very quickly, implying that by the time the star was born, insufficient time had passed for neutron stars to form, spiral together, and collide. That reasoning left the hypernova as the only remaining explanation, and the star became one of the flagship cases for these ultra-powerful explosions.

Aggarwal and Schoenrich presented a different picture. They noted the star’s exceptional speed and the fact that it orbits the Milky Way in the opposite direction to most other stars, characteristics typical of stars captured from smaller galaxies. They argued that the star most likely originated in a small dwarf galaxy, where metal-poor stars can form much later, many millions of years after the Big Bang, before being pulled into the Milky Way. In tiny galaxies, star formation proceeds more slowly, and when stars do explode, most of the ejected material is blasted out into intergalactic space, keeping the galaxy’s gas metal-poor for far longer than in a large galaxy like our own. Using their mathematical model of asymmetric explosions, the researchers concluded that this star most likely formed from a neutron star merger combined with a single ordinary supernova, removing the need to invoke a hypernova at all.

The implications of the two studies extend well into the future of galactic chemistry modelling. If supernova ejecta are not well mixed, then models of how chemicals evolve in galaxies, and of how the interstellar medium behaves, need to incorporate this incompleteness explicitly. Chemical evolution models that treat each supernova as delivering a uniform, averaged blend of elements may misinterpret unusual abundance patterns as signs of exotic events, when they are in fact the natural consequence of patchy, directional ejecta. The UCL team’s work suggests that astronomers should be cautious before attributing strange stellar chemistries to the most extreme explosions imaginable, and that the humble supernova, viewed through a more realistic lens, may be capable of far more chemical variety than previously assumed. Whether genuine hypernovae exist at all remains an open question, but according to this new analysis, the evidence once marshalled in their favour has suddenly gone.

Subject of Research: The chemical signatures of hypernovae and supernovae in metal-poor stars

Article Title: Ultra-powerful star explosions might not have occurred after all

Article References: Ultra-powerful star explosions might not have occurred after all. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: hypernova, supernova, stellar nucleosynthesis, metal-poor stars, neutron star merger, Milky Way halo, dwarf galaxy, chemical evolution, interstellar medium, UCL, Monthly Notices of the Royal Astronomical Society, stellar explosions

The diffuse Galactic emission is the fog of high-energy astrophysics. It arises when cosmic rays, protons and electrons accelerated by sources such as supernova remnants, collide with interstellar gas, dust grains, magnetic fields, and the sea of starlight that permeates the Galaxy. Protons smashing into gas produce neutral pions that decay into gamma rays, while electrons scatter low-energy photons up to gamma-ray energies through inverse Compton processes and radiate synchrotron light in magnetic fields. Disentangling these components is notoriously difficult because their contributions overlap and depend on poorly constrained distributions of gas, radiation, and cosmic rays. Traditional models, such as the widely used GALPROP code, solve the cosmic-ray transport equation with detailed parameterizations, but they require extensive tuning and still struggle in the structurally complex inner Galaxy.

The research team, led by Xi Liu and colleagues at Sun Yat-sen University, took a deliberately different route. Rather than assuming a physical model from first principles, they trained two standard machine learning algorithms, Random Forest and K-Nearest Neighbors regression, to learn a direct mapping between Planck’s nine all-sky foreground maps, spanning 30 to 857 gigahertz, and the gamma-ray intensity predicted by the Fermi-LAT interstellar emission model across 28 energy bins from 50 megaelectronvolts to 814 gigaelectronvolts. Each Planck band traces a different physical regime: the lowest frequencies are dominated by synchrotron radiation from cosmic-ray electrons, intermediate bands mix free-free and spinning-dust emission, and the highest frequencies are dominated by thermal dust, a reliable tracer of interstellar gas.

The results are striking. In the 0.1 to 10 gigaelectronvolt range, where the diffuse gamma-ray sky is brightest, the models achieve coefficients of determination above 0.90, peaking at 0.96 at 687 megaelectronvolts. The predicted full-sky maps reproduce the major morphological features of the gamma-ray Galaxy, and the extracted spectral energy distribution of the inner Galactic plane matches the reference model closely. The team verified that the learned relationship is not a local accident: when a model trained on one Galactic hemisphere was tested on the opposite side, performance varied by no more than 6 percent regardless of how large a training region was used, indicating that the microwave-to-gamma-ray connection reflects a robust, large-scale property of the Milky Way rather than a statistical fluke of nearby structures.

Perhaps the most physically revealing result comes from asking which Planck frequencies matter most. When the high-frequency bands, dominated by thermal dust emission, were used alone, they achieved predictive power nearly identical to that of all nine bands in the 0.1 to 10 gigaelectronvolt range, while the low-frequency synchrotron channels performed markedly worse. This is exactly what one would expect if the diffuse gamma-ray emission at these energies is hadronically dominated: gamma-ray intensity is proportional to the product of gas column density and cosmic-ray density, and dust is an excellent proxy for the gas that serves as the target material. Above 10 gigaelectronvolts, however, the predictive power of the low-frequency bands rises, consistent with a growing leptonic contribution in which the synchrotron-emitting electron population also produces gamma rays through inverse Compton scattering. The analysis thus provides an independent, data-driven confirmation of the standard picture of Galactic gamma-ray production.

The residual maps, showing where the learned mapping breaks down, read like a tour of the Galaxy’s most enigmatic structures. Positive residuals appear toward the Magellanic Clouds and Centaurus A, where abundant gas and dust in nearby systems tempt the model to predict hadronic gamma-ray emission that the reference maps treat separately. Warm ionized gas complexes such as Barnard’s Loop in Orion and the Gum Nebula also show over-predictions, because free-free emission is subdominant across the Planck bands. On the other side, negative residuals trace Loop I and the North Polar Spur, hinting at inverse-Compton-dominated radiation fields or enhanced electron populations that submillimeter dust tracers cannot capture. Around the Galactic Center, the residuals mirror the shape of the Fermi Bubbles, the giant gamma-ray lobes whose faint, hard-spectrum microwave counterpart is too subtle for the model to disentangle without explicit priors.

The team also uncovered a curious hemispheric asymmetry: at high Galactic latitudes, the model systematically under-predicts the north while over-predicting the south by 10 to 30 percent on average. Because the discrepancy is spatially coherent rather than random, it may point to genuine physical differences between the two Galactic hemispheres, such as variations in cosmic-ray density or radiation fields. A segmented analysis of the inner Galaxy further showed that the asymmetry is concentrated within 30 degrees of the Galactic Center, likely linked to the Galactic bar and spiral-arm tangents, while the outer disk behaves with remarkable uniformity.

In a direct head-to-head comparison at roughly 4.3 gigaelectronvolts, the machine learning approach outperformed a GALPROP simulation tuned to the latest AMS-02 cosmic-ray data in the inner disk and Galactic Center region, achieving a coefficient of determination of 0.9465 and a mean absolute relative error of 14.7 percent, compared with 0.8624 and 22.3 percent for the physical model. GALPROP, by contrast, held a slight edge in the more homogeneous outer disk and in parts of the halo, confirming its reliability under simpler conditions. The authors emphasize that the two approaches are complementary: the data-driven model absorbs non-linear multi-frequency correlations that static gas distributions miss, while the physical simulation remains indispensable for interpretation and extrapolation.

An appendix to the study adds a compelling validation using raw Fermi-LAT photon counts rather than the smoothed emission model. As exposure accumulated over 40 weeks of mission data, the predictive performance climbed steadily and stabilized, and the negative residuals in the predicted maps turned out to coincide with 86 known gamma-ray point sources cataloged by Fermi, including pulsars, blazars, and a radio galaxy. Because the machine learning model learned the diffuse background purely from gas and dust maps without ever seeing point-source information, any compact gamma-ray emitter naturally surfaces as a localized deficit, offering an independent check on standard source-detection pipelines.

Beyond its immediate results, the work positions machine learning as a physically interpretable instrument rather than a black box. The learned mapping encodes real relationships between interstellar matter, radiation fields, and cosmic-ray processes, and its failures mark precisely the regions where conventional templates are incomplete or biased. The authors propose extending the framework with low-frequency radio surveys such as the Haslam 408-megahertz map, polarized microwave channels, and hydrogen-alpha data on ionized gas, and eventually integrating X-ray, ultra-high-energy gamma-ray, and neutrino observations. Such a multi-messenger, data-driven baseline could help separate standard emission from exotic components, sharpen constraints on cosmic-ray propagation, and, in an era when IceCube has detected Galactic neutrinos and LHAASO has measured PeV-scale diffuse emission, provide the empirical foundation that next-generation high-energy astrophysics will demand.

Subject of Research: Data-driven machine learning modeling of Galactic diffuse gamma-ray emission using multi-wavelength Planck observations

Article Title: Data-driven modeling of Galactic diffuse emission with multi-wavelength observations

Article References: Data-driven modeling of Galactic diffuse emission with multi-wavelength observations. (n.d.). https://doi.org/10.1140/epjc/s10052-026-16408-2

Image Credits: AI Generated

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

Keywords: Galactic diffuse emission, machine learning, Planck satellite, Fermi-LAT, cosmic rays, gamma-ray astronomy, interstellar medium, inverse Compton, hadronic emission, GALPROP, Fermi Bubbles, multi-messenger astrophysics

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 controversy began in March 2025, when a group led by G. Kocheril reported measurements suggesting that cyclic C6H5+ is essentially inert toward the key molecules it would need to consume in order to build larger aromatic structures, including molecular hydrogen and acetylene, C2H2. If phenylium truly refused to react with acetylene, the ion could not serve as the gateway to naphthalene-like structures and the wider family of polycyclic aromatic hydrocarbons, or PAHs, that astronomers believe account for a substantial fraction of the carbon in the galaxy. The 2025 paper went so far as to argue that bottom-up interstellar aromatic ring formation effectively terminates at C6H5+, a conclusion that threatened to upend decades of astrochemical modelling.

That conclusion mattered because the formation of the first aromatic ring is widely regarded as the critical bottleneck in PAH growth. Once a single six-membered carbon ring exists, subsequent attachment of acetylene units can, in principle, build fused ring systems that grow into the large aromatic molecules whose infrared emission signatures pervade the Milky Way and beyond. Models of dark molecular clouds and ionized regions have long included phenylium-mediated pathways as a plausible route across this bottleneck. Removing that route would leave astrochemists scrambling to explain how the observed abundance of aromatic material could arise at all.

The new study, led by Jean-Christophe Loison of the Institute of Molecular Sciences at CNRS and the University of Bordeaux, together with Corentin Rossi and Ugo Jacovella of Université Paris-Saclay and colleagues at the Institute of Physical Chemistry and the SOLEIL synchrotron, set out to test the claimed unreactivity directly. Writing as a Matters Arising contribution in Nature Astronomy, the team combined laboratory measurements with high-level quantum chemical calculations and astrochemical modelling to re-examine how phenylium behaves when it encounters acetylene under conditions mimicking those in interstellar clouds, where temperatures hover near ten kelvin and densities are vanishingly low by terrestrial standards.

Their central result is striking: phenylium reacts efficiently with acetylene through a barrierless mechanism. In chemical kinetics, a barrierless reaction is one whose potential energy surface descends smoothly from reactants to products without any energetic hill that the colliding species must climb. Such reactions proceed at essentially every collision, governed only by how fast the two partners can find each other, which makes them extraordinarily effective in cold environments where molecules lack the thermal energy to overcome activation barriers. The team mapped the energy pathways along the reaction coordinate, showing that the association of C6H5+ with C2H2 leads to a stabilized intermediate that can proceed onward to products that incorporate the two-carbon unit into the growing carbon skeleton — precisely the chemistry needed to advance from one aromatic ring toward two.

A key technical strength of the work lies in how the phenylium ions were produced and characterized. The C6H5+ cation is a peculiar species: it possesses a low-lying triplet electronic state close in energy to its singlet ground state, and it can exist in both cyclic, ring-retaining forms and acyclic isomers with rearranged carbon skeletons. Earlier gas-phase studies, dating back to ion cyclotron resonance experiments in the 1970s and flow-tube measurements of acetylene ion chemistry in the 1980s, had already hinted at rich ion–molecule chemistry involving phenylium, but the identity and internal energy of the ions under study have often been uncertain. In the new experiments, the researchers exploited tunable synchrotron radiation to ionize a suitable precursor and measured the photon-energy dependence of C6H5+ production, allowing them to disentangle which isomers were being formed and to probe the reactivity of the cyclic cation specifically.

By monitoring how the phenylium signal decayed in the presence of acetylene, and by measuring the appearance of product ions as a function of photon energy, the team could extract rate behaviour that contradicts the earlier report of inertness. The measurements indicate that the cyclic cation, when prepared cleanly, participates in fast ion–molecule reactions with acetylene, consistent with a barrierless entrance channel predicted by the accompanying quantum chemical calculations. Loison performed the electronic structure computations that characterize the potential energy surface, identifying the intermediates and transition states that connect the initial association complex to the final products, while also running astrochemical models to assess what the revised rates mean for aromatic chemistry in interstellar clouds.

The theoretical picture also helps explain how the earlier discrepancy could have arisen. Phenylium ions are prone to isomerization and to fragmentation when produced by energetic ionization processes, and different production methods can yield mixtures of cyclic and acyclic C6H5+ with different internal energies. Acyclic isomers are known to display distinct reactivity patterns, and contamination of an ion population by unreactive or differently reactive isomers can mask the true behaviour of the cyclic cation. The photon-energy-resolved approach adopted by the French team provides a way to control for this, and their data suggest that earlier conclusions about phenylium unreactivity may have been confounded by such isomeric and energetic effects rather than reflecting an intrinsic property of the aromatic cation.

The astrophysical implications are immediate. In the cold gas of dark molecular clouds, where carbon-chain chemistry thrives and where observations have revealed surprising abundances of aromatic molecules, the availability of a fast, barrierless route from the first aromatic ring to larger structures changes the predicted chemistry substantially. Sensitivity analyses published in 2024 by Byrne and colleagues had already shown that models of aromatic chemistry in dark clouds are highly sensitive to the assumed gas-phase kinetics, meaning that a single rate coefficient can swing predicted abundances of benzene derivatives and PAH precursors by orders of magnitude. Reinstating efficient C6H5+ + C2H2 reactivity therefore reopens a pathway that recent models had been forced to close, and it strengthens the case that ion–molecule chemistry, rather than chemistry on the surfaces of dust grains alone, can seed the aromatic inventory of the interstellar medium.

The episode is also a reminder of how demanding laboratory astrochemistry can be. Ions as short-lived and isomer-sensitive as C6H5+ must be generated, transported and interrogated with exquisite control, and small differences in experimental design can produce apparently conflicting answers. The new work does not merely assert that the earlier measurements were wrong; it provides an independent, energy-resolved dataset, openly archived on Zenodo along with the data-reduction scripts, so that other groups can scrutinize and extend the analysis. Combined with the theoretical characterization of the reaction’s potential energy surface, the study offers a coherent physical picture: the phenylium cation, far from terminating bottom-up aromatic growth, is an efficient springboard from the first ring toward the polycyclic structures that glow across the galaxy. As infrared observatories continue to uncover aromatic signatures in ever more exotic environments, from dark clouds to protoplanetary disks, the rehabilitation of phenylium ensures that models of cosmic aromatic chemistry have a firmer — and faster — foundation on which to build.

Subject of Research: Laboratory and theoretical study of phenylium cation reactivity with acetylene under interstellar conditions

Article Title: Evidence for phenylium reactivity under interstellar-relevant conditions

Article References: Loison, J.-C., Rossi, C., Solem, N., Thissen, R., Romanzin, C., Alcaraz, C., & Jacovella, U. (2026). Evidence for phenylium reactivity under interstellar-relevant conditions. Nature Astronomy. https://doi.org/10.1038/s41550-026-02973-9

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02973-9

Keywords: phenylium, interstellar medium, polycyclic aromatic hydrocarbons, acetylene, ion–molecule reactions, laboratory astrophysics, aromatic ring formation, reaction kinetics, molecular clouds, astrochemistry, quantum chemistry, Nature Astronomy

The dispute began when a comment by J.-C. Loison and colleagues, published alongside the reply, presented results from an ion–molecule reaction experiment in which phenylium was generated by photodissociation of a precursor molecule. In that experiment, the resulting C6H5+ ions were found to react readily with acetylene, C2H2. According to the comment, this directly contradicts the Colorado group’s earlier publication, which reported that phenylium produced through sequential bottom-up ion–molecule reactions with acetylene was unreactive toward both acetylene and molecular hydrogen. The comment’s authors also calculated, using density functional theory, a barrierless pathway leading from phenylium and acetylene to a C8H7+ adduct, reinforcing their claim that the reaction should occur. They further argued that because their measurements operate in the single-collision regime, their conditions are highly relevant to the cold, collision-poor environment of the interstellar medium.

The stakes of this debate are considerable. Polycyclic aromatic hydrocarbons are thought to account for a substantial fraction of the carbon in the universe, and understanding how the first aromatic ring forms and grows from simple molecules such as acetylene is a central problem in astrochemistry. If phenylium reacts efficiently with acetylene, ring growth can continue, building larger hydrocarbons that eventually become the seeds of interstellar aromatic clouds. If it does not, the bottom-up assembly of aromatic rings stalls at the six-carbon stage, and astrochemical models must be revised accordingly. Which experiment is right changes how scientists simulate the chemistry of dark molecular clouds and circumstellar envelopes.

In their reply, the Colorado researchers make a striking concession that reframes the entire controversy: they are not surprised by the data presented by Loison and colleagues. Decades of similar experiments, they note, have consistently demonstrated that a reactive form of C6H5+ is produced through methods like photodissociation and dissociative ionization. That behavior stands in sharp contrast to observations of an unreactive species formed through ion–molecule chemistry, a pattern documented in the literature going back to the 1980s. The discrepancy between reactive and unreactive outcomes has appeared repeatedly across different laboratories, techniques and generations of instrumentation, making it one of the longest-running puzzles in gas-phase ion chemistry.

Historically, the debate over these divergent results has centered on molecular structure. When two experiments produce ions with the same chemical formula but different reactivity, the simplest explanation is that the experiments are producing different structural isomers. In this case, the acyclic form of C6H5+ would be expected to react with acetylene, while the ring-shaped phenylium cation might resist the reaction. Loison and colleagues adopted precisely this line of reasoning, suggesting that the simplest explanation for the discrepancy is that the Colorado experiments have produced a higher-energy, acyclic C6H5+ isomer instead of the true phenylium cation. Under this interpretation, the unreactivity reported by the Colorado team would be an artifact of making the wrong molecule.

The reply firmly rejects that interpretation. Kocheril, Zagorec-Marks and Lewandowski state their position plainly: they believe the structure of the C6H5+ produced in their experiments is the same as that produced by Loison and colleagues. What differs, they argue, is the internal energy of the ions. Two ionic species can share an identical atomic connectivity while carrying very different amounts of vibrational and electronic excitation, and that internal energy can govern whether a barrierless association reaction proceeds on the timescales probed by an experiment. An ion born with substantial internal energy can find reaction pathways that a cold, relaxed counterpart cannot access, even if both species are structurally indistinguishable.

This energy-based explanation carries particular weight because of how the ions are made. In the Colorado experiments, phenylium is built from the bottom up, assembled through a sequence of ion–molecule reactions with acetylene at low temperature. Such gentle assembly tends to leave ions in comparatively low internal energy states. In the commenting authors’ experiments, the ion is produced through photodissociation of a precursor, a process that can deposit significant energy into the fragment. The reply includes a schematic representation of the energetics of C6H5+ produced through dissociative ionization, underscoring the authors’ argument that the production method itself dictates the internal energy content and therefore the observed reactivity. The same ion, they contend, can behave differently depending on how hot it is when it meets a collision partner.

The reply also addresses the claim that single-collision experiments are the most faithful proxies for interstellar conditions. While the single-collision regime does mimic the extreme isolation of the interstellar medium, where particles collide rarely, the Colorado authors point out that decades of experimental history show reactive C6H5+ arising from precisely the kind of photodissociation-based production used in the comment. If that production route systematically yields an internally excited ion, then reactivity measured under those conditions may reflect the excess energy of the ion rather than the intrinsic low-temperature chemical behavior of relaxed phenylium in space. The unreactive behavior observed for ions assembled through ion–molecule reactions, by contrast, may better represent the quiescent chemistry of cold molecular clouds. Resolving which regime matters most for astrochemical models will require disentangling internal energy from structural identity, a task that demands careful spectroscopic characterization of the ions involved.

The authors point to a growing body of independent evidence relevant to this question, including recent photoelectron spectroscopy work that directly observed the fundamental arylium species and measured its singlet–triplet gap, as well as theoretical studies of the isomers of C6H5+ and their formation pathways in the interstellar medium. For readers following the debate, the exchange is a vivid reminder that in laboratory astrophysics, how you make a molecule can matter as much as what molecule you make. The Colorado team refers readers to their original publication for a full account of why they believe they have produced the phenylium isomer, and the exchange is likely to prompt follow-up experiments designed to measure the internal energy distributions of C6H5+ ions directly. Until then, the fate of aromatic ring growth in the coldest corners of the galaxy hangs on a subtle question of molecular temperature.

Subject of Research: Laboratory study of the internal-energy-dependent reactivity of the phenylium cation C6H5+ and its role in interstellar aromatic ring formation

Article Title: Reply to: Evidence for phenylium reactivity under interstellar-relevant conditions

Article References: Kocheril, G. S., Zagorec-Marks, C., & Lewandowski, H. J. (2026). Reply to: Evidence for phenylium reactivity under interstellar-relevant conditions. Nature Astronomy. https://doi.org/10.1038/s41550-026-02972-w

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02972-w

Keywords: phenylium, astrochemistry, interstellar medium, ion–molecule reactions, aromatic ring formation, acetylene, C6H5+, laboratory astrophysics, internal energy, Nature Astronomy, isomers, molecular clouds

Turbulence is not a cosmetic feature of galaxies; it is a controlling one. The random, supersonic motions of interstellar gas set the thickness of galactic disks, regulate how molecular clouds collapse into new stars, and determine how efficiently matter converts into luminous stellar populations. Without a continuous supply of energy, turbulence in a galaxy like Andromeda would decay on timescales of only a few tens of millions of years, far shorter than the age of the disk itself. Something must constantly re-energize the gas, and supernovae have topped the list of candidate drivers since the earliest theoretical treatments of the multiphase interstellar medium. Alternative mechanisms, including gravitational instabilities in the rotating disk and the magnetorotational instability, have also been proposed, and disentangling their contributions observationally has proven exceptionally difficult.

The key to the new result lies in superbubbles, the gigantic cavities that clusters of massive stars carve out of the neutral hydrogen gas around them. When a group of hot, short-lived stars forms together, the combined winds and subsequent core-collapse supernovae of its members blow a common expanding shell into the surrounding medium. These shells, sometimes hundreds of light-years across, act as fossil records of clustered stellar feedback: their sizes, expansion velocities, and ambient gas densities encode how much kinetic energy the parent star cluster deposited into the disk. Because supernovae in galaxies tend to occur in clusters rather than in isolation, superbubbles are the natural conduits through which stellar feedback feeds galactic-scale turbulence, making a complete census of them the decisive observational test.

Until now, no such dynamically complete, galaxy-wide census existed. Earlier surveys of hydrogen holes and shells in Andromeda and other nearby galaxies, dating back to work in the 1980s, identified cavities in the gas but lacked the sensitivity and velocity resolution to measure expansion reliably across a whole disk, leaving the energy budget unsettled. The new study overcame this limitation by combining the extraordinary sensitivity of FAST, the world’s largest single-dish radio telescope, with the fine angular resolution of the Jansky Very Large Array. The combination was essential: single-dish data recover the diffuse, large-scale hydrogen emission that interferometers miss, while interferometric data resolve the fine structure of individual shells. An image-fusion technique merged the two data sets into a single, dynamically complete hydrogen data cube of Andromeda, capturing structures on all relevant scales simultaneously.

From this combined data set, the team identified 118 superbubbles distributed across the entire disk of M31, each appearing as a coherent shell in both space and velocity. The dynamical ages of the shells extend up to roughly 40 million years, a range that matches theoretical expectations for the duration of supernova activity within a star cluster, since the most massive members explode within a few million years while lower-mass stars detonate tens of millions of years later. This consistency between the observed age distribution and stellar-evolution models strengthens the interpretation that the shells are indeed powered by clustered supernovae rather than by other processes such as infalling clouds or galactic-scale instabilities.

The crucial quantitative step was an energy accounting performed two independent ways. First, from the measured sizes, expansion velocities, and surrounding gas densities of the superbubbles, the researchers calculated the rate at which supernovae inject kinetic energy into the neutral medium. The inferred injection rates span 10^49 to 10^51.5 erg per cubic kiloparsec per million years. Second, from the same hydrogen data, they measured the turbulent velocity dispersion of the gas and derived how quickly turbulent energy dissipates at each location in the disk, using established scaling relations for supersonic, magnetized turbulence. If supernovae truly power the turbulence, these two independently determined rates should agree, both in total magnitude and in how they vary from place to place across the galaxy.

They do agree, and remarkably well. The kinetic-energy-injection rates inferred from the superbubble population closely match the turbulence dissipation rates derived from the gas kinematics, not only in overall magnitude but also in their spatial distribution across the Andromeda disk. Regions where the shells inject more energy are precisely the regions where the gas exhibits the strongest turbulent motions. This point-by-point correspondence is far more constraining than a global average, because it rules out a coincidence in which supernovae supply the right amount of energy somewhere in the galaxy while a different mechanism actually drives the local turbulence. The result demonstrates that clustered supernova feedback alone is sufficient to sustain galactic-scale turbulence in a giant spiral galaxy.

The implications extend well beyond Andromeda. Because M31 is the nearest giant spiral and a close analogue of the Milky Way, the result provides the strongest direct evidence to date that our own galaxy’s turbulent interstellar medium is likewise maintained by the death throes of massive stars. Turbulence, in turn, feeds back into star formation: it both prevents gas from collapsing too quickly into stars and concentrates density enhancements that seed new star-forming clouds. A galaxy’s structure and evolutionary trajectory therefore depend on this feedback loop, and models of galaxy formation and evolution can now anchor their prescriptions for supernova-driven turbulence to a directly measured, observationally verified energy budget rather than to theoretical assumption alone.

The study also showcases the power of combining complementary radio facilities. FAST’s collecting area delivers sensitivity to faint, extended hydrogen emission at a level no other instrument can reach, while the Very Large Array contributes the sub-arcminute resolution needed to resolve individual shells hundreds of parsecs away in a neighboring galaxy. The data, the 118-object superbubble catalogue, and the analysis code have been released publicly, allowing other researchers to scrutinize the classification, refine the energy estimates, and extend the method to additional galaxies. As similar combined surveys target more spirals, astronomers will be able to test whether supernova-driven turbulence dominates universally or whether gravitational and magnetic mechanisms take over in particular environments, such as low-star-forming outer disks or violently interacting systems.

For decades, the image of galaxies as serene, slowly rotating pinwheels has coexisted with the reality that their gas is in constant, violent motion, churned by forces whose origin remained unproven. With 118 superbubbles now mapped across Andromeda and their energy output shown to balance the turbulent dissipation of the entire disk, that origin is no longer a hypothesis but a measurement. The explosions of massive stars, it turns out, are not merely the spectacular endings of stellar lives; they are the beating heart that keeps entire galaxies stirred, structured, and capable of forming the next generation of stars.

Subject of Research: Observational evidence that clustered supernova feedback, traced by hydrogen superbubbles, sustains galactic-scale turbulence in the Andromeda galaxy.

Article Title: Supernova origin of galactic turbulence revealed by superbubbles

Article References: Meng, F., Tsai, C.-W., Wu, J., Jiao, S., Mac Low, M.-M., Zhang, Z.-Y., Saintonge, A., Li, H., Li, Z., Wang, J., Wang, L., Xu, H., Yang, Y., Zhang, K., Li, R., & Li, D. (2026). Supernova origin of galactic turbulence revealed by superbubbles. Nature Astronomy. https://doi.org/10.1038/s41550-026-02981-9

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02981-9

Keywords: supernovae, superbubbles, galactic turbulence, Andromeda galaxy, M31, neutral hydrogen, FAST telescope, Jansky Very Large Array, interstellar medium, star formation, stellar feedback, Nature Astronomy