Peer Revenue Publish or profit

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Thursday 5 March

Astrophysics of Galaxies

The Age of the Universe with Globular Clusters IV: Multiple Stellar Populations

We revisit the determination of the age of the Universe from galactic globular clusters, extending previous analyses by explicitly accounting for the presence of multiple stellar populations within each cluster. Using high–quality Hubble Space Telescope color–magnitude diagrams for 69 globular clusters, we relax the standard single–population assumption, and model two stellar populations with independent ages, metallicities, helium abundances, and population fractions. The inference is performed using the full color–magnitude diagram morphology, an explicit treatment of field contamination, and a hierarchical framework that propagates non–Gaussian age posteriors. Allowing for multiple stellar populations has a negligible impact on globular cluster age estimates. The ages of the oldest populations remain fully consistent with those obtained under the single–population assumption, with differences at the $0.6\sigma$ level. Restricting to the metal–poor subsample ([Fe/H] $< -1.5$), we infer a dominant old component with mean age $t_{\rm GC}=13.61\pm0.25\,\mathrm{(stat)}\,\pm0.23 \mathrm{(sys)}\,\mathrm{Gyr}$. Adopting a conservative delay between the Big Bang and the formation of the first globular clusters, we obtain an age of the Universe of $t_{\rm U}=13.81\pm0.25\,\mathrm{(stat)}\,\pm0.23 \mathrm{(sys)}\,\mathrm{Gyr}$. In addition to age constraints, our analysis yields simultaneous measurements of metallicity and helium content for the different populations, including constraints on helium enrichment and population fractions which are consistent with independent determinations from the literature. These results demonstrate that globular–cluster–based cosmic chronometry is robust to stellar population complexity, reinforcing its role as a precise and largely cosmological model–independent probe of the age of the Universe.
Read the full abstract

All authors David Valcin, Raul Jimenez, Carmela Lardo, Uroš Seljak, Licia Verde

Cosmology and Nongalactic Astrophysics

STOchastic LAttice Simulation of hybrid inflation

We investigate the spatial profile of the curvature perturbation generated in multi-waterfall hybrid inflation models, which are known to produce various topological defects. Using the lattice simulation code STOchastic LAttice Simulation, based on the stochastic formalism of inflation, we analyse six cases by varying the number of waterfall fields $n$ and the functional form of the inflaton potential (“Quadratic” and “Cubic” cases). Our statistical analysis shows that the probability density functions (PDFs) and power spectra are broadly consistent with the so-called stochastic-$\delta N$ algorithm. The “Cubic” case also exhibits a characteristic upper bound in the PDF, as discovered in our previous work, that suppresses \acl{PBH} formation while potentially affecting halo formation. Furthermore, we employ the Euler characteristic as a topological diagnostic tool to identify the structures of the waterfall fields as well as the curvature perturbation. We find that the topological defects, such as domain walls ($n=1$), cosmic strings ($n=2$), and monopoles ($n=3$), are reconnected during inflation into finer structures by the stochastic noise, making their correlation lengths much smaller than the Hubble scale at the critical point of the waterfall phase transition counterintuitively. The Euler characteristic also implies global structures of the curvature perturbation for $n=1$, though we do not conclude if they are due to the domain wall, because neither the strings ($n=2$) nor monopoles ($n=3$) leave such structures. The global structures of the curvature perturbation will provide a novel probe for the physics of the early universe.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Probing Dark Energy on the Moon

The effective field theory (EFT) of cosmic acceleration provides a model-independent framework for describing dark energy and modified gravity, yet many of its defining operators remain weakly constrained by existing observations. We show that measurements of horizon-scale metric fluctuations with a lunar laser interferometer can directly probe the kinetic sector of the EFT of dark energy, enabling constraints on operators governing scalar perturbation dynamics rather than only the background expansion history. In particular, we demonstrate sensitivity to the EFT kinetic coefficient $M_2^4$ and the associated sound speed of dark energy, $c_s^2$. This establishes a qualitatively new observational handle on the microphysical consistency conditions of late-time acceleration models, allowing broad regions of EFT parameter space to be probed, constrained, or potentially discovered.
Read the full abstract

Instrumentation and Methods for Astrophysics

A New Method for Identifying Contaminating Sources and Locating Target Sources through the Cross-Arm Features of Micro Pore Optics

The Pathfinder of the Type-A satellites in the Chasing All Transients Constellation Hunters (CATCH) space mission is equipped with Micro-Pore Optics (MPOs) and four single-pixel Silicon Drift Detectors (SDDs). Due to the lack of position resolution in an individual SDD, we propose a new method based on the cross-arms in the point spread function (PSF) of MPOs to enhance the satellite's capability in identifying contaminating sources and locating target sources. By placing one detector on each of the horizontal and vertical cross-arms on the focal plane, we can use the changes in the relative counts on the cross-arms detectors to deduce the location of the source. Simulated observations demonstrate that, for a target source with a flux of 1 Crab and an exposure time of 200 s, the cross-arms detectors can identify contaminating source with the same flux level at an off-axis angle larger than 8', and improve positioning accuracy to 6'. Furthermore, we extend the simulation study to CATCH Type-A, which plans to use an SDD array. In situations where sources exhibit the same flux of 1 Crab and the exposure time is merely 1 s, a 16x16 SDD array is capable of identifying contaminating source with an off-axis angle greater than 2.4' and can achieve a positioning precision of 1.8'.
Read the full abstract

All authors Yiming Huang, Lian Tao, Jin-Yuan Liao, Shuang-Nan Zhang, Stéphane Schanne, Bertrand Cordier, Shaolin Xiong, Juan Zhang, Zhengwei Li, Qian-Qing Yin, Xiangyang Wen, Sheng Yang, Min Gao, Donghua Zhao, Xiang Ma, Yue Huang, Liang Zhang, Liming Song

High Energy Astrophysical Phenomena

Resolving diffusion signatures in distant pulsar halos with current and future experiments

Pulsar halos provide a unique probe of cosmic-ray propagation in the vicinity of pulsars and have important implications for our understanding of particle diffusion in the interstellar medium. However, the number of firmly identified pulsar halos remains limited. One of the main challenges is the difficulty in unambiguously confirming halo candidates through precise morphological measurements with current $\gamma$-ray observations. In this work, we investigate the prospects for identifying pulsar halo candidates through morphological discrimination using simulations of two advanced $\gamma$-ray experiments: LHAASO-KM2A and the Cherenkov Telescope Array (CTA). Using mock observations with realistic instrumental responses, we assess the ability of each experiment to distinguish diffusion-based halo morphologies from alternative simplified spatial models. Our results show that both increased photon statistics and improved angular resolution significantly enhance the power of morphological discrimination. In particular, CTA benefits from its superior angular resolution, while LHAASO-KM2A gains sensitivity from its large effective area at the highest energies. These results indicate that future $\gamma$-ray observations have the potential to expand the sample of pulsar halos and provide further insights into cosmic-ray transport around pulsars.
Read the full abstract

All authors Yong-Jian Wei, En-Sheng Chen, Kun Fang, Xiao-Jun Bi

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

Search for continuous gravitational waves from neutron stars in five globular clusters in the first part of the fourth LIGO-Virgo-KAGRA observing run

We present the results of directed searches for continuous gravitational waves from unknown neutron stars in five Milky Way globular clusters. We carry out these searches in the LIGO data from the first eight months of the fourth LIGO-Virgo-KAGRA observing run using the WEAVE semicoherent program, which sums matched-filter detection-statistic values over many time segments spanning the observation period. No gravitational wave signal is detected in the search band of 20-475 Hz. Injections of simulated continuous wave signals in the data indicate that we achieve the most sensitive results to date across most of the explored parameter space volume, obtaining median 95% confidence level upper limits as low as $\sim 4.2 \times 10^{-26}$ near 282 Hz for NGC 6397. We also derive upper limits on neutron star ellipticity and $r$-mode amplitudes, reaching $\lesssim 10^{-5}$ and $\lesssim 10^{-3}$, respectively, at frequencies above 200 Hz.
Read the full abstract

All authors Damon H. T. Cheung, Keith Riles, Rafel Amengual, Preet Baxi, Alicia Calafat, Anamaria Effler, Tabata Aira Ferreira, Evan Goetz, David Keitel, Tom Kimpson, Alan M. Knee, Joan-Rene Merou, Quynh Lan Nguyen, Joseph O'Leary, Ornella J. Piccinni, Alicia M. Sintes, Rebekah Slocum, Karl Wette

Astrophysics of Galaxies

Accelerating massive galaxy formation with primordial black hole seed nuclei

If massive primordial black holes (PBHs) exist and constitute a fraction of the dark matter, they can dramatically catalyze galaxy formation. By acting as pre-existing, high-density seeds, they can shorten the galaxy assembly time to as little as 100 Myr for up to 10^8 solar mass PBH seeds, allowing for the rapid formation of host halos. Furthermore, low surface brightness or diffuse galaxies may represent a natural outcome of this process, perhaps as the residue of halos seeded by smaller PBHs that failed to accrete a major baryonic component.
Read the full abstract

High Energy Astrophysical Phenomena

Could the interaction of jet and SN ejecta be the cause of X-ray knots observed in a radio galaxy?

We investigate the interaction between relativistic jets and supernova (SN) ejecta as a potential origin of X-ray knots in radio galaxies, employing knot A in M 87 as a test case. By modeling the dynamical evolution of the interaction, we evaluate this scenario based on particle acceleration efficiency and spatial morphology. Our modeling indicates that the ejecta shock expands to only ~ 30 pc, which is inconsistent with the observed spatial scale of knot A (~ 60 pc). In contrast, the jet shock can successfully reproduce the observed scale after approximately 3000 yr, with the ejecta being accelerated to a bulk velocity of \beta~ 0.43. We fit the multi-wavelength spectral energy distribution (SED) using a one-zone leptonic framework, attributing the X-rays to synchrotron radiation from electrons accelerated up to~1 PeV at the jet shock. The derived magnetic field is approximately 70 uG in the SN ejecta rest frame, which is significantly below the equipartition value. Protons may be accelerated up to ~ EeV, supporting the hypothesis that the jets of radio galaxies (RGs) may be the potential site for ultra-high-energy cosmic-ray (UHECR) acceleration within the framework of the jet-ejecta interaction.
Read the full abstract

All authors Jia-Chun He, Xiao-Na Sun, Hao-Qiang Zhang, Yun-Feng Liang, Hai-Ming Zhang, Da-Bin Lin, En-Wei Liang

Astrophysics of Galaxies

Early Planet Formation in Embedded Disks (eDisk). XVIII. Indication of a possible spiral structure in the dust-continuum emission of the protostellar disk around IRAS 16544-1604 in CB 68

We performed numerical simulations along with radiative transfer calculations to reproduce an intriguing asymmetric shoulder feature in the dust-continuum emission of the protostellar disk around one of the eDisk targets, the Class 0 protostar IRAS 16544-1604 in CB 68. This is our first attempt to bridge the theoretical works of protostellar disk evolution and the eDisk observations. We found that while our hydrodynamic simulations form spiral structures caused by gravitational instability, they become less discernible after the disk is inclined and convolved with the telescope beam. The widths of the spiral structure as obtained by our numerical simulations are ~0.1-0.8 times the eDisk beam size of 4.5 au. Our modeling effor implies that the apparent absence of spiral features in the eDisk observations does not necessarily indicate the real absence of internal substructures and gravitational instability. We also found that the asymmetric shoulder structure of the continuum profile along the major axis appears when the disk is massive enough with a Toomre parameter Q~1. This mechanism offers a potential explanation for the observed, asymmetric shoulder features in the disks surrounding IRAS 16544-1604 and the other eDisk sources.
Read the full abstract

All authors Sanemichi Z. Takahashi, Shigehisa Takakuwa, Ryosuke Nakanishi, Yusuke Tsukamoto, Kazuya Saigo, Miyu Kido, Nagayoshi Ohashi, Zhi-Yun Li, Leslie W. Looney, Zhe-Yu Daniel Lin, Mayank Narang, Kengo Tomida, John J. Tobin, Jes K. Jørgensen

Instrumentation and Methods for Astrophysics

AstroInspect: a web-based system to organize, assess, and visually inspect astronomical objects

The rapid growth of imaging and spectroscopic surveys has intensified the need for efficient tools that support visual inspection, a practice that remains essential for tasks such as classification, catalog refinement, and validation of automated methods. Existing solutions, however, often require the use of multiple platforms and complex workflows to integrate heterogeneous data. To address this challenge, we present the first release of the AstroInspect (https://astroinspect.github.io), a web-based system which ensures seamless access to several astronomical resources. The system provides an intuitive graphical user interface (GUI) through which users can upload catalogs of objects defined by celestial coordinates. AstroInspect automatically enriches these catalogs with complementary information, including imaging, spectroscopic, and photometric data retrieved in real time from surveys such as the Sloan Digital Sky Survey (SDSS), the Legacy Surveys (LS), and the Southern Photometric Local Universe Survey (S-PLUS). As an example of its scientific utility, we used AstroInspect to identify H$\alpha$ emission-line galaxies within a 7 deg radius in the direction of the Hydra I cluster (also known as Abell 1060) by visual inspection. Using a candidate set of 981 galaxies selected from S-PLUS photometric data, we produced a catalog of 80 galaxies with confirmed H$\alpha$ emission. These results highlight the potential of AstroInspect to support efficient visual inspection workflows.
Read the full abstract

All authors Natanael M. Cardoso, Claudia Mendes de Oliveira, Angela C. Krabbe, Analia V. Smith Castelli, Gustavo B. Oliveira Schwarz, Lilianne Nakazono, Ricardo Demarco, Maiara S. Carvalho, William Schoenell, Tiago Ribeiro, Antonio Kanaan, Antonio M. Saraiva

Astrophysics of Galaxies

3D Rotation of the Open Cluster NGC 2516

We have combined Gaia astrometry with Gaia-ESO Survey radial velocities to measure the 3D rotation of the open cluster NGC 2516. We compiled a sample of 376 members with astrometry and spectroscopy and use these to determine a distance to the cluster of 406.4 +/- 0.8 pc, which we then use to infer the 3D positions and velocities of all stars in the cluster using a Bayesian model. We identify the axis of maximum cluster rotation and measure a median rotational velocity of 0.11 +/- 0.02 km/s. We find the axis of maximum cluster rotation to be 52 +/- 22 degrees to the plane of our galaxy. We compare this rotation rate to measurements of cluster rotation in other open clusters and find that it is inconsistent with the expected dependences on cluster age and mass.
Read the full abstract

Wednesday 4 March

High Energy Astrophysical Phenomena

X-ray Doppler tomography of Fe K$\alpha$ emission in a low-mass X-ray binary 4U 1822-371 - a localized reflector at the accretion stream-disk overflow

We present the X-ray Doppler tomography of the Fe K$\alpha$ (6.4 keV) fluorescence line of the low-mass X-ray binary 4U 1822-371 obtained with XRISM. Eleven orbits of this short period (5.57 hr) binary were covered. The Doppler shift of the line shows clear modulation with the orbital period, motivating us to apply the Doppler tomography in the X-ray band for the first time. The resulting velocity map reveals a compact feature at ($v_{\mathrm{x}}$, $v_{\mathrm{y}}$) $\sim$ ($-$550, $+$125) km s$^{-1}$. This is inconsistent with the emission originating from a symmetric accretion disk, an extended corona around the neutron star, or the surface of the neutron or companion star. Instead, it suggests that the emission originates from the accretion stream-disk overflow. Remarkably, the Fe K$\alpha$ velocity map closely resembles that of the O VI 3811 Å, indicating that both X-ray and optical lines arise from the same site irradiated by the central X-ray source. These results provide the first velocity-resolved X-ray map of the fluorescent line, directly localizing the major reflector in an X-ray binary and establishing X-ray Doppler tomography as a new probe of the structures of accreting systems.
Read the full abstract

Instrumentation and Methods for Astrophysics

PulSKASim: A Pulsar Simulator for SKA-Scale Interferometric Observations

Accurate simulation of pulsar flux variability is critical for testing Square Kilometre Array (SKA) interferometric pipelines. However, most existing simulators neglect the effects of integration time and related observational parameters, limiting their realism and utility for interferometric end-to-end testing. To address these shortcomings, we develop a Pulsar Simulator for SKA-scale interferometric observations (PulSKASim), which models pulsar flux evolution across pulsar period, maximum flux, duty cycle, and noise, accounting for integration time, sampling, and observation duration, and naturally models flux smoothing that arises from finite integration within each dump time. PulSKASim generates synthetic measurement sets using functions in simulators for radio interferometers, such as OSKAR and Pyuvsim, where each snapshot contains pulsars with controlled flux levels, enabling realistic per-time-slot experiments. This simulator allows for detailed assessment of calibration, imaging, and detection pipelines under realistic SKA-like conditions, bridging pulsar variability modelling with interferometric simulation in a way not achievable by existing tools.
Read the full abstract

Instrumentation and Methods for Astrophysics

NEXUS: Quick Release Notes

NEXUS is a JWST Multi-Cycle (Cycles 3-5) GO Treasury imaging and spectroscopic survey around the North Ecliptic Pole during 2024-2028. It contains two overlapping tiers in depth and area coverage. The Wide tier ($\sim 400~{\rm arcmin}^2$) performs NIRCam/WFSS 2.4-5 $\mu$m grism spectroscopy with three annual epochs over 3 years (final spectral continuum ${\rm S/N/pixel>3}$ at F444W $<22.2$), accompanied by NIRCam multi-band imaging in F090W, F115W, F150W, F200W, F356W and F444W. The Deep tier ($\sim 50~{\rm arcmin}^2$) performs high-multiplexing NIRSpec 0.54-5.5 $\mu$m MOS/PRISM spectroscopy for ~10,000 targets in total, over 18 epochs with a 2-month cadence, along with F200W+F444W NIRCam imaging for each epoch. Parallel imaging observations with MIRI and additional NIRCam filters are also performed within the Wide and Deep tiers. The primary data covering the Deep tier (including NIRCam imaging, NIRSpec/MSA spectra, and vetted MSA spectroscopic redshifts) are released in regular Quick Data Releases to facilitate follow-up studies. This evolving document describes the MSA targeting information and observing status for each of the 18 Deep epochs, which started in May 2025 and continue on the regular 2-month cadence. We also describe the content and caveats of the quick release data and report selected cases of diverse scientific interests.
Read the full abstract

All authors Ming-Yang Zhuang, Yue Shen, Junyao Li, Zhiwei Pan, Lei Hu, Adam J. Burgasser, David A. Coulter, Jenny E. Greene, Feige Wang

2603.04586 · 4 Mar · This is an arXiv-only document that describes the quick data releases from the NEXUS-JWST program through 2028. This document will be updated regularly to follow the cadenced observations from NEXUS (with a delay of ~2 months after each epoch). Data products are available at https://ariel.astro.illinois.edu/nexus/qdr/. Please send comments and suggestions to M. Zhuang and Y. Shen

Astrophysics of Galaxies

The First Detection of Forbidden Emission Lines at the Outskirts of the AGN Broad Line Region?

Double-peaked (DP) broad emission line profiles in active galactic nuclei (AGNs) are often interpreted as signatures of rotating disk-like structures in the broad-line region (BLR) and are commonly observed in low-luminosity AGNs using recombination lines. We use optical spectroscopy to investigate the origin of double-peaked broad emission line profiles observed not only in hydrogen recombination lines but also in forbidden transitions in the LINER galaxy IC 1459. We detected DP emission in all strong optical lines except for the [S II] doublet, which has the lowest critical density among all the lines. We successfully fitted the DP broad profiles using a disk-like BLR model, assuming a circular accretion disk with an inclination of approximately 35 degrees and internal turbulence of about 500 km/s, confined within a maximum radius of 9.6 (+4.8, -1.1) light-years. We estimate a full width at half maximum of the DP profiles of about 3300 km/s. Our results provide new insights into the structure of the BLR, indicating that forbidden emission lines can be produced in lower-density regions near the outskirts of the BLR.
Read the full abstract

All authors Kelly F. Heckler, Rogemar A. Riffel, Daniel Marsango, Tiago V. Ricci, Angela C. Krabbe, Oli L. Dors

Solar and Stellar Astrophysics

The Age of the R127 & R128 Clusters: Implications for the LBV

We infer the age of the R127 and R128 clusters in the Large Magellanic Cloud (LMC) using Strömgren photometry from the literature and the age-dating algorithm, Stellar Ages. Analysis using single-star evolutionary models shows a substantial discrepancy between the relative numbers of bright blue stars and lower-mass stars as compared to expectations from a Salpeter mass function, and yields a younger age for the brightest blue stars than for the rest of the cluster. This inconsistency reflects an emerging trend among young clusters in the Local Group. In general, the resolution may be binary evolution or very rapid rotation, although in the specific case of the R127 and R128 clusters, unknown incompleteness in the data may also affect the relative numbers of low- and high-mass stars. The discrepancy grows toward fainter magnitudes, suggesting that the dataset is likely incomplete. However, when the five brightest stars are excluded, the observed and expected counts become consistent, demonstrating that the brightest stars are peculiar. These findings have direct implications for the luminous blue variable (LBV) R127, which is the only confirmed LBV in the LMC located within a young stellar cluster. LBVs have traditionally been considered products of single-star evolution, although there is growing recognition that binary interactions may play a critical role in their evolution. A more complete dataset, particularly deeper imaging with the Hubble Space Telescope, is needed to confirm whether the apparent absence of coeval stars arises solely from observational incompleteness or the broader trend of inconsistency in young cluster modeling.
Read the full abstract

All authors Mojgan Aghakhanloo, Jeremiah W. Murphy, Nathan Smith, Joseph Guzman

High Energy Astrophysical Phenomena

Repeating Nuclear Transients from Repeating Partial Tidal Disruption Events

Extragalactic nuclear transients that exhibit repeating outbursts can be modeled as the repeated dynamical interaction between bound stars and supermassive black holes (SMBHs). A subset of these transients, with recurrence timescales of months-to-years, have been explained as accretion flares from the repeated tidal stripping of a star by an SMBH, in a repeating partial tidal disruption event (rpTDE). We outline the scope of the rpTDE model and discuss hydrodynamical simulations and analytical predictions for the stability of stars undergoing repeated mass loss, and the long-term evolution of these flares as a function of stellar type and orbital parameters. Our findings demonstrate that high-mass and centrally concentrated stars undergo negligible changes in structure in response to small amounts ($\sim 1-10\% M_\star$) of mass loss, and can survive many mass-stripping encounters with an SMBH. Contrarily, low-mass and less evolved stars are unstable to mass loss, and would be destroyed within a few orbits. We discuss the implications of these results for constraining the stellar type and orbital parameters of observed sources, such as ASASSN-14ko, for which $\gtrsim 20$ flares have been observed, and AT2020vdq, which exhibits a second flare that is brighter than its primary outburst.
Read the full abstract

All authors Ananya Bandopadhyay, Eric R. Coughlin, Julia Fancher, C. J. Nixon, Dheeraj R. Pasham

Earth and Planetary Astrophysics

A Comparative Study of the Streaming Instability: Unstratified Models with Marginally Coupled Grains

The streaming instability is a leading mechanism for concentrating solids and initiating planetesimal formation in protoplanetary disks. Although numerous studies have explored its linear growth, nonlinear evolution, and implications for planet formation, the diversity of numerical methods and dust treatments used across the literature has made it difficult to assess which features of the instability are physically robust and which arise from code-dependent choices. We present the first systematic comparison of seven hydrodynamic codes–spanning finite-volume and finite-difference schemes and modeling dust either as Lagrangian particles or as a pressureless fluid–applied to the unstratified streaming instability with a dimensionless stopping time of unity. All codes reproduce the characteristic sequence of exponential growth, filament formation, and turbulent saturation, demonstrating broad agreement in the qualitative behavior of the instability. Quantitatively, however, the dust model remains the dominant source of variation at moderate resolution: particle-based simulations reach higher peak densities and exhibit broader high-density tails than fluid-based models at $512^2$ resolution, although increasing the number of particles brings their initial maximum density evolution into close agreement with that of dust-fluid models. At $1024^2$, these differences diminish substantially, indicating better agreement of the saturated-state statistics across dust treatments. In terms of computational performance, most particle implementations suffer from imbalanced parallelized loads, while execution on a GPU is at least two to three times more energy efficient and scales better at higher resolutions than on CPUs. Given the intrinsic stochasticity of this nonlinear system, only statistical diagnostics remain meaningful across codes.
Read the full abstract

All authors Stanley A. Baronett, Wladimir Lyra, Hossam Aly, Olivia Brouillette, Daniel Carrera, Victoria I. De Cun, Linn E. J. Eriksson, Mario Flock, Pinghui Huang, Leonardo Krapp, Geoffroy Lesur, Rixin Li, Shengtai Li, Jeonghoon Lim, Sijme-Jan Paardekooper, David G. Rea, Debanjan Sengupta, Jacob B. Simon, Prakruti Sudarshan, Orkan M. Umurhan, Chao-Chin Yang, Andrew N. Youdin

Solar and Stellar Astrophysics

Energetics and Emission in a Simulated Solar Flare Initialised by a Non-Force Free Magnetic Field

Solar flare simulations are commonly initialised using non-linear force free field (NLFF) extrapolations derived from photospheric vector magnetograms. However, the force free assumption neglects plasma forces and may limit the available free magnetic energy. In this work, we perform a controlled comparison of two three-dimensional resistive magnetohydrodynamic simulations of the X2.1-class flare that occurred on 2011 September 06 in NOAA Active Region 11283. The simulations differ only in their initial magnetic configuration: one is based on a conventional NLFF extrapolation, while the other employs a non-force free extrapolation. Both models are evolved in an identical stratified atmosphere using the same numerical framework, enabling direct assessment of how the initial magnetic assumptions influence flare dynamics and energetics. We find that the non-force free model undergoes more extensive magnetic restructuring and releases approximately twice as much magnetic energy ($\approx4.4 \times 10^{31}$ erg) as the NLFF case ($\approx2.3 \times 10^{31}$ erg), bringing the energy budget into closer agreement with expectations for X-class flares. Synthetic extreme ultraviolet emission in the 94A channel is computed for both simulations and compared with observations from the Solar Dynamics Observatory. The non-force free model produces a brighter and more spatially extended emission structure that more closely resembles the observed flare morphology and light curve. These results demonstrate that assumptions made in constructing the pre-flare coronal magnetic field can significantly affect flare energetics and observable signatures, and suggest that non-force free extrapolations provide a promising pathway toward more realistic data-constrained flare modelling.
Read the full abstract

All authors W. Bate, M. Gordovskyy, A. Prasad, A. S. Brun, A. Strugarek, M. V. Sieyra, P. Browning, S. Inoue, K. Matsumoto, A. Roddanavar

Solar and Stellar Astrophysics

EBLM XVII - Tidal Synchronization and Circularization in Tight Stellar Binaries

Tidal interactions in close stellar binaries are central to their orbital and rotational evolution, making observational tests of theoretical predictions essential for our understanding of the evolution of these, as well as close exoplanetary systems. Such tests require precise measurements of the orbital eccentricity and stellar rotation. The EBLM (Eclipsing Binary Low Mass) survey delivers a homogeneous sample of eclipsing binaries, composed of F/G/K primaries and M-dwarf (or low-mass K-dwarf) secondaries. We analyze 68 unequal mass binaries ($0.1 \leq q \leq 0.6$, where $q$ is the mass ratio), with measurable primary star rotation rates from TESS, and over a decade of radial velocity observations. This sample probes the critical regime where tidal effects are expected to transition between being efficient and inefficient. We find that ~75% of our sample has circularized, with eccentric systems confined to $P_{\rm orb} \gtrsim 3$ days, with modest eccentricities (e < 0.25). Roughly ~78% of our sample is synchronized, with nearly all binaries within a 3-day orbital period residing in a well-defined "synchronization zone". Beyond this, a minority of asynchronous systems persist, which cannot be easily explained by our application of current tidal mechanisms or by differential rotation.
Read the full abstract

All authors Ritika Sethi, David V. Martin, Adrian Barker, Pierre F. L. Maxted, Amaury H.M.J. Triaud, Vedad Kunovac, Wata Tubthong, Alison Duck, François Bouchy, Stéphane Udry

Cosmology and Nongalactic Astrophysics

Beyond21: A Global Framework for Cosmic Dawn and Reionization Within and Beyond the Standard Model

Observations of the Cosmic Dawn (CD) and Epoch of Reionization (EoR) are steadily improving, opening new opportunities to study early galaxies through complementary probes. To enable consistent interpretation of these observations, we present Beyond21, a fully open-source Python package that implements flexible prescriptions for Pop II and Pop III star formation and computes the resulting radiation backgrounds and their impact on the intergalactic medium. From this coupled evolution, Beyond21 predicts the global 21-cm signal, UV luminosity functions (UVLFs), the ionization history, and the contribution to the observed cosmic X-ray background (CXB) within a single, self-consistent pipeline. A full global evolution run executes in $\sim0.1 \ {\rm s}$ on a single CPU core, enabling broad, high-resolution parameter exploration. The modular architecture facilitates straightforward modification of astrophysical prescriptions and the incorporation of new physics. As an illustrative example, we implement a scenario in which a small fraction of dark matter is millicharged, leading to baryon cooling through elastic interactions.
Read the full abstract

Astrophysics of Galaxies

Aromatic Species in the Molecular Universe

Interstellar polycyclic aromatic hydrocarbon (PAHs) are an important component of the interstellar medium of galaxies, containing some 10 percent of the elemental carbon. Their vibrational emission dominates the mid-infrared spectra of galactic and extragalactic objects. PAHs control the heating of interstellar neutral gas and the charge balance of molecular clouds. PAHs are formed in the outflows from late type stars through chemical processes akin to those in sooting flames and then further processed in the interstellar medium by UV photolysis and strong shock waves. PAHs are also formed through ion molecule reactions and neutral radical reactions in dense cloud cores. The James Webb Space Telescope has provided a wealth of high-quality spectra that have provided new insights in the characteristics of the interstellar PAH family. Their analysis is supported by dedicated laboratory and quantum chemistry studies, feeding into detailed molecular physics models relevant to astronomical environments. Laboratory studies have also provided deeper insight in the chemical evolution of PAHs in the interstellar medium. This paper will review progress in the field and chart its future.
Read the full abstract

Instrumentation and Methods for Astrophysics

A Fast Generative Framework for High-dimensional Posterior Sampling: Application to CMB Delensing

We introduce a deep generative framework for high-dimensional Bayesian inference that enables efficient posterior sampling. As telescopes and simulations rapidly expand the volume and resolution of astrophysical data, fast simulation-based inference methods are increasingly needed to extract scientific insights. While diffusion-based approaches offer high-quality generative capabilities, they are hindered by slow sampling speeds. Our method performs posterior sampling an order of magnitude faster than a diffusion baseline. Applied to the problem of CMB delensing, it successfully recovers the unlensed CMB power spectrum from simulated observations. The model also remains robust to shifts in cosmological parameters, demonstrating its potential for out-of-distribution generalization and application to observational cosmological data.
Read the full abstract

High Energy Astrophysical Phenomena

Geometric masking in AGN jets and its implications for unification and blazar physics

I explore the theoretical implications of the Geometric Masking scenario for AGN Unification and the blazar sequence. In PKS 2155$-$304, GeV spectral hardening appears locked to its $\sim 1.7$ yr quasi-periodic oscillation trough, suggesting high-flux states are dominated by a soft, geometrically amplified envelope that masks an underlying hard core. While this foundational mechanism builds on observations of coexisting emission regions (Madero & Domínguez 2026) and phase-locked spectral hardening (Domínguez et al. 2026), here I propose hypothetical extensions that fall outside the scope of those works. I propose that Doppler boosting preferentially enhances this soft component in closely aligned jets. This visibility bias could extend AGN Unification and introduce a geometric modulation layer to the blazar sequence. Under this speculative framework, low-flux states represent windows of geometric transparency, and misaligned systems act as naturally unmasked laboratories. Consequently, the intrinsic duty cycle of extreme acceleration in AGN jets may be higher than inferred from flux-selected observations.
Read the full abstract

High Energy Astrophysical Phenomena

Spectral Hardening Revealed by Geometric De-boosting in the Masked Jet of PKS 2155$-$304

Blazar gamma-ray variability is predominantly stochastic and well described by red-noise processes. However, a subset of sources shows quasi-periodic oscillations (QPOs) on year-long timescales, whose physical origin remains debated. In high-synchrotron-peaked (HSP) blazars, departures from a single power-law gamma-ray spectrum, manifested as high-energy upturns in the GeV band, may probe emission mechanisms and the intrinsic duty cycle. We investigate the link between the 1.7 yr gamma-ray QPO in PKS 2155-304 and an exceptional spectral hardening event identified in the Fermi-LAT HSP blazar population. We analyze 17.4 years of Fermi-LAT data using 30-day binning, applying Singular Spectrum Analysis to mitigate red-noise effects and a Moving Block Bootstrap approach to quantify the correlation between photon flux and photon index. We find a statistically significant softer-when-brighter chromatic trend, supporting a geometric origin of the flux modulation. The spectral hardening event is phase-locked to the QPO trough, implying that the hardening signature is detectable only when geometrically boosted soft emission is suppressed at the flux minimum. We propose a Geometric Masking scenario in which jet geometry regulates the visibility of acceleration processes. These results favor a two-component jet structure and suggest that spectral hardening during low-flux states, even in non-periodic sources, may reveal jet physics otherwise obscured by relativistic amplification.
Read the full abstract

High Energy Astrophysical Phenomena

Coexistence of Chromatic Flares and an Achromatic QPO in the Gamma-ray Blazar PG 1553+113

The physical origin of quasi-periodic oscillations (QPOs) in blazars remains debated, with geometric and plasma-driven scenarios as the main competing interpretations. Discriminating between them requires probing variability beyond flux periodicity. We study the spectral evolution of the BL Lac object PG 1553+113 over its 2.2-year gamma-ray QPO cycle to constrain the mechanism driving the oscillation. In particular, we test whether the variability is chromatic (coupled to spectral changes) or achromatic (independent of spectral shape), allowing us to distinguish between plasma-driven and geometric scenarios. We analyze 17 years of Fermi-LAT data (2008-2025) with 30-day binning. To mitigate red-noise effects, we apply Singular Spectrum Analysis (SSA) to remove slow baseline trends and use a Block Bootstrap approach to quantify correlations between photon flux and photon index while preserving temporal dependence. We find a robust softer-when-brighter chromatic trend, atypical for high-synchrotron-peaked blazars such as PG 1553+113 and which, based on our analysis, physically corresponds to softer-when-flaring episodes, that persists after detrending and accounting for temporal autocorrelation. In contrast, no significant correlation is detected between the photon index and the QPO phase, indicating that the periodic modulation is effectively achromatic. The coexistence of plasma-driven chromatic flares and an achromatic QPO disfavors scenarios in which the periodicity arises from intrinsic jet processes. Instead, the results support a geometric origin for the QPO modulation, such as jet precession, where Doppler-factor variations modulate the flux without altering the intrinsic particle energy distribution.
Read the full abstract

Astrophysics of Galaxies

Detection of non-thermal radio emission components from the Orion Nebula: stellar jets, cloud collision or feedback from stellar winds?

The Orion Nebula is the closest high-mass star-forming region, making it an ideal laboratory to investigate physical processes in complex star-forming environments. At radio frequencies, the dominant emission mechanisms are thermal bremsstrahlung and non-thermal synchrotron. HII regions typically emit thermal radiation tracing the ionised gas; however, detecting and characterising non-thermal emission can provide insights into magnetic fields and the energy distribution of relativistic particles in star-forming regions. We have utilised the upgraded Giant Metrewave Radio Telescope (uGMRT) to study radio emission in the Extended Orion Nebula (EON) region. We present results from wide-band interferometric observations using uGMRT bands 3 and 4, probing a frequency range not covered by other sensitive radio interferometers. We produced deep continuum images with RMS noise levels of $\sim400\,\mu$Jy~beam$^{-1}$ in band 3 and $\sim200\,\mu$Jy~beam$^{-1}$ in band 4. We further generated in-band and broad-band spectral index maps using these images. To establish the robustness of the spectral index measurements, we conducted a detailed analysis using simulated uGMRT data. From the continuum spectral index analysis, we report the unambiguous presence of non-thermal radio emission in the EON region. To investigate its plausible origin, we correlated our results with multiwavelength observations, identifying a strong association between non-thermal emission and outflows from young stellar objects, while also exploring alternative explanations. In future, reliable broad-band radio spectral index measurements, together with dedicated multiwavelength observations, will be invaluable for resolving the origin of non-thermal emission in the Orion Nebula and other star-forming regions.
Read the full abstract

All authors Md Rashid, Nirupam Roy, Prasun Dutta, Jagadheep D. Pandian, Sarita Vig, Srijita Pal, Arnab Chakraborty, Samir Choudhuri

Earth and Planetary Astrophysics

NASA's Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA's Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible photometry and NIR spectroscopy of exoplanet host stars. Here we present an end-to-end simulation study quantifying Pandora's ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of Pandora spectrophotometry recover photospheric temperatures with typical uncertainties of ${\approx}30$ K, with no significant bias. Models with two spectral components (i.e., quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ${\approx}0.3$%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of $10^2{-}10^3$ ppm are reduced to ${\lesssim}10$ ppm, well below Pandora's expected transmission spectroscopy precision (30$-$100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the $10^3$ ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar-planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.
Read the full abstract

All authors Benjamin V. Rackham, Aishwarya R. Iyer, Dániel Apai, Peter McGill, Yoav Rotman, Knicole D. Colón, Brett M. Morris, Emily A. Gilbert, Elisa V. Quintana, Jessie L. Dotson, Thomas Barclay, Pete Supsinskas, Jordan Karburn, Christina Hedges, Jason F. Rowe, David R. Ciardi, Jessie L. Christiansen, Trevor O. Foote, Thomas P. Greene, Kelsey Hoffman, Rae Holcomb, Aurora Y. Kesseli, Veselin B. Kostov, Nikole K. Lewis, James P. Mason, Gregory Mosby Jr., Susan E. Mullally, Joshua E. Schlieder, Megan Weiner Mansfield, Luis Welbanks, Allison Youngblood

Astrophysics of Galaxies

MXDFz4.4: A LyC emitter 250Myr after the epoch of reionization and a first test of Ly-alpha morphology as a tracer of LyC escape at high redshift

Assessing the contribution of ionizing sources to cosmic reionization is a central goal of extragalactic astrophysics. Understanding and quantifying ionizing escape remains challenging near the epoch of reionization. We present the highest-redshift Lyman continuum (LyC) emitter detected to date, MXDFz4.4 at z=4.442 in the MUSE eXtremely Deep Field, observed only ~0.25Gyr after the end of reionization. A high confidence Ly-alpha line confirms the redshift. LyC flux is detected at 5.3sigma in the F435W filter with a flux of 4.2+/-0.8nJy, corresponding to a flux measurement at 5.2sigma. After correcting for the intrinsic production of LyC photons and the IGM opacity at z=4.44, we derive high escape fractions, ranging from 50 - 100%. We apply established low-redshift tracers of LyC escape and, for the first time at high redshift, promising Ly-alpha morphological tracers such as the halo fraction. SED fitting indicates the presence of a recent burst of star formation; we explore its impact on the production and escape of ionizing photons. Ly-alpha-based tracers of LyC escape reveal a complex scenario in which the recent burst strong influences LyC production and escape, combined with a more evolved stellar population. This interpretation is supported by UV diagnostics, including the star formation rate surface density and sSFR. Our results provide cautious support for the Ly-alpha halo fraction as a LyC escape tracer at high redshift. Considering the burst-driven enhancement in LyC production and escape, we conclude that stochastic star formation in the early Universe likely plays a significant role in the contribution of galaxies to cosmic reionization.
Read the full abstract

All authors Ilias Goovaerts, Marc Rafelski, Alexander Beckett, Grecco Oyarzùn, Annalisa Citro, Farhanul Hasan, Kalina V Nedkova, Calum Hawcroft, Anton M Koekemoer, Mitchell Revalski, Matthew J Hayes, Claudia Scarlata, Ray A Lucas, Norman A Grogin, David V Stark, Paolo Suin, Nor Pirzkal

Instrumentation and Methods for Astrophysics · Cross-listed from cs.LG

Augmenting representations with scientific papers

Astronomers have acquired vast repositories of multimodal data, including images, spectra, and time series, complemented by decades of literature that analyzes astrophysical sources. Still, these data sources are rarely systematically integrated. This work introduces a contrastive learning framework designed to align X-ray spectra with domain knowledge extracted from scientific literature, facilitating the development of shared multimodal representations. Establishing this connection is inherently complex, as scientific texts encompass a broader and more diverse physical context than spectra. We propose a contrastive pipeline that achieves a 20% Recall@1% when retrieving texts from spectra, proving that a meaningful alignment between these modalities is not only possible but capable of accelerating the interpretation of rare or poorly understood sources. Furthermore, the resulting shared latent space effectively encodes physically significant information. By fusing spectral and textual data, we improve the estimation of 20 physical variables by 16-18% over unimodal spectral baselines. Our results indicate that a Mixture of Experts (MoE) strategy, which leverages both unimodal and shared representations, yields superior performance. Finally, outlier analysis within the multimodal latent space identifies high-priority targets for follow-up investigation, including a candidate pulsating ULX (PULX) and a gravitational lens system. Importantly, this framework can be extended to other scientific domains where aligning observational data with existing literature is possible.
Read the full abstract

All authors Nicolò Oreste Pinciroli Vago, Rocco Di Tella, Carolina Cuesta-Lázaro, Michael J. Smith, Cecilia Garraffo, Rafael Martínez-Galarza