Peer Revenue An exchange of ideas

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Friday 17 April

Solar and Stellar Astrophysics

Revisiting the Galactic age-metallicity relation from wide white dwarf-main-sequence binaries

The age-metallicity relation (AMR) is a fundamental observational constraint for understanding the chemical evolution of the Galaxy. As reliable cosmochronometers, white dwarfs in binary systems with main sequence companions (WD+MS binaries) provide excellent laboratories to study this relation, since both components are expected to be coeval. We construct a sample of widely separated WD+MS binaries using data from the third data release of the Gaia mission in order to investigate the AMR of the Galactic disk. The sample is identified using photometric measurements and parallaxes of both components. White dwarf ages are derived by interpolating their Gaia absolute G magnitudes and BP-RP colours within state-of-the art white dwarf evolutionary sequences. We compile publicly available [Fe/H] abundances for the main sequence companions from the literature and combine them using different statistical approaches to obtain representative metallicity values for each system. We derive the AMR from several sub-samples of WD+MS that use independent measurements of [Fe/H] abundances and consistently find a large dispersion in [Fe/H] at all ages. This behaviour agrees with previous determinations of the AMR based on both WD+MS binaries and samples of isolated stars. Our results reinforce the observational evidence that the AMR in the Galactic disk exhibits substantial intrinsic scatter, likely reflecting the combined effects of multiple mechanisms such as radial migration, inhomogeneous chemical enrichment, and variations in the star formation history.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Inferring Halo Mass and Scale Radius of Galaxy Clusters Using Convolutional Neural Networks and Uchuu-UniverseMachine Catalogs

We investigate the ability of machine learning to infer the virial mass ($M_{\rm vir}$) and the scale radius ($r_{\rm s}$) of galaxy clusters from their observables. Using the Uchuu–UniverseMachine galaxy catalog at $z=0.093$, we generate mock cluster observations that include interlopers, and we encode each cluster as an image representing the two-dimensional joint probability distribution of member galaxies' projected position and line-of-sight velocity. We train two architectures: a baseline convolutional neural network (CNNb) following a previous approach, and an extended model (CNNr) that appends richness as an additional scalar input. We further compare the performance of networks trained on the all cluster sample and on a dynamically relaxed subsample. Across the test ranges $10^{13.7}\leq M_{\rm vir}\leq10^{15.3}$ Msun/h and $10^{1.7}\leq r_{\rm s}\leq10^{2.7}$ kpc/h, all configurations yield nearly unbiased absolute median residuals (within 0.01 dex). For the halo mass, adding richness narrows the residual distribution, reducing the standard deviation from 0.133 to 0.122 dex for the all sample, and from 0.124 to 0.111 dex for the relaxed sample. For the scale radius, restricting the training to relaxed clusters improves the performance more than adding richness. The standard deviation decreases from 0.180 to 0.154 dex for CNNb and from 0.175 to 0.148 dex for CNNr, while the inclusion of richness yields only a modest improvement of 0.005 dex. These results demonstrate that machine learning is a powerful tool to infer the mass and internal mass distribution of clusters, providing a new window for cosmological inferences and understanding galaxy formation processes.
Read the full abstract

All authors Hirobumi Tominaga, Asuka Nakamura, Tomoaki Ishiyama, Mohamed H. Abdullah

Astrophysics of Galaxies

The Reaction between Atomic Carbon and Molecular Nitrogen as a Source of Cyanamide and Carbodiimide on Interstellar Ices

Reactions occurring on the ice-covered surfaces of interstellar dust grains are considered to be among the most important sources of complex species in the interstellar medium. Despite this, molecules such as cyanamide, NH2CN, are largely underpredicted by current astrochemical models suggesting that the network of reactions currently used to describe this species and its tautomer carbodiimide, HNCNH, are incomplete. Here, we performed a theoretical investigation of the reaction of ground state atomic carbon C(3P) with molecular nitrogen N2 in both the gas-phase and on the surface of amorphous solid water (ASW) clusters to examine its potential importance in the formation of NH2CN and HNCNH. We show that the reaction of gas-phase C-atoms with N2 molecules already present on the ASW surface results in the barrierless formation of CNN. Following exothermic hydrogenation reactions, the N-N bond of the C-N-N bearing intermediates is broken allowing the formation of molecules with N-C-N backbones through cyclic intermediates over low barriers. To test the importance of these processes to NH2CN and HNCNH formation, these reactions were included in a three-phase astrochemical model of low-mass protostellar evolution employing a reaction network that was updated to better describe the formation and destruction pathways of related small nitrogen bearing molecules. These simulations demonstrate that the ice surface reaction between C and N2 represents by far the dominant source of NH2CN and HNCNH in protostellar environments and in dense clouds.
Read the full abstract

Earth and Planetary Astrophysics

Beyond the mass-radius plane: Integrated radiative-convective and interior structure simulations of the exoplanet continuum

Static structure models, which map mass-radius constraints to bulk planet composition, are frequently used to categorise exoplanets due to their computational efficiency and the high-level insight they offer into planetary properties. However, static structure models typically have simplified atmospheric treatments, which may introduce systematic biases when interpreting the structures - and therefore the climates - of sub-Neptunes and super-Earths. We present a framework for recovering exoplanet properties using static structure models that accounts for necessary physical-chemical complexity in their atmospheres. We produce a comprehensive library of 504,000 exoplanet simulations that unify deep planetary interior structure with radiative-convective-chemical climate calculations. From these models we demonstrate that a planet's envelope mass fraction - a critical parameter to infer - is frequently degenerate with its instellation flux and atmospheric metallicity, and sensitive to the treatment of gravitational acceleration at the mbar level. Such uncertainties have significant implications for inferring planetary processes, as our modelling shows that habitable-zone sub-Neptunes readily host supercritical surfaces or deep magma oceans, despite their temperate irradiation regime. To marginalise over these uncertainties, we introduce a Bayesian retrieval tool that uses our library of self-consistent models. By applying this Bayesian approach to case-studies of pi Men c and TOI-421 b, we show that robust physical interpretations are achievable through whole-planet mass-radius retrievals. While new data from JWST, Ariel, and PLATO will expand our observational horizon, physically-consistent modelling provides the means to transition from categorical interpretations toward a comprehensive picture of the exoplanet continuum.
Read the full abstract

All authors Harrison Nicholls, Oliver Shorttle, Tim Lichtenberg, Flavia Pascal

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

Gravitational-wave astronomy requires population-informed parameter estimation

Gravitational-wave events are interpreted in terms of Bayesian posteriors for their source properties inferred under unphysical reference priors. Though these parameter estimates are important intermediate data products for downstream analyses, we demonstrate that they are generically biased and therefore should not be used for astrophysical interpretation directly, as is common. Hierarchical parameter estimation is the solution, as joint analysis of the entire catalog of observations reduces statistical uncertainties and actually informs the correct prior, with population-informed event parameters now appropriate for astrophysical interpretation. As an example, we show how the most extreme measurements from a catalog can be derived and used to identify exceptional events from previous and ongoing observing runs, pointing out they are more informative about the population than any individual event. Using LIGO-Virgo-KAGRA data, we thus demonstrate that population inference is not optional to interpret gravitational-wave observations.
Read the full abstract

Solar and Stellar Astrophysics

PISP: Projected-Space Inference of Stellar Parameters

To improve the accuracy and efficiency of high-dimensional stellar parameter inference in large spectroscopic datasets, we propose a projection-assisted parameter-inference framework – Projected-Space Inference of Stellar Parameters (PISP). PISP constructs an orthonormal basis and optimizes in the projected space, reducing the impact of parameter correlations on inference. The basis is constructed using either principal component analysis (PCA) or the active-subspace (AS) method and is combined with two inference strategies – Non-L1, which optimizes the projection coefficients for a user-specified projected dimensionality, and L1, which introduces L1 regularization in the full projected space to adaptively select projection directions – yielding four strategies: PCA-Non-L1, AS-Non-L1, PCA-L1, and AS-L1. For different computational scenarios, we implement two versions: PISP-CurveFit for fast single-spectrum inference and PISP-Adam for large-scale GPU-parallel inference. Using a fully connected neural network and a residual network as spectral emulators, we evaluate PISP on Kurucz synthetic spectra and on $722{,}896$ APOGEE DR$17$ observed spectra. Compared to the baseline strategy, PISP improves inference accuracy for multiple parameters across all emulator-optimizer combinations. In synthetic data, PCA-L1 performs best, reducing the standard deviation of differences ($\sigma(\Delta)$) by at least $0.01$ dex for $12$ of $20$ elemental abundances, with [N/H], [O/H], [Na/H], [Co/H], [P/H], [V/H], [Cu/H] showing $0.05$–$0.72$ dex reductions. In observed data, PCA-Non-L1 reduces $\sigma(\Delta)$ by $>30$ K for effective temperature and by at least $0.01$ dex for $9$ of $17$ elemental abundances, with [O/H], [Na/H], [V/H] showing $0.05$–$0.20$ dex reductions, while achieving a $\sim$$4\times$ efficiency gain, slightly outperforming PCA-L1.
Read the full abstract

All authors Jun-Chao Liang, Yin-Bi Li, A-Li Luo, Shuo Li, Xiao-Xiao Ma, Hai-Ling Lu, Shu-Guo Ma, Ming-Hui Jia, Shuo Ye, Hao Zeng, Ke-Fei Wu, Zhi-Hua Zhong, Xiao Kong, Li-Li Wang, Hugh R. A. Jones

Earth and Planetary Astrophysics

Benchmarking Two Chemical Networks used in General Circulation Models of Hot Jupiters

Chemical kinetics is becoming an increasingly vital component of hot Jupiter general circulation models (GCMs). Here we simulate the hot Jupiter WASP-96b using two chemical networks, a reduced chemical network frequently used in the GCM literature (which we refer to as V19) and a more recent effective network making use of tables of net reactions (MiniCHEM), coupled to the same GCM in order to provide a robust benchmark. We find a numerical escape criterion used by the Unified Model chemical kinetics solver to stop integration for the duration of the chemical timestep, independent of the chemical network, results in artificial quenching, overestimating of HCN, CH$_4$, and NH$_3$ abundances by factors of 1.5 to 3. With this criterion disabled, agreement between the two networks is improved, except for HCN and NH$_3$, where different reaction rates and included species results in lower abundances in the V19 network. While many rates differ between the networks, the lower quenched NH$_3$ abundances in the V19 simulations are, in particular, due to the choice of NH$_2$ + NH$_3$ $\rightarrow$ N$_2$H$_3$ + H$_2$ reaction rate, which is poorly constrained in the literature. This reaction also impacts the quenching of HCN, which is additionally affected by the lack of CH$_2$NH$_2$ in the V19 network. While there are reasons to favour the MiniCHEM HCN and NH$_3$ abundances, ultimately, improved experimental and theoretical determination of reaction rates are needed to address the uncertainties and better characterize the quenching behaviour.
Read the full abstract

All authors D. A. Christie, M. Zamyatina, E. Hébrard, T. M. Evans-Soma, N. J. Mayne, E. K. H. Lee, S.-M. Tsai, D. E. Sergeev, R. Veillet, K. Kohary

Earth and Planetary Astrophysics

A gem system with a lava world and a habitable zone sub-Neptune orbiting TOI-1752

The Transiting Exoplanet Survey Satellite (TESS) has delivered a large number of transiting planet candidates around nearby stars by identifying periodic decreases in stellar brightness. Establishing the planetary nature of these signals and determining their fundamental properties is a necessary step toward detailed studies of their internal structure, atmospheres, and formation pathways. In this work, we investigate the planetary nature of the TOI-1752 system (M1 V, $103.02\pm0.34$ pc), which hosts two TESS candidates: TOI-1752 b, a short-period object consistent with a lava-world scenario, and TOI-1752 c, a sub-Neptune-size planet candidate located in the optimistic habitable zone. We obtained ground-based multi-color photometric follow-up observations of TOI-1752, which we combined with TESS photometry to assess the nature of both signals. We performed a formal statistical validation using the TRICERATOPS framework, while independently vetting the candidates with the neural-network-based classifier WATSON-Net, which provides a machine-learning assessment of their planetary likelihood based on light-curve morphology, centroid diagnostics, and auxiliary vetting features. We validate TOI-1752 b as a bona fide planet with a radius of $1.69\pm0.07 R_{\oplus}$ and an orbital period of $0.935186^{+0.000001}_{-0.000002}$ days, and TOI-1752 c with a radius of $2.29^{+0.13}_{-0.14} R_{\oplus}$ and an orbital period of $32.7144\pm0.0004$ days. The combined analysis confirms TOI-1752 as a new planetary system, places TOI-1752 c within the optimistic habitable zone of its host star, and identifies TOI-1752 b as a promising target for atmospheric characterization, with an estimated emission spectroscopy metric (ESM) of up to $\sim8$.
Read the full abstract

All authors A. Peláez-Torres, F. J. Pozuelos, G. Morello, M. Dévora-Pajares, K. Barkaoui, L. Gkouvelis, E. Pallé, K. A. Collins, B. V. Rackham, S. Geraldía-González, M. Centenera-Merino, R. Varas, E. Esparza-Borges, Z. Parlapani, J. Flores, J. Aceituno, P. J. Amado, A. Burdanov, Y. Calatayud-Borras, D. R. Ciardi, B.-O. Demory, T. Gan, S. Giacalone, M. Gillon, Y. Gómez Maqueo Chew, K. Kawauchi, A. Khandelwal, J. Korth, M. Lendl, J. P. de Leon, J. Livingston, N. Morales, F. Murgas, N. Narita, J. L. Ortiz, H. Parviainen, M. Pichardo Marcano, I. Plauchu-Frayn, D. Queloz, D. Rapetti, J. Saito, A. Sánchez-López, A. B. Savel, R. P. Schwarz, U. Schroffenegger, M. Serra-Ricart, C. Stockdale, A. H. M. J. Triaud, J. de Wit, F. Zong Lang

Earth and Planetary Astrophysics

Asteroid modelling by starlight diffraction: The shape of Dimorphos, the satellite of (65803) Didymos

The DART spacecraft impacted Dimorphos, the satellite of (65803) Didymos, in September 2022. Evidence of crater formation and possible global reshaping has been obtained indirectly from spacecraft and ground-based data. Since the impact, several stellar occultations by Didymos have been observed, but only one in particular, on January 21, 2023, can provide useful constraints on the size and shape of Dimorphos. We modelled the diffraction signatures recorded on multiple occultation chords to constrain the projected shape and size of Dimorphos, assuming an ellipsoidal model. This is the first time diffraction observed simultaneously on several chords of a single event has been used for such a purpose. The projected shape at the epoch of the event is well constrained and consistent with recent post-DART determinations. When extended to a full three-dimensional ellipsoidal solution, the result remains compatible with previous studies, suggesting an equatorially elongated post-impact shape.
Read the full abstract

All authors P. Tanga, K. Tsiganis, D. Souami, R. Anderson, E. Barbotin, A. Cazaux, F. Colas, J. Hanuš, F. Marchis, J-L. Dauvergne, G. Langin, A. Leroy, B. Lott, A. Manna, L. Rousselot, A. Siakas, S. Sposetti, Ch. Vigna, F. Weber, A. Wünche

High Energy Astrophysical Phenomena

Induced Scattering of Strong Waves in Pair Plasmas

We study induced (stimulated) scattering of linearly polarized, strong electromagnetic waves in pair plasmas, which is crucial for understanding the propagation of fast radio bursts (FRBs). Magnetars are the most promising progenitors of FRBs, and FRBs propagate through the magnetar wind and successfully escape before being significantly scattered. We revisit the steady-state solution of linearly polarized electromagnetic waves in pair plasmas with arbitrary amplitude, and demonstrate that the nonlinearity is characterized by the nonlinearity parameter $a_0\omega_{pe}/\omega_0$ rather than the dimensionless amplitude $a_0$, where $\omega_{pe}$ is the electron plasma frequency and $\omega_0$ is the wave frequency. We follow the time evolution of the steady-state solution for the linear regime $a_0\omega_{pe}/\omega_0 \ll 1$ by performing one-dimensional particle-in-cell simulations, and show that the conventional linear analysis of induced scattering assuming $a_0 \ll 1$ is applicable even for $a_0 > 1$ when the Lorentz boost due to the plasma motion in the incident wave is considered. The saturation level is controlled by $a_0\omega_0/\omega_{pe}$, which corresponds to the ratio of the wave energy to the plasma energy, and the incident wave is hardly scattered for $a_0\omega_0/\omega_{pe} \gg 1$. We discuss the application of our results to FRBs.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Inverting Fisher biases for fast systematics exploration

Upcoming cosmological surveys will achieve increasingly precise constraints in cosmological parameter estimation. To guarantee the robustness of cosmological analyses, it is essential to account for and model systematic effects that can bias cosmological constraints, shifting the best fit parameters away from their fiducial values. It is possible to approximately infer the biases that un-modelled systematic effects might introduce in cosmological parameter estimation by means of the Fisher matrix formalism. In this paper, we introduce a new application of this formalism, where by inverting the process, we investigate whether a specific missing or mis-modelled systematic effect can explain away a given tension between two different probes or experiments. We showcase the proposed methodology by examining two representative systematics: galaxy intrinsic alignments and baryonic feedback. As the method is agnostic to the systematic effect and can be applied to a wider range of scenarios, we discuss more possible future applications. While the proposed approach is accurate in the limit of small offsets in the cosmological parameters, where the likelihood can be considered linear in both the cosmological parameters and the systematics, in practice, the region of validity depends on the systematic effect. In general, even beyond this region, the approach still provides a useful test that helps indicate the magnitude and direction of potential biases from systematic effects in data.
Read the full abstract

Instrumentation and Methods for Astrophysics · Cross-listed from gr-qc

Mechanical Long Baseline Differential Gradiometers as Low Frequency Gravitational Wave Detectors

We present a new differential mechanical gradiometer for the detection of low-frequency Gravitational Waves. The frequency range is 0.05 to 1 Hz, a frequency gap not covered either by future space-based detectors such as LISA or by ground-based observatories such as Einstein Telescope or Cosmic Explorer. The proposed detection principle is similar to antennas based on torsion pendulums but solves the problem of physical confinement of these antennas by operating vertically and by having a counterweight at one end of each bar and a mass suspended from a long wire at the other. With this configuration, we enlarge the gravitational force acting on the system without changing the moment of inertia of the system, so that we move from a signal $\Delta \theta$ of the order of $\Delta \theta = h$, where h is the amplitude of the gravitational wave, to a signal of the order $\Delta \theta = h\frac{L}{D}$, where D is the length of the arm and L is the length of the wire suspending the test mass. This configuration is a further evolution of the recent development of tiltmeters and balances with double suspended arms and interferometric read-out, where the main working principles are already tested. The expected sensitivity will be discussed with respect to the proposed parameters and the present technology.
Read the full abstract

All authors Enrico Calloni, Annalisa Allocca, Antonino Chiummo, Rosario De Rosa, Luciano Errico, Marina Esposito, Edoardo Imparato, Bruno Mantice, Luigi Rosa, Paolo Ruggi, Alessandra Ruggiero, Valeria Sequino, Daniela Stornaiuolo, Vittorio Tortorella, Lucia Trozzo

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

Geometrically Regular Black Holes with Hedgehog Scalar Hair

We study a simple theory based on general relativity, minimally coupled to a constrained scalar triplet and to an auxiliary non-propagating three-form sector. Within a spherically symmetric hedgehog ansatz, the theory admits a continuous exact family of asymptotically flat geometrically regular black holes. For a simple choice of kinetic function, the solutions possess a de Sitter core and approach Schwarzschild with the first correction appearing only at order $r^{-4}$. We analyse their horizon structure, thermodynamics, and main strong-field properties. The black holes carry topological scalar hair and a continuous secondary parameter, but no scalar charge. The regularity established here is geometric: the curvature invariants remain finite, although the matter sector is not completely smooth at the centre.
Read the full abstract

High Energy Astrophysical Phenomena

A universal relationship between the variability timescale and black hole mass in black hole jetted and non-jetted accreting systems

A long-term variability study spanning a range of black hole mass systems, from microquasars hosting stellar-mass black holes to active galactic nuclei (AGNs) harboring supermassive black holes, provides new insights into the physics of relativistic jets. In this work, we investigate the optical variability of both jetted and nonjetted AGNs. We apply a stochastic process known as the Damped Random Walk (DRW) to model light curves from the Zwicky Transient Facility (ZTF) DR23. Our results show that the mass-scaled characteristic timescales across the black hole mass exhibit a linear relationship with a slope of 0.35-0.50. This analysis confirms a previously observed correlation between the damping timescales and black hole mass and extends it by incorporating 125 newly identified non-jetted AGNs selected from the Burst Alert Telescope (BAT) AGN catalogue. The derived slope of the relation between the damping timescales and black hole mass aligns with recent theoretical predictions, supporting the presence of a universal accretion mechanism in AGNs across different mass scales. Furthermore, our findings suggest a novel implication: the properties and production mechanisms of relativistic jets may be largely independent of black hole mass.
Read the full abstract

All authors Yongyun Chen, Qiusheng Gu, Junhui Fan, Dingrong Xiong, Xiaoling Yu, Xiaogu Zhong, Xiaotong Guo

Astrophysics of Galaxies

A New Cloud-Cloud Collision Source N68 toward the G35 Molecular Cloud Complex

Bubble N68 in the G35 complex shows clear cloud-cloud collision (CCC) signatures. Its semi-ring-like morphology harbors many significant massive star formation tracers: 6 HII regions, 4 6.7 GHz masers, 5 Midcourse Space Experiment sources, 9 radio peaks, and nearly 10 O/B-type stars. We also identified 163 young stellar objects (45 Class I, 5 Flat, 113 Class II), indicating active star formation toward N68. Our molecular study with CO reveals two distinct molecular clouds (N68a: 47-56 km s$^{-1}$; N68b: 56-64 km s$^{-1}$), with broad bridge features and complementary distributions at their borders, indicating an ongoing CCC. Star formation in N68 is collectively driven by collect-and-collapse (CC), radiation-driven implosion (RDI), and CCC mechanisms. However, compared with the CC and RDI mechanisms, the CCC mechanism does not enhance the star formation efficiency; instead, it tends to trigger the formation of massive stars. N68, along with bubbles N65 and N61, constructs a $\sim100$ pc scale CCC system in the G35 complex.
Read the full abstract

Solar and Stellar Astrophysics

Radiation-driven stellar winds at the fast-slow transition: new hydrodynamic solutions

Radiation-driven winds of massive stars can be described within the modified CAK theory, which parametrises the radiation force through three key quantities: $\alpha$, $\delta$, and $k$. Different combinations of these parameters, together with rotation, result in three types of stationary solutions, namely fast (or classical), $\delta$-slow, and $\Omega$-slow solutions. The primary objective of this work is to model radiation-driven winds inside the gap region between the fast and $\delta$-slow regimes, where stationary solutions have proven elusive. In addition, we compute synthetic line profiles of H I, He I, and Si IV to illustrate the morphology of different wind regimes. We employ the time-dependent hydrodynamic code ZEUS-3D, capable of obtaining stationary solutions by progressing through an initial solution. Then we compute the line profiles solving the transfer equation for an expanding atmosphere, assuming spherical symmetry in the comoving frame, under non-local thermodynamic equilibrium (NLTE) conditions. We found new stationary solutions in the gap region, alongside their corresponding line profiles, for a typical B supergiant star model. In this model, the new solutions are stable, and some of them present a kink in the velocity profile at a fixed distance from the star, depending on the $\delta$ value. Perturbations in the wind ionisation may trigger transitions between different hydrodynamic regimes and offer a plausible explanation for structured and variable winds. A systematic investigation of these effects will be the subject of future work. Furthermore, we investigate the resulting line profiles from different hydrodynamic solutions and compare them with those predicted by a velocity profile given by a $\beta$-law using the same global wind parameters.
Read the full abstract

All authors M.C. Fernandez, R.O.J. Venero, L.S. Cidale, I. Araya, M. Curé

High Energy Astrophysical Phenomena

Circular polarization images of Sgr A* for different magnetic field geometries

Sgr A* exhibits a persistent negative circular polarization (CP) at 230\,GHz, offering a powerful probe of the magnetic field geometry in its accretion flow. Using a stationary semi-analytic radiatively inefficient accretion flow (RIAF) model in Kerr spacetime with polarized radiative transfer, we systematically analyze CP images for six poloidal magnetic field configurations across varying black hole spins, inclinations, and field polarities. We find that CP production is dominated by Faraday conversion in radial, parabolic, quadrupole, and combined geometries, but by intrinsic emission in dipole and vertical fields. The radial and parabolic configurations produce the polarity-invariant net CP, while dipole and vertical fields yield the polarity-sensitive one. As the accretion disk is prograde with respect to the black hole spin, the CP production across all six field geometries is found to be lower at high spin case, while the situation is more complicated in the retrograde case. Moreover, the net CP observed from edge-on views $V_{\rm net} \approx 0$ except for the quadrupole geometry. Comparing with ALMA data, the reversed-field model is excluded at high inclinations and then the magnetic field geometry of Sgr A* is constrained.
Read the full abstract

Astrophysics of Galaxies

Distances to molecular clouds in the Galactic longitude l=10-20 deg from the MWISP 12CO 1-0 survey

We present distances to 56 molecular clouds within $10\degr \leq l \leq 20\degr$ and $|b| \leq 5.25\degr$ from the Milky Way Imaging Scroll Painting (MWISP) $^{12}$CO survey, 47 of which are first-time determinations. The molecular clouds were identified using the DBSCAN algorithm, and their distances were measured with the model-calibrated color-distance method using $J-K{_s}$ colors and the distances provided by 2MASS and Gaia EDR3. The distances range from $\sim$275 pc to $\sim$2118 pc. We also derived the physical properties of molecular clouds and found a moderate correlation between the dust extinction and the $^{12}$CO integrated intensity.
Read the full abstract

All authors Juan Mei, Zhiwei Chen, Min Fang, Miaomiao Zhang, Shiyu Zhang, Zhibo Jiang

Astrophysics of Galaxies

The Impact of the Bar on Dense Gas and Star Formation in M83

Although the relationship between molecular gas content and the star formation rate (SFR) has been extensively studied in nearby galaxies, it remains controversial whether the star formation efficiency (SFE) depends on galactic structure. In particular, whether the SFE is suppressed in the bar region compared with other structures, and the physical origin of this suppression, remain poorly understood. In this study, we investigate variations in the SFE and its physical drivers in the bar region of the nearby spiral galaxy M83, using multi-wavelength observations toward the bar and spiral arm regions on a scale of 200 pc observed with ALMA. From the molecular gas surface density derived from $^{12}$CO($J=2-1$) and $^{13}$CO($J=1-0$), the dense molecular gas surface density derived from HCN ($J = 1-0$), and the star formation rate surface density determined from extinction-corrected H$\alpha$, we find that the SFEs in the bar region are roughly a factor of two lower than those in the spiral arm, indicating that the SFE is systematically suppressed in the bar. Moreover, we find that the SFEs of dense gas are lower in the bar than in the arm by a factor of about 0.35. These results suggest that not only the efficiency of converting bulk molecular gas into stars reduced in the bar, but the efficiency of star formation from dense gas is also lower. In addition, the CO line widths are systematically larger in the bar region and exhibit a negative correlation with both the SFE and the dense-gas SFE, consistent with the interpretation that enhanced turbulent motions hinder star formation. Although the analysis is limited to small regions of M83, our results suggest that the suppression of the SFE is related to large-scale dynamical effects on the molecular gas, such as strong shocks induced by cloud-cloud collisions and/or shear, both driven by non-circular motions in the bar.
Read the full abstract

Thursday 16 April

Astrophysics of Galaxies

Predicting Supermassive Black Hole-Host Mass Offsets from Broadband Photometry Across Cosmological Simulations with Forecasts for LSST

The possibility of over-massive black holes suggested by James Webb Space Telescope photometric discoveries of 'little red dots', may disfavor light supermassive black hole (SMBH) seeds. However, what should constitute the mass (range) of 'heavy' seeds remains relatively unconstrained. Moreover, Vera C Rubin Observatory's Legacy Survey of Space and Time will photometrically characterize galaxies without direct black hole mass measurements. We forward-model the SIMBA, IllustrisTNG, and EAGLE cosmological simulations into the photometric bands of LSST to train an ensemble machine learning classifier. Our framework achieves $91\%$–$94\%$ accuracy across SIMBA and IllustrisTNG in distinguishing between over-massive and under-massive SMBH growth regimes under LSST magnitude limits, using only broadband photometry. Furthermore, cross-simulation transfer experiments (training on one cosmological simulation and evaluating on another using rank-normalized features) achieve $83\%$–$89\%$ accuracy. This suggests the relative photometric ordering of growth regimes is largely preserved even across fundamentally different sub-grid SMBH feedback prescriptions. Signal decomposition shows our classification is driven by host galaxy colors ($82\%$–$87\%$ accuracy) and, relatedly, the accretion-state's spectral energy distribution shape as opposed to an inversion of our forward model's analytical luminosity prescription. Given that the evaluated simulations employ heavy seed prescriptions ($\geq 10^{4}~M_\odot$), our methodology establishes a validated baseline for classifying post-seeding growth regimes.
Read the full abstract

Cosmology and Nongalactic Astrophysics

AI–Assisted Exploration: DHOST Theories without Quantum Ghosts

We ask whether a local symmetry can organize both classical degeneracy and perturbative stability in DHOST theories. An AI-assisted search finds a candidate in a constant-disformal image of Einstein gravity. We identify it as a field-dependent diffeomorphism plus a vertical translation of a spectator scalar and extend it to every regular first-derivative map $\widetilde g_{\mu\nu}=C(\phi,X)g_{\mu\nu}-D(\phi,X)\nabla_\mu\phi\nabla_\nu\phi$, which pulls the spectator translation back to an exact local gauge redundancy. We derive its closed generator and show that $\mathcal J=C-XC_X+X^2D_X$ controls kinetic invertibility, the generator denominator and the metric-map determinant. This regular ($\mathcal J\neq0$) construction is distinct from non-invertible mimetic symmetry. We prove that a same-field sector $K(\phi,\widetilde X)$ preserves the local shift exactly when its scalar density is variationally trivial. For one scalar on a generic branch this requires constant $K$; nonconstant $K(\widetilde X)$ leaves only a global shift. The multi-field extension has local group $\mathbb R^n$ and admits determinant-type topological exceptions for $n\geq d$. Under standard local BV hypotheses, the spectator, its ghost and their antifields form a contractible quartet, so the spectator fibre adds no independent local counterterm, gauge deformation or anomaly class. The constant-map Einstein image is a special two-mode quartic-Horndeski subclass, while regular nonconstant $K$ restores one scalar. Finally, the complete canonical Einstein–scalar one-loop divergence shows that $(\widetilde\Box\phi)^2$ is equation-of-motion redundant at first loop order. The essential $\widetilde X^2$ pole pulls back to a first-derivative counterterm throughout the regular orbit and introduces no additional perturbative mode at $O(\hbar)$ within the stated hyperbolic EFT domain.
Read the full abstract

Astrophysics of Galaxies

The ubiquity of turbulence in the expanding kinematics of the ionized shells of Galactic planetary nebulae

We present an analysis of the residual velocities from a sample of 105 Galactic planetary nebulae (PNe), the largest done to date on this subject. The analysis has been carried out with long-slit, high dispersion echelle spectra. The data were drawn from the San Pedro Mártir Kinematic Catalogue of Galactic Planetary Nebulae. The residual velocity is identified with turbulence in the plasma and is derived by decomposing the emission line profiles into their structural contributors. Turbulence seems pervasive throughout all the PNe in the sample. We find the values for the residual velocities in the sample to be either transonic or slightly supersonic in the ionized environment. When residual velocities of [N II], [O III] and He II in the same PNe are compared, there is a tendency for the residual velocities of the higher ionized ion to be larger by about 5-10 km s$^{-1}$, indicating that the turbulent structure is larger in the inner zones of the PN. We find in general no clear correlation between the residual velocities and other nebular parameters such as morphology, global expansion velocities, ionization degree and binary cores. The only exception is the case of PNe with H-poor ([WR]-type) central stars, where we confirm previous results that have consistently shown higher residual velocities for this group of PNe as compared to those with H-rich central star atmospheres. Turbulence seems to be a localised, random, dissipative process occurring in the inner sections of the shell and may affect its early evolution.
Read the full abstract

Earth and Planetary Astrophysics

ExoNet: Calibrated Multimodal Deep Learning for TESS Exoplanet Candidate Vetting using Phase-Folded Light Curves, Stellar Parameters, and Multi-Head Attention

The discovery of exoplanets at scale has become one of the defining data science challenges in modern astrophysics. NASA's Transiting Exoplanet Survey Satellite (TESS) had catalogued over 7,800 planet candidates by early 2026, yet confirmation stands at fewer than 720. This paper introduces ExoNet, a multimodal deep learning framework that jointly processes phase-folded global and local light curve views alongside stellar parameter features through a calibrated late-fusion architecture combining 1D Convolutional Neural Networks, 8-head Multi-Head Attention over temporal feature maps, and a residual fusion head with post-hoc Temperature Scaling calibration. Trained on 7,585 labeled Kepler Objects of Interest, ExoNet achieves Test AUC = 0.9549 and 86.3% accuracy. Applied to 4,720 verified unconfirmed TESS Planet Candidates with TOI-TIC cross-identification verified against the NASA Exoplanet Archive, the model yields 1,754 high-confidence signals, 52 habitable-zone candidates, and six Earth-sized habitable-zone targets below 1.6 Earth radii. TOI-5728.01 and TOI-6716.01 emerge as the most Earth-like unconfirmed candidates. Full ablation confirms each modality improves AUC. Code and catalog are openly released.
Read the full abstract

High Energy Astrophysical Phenomena

On the contribution of the bow shock pulsar wind nebula PSR J0437-4715 to the observed fluxes of GeV-TeV positrons and antiprotons

The orbital observatories PAMELA and AMS-02 have detected a significant excess in the cosmic ray (CR) positron flux at energies above several tens of GeV. The measured values exceed those expected in models of secondary origin of positrons due to inelastic collisions of CR nuclei with the interstellar matter. This excess may be due to the annihilation or decay of hypothetical dark matter particles or, alternatively, to the contribution of the nearest pulsars. The measured positron-to-antiproton flux ratio is virtually independent of energy from 60 to 400 GeV. This can be understood if the observed local spectra of positrons and antiprotons in the range of tens to hundreds of GeV are formed by the same source. Fast-moving pulsars form pulsar wind nebulae with bow shocks (BSPWNe), which accelerate both the positrons and electrons of the pulsar wind and the hadrons and leptons of galactic CRs from the interstellar medium via the Fermi acceleration mechanism in colliding flows. Such a system can produce identical particle spectra regardless of the site of particle injection. The nearest to the Earth millisecond pulsar PSR J0437-4715 forms a pulsar wind nebula (PWN) with a bow shock observable in optical and ultraviolet wavelengths. This BSPWN is a possible candidate for the main near-Earth “factory” of antiparticles along with the Geminga PWN. Considering PSR J0437-4715, we provide the Monte Carlo simulations of particle acceleration in its BSPWN and the analytical model of anisotropic diffusion in the local interstellar medium. We show that this pulsar's contribution can explain the observed positron flux from 30 GeV to 1 TeV, and simultaneously the antiproton flux at hundreds of GeV with an almost energy-independent positron-to-antiproton flux ratio using $\sim 25\%$ of the PSR J0437-4715 pulsar wind power to accelerate positrons and electrons and re-accelerate antiprotons.
Read the full abstract

Instrumentation and Methods for Astrophysics

Sequential and non-sequential Zemax Dynamic Link Libraries for generating image slicer integral field units

In astronomy, image slicer integral field units (IFUs) are often used in integral field spectrographs to simultaneously record spatial and spectral information. The majority of astronomical instruments, including integral field spectrographs, are designed using the Zemax OpticStudio optical design software. Modeling an image slicer IFU in Zemax traditionally requires using many separate configurations, which is slow, cannot accurately model diffraction, and can prevent one from fully describing their instrument within a single file. This paper presents the implementation of sequential and non-sequential Dynamic Link Libraries (DLLs) that efficiently model image slicer IFUs with a known design. The parameters used to manipulate the surfaces are chosen to match fabrication processes. The DLLs identically reproduce natively transformed surfaces in Zemax and have also been used to replicate the design of SPECTRE, a facility spectrograph with a 36-slice image slicer for the NASA Infrared Telescope Facility. The DLLs also work in transmission and can be used in other applications that require modeling a nearly arbitrary grid of surfaces. In the future, this work may facilitate the creation of an optical design tool for IFUs that starts from basic system requirements.
Read the full abstract

Solar and Stellar Astrophysics

The Effect of External Photoevaporation on the Disk Fraction in M17

A major obstacle to improving models of planet formation is understanding how the local environment influences the lifetime of the disks in which they form. The spread in observed disk lifetimes is caused by effects both observational (e.g., target selection, survey sensitivity) and physical (e.g., disk destruction by internal and external photoevaporation); however, the degree to which each plays a role remains poorly constrained. Isolating the impact of external photoevaporation on the disk lifetime benefits from the inclusion of low-mass ($\lesssim0.5$ M$_{\odot}$) YSOs, for which this effect is most predominant. In this work, we measure the inner disk fraction from JHK excess in the ~6000 M$_{\odot}$, ~1 Myr-old star-forming region M17. Using VLT/HAWK-I, we perform a deep photometric survey of an ~8$^{\prime}\times$8$^{\prime}$ field towards the region. The ~4 times greater sensitivity and ~2-3 times higher resolution than previous surveys of M17 reveal 10,339 sources. We select cluster members using the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) catalog and find a disk fraction of 28$\pm$2%: the first X-ray-selected disk fraction measurement in M17 to include low-mass YSOs, and only the second such measurement in any high-mass star-forming region. After correcting for observational biases, we find no correlation between disk fraction and incident UV flux within M17, likely due to dynamical mixing within the region. However, when compared to other regions of similar age, we find lower disk fractions in regions with higher UV fields, suggesting that external photoevaporation decreases the average disk lifetime.
Read the full abstract

All authors Samuel Millstone, Megan Reiter, Morten Andersen, Thomas J. Haworth, Dominika Itrich, Anna McLeod, Richard J. Parker, Andrew Winter, Peter Zeidler

High Energy Astrophysical Phenomena

Refined Constraints on the Hard X-ray Polarization of the Crab Pulsar and Nebula Derived from an Extended XL-Calibur Dataset

We present updated hard X-ray polarization measurements of the Crab pulsar and nebula obtained with the balloon-borne polarimeter XL-Calibur in the ~19-64 keV energy range. During the flight, intermittent GPS-failure resulted in poorly constrained timing for ~38% of the Crab dataset. By implementing a new phase-recovery method that reconstructs timing during extended GPS-off intervals, phase tag data is recovered for ~95% of the GPS-off dataset, increasing the precision of the phase-resolved analysis. Phase-information for the data is recovered by using the Crab pulsar, with its 33 ms period, as an external timing source. Using a Markov-Chain Monte-Carlo framework to jointly fit phase offsets and frequency derivatives, sufficient phase accuracy is achieved, across multiple periods without GPS for a phase-resolved analysis. This enables inclusion of nearly the full dataset in the polarization study. The polarization degree of the nebular emission is found to be (27.7${\pm}$4.9)% at a polarization angle of 127.2{\deg}${\pm}$5.1{\deg} confirming previous XL-Calibur results and remaining aligned with the Crab's spin axis, consistent with synchrotron emission from the inner nebula. Phase-resolved measurements show that the off-pulse and bridge intervals exhibit a strong polarization, while the pulsar peaks, although weakly constrained, remain in agreement with the softer-energy trends of IXPE. These findings reinforce a scenario in which hard X-ray emission arises primarily in the nebular torus and wind regions. The successful recovery of precise phase tagging from GPS-off data demonstrates the capacity to use the pulsar as an external clock even in the case of sparsely populated data.
Read the full abstract

All authors Matthew G. Baring, Jacob Casey, Sohee Chun, Ephraim Gau, Tomohiro Hakamata, Kun Hu, Daiki Ishi, Fabian Kislat, Mózsi Kiss, Merlin Kole, Henric Krawczynski, Haruki Kuramoto, Lindsey Lisalda, Bingkun Liu, Yoshitomo Maeda, Hironori Matsumoto, Shravan Vengalil Menon, Takuya Miyazawa, Kaito Murakami, Takashi Okajima, Mark Pearce, Brian Rauch, Kentaro Shirahama, Sean Spooner, Hiromitsu Takahashi, Sayana Takatsuka, Yuusuke Uchida, Varun, Andrew Thomas West

High Energy Astrophysical Phenomena

Joint Detection and Characterization of the Standing Accretion Shock Instability for Core-Collapse Supernovae with cWB XP

The most sensitive to-date multimessenger detection of the standing accretion shock instability in real interferometric data is presented, which quantitatively identifies the presence of the SASI in core-collapse supernovae using neutrino and gravitational-wave (GW) signals. In the GW channel, the coherent WaveBurst (cWB) software on its version XP is implemented, among with real LIGO data from the O3 and O4 observing runs. With this, a more accurate estimation of parameters, such as the central frequency and signal duration, is obtained for both sets of data. The SASI identification probability versus false alarm rates is presented in the form of Receiver Operating Characteristic (ROC) curves. For O3, the new study for the combined GW and neutrino detection condition, labeled as $x + y$, shows an identification probability (previous best results from Lin et al. [1]) of 1 (1), 0.90 (0.70) and 0.37 (0.34) at 1, 5 and 10 kpc for a false identification probability of 0.10. On the other hand, using O4 shows that the GW channel by itself is sensitive enough to provide almost perfect identification probability scores, with identification probability values of 1, 0.99 and 0.97 for a false identification probability of 0.01 at 1, 5 and 10 kpc, respectively.
Read the full abstract

All authors Vicente Sierra, Zidu Lin, Michelle Zanolin, Claudia Moreno, Javier M. Antelis, Marek J. Szczepańczyk

High Energy Astrophysical Phenomena

Detectability of Gravitationally Lensed Kilonovae in the Rubin LSST

Identification and characterisation of Kilonovae (KNe) can be instrumental in improving our understanding of cosmology and astrophysics. However, their detection poses unique challenges due to rarity and faintness. Upcoming telescopes, with their deep imaging capabilities and wide field-of-views, will provide a unique opportunity to observe these rare and faint transients. The Rubin Legacy Survey of Space and Time (LSST) will generate a deluge of data, making it essential to deploy fast, efficient methods for identifying genuine KNe, especially when they are gravitationally lensed. To address this, we simulate realistic populations of both unlensed and lensed KNe in the six LSST bands. Comparing with the Type Ia Supernovae, we find that the KNe color evolution is more rapid and the two separate out when their colors are compared at two epochs. Since the mergers of compact binaries are probable progenitors of KNe, the KNe properties may be affected by the delay time distribution (DTD) of the mergers, which is dictated by the minimum delay time ($\tau$) and power-law slope. For longer $\tau$ and shallower slopes, we find an increased rate of detectable KNe in LSST. We generate the first statistically realistic lensed KNe population for different DTDs and find that the rate of detectable lensed KNe increases for DTDs with longer $\tau$ for a fixed slope. We further note that an AT2017gfo-like event at a redshift of 0.5~(1.0) needs magnification of at least 5~(44) to be detectable in LSST.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Probing the large-scale structure with 21cm-galaxy cross-bispectrum: Estimates from simulations and forecasts for upcoming cosmological surveys

The redshifted 21cm signal from the post-reionization epoch is highly non-Gaussian; thus, higher-order statistics, such as the bispectrum, are required to extract this non-Gaussian information. However, high signal-to-noise ratio (SNR) detection of the 21cm auto-bispectrum will be hindered by the presence of residual systematics. Cross-correlating the 21cm signal with galaxies offers a promising path to suppress this uncertainty from residual systematics and potentially increase the SNR. We present a comprehensive analysis of the HI-galaxy cross-bispectrum using the predictions of theoretical galaxy evolution models defined on large cosmological volumes. Our analysis includes the cross-bispectrum for different triangle sizes and shapes, as well as for different combinations of the HI and galaxy fields. We forecast the detectability of the 21cm-galaxy cross-bispectrum at redshift $z\approx1$ with the Euclid-like galaxy survey and SKA-Mid observations in both interferometric and single-dish modes of the survey. We find that the 21cm-galaxy cross-bispectrum shows enhanced detectability compared to the 21cm auto-bispectrum for all unique triangles in the interferometric mode of observations. We forecast a 10$\sigma$ detection of the cross-bispectrum for squeezed-limit triangles and a 100$\sigma$ detection for all shapes combined for scales $0.2~\text{Mpc}^{-1}\leq k_1 \leq 0.9~\text{Mpc}^{-1}$ with 100 hours of SKA-Mid observations per pointing. However, the detectability of the cross-bispectrum for large scales ($k_1 < 0.1~\text{Mpc}^{-1}$), which is accessible with the single-dish mode of the survey, is limited by cosmic variance. Additionally, the signal loss due to foreground removal further suppresses the detectability. Our analysis presents a first step toward an end-to-end analysis pipeline for the future 21cm-galaxy cross-bispectrum observations.
Read the full abstract

All authors Leon Noble, Suman Majumdar, Matteo Viel, Fabio Fontanot, Gabriella De Lucia, Abinash Kumar Shaw, Marta Spinelli, Mohd Kamran, Lizhi Xie, Michaela Hirschmann