Peer Revenue Returns subject to peer review

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics

astro-ph Watch Astrophysics

Tuesday 24 March

Solar and Stellar Astrophysics

Magnetic flux systems involved in the May 2024 solar energetic events from AR 13664 inferred through coronal dimmings

Amaia Razquin, Karin Dissauer, Astrid M. Veronig et al. (+1 more)

Coronal dimmings are transient depletions of coronal plasma observed in extreme ultraviolet and soft X-rays, interpreted as low-corona signatures of coronal mass ejections (CMEs). Their evolution is closely linked to CME dynamics, flare reconnection, and large-scale coronal magnetic restructuring. In May 2024, active region (AR) 13664 produced 66 > M-class flares and multiple fast CMEs that caused the strongest geomagnetic storm since 2003. We analysed 16 on-disc dimmings observed in SDO/AIA 211Ådata between May 1-14, 2024. We extracted dimmings using logarithmic base-ratio thresholding, and derived magnetic properties from SDO/HMI line-of-sight magnetograms. We identified flare ribbons in AIA 1600Ådata via adaptive thresholding, and computed reconnection fluxes from radial magnetic field maps. We examined the magnetic flux systems involved in the eruptions using PFSS and NLFF extrapolations. We found correlations between dimming and flare ribbon properties stronger than previously reported. The dimming morphology evolved systematically, with southward expansion before May 9 and northward afterward, coinciding with a shift in flare ribbon locations between two mayor east-west polarity inversion lines (PILs). These changes imply the presence of two distinct magnetic domains. The PFSS extrapolations showed that southward (northward) dimmings are mainly strapping flux dimmings with magnetic field lines vaulting above the southern (northern) PIL. The final extent of the dimmings was then given by the exterior flux involved in the eruption via stretching and reconnection. One event revealed an extended quiet-Sun dimming likely triggered by EUV wave-induced field opening.
Read the full abstract

All authors Amaia Razquin, Karin Dissauer, Astrid M. Veronig, Graham Barnes

High Energy Astrophysical Phenomena · Cross-listed from gr-qc

Radial Oscillations of Viscous Stars

Lennox S. Keeble and Jaime Redondo-Yuste

Oscillation modes of neutron stars, a key target for third-generation gravitational wave detectors, encode key information about their constituent nuclear matter. In this work, we study the effect of viscosity on oscillations of cold, polytropic, spherically symmetric neutron stars. We focus on purely radial oscillations and work perturbatively to linear order within two hydrodynamic frameworks: the acausal covariant generalization of the Navier-Stokes equations proposed by Eckart, and the causal generalization formulated by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). We find that viscosity damps the radial modes on millisecond timescales and induces fractional shifts in the oscillation frequency which increase both with the compactness and viscosity of the star, reaching up to the percent level for the fundamental mode with bulk viscosities $\zeta\sim10^{30}\mathrm{g}/\mathrm{cm}/\mathrm{s}$. For more viscous stars, the oscillation frequency decreases, becoming zero (i.e., an overdamped mode) for $\zeta\gtrsim10^{31}\mathrm{g}/\mathrm{cm}/\mathrm{s}$. We also study the linear threshold of gravitational collapse. Consistent with recent analytic results in the zero heat conductivity limit, we find that viscosity in Eckart theory cannot stabilize an unstable inviscid star. We provide numerical evidence that viscosity in BDNK theory is similarly unable to prevent gravitational collapse, but it slightly modifies the threshold of collapse. Overall, our results advance our understanding of the impact of viscosity on the oscillation modes of neutron stars, a key component of viscous asteroseismology with next-generation gravitational wave detectors.
Read the full abstract

Astrophysics of Galaxies

SPHEREx mapping of diffuse PAH and H II emission in the Galactic plane

Giulia Murgia, Ari J. Cukierman, Brandon S. Hensley et al. (+29 more)

We present preliminary SPHEREx maps of diffuse Galactic emission tracing polycyclic aromatic hydrocarbons (PAHs) and ionized hydrogen gas, and we study their relationship across the Galactic plane. Since its launch in early 2025, the SPHEREx space telescope has been conducting an all-sky near-infrared spectral survey from 0.75 to 5.0 microns. We produce a large-scale map of the 3.3-micron PAH emission feature, which is bright and detectable throughout the Galactic plane, and find a strong correlation with the thermal dust radiance measured by Planck. We also trace ionized hydrogen gas by producing a map of Brackett-alpha emission at 4.05 microns. By combining the two maps, we identify extended shells of PAH emission associated with photodissociation regions surrounding ionized gas. We construct a PAH abundance map and find a significant anticorrelation between PAH abundance and ionized hydrogen, indicating systematic PAH depletion within ionized gas regions across the Galactic plane and demonstrating that ionizing radiation is a dominant driver of PAH abundance variations. These early SPHEREx results provide a large-scale view of PAHs and ionized hydrogen and preview the capability of the mission to map diffuse emission in the interstellar medium.
Read the full abstract

All authors Giulia Murgia, Ari J. Cukierman, Brandon S. Hensley, Matthew L. N. Ashby, James J. Bock, Tzu-Ching Chang, Shuang-Shuang Chen, Yun-Ting Cheng, Yi-Kuan Chiang, Asantha Cooray, Brendan P. Crill, Olivier Doré, C. Darren Dowell, Andreas L. Faisst, Joseph L. Hora, Howard Hui, Miju Kang, Jae Hwan Kang, Phil M. Korngut, Dennis Lee, Jeong-Eun Lee, Bomee Lee, Carey M. Lisse, Daniel C. Masters, Gary J. Melnick, Mary H. Minasyan, Chi H. Nguyen, Roberta Paladini, Volker Tolls, Robin Y. Wen, Michael W. Werner, Michael Zemcov

Cosmology and Nongalactic Astrophysics

A velocity-dependent two-scale model for cosmic string networks with small-scale structure

Teresa O. Miranda and Lara Sousa

We develop a semi-analytical model to describe the cosmological evolution of networks of cosmic strings with small-scale structure, by extending the velocity-dependent one-scale model to include an additional lengthscale describing the typical interkink density. We study the impact of the different physical processes involved in the production and removal of small-scale structure from cosmic strings on the attainment of a full linear scaling regime, in which the characteristic lengths of the network and of small-scale structure evolve proportionally to physical time and the root-mean-squared velocity of the network remains constant. We find, using this novel velocity-dependent two-scale model, that quite generally small-scale structure does not prevent the attainment of a linear scaling regime since, even if not enough kinks are carried away when loops are chopped from the network, gravitational backreaction is generally enough to ensure that the interkink density scales. We find, however, that this regime is characterized by a smaller energy density and root-mean-squared velocity when compared to strings without small-scale structure and that this reduction may be significant when scaling is maintained by gravitational backreaction. In this case, we also find that, before reaching full scaling, the network should evolve in a transient quasi-scaling regime, in which its evolution is very similar to that of cosmic strings without small-scale structure.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Reducing cosmological degeneracies by combining multiple classes of LISA gravitational-wave standard sirens

Danny Laghi, Nicola Tamanini, Alberto Sesana et al. (+6 more)

We present the first joint gravitational-wave cosmological inference with LISA extreme mass-ratio inspirals at $z\lesssim1$ (galaxy redshifts) and massive black hole binaries at $z\gtrsim1$ (electromagnetic counterparts). Combining these standard sirens reduces cosmological degeneracies and yields competitive constraints on the Hubble constant $H_0$ and the dark-energy equation-of-state parameter $w_0$. This highlights LISA's potential for late-time cosmology across a broad redshift range with systematics distinct from electromagnetic distance indicators.
Read the full abstract

All authors Danny Laghi, Nicola Tamanini, Alberto Sesana, Jonathan Gair, Enrico Barausse, Chiara Caprini, Walter Del Pozzo, Alberto Mangiagli, Sylvain Marsat

Astrophysics of Galaxies

Revealing the nature of the starburst galaxies in the $z=2.4$ overdensity HATLAS J0849

Melanie Kaasinen, Francesca Rizzo, Francesco Valentino et al. (+4 more)

Today's most massive ellipticals are proposed to originate from starbursting galaxies in $z\gtrsim2$ overdensities. To discern what triggers these starbursts, and their $z=0$ descendants, we performed a detailed case study of five gas-rich galaxies in the $z=2.41$ overdensity, HATLAS J084933.4+021443. Using 0.15" resolution CO(4-3), [C I] 1-0, and dust-continuum observations, we characterised their cold gas morphology and kinematics. We find two rotating discs, W and C, both exhibiting non-axisymmetric radial gas motions (consistent with bars). Of the two extreme starbursts, W is a lopsided, rotation-dominated disc with a rotation velocity of $\sim520$ km s$^{-1}$, whereas T is most likely a late-stage merger. Combined with recent studies, we find that $\gtrsim42\%$ of gas-rich, massive starbursts in overdensities are rotation-dominated discs, a fraction not yet systematically reproduced by galaxy evolution models. Beyond $z=1$, disc galaxies with rotation velocities of $>400$ km s$^{-1}$ reside almost exclusively in overdensities, consistent with early mass assembly in dense environments. By comparing to local early-type galaxies with cold gas discs, we confirm that these systems already reside in halos comparable to the most massive $z\sim0$ ellipticals at the centres of groups and clusters. Despite their extreme star-formation rates, these discs lie on the same $\sigma-$SFR locus as lower-SFR field galaxies, implying that stellar feedback remains the dominant turbulence driver. We postulate that this is because inflowing gas is effectively transported through ordered streaming, such that only a small fraction of kinetic energy feeds disc-wide turbulence.
Read the full abstract

All authors Melanie Kaasinen, Francesca Rizzo, Francesco Valentino, Cecilia Bacchini, Jianhang Chen, Takafumi Tsukui, Aristeidis Amvrosiadis

Astrophysics of Galaxies

A hydrodynamical CLONE of the Virgo cluster: II. Confronting observed and synthetic galaxy population twins in a dense environment

Jenny G. Sorce, Sean L. McGee, Yohan Dubois et al. (+3 more)

Galaxy clusters offer powerful laboratories for studying galaxy evolution in dense environments. In this context, the Constrained LOcal and Nesting Environment (CLONE) project provides a zoom-in hydrodynamical simulation of the Virgo cluster, including active galactic nucleus and supernova feedback, with a resolution down to 350~pc, designed to mirror Virgo's observed properties. Previous work showed that this replica and Virgo share the same history, mass, and luminosity distributions including the central M87. This study examines several observational relations extending to lower stellar masses than previous synthetic-population studies: star formation density, (specific) star formation rate, metallicity, and quenched fraction of galaxies as a function of stellar mass and cluster-centric distance. Despite the slightly low metallicity and the sufficiently high quenched fraction, simulated galaxies reproduce key observational trends even without averaging or accounting for observational uncertainties, aside from the consideration of projection effects: At fixed stellar mass, cluster galaxies form fewer stars than field counterparts; Most galaxies are quenched, except for intermediate-mass or isolated galaxies; Low-mass galaxies are highly quenched, thus implying a sharp metallicity drop, and low metallicity does not imply youth; Quenching occurs earlier for the most massive and the smallest galaxies than for those of intermediate mass, at least until they enter the cluster; Quenched galaxies have undergone dark matter stripping; Gas depletion drives quenching, especially in low-mass galaxies and the farther from the cluster center they are. Overall, the synthetic population jointly reproduces multiple observational trends, making it a valuable tool to probe processes from jellyfish galaxies to cluster-core gas dynamics.[Abridged]
Read the full abstract

All authors Jenny G. Sorce, Sean L. McGee, Yohan Dubois, Jérémy Blaizot, Alexander Knebe, Gustavo Yepes

Astrophysics of Galaxies

The impact of disc disruption on Milky Way satellite counts

Mark R. Lovell, Alexander H. Riley and Isabel Santos-Santos

Estimates for the total number of Milky Way (MW) satellites are often generated from a combination of the observed number of satellites in surveys, adjustments for the completeness of those surveys, and theoretical expectations from halo assembly modelling. One of the features of this modelling is disruption by the MW stellar disc. We examine the effect of degrees of disc disruption on inferred satellite counts, by means of an N-body simulation of a MW-mass halo plus a toy model for this disruption. We use a fictional all-sky survey to show that high resilience to disc disruption predicts small populations of satellites that are radially very concentrated around the central galaxy and are hosted by massive subhaloes, while low resilience predicts many more satellites with a less concentrated radial distribution and hosted within less massive subhaloes. We show that the most massive subhaloes are particularly susceptible to disruption due to their radial orbits, and in their putative absence galaxy formation must occur in lower mass haloes that have a shallower radial number density profile. We then demonstrate this phenomenon for a combination of the Pan-STARRS and DES surveys. It is therefore necessary to account for uncertainty in the disc disruption radius when making predictions for MW satellite distributions.
Read the full abstract

High Energy Astrophysical Phenomena

The Environments of Luminous Fast Blue Optical Transients: Evidence for a Compact Object and Wolf-Rayet Star Merger Origin

Anya E. Nugent, V. Ashley Villar, Brian D. Metzger et al. (+5 more)

We present a comprehensive analysis of the host galaxies of 11 luminous fast blue optical transients (LFBOTs). We model new and archival host photometry and spectroscopy with Prospector. We determine that all LFBOT hosts are actively star-forming with recent bursts of star formation and have a median stellar mass of $\log(M_*/M_\odot)=9.61^{+0.74}_{-1.61}$, present-day star formation rate SFR=$0.99^{+14.85}_{-0.95}$~$M_\odot$~yr$^{-1}$, and gas-phase oxygen abundance metallicity 12+log(O/H)=$8.59^{+0.18}_{-0.22}$. To contextualize these results, we compare them to the host properties of Hydrogen-poor superluminous supernovae (SLSNe-I), several core-collapse supernova subtypes (CCSN; SNe Ibc, II, and Ibn) and long gamma-ray bursts (LGRBs). We find that LFBOT hosts are more star-forming than CCSN hosts, but less star-forming than SLSN-I hosts. We further show that LFBOT hosts are more metal-poor than SN Ibc and II hosts, but more metal-rich than SLSN-I and LGRB hosts. Finally, we find that, similar to SLSNe-I and unlike CCSNe and LGRBs, a large fraction of LFBOTs occur in their hosts' faintest pixel or outside their host galaxy's light. Our results indicate that LFBOTs have a massive stellar origin that do not trace active star-forming regions within their hosts and have a weaker metallicity-dependence than other extreme transients. For these reasons, we favor a compact-object and Wolf-Rayet star merger progenitor scenario over other previously proposed models, such as tidal disruption events and failed or successful CCSN. Future discoveries of LFBOTs with the Rubin observatory will help to increase their sample size and place firmer constraints on their environments and progenitors.
Read the full abstract

All authors Anya E. Nugent, V. Ashley Villar, Brian D. Metzger, Christopher L. Fryer, Eric Burns, Alexa Gordon, Danielle Frostig, Yuxin Dong

Astrophysics of Galaxies

Direct Evidence for Stellar Initial Mass Function Variation in the Milky Way

Charles L. Steinhardt, Carter Meyerhoff and Alexander J. Luening

Because direct measurements require resolved stellar populations including low-mass stars, determining the stellar initial mass function (IMF) has been a historically difficult problem even within our own Galaxy and impossible everywhere else. As a result, even though it is predicted that the IMF should vary depending upon the properties of each individual star-forming molecular cloud, it is standard to assume a Universal IMF. Using recent observations from Gaia, it is now possible to test for IMF variation using resolved stellar populations in open clusters and a parameterization that separates properties of the IMF from subsequent dynamical evolution. Here, we show that the IMF is not Universal but instead varies across individual Galactic stellar populations, reflecting evolution in the average conditions of molecular clouds over cosmic time. This evolution is consistent with the predictions of a simple astrophysical model in which the IMF is environmentally-dependent and the Milky Way reflects typical galactic behavior in recent cosmic history. Thus, observational evidence now agrees with long-standing theoretical and numerical predictions.
Read the full abstract

Astrophysics of Galaxies

Deviations from the radial acceleration relation in the central galaxies of clusters, subclusters, and groups

Michal Bílek, Florent Renaud and Srdjan Samurović

Most galaxies closely follow the radial acceleration relation (RAR), which tightly links the observed accelerations to those predicted by Newtonian gravity from visible baryonic matter. Galaxy clusters, however, deviate from this relation. Several explanations have been proposed. Some of them predict that even some galaxies in clusters should deviate, but this hypothesis remains largely untested. We test it here by analyzing acceleration profiles for 17 early-type galaxies, derived from Jeans modeling of their globular cluster systems in our older work. Our sample spans central galaxies in clusters and groups, non-central galaxies, isolated ones, and-uniquelly for this paper-centrals in galactic subclusters, which are smaller clusters being accreted by larger ones. We compare these profiles to the standard RAR for non-cluster galaxies and its counterpart for clusters. We find that isolated and non-central galaxies adhere to the standard RAR. In contrast, central galaxies of clusters, subclusters, and groups exhibit enhanced accelerations in most cases, tracing instead the cluster acceleration behavior either partly or fully. The radius at which divergence from the standard RAR begins tends to decrease with increasing group mass. These findings imply that if cluster fields depart from the standard RAR due to undetected material, it must be dynamically cold and collisionless, such as non-baryonic cold dark matter, but also compact clouds of cold gas.
Read the full abstract

Astrophysics of Galaxies

oMEGACat. X. Shedding light on the disrupted dwarf galaxy of Omega Centauri

Stefano Souza, Nadine Neumayer, Anil C. Seth et al. (+15 more)

Omega Centauri ($\omega\,$Cen) is the most massive and chemically complex star cluster in the Milky Way and is widely regarded as the surviving nuclear star cluster of an accreted dwarf galaxy. However, its parent host remains uncertain. Here, we investigate a scenario in which Sequoia, Thamnos, and Gaia–Enceladus (GE) are debris from a single disrupted progenitor, the $\omega\,$Dwarf, whose nucleus survives today as $\omega\,$Cen. Using APOGEE and GALAH abundances together with Gaia astrometry, we reconstruct the chemical structure across this progenitor adopting orbital energy as a proxy for pre-merger radius. We find that the chemically evolved (younger Al-N-He-rich) population is strongly concentrated toward the inner regions, representing a population formed after/during the merger, while the primordial population represents a dwarf-galaxy-like population, supporting a common dwarf-galaxy origin for its components. The metallicity profile shows an inverted U-shaped gradient similar to those observed in present-day nucleated dwarf galaxies. At the same time, the inner regions ($\omega\,$Cen+Thamnos) are more $\alpha$-enhanced than the outskirts, pointing to shorter and more efficient star formation and indicating that the nucleus may have assembled through the merger of inspiraling globular clusters. Neutron-capture abundances reveal a Eu-rich, r-process-dominated outskirts and inner regions enhanced in [Ba/Eu] and [La/Eu], requiring delayed enrichment and more complex chemical evolution. Finally, our analysis shows that Sequoia and Thamnos naturally fit an outside-in stripping sequence around $\omega\,$Cen, whereas the connection with GE remains unsure.
Read the full abstract

All authors Stefano Souza, Nadine Neumayer, Anil C. Seth, Zixian Wang, Callie Clontz, Maximilian Häberle, Maria S. Nitschai, Peter J. Smith, Tadafumi Matsuno, Guillaume Guiglion, Anja Feldmeier-Krause, Nikolay Kacharov, Glenn van de Ven, Jiadong Li, Mattia Libralato, Andrea Bellini, Antonino P. Milone, Mayte Alfaro-Cuello

Astrophysics of Galaxies

Phase spirals across galactic disks I: Exploring dynamical influences on winding

Kiyan Tavangar, Kathryn V. Johnston, Jason A.S. Hunt et al. (+4 more)

The vertical phase-space spirals in the Milky Way are clear evidence of disequilibrium. However, they are challenging to study because phase mixing signals evolve under the influence of many different dynamical processes and can be driven by many sources of disequilibrium. We characterize phase spirals in two simulations – one test particle and one N-body – with basis function expansions, using these to derive winding times ($T_{\rm fit}$). We find that phase spirals in the test particle simulation wind up as expected from pure phase mixing theory while those in the self-consistent simulation do not. Specifically, in the N-body simulation we find that (i) the onset of winding is delayed, (ii) the winding rate is slowed, and (iii) the rate of winding oscillates with time. The extent of these effects depends on the azimuthal action $J_\phi$ of the phase spiral region. We build some physical intuition for these effects through 1-D toy models which follow a group of co-moving stars traveling through several different evolving potentials. We find that phase spiral winding can be delayed until the group no longer moves coherently with the midplane of the (perturbed) potential and oscillates with time as the group experiences (e.g.) a breathing mode traveling through the disk. Rates of winding change as the vertical structure of the disk evolves. The modifications to winding are strongest in the inner galaxy where the disk potential dominates. We conclude that in the Milky Way, all calculations of the winding time should be interpreted as lower limits and that the most trustworthy winding times are likely in the outer disk.
Read the full abstract

All authors Kiyan Tavangar, Kathryn V. Johnston, Jason A.S. Hunt, Axel Widmark, Chris Hamilton, Michael S. Petersen, Martin D. Weinberg

Cosmology and Nongalactic Astrophysics

Investigating the radio emission in the Perseus cluster with LOFAR sub-80 MHz LBA

C. Groeneveld, R. J. van Weeren, M.-L. Gendron-Marsolais et al. (+7 more)

The Perseus cluster is a nearby cool-core galaxy cluster that hosts an archetypal radio mini-halo. Recent Low Frequency Array (LOFAR) High Band Antenna (HBA) observations at 120 - 168 MHz have revealed the presence of a giant radio halo within the cluster with a size of 1.1 Mpc enveloping the mini-halo. By exploring the spectral properties of the radio emission at low frequencies, we can gain deeper insights into the nature of this emission and improve our understanding of its origin. Here we present LOFAR Low Band Antenna (LBA) images of the cluster between 30.0 - 57.7 MHz, with a resolution of 19.2” x 15.0” and a r.m.s. noise of 3.7 mJy/beam . In our images, we detect both the mini-halo and giant radio halo. We measured the spectral indices between 44 and 144 MHz of the mini-halo and giant radio halo to be -1.34 +- 0.10, and -1.01 +- 0.11, respectively. An alternative and more direct measurement of the spectrum of the giant radio halo results in a spectral index of -1.28 +- 0.15. The discrepancy between both values is caused by the poor ionospheric conditions. In addition, we study two X-ray 'ghost cavities' in the cluster. These cavities are thought to have been produced by an older outburst from the central AGN 3C 84. We measure a spectral index between 44 and 144 MHz for the radio plasma in these cavities of -1.86 +- 0.12 and -1.90 +- 0.12 for the northwest and southern ghost cavities, respectively. Furthermore, by including VLA 352 MHz data, we find that the spectrum steepens at higher frequencies. These results are consistent with the ghost cavities being filled with old and aged radio plasma. We also detect the tailed radio galaxies NGC 1265 and IC 310. In our analysis, these sources show signs of spectral steepening along their tails.
Read the full abstract

All authors C. Groeneveld, R. J. van Weeren, M.-L. Gendron-Marsolais, E. Osinga, A. Botteon, F. de Gasperin, M. Cianfaglione, G. di Gennaro, G. Brunetti, R. Cassano

Cosmology and Nongalactic Astrophysics

Derivation of the Kompaneets equation using the boost operator approach

Alex Hoey, Jacob Long and Jens Chluba

The repeated scattering of photons by thermal electrons at low temperatures is described by the Kompaneets equation and its generalized forms that include anisotropies and higher order temperature corrections. In this work, we use the boost operator approach to derive the related expressions in a transparent way that showcases the generality of the formalism and its application to radiative transfer problems. We consider the simplest form of the Kompaneets equation for the scattering in isotropic media at the leading order in the electron temperature and then include anisotropies in the photon field, reproducing previously obtained expressions for the evolution equations. For this we use expressions for the scattering operator in the electron rest frame up to first order in the electron recoil, O(h nu/m_e c^2), but then work at all orders in the electron momentum, p, as easily obtained with the boost operator approach. This shows how specific transformation rules can be formulated that allow simplification of the otherwise cumbersome and repetitive calculations. We also confirm the expressions for higher order temperature corrections in isotropic media, highlighting the validity of the approach presented here. As part of the derivation, we find expressions for the boost operator in general boost directions which we believe will also be useful in other applications of the formalism.
Read the full abstract

Earth and Planetary Astrophysics · Cross-listed from physics.hist-ph

Mars in the Australian Press, 1875-1899. 1. Interpretation, Authority and Planetary Science

Richard de Grijs

[Abridged] In the late nineteenth century, Mars emerged as one of the most intensively reported astronomical objects in the popular press, driven by favourable oppositions, improved telescopic capabilities and growing speculation regarding planetary habitability. I examine how Mars was interpreted in Australian newspapers between the 1870s and 1899, focusing on the ways in which astronomical knowledge was framed, contextualised and debated within a colonial media environment. Drawing on a large collection of digitised newspaper articles, I analyse how observational authority, instrumental credibility and individual expertise were harnessed in press reporting. The paper situates Australian Mars coverage within a global network of scientific communication dominated by metropolitan centres in Europe and North America, while highlighting the distinctive role played by southern-hemisphere visibility. Australian observatories and observers were frequently positioned as contributors of confirmatory observation rather than interpretive leadership, reinforcing a pattern of locally grounded but internationally oriented scientific engagement. The analysis traces a shift from early emphasis on disciplined observation and measurement to later periods characterised by contested interpretations, particularly surrounding the so-called Martian "canals" and the speculative claims advanced by personalities such as Percival Lowell in the USA. By examining how newspapers mediated between observational astronomy, engineering analogies and popular imagination, this study contributes to a broader understanding of how planetary science entered public discourse beyond metropolitan centres. In doing so, it underscores the active role of colonial newspapers in shaping scientific meaning and situates Australian Mars reporting within the wider history of nineteenth-century astronomical culture.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Solving the Cosmic Coincidence Problem: The Locally Pumped Dark Energy Model

Carlo R. Contaldi and Mauro Pieroni

We propose the Locally Pumped Dark Energy (LPDE) mechanism in which cosmic acceleration is triggered by the emergence of non-linear dark matter structure. In an effective-field-theory description, coarse-graining over the density contrast profile, whose short-wavelength modes grow during halo formation, induces a shift in the local equilibrium point of a second, sufficiently heavy scalar field $\chi$. At early times, the pump mechanism is negligible and $\chi$ remains fixed at the origin, contributing no DE. As structures form, the equilibrium value of $\chi$ is locally displaced within halos, generating a vacuum energy whose global contribution, in a mean-field picture, is controlled by the halo volume filling factor. If the $\chi$ field is sufficiently heavy, with a Compton wavelength limited by halo scales, its response is localised, and spatial gradients are exponentially suppressed on large scales. After volume-averaging over the halo population, the resulting contribution on large scales behaves as a homogeneous DE component. Using the halo mass function of a fiducial $\Lambda$CDM cosmology, we show that vacuum-energy domination generically emerges at $z\sim\mathcal{O}(1)$, naturally correlating cosmic acceleration with structure formation. For reference, we present an explicit realisation of such a mechanism and show that, by naturally featuring a transient acceleration epoch, it can be in excellent agreement with the most recent cosmological data, including the Dark Energy Spectroscopic Instrument (DESI).
Read the full abstract

Cosmology and Nongalactic Astrophysics

Microlensing by Cluster of Primordial Black Holes

K.A. Toshchenko, P.V. Baklanov, K.M. Belotsky et al. (+1 more)

Numerous microlensing survey programs have constrained the possibility of dark matter existing in the form of compact objects within the Galactic halo. These constraints on the dark matter fraction were derived under the assumption of isolated, widely separated objects. This work investigates microlensing by primordial black holes (PBHs) organized into clusters. In this scenario, it is necessary to account for both the influence of neighboring PBHs and the collective gravitational potential of the entire cluster, which significantly complicates the microlensing light curve. Events exhibiting such complex light curves elude detection in observational experiments such as MACHO, EROS, OGLE, POINT-AGAPE, and HSC. It is demonstrated that a significant fraction of PBH dark matter (up to 93% for the models studied) remains undetected in these observational data. However, for all considered cluster models, a substantial population of PBHs still behaves as isolated lenses. Consequently, the clustering of PBHs does not completely eliminate the microlensing constraints on the PBH contribution to dark matter.
Read the full abstract

All authors K.A. Toshchenko, P.V. Baklanov, K.M. Belotsky, S.I. Blinnikov

High Energy Astrophysical Phenomena

Investigating peculiar prompt emission properties of the multi-Peaked GRB 250129A

Ankur Ghosh, Tamador K. M. Aldowma and Soebur Razzaque

We present a high-energy spectral analysis of GRB 250129A, which was triggered by the Swift-BAT. The burst exhibits a complex, multi-peaked temporal structure characterized by two distinct emission episodes, with the main peak occurring approximately 180 seconds after the BAT trigger. The time-integrated spectral analysis in the 15 - 150 keV energy range indicates that a broken power-law (BPL) model provides the best fit, signifying a non thermal origin of the prompt emission. A time resolved spectral analysis, performed using the Bayesian block technique, shows that the intervals around the main emission peak are well described by the BPL model, while the fits for low count intervals remain less constrained. An evident intensity tracking behavior is observed between the flux and the spectral peak energy (Ep). Furthermore, both the Amati relation and hardness - intensity correlation suggest that GRB~250129A occupies an intermediate regime, acting as a bridge between long and ultra long GRBs.
Read the full abstract

Cosmology and Nongalactic Astrophysics

An Opacity-Free Test of the Cosmic Distance Duality Relation Using Strongly Lensed Gravitational Wave Signals with Space-Based Detector Networks

Yong Yuan, Minghui Du, Benyang Zhu et al. (+4 more)

The cosmic distance duality relation (CDDR), expressed as $d_L(z) = (1+z)^2 D_A(z)$, is a fundamental relation in modern cosmology. In this work, we apply a method to test the CDDR using simulated strongly lensed gravitational-wave (SLGW) signals from massive binary black holes (MBBH) as observed by proposed space-based detector networks. Our analysis is conducted under the point-mass lens model, considering the strong lensing scenario that produces two images. We generate 90 days of simulated SLGW data for 10 events based on the Population III stellar formation model, with source redshifts in the range $z_s \in [2,6]$ and lens redshifts in $z_L \in [0.2,1]$. The deviation of CDDR is parameterized by $\eta_1(z) = 1 + \eta_0 z$ and $\eta_2(z) = 1 + \eta_0 z/(1+z)$, and we incorporate the deviation parameter $\eta_0$ directly into the waveform model. Parameter estimation is performed within a Bayesian statistical framework, combining simulated data from both Taiji and LISA. For a single lensed event, the joint Taiji+LISA analysis improves the measurement precision of $\eta_0$ by roughly a factor of two compared with Taiji-only observations. By combining 10 simulated events, the population-level constraints on $\eta_0$, quantified by the half width of the $95\%$ credible interval, reach approximately $2.61\times10^{-4}$ ($1.72\times10^{-4}$) for the $\eta_1(z)$ parameterization and $1.22\times10^{-3}$ ($6.86\times10^{-4}$) for $\eta_2(z)$ in the Taiji-only (Taiji+LISA) scenario, respectively. The inferred values of $\eta_0$ remain consistent with $\eta_0 = 0$ within the estimated uncertainties, with no statistically significant evidence for deviations from the CDDR at the achieved precision. These results demonstrate the significant advantage of joint space-based observations for high-precision tests of the CDDR.
Read the full abstract

All authors Yong Yuan, Minghui Du, Benyang Zhu, Xin-yi Lin, Wen-Fan Feng, Peng Xu, Xilong Fan

High Energy Astrophysical Phenomena

Photometric and late-time spectropolarimetric observations of GRB 250129A afterglow

A. Ghosh, S. Razzaque, J. Barnard et al. (+15 more)

Gamma-Ray Burst (GRB) afterglows arise from the interaction of relativistic ejecta with the circumburst medium and are observed across the electromagnetic spectrum. Afterglow polarisation is expected at early and late phases depending on the presence of reverse shocks (RS) and the observer's viewing geometry relative to the jet. Polarimetric observations of GRB afterglows provide a unique diagnostic tool to probe the geometry and structure of magnetic fields in the emitting region, which cannot be inferred from photometric or spectroscopic data alone. We report late-time (~19 hours post-burst) spectropolarimetric observations of GRB 250129A using the Southern African Large Telescope (SALT). The data reveal a hint of linear polarisation, with no evidence for rotation in the polarisation angle across wavelengths. Polarisation is typically expected during the early afterglow (<100 s) when the RS dominates. However, multi-wavelength modelling shows no indication of RS contribution at late times. Modelling incorporating both forward shock (FS) and RS components confirms that the RS fades rapidly after ~100 s. The afterglow emission is best explained by an off-axis viewing geometry of a jet with a Gaussian core and wings evolving in a uniform density environment. GRB 250129A thus provides rare observational evidence linking late-time polarisation to jet geometry and structure.
Read the full abstract

All authors A. Ghosh, S. Razzaque, J. Barnard, J. C. Joshi, R. Gupta, D. A. H. Buckley, B. van Soelen, N. Dukiya, A. Gupta, A. S. Moskvitin, J. Cooper, S. Chandra, K. M. Jayasurya, K. Misra, N. Rawat, L. Resmi, O. I. Spiridonova, R. I. Uklein

Instrumentation and Methods for Astrophysics

Data-calibrated point spread function prediction: General description of the method and demonstration on MUSE-NFM

Arseniy Kuznetsov, Benoit Neichel, Sylvain Oberti et al. (+1 more)

Precise knowledge of the point spread function (PSF) underpins many data analysis steps in astronomy, from photometry and astrometry to source de-blending and deconvolution. In adaptive optics (AO) observations, however, the PSF is highly variable with wavelength, field position, and observing conditions, making it difficult to model. Traditional PSF reconstruction (PSF-R) requires full AO telemetry and complex infrastructures, limiting its routine use, especially for tomographic systems. We present a practical framework for fast, accurate, and data-calibrated PSF modeling that captures the spatial and spectral variability of AO-corrected PSFs without relying on complete AO telemetry. Our approach builds on a Fourier-based PSF model inspired by astro-TIPTOP. As inputs, our model uses only a compact set of physically meaningful parameters retrievable from the ESO archive. A lightweight neural network corrects these inputs to achieve the best match with real data. It is trained end to end with the PSF model, allowing it to learn any miscalibrations directly from on-sky data. The framework achieves high accuracy on on-sky data. On a test set of MUSE-NFM standard stars, it yields median errors of 13.5% in the Strehl ratio and 10.9% in the core full width at half maximum (FWHM). In crowded MUSE-NFM observations of $\omega$ Centauri, the method predicts dozens of off-axis, wavelength-dependent PSFs with a Strehl error of <5% and a FWHM error of 4.6%, enabling source separation without per-star PSF extraction. Our compact, physics-informed, and data-calibrated model delivers accurate, polychromatic, and field-varying PSFs without relying on full AO telemetry. While demonstrated on MUSE-NFM, the method is still transferable to other AO-assisted instruments.
Read the full abstract

All authors Arseniy Kuznetsov, Benoit Neichel, Sylvain Oberti, Thierry Fusco

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Dark Matter Detection through Rydberg Atom Transducer

J. F. Chen, Haokun Fu, Christina Gao et al. (+5 more)

Ultralight bosonic dark matter with masses in the meV range, corresponding to terahertz (THz) Compton frequencies, remains largely unexplored due to the difficulty of achieving both efficient signal conversion and single-photon-sensitive detection at THz frequencies. We propose a hybrid detection architecture that integrates a dielectric haloscope, Rydberg-atom transducer, and superconducting nanowire single-photon detection within a unified cryogenic platform operating at $\lesssim 1\,\text{K}$. The dielectric haloscope converts dark matter into THz photons via phase-matched resonant enhancement, achieving form factors $C \sim 0.4$ and loaded quality factors $Q_L \sim 10^4$. A cold $^{87}$Rb ensemble then coherently up-converts the THz signal to the optical domain through six-wave mixing among Rydberg states. The intrinsic directionality and narrow bandwidth ($\Delta\nu_{\mathrm{atomic}} \sim 1\,\text{MHz}$) of this process provide extra suppression of isotropic thermal backgrounds. With 10 days of integration at $0.3\,\text{K}$, we project sensitivity to the axion-photon coupling $g_{a\gamma\gamma} \sim 10^{-13}\,\mathrm{GeV}^{-1}$ at $m_a \sim 0.4\,\text{meV}$, reaching the QCD axion band and opening the THz window for searches of both axion and dark photon dark matter.
Read the full abstract

All authors J. F. Chen, Haokun Fu, Christina Gao, Jing Shu, Geng-Bo Wu, Peiran Yin, Yi-Ming Zhong, Ying Zuo

Solar and Stellar Astrophysics

Signatures of localised particle acceleration at a global coronal shock wave

C. Cuddy, D. M. Long, M. Nedal et al. (+2 more)

Extreme ultraviolet (EUV) waves are global waves in the solar corona which can accelerate particles. The efficiency of the acceleration depends on local plasma characteristics e.g. Alfvén speed and the geometry of the magnetic field. This shock-driven particle acceleration can produce radio signatures such as Type II radio bursts and herringbone emission. Here we investigate signatures of particle acceleration by a weak coronal shock on 10 March 2024. In particular, we combine EUV images with radio imaging and spectral observations to determine how and where this weak shock could accelerate energetic particles. A potential field source surface extrapolation was used to examine the pre-eruption ambient magnetic field while the evolution of the global wave was probed using running difference and base difference EUV images. The EUV images enabled the speed and Alfvén Mach number of the EUV wave to be characterised. The combination of radio images and dynamic spectra provide evidence of beams of shock-accelerated electrons localised to a dimming region at the time the EUV wave passes through it. The speeds and energies of these electrons were estimated from the drift rates of their herringbones. The EUV wave initially propagated West, channelled by loop systems, before changing direction northward. From the EUV intensity jump at the wavefront, the Alfvén Mach number was estimated to be approximately 1.005 at the time that the herringbones were produced. The herringbone drift rates revealed accelerated electron energies of 75-122 keV, using Newkirk density models with scaling factors of 1.3-2.6. These observations suggest that the weak lateral shock impacted quasi-perpendicular open field in a dimming region, enabling localised particle acceleration. This indicates that the geometry of the ambient magnetic field relative to the shock strongly governs where particles can be accelerated.
Read the full abstract

All authors C. Cuddy, D. M. Long, M. Nedal, S. Bhunia, P. T. Gallagher

Astrophysics of Galaxies

3D CMZ V: A new orbital model of our Galaxy's Center, informed by data across the electromagnetic spectrum

Dani R. Lipman, Cara Battersby, Daniel Walker et al. (+9 more)

The 3D structure of The Milky Way's Central Molecular Zone (CMZ) informs our understanding of star formation cycles, black hole accretion, and the evolution of galactic nuclei. However, a comprehensive 3D model has remained elusive, as no singular dataset nor theory contains the requisite information to describe the orbital motion of the gas. We implement a Bayesian framework to flexibly combine datasets across the electromagnetic spectrum for molecular clouds in our CMZ catalog. We develop near/far metrics for each dataset, including dust extinction, absorption, stellar densities, X-ray echoes, and proper motions; and report a posterior positional probability density function (PPDF) for each cloud. We then use the posterior PPDF distributions for all CMZ clouds to search for a best fitting x$_2$ orbit. We find that no single orbit is a perfect fit, but the structure can overall be represented by nested x$_2$ orbits, with major axes ranging from about $72 < a < 146$ pc. We also present projected line of sight distance estimates for all 31 clouds in the catalog. Our results highlight asymmetries along the line of sight, with most clouds lying on the near side of the Galactic Center, and agree overall with current near/far assumptions for most CMZ clouds, including those in the Sgr A region, which may be much closer to the center. We conclude that the CMZ can be well-described by x$_2$ orbital families, and that the overall gas distribution is more complex than a single closed or open elliptical orbit.
Read the full abstract

All authors Dani R. Lipman, Cara Battersby, Daniel Walker, Maïca Clavel, B.L. DuBois, Adam Ginsburg, Jonathan D. Henshaw, Ralf S. Klessen, Elisabeth A.C. Mills, Francisco Nogueras-Lara, Mattia C. Sormani, Robin G. Tress

Solar and Stellar Astrophysics

Formation of spirals in early stage protoplanetary discs

Marc Van den Bossche and Oliver Gressel

Class II protoplanetary discs feature numerous non-axisymmetric substructures like spirals and the underlying mechanisms for their formation are still highly debated. Coincidentally, early stage, massive discs are subject to the gravitational instability that causes them to collapse into denser substructures. However, like for most instabilities, real systems usually remain marginally stable, here with Toomre parameter $Q \gtrsim 1$. We study how the self-gravity of the gas triggers the growth of spiral structures in the disc. We specifically focus on discs that are considered stable, that is, with respect to the gravitational instability (with $Q > 1$), as these discs remain unstable to non-axisymmetric perturbations like spirals. After a linear stability analysis, we produce high-resolution 2D shearing sheet simulations with the GPU-accelerated code \idefix of self-gravitating discs. We probe different initial densities and thermodynamical models of Toomre-stable discs. The initial transient growth of the spiral wave matches the linear theory provided we take into account the time dependency of the amplification. The spirals are then rapidly non-linearly amplified with growth rate $\approx 10$ orbital time scale. After this time spiral large scale mode are amplified up to 1000 times more than linear theory predicts. At later times, low density discs reach a weak gravito-turbulent state with $\alpha\approx 10^{-3}$ and discs with higher density undergo runaway collapse of the spiral arms. All discs exhibit dominant large-scale spirals.
Read the full abstract

Solar and Stellar Astrophysics

The Lyman-alpha Emission in Solar Flares. II. A Statistical Study on Its Relationship with the White-light plus Soft X-ray Emission

De-Chao Song, Ying Li, Qiao Li et al. (+1 more)

The hydrogen \lya\ line and the white-light (WL) continuum are two key diagnostics of energy transport in the lower atmosphere during solar flares, yet their relationship remains poorly understood. Here we present a statistical analysis of 69 white-light flares (WLFs) to investigate the relationships among the \lya, soft X-ray (SXR), and WL continuum emissions using the data from GOES and the Helioseismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory. We find that the \lya\ contrast in these WLFs ranges 0.8–28.5% with a mean value of 7.0%. Positive power-law relationships exist among peak enhancements in SXR, \lya, and WL. For most events, the \lya\ peak is nearly co-temporal with the peak of SXR time derivative, whereas the WL peak is either co-temporal with or lags those of \lya\ and SXR derivative. The \lya\ and WL rise times are similar ($\sim$3–4 min) and correlated. We also find that the radiated energy in \lya\ and HMI narrow-band WL has a positive power-law relationship with duration. In particular, the power-law index for the narrow-band WL is very close to 1/3 as predicted by magnetic reconnection theory. On average, the radiated energies in GOES \lya\ and SXR bands are approximately three orders of magnitude greater than the energy emitted in the continuum near 6173 Å\ with a bandwidth of 1 Å. Our findings provide new constraints on lower-atmosphere energy transport in solar flares and can serve as valuable references for modelling and interpreting the flares on solar-type stars.
Read the full abstract

All authors De-Chao Song, Ying Li, Qiao Li, Xiaofeng Liu

Solar and Stellar Astrophysics

Upper Limits to Long-Term Variability of Solar-Type Stars from Observations of the Open Cluster M67

Steven R. Spangler

Variations in the luminosity of the Sun on timescales of thousands to millions of years could potentially be responsible for terrestrial climate variations in the Phanerozooic geological period (last 540 million years). In this paper, I consider a method that utilizes observations of an open star cluster with approximately the age of the Sun, specifically M67, with data taken from Geller (2015). The idea is to measure the width of the main sequence in the region of solar-type stars, here generously defined to be about spectral class G0 - K1. This width gives an estimate of the dispersion in absolute magnitude of nominally solar-type stars. The sample used consists of 170 solar-type main sequence stars which are not known to be binaries. With this sample, I form an empirical measurement of the width of the main sequence, which is compared with a theoretical expression from Spangler (2025). The measured spread is consistent with a value of sigma, the normalized Gaussian variability of the primary and (if present) secondary star, in the range of 0.100 - 0.135. However, the expected value of sigma from purely photometric noise is estimated as 0.101. Generous upper limits to the intrinsic variability contribution to the inferred width are sigma in the range 0.058-0.089. These limits are not totally devoid of interest in a paleoclimatic context. However, major improvements in the technique are possible with the use of existing data sets from space astronomy missions such as Gaia and Kepler.
Read the full abstract

Cosmology and Nongalactic Astrophysics

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

Mohit Yadav, Pankaj Chavan and Tapomoy Guha Sarkar

Pure cosmographic reconstructions provide a model-agnostic approach towards constraining cosmic evolution.We develop a hybrid phenomenological semi-cosmographic framework that adopts a Padé-rational fraction parametrization of the luminosity distance, and incorporates a phenomenologically motivated two-body decaying dark matter (DDM) sector. In this framework, the effective residual dark-energy equation of state is obtained by subtracting the contributions of DDM, baryons, and radiation from the Padé approximation to the expansion history under the assumptions of flat FLRW cosmology and the standard Friedmann equations. Cosmological observables computed from these semi-cosmographic quantities are fitted to the data to reconstruct cosmological evolution. We present two separate analyses. First, we perform a data fit to the current DESI DR2 BAO measurements. Second, we report a forecast based on an optimistic mock 21\,cm intensity-mapping power spectrum at $z\simeq 1.75$ for a SKA1-MID-like instrument, used jointly with the DESI DR2 BAO data. The forecast serves as proof of concept to demonstrate a possible improvement in model constraints when measurements of the 21\,cm power spectrum complement the BAO background observables. Posterior constraints on the Padé and DDM parameters are obtained using Markov Chain Monte Carlo (MCMC) analysis, and we present the corresponding equations of state of the massive daughter and the effective reconstructed dark energy.
Read the full abstract

Earth and Planetary Astrophysics

Perihelion observations of interstellar comet 3I/ATLAS with the IRAM 30-m telescope

N. Biver, D. Bockelée-Morvan, R. Moreno et al. (+6 more)

3I/ATLAS is the third interstellar comet identified as passing through the Solar System. Its high outgassing activity and favourable perihelion passage on October 29, 2025 UT provided an excellent opportunity to investigate the composition of its coma gases through millimeter spectroscopy. We present observations undertaken with the IRAM 30-m telescope on November 1–3, 2025 at an heliocentric distance of 1.36–1.37 au. Lines of HCN, CH$_3$OH, CO, and H$_2$CO are well detected, and $\sim$4$\sigma$ detections are obtained for CS and CH$_3$CN. The search for H$_2$S was unsuccessful. Abundances of CO, H$_2$CO, CH$_3$OH, and CH$_3$CN relative to HCN are in the upper ranges of values measured in Solar System comets. The sulfur-to-carbon abundance ratio in 3I/ATLAS's coma is at most the minimum value observed in comets. The unusually low expansion velocity of coma gases suggests a near-nucleus gas flow driven by heavy molecules such as CO$_2$, and/or a large fraction of the gaseous production coming from subliming icy grains.
Read the full abstract

All authors N. Biver, D. Bockelée-Morvan, R. Moreno, J. Crovisier, G. Paubert, V. Zakharov, J. Boissier, M. A. Cordiner, N. X. Roth