Peer Revenue Impact factor: volatile

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics

Friday 28 August

Solar and Stellar Astrophysics

A decade of radial-velocity monitoring of rho Leo: moment analysis and periodic variability

V.A. Checha, A. Aret, I. Kolka et al. (+6 more)

We investigate the origin of long-term spectroscopic and photometric variability in the blue supergiant rho Leo, with particular emphasis on distinguishing between intrinsic pulsations and variability induced by a possible companion. Our analysis is based on an 11.5-year spectroscopic time series obtained at Tartu Observatory, complemented by high-cadence, high-resolution spectroscopy from the Hertzsprung SONG telescope and space-based photometry from TESS. We studied line-profile variability using normalised moments of the HeI 6678, HeI 5875, and SiIII 4552 lines. Periodic signals were identified using the generalised Lomb-Scargle periodogram with iterative pre-whitening, and their temporal stability was examined with the weighted wavelet Z-transform. We detect a persistent periodic signal at P = 16.46 d in the first and third moments, present throughout the full observing interval, with a radial-velocity amplitude of 3.9 km/s. This signal is also present in the SONG data and is visible in multiple spectral lines, indicating a global origin. Photometric observations reveal a dominant variability timescale near 33 d, approximately twice the spectroscopic period. The stable 16.46-day period, present throughout the entire observing interval, most likely results from non-radial pulsations of the supergiant. A binary origin of the signal is not excluded, but distinguishing between these scenarios is complicated by the supergiant's complex variability pattern.
Read the full abstract

All authors V.A. Checha, A. Aret, I. Kolka, T. Liimets, V. Mitrokhina, A. Kasikov, T. Eenmäe, S. P. D. Borthakur, H. Ramler

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

Devoured by a Hairy Gargantua: Probing Massive Scalar Charges with Non-minimal Curvature Coupling with Extreme-Mass-Ratio Inspirals

Leif Lui, Adrian Ka-Wai Chung and Alejandro Torres-Orjuela

Ultralight scalar fields nonminimally coupled to curvature can endow rotating black holes with stationary hair and alter nearby orbits. Using recently constructed hairy black-hole spacetimes, we model the transition to plunge of an extreme-mass-ratio inspiral. For scalar masses $\mu M=0.01$ and $0.1$ and dimensionless curvature coupling $\zeta=10^{-3},10^{-4}$, we find $\mathcal{O}(10)$ rad of gravitational-wave dephasing relative to general relativity over $\mathcal{O}(10^3)$ orbits. A multimode Fisher forecast at $\mu M=0.2$, where neighbouring hair solutions permit a numerical $\mu$-derivative, suggests that $\zeta$ and $\mu$ could be measured to precisions of $\simeq0.2\%$ and $0.6$–$1.1\%$, respectively, for a source with signal-to-noise ratio of about 80 in LISA.
Read the full abstract

High Energy Astrophysical Phenomena

The Geometry of Ultra-Fast Outflows Probed by Soft X-ray Variability in PDS 456

Riki Sato, Kouichi Hagino, Toshiya Iwata et al. (+20 more)

Constraining the location and geometry of ultra-fast outflows (UFOs) is essential for identifying where they are launched and how they are accelerated. We investigate the soft X-ray variability of the luminous quasar PDS 456 using simultaneous March 2024 observations with XRISM/Xtend and NuSTAR. A model-independent comparison between the flare and quiescent phases reveals spectral variability around 1 keV in the rest frame, while the hard X-ray spectral shape remains nearly unchanged. Broadband spectral fitting shows that the soft X-ray structure is well described by a partial-covering low-ionization UFO with $\log (\xi/(\mathrm{erg~cm~s^{-1}})) \simeq 3.1$ and $v_{\mathrm{out}} \simeq 0.30c$. Time-sliced spectral analysis further reveals significant variability in the covering fraction of this absorber. Interpreting this variability as transverse motion across the X-ray source, we constrain the crossing velocity to be $v_{\mathrm{cross}} \lesssim 5 \times 10^{-3}c$ and derive a lower limit on the absorber distance of $r \gtrsim 4 \times 10^3 R_{\mathrm{g}}$. This location is substantially farther out than the high-ionization UFO previously inferred at ~200-600 $R_{\mathrm{g}}$, while the two phases have comparable outflow velocities. The resulting velocity-distance structure disfavors a self-similar magnetocentrifugal wind and instead suggests either radiation-pressure acceleration following a Castor-Abbott-Klein-like velocity law or compact magnetic acceleration through magnetic reconnection. These results demonstrate that soft X-ray partial-covering variability can provide a geometrical probe of UFOs and directly connect spectral variability to wind acceleration.
Read the full abstract

All authors Riki Sato, Kouichi Hagino, Toshiya Iwata, Ehud Behar, Valentina Braito, Chris Done, Misaki Mizumoto, James N. Reeves, Rozenn Boissay-Malaquin, Pierpaolo Condó, Luigi C. Gallo, Adam G. Gonzalez, Alfredo Luminari, Aiko Miyamoto, Ryuki Mizukawa, Hirokazu Odaka, Atsushi Tanimoto, Makoto Tashiro, Francesco Tombesi, Yerong Xu, Satoshi Yamada, Tahir Yaqoob, Aya Bamba

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

Dynamical friction on black holes in scalar field environment

Adith Praveen and Dawood Kothawala

Astrophysical black holes exist within non-vacuum environments, and their motion through these environments generically result in a force on these black holes through processes such as dynamical friction and Bondi accretion. We explore these forces numerically for black holes moving through a complex scalar field of mass $m$ with constant acceleration $a$, while also revisiting the constant velocity case considered in the existing literature. The former is modeled by the C-metric, while we use Painlevé-Gullstrand metric with an additional divergence and vorticity free velocity field to mimic the constant velocity motion. Our simulations reveal several novel and interesting aspects of the drag force in both cases. For the constant velocity case, the force saturates at late times, with a value that depends on velocity in a manner distinct from known dependence. For the constant acceleration case, the force increases to a maximum value $F_\star$ and then reduces drastically, with $F_\star \approx 0.3~ m a$. Moreover, the density wake in this case shows revivals separated by decreasing time intervals.
Read the full abstract

Earth and Planetary Astrophysics

Cooling-regulated gas accretion onto gap-opening planets

Ayumu Kuwahara, Michiel Lambrechts, Minghao Zhang et al. (+2 more)

Gas accretion onto forming planets controls the final masses of giant planets and provides observable signatures of ongoing formation. How this process depends on the cooling properties of these newly attracted gas remains poorly constrained. We present long-term, three-dimensional global hydrodynamical simulations to quantify gas accretion onto gap-opening planets in the mass range between 1 and 3 Jupiter masses. We systematically vary the cooling time, $\beta$, from near-isothermal ($\beta=10^{-2}$ in units of orbital time) to near-adiabatic ($\beta=10^{2}$), and follow the evolution until a quasi-steady state is reached. Our simulations show that the gas accretion rate decreases monotonically with increasing $\beta$, as $\dot{M}_{\rm acc}\propto\beta^{-0.18}$, reaching values at $\beta=10^2$ that are approximately an order of magnitude lower than locally isothermal predictions, largely independent of planet mass. The reduction in accretion is traced to thermodynamic restructuring of the circumplanetary region: inefficient cooling weakens shocks, narrows the accretion bands feeding the circumplanetary disk. Our results imply that thermodynamic effects should be taken into account when interpreting observed accretion rates of young planets, and may introduce systematic uncertainties in commonly used locally isothermal assumptions.
Read the full abstract

All authors Ayumu Kuwahara, Michiel Lambrechts, Minghao Zhang, Ruobing Dong, Hideko Nomura

High Energy Astrophysical Phenomena

Hyperon bulk viscosity effects in neutron-star inspirals

Elena M. Kantor, Mikhail E. Gusakov and Kirill Y. Kraav

The paper revisits the role of hyperon bulk viscosity during neutron-star inspiral. We find that hyperon bulk viscosity exerts only a minor influence on the gravitational-wave phase, causing a phase shift of $\sim 10^{-3}\,\rm rad$. This shift is too small to be detected by existing gravitational-wave observatories, and it seems unlikely that next-generation detectors will be able to isolate this effect from other factors. However, our analysis indicates that hyperon bulk viscosity can significantly heat the hyperon core of a neutron star, raising its temperature to approximately $5 \times 10^8\,\rm K$.
Read the full abstract

High Energy Astrophysical Phenomena

From Inspiral to Expansion: The Wake-Driven Torque on Binary Black Holes in Gaseous Medium

Jixuan Yang, Lile Wang, Xinyu Li et al. (+1 more)

Binary black holes (BBHs) in gaseous medium, such as active galactic nucleus (AGN) disks, are important gravitational-wave sources, yet the gas-driven torque that governs their orbital evolution remains to be fully understood. Most existing studies of BBHs in gas often approximate the net torque as the sum of independent dynamical friction (DF) forces exerted on each black hole by its own wake, neglecting the mutual gravitational coupling between the two wakes. We perform three-dimensional hydrodynamic simulations of circular BBHs in a uniform flow and find the torque is determined by the wake-wake interactions. The net torque is controlled by a single parameter $\eta \equiv v_g/v_o$, the ratio of the gas flow velocity to the binary orbital velocity. At small $\eta$, the wakes merge into a single overdense envelope and the time-averaged torque is negative; as $\eta$ increases, the wakes separate and the torque becomes positive. In AGN disks, capture-channel binaries in a disk model naturally produce $\eta$ in the positive-torque regime; the expansion timescale is comparable to or shorter than the disk lifetime, suggesting that gas-driven expansion can compete with gravitational-wave inspiral and suppress the capture-channel merger rate.
Read the full abstract

All authors Jixuan Yang, Lile Wang, Xinyu Li, Rixin Li

Cosmology and Nongalactic Astrophysics

Learning the averaged history of an inhomogeneous universe from its present day density field

Jonas Broe Bendtsen and Sofie Marie Koksbang

In this work, we investigate whether machine learning can be used to infer averaged cosmological quantities directly from the present day matter density distribution. Using an implementation of the simplified silent universe framework, we generate a dataset consisting of 92610 independent relativistic, simplified cosmological simulations spanning a range of initial conditions and average cosmological parameters. We train a convolutional neural network to take the present-time matter density field as input and predict the averaged matter, cosmological constant, curvature, and kinematical backreaction density parameters as well as the Hubble parameter, at both initial and present time. The network achieves coefficients of determination exceeding 0.9 for all predicted quantities, with error distributions showing that only a small fraction of predictions reach percent-level errors or above. Although our use of the simplified silent universe approximation precludes direct application of the trained model to observational data, the results provide a proof-of-principle that neural networks can successfully recover averaged cosmological properties, including backreaction, at different epochs, simply from the late-time matter distribution. More broadly, our findings demonstrate that information about the averaged cosmological history of a universe is encoded in its present-time density field and can be extracted using machine learning techniques. This opens the possibility of applying similar approaches to more realistic cosmological simulations and observational probes such as N-body simulations and weak-lensing maps, ultimately providing a new avenue for constraining the large-scale evolution of the Universe.
Read the full abstract

High Energy Astrophysical Phenomena

Rapid Variability and Broadband Spectral Modeling in the Flaring Activity of BL Lacertae

Xin Chang, Dingrong Xiong, Chenxu Liu et al. (+17 more)

We report a multi-wavelength study of two flaring episodes of the blazar BL Lacertae during MJD 60500-60800 (9 July 2024 - 5 May 2025). The source reached a daily-averaged $\gamma$-ray flux of $(1.03 \pm 0.05) \times 10^{-5} \, \mathrm{ph \, cm^{-2} \, s^{-1}}$ ($E > 100$ MeV) on MJD 60588 (5 October 2024). Using orbit-binned data from the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope, we identify a minimum flux halving timescale of $\tau = 1.33 \pm 0.29$ hr. This constrains the upper limit on the $\gamma$-ray emitting region size to $R \le 2.0 \times 10^{15}$ cm, as well as its distance from the central supermassive black hole to $R_\mathrm{H} \le 5.9 \times 10^{16}$ cm, assuming a Doppler factor of $\delta = 14.8$ derived from the spectral energy distribution (SED) modeling. We find tentative evidence for sub-minute $\gamma$-ray variability with a minimum doubling time of $0.7 \pm 0.2$ min ($p$-value = 0.03). This may originate from an extremely compact region with a size of $R \le 1.8 \times 10^{13}$ cm, suggesting that the emission arises from magnetohydrodynamic substructures, such as plasmoids within a magnetic reconnection zone. Spectral analysis reveals a significant “softer-when-brighter” trend ($r = 0.96, p = 4.5 \times 10^{-4}$) during the minute-scale flare peaks, indicating a complex interplay between particle acceleration and radiative cooling. The SED is reproduced using a one-zone leptonic model, in which synchrotron self-Compton (SSC) and external Compton (EC) scattering effectively account for the high-energy emissions. The reduced magnetic field strengths and hard electron injection spectral indices observed during the flaring states suggest enhanced particle acceleration efficiency, possibly associated with relativistic magnetic reconnection.
Read the full abstract

All authors Xin Chang, Dingrong Xiong, Chenxu Liu, Rui Xue, Tingfeng Yi, Jia Zhang, Yu Pan, Xingzhu Zou, Xinlei Chen, YeHao Cheng, Yuanpei Yang, Jinghua Zhang, Xiangkun Liu, Yuan Fang, Guowang Du, Tao Wang, Xufeng Zhu, Zhongxiang Wang, Sarira Sahu, Xiaowei Liu

Astrophysics of Galaxies

Dense Cores in the Vicinity of an HII Region

Ruofei Zhang, Xing Lu, Jingwen Wu et al. (+23 more)

Massive stars strongly influence their surroundings through radiative and mechanical feedback, but its effects on dense gas structures at sub-pc scales remain poorly constrained. We investigate how feedback from a newly formed massive star affects dense cores in the filamentary molecular cloud IRAS 18530+0215. We analyze ALMA Band 6 observations of 1.3 mm dust continuum and DCN, N$_2$D$^+$, and $^{13}$CS line emission, together with VLA K-band continuum and NH$_3$ observations. Dense cores are identified with astrodendro, and their temperatures, masses, velocity dispersions, and virial parameters are derived. The dynamical state of the ultra-compact H II region is examined through energy and pressure estimates. The H II region has a radius of $\sim$0.1 pc and an expansion velocity of $\sim$2.5 km s$^{-1}$, corresponding to a shell dynamical age of $\sim$0.06 Myr. DCN and $^{13}$CS cores are concentrated near the H II region, whereas N$_2$D$^+$ cores preferentially lie farther away. Core temperatures and velocity dispersions decrease with projected distance from the H II region. Virial parameters increase within the inner $\sim$0.3 pc but decline sharply beyond this scale, while core masses show no significant trend with distance. Strong star formation signatures are found at $\sim$0.2 pc, whereas more distant regions still host quiescent, cold dense cores. The compact H II region appears trapped or choked within $\sim$0.1 pc, while its feedback extends to at least $\sim$0.3 pc. Within this region, feedback enhances core velocity dispersions, gas temperatures, and virial parameters, with no evidence that it promotes the formation of more massive dense cores.
Read the full abstract

All authors Ruofei Zhang, Xing Lu, Jingwen Wu, Sihan Jiao, Hauyu Baobab Liu, Guang-Xing Li, Roberto Galván-Madrid, Aiyuan Yang, Siju Zhang, Shanghuo Li, Fengwei Xu, Xindi Tang, Yu Cheng, Weiyuan Zhang, Andrés E. Guzmán, Yuxin Lin, Yuhua Liu, Qizhou Zhang, Patricio Sanhueza, Ke Wang, Siyi Feng, Linjing Feng, Fangyuan Deng, Hao Ruan, Yuanzhen Xiong, Yuxiang Liu

Earth and Planetary Astrophysics

Modeling Carbon-Based Planets with MAGRATHEA

Robert Royer III, David R. Rice and Jason H. Steffen

A planet's interior structure not only influences its habitability, but it also provides information about the planet's formation and history. In the quest to characterize the variety of planetary discoveries, the most common practice is to rely on our understanding of the Earth's differentiated interior structure, with an iron core and silicate mantle as a starting point. However, recent work by arXiv:2408.07761 [astro-ph.EP] shows that these assumptions may not be true for all planetary systems as the composition of the condensed material can change drastically towards a carbon mantle instead of a silicate one, once the C/O ratio of the system approaches and exceeds approximately 0.9. In an effort to more fully understand the range of possible planetary bodies, we added the ability to model carbon-based mantles to our open-source planetary interior solver, MAGRATHEA. The newly added carbon mantle includes three relevant allotropes of carbon: graphite, diamond, and BC8 (body centered cubic with eight atoms per cell).
Read the full abstract

Solar and Stellar Astrophysics

Dynamically driven collapse of a thermally stable accretion disk in a nova-like system

Wentao Li, Jie Lin, Tinggui Wang et al. (+8 more)

While standard disk instability theory predicts thermal-steady, outburst-free disks in nova-like variables due to their high mass transfer rates, mass transfer variations are prevailingly invoked to explain the occasional flaring/fading phenomena observed in these systems. Here we report the observational evidence for a dynamical collapse of a thermally stable accretion disk, captured serendipitously by TESS during a fading episode of the VY Scl-type nova-like system MASTER OT J072703.91-631952.8, and traced by the emergence of an unusual negative superhump that evolves toward its orbital frequency. Simultaneously, the system underwent an anomalous eruptive event featuring a remarkably symmetric 45-day light-curve profile, indicative of a mild energy-release mechanism fundamentally distinct from documented eruptive events in cataclysmic variables. Notably, the concurrence of the eruption and the disk collapse is difficult to reconcile with the paradigm of mass transfer variations: the eruption implies enhanced mass transfer, whereas the disk collapse indicates a decline. The morphological disk evolution over $\sim 700$ days, characterized by both negative and positive superhumps, indicates a transition from a circular disk to an eccentric disk, followed by a tilted state and ultimately a minimal disk configuration. This evolutionary sequence provides evidence for a previously unrecognized dynamics-driven cycle operating in the thermally stable accretion disk of a VY Scl-type nova-like star.
Read the full abstract

All authors Wentao Li, Jie Lin, Tinggui Wang, Yongkang Sun, Chengyuan Wu, Heran Xiong, Krystian Iłkiewicz, Luca Casagrande, Michael Bessell, Ning Jiang, Chengyi Wang

Cosmology and Nongalactic Astrophysics

A Unified Tracer Analysis of DESI DR2 Baryon Acoustic Oscillations

N. Sanders, H. Seo, M. Rashkovetskyi et al. (+65 more)

We improve upon previous efforts to optimally combine overlapping galaxy samples in the DESI baryon acoustic oscillation analysis. By weighting each galaxy by its linear bias, overlapping galaxies are combined into a single, unified catalog, naturally avoiding double counting of cosmic volume and including all auto- and cross- information at the catalog level. Improvements over the previous effort include the addition of QSO out to $z=1.6$ to account for all overlapping DR2 tracers and redshift-dependent bias treatment to improve reconstruction. We report distance measurements using this unified tracer, and find them to be highly consistent with the baseline DR2 BAO analysis. We also test for tracer-dependent systematics within the DESI data, and find no evidence of tracer-dependent systematics within $0.8<z<1.6$. Finally, we take advantage of the unified tracer to rebin the analysis in redshift in order to more finely resolve the redshift-to-distance relation. Dynamical dark energy results on this finer redshift binning indicate that there is no missed feature in the expansion history in the redshifts $0.8<z<1.6$, and reproduces DESI's preference for an evolving dark energy equation of state.
Read the full abstract

All authors N. Sanders, H. Seo, M. Rashkovetskyi, U. Andrade, D. Valcin, E. Paillas, P. McDonald, J. Aguilar, S. Ahlen, O. Alves, E. Armengaud, A. Aviles, F. Beutler, D. Bianchi, D. Brooks, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, S. Cole, A. Cuceu, A. de la Macorra, Biprateep Dey, Z. Ding, P. Doel, S. Ferraro, A. Font-Ribera, D. Forero-Sánchez, J. E. Forero-Romero, E. Gaztañaga, Satya Gontcho A Gontcho, A. X. Gonzalez-Morales, G. Gutierrez, C. Hahn, H. K. Herrera-Alcantar, K. Honscheid, S. Juneau, T. Karim, D. Kirkby, A. Kremin, O. Lahav, C. Lamman, M. Landriau, L. Le Guillou, M. Manera, A. Meisner, R. Miquel, J. Moustakas, A. Muñoz-Gutiérrez, S. Nadathur, J. A. Newman, H. E. Noriega, N. Palanque-Delabrouille, W. J. Percival, F. Prada, I. Pérez-Ràfols, C. Ravoux, A. J. Ross, G. Rossi, L. Samushia, E. Sanchez, C. Saulder, D. Schlegel, M. Schubnell, J. Silber, M. Siudek, G. Tarlé, M. Vargas-Magaña, B. A. Weaver

Instrumentation and Methods for Astrophysics · Cross-listed from hep-ph

Future of Artificial Intelligence for Science in Japan 2024 Community Report

Yoshitaka Itow, Jia Liu, Hirokazu Maesaka et al. (+16 more)

This white paper summarizes scientific challenges and AI/ML research opportunities identified through the FAIRS Japan 2024 unconference process. The discussion focuses on three major physics domains: accelerator physics, cosmology and astrophysics, and neutrino physics. Although each domain has distinct scientific goals and experimental constraints, several common technical themes emerge: high-dimensional reconstruction, fast and accurate simulation, uncertainty propagation, simulation-to-data mismatch, anomaly detection, real-time decision-making, and shared infrastructure.
Read the full abstract

All authors Yoshitaka Itow, Jia Liu, Hirokazu Maesaka, Vinicius Mikuni, Nhat-Minh Nguyen, Hironao Miyatake, Atsushi J. Nishizawa, Patrick de Perio, Daniel Ratner, Kazuhiro Terao, Leander Thiele, Omar Alterkait, Francois Drielsma, Rocio Garcia, Masako Iwasaki, Ahsani Hafizhu Shali, Federica Tarsitano, Takahiro Terada, Junjie Xia

Solar and Stellar Astrophysics

Quiescent and Interacting Compact Remnant Binary Populations in the Space Telescope Era

Erin Eastep, Jan J. Eldridge and J. C. Bray

Since the discovery of Gaia BH1, interest in quiescent binary systems has grown. Two additional black-hole binaries have since been discovered alongside several candidate neutron-star binaries, yet the formation pathways of these systems remain poorly understood. Theoretical studies of these systems, together with the populations into which they evolve, such as X-ray binaries (XRBs), can help constrain their formation channels. We present predictions for the Galactic population of quiescent binary systems and their later evolutionary phases, including XRBs and post-interaction quiescent binaries. We combine stellar models from the Binary Population and Spectral Synthesis code (BPASS) with a synthetic Milky Way generated using Feedback In Realistic Environments (FIRE) simulations. We compare our results with existing observations to assess how well the synthetic population reproduces the observations. We find generally good agreement, with observed donor masses, compact-object masses, and orbital periods largely reproduced, except for the Gaia BH systems. Our results also suggest that some accreting compact remnants with low-mass donors and wide orbital periods may not be X-ray bright. Remaining discrepancies indicate that the BPASS initial parameter space should be expanded, the common-envelope evolution prescription revised, and magnetic braking included in low-mass donor models. We also predict the number of systems detectable in future Gaia data releases and by the Nancy Grace Roman Space Telescope (hereafter referred to as Roman). In particular, the Galactic Bulge Time Domain Survey is expected to yield a significant population of quiescent black-hole binaries. We predict that ~200-600 candidate quiescent binary systems will be detected in the Galactic Bulge Time Domain Survey, while ~2500-4000 candidate systems will be identified in Gaia DR5.
Read the full abstract

Thursday 27 August

High Energy Astrophysical Phenomena

Low-energy antinuclei measurements for background-free indirect dark matter searches and PBH signatures

Tsuguo Aramaki, Ilias Cholis, Philip von Doetinchem et al. (+7 more)

Light low-energy cosmic ray antinuclei constitute powerful probes for detecting dark matter or other sources of new physics in our Galaxy. This ASTRA Initiative Mission Concept reviews the theoretical motivation, highlighting the transformative potential and current experimental status. Motivated by exciting tentative results, it makes a clear case for the need for the next flagship mission to confirm or refute claims of cosmic antinuclei as signs of new physics.
Read the full abstract

All authors Tsuguo Aramaki, Ilias Cholis, Philip von Doetinchem, Fiorenza Donato, Priya Ghosh, Dan Hooper, Timothy Linden, Keith McBride, Stefano Profumo, Stephanie Wissel

Instrumentation and Methods for Astrophysics

JUG: JAX-based Unified pulsar timinG

Matthew T. Miles, Stephen R. Taylor, Matthew Bailes et al. (+14 more)

We present JUG (JAX-based Unified pulsar timinG), a JAX-based, fully independent pulsar timing package emphasising speed and ease of use, designed to confidently handle the increasingly large and complex pulsar timing array datasets that are being created in the pulsar timing field. JUG implements the entire pulsar timing pipeline itself, from data handling and clock corrections through to the timing model and fitting, without relying on other timing software. It enables Pythonic programming at the speed of compiled code, is GPU-capable, and can be operated via a Python API or an interactive GUI. A user of JUG can interactively explore data, fit timing models with complex stochastic noise, model deterministic signals such as continuous gravitational waves, and obtain accurate point estimates of the parameters of the stochastic processes present, thereby bridging frequentist timing and Bayesian noise analysis. JUG is faster than PINT by more than fifty times and is comparably fast to Tempo2, can handle millions of arrival times, agrees with PINT at the picosecond level, and can reliably recover known timing model and noise parameter values. In this paper we describe its design, performance, and validation, and demonstrate its advantages for pulsar timing data analysis.
Read the full abstract

All authors Matthew T. Miles, Stephen R. Taylor, Matthew Bailes, Aurelien Chalumeau, H. Thankful Cromartie, Kyle A. Gersbach, Rutger van Haasteren, Michael J. Keith, Nima Laal, Michael T. Lam, Kuo Liu, Aditya Parthasarathy, Scott M. Ransom, Daniel J. Reardon, Ryan M. Shannon, David C. Wright, Andrew Zic

Astrophysics of Galaxies

High Dispersion Measure Fast Radio Bursts from Young, Dust-Poor H II Regions

Eliza C. Diggins, Calvin Leung and Wenbin Lu

Several repeating fast radio burst (FRB) sources exhibit host-frame dispersion measure (DM) excesses far exceeding those expected from the interstellar medium in their host galaxies, implying a substantial contribution from the sub-kiloparsec-scale local environment. These excesses are often attributed to an expanding supernova remnant; however, this canonical interpretation does not explain their apparent preference for low-metallicity dwarf galaxy hosts, suggesting that the mechanism producing the local DM is coupled to the properties of the host galaxy. In this work, we investigate whether the H II regions generated by young, massive stellar associations can simultaneously explain these DM excesses and their preferential association with low-metallicity dwarf galaxies. We develop a semi-analytic model for FRB-emitting neutron stars in cluster-forming molecular clouds and perform idealized one-dimensional radiation-hydrodynamic simulations using time-dependent ionizing luminosities from pySTARBURST99. Our findings indicate that high-mass $(M_\star \sim 10^5-10^6\;M_\odot)$, young ($\lesssim 10\;{\rm Myr}$) stellar associations in low-metallicity environments (LMC-like or lower) can produce ${\rm DM}_{\rm HII} \gtrsim 10^3\;{\rm pc\;cm^{-3}}$ for FRBs that form early ($\lesssim 5$ Myr) and remain embedded in the ionized gas. The low-mass/metallicity preference arises from dust: absorption of the Lyman continuum suppresses the ionized column at higher metallicity, making extreme local DMs preferentially generated in dust-poor, low-metallicity environments. At cosmological distances, selection of a metal-rich, massive, and quiescent host galaxy sample where such extreme local contributions are likely to be suppressed may yield a cleaner FRB sample for using extragalactic DMs as tracers of the intervening intergalactic and circumgalactic media.
Read the full abstract

Instrumentation and Methods for Astrophysics

PolPy: A universal tool for X-ray/gamma-ray polarimetry using common data format standards

Sujay Mate, Hancheng Li, Utkarsh Pathak et al. (+4 more)

X-ray/gamma-ray polarimetry has been a rapidly developing field in recent years. However, the lack of standardized analysis frameworks has hindered cross-instrument comparisons and joint multi-mission studies. To address this, we present PolPy, a universal tool designed for polarimetry of high-energy astrophysical transients. Built as a plugin for the Multi-Mission Maximum Likelihood (3ML) framework, PolPy introduces standardized data and response formats, common reference frames based on IAU conventions for defining polarization angles, and a robust likelihood-based template-matching algorithm. The article briefly describes these standards and presents tests and results performed to validate the software using simulated gamma-ray burst (GRB) injections across two different mission mass models: POLAR and Daksha. Our verification tests demonstrate that PolPy accurately recovers the injected polarization angles and fractions without systematic bias. Furthermore, we show that performing simultaneous joint fits across multiple detectors and instruments can significantly improve parameter constraints compared to single-instrument analyses. By unifying data formats and methodology, PolPy enables reliable joint polarimetric analysis for current and future high-energy polarimeters.
Read the full abstract

All authors Sujay Mate, Hancheng Li, Utkarsh Pathak, Yashowardhan Rai, Nicolas De Angelis, Varun Bhalerao, Merlin Kole

Astrophysics of Galaxies

Near-Infrared Spectroscopy of Low-Redshift Palomar-Green Quasars and Composite Spectra

Amy Cavanaugh, Michael Brotherton, Jacob McLane et al. (+1 more)

The Palomar-Green (PG) Quasars have been used extensively to establish the spectral properties of luminous active galactic nuclei (AGNs) at low redshifts in many wavebands, but not the near-infrared (NIR). We present NIR spectra of all 87 PG quasars with z < 0.5 observed with the Apache Point Observatory 3.5-m telescope and the TripleSpec spectrograph in standard cross-dispersed mode. Here we employ our data set to construct a NIR composite quasar spectrum spanning rest-frame wavelengths from approximately 8800-22000 angstroms reaching a signal-to-noise ratio of several hundred in the H-band for 4 angstrom pixels. We characterize the continuum and emission-line properties of this composite. Additionally, we construct composite spectra as a function of luminosity and Eddington ratio, finding some differences, not entirely unexpected. In particular, more luminous composite NIR spectra grow increasingly red at long wavelengths, consistent with relatively stronger dust emission at higher temperatures.
Read the full abstract

All authors Amy Cavanaugh, Michael Brotherton, Jacob McLane, Richard Green

High Energy Astrophysical Phenomena

The Impact of Errors in the Shape Function of Rotating Neutron Stars in the Oblate Schwarzschild Approximation

John Ming Ngo, Charlee Amason and Sharon M. Morsink

Measurements of thermal X-ray flux emitted by neutron stars can constrain their radii and the equation of state of cold dense matter. Accurate flux modeling of rapidly rotating neutron stars requires an approximation for their oblate surface shape, usually given by a shape function. We examine the accuracy of five shape functions by comparing them against numerical relativistic stellar surfaces for seven representative equations of state over a broad range of compactness and spin. We characterize their errors in the surface radius and their derivative with respect to colatitude, test the extent of quasi-universality, and propagate these errors into the differential solid angle d{\Omega}, a geometric factor underlying the observed flux, under the oblate Schwarzschild approximation. We find that the shape functions produce distinct error patterns associated with their forms and parameter ranges. These errors remain structured in compactness, spin, and colatitude. Simple polynomial corrections substantially reduce the radius errors of several shape functions, with most corrected errors within approximately 0.5%. We show that errors in d{\Omega} may be understood as coming from an area contribution, which scales with twice the relative radius error, and a projection contribution controlled primarily by the error in the surface derivative. The projection contribution grows towards the limb, where small surface errors can produce much larger local errors in d{\Omega}. We provide a set of guidelines for the construction of improved quasi-universal shape functions that can be used in future radius estimation efforts. Our recommendations include that both the shape and its derivative be modeled accurately, that more accurate predictions of the polar radius be developed, and that the shape function be calibrated over the same compactness-spin range that is to be tested.
Read the full abstract

Astrophysics of Galaxies

Overmassive supermassive black holes in SDSS close galaxy pairs

Maggie C. Huber, Julia M. Comerford and Joseph Simon

Supermassive black holes (SMBHs) and their host-galaxies coevolve through channels such as hierarchical merging and AGN feedback, establishing tight scaling relations between SMBH masses and properties of the host-galaxy. These scaling relations, particularly those with host-galaxy bulge mass (M${_\text{b}}$) and stellar velocity dispersion ($\sigma$), are widely used to estimate SMBH masses. Galaxy mergers drive gas into the centers of galaxies, enhancing SMBH growth and stellar bulge growth but not necessarily in tandem, so it is necessary to determine if galaxies undergoing a merger still follow SMBH-host-galaxy scaling relations. In this study, we explore scaling relations in galaxy mergers by taking advantage of the well-established value-added catalogs in the Sloan Digital Sky Survey. These catalogs provide close galaxy pairs, AGN broad lines for SMBH mass measurements, and bulge-disk decomposed stellar masses. We compare SMBH mass estimates from M$_\bullet-$M${_\text{b}}$ to SMBH mass estimates from AGN broad lines and find that SMBH mass growth outpaces the bulge for AGN residing in the secondary (less massive) bulge and for smaller physical separations between galaxies in a pair. Using our full sample of close galaxy pairs, we find that M$_\bullet-$M${_\text{b}}$ and M$_\bullet- \sigma$ predict significantly different fractions of major and minor black hole mergers with M$_\bullet-$M${_\text{b}}$ predicting black hole mass ratios closer to 1:1 and $6-20\%$ more major mergers overall. These results have major implications, including for predictions of astrophysical gravitational-waves and high-redshift overmassive SMBHs.
Read the full abstract

High Energy Astrophysical Phenomena

Magnetospheric birefringence breaks the rotating-vector model and depolarises X-ray pulsars

Alexander A. Mushtukov, Sergey S. Tsygankov, Silvia Zane et al. (+4 more)

Accreting X-ray pulsars have long been expected to show large linear polarisation below the cyclotron resonance. In the standard picture, this follows from two facts: radiative opacities in the strongly magnetised neutron star environment depend strongly on polarisation mode, and the emission comes from compact polar regions. The emerging radiation should therefore be locally polarised and add almost coherently, with a polarisation angle tracing the projected magnetic field as in the rotating-vector model. Recent X-ray polarimetric observations, instead, show modest polarisation fractions and position angle swings that often depart from this simple picture. Here we show that propagation through the magnetosphere of an accreting neutron star can naturally produce both effects. Using general-relativistic ray tracing and Stokes-parameter transport, we find that photons emitted from compact polar regions acquire ray-dependent polarisation angle offsets in two ways. First, nearly field-aligned trajectories can cross local non-adiabatic pockets inside the QED birefringent region. Second, photons that intersect the plasma-loaded accretion flow accumulate additional phase shifts between plasma normal modes. Both effects broaden the distribution of polarisation angles across the ray bundle and reduce the hotspot-integrated linear polarisation, while also producing systematic deviations from rotating-vector-model expectations. Magnetospheric propagation should therefore be treated as part of the polarimetric transfer problem in accreting X-ray pulsars, and may help explain the observed low polarisation fractions and anomalous angle swings.
Read the full abstract

All authors Alexander A. Mushtukov, Sergey S. Tsygankov, Silvia Zane, Nabil Brice, Ruth M. E. Kelly, Roberto Taverna, Roberto Turolla

High Energy Astrophysical Phenomena

Spin-disc misalignment drives periodic accretion and pulse profile asymmetry in X-ray pulsars

A. A. Mushtukov, V. Ganesh, S. S. Tsygankov et al. (+2 more)

We investigate magnetospheric accretion in disc-fed X-ray pulsars assuming that the neutron star spin axis is not perpendicular to the disc plane. We focus on X-ray pulsars, where a geometrically thin disc is truncated far from the stellar surface the channelled part of the accretion flow is expected to be guided by the large-scale dipolar magnetic field after coupling to it near the disc-magnetosphere boundary. Using numerical simulations of plasma motion from the inner disc edge to the neutron star surface, we show that a finite inclination between the disc normal and the stellar spin axis leads to periodic modulation of the mass accretion rate onto the magnetic poles even for a steady mass supply through the disc. This purely geometrical effect arises because stellar rotation changes the orientation of the magnetosphere relative to the disc, producing phase-dependent mass loading of magnetic field lines. The amplitude and shape of the modulation are determined by the system geometry and by the ratio of the stellar spin period to the flow time through the magnetosphere. The resulting variability affects the structure and luminosity of emitting regions near the neutron star surface and leads to asymmetric X-ray pulse profiles. Even without intrinsic asymmetries of the emission regions, this mechanism breaks the time-reversal symmetry expected for stationary accretion and naturally contributes to the observed asymmetry of pulse profiles and phase-resolved spectral features. The effect may also be relevant for ULX pulsars, where intrinsic accretion rate modulation can help preserve strong pulsations in the presence of geometric beaming.
Read the full abstract

All authors A. A. Mushtukov, V. Ganesh, S. S. Tsygankov, S. Portegies Zwart, A. Palyam

Solar and Stellar Astrophysics

A 500 pc volume-limited sample of hot subluminous stars III. The short-period binary population

H. Dawson, M. Dorsch, J. Munday et al. (+16 more)

Hot subdwarf stars of spectral types O and B (sdO/B) in binaries form as products of substantial mass loss near the tip of the red giant branch and offer powerful constraints on binary-star evolutionary models. However, details of some formation channels remain missing. We present a comprehensive analysis of the short-period binary population in the 500 pc volume-limited sample of hot subdwarfs, enabled by precise Gaia DR3 parallaxes. Besides 45 known binaries within 500 pc, this work identifies 50 new single-lined radial-velocity variable systems, 34 with orbital solutions with periods between 0.07 and 22 days. We derive an overall short-period binary fraction of 34.7 (+2.8/-2.9)% for the full sample of 301 hot subdwarfs within 500 pc, corrected for orbital inclination and detection efficiency. The newly solved binaries predominantly occupy the 1-20 day range, a parameter space previously under-represented in the literature. We identify five new reflection-effect systems, three ellipsoidal modulation systems, a newly solved HW Vir system, and a new triple candidate. The sdB and sdO binary fractions are similar, while the sdOB class shows a fraction roughly half that of sdB/sdO stars, with periods over one day, suggesting a different evolutionary pathway. Underluminous hot subdwarfs below the canonical extreme horizontal branch show a binary fraction comparable to other sdB/sdO stars but host only white dwarf companions, which are more massive and mainly found at periods under one day. We provide the first volume-complete birthrate estimates for sub-populations, determining a Galactic merger rate of 2.5+-1.5x10^-5 yr^-1 for sdO/B binaries, explaining no more than ~12% of the eHe-sdO population. We also derive birthrates for two SN Ia progenitor channels, finding hot subdwarf binaries could account for up to 2.5 (+0.7/-0.5)% of the observed Galactic SN Ia rate.
Read the full abstract

All authors H. Dawson, M. Dorsch, J. Munday, S. Geier, F. Mattig, M. Pritzkuleit, D. Benitez-Palacios, M. Vučković, K. Deshmukh, A. Bhat, U. Heber, I. Pelisoli, R. Raddi, P. Fernandez-Schlosser, A. Durán-Reyes, E. Arancibia-Rojas, A. Bobrick, V. Schaffenroth, G. T. Jones

Instrumentation and Methods for Astrophysics · Cross-listed from eess.SY

Horizon-Dependent Tube MPC for Elliptical-Orbit Rendezvous Under Mass Uncertainty

Omer Burak Iskender and Keck-Voon Ling

A spacecraft closing on a target from three hundred kilometres to contact flies one guidance law across five orders of magnitude of range, and a controller that is provably safe at close range can lose that guarantee completely at long range while continuing to fly as though nothing were wrong. This paper derives the range at which the guarantee lapses and uses it as a design rule. The bound compares the prediction model's own linearisation error against the disturbance set the controller was built to reject, and needs only the sampling period, the orbit and that disturbance bound, so it can be evaluated before any simulation. On a Mars Sample Return approach it disqualifies the homing phase, where most of the propellant is spent, and clears the other two. Re-posing the disqualified phase in relative orbital elements restores the guarantee; re-posing a phase the rule already clears, in a frame two orders of magnitude more accurate, changes propellant by under a tenth of one per cent, and it is that second prediction that makes the rule falsifiable rather than descriptive. The constraint tightening also ties the prediction horizon to feasibility, so the horizon search limit becomes a mission parameter rather than a solver setting. Against a reimplementation of a published benchmark that reproduces its propellant to within one per cent, over five hundred dispersed Monte Carlo transfers per case on matched seeds, the controller saves 29% of the propellant on a circular target orbit and 40% on an eccentric one, docking inside the 0.20 m capture requirement on essentially every draw at a median miss near 5 cm. Two results run the other way: the saving is bought with time of flight and computation, and it comes from what the guarantee demanded of the terminal condition rather than from better disturbance rejection. Recursive feasibility and asymptotic stability are not claimed.
Read the full abstract

2608.27659 · 27 Aug · 19 pages, 12 figures, 8 tables. Accepted for presentation at the 77th International Astronautical Congress (IAC 2026), Antalya, Turkiye, 5-9 October 2026; paper IAC-26,C1,4,9,x108305. Accepted version (authors' preprint), not the version of record. Typeset in a plain article class rather than the IAF conference template

Abstract PDF $BNCH ×1

Instrumentation and Methods for Astrophysics

Toward a realistic test of black-hole movie correlations as a probe of extreme lensing

Barbora Bezděková and Shahar Hadar

In an optically thin setting, gravitational lensing around a black hole creates multiple delayed images of any localized source, appearing at different positions on the observer's screen. Even if particular sources cannot be singled out, this effect may be revealed in a black-hole movie (time-dependent image) by examining the two-point spatiotemporal correlation function of intensity fluctuations in the movie. Recently, the correlation function derived from a movie obtained by ray tracing of a general-relativistic magnetohydrodynamic simulation demonstrated the plausibility of this idea, even when the movie was blurred to angular resolutions corresponding to next-generation terrestrial very-long-baseline interferometric arrays. In this work, we further develop this method and study how some possible corruptions of the movie, motivated by real observational limitations, affect the ability to identify the extreme-lensing signatures by correlations. The effects are compared for both the original and blurred movies. We discuss the implications of our analysis for upcoming observations of M87* and Sgr A*, the two main targets of the Event Horizon Telescope and its future extensions. Our case study suggests that the types of data corruption we have explored, in themselves, do not preclude an identification of extreme-lensing effects in correlation measurements. We use our results to roughly estimate that, most importantly, for a successful correlation measurement, M87* may need to be continuously monitored for $\gtrsim 2$ years, while Sgr A* may need to be measured at a cadence $\gtrsim 1$ frame per minute. Finally, we outline a new analysis method that reduces the data in the complete correlation function—a large dataset defined on a five-dimensional configuration space—to a two-dimensional histogram that could be useful for black-hole parameter inference.
Read the full abstract

Solar and Stellar Astrophysics

Spectral Synthesis for Classical Ae Stars: H${\alpha}$ and CaII

R. Anusha and T. A. A. Sigut

Classical Ae stars (CAe) are main sequence, A-type stars with H${\alpha}$ emission but no signature of dust; they are thought to be the cool extension of the classical Be (CBe) stars to lower stellar masses. We investigate the use of H${\alpha}$ and the CaII IR triplet observed at low spectroscopic resolution ($R=1800$) as diagnostics for the circumstellar disks of the CAe stars. Using a sample of 159 CAe stars from the LAMOST survey (DR5), we find that circumstellar disk masses and star/disk viewing inclination angles can be reliably determined from low resolution H${\alpha}$ profiles. The CAe disk mass distribution agrees well with the previous range found for late-type CBe stars, although there is an apparent increase in disk mass for the latest CAe spectral sub-types (A3 and A4) which we attribute to selection effects. Over the spectral range A2 through B6, the average disk-to-star mass ratio for the CAe and CBe stars is consistent with $\log_{10}(M_{\rm disk}/M_*)=-9.1$. A minority of the CAe sample ($\sim$16%) also show emission in the CaII IR triplet. The observed emission is well within the range predicted by the disk models, and there is statistical evidence that the CaII emitters have systematically larger H${\alpha}$ disk masses.
Read the full abstract

Astrophysics of Galaxies

Cluster Ages to Reconstruct the Milky Way Assembly (CARMA). V. The chronological merger tree of the Milky Way

Fernando Aguado-Agelet, Chiara Zerbinati, Davide Massari et al. (+9 more)

We present a new age determination of 24 globular clusters (GCs) dynamically associated with the main accretion events experienced by the Milky Way (MW), as part of the Cluster Ages to Reconstruct the Milky Way Assembly (CARMA) project's effort to trace the Galaxy's assembly history. We used deep and homogeneous archival Hubble Space Telescope data, and applied the CARMA isochrone-fitting code to derive homogeneous estimates of age, metallicity, reddening, and distance modulus for systems dynamically associated with Gaia-Sausage-Enceladus (GSE), the Sagittarius dwarf galaxy (Sag), the Helmi streams (H99), and the Sequoia galaxy (Seq). These 24 new determinations are supplemented by 11 previously studied GSE clusters to construct the complete age-metallicity relation (AMR) of the GSE system. We find that each progenitor system describes a well-defined AMR, with a distinct slope and extent reflecting its individual star-formation efficiency and chemical enrichment history. By fitting analytical AMR models within a Markov chain Monte Carlo framework, we quantify the stellar mass and accretion time for each progenitor galaxy. This results in the first detailed merger tree obtained from strictly homogeneous chronological information, according to which the Low-energy-Kraken-Heracles (LKH) system is the first merger experienced by the MW that brought GCs in, followed by Sequoia, H99, GSE, and finally Sgr. The most significant events in terms of stellar mass are LKH, GSE, and Sgr, which together contribute a total of approximately 2.5 x 10^9 solar masses. This corresponds to more than 95% of the stellar mass accreted by the MW from mergers massive enough to host GCs.
Read the full abstract

All authors Fernando Aguado-Agelet, Chiara Zerbinati, Davide Massari, Matteo Monelli, Santi Cassisi, Cristiano Fanelli, Carme Gallart, Edoardo Ceccarelli, Yllari Kay González Koda, Tomás Ruiz-Lara, Sara Saracino, Maurizio Salaris

Astrophysics of Galaxies

Why is there cold gas inside the Local Bubble?

Tim-Eric Rathjen and Jeffrey L. Linsky

The Local Bubble (LB) extends approximately 100–200~pc from the Sun and contains $10^6$~K plasma and warm photoionised gas. Yet the Local Leo Cold Cloud (LLCC) and the Local Ribbon of Cold Clouds lie only 11–24~pc away, well inside the LB, at temperatures near 20~K. How such cold material forms and survives in this hotter environment remains uncertain, and explanations based on colliding warm clouds are difficult to test observationally. We investigate whether LLCC-like gas can assemble locally in a feedback-driven environment or must be imported as pre-existing cold material, and whether it can survive for several Myr. We use magnetohydrodynamic simulations from the SILCC Project including self-consistent star formation, stellar winds and supernovae, ionising and far-ultraviolet radiation, non-equilibrium chemistry, and cosmic-ray transport. We identify cold diffuse gas with $T<100$~K and $0<n_\mathrm{H_2}\leq2$~cm$^{-3}$ and use Lagrangian tracers to reconstruct its assembly, cooling, survival, and dispersal. The simulations produce cold diffuse gas in warm–hot, feedback-driven surroundings, where non-thermal pressure support is crucial. Cosmic-ray pressure remains nearly continuous across cold-gas interfaces and establishes a pressure floor, while magnetic pressure partly offsets the external thermal-pressure excess. The structures persist in the cold diffuse phase for a median of 2.5~Myr, and most of their material remains cold after they lose spatial coherence. The gas can assemble locally from nearby cold and thermally unstable warm material, developing filamentary or sheet-like morphologies. Cooling to $\sim20$~K is driven not by adiabatic expansion but by non-adiabatic cooling. The LLCC therefore need not have entered the LB as an already-formed cold cloud. LLCC-like gas can instead assemble and cool non-adiabatically in situ.
Read the full abstract