Peer Revenue Fully funded

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Saturday 4 April

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

What spectators do during inflation

The inflaton equation of motion including one loop radiative corrections from spectator fields is obtained. We consider a massless scalar conformally coupled to gravity and a massless fermion Yukawa coupled to the inflaton as models for spectators that do not feature gravitational particle production, their production during slow roll is solely a consequence of their coupling to the inflaton. The one-loop self energy and the fully renormalized equation of motion of the inflaton are obtained and solved explicitly for an inflaton potential $m^2\varphi^2/2$. The solution features Sudakov-type logarithmic secular terms, which are resumed via the dynamical renormalization group and compared to the solutions with a phenomenological friction term. During $N_e$ e-folds of slow roll inflation the inflaton evolves as $\varphi^{(0)}_{Isr}(t)\,e^{\frac{m^2\Gamma}{9H^3}\,N_e(t)}$ for the phenomenological friction term $\Gamma$ and $\varphi^{(0)}_{Isr}(t)\,e^{\Upsilon N^2_e}$ with $\Upsilon = -\frac{\lambda^2}{24\pi^2 H^2} ; \frac{y^2_R}{12\pi^2}$ for the radiative corrections from bosonic and fermionic spectators respectively where $\varphi^{(0)}_{Isr}(t)$ is the slow roll solution in absence of interactions, showing that a phenomenological friction term is not reliable. A generalization of the optical theorem to a finite time domain and cosmological expansion is introduced to obtain the distribution function $f(k,t)$ and total number of spectators produced during slow roll. $f(k,t)$ is peaked at superhorizon scales and the total number of particles grows $\propto e^{3N_e}$. A non-perturbative mean field theory is introduced to describe the self-consistent evolution of the inflaton coupled to spectators, its linearized version reproduces the self-energy, the inflaton equation of motion and the results on particle production.
Read the full abstract

High Energy Astrophysical Phenomena

Characterizing Gamma-Radio Delayed Flaring Activity from Blazars

Flaring activity from the jets of active galactic nuclei has been studied for several decades, closely related to the loading and evolution of the jet. In this work, we focus on the sub-hundred parsec jet region, well traced by non-thermal radio and gamma-ray emission. Only in recent years have light curves capturing the decade-long behavior of such sources become available for a large ensemble of objects. While previous studies have focused on a direct correlation or few-month lag between gamma-ray and radio activity, recent neutrino-bright blazars observed by the IceCube Neutrino Observatory present multi-year delays between initial gamma-ray activity and subsequent radio flares. In this work, we search for similar-timescale correlations between Fermi-LAT gamma-ray data and RATAN-600 radio data from ~100 blazars. We consider two gamma-ray bands, 100 MeV-1 GeV and 1 GeV-500 GeV, as well as the integral band, to compare correlations between potential opaque and unabsorbed regions of the jet. Gaussian process modeling is used for smooth light curve function prediction. We also analyze morphological AGN core data from the MOJAVE survey, forming a sub-selection to better illustrate potential dependence on location. In the broader selection, several sources exhibit delayed flares on the order of 1-3 years. In the stacked analysis, we find the highest correlation for a radio delay on the order of 180 days. The stacked correlation resulting from the MOJAVE sub-selection corresponds to a slightly smaller time lag. Delayed radio flares or extended radio emission appear to be notable features within the general blazar population.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Beyond $f(\phi)\mathcal{G}$: Gauss–Bonnet inflation with $\mu(\phi,X)$

Gauss–Bonnet inflation typically affects the dynamics over an extended portion of the trajectory, making it difficult to isolate a controlled imprint at CMB scales. We consider a trajectory-selective coupling \(\mu(\phi,X)\) that gates the Gauss–Bonnet sector in phase space, enabling the higher-curvature contribution to be localized within a finite e-fold window while remaining negligible elsewhere. We identify stable inflationary solutions consistent with this localization and enforce standard ghost and gradient stability conditions for both scalar and tensor perturbations. For these viable backgrounds we compute pivot-scale observables and examine their dependence on the overall Gauss–Bonnet strength and on the kinetic gating. The framework offers a controlled route for realizing localized higher-curvature effects with predictable consequences for CMB-scale measurements.
Read the full abstract

Solar and Stellar Astrophysics

Interaction between Winds from Weak-lined T Tauri Stars with Exoplanetary Magnetospheres

T Tauri stars, in more advanced stages of evolution, during the final accretion phase of stellar formation, exhibit intense stellar winds and surface magnetic fields with intensities around a kilogauss. With the growing interest in the search for rocky exoplanets with Earth-like dimensions, it is essential to deepen our understanding of the interaction between stellar winds and planetary magnetospheres. We investigated the interaction between stellar winds from 46 weak-lined T Tauri stars (WTTSs) and the magnetospheric protection of Earth-like planets located within their habitable zones. We employ two distinct stellar wind models, nonmagnetized and magnetized with both constant and resonant Alfven wave damping, to evaluate the pressure balance between the stellar wind and the planetary magnetic field. Our results show that the strong wind dynamic and magnetic pressures characteristic of WTTSs lead to systematically compressed planetary magnetospheres, significantly smaller than that of the present-day Earth. The analysis further indicates that planetary magnetospheric sizes increase with stellar age, following the decay of stellar magnetic activity, in agreement with previous findings for solar-type stars.
Read the full abstract

All authors Y. F. Tamburus, N. F. S. Andrade, G. R. C. Sampaio, V. Jatenco-Pereira

Astrophysics of Galaxies

Dynamical clock of the Helmi stream – Analysis of the clumping of stars in the orbital frequency-space

Reconstructing the assembly history of the Milky Way requires precise constraints on the dynamical age of its merger remnants – the time elapsed since a progenitor satellite was disrupted by the Galactic tidal force. We present a new framework to derive this dynamical age for disrupted stellar systems by extending the Fourier analysis of the orbital frequency distribution proposed by Gomez and Helmi. To overcome the smearing of frequency-space structures caused by observational noise, we introduce the Greedy Optimistic Clustering algorithm. This method allows for an optimistic exploration of the density contrasts in the orbital frequency space by taking into account the observational uncertainty in the data, effectively sharpening the signal required for age estimation. By applying this method to the Helmi stream, we derive a dynamical age of $6.8 \pm 0.8$ Gyr. Our derived accretion epoch is consistent with the observed kinematic properties of the Helmi stream. In particular, the marked asymmetry in the vertical velocity distribution – where approximately two-thirds of the stars have negative $v_z$ in the solar neighborhood – supports a relatively recent arrival. This suggests that the progenitor of the Helmi stream was accreted during an epoch of Galactic growth distinct from the much earlier Gaia-Sausage-Enceladus merger ($\sim 10$ Gyr ago). We validate our methodology using error-added mock simulations, demonstrating the reliability of our approach. Our results establish the Greedy Optimistic Clustering framework as a powerful chronometric tool for reconstructing the hierarchical assembly of the Milky Way using current and future high-precision astrometric datasets.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Dispersion Measure Distribution of Unlocalized Fast Radio Bursts as a Probe of the Hubble Constant

We present constraints on the Hubble constant ($H_0$) derived from the observed dispersion measure (DM) distribution of unlocalized fast radio bursts (FRBs). While localized FRBs with redshift measurements have been used to investigate the Hubble tension, their sample remains limited. Here we demonstrate that unlocalized FRBs—which are far more numerous—can independently constrain $H_0$ without requiring redshift information, as cosmic expansion imprints itself on their DM distribution. Analyzing a selected sample of 2124 unlocalized FRBs from the CHIME Catalog II, we obtain $H_0 = 69^{+17}_{-15}~\mathrm{km\,s^{-1}\,Mpc^{-1}}$ at the $1\sigma$ confidence level, corresponding to an uncertainty of about 22%. Disentangling the parametric degeneracy among $H_0$, the FRB spectral index $\alpha$, and the characteristic cutoff energy $E_*$ of the FRB energy distribution would reduce the fractional uncertainty in $H_0$ to 9%. This work constitutes the first $H_0$ measurement derived solely from the DM distribution of unlocalized FRBs, highlighting their potential as a new cosmological probe. Future joint analyses with localized FRBs promise even tighter constraints.
Read the full abstract

All authors Yang Liu, Jun-Jie Wei, Puxun Wu, Xue-Feng Wu

Solar and Stellar Astrophysics

Discovery of a Low-Mass Companion to the Accelerating Star HIP 53005 with Strongly Conflicting Mass Estimates

We present the discovery of a low-mass companion located at $\rho$ $\sim$ 0\farcs{}85 ($r_{\rm proj} \approx 62~au$) from the early-type 1.2 Gyr-old star HIP 53005 using direct imaging data from the Subaru and Keck Telescopes and astrometry from the Hipparcos-Gaia Catalog of Accelerations. The companion, HIP 53005 C, is a component of a multiple system also including a $\approx$ 12\farcs{}4-separation M dwarf companion inducing a negligible proper motion acceleration. HIP~53005 C's position on color-magnitude diagrams, the fit of its spectral energy distribution to atmosphere models, and its location on an empirical mass-magnitude diagram all suggest that it lies at the M/L transition and near the hydrogen-burning limit ($\sim80~M_{\rm Jup}$). However, our orbital fitting combining direct-imaging relative astrometry with proper motion acceleration favors a much higher dynamical mass of $\sim185\ M_{\rm Jup}$. An additional unseen, more closely-orbiting companion below the detection limit (at $\rho\lesssim0\farcs2$)) may explain this discrepancy. Alternatively, HIP~53005C could be a low-mass binary like Gliese~229Bab, making this system an intriguing laboratory for studying multiple star formation.
Read the full abstract

All authors Taichi Uyama, Thayne Currie, Jerry W. Xuan, Robert De Rosa, Masayuki Kuzuhara, Minghan Chen, Vito Squicciarini, Charles Beichman, Timothy D. Brandt, Vincent Deo, Olivier Guyon, Teruyuki Hirano, Markus Janson, Michael C. Liu, Dimitri Mawet, Julien Lozi, Stevanus Nugroho, Motohide Tamura, Sebastien Vievard, Danielle Bovie, Yasunori Hori, Hajime Kawahara, Takayuki Kotani, Yiting Li, Jason Wang

Astrophysics of Galaxies

Comparative Bayesian SED Fitting of PEARLSDG

The initial distance to PEARLSDG estimated from the Tip of the Red Giant Branch suggested it was an exotic isolated quiescent dwarf galaxy. We combine recent and archival Hectospec spectroscopy to place it at $z = 0.02843 \pm 0.00012$ ($D \approx 124$\,Mpc) within a galaxy group, revising the distance from 30\,Mpc to $\sim$124\,Mpc. We then carry out Prospector SED fitting using parametric and non-parametric star-formation histories sampled with dynesty, nautilus, and emcee, recovering metallicity $\log(Z/Z_\odot) = -0.44^{+0.35}_{-0.06}$, stellar mass $\log_{10}(M_*/M_\odot) = 9.25^{+0.02}_{-3.73}$, and dust attenuation $\hat{\tau}_V = 0.67^{+0.02}_{-0.05}$. The updated metallicity places PEARLSDG squarely on the standard mass–metallicity relation, resolving its former outlier status, with its quenched star-formation history consistent with environmental quenching in a group setting.
Read the full abstract

All authors Bar{\i}ş Bayraktar, Timothy Carleton, Rogier A. Windhorst, S. P. Willner, Christopher N. A. Willmer

Cosmology and Nongalactic Astrophysics

Is the $w_0w_a$CDM cosmological parameterization evidence for dark energy dynamics partially caused by the excess smoothing of Planck PR4 CMB anisotropy data?

We study the performance of the flat $\Lambda$CDM model and the dynamical dark energy parameterizations $w_0$CDM and $w_0w_a$CDM, in which the dark energy (DE) equation of state is either constant ($w=w_0$) or redshift-dependent [$w(z)=w_0+w_a z/(1+z)$], without and with a varying CMB lensing consistency parameter $A_L$, using combinations of Planck PR4 CMB data (PR4 and lensing), and a compilation of non-CMB data composed of baryon acoustic oscillation (BAO) data that do not include DESI BAO data, Pantheon+ type Ia supernova observations, Hubble parameter measurements $H(z)$, and growth rate $f\sigma_8$ data. We also compare results from earlier Planck PR3 data with those obtained using PR4 data in order to assess the stability of cosmological constraints. For the largest data combinations, PR3/PR4+lensing+non-CMB, the cosmological parameters inferred from PR3 and PR4 data are consistent, almost all differing by $1\sigma$ or less. For the $\Lambda$CDM$+A_L$ model, we have $A_L=1.087 \pm 0.035$ for PR3 and $A_L=1.053 \pm 0.034$ ($1.6\sigma$ above unity) for PR4, which indicates that the CMB lensing anomaly is reduced when PR4 data are used. For the $w_0 w_a$CDM parameterization, we find $w_0 = -0.863\pm0.060$ (quintessence-like) and $w_0+w_a=-1.37^{+0.19}_{-0.17}$ (phantom-like), suggesting that the current observations favor dynamical DE over a cosmological constant at about $1.8\sigma$. For the $w_0w_a$CDM$+A_L$ parameterization, we find $w_0=-0.877\pm 0.060$ and $w_0 + w_a =-1.29_{-0.17}^{+0.20}$, corresponding to a preference for dynamical DE over a cosmological constant of about $1.5\sigma$ and with $A_L = 1.042 \pm 0.037$ exceeding unity at $1.1\sigma$. These results indicate that while the PR4 data mildly favor a time-evolving DE, part of this preference may be associated with possible residual excess smoothing present in the Planck PR4 CMB anisotropy spectra (abridged).
Read the full abstract

Astrophysics of Galaxies

Investigating Extended Main-Sequence Turnoffs in Galactic Open Clusters

The extended main sequence (eMS) and extended main sequence turnoff (eMSTO) phenomena have been observed in some young and intermediate-age star clusters in the Milky Way and in the Magellanic Clouds. In this study, we conduct a survey of 53 galactic open clusters (OCs) to investigate the roles of stellar rotation, differential extinction, and cluster properties in the emergence of eMS and eMSTO. The projected rotational velocities are taken from the Gaia ESO spectroscopic survey and the Gaia DR3 line-broadening velocities. Stellar members of each OC are identified using the ML-MOC algorithm with Gaia DR3 astrometry. We divide clusters into four classes based on the color-rotation distribution, extinction, and MSTO morphology and report 14 clusters ($A_{\rm v} \lesssim 0.15$ mag, Class I) that exhibit split MS with fast and slow rotators populating the redder and bluer parts of MSTO. For the remaining clusters, differential extinction hampers the color-rotation distinction and also inflates MSTO width and therefore introduces a systematic offset in the MSTO-age relation. We also quantify the fraction of slow rotators among MSTO stars, finding a median value of $f_{\rm slow\, rot}^{v \sin i<100} \approx 0.41$ and the fraction reaching the spin-down limit, $f_{\rm slow\, rot}^{v \sin i<30}$, is $ \approx 0.08$. We find no statistically significant correlation between $f_{\rm slow\, rot}$ and either the binary fraction or cluster age.
Read the full abstract

Astrophysics of Galaxies

LensAgent: A Self Evolving Agent for Autonomous Physical Inference of Sub-galactic Structure

Probing dark matter distribution on sub-galactic scales is essential for testing the Cold Dark Matter ($\Lambda$CDM) paradigm. Strong gravitational lensing, as one of the most powerful approach by far, provides a direct, purely gravitational probe of these substructures. However, extracting cosmological constraints is severely bottlenecked by the mass-sheet degeneracy (MSD) and the unscalable nature of manual and neural-network modeling. Here, we introduce LensAgent, a pioneering training-free, large language model (LLM)-driven agentic framework for the autonomous physical inference of mass distributions. Operating as an autonomous scientific agent, LensAgent couples high-level logical reasoning with deterministic physical modeling tools, demonstarting successful reconstruction of mass distribution in SLACS Grade A strong lensing systems. This self-evolving architecture enables the robust extraction of sub-galactic substructures at scale, unlocking the cosmological potential of upcoming wide-field surveys such as the Rubin Observatory (LSST) and Euclid.
Read the full abstract

All authors Xiaotang Feng, Zihan Wang, Zilang Shu, Jean-Paul Kneib, Philip Torr

Cosmology and Nongalactic Astrophysics

Lattice simulations of scalar-induced gravitational waves from inflation

Scalar-induced gravitational waves (SIGWs) provide a powerful probe of inflationary dynamics on scales far smaller than those accessible to the cosmic microwave background and large-scale structure. In scenarios with a transient ultra-slow-roll (USR) phase, the curvature power spectrum can be strongly enhanced on small scales, potentially generating an observable stochastic GW background. In this regime, scalar dynamics during inflation can become nonlinear, challenging the validity of standard perturbative predictions. Existing semi-analytical calculations of SIGWs rely on the linear evolution of inflation fluctuations. In this work, we compute SIGWs from USR inflation using lattice simulations. We evolve the inflaton field nonlinearly during inflation and extract the curvature perturbation nonperturbatively, then simulate its post-reheating horizon re-entry by evolving the Newtonian potential linearly while retaining the full non-Gaussian structure of the initial conditions for the primordial fluctuations in the tensor source. For moderate non-Gaussianity, the semi-analytical prediction captures the correct order of magnitude of the GW signal but receives important corrections. When inflationary non-Gaussianities are large, it can fail dramatically in both amplitude and spectral shape, independently of the overall size of the tensor power spectrum. Our results show that reliable predictions of SIGWs in such scenarios require nonperturbative control of the inflationary scalar dynamics. The code used for this work is available at https://github.com/caravangelo/inflation-easy.git.
Read the full abstract

Solar and Stellar Astrophysics

Preferential Appearance of H$\alpha$ Moreton Waves along Photospheric Magnetic Network Boundaries

Moreton waves are rare chromospheric signatures of large-scale coronal disturbances, often associated with big flares and coronal mass ejections (CMEs). Using high-cadence, full-disk H$\alpha$ spectroscopic observations from CHASE, together with the EUV data from SDO/AIA and magnetograms from SDO/HMI, we analyzed a coronal EUV wave and an H$\alpha$ Moreton wave event associated with a filament eruption on 2024 July 29. The Moreton wave fronts are roughly cospatial with the fast-mode coronal EUV wave fronts, which propagate with a speed of $\sim$600 km s$^{-1}$. By comparing the Moreton wave fronts with photospheric features, we found that they preferentially appear along photospheric supergranule boundaries characterized by 1600 Å\ bright ridges, concentrated magnetic fields, and convective downflows. It is shown that the H$\alpha$ line profiles at the Moreton wave fronts are systematically redshifted. Gaussian fit yields a systematic downward Doppler velocity of 1.73 km s$^{-1}$. Using the bisector method, we further derived height-dependent Doppler velocities in the chromosphere. While there is an expected tendency for the downward velocity to decrease from 4.12 km s$^{-1}$ in the upper chromosphere to 1.60 km s$^{-1}$ in the lower chromosphere, it is intriguing to see an unexpected velocity enhancement in the lower chromosphere. We conjecture that coronal fast-mode MHD waves experience mode-conversion to slow-mode waves, which propagate along magnetic field lines of the magnetic canopy, forming preferential appearance of Moreton waves at magnetic networks, where the convective downflow contributes to the velocity enhancement in the lower chromosphere.
Read the full abstract

All authors C. H. Zhai, Y. W. Ni, J. H. Guo, P. F. Chen

Earth and Planetary Astrophysics

Influence of CO versus CH$_4$ on organic haze formation in atmospheres of diverse terrestrial exoplanets

Context. Terrestrial exoplanets are expected to host secondary, high-metallicity atmospheres derived from outgassing of volatiles such as N2, CO2, H2O, CH4, and CO. Photochemical organic hazes are likely to form in such environments, significantly affecting atmospheric observations and planetary habitability. Aims. We investigate haze formation in representative terrestrial exoplanet atmospheres and assess how CH4 versus CO as the primary carbon source affects haze production rates, particle properties, and chemical complexity. Methods. We performed six laboratory simulations by exposing gas mixtures at a few mbar to glow discharge at 300 K. Each atmosphere contained 75% N2, CO2, or H2O, 10% of each of the other two gases, and 5% CH4 or CO. Gas-phase products were analyzed with a residual gas analyzer, and solid products were characterized by production rate, particle density, atomic force microscopy, Fourier-transform infrared spectroscopy, and very high-resolution mass spectrometry. Results. CH4 experiments produced more diverse gas-phase species and much higher haze yields than the corresponding CO experiments. CO-derived hazes showed a narrow particle size range of 10-80 nm, whereas CH4-derived hazes were denser and chemically more complex. The identified molecular formulas suggest growth pathways linked to gaseous precursors such as HCN, CH2O, and C2H4. Conclusions. The atmospheric redox state critically controls haze formation in simulated terrestrial exoplanet atmospheres. CH4 is significantly more effective than CO in initiating organic growth, leading to higher haze production rates and greater chemical complexity. These results provide useful constraints for exoplanet atmospheric modeling and spectral interpretation, and further support the possibility that reducing atmospheres may facilitate prebiotic organic chemistry relevant to the emergence of life.
Read the full abstract

All authors Sai Wang, Zhengbo Yang, Haixin Li, Chao He, Yingjian Wang, Xiaoou Luo, Yu Liu, Sarah M. Horst, Sarah E. Moran, Veronique Vuitton, Laurene Flandinet, Patricia McGuiggan

Astrophysics of Galaxies

SPURS: Evidence for Clumpy Neutral Envelopes and Ionized IGM Surrounding Little Red Dots in Abell 2744 from Ultra-Deep Rest-UV Spectroscopy

Rest-frame ultraviolet (UV) spectra of Little Red Dots (LRDs) often show Ly$\alpha$ emission. Along with broad Balmer emission, LRDs are expected to produce broad Ly$\alpha$ emission. However, the large column density of neutral gas invoked to explain the Balmer break should significantly redshift and further broaden the Ly$\alpha$ line, making it challenging to detect without sensitive, moderate-resolution spectra. We present ultra-deep (29 hours) G140M JWST/NIRSpec observations covering the rest-UV of two LRDs in Abell2744 from the SPURS Cycle 4 Large Program. One of our targets is Abell2744-QSO1, a gravitationally-lensed LRD at $z=7.04$ with faint UV emission (M$_{\rm UV}=-16.9$), and the other source (UNCOVER-2476) is newly-confirmed at $z=4.02$ with a very bright UV continuum (M$_{\rm UV}=-19.6$). We find that Abell2744-QSO1 has a broad Ly$\alpha$ profile, along with narrow CIV, FeII$\lambda1786$, and OI$\lambda1302$ emission. The Ly$\alpha$ profile suggests an origin similar to the broad H$\alpha$, but the line is considerably less redshifted than expected from existing dense gas models. We show that the line profile can be explained if the dense neutral gas is clumpy, allowing Ly$\alpha$ to escape by scattering off of the clump surfaces. We find that UNCOVER-2476 has narrow [NeIV] emission, indicating either a hard radiation field or shocks. We confirm two close neighbors with Ly$\alpha$ emission around Abell2744-QSO1, indicating it traces a dense environment that may have ionized its surrounding IGM. We suggest that LRDs may preferentially trace bubbles carved by their dense environments, contributing to the prevalence of Ly$\alpha$ in the population.
Read the full abstract

All authors Mengtao Tang, Daniel P. Stark, Charlotte A. Mason, Zuyi Chen, Harley Katz, Max Gronke, Lukas J. Furtak, Seok-Jun Chang, Jorryt Matthee, Lily Whitler, Adi Zitrin, Ryan Endsley, Viola Gelli, Tamojeet Roychowdhury, Peter Senchyna, Michael W. Topping, Meng Zhang

Friday 3 April

Instrumentation and Methods for Astrophysics

FluxMC: Rapid and High-Fidelity Inference for Space-Based Gravitational-Wave Observations

Bayesian inference in the physical sciences faces a fundamental challenge: the imperative for high-fidelity physical modeling often clashes with the intrinsic limitations of stochastic sampling algorithms. Complex, high-dimensional parameter spaces expose the universal vulnerability of conventional methods, e.g., Markov Chain Monte Carlo (MCMC), which struggle with the prohibitive costs of likelihood evaluations and the risk of entrapment in local optima. To resolve this impasse, we introduce FluxMC (Flow-guided Unbiased eXploration Monte Carlo), a machine learning-enhanced framework designed to shift the inference paradigm from blind local search to globally guided transport. It integrates Flow Matching with Parallel Tempering MCMC, effectively combining the global foresight of generative AI with the rigorous asymptotic convergence and local robustness of temperature-based sampling. We showcase the efficacy of this framework through the lens of space-based gravitational-wave (GW) astronomy – a field representing the frontier of challenging parameter inversion. In the analysis of massive black hole binaries using high-fidelity waveforms (IMRPhenomHM), FluxMC achieves robust convergence in under five hours, whereas traditional Parallel Tempering MCMC fails to converge even after hundreds of hours, yielding high Jensen-Shannon divergences (JSD) of $O(10^{-1})$. Our method reduces the distributional error by two to three orders of magnitude. Furthermore, for computationally efficient models (IMRPhenomD), it eliminates systematic biases caused by local-optima entrapment. Ultimately, FluxMC removes the necessity to compromise between model accuracy and analysis speed, establishing a new computational foundation where scientific discovery is limited only by observational data quality, not by algorithmic capacity.
Read the full abstract

All authors Bo Liang, Chang Liu, Hanlin Song, Tianyu Zhao, Minghui Du, He Wang, Haohao Gu, Sensen He, Yuxiang Xu, Wei-Liang Qian, Li-e Qiang, Peng Xu, Ziren Luo, Mingming Sun

Earth and Planetary Astrophysics · Cross-listed from math.DS

Effective stability estimates close to resonances with applications to rotational dynamics

We consider nearly-integrable Hamiltonian systems defined over a non-resonant domain. In the neighborhood of resonances, we use Nekhoroshev-like estimates to provide effective stability bounds for the action variables over long time. The applicability conditions of these estimates allow some freedom in the choice of parameters. Hence, we develop an optimization algorithm for choosing parameters that maximize the stability time. To further improve the stability estimates, we use perturbation theory to reduce the norm of the perturbing function. We implement this procedure (effective stability estimates and perturbation theory) to analyze the stability of sequences of irrational (Diophantine) frequencies converging to frequencies corresponding to resonances. We consider two applications to models describing problems of rotational dynamics in Celestial Mechanics: the spin-orbit problem, described by a 1D time-dependent Hamiltonian, and the spin-spin-orbit model, described by a 2D time-dependent Hamiltonian. We show stability results for orbits close to the main resonances associated with such models.
Read the full abstract

Solar and Stellar Astrophysics

Prediction of Magnetic Flux Evolution During Solar Active Region Emergence using Long Short-Term Memory Networks

Solar active regions (ARs) are the primary drivers of space weather events, making their early prediction crucial for operational forecasting systems. We develop machine learning models capable of predicting the evolution of magnetic flux during AR emergence using 1D time series of the continuum intensity and solar oscillation power maps for 53 active regions and their surrounding quiet-Sun areas. Each observable is sampled over a fixed 30.66{\deg}x30.66{\deg} field of view. These observations capture the temporal evolution of each active region and serve as inputs for training and validation of our MagFluxLSTM and MagFluxEnc-Dec models. The MagFluxLSTM architecture implements a single-stage standard Long-Short Term Memory (LSTM) network. MagFluxEnc-Dec represents an LSTM encoder-decoder with teacher forcing. To test and evaluate the models' performance, we use the continuum intensity and oscillation power maps (calculated for several frequency bands from Doppler velocity) as input to predict the magnetic flux. Among the top 100 hyperparameter configurations ranked by validation derivative RMSE, 98% correspond to MagFluxLSTM, compared to only 2% for MagFluxEnc-Dec. Thus, although the MagFluxEnc-Dec architecture has higher model complexity, it leads to poorer generalization to ARs outside the training set and less stable training than the simpler MagFluxLSTM, which can predict magnetic flux emergence 3-10 hours in advance within a 12-hour prediction window in both experimental and operational-type settings for the 5 testing active regions.
Read the full abstract

All authors Eren Dogan, Spiridon Kasapis, Sarang Patil, Jonas Tirona, John Stefan, Irina Kitiashvili, Mengjia Xu, Alexander Kosovichev

Astrophysics of Galaxies

The morphologies of present-day galaxies in the COLIBRE simulations

The diversity of galaxy morphologies and their relations with galaxy and halo properties is fundamental to understanding galaxy formation. Cosmological simulations of representative volumes can help disentangle the origin of observed correlations, but most suffer from two main limitations that affect morphologies: an over-pressurised interstellar medium and spurious interactions between stellar and dark matter particles. We present an overview of galaxy morphologies in the COLIBRE simulations, which address these limitations and reproduce many observed galaxy scaling relations. To quantify galaxy morphology, we use four (strongly-correlated) theory-space metrics, three kinematic and one spatial. We explore how different choices and limitations affect these indicators, including luminosity- versus mass-weighting, aperture size and shot noise. Overall, we find good convergence in present-day morphologies across two orders of magnitude in mass resolution. COLIBRE predicts that kinematic morphology correlates strongly with stellar mass and colour, and that galaxies with stellar masses of $\approx(1-2)\times 10^{10}\,\mathrm{M}_{\odot}$ tend to be the most rotationally-dominated. At fixed stellar mass, the morphology of central galaxies correlates weakly with the properties of their host halo. Morphology correlates more strongly with internal galaxy properties, with more disky galaxies being more gas-rich, having higher star formation rates and exhibiting younger and more extended stellar populations. Other properties, like the mass of the most massive black hole, the fraction of stars that are accreted and stellar metallicity, also correlate with morphology, but with correlation strengths sensitive to the stellar mass of the galaxy and whether it is a central or satellite.
Read the full abstract

All authors Victor J. Forouhar Moreno, Joop Schaye, Matthieu Schaller, Aaron Ludlow, Robert J. McGibbon, Alejandro Benítez-Llambay, Evgenii Chaikin, Carlos S. Frenk, Filip Huško, Sylvia Ploeckinger, Alexander J. Richings, James W. Trayford

Instrumentation and Methods for Astrophysics · Cross-listed from physics.optics

Silicon Photonics-based Heterodyne Interferometric Imager for free-space imaging

This paper reports on the design, fabrication, and demonstration of a silicon photonics based heterodyne interferometric imaging system. The photonic integrated circuit (PIC) can perform one-dimensional spectroscopy for unique input spectrums using a single baseline within its 91 available baselines. The PIC uses polarization diversifying gratings to separate incoming light into two distinct polarizations, an on-chip 2x4 optical hybrid, and a strong local oscillator (LO) to perform the heterodyne measurements. The optical hybrids combine the input signals with the LO and splitting them into 2 components pairs for phase sensitive measurements. Furthermore, the PIC can perform 2-D image reconstruction by combining many baseline pairs to measure the visibility of a simple target. These demonstrations show the PIC's capabilities for 1-D spectroscopy and 2-D imaging applications.
Read the full abstract

All authors Humphry Chen, Mingye Fu, Shun-Hung Lee, Shelbe Timothy, Lawrence Shing, Gopal Vasudevan, Tony Kowalczyk, Neal Hurlburt, Sung-Joo Ben Yoo

High Energy Astrophysical Phenomena

Milliarcsecond-scale spectrum of the persistent radio source associated with FRB 20190417A and constraints for FRB 20181030A

We aim to confirm the compact nature and constrain the radio spectra of candidate persistent radio sources (PRSs) associated with repeating fast radio bursts (FRBs). We performed European VLBI Network (EVN) observations at 5 and 8 GHz targeting two candidates identified in a recent VLA survey. We measured flux densities and upper limits at milliarcsecond resolution and combined them with published VLBI data at lower frequencies to derive spectral constraints. We detect a compact source associated with FRB 20190417A at 5 GHz with a flux density of $150\pm45$ uJy, while no detection is obtained at 8 GHz. The source is unresolved and has a brightness temperature $T_{\rm b}>10^{5}$ K, confirming its non-thermal nature. Combining our measurement with VLBI data at 1.4 GHz, we derive a spectral index $\alpha = -0.19 \pm 0.29$, consistent with a nearly flat spectrum. This makes FRB 20190417A only the second PRS with a spectral index constrained using VLBI data. The inferred luminosity places the source on the proposed $L_{\nu}$-|RM| relation. Including this source yields a scatter of $\sigma_\Delta = 0.65$, corresponding to $\hat{\alpha}|\epsilon| = 1.5 \pm 0.7$, consistent with forward shocks in the free-expansion phase or young pulsar wind nebulae. For the candidate PRS associated with FRB 20181030A, we report upper limits of 80 uJy at 5 GHz and 150 uJy at 8 GHz, corresponding to $L_{5\,\mathrm{GHz}} \lesssim 3.8 \times 10^{25}\ {\rm erg\ s^{-1}\ Hz^{-1}}$, and implying a steep spectral index ($\alpha \lesssim -1.2$) if the VLA emission arises from a compact component. Our results highlight the importance of VLBI in isolating compact emission from FRB engines and provide one of the few spectral constraints for PRSs at milliarcsecond resolution. The consistency of FRB 20190417A with the $L_{\nu}$-|RM| relation supports a nebular origin for the persistent emission.
Read the full abstract

All authors G. Bruni, L. Piro, Y.-P. Yang, L. Nicastro, A. Rossi, E. Palazzi, E. Maiorano, S. Savaglio, B. Zhang

Instrumentation and Methods for Astrophysics

A Foundation for Gravitational-Wave Population Inference within the LISA Global Fit

Population inference in gravitational-wave astronomy allows us to connect individual detections to the astrophysics of compact objects and their environments. Current approaches employed for population inference with LIGO-Virgo-KAGRA data approximate evaluation of the hierarchical population likelihood via post-processing of individual-event posteriors. However, the case of the Laser Interferometer Space Antenna (LISA) will be more complex for two main reasons: the transdimensional "global fit" approach to LISA data analysis which models all signals and noise simultaneously, and the presence of both individually-resolved signals and the unresolved stochastic “Galactic foreground" arising from the Galactic binary population, which induces a circular dependence between the resolved and unresolved systems and our ability to detect the former. These challenges are not without opportunity; LISA's data will contain every mHz compact binary in the Milky Way – either individually or within the Galactic foreground – with great potential for Galactic and stellar astrophysics. We therefore propose an alternative approach: direct evaluation of the full hierarchical population likelihood within the LISA global fit. We develop a statistical formalism for joint inference of individually-resolved gravitational-wave sources, an unresolved stochastic foreground, and a shared, underlying astrophysical population, present PELARGIR, a prototype GPU-accelerated population inference module for the LISA global fit, demonstrate the formalism and PELARGIR via a toy model analysis, and lay out a roadmap towards an astrophysically-motivated LISA global fit with embedded population inference. While we apply the formalism here to the population of LISA Galactic binaries, it is applicable across the gravitational-wave spectrum with use cases in pulsar timing and next-generation terrestrial observatories.
Read the full abstract

All authors Alexander W. Criswell, Sharan Banagiri, Vera Delfavero, Maria Jose Bustamante-Rosell, Stephen R. Taylor, Robert Rosati

Astrophysics of Galaxies

The Role of Baryonic and Dark Matter in Bar Kinematics

Simulations predict that bars in galaxies should slow down over time. This is often attributed to the exchange of angular momentum between the bar and other regions of the galaxy, such as the outer disc and dark matter halo, which implies that galaxies with a more massive halo or disc should be able to slow down the bar more efficiently. However, observational evidence for this process has been limited. In this work, we provide observational support for the slowing down of bars as predicted by simulations. We combine bar kinematics measurements obtained with the Tremaine-Weinberg method and host galaxy mass estimates derived from Jeans anisotropic models for a sample of 30 galaxies from the MaNGA survey. We find a statistically significant anti-correlation (>4sigma) between the bar pattern speed and both the stellar and total dynamical mass, which suggests that the slowest bars reside in the most massive galaxies. However, while the slope of the best-fit line between the pattern speed and dark matter mass is negative, it is not statistically significant (2.43sigma). We also find that bars with lower pattern speeds have more extended NFW dark matter profiles with lower central densities. Additionally, we find statistically significant correlations (>3sigma) between the corotation radius and the stellar mass, dark matter mass, and total dynamical mass. Finally, we find no significant correlations that involve the dark matter fraction or R, likely due to the inherent challenges associated with measuring these specific parameters accurately.
Read the full abstract

All authors Tobias Géron, Kai Zhu, Steph Campbell, Chris Lintott, R. J. Smethurst, Behzad Tahmasebzadeh

Solar and Stellar Astrophysics

Deep Adaptive Optics Imaging Rules Out a Helium Star Companion to PSR J1928+1815

PSR J1928+1815 is a 10.55 ms millisecond pulsar in a 3.6 hr orbit with a massive ($1.0$-$1.6\,M_{\odot}$) companion that produces extended radio eclipses. The companion, proposed to be a stripped helium star, is undetected in optical and infrared surveys. We present deep near-infrared imaging using Keck/NIRC2 with laser guide star adaptive optics. No source is detected at the pulsar position down to a $5\sigma$ limit of $K_s \approx 21.3$. Using stripped-star atmosphere models and conservative extinction estimates, we show that any plausible helium star companion would have been detected, ruling out this interpretation. A massive white dwarf (WD) companion remains consistent with the non-detection. We consider two possible origins for the eclipses: (1) absorption in a wind driven by a young, hot WD, and (2) material ablated from the WD by the pulsar. The former can naturally arise following Case BB mass transfer, which produces $\sim 1.2\,M_\odot$ WDs capable of sustaining winds of $\dot{M} \gtrsim 10^{-12}$-$10^{-13}\,M_\odot\,{\rm yr}^{-1}$ for $\sim 10^4$-$10^5$ yr, sufficient to obscure the pulsar at GHz frequencies. The latter requires efficient coupling of the pulsar's spin-down luminosity to the companion to drive the needed mass loss, which may be difficult to achieve. If the eclipse is powered by a WD wind, the system is likely observed in a short-lived phase; alternatively, if the companion is an older WD, the origin of the eclipsing material remains unclear. The apparent uniqueness of PSR J1928+1815 is consistent with a short detectability lifetime, though formation rate estimates remain uncertain.
Read the full abstract

All authors Pranav Nagarajan, Kareem El-Badry, Jim Fuller, Yunlang Guo, Thomas M. Tauris

Cosmology and Nongalactic Astrophysics

Big Bang Nucleosynthesis Constraints on the CCC+TL Cosmology

We investigate whether Big Bang nucleosynthesis (BBN) remains compatible with the Covarying Coupling Constants plus Tired Light (CCC+TL) cosmology. In this framework, only quantities with explicit length dimensionality covary through a universal scaling function $f \left( z \right)$, while dimensionless constants and dimensionless ratios remain invariant. At the redshifts $z$ relevant to BBN, $f \left( z \right )$ approaches a constant plateau $f_{\text{max}} \left( z \right) \simeq 3$, and the tired-light contribution is negligible, so the early-time dynamics reduce to a global rescaling of dimensioned quantities. In particular, the Hubble expansion rate $H$ at fixed temperature $T$ satisfies $H_{\text{CTL}} \left( T \right) = f^{-1}_{\text{max}} H_{\Lambda\text{CDM}}\left( T\right)$, implying a longer cooling time $\Delta t$ between weak freeze-out and the onset of nucleosynthesis by the same factor (CCC+TL labeled as $\textit{CTL}$). We find that BBN predictions are preserved provided the relevant interaction rates $\Gamma$ and decay rates governing the neutron lifetime ${\tau}_n$ share the same plateau scaling as $H$, so that governing combinations such as $\Gamma\text{/}H$ and $\text{exp} \left( -\Delta t \text{/} {\tau}_n \right)$ remain invariant. Implementing these plateau rescalings in the Kawano/NUC123 network (via a single control parameter $\texttt{fctl} \equiv f_{\text{max}}$) yields identical light-element abundances for $\texttt{fctl}= 1$ ($\Lambda$CDM) and $\texttt{fctl} = 3\left( \text{CCC+TL} \right)$ to within $10^{-3} - 10^{-4}$ level, consistent with numerical rounding. We also illustrate that adopting the lower late-time CCC+TL baryon density from the Pantheon+ data fit can reduce the ${}^7\text{Li}$ discrepancy but simultaneously increases D/H, implying that BBN alone does not select between the late-time baryon-density inferences considered here.
Read the full abstract

Astrophysics of Galaxies

Silicon, sulfur and iron in the interstellar medium: a high-resolution X-ray spectral study of GX 340+0

High-resolution X-ray spectroscopy provides a powerful probe of the interstellar medium (ISM), giving direct access to the composition and physical state of dust grains and atomic species in dense environments. We present a study of the gas and dust along the line of sight to the bright low-mass X-ray binary GX 340+0, which samples higher-density gases in the ISM. Using deep Chandra/HETG spectra, we characterize X-ray absorption fine structure from dust, gas absorption lines, and the optical depths of the Si, S, and Fe K-edges. By fitting these three edges simultaneously, we reduce degeneracies in the dust composition and find that amorphous olivine dominates the fractional contribution among the dust columns ($\sim$65%), followed by metallic iron ($\sim$19%), iron sulfides (pyrrhotite and troilite; $\sim$10%), and fayalite ($\sim$5%), with the remaining species contributing only a few percent in total. From the inferred stoichiometry, we estimate that $\sim$74% of Fe is associated with silicates, $\sim$8% with sulfides, and $\sim$18% with metallic iron, suggesting that Fe is predominantly incorporated in iron rich silicate grains along this sightline. We detect S ii absorption and infer a sulfur dust fraction of $\sim$35%. We also detect absorption structure near the Ca and Ar K edges, highlighting the need for improved atomic photoabsorption data. The Chandra/HETG spectral resolution remains essential to disentangle dust and gas contributions at the Si and S K edges, providing a benchmark for dust characterization in high-density regions in the ISM.
Read the full abstract

All authors Daniele Rogantini, Claude Canizares, Elisa Costantini, Liyi Gu, Missagh Mehdipour, Ioanna Psaradaki, Norbert S. Schulz, Sascha T. Zeegers

Instrumentation and Methods for Astrophysics

GOPREAUX I: Open-source Code and Data to Model Multi-wavelength Emission of Extragalactic Transients using Gaussian Processes

Contemporary all-sky surveys have observed thousands of extragalactic transients in the nearby universe, and upcoming surveys will discover exponentially more at higher redshifts. With these large samples, population-level analysis of the photometric behavior of different transient classes is now possible, allowing for photometric classification and physical parameter inference from relatively sparse individual light curves. To enable such studies, we introduce Gaussian process Optimized Photometric Regression of Extragalactic Archival Ultraviolet-infrared eXplosions, a.k.a GOPREAUX–a Python package for Gaussian Process Regression of multi-wavelength transient photometry. Our modeling is unique in that it interpolates transient emission across phase and wavelength in a non-parametric, data-driven way. This allows for predictions of light curves and spectra at higher redshifts, where the rest-frame ultraviolet (UV) emission is redshifted into the observer-frame optical or infrared. To this end, we aggregate a sample of almost 1,300 transients observed in the UV and optical with the Neil Gehrels Swift Telescope, complemented with additional optical and infrared coverage from surveys such as ZTF and open-source data releases. Our sample includes 275 Type II SNe, 172 stripped-envelope SNe, 72 superluminous SNe, and 58 tidal disruption events, among other classes. Our code and reduced photometry–comprising over 146,000 photometric observations–are available as open-source software and data products. Here we discuss our sample criteria, data reduction and modeling methodologies, the multi-wavelength light curves and spectral templates produced by our models, and the future directions in photometric classification and physical parameter inference this code and data repository enables.
Read the full abstract

All authors C. Pellegrino, T. A. Pritchard, M. Modjaz, A. Crawford, S. Khakpash, F. Bianco

Astrophysics of Galaxies

Holes in the BH$^\star$? AGN signatures in the FUV spectrum of a black-hole dominated Little Red Dot at $z=7.04$

It has been suggested that "Little Red Dots" (LRDs) might be accreting black holes enshrouded by dense gas in a nearly closed geometry, which completely covers the central black hole, leading to an atmosphere-like structure known as the "black-hole star" ($\rm BH^\star$). We test this scenario by analysing new JWST spectroscopy in the far ultraviolet (FUV, rest-frame) of the prototypical LRD Abell2744-QSO1, at $z=7.04$. We found the presence of broad Ly$\alpha$ emission with an FWHM of $\sim 1000$ km/s, and detections of OI, CIV, and/or FeII emission lines. The NIRCam imaging and NIRSpec slit images indicate that the low-velocity component ($v\lesssim 200$ km/s) of Ly$\alpha$ is likely spatially extended, but the high-velocity component ($v\gtrsim 200$ km/s) of Ly$\alpha$ remains unresolved. Based on the multi-component kinematics and flux of Ly$\alpha$ relative to Balmer lines, we conclude that the observed line profile is unlikely to be broadened by subsequent resonant scattering through the interstellar medium. This suggests that the high-velocity component of Ly$\alpha$ originates in the broad-line region, although resonant scattering in the dense gas likely makes Ly$\alpha$ broader than H$\alpha$ as observed. The nebular features of this LRD indicate that there is at least one relatively optically thin direction where Ly$\alpha$ can escape from the broad-line region (BLR). We also found indications that photons from the BLR are powering fluorescence of FeII and OI on a larger physical scale. The FUV features thus challenge the fully-covered geometry interpretation and suggest that there are "holes" in the $\rm BH^\star$, or the absorbing medium is simply clumpy.
Read the full abstract

All authors Xihan Ji, Gabriele Pezzulli, Francesco D'Eugenio, Roberto Maiolino, Stefano Carniani, Sandro Tacchella, Gareth Jones, Aaron Smith, Joris Witstok, Andrew C. Fabian, Sophia Geris, Anishya Harshan, Yuki Isobe, Lucy R. Ivey, Ignas Juodžbalis, Robert Pascalau, Jan Scholtz, Callum Witten

Earth and Planetary Astrophysics

The PLATO Input Catalogue of targets (tPIC) for the first Long Pointing Field

The ESA PLAnetary Transits and Oscillations of stars (PLATO) mission is designed to detect terrestrial planets in the habitable zones of solar-type stars. Owing to telemetry constraints, the selection of PLATO targets must be performed in advance. We present the first public release of the PLATO Input Catalog of targets (tPIC2.2), which provides the list of stars that will be observed during PLATO's first Long-duration Observation Phase field at South (LOPS2) as part of its core science program. We used astrometric and photometric data from Gaia data release 3 (DR3), together with three-dimensional maps of the local interstellar medium, to identify stars belonging to the PLATO stellar samples in accordance with the mission requirements. The tPIC comprises 217,741 stars, including 202,315 FGK dwarfs and subgiants, 15,037 M dwarfs, and 789 known planet-host stars. The median distances of the samples are 512 pc for FGK stars and 133 pc for M dwarfs. We estimated interstellar reddening for almost all targets and developed an algorithm to infer fundamental stellar parameters (effective temperature, radius, and mass) homogeneously from astrometric and photometric observables. The tPIC fulfills all the science requirements of the PLATO mission. It also includes a list of stars that host known exoplanets, confirmed or still candidate, located within the LOPS2 field.
Read the full abstract

All authors M. Montalto, G. Piotto, P. M. Marrese, L. Prisinzano, S. Marinoni, V. Granata, J. Cabrera, V. Nascimbeni, S. Desidera, V. Adibekyan, S. Ortolani, E. Alei, C. Aerts, G. Altavilla, K. Belkacem, S. Benatti, A. Börner, M. Deleuil, M. Fabrizio, L. Gizon, M. J. Goupil, M. Günther, A. M. Heras, D. Magrin, L. Malavolta, J. M. Mas-Hesse, H. P. Osborn, I. Pagano, C. Paproth, D. Pollacco, R. Ragazzoni, G. Ramsay, H. Rauer, S. Udry

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

Multimessenger Signatures of Tilted, Self-Gravitating, Black Hole Disks

We perform fully relativistic GRMHD simulations of magnetized, self-gravitating black hole-disk (BHD) systems in which the black hole spin is misaligned with the disk angular momentum. Massive disks (disk to BH mass ratios of $16-28\%$) around rapidly rotating black holes ($\chi\lesssim 0.97$) develop a nonaxisymmetric instability for tilt angles from $0^\circ$ to $180^\circ$. Magnetic stresses damp, but do not completely suppress, the nonaxisymmetric instability, and corresponding gravitational wave (GW) emission, in aligned systems, while they enhance it in antialigned BHDs: MRI-driven turbulence enhances angular momentum transport and accelerates nonlinear instability evolution in misaligned configurations. All models launch magnetically driven jets consistent with the Blandford-Znajek (BZ) mechanism, with collimation depending on spin orientation. The GWs reflect strong nonaxisymmetric structure from a persistent $m=1$ mode. The coupling between fast MRI and the slower nonaxisymmetric instability growth governs the outcome, with tilt controlling how MRI modifies the global mode. These simulations provide the first self-consistent GRMHD treatment of tilted, self-gravitating BHD systems and support their role as multimessenger sources.
Read the full abstract